# Impact of karst areas on runoff generation, lateral flow and interbasin groundwater flow at the storm-event timescale

Martin Le Mesnil, Roger Moussa, Jean-Baptiste Charlier, Yvan Caballero

<figure><img src="image_1.png" /></figure>

## To cite this version:

Martin Le Mesnil, Roger Moussa, Jean-Baptiste Charlier, Yvan Caballero. Impact of karst areas on runoff generation, lateral flow and interbasin groundwater flow at the storm-event timescale. Hydrology and Earth System Sciences, 2021, 25 (3), pp.1259 - 1282. 10.5194/hess-25-1259-2021. hal-03169475

HAL Id: hal-03169475

https://brgm.hal.science/hal-03169475

Submitted on 15 Mar 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

# Impact of karst areas on runoff generation, lateral flow and interbasin groundwater flow at the storm-event timescale

Martin Le Mesnil<fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Roger Moussa<fnref n="1" />, Jean-Baptiste Charlier<fnref n="2" />,<fnref n="3" />, and Yvan Caballero<fnref n="2" />,<fnref n="3" />

**Correspondence:** Martin Le Mesnil (martin.le-mesnil@inrae.fr)

Received: 15 May 2020 – Discussion started: 1 July 2020

Revised: 1 February 2021 – Accepted: 7 February 2021 – Published: 12 March 2021

**Abstract.** Karst development influences the hydrological response of catchments. However, such an impact is poorly documented and even less quantified, especially over short scales of space and time. The aim of this article is thus to define karst influence on the different hydrological processes driving runoff generation, including interbasin groundwater flow (IGF) for elementary catchments at the storm-event timescale. IGFs are estimated at the scale of the river reach, by comparing inlet and outlet flows as well as the effective rainfall from the topographic elementary catchment. Three types of storm-event descriptors (characterizing water balance, hydrograph shape and lateral exchanges) were calculated for the 20 most important storm events of 108 stations in three French regions (Cévennes Mountains, Jura Mountains and Normandy), representative of different karst settings. These descriptors were compared and analysed according to catchment geology (karst, non-karst or mixed) and seasonality in order to explore the specific impact of karst areas on water balance, hydrograph shape, lateral exchanges and hydrogeological basin area. A statistical approach showed that, despite the variations with study areas, karst promotes (i) higher water infiltration from rivers during storm events, (ii) increased characteristic flood times and peak-flow attenuation, and (iii) lateral outflow. These influences are interpreted as mainly due to IGF loss that can be significant at the storm-event scale, representing around 50 % of discharge and 20 % of rainfall in the intermediate catchment. The spatial variability of such effects is also linked to contrasting lithology and karst occurrence. Our work thus provides a generic framework for assessing karst impact on the hydrological response of catchments to storm events; moreover, it can anal-

yse flood-event characteristics in various hydro-climatic settings and can help with testing the influence of other physiographic parameters on runoff generation.

## 1 Introduction

Understanding runoff generation requires a good knowledge of the different processes involved in catchment response to rainfall events, i.e. how precipitation is converted into underground, subsurface or surface flow. These processes are affected by several factors, such as the thickness and nature of soil, land use, hydrological conditions, or geology. While most of these can be documented or measured, it may be difficult to define the role of geology in a comprehensive way, especially when underground drainage networks are involved, as in karst areas. Karstification, the result of carbonate-rock dissolution, promotes infiltration and groundwater flow through enlarged fissures and voids (Bakalowicz, 2005), locally dramatically reducing drainage density on the surface and thus affecting the hydrological response of a catchment. These groundwater flows can be considered interbasin groundwater flows (IGF) if they do not return to surface within the considered catchment. IGF amplitude and sign can be linked to the area of the hydrogeological basin, which can be different from the area of the topographic one.

Karst impacts on flood processes are mostly documented through case studies. As an example, Zanon et al. (2010) showed that during a flash flood in 2007 in Slovenia, a karst area reduced flooding, which was more important in a non-karst neighbouring zone, receiving less precipitation. Like-

<footnote n="1">¹LISAH, Univ. Montpellier, INRAE, IRD, Institut Agro, Montpellier, France</footnote>

<footnote n="2">²BRGM, Univ. Montpellier, Montpellier, France</footnote>

<footnote n="3">³G-eau, INRAE, CIRAD, IRD, AgroParisTech, Institut Agro, BRGM, Montpellier, France</footnote>
wise, Delrieu et al. (2005) observed, for an exceptional storm event in 2002, lower runoff coefficient values for the karstic catchment compared to the hard-rock catchment in the eastern zone of the Cévennes Mountains. De Waele et al. (2010) and Charlier et al. (2015, 2019) determined that, depending on the location on the river profile, karst areas could result in streamflow losses or gains due to the high spatial variability of the hydrogeological karst features. Other frequently described processes are groundwater rising, leading to reduced infiltration and important surface runoff (López-Chicano et al., 2002; Bonacci et al., 2006), and backflooding or sinkhole flooding due to conduit constriction (Maréchal et al., 2008; Bailly-Comte et al., 2009).

The diversity of observed processes during storm events in karst catchments does not allow for drawing a straightforward analysis on the control of karst in flood runoff generation. For the purpose of understanding the mechanisms involved in this control, there is a need for regionalized studies, covering a large-scale analysis of karst impact over short time periods when the catchment reacts after storm events. It is reasonable to think that karst can alternately increase or decrease storm impacts, depending on its capacity to infiltrate precipitation or to release stored water, i.e. depending on the direction of IGF it promotes. Despite the early conceptualization of IGF (Eakin, 1966), its major role in karst hydrological processes is tackled by very few studies (Le Moine et al., 2007; Lebecherel et al., 2013). Some authors tried to improve model capacities to reproduce karst-based IGF, such as Nguyen et al. (2020) with SWAT, Le Moine et al. (2008) with GR4J or Scanlon et al. (2003) comparing a distributed and a lumped model. Nevertheless, those studies that are dedicated to the improvement of model performance are not devoted to describe and understand all flood components in karst catchments.

On the one hand, most studies including karst system descriptors are based on a purely hydrogeological point of view and are very integrative, as they tend to characterize the karst aquifer as a whole, by analysing daily spring discharge. Gárfias-Soliz et al. (2009) found that system memory, response time and mean input-output delay are relevant indicators for karstification, in addition to a necessary consideration of the structural complexity and heterogeneity of the lithology. Hartmann et al. (2013), using 10 system signatures, performed a model parameter sensitivity analysis to investigate their links with hydrological processes on five European and Middle Eastern karst sites. Basha et al. (2020) proposed six recession curve equations for the classification of karst aquifers, depending on their flow characteristics. On the other hand, some studies accounting for a spatialization of catchments focus on low-flow issues and surface-water-groundwater interaction (e.g. Covino et al., 2011; Mallard et al., 2014). Moreover, most regionalization studies tend to spatialize annual indices (Sivapalan et al., 2011) or model parameters (Parajka et al., 2005; Oudin et al., 2008), and they

usually exclude catchments with identified IGF (Merz and Blöschl, 2004) or karst areas (e.g. Laaha and Blöschl, 2006).

Regional spatial analyses need to be based on reliable data at the highest resolution available. For this purpose, the scale of the elementary catchment – i.e. subdivision of a basin following available gauging stations – appears to be the best resolution for long-term monitoring. Elementary catchment can be either the drained area of a headwater catchment controlled by a gauging station or the drained area between two gauging stations (intermediate catchment). When considering surface and groundwater components, the delineation method of elementary catchments is questionable (topographic versus hydrogeological boundaries). Despite the importance of groundwater processes in karst areas, topographic catchment delineation remains a more robust reference for several methodological reasons. First, IGF can be defined as groundwater flow crossing topographic divides, as this concept emerged with the evidence of certain groundwater systems extending beyond the limits of valleys (Eakin, 1966). A perfectly delineated groundwater basin would then show IGF equal to zero. For this reason, studies related to IGF often use the topographic catchment spatial reference (Genereux et al., 2005; Schaller and Fan, 2009; Bouaziz et al., 2018; Nguyen et al., 2020; see also a synthesis in Fan, 2019). Second, although groundwater contributes to flood flow in karst catchments, the surface runoff component should not be completely discarded. Exclusive consideration of hydrogeological catchments could thus lead to wrong surface contribution assessment depending on their surface drainage network. Third, as some groundwater flows are aligned with the main surface drainage axis, hydrogeological catchments would encompass the whole river, making it impossible to study the spatial variability of parameters along the river, at the elementary catchment scale. Finally, topographic delineation is reliable and easily reproducible, while groundwater delineation is characterized by a strong uncertainty and variability in karst areas.

This article aims at providing a new methodology to characterize the spatial variability of karst influence on hydrological processes affecting runoff generation, including IGF, at the storm-event timescale. IGFs are estimated at the scale of the river reach, by comparing inlet and outlet flows as well as the effective rainfall from the topographic elementary catchment. The present study complements the previous work by Le Mesnil et al. (2020), which described the role of karst areas using annual water-budget indicators at the elementary catchment scale. Here, descriptors are calculated for major storm events at 108 stations in three areas in France (Cévennes Mountains, Jura Mountains and Normandy) with different karst settings. The descriptors are of three types: water balance, hydrograph shape and lateral exchange. Water-balance descriptors are obtained from discharge and precipitation depths, analysing the respective importance of the different flows during storm events. They help with understanding how catchments transform precipitation
into surface and underground flow. Such descriptors are of great interest to assess the spatial variability of catchment hydrological response (Sivapalan et al., 2011). Hydrograph-shape descriptors are derived from catchment hydrographs recorded at inlet and outlet stations during storm events. They describe hydrological processes (Raghunath, 2006) and, when analysed on successive stations, help with characterizing flood-wave routing. Lateral-exchange descriptors are based on lateral hydrographs, simulated with the diffusive wave equation (DWE; Moussa, 1996) applied between two gauging stations (for intermediate catchments only). They provide information on lateral inflow and outflow at the elementary catchment scale. Lateral flow is mainly the result of IGF, effective rainfall, variations in aquifer storage and overbank phenomena. A particular analysis is performed on IGF, which are expressed in depth as well as in theoretical catchment area variation, and seasonal influence is discussed.

The three types of descriptors are compared and analysed according to the catchment geology type (classified as “karst”, “non-karst” or “mixed”), in order to explore its impact on runoff generation processes. The paper thus provides a framework for assessing the impact of a given physiographical parameter on the hydrological response of catchments to storm events.

## 2 Methodology

### 2.1 General methodology

We calculated 15 descriptors (five for each of the three types) of catchment response to storm events and assessed their variability as a function of karst occurrence. To this end, we grouped elementary catchments into three different geology types, based on relative areas of their main geological formations. Catchments underlain by only karstified rock are in the karst group (K), whereas catchments containing only non-karstified rock are in the non-karst group (NK). Any catchment with a combination of both karstified (> 10%) and non-karstified rocks is in the mixed group (M). The karstifiable nature of rocks was assessed with the BDLISA (Base de Donnée des Limites des Systèmes Aquifères) database (Sect. 3.1.2). Descriptors were calculated for the 20 strongest storm events for all catchments (Sect. 2.2). The values obtained for each group (K, M, NK) were compared with each other to assess karst influence on them. Then, for K catchments only, the three study areas (Cévennes, C; Jura, J; Normandy, N; see Sect. 3.2 for area descriptions) were compared to assess the area-specific nature of karst influence. Our framework was kept as generic as possible to propose an approach that would be easily adaptable to another investigation field. Consequently, the karst-catchment classification can be replaced by any other physiographical typology.

Descriptors are complementary but not necessarily independent from each other. They are of three types and were

chosen to provide relevant information on different processes:

- Water-balance descriptors show the respective importance of the different flows occurring during storm events and allow for understanding how catchments transform rainfall into surface and underground flow. Such descriptors are volumes (expressed as depth [mm]) or volume ratios. They can be calculated at headwater-catchment or intermediate-catchment outlets, with the latter involving subtracting inlet flow from outlet flow. This was applied to the 108 elementary catchments.

- Hydrograph-shape descriptors describe the dynamics of storm events and flood-wave routing. They combine peak-flow variation, characteristic times and flood-wave celerity. Characteristic times can be obtained from any measured hydrograph, whereas peak-flow variation and celerity are evaluated between inlet and outlet hydrographs (for intermediate catchments only), considering catchments with only one inlet station and some intermediate catchments having several inlets when covering the confluence of streams. This was applied to all 108 elementary catchments for characteristic times and to the 36 intermediate catchments with only one inlet regarding peak-flow variation and celerity.

- Lateral-exchange descriptors provide information on the dynamics of lateral inflow and outflow affecting an elementary catchment reach, as well as on the respective contributions of channel diffusivity and lateral exchanges to peak-flow variations. This analysis is based on lateral hydrographs, simulated with the DWE, that are applied to intermediate catchments with one inlet. Lateral exchanges may be a combination of effective rainfall ($P_{eff}$) over the elementary catchment, IGF, aquifer storage variation and overbank phenomena. Overbank flow is considered a minor process assuming that the overflow water returns to the river after a relatively short time during the recession. Thus, a water balance including $P_{eff}$ allows for discussing the importance of IGF, along with aquifer storage variation. This was applied to the 36 intermediate catchments with only one inlet for lateral exchange dynamics and to all 108 elementary catchments for IGF assessment.

### 2.2 Flow assessment at the elementary catchment scale

The 15 descriptors were calculated on elementary catchments for 20 selected strongest storm events. For the sake of representativeness, this selection was based on both rainfall and streamflow records. From the available data time series spanning several decades (see Sect. 3.1 for more details), the 10 strongest precipitation and 10 strongest streamflow events were extracted. Care was taken not to select the same events via the two extraction methods.
Streamflow is measured at several gauging stations along a given river, defining elementary catchments. For the most upstream station of a river, the elementary catchment corresponds to the ordinary topographic catchment. Otherwise, the elementary catchment is an intermediate one, covering the portion of the basin drained between two gauging stations (Fig. 1). In the case of intermediate catchments, streamflow is calculated following Eq. (1) as the difference between outlet flow ($Q_O$) and inlet flow ($Q_I$). A hydrograph decomposition is operated for separating the quick- (storm flow $Q_S$) and slow-flow (underground flow $Q_U$) components, as shown in Eqs. (2) and (3) (see Appendix A for more details). The $Q_S$ and $Q_U$ components of the intermediate catchment are also obtained from the difference between inlet and outlet flow components as shown in Eqs. (4) and (5). For these intermediate catchments, the $Q$, $Q_S$ and $Q_U$ values can thus be negative; note also that the associated uncertainties are twice higher than those of measured flows ($Q_O$ and $Q_I$).

$$Q = Q_O - Q_I \qquad (1)$$

$$Q_I = Q_{I, U} + Q_{I, S} \qquad (2)$$

$$Q_O = Q_{O, U} + Q_{O, S} \qquad (3)$$

$$Q_S = Q_{O, S} - Q_{I, S} \qquad (4)$$

$$Q_U = Q_{O, U} - Q_{I, U} \qquad (5)$$

In addition to inlet and outlet flows, a lateral flow ($Q_L$) is calculated. This flow is mainly a combination of IGF, effective rainfall ($P_{eff}$) over the elementary catchment, aquifer storage variation and overbank phenomena (Eq. 6). Overbank phenomena are neglected here (see Sect. 2.1). The $Q_L$ hydrograph is simulated using an inverse model, solving the DWE accounting for lateral flow between inlet and outlet stations, assuming a uniform lateral distribution of exchanges along the river reach (see Appendix B for more details).

$$Q_L = \text{IGF} + P_{\text{eff}} + \delta \qquad (6)$$

Finally, after an estimation of $P_{eff}$ (Appendix C), all measurable inflow ($Q_I$ and $P_{eff}$) and outflow ($Q_O$) are known. Hence, IGF values (groundwater flowing inside or outside the elementary catchments delimited by topographic boundaries) can be calculated for each storm event, along with the aquifer storage (noted $\delta$). Figure 1 shows all considered flows and their corresponding hydrographs at the spatial scale of an elementary catchment and the timescale of a storm event.

## 2.3 Descriptors

### 2.3.1 Water-balance descriptors

For each storm event, volumes of the different flows in an elementary catchment are calculated (Fig. 1). Total discharge volume, noted $V_O$, is calculated as

$$V_Q = V_O - V_I = (V_{O, S} + V_{O, U}) - (V_{I, S} + V_{I, U}), \qquad (7)$$

with $V_I$ and $V_O$ being the volumes of total streamflow (expressed as water depths in mm) and $V_{I, S}$, $V_{I, U}$, $V_{O, S}$ and $V_{O, U}$ being those of the quick- and slow-flow components of inlet $Q_I$ and outlet $Q_O$, respectively. $V_P$ is the volume of precipitation falling on the elementary catchment during the storm event. With the values being dependent on the catchment surface, they are normalized by the topographic area (A) of the considered catchment and expressed in water depth [mm]. From these volumes, five water-balance descriptors are calculated (see Appendix D for equations).

- $V_S$, event storm flow [mm];

- $V_U$, event base flow [mm];

- $R_C$, event runoff coefficient [dimensionless], calculated as total event streamflow divided by event rainfall;

- $R_C$, s, event storm runoff coefficient [dimensionless], calculated as event storm flow divided by event rainfall;

- $S_C$, event storm-flow coefficient [dimensionless], calculated as event storm flow divided by total event streamflow.

### 2.3.2 Hydrograph-shape descriptors

For each storm event, five more descriptors were calculated, characterizing hydrograph morphology and storm-event dynamics. Characteristic times are of great interest in storm hydrology and constitute a widely used framework, allowing convenient catchment and event comparisons (Bell and Om Kar, 1969). Here, we use three of them. The time constant of the rising limb $T_{Ri}$ corresponds to the time duration needed for streamflow to increase from half of peak flow (0.5$Q_{O,X}$) to peak flow ($Q_{O,X}$); $T_{08}$ is the time duration for $Q_O(t)>0.8Q_{O,X}$, important in terms of operational flood management. $T_{Re}$ is the characteristic recession time, obtained from the linearization method of Maillet (1905), which approximates the recession curve by an exponential function as shown in Eq. (8).

$$Q_{\text{Maillet}}(t) = q \cdot e^{-\beta t} = q \cdot e^{-t/T_{\text{Re}}}, \qquad (8)$$

with $\beta$ the recession coefficient and $q$ the discharge at the beginning of the recession phase. The $T_{Re}$ time constant is the duration needed for streamflow to decrease by a factor $e \approx 2.7$ during recession.

In the case of intermediate catchments, the peak-flow variation ($\Delta$) is calculated as the difference between outlet station peak flow ($Q_{O,X}$) and inlet station peak flow ($Q_{I,X}$), normalized by rainfall. A parameter $C_G$ is also calculated, used as the celerity for applying the DWE. $C_G$ is equal to the river reach length ($l$) divided by the elapsed time ($T_G$) between $Q_I$ and $Q_O$ gravity centres ($G_I$ and $G_O$). The five hydrograph-shape descriptors (see Appendix D for formulae) are the following.
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Main hydrological flow types at the elementary catchment and storm-event scales, with corresponding measured and simulated hydrographs. Flow definitions are given in Eqs. (1) to (5). The corresponding volumes, integrated for the storm event, are noted *V* with the same indices (see Appendix D for symbols). Lateral flow *Q*<sub>L</sub> is a combination of IGF, effective rainfall (1), aquifer storage variation (2) and overbank phenomena (3, neglected here).</figcaption></figure>

- $\Delta/V_p$, peak-flow variation normalized by rainfall [$s^{-1}$];

- $T_{Ri}$, event time duration of rise [h];

- $T_{08}$, event time duration to 80% of peak flow [h];

- $T_{Re}$, event characteristic recession time as per the approximation by Maillet (1905) [h];

- $C_G$, celerity based on elapsed time between hydrograph gravity centres [$m s^{-1}$].

### 2.3.3 Lateral-exchange descriptors

For catchments with one inlet station, a lateral-flow hydrograph is simulated, using the solution of the inverse problem of the DWE assuming that the lateral flow is uniformly distributed along the channel (Moussa, 1996; see Appendix B). The DWE has two free parameters: the celerity $C$ [$m s^{-1}$] and the diffusivity $D$ [$m^2 s^{-1}$] of the flood wave. First, $Q_I$ is routed using the DWE without lateral flow for given values of $C$ and $D$, simulating a theoretical outflow without lateral exchange (noted $Q_{I,R}$ in Fig. 3). The inverse problem supposes known $Q_I$, $Q_O$, $C$ and $D$ values, simulating the lateral flow $Q_L$. Sensitivity analysis of the DWE to the two parameters is largely available in the literature, showing that the DWE is more sensitive to parameter $C$ than $D$ (Moussa and Bocquillon, 1996; Cholet et al., 2017; Charlier et al., 2019). Here, $C$ was assumed to be equal to $C_G$, and we used a fixed value for $D = 10000 m^2 s^{-1}$ (for medium size catchments of 100 to 500 km$^2$; Moussa and Bocquillon, 1996; Todini, 1996) in a matter of parsimony, which is acceptable as the model is much more sensitive to $C$ than to $D$.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Storm hydrograph with characteristic times and discharge values. $Q_O$, $Q_I$, and $Q_{I,R}$ are peak flows of outlet and inlet station hydrographs $Q_O$ and $Q_I$, respectively. $T_G$ is the elapsed time between the corresponding gravity centres $G_O$ and $G_I$.</figcaption></figure>

The DWE solution was validated experimentally under controlled conditions (Moussa and Majdalani, 2019), and it has also been implemented on natural karst catchments (Charlier et al., 2015, 2019; Cholet et al., 2017).

From the simulated lateral hydrograph, five descriptors were calculated. First, the peak-flow-variation descriptor $\Delta$ was decomposed into two components ($\Delta_D$ and $\Delta_L$), which are the respective contributions of channel diffusivity and lateral exchanges (Charlier et al., 2019). $\Delta_D$ and $\Delta_L$ are also normalized by $V_p$. In the case of a reach without lateral ex-
<figure><img src="image_7.png" /><figcaption>Figure 3. Theoretical examples of simulated lateral hydrographs $Q_L$, with corresponding values of peak-flow variation $\Delta$ and its two components $\Delta_D$ and $\Delta_L$. Modified from Charlier et al. (2019). $Q_{O,X}$ and $Q_{I,X}$ are peak flows of the outlet and inlet station hydrographs $Q_O$ and $Q_I$, respectively. $Q_{I,R}$ is the theoretical outlet flow without lateral exchange (routed inflow by DWE, with C and D parameters).</figcaption></figure>

changes, $\Delta$ is negative and only due to channel diffusivity (equal to $\Delta_D$, Fig. 3a). In the case of a gaining catchment, $\Delta$ can be positive or negative, with compensating contributions of channel diffusivity and lateral inflow (Fig. 3b). In the case of a losing catchment, $\Delta$ is negative, with cumulated contributions of channel diffusivity and lateral outflow (Fig. 3c). Second, volumes of lateral inflow and outflow ($V_{L+}$ and $V_{L-}$) were calculated from the lateral hydrograph and normalized by event rainfall. Note that $Q_L$ can be successively positive and negative within a single storm event (Fig. 1).

As simulated lateral hydrographs are based on the difference between discharge volumes at the inlet and outlet of catchments, they also integrate precipitation on the elementary catchment during storm events. In order to focus on groundwater exchanges and suppress the influence of effective rainfall ($P_{eff}$), this latter has to be estimated. It was done using three different modelling approaches, providing a range for this component characterized by major uncertainties. The three approaches respectively based on the methods of Thornthwaite (1948), Dingman (2015), and the GR model (Edijatno et al., 1999), are described in Appendix C. The final value is the mean of those obtained from the three methods.

The remaining lateral-exchange term combines IGF and potential aquifer storage variation ($\delta$). $\delta$ can either be positive, corresponding to an aquifer recharge, or negative, corresponding to an aquifer draining. Aquifer draining is less likely during storm events, as important rainfalls generally occur. In case of an aquifer recharge, and as our analysis is performed on the whole storm-event period (including the entire recession), a substantial part of the infiltrated water should be released, either inside the considered elementary catchment or outside of it. In the first case, the considered amount of water is accounted for in the $V_O$ term and will not

influence the lateral-exchange term. In the second case, the released water actually constitutes an IGF. For these reasons, we use the term IGF* for the obtained water-balance term. For each storm event, IGF*, combining IGF and potential aquifer storage variation ($\delta$), is calculated as in Eq. (9):

$$IGF^* = IGF + \delta = V_Q - V_{P_{eff}} = V_O - V_I - V_{P_{eff}} \quad (9)$$

The five lateral exchange descriptors (see Appendix D for formulae) are the following:

- $\Delta_D/V_P$, part of peak-flow variation due to channel diffusivity, normalized by rainfall [$s^{-1}$];

- $\Delta_L/V_P$, part of peak-flow variation due to lateral exchanges, normalized by rainfall [$s^{-1}$];

- $V_{L+}/V_P$, event lateral inflow, normalized by rainfall [dimensionless];

- $V_{L-}/V_P$, event lateral outflow, normalized by rainfall [dimensionless];

- $IGF^*/V_P$, event interbasin groundwater flow including potential aquifer storage variation ($\delta$), normalized by rainfall [dimensionless].

## 2.4 Statistical approach

Once all descriptors have been calculated, a statistical analysis is performed for comparative purposes. For each descriptor, the obtained values are grouped in different samples, by (i) geology type (K, M, NK) and (ii) study area (C, J, N) for karst catchments only (K geology type). This allows for characterizing the impact of karst areas on the hydrological response and provides additional information on regional specifics of this impact.
The results are presented as boxplots to discuss how the distribution of the descriptors varies for all samples. Then, statistical tests assess the significance of the results. Twin-sample *t* tests are performed successively on K vs. NK, M vs. NK and M vs. K catchments and – only for K catchments – on C vs. J, C vs. N and J vs. N. Since it is not assumed that the two data samples are from populations with equal variances, the test statistics under the null hypothesis have an approximate Student's *t* distribution with a number of degrees of freedom given by Satterthwaite's approximation (Satterthwaite, 1946). This test provides a decision for the null hypothesis that the data in paired tested samples come from independent random samples from normal distributions with equal means and equal but unknown variances. The result is 1 if the test rejects the null hypothesis at the 5% significance level and 0 otherwise.

# 3 Data sets and study areas

## 3.1 Data sets

### 3.1.1 Temporal data

Temporal data used in this paper are the following:

- Hourly streamflow data from the French public streamflow database "Banque Hydro" (http://www.hydro.eaufrance.fr/, last access: 10 April 2020), managed by the French Regional Environment Directions (Direction Régionale de l'Environnement, de l'Aménagement et du Logement, DREAL). These data are produced by interpolation of water-depth measurements at a variable infra-hourly time step and converted into streamflow values.

- Hourly rainfall data are from Comephore (Tabary, 2007) covering 1997 to 2005 and Antilope (Champeaux et al., 2011) from 2006 to 2018. Both data sets are measurement reanalyses edited by the French public meteorological service Météo-France (http://www.meteofrance.fr/, last access: 3 February 2020). They have a 1 km² spatial resolution and consist of radar rainfall measurements, calibrated to fit the data from surface precipitation gauges.

- Daily potential evapotranspiration depths are from "Safran" (Système d'Analyse Fournissant des Renseignements Atmosphériques à la Neige; Vidal et al., 2010), edited by the French meteorological service (Météo-France). They are used for estimating effective rainfall (Appendix C).

Hourly rainfall data are available from 1997 onward, and hourly streamflow data periods vary depending on the catchments. Table 1 shows periods of availability of both data sets

for 11 hydrographic catchments in three areas of France, covering a total area of almost 25 000 km². Attention was paid to use only validated streamflow data from gauging stations with a hydrological significance, i.e. not influenced by human activities such as damming or pumping. Nevertheless, a 10% uncertainty is associated with streamflow data, which can be higher during extreme storm events due to uncertain rating equations or measuring ranges.

### 3.1.2 Spatial data

Spatial data used in this paper are the following:

- Boundaries of topographic catchments, from the French National Watershed Database (Base Nationale des Bassins Versants, BNBV). It was edited by the French Central Service for Hydrometeorology and Support on Flood Forecasting (Service Central d'Hydrométéorologie et d'Appui à la Prévision des Inondations, SCHAPI) and the French Research Institute for Agriculture, Food, and the Environment (Institut National de Recherche en Agriculture, Alimentation et Environnement, INRAE).

- BDLISA database (https://bdlisa.eaufrance.fr/, last access: 3 April 2020), which describes the properties of hydrogeological entities and aquifers in France, edited by the French Geological Survey (Bureau de Recherches Géologiques et Minières, BRGM).

- Map of available soil-water capacity from INRAE (Le Bas, 2018), used for estimating effective rainfall (Appendix C).

## 3.2 Study areas

### 3.2.1 General settings

The previously described methodology was applied to three areas in France, including 11 hydrographic basins and representing a total area of 25 000 km². The three areas are totally or partially karstified, with different geological and hydrometeorological settings. From south to north, they belong to the regions of Cévennes Mountains, the Jura Mountains and Normandy. Figure 4 shows the locations of the 11 hydrographic basins, on the map of French karst aquifers (Fig. 4a), and the gauging station networks with geology type of the elementary catchments on maps Fig. 4b, c and d.

In the Cévennes Mountains (Fig. 4b), six hydrographic basins were studied, including 51 gauging stations. They are mostly so-called binary karst basins, with head catchments on hard rock receiving around 1500 mm yr⁻¹ precipitation, and median and downstream parts underlain by limestone plateaus, with around 1000 mm yr⁻¹ rainfall.

The Jura Mountains region corresponds to the Doubs river basin, a few kilometres upstream from its confluence with
<table><caption>**Table 1.** Studied catchments and associated available data.</caption><thead><tr><th></th><th>Study<br>zone area<br>[km²]</th><th>Gauging<br>stations</th><th>Median<br>gauged<br>area [km²]</th><th>Time<br>series<br>length</th></tr></thead><tbody><tr><td>Ardèche</td><td>2257</td><td>9</td><td>193</td><td>1996–2018</td></tr><tr><td>Cèze</td><td>1048</td><td>6</td><td>192</td><td>2002–2018</td></tr><tr><td>Gardons</td><td>1853</td><td>10</td><td>137</td><td>2008–2018</td></tr><tr><td>Vidourle</td><td>772</td><td>4</td><td>182</td><td>2009–2018</td></tr><tr><td>Hérault</td><td>2203</td><td>9</td><td>223</td><td>2007–2018</td></tr><tr><td>Tarn</td><td>2145</td><td>13</td><td>64</td><td>1984–2018</td></tr><tr><td>Total Cévennes</td><td>10 300</td><td>51</td><td>172</td><td>13 years (median)</td></tr><tr><td>Doubs</td><td>7400</td><td>39</td><td>121</td><td>1998–2018</td></tr><tr><td>Total Jura</td><td>7400</td><td>39</td><td>121</td><td>21 years</td></tr><tr><td>Iton</td><td>1048</td><td>2</td><td>524</td><td>1999–2018</td></tr><tr><td>Risle</td><td>1803</td><td>5</td><td>84</td><td>2001–2018</td></tr><tr><td>Touques</td><td>800</td><td>5</td><td>106</td><td>2010–2018</td></tr><tr><td>Dives</td><td>879</td><td>6</td><td>113</td><td>2009–2018</td></tr><tr><td>Total Normandy</td><td>4500</td><td>18</td><td>110</td><td>12 years (median)</td></tr><tr><td>Total all basins</td><td>22 200</td><td>108</td><td>175</td><td>14 years (median)</td></tr></tbody></table>

the Saône river, which includes 39 gauging stations. Out-crops mostly consist of extensively karstified Jurassic limestone and marl, except in the extreme northern and western parts of the region. Precipitation follows a strong elevation gradient, with annual values ranging from 1700 mm on up-stream catchments at heights of 1400 m (a.s.l) to 1200 mm at the outlet at an elevation of around 200 m (a.s.l).

In Normandy, four hydrographic basins covered 18 gauging stations. The two eastern basins are tributaries of the Seine, with the other two being coastal basins. The climate is maritime, and annual rainfall ranges from 700 to 1000 mm. Rivers in the eastern part of the area drain chalky limestone with karst covered by clay. The midwestern zone is underlain by Jurassic limestone, corresponding to the western border of the Paris Basin.

### 3.2.2 Hydrometeorology and physiography of study areas

The surface area of elementary catchments depends upon the location of gauging stations. They have similar ranges for each geology type and study areas (Fig. 5a), and the potential bias induced by scale effect is limited. Figure 5b shows the strong variability of reach mean slope with geology type. This can be explained by the morphology of karst plateaus that are prone to intense erosion, forming canyons with low slopes controlled by the base level. Moreover, higher-elevation ground is mostly underlain by hard-rock (i.e. non-karst) terrains. Contrasting slopes should thus not be seen as a bias but as the result of an intrinsic characteristic of limestone areas, coinciding with karst occurrence. Reach mean slope on karst catchments has similar values for the three

study areas, with a maximum variation in median slope of 5‰ from Normandy to the Cévennes.

Figure 5c and d present the distribution of precipitation and discharge depth for the 20 strongest storm events (see Sect. 2.1) in the 108 elementary catchments, grouped by geology type (K: karst, M: mixed, NK: non-karst) and by study area (C: Cévennes, J: Jura, N: Normandy) for K catchments. A major contrast in precipitation depth is highlighted, with a median value of 100 mm per event for the Cévennes catchments, which is higher than the maximum recorded value for the Normandy catchments (80 mm). Jura Mountains catchments have an intermediate position, with a median rainfall-event depth of around 50 mm. A similar variation is observed for geology types, with median-event rainfall depth increasing from 45 mm on K catchments to 105 mm on NK catchments. This is partly due to the Cévennes upstream catchments of non-karst hard rock receiving intense rainfall. Nevertheless, this has no major influence on the descriptor values since they are normalized by rainfall. The same trend is seen for flow depths, with median values of 20, 10 and 2 mm and 50, 15 and 5 mm for C, J and N and NK, M and K catchments, respectively.

## 4 Results

### 4.1 Water-balance descriptors

The distribution of water-balance descriptors is presented in Fig. 6 (top row) (i) by geology type (K: karst, M: mixed, NK: non-karst) and (ii) by location for karst catchments only (C: Cévennes, J: Jura, N: Normandy). Storm-flow and base-flow depths ($V_S$ and $V_U$) are influenced by geology type, with
<figure><img src="image_8.png" /><figcaption>**Figure 4.** (a) Location of the studied areas with karst regions in France (from BDLISA database). (b, c, d) Distribution of gauging stations and geology type of elementary catchments. Red circles indicate gauging stations in intermediate catchments; green circles indicate gauging stations in headwater catchments.</figcaption></figure>

both decreasing by a factor 10 from NK to K reaches. $V_P$ may have an influence on these results but cannot explain all of this correlation since it only varies by a factor 2. $V_S$ and $V_U$ being calculated from differences between inlet and outlet runoff (see Appendix D) show that karstified catchments generally produce less streamflow than other ones, regarding both quick- and slow-flow components. K and M catchments have their first $V_S$ and $V_U$ quartiles close to zero, highlighting the major streamflow losses along their reaches. $V_S$ and $V_U$ values also vary with location, with $V_S$ marking the typical rainfall intensity of each climatic region (Cévennes, Jura and Normandy in descending order; Fig. 5a), whereas $V_U$ seems to be more influenced by baseflow index information that is lower in Cévennes, probably due to a higher number of NK catchments.

Runoff coefficients are also significantly influenced by geology type, $R_C$ values dropping from 0.42 on NK to 0.14

and 0.19 on K and M, respectively. This trend is observed in a comparable way in the three areas, all showing $R_C$ values between 0.1 and 0.2 for their K reaches. Storm-runoff coefficients show a similar behaviour, with even more contrasting values between K and NK reaches ($R_{C,S}$ being fivefold higher in NK compared to threefold for $R_C$). Storm-flow coefficient $S_C$ corresponds to the part of streamflow due to quick runoff. It decreases from 0.45 to 0.3 in K reaches and also has a geographic pattern, with Cévennes K catchments having 50% of quick-flow components compared to 25% for Jura and Normandy K catchments. This correlates with lower $V_U$ values in the Cévennes K catchments, even with higher rainfall depths.

Water-balance descriptors globally show that karst areas promote more infiltration (lower $R_C$), with the infiltrated water not being released through streamflow at the storm-event timescale (lower $S_C$).
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Elementary catchment area (a); reach mean slope (b) distribution of the 108 gauging stations, grouped by geology type (K: karst, M: mixed, NK: non-karst) and by study area (C: Cévennes, J: Jura, N: Normandy) for K catchments. Precipitation (c) and runoff (d) depth distribution for the 20 selected storm events on the 108 gauging stations, grouped by geology type and study area for K catchments. Values beyond dashed lines are represented on the lines.</figcaption></figure>

## 4.2 Hydrograph-shape descriptors

Distributions of hydrograph-shape descriptors are shown in Fig. 6 (middle row) (i) by geology type and (ii) by study area for karst catchments only. Median values of peak-flow variation normalized by rainfall ($\Delta/V_p$) vary from $4.6 \times 10^{-6}$ s$^{-1}$ for NK catchments to $1.9 \times 10^{-6}$ and $0.5 \times 10^{-6}$ s$^{-1}$ for M and K catchments, respectively. It shows that the larger the karst outcrops, the smaller the peak-flow amplification. This trend was noted in all study areas, even though Cévennes K catchments show more peak-flow amplification than others (and especially a large $\Delta$ variability).

Karst catchments have a median rise duration $T_{Ri}$ of 11 h, whereas M and NK rises take 6 and 4 h, respectively. Stream-flow increase is thus considerably slower in karst reaches. Median durations at 80% of peak flow, $T_{08}$, were also correlated with geology type, with values of 11, 6 and 3 h for K, M and NK catchments, respectively, showing that the slower discharge increase is associated with a buffered peak flow. The median value of the recession time constant ($T_{Re}$) is 44 h for K catchments and 28 h for M and NK ones. Typical recessions in karst reaches are thus 50% slower than in others. The characteristic times $T_{Ri}$, $T_{08}$ and $T_{Re}$ are influenced by karst in different ways depending on study area. For each of them, Cévennes catchments show shorter durations followed by Normandy with durations close to the median values of the whole K-catchments samples and by the Jura.

The flood-wave celerity ($C_G$), obtained by comparing inlet and outlet hydrograph gravity centres, is constant around $1.3 \text{ m s}^{-1}$, regardless of geology type. This means that, in karst catchments, storm events tend to have longer $T_{Ri}$, $T_{08}$ and $T_{Re}$ values, without slowing down of flood-wave routing, corresponding to an increase in the diffusivity of the flood wave at constant celerity. However, celerity shows a regional pattern, consistent with mean regional reach slopes (Fig. 5), with values increasing from $0.8 \text{ m s}^{-1}$ in Normandy to $1.6 \text{ m s}^{-1}$ in the Cévennes, with Jura catchments showing intermediate values.

Figure 7 shows the peak-flow evolution towards catchments with karst. K catchments globally align on the first bisector, showing in most cases low peak-flow amplification. NK catchments mostly lie above the equation line of $y = 1.5x$, showing an increase of at least 50% in specific peak flow, meaning that lateral inflow is higher in those catchments. M catchments have an intermediate hydrologic response.

To summarize, hydrograph-shape descriptors globally show that during the strongest storm events, karst areas tend to decrease peak-flow amplification and increase characteristic flood times, without impacting flood-wave celerity. It should be noted that this general pattern is associated with a large variability of hydrological response in K catchments with locally contrasting behaviour.
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Distribution of hydrological descriptors, grouped by geology type (K: karst, M: mixed, NK: non-karst) and by location for karst catchments (C: Cévennes, J: Jura, N: Normandy). First row: water-balance descriptors; second row: hydrograph-shape descriptors; third row: lateral-exchange descriptors. Red dotted lines show the median value of karst catchments for the whole sample. Values beyond black dashed lines are represented on the lines.</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Figure 7.** Variation of peak flow ($Q_X$, see Appendix D) from inlet to outlet stations, with representation of geology type.</figcaption></figure>

## 4.3 Lateral-exchange descriptors

The distribution of lateral-exchange descriptors is shown in Fig. 6 (bottom row). The component of peak-flow variation due to channel diffusivity ($\Delta_D/V_p$) is quite stable with geology type, with values around $-4 \times 10^{-7} s^{-1}$. The component of peak-flow variation due to channel lateral exchanges ($\Delta_L/V_p$) correlation with geology type is more significant, with median values increasing from $1 \times 10^{-6}$ to $2.5 \times 10^{-6}$ and $5 \times 10^{-6} s^{-1}$ for K, M and NK reaches, respectively. This shows that the strongest lateral inflow in intermediate catchments occurs in NK areas, as expected from the $V_p$ and $V_Q$ variability (Fig. 5), as well as following $Q_O, X$ versus $Q_I, X$ (Fig. 7). It also highlights the compensation of diffusivity peak-flow attenuation by lateral inflow, which is stronger in NK catchments. Lateral exchanges are negative for the first quartile of events on K reaches compared to very few events for NK reaches. Karst influence on $\Delta_D/V_p$ and $\Delta_L/V_p$ is similar in Jura and Normandy catchments, whereas Cévennes catchments show higher peak-flow variations in both components, leading to a strong variability in $\Delta$ values (Fig. 6).

This trend is confirmed by the volume of lateral inflows ($V_{L+}/V_p$) values that decreases with karst occurrence. Volume of lateral outflows ($V_{L-}/V_p$) distribution also varies with geology type, with lateral outflow being reduced in NK catchments. Finally, both descriptors undergo similar karst influence in the three study areas, with only Cévennes K catchments having higher lateral outflow.

Analysis of the simulated lateral hydrographs shows that the weak peak-flow amplification of karst reaches is mostly due to a low exchange component $\Delta_L$, associated with a peak-flow attenuation caused by diffusivity $\Delta_D$. This is confirmed by lateral $V_{L+}/V_p$ inflow that is 3 times lower in K catchments than in NK catchments.

Figure 6 also shows the distribution of $IGF^*/V_p$ ($IGF^* = IGF + \delta$), which indicates that storm events over NK catchments are mostly characterized by incoming IGFs (>75% of

all events). K and M catchments are subject to more balanced incoming and outgoing IGFs, with a median $IGF^*/V_p$ value of 0.06 against 0.2 for NK. Absolute $IGF^*$ depths for K and M catchments are mostly <50 mm per storm event compared to 100 mm or more in NK catchments (results not shown).

$IGF^*/V_p$ distribution highlights some differences for K catchments in each study area. Storm events in Jura K catchments mostly show streamflow excess (>75% of all events), with a mean $IGF^*/V_p$ value of 0.1. Cévennes and Normandy K catchments have a more balanced behaviour, with median $IGF^*/V_p$ values close to zero, and storm events resulting in streamflow excess or deficit in comparable proportions.

## 4.4 Statistical tests

Table 2 shows for each descriptor the results of *t* tests for all combinations of paired samples (Sect. 2.5). The quick- and slow-flow components, $V_S$ and $V_U$, both show statistically significant variations with geology type (increasing from K to NK) and with study area (Fig. 6). Despite visible trends in $R_C, R_{C,S}$ and $S_C$ variation with geology type (Fig. 6), this correlation has a low statistical significance, with only $R_{C,S}$ being discriminant between K and NK catchments. $R_C$ and $R_{C,S}$ on K catchments show significant variation with study area, which was not verified for $S_C$.

Storm-hydrograph shape is strongly related to geology type, with peak-flow amplification reduced and locally attenuated. The longer characteristic times for karst catchments (Fig. 6) have a statistical significance regarding characteristic times. Only flood-wave celerity is not affected by geology type. Karst influence is area specific for $T_{Ri}$ and $T_{08}$ characteristic times (increasing order of characteristic times is C to N to J).

Lateral-exchange descriptors show that peak-flow amplification in karst reaches is limited, mostly due to its exchange component $\Delta_L$ being reduced (Fig. 6). This is confirmed by the volumes of lateral exchange, with $V_{L-}/V_p$ showing that K and M catchments are more prone to lateral outflow than NK catchments and $V_{L+}/V_p$ showing that NK catchments are more prone to lateral inflow than K and M catchments. $IGF^*/V_p$ also seems related to geology type. Table 2 shows that the statistical significance of these relationships is not systematic, M catchments being generally different from NK ones. Only the $IGF^*/V_p$ descriptor does not show a statistical distinction between all samples.
<table><caption>**Table 2.** Synthesis of *t* test decisions for the null hypothesis that data in paired tested samples come from independent random samples from normal distributions with equal means and equal but unknown variances. The result is 1 if the test rejects the null hypothesis at the 5% significance level and 0 otherwise. The test is performed (i) in the three study areas for all combinations of paired samples by geology type (K: karst, M: mixed, NK: non-karst) and (ii) for all combinations of karst-catchment paired samples by study area (C: Cévennes, J: Jura, N: Normandy).</caption><thead><tr><th>Descriptor</th><th>K-NK</th><th>M-NK</th><th>M-K</th><th>C-J</th><th>C-N</th><th>J-N</th></tr></thead><tbody><tr><td colspan="7">Water-balance descriptors</td></tr><tr><td>V<sub>S</sub></td><td>1<sup>a</sup></td><td>1<sup>a</sup></td><td>1<sup>a</sup></td><td>1<sup>c</sup></td><td>1<sup>a</sup></td><td>1<sup>a</sup></td></tr><tr><td>V<sub>U</sub></td><td>1<sup>a</sup></td><td>1<sup>a</sup></td><td>1<sup>a</sup></td><td>1<sup>a</sup></td><td>1<sup>a</sup></td><td>1<sup>a</sup></td></tr><tr><td>R<sub>C</sub></td><td>0</td><td>0</td><td>0</td><td>1<sup>a</sup></td><td>1<sup>b</sup></td><td>1<sup>b</sup></td></tr><tr><td>R<sub>C, S</sub></td><td>1<sup>c</sup></td><td>0</td><td>0</td><td>1<sup>a</sup></td><td>1<sup>b</sup></td><td>1<sup>b</sup></td></tr><tr><td>S<sub>C</sub></td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td colspan="7">Hydrograph-shape descriptors</td></tr><tr><td>Δ/V<sub>P</sub></td><td>0</td><td>1<sup>a</sup></td><td>0</td><td>1<sup>b</sup></td><td>1<sup>c</sup></td><td>0</td></tr><tr><td>T<sub>Ri</sub></td><td>1<sup>a</sup></td><td>1<sup>b</sup></td><td>1<sup>a</sup></td><td>1<sup>a</sup></td><td>1<sup>a</sup></td><td>1<sup>a</sup></td></tr><tr><td>T<sub>08</sub></td><td>1<sup>a</sup></td><td>1<sup>a</sup></td><td>1<sup>a</sup></td><td>1<sup>a</sup></td><td>1<sup>a</sup></td><td>1<sup>a</sup></td></tr><tr><td>T<sub>Re</sub></td><td>1<sup>b</sup></td><td>0</td><td>1<sup>c</sup></td><td>0</td><td>0</td><td>0</td></tr><tr><td>C<sub>G</sub></td><td>0</td><td>0</td><td>0</td><td>0</td><td>1<sup>c</sup></td><td>0</td></tr><tr><td colspan="7">Lateral-exchange descriptors</td></tr><tr><td>Δ<sub>D</sub>/V<sub>P</sub></td><td>0</td><td>1<sup>c</sup></td><td>1<sup>c</sup></td><td>1<sup>a</sup></td><td>1<sup>b</sup></td><td>0</td></tr><tr><td>Δ<sub>L</sub>/V<sub>P</sub></td><td>0</td><td>1<sup>b</sup></td><td>0</td><td>1<sup>a</sup></td><td>1<sup>c</sup></td><td>0</td></tr><tr><td>V<sub>L+</sub>/V<sub>P</sub></td><td>0</td><td>1<sup>c</sup></td><td>0</td><td>0</td><td>1<sup>b</sup></td><td>1<sup>c</sup></td></tr><tr><td>V<sub>L-</sub>/V<sub>P</sub></td><td>0</td><td>1<sup>b</sup></td><td>0</td><td>1<sup>b</sup></td><td>0</td><td>0</td></tr><tr><td>IGF*/V<sub>P</sub></td><td>0</td><td>0</td><td>0</td><td>1<sup>a</sup></td><td>0</td><td>0</td></tr></tbody></table>

p values are indicated with letters:

a = *p* < 0.001, b = 0.001 < *p* < 0.01 and c = 0.01 < *p* < 0.05.

# 5 Discussion

## 5.1 Factors driving flood processes

The analysis of descriptor distribution (Fig. 6) and its statistical significance (Table 2) shows that several factors affect flood processes at the elementary catchment scale. Water-balance descriptors show that quick- and slow-flow depths, $V_S$ and $V_U$, vary with geology type in a statistically significant way. This is partly explained by the role of karst in how catchments transform precipitation to runoff ($R_C$ and $R_{C, S}$) and on the relative proportions of quick and slow flow ($S_C$). Nevertheless, only the $R_{C, S}$ test is statistically significant according to geology type. This highlights the role of climatic influence, as $V_S$ and $V_U$ also depend on rainfall depth $V_P$ that varies with study area, but their range is clearly higher than that of $V_P$.

In addition, hydrograph descriptors show a reduction in peak-flow amplification on K catchments in all study areas, associated with an increase of characteristic times. This trend is consistent with previous work showing peak-flow attenuation

in karst areas (De Waele et al., 2010; Charlier et al., 2019). The influence on characteristic times is area specific, with Jura catchments having greater inertia than Cévennes ones. This area-dependent nature of characteristic flood times may be partly explained by different karst settings and occurrence, but it is also linked with rainfall patterns. Cévennes catchments have a Mediterranean climate, with typical intense and short storm events triggering flash floods (Marty et al., 2013). The inertial events of Jura catchments can be explained by their great length compared to other areas: Cévennes and Normandy rivers hardly reach 150 km before reaching the sea or the Rhône, but the downstream Jura station is 400 km from the source along the Doubs river.

Finally, lateral-exchange descriptors show contrasting exchanges between different geology types, with more lateral inflow in NK catchments and more lateral outflow in K catchments. This agrees with previous studies on storm events in karst catchments highlighting river losses in karst reaches (Delrieu et al., 2005; Perrin and Tournoud, 2009; Bailly-Comte et al., 2012; Charlier et al., 2019). This trend companies a study-area variation, where the main aspect is a higher variability of exchanges for Cévennes catchments and a lower one for Normandy. This is probably related to the need of harmonizing spatial scales between gain-loss processes and gauging networks. In fact, earlier work (Toth, 1963; Schaller and Fan, 2009; Bouaziz et al., 2018; Fan, 2019) shows that the size of the investigated catchment affects the IGF importance, with greater areas being more likely to be self-contained. In the case of our study areas, as the median gauged areas are similar for each of them (Table 1), it might be explained by the generally thick soil and epikarst in Normandy, which is very reduced in the mountainous and Causses areas of Cévennes. Indeed, soil and epikarst are more likely to promote subsurface flow with closer zones of gains and losses, whereas exposed karst drains, as in the Cévennes, enlarge the spatial scale of IGF processes by connecting river losses to trans-catchment karst aquifers.

Hydrological processes have been shown to be influenced by physiographic parameters such as karst occurrence. Nevertheless, other drivers can control catchments hydrological response to storm events. As an example, Fig. 8 shows IGF* depths distribution for the three studied sites, according to seasons. Seasons have been selected in order to reflect the main periods of the hydrological year, and to be suitable for the three sites. Throughout the hydrological year (from September to August), median IGF* depth is continuously decreasing for all sites towards zero. This is probably linked to the hydrological conditions of catchments, in particular the saturation state of aquifers. Indeed, low-water-table periods (April to August) are more likely to limit IGFs. In the case of Jura catchments, a majority of events show IGF gains during high-flow periods (fall and winter), whereas a majority of events promote IGF losses during low-flow periods. This reversal of IGF direction, in the particular case of karst catch-
<figure><img src="image_12.png" /><figcaption>**Figure 8.** IGF* depth distributions for the three studied sites, according to seasons.</figcaption></figure>

ments, can also be favoured by the complex organization of
the underground conduit networks. In fact, connection or dis-
connection of the main networks could temporally modify
the hydrogeological catchment boundaries with threshold ef-
fects (Charlier et al., 2012; Bonacci, 2015). For those rea-
sons, spatial and temporal variability of karst catchments be-
haviour are still challenging to characterize and predict. A
comprehensive understanding of flood processes would thus
imply accurate and continuous monitoring of climatic, hy-
drological and hydrogeological variables on multiple catch-
ments covering various physiographical settings.

## 5.2 Regional patterns and typology

Figure 6 shows that storm events in K catchments in the Cévennes are mostly associated with losing IGF\*s (that can include storage variation), whereas Jura events are mostly associated with gaining IGF\* values. This can be linked to the typology of karst settings. The Cévennes region is characterized by binary karst systems, with large upstream NK (hard-rock) terrains and downstream limestone plateaus where karst influence occurs. In such a geological setting, karst areas mostly play a role in flood attenuation as they lie downstream of reliefs with intense rainfall events and high runoff coefficients (see the example of the Tarn river in Charlier et al., 2015). The Jura Mountains region is regionally much more homogeneous, with widespread karst formation affecting the limestone plateaus, and few areas covered with Quaternary deposits. In this setting, karst can alternately promote streamflow capture (attenuation) and generation (amplification), depending on the location of river losses and the interaction between surface water and groundwater (Le Mesnil et al., 2020).

To discuss the spatial variability of IGFs, the median
storm-event IGF* value [mm] is shown for each elementary

catchment of the three studied areas in Fig. 9d to f, along
with the geological maps (Fig. 9a to c). In the Cévennes
(Fig. 9, left row), downstream karstified parts can be divided
into two zones defined by a different lithology. The eastern
part is mostly underlain by Cretaceous limestone and Ceno-
zoic formations, whereas the western part is mostly com-
posed of (older) Jurassic limestone. Despite IGF* account-
ing for the potential storage variation (IGF* = IGF + δ), all
eastern Cretaceous catchments have negative median IGF*
values, meaning that almost certainly groundwater flows out
through karst aquifers during storm events, without being re-
covered in the investigated area during the considered event.
The western Jurassic catchments have a less significant neg-
ative IGF*, with most values comprised between -1.45 and
3.5 mm. Through the lithology, this highlights the local role
of karst occurrence, which is superimposed on the regional
role of general karst settings.

In the Jura Mountains along the main Doubs river (Fig. 9,
second row), almost all elementary catchments have neg-
ative IGF* during storm events, contrary to the tributary
catchments that mostly have positive IGF*. This is coherent
with the well-known Doubs river capture, at least for
the south-eastern catchments. The tributaries being upstream
and higher above the base level, river loss into karst aquifers
is less important or (partly) recovered within a short dis-
tance. Storm-event IGF mapping highlights the already well-
documented zone of Doubs losses feeding the Loue catch-
ment (e.g. Charlier et al., 2014).

Normandy catchments (Fig. 9, third row) show less im-
portant IGF* depths, most catchments having values close to
zero. This is the result of a too widespread gauging network
compared to the spatial scale of IGF processes here. River-
capture phenomena are well known in some of the catch-
ments (e.g. Charlier et al., 2019), but the gaining and los-
ing zones mostly seem to fall in the same catchment due to
gauging station locations (Sect. 5.1).

Figure 9g to i show the percentage of area variation of the elementary catchments, from the topographic catchment area to the theoretical area of the hydrogeologically active catchment (i.e. without excess or deficit in water balance). The area of such catchments is calculated from the IGF* and Vp depths. Assuming a spatially homogeneous rainfall intensity, the area variation corresponds to the surface to be added or withdrawn, in regards to the topographic elementary catchment, to obtain the median IGF* value. These variations reflect the incoming or outgoing IGF processes, respectively. For Cévennes and Jura catchments, the values of area variation are of similar magnitudes, ranging from approximately -7% to +90%. Regarding Normandy catchments, area variations are slightly lower, with values ranging from -20% to +32%. This shows that IGF can be an important term of the water balance during storm events and has to be considered as such in order to better understand flow processes during floods.
<figure><img src="image_13.png" /><figcaption>**Figure 9.** Main geological features in the three studied areas (a, b, c), median storm-event IGF* values (including potential aquifer storage variation) (d, e, f), and area variation from actual topographic to theoretical hydrogeologically active catchment [in %] (g, h, i) and in square kilometres [km²] (j, k, l).</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>**Figure 10.** Relationships between multi-annual water-balance-derived IGF<sub>A</sub> and storm-event-derived IGF* (including potential aquifer storage variation) depths. IGF<sub>A</sub> depth is normalized by event duration for consistent comparisons. First, all 108 elementary catchments were plotted (a); then catchments were filtered to retain only α values between 0.5 and 1; these were plotted and grouped by geology type (b) and study area (c).</figcaption></figure>

Figure 9j to l show the same elementary catchment area variations expressed in square kilometres [km²]. These values are correlated to the size of topographic catchments. They allow for a concrete representation of the recharge areas located outside of the surface catchments. As this study is based on storm events only, and considering that the studied catchments extend further downstream, the sum of theoretical active catchments (which is a first approximation of the hydrogeological catchment) do not match with the topographic catchment area. Nevertheless, a previous study by Le Mesnil et al. (2020) showed that, along the Doubs river, annual IGF takes part in the water balance that progressively decreases from spring to outlet and tends towards zero.

## 5.3 IGF: annual vs. event scale and impacts on quick- and slow-flow components

In this section, we compare a major missing term of the water budget, alternately estimated by two approaches. First, at the annual scale, the multi-annual IGF (noted IGF<sub>A</sub>) is calculated by Le Mesnil et al. (2020), under the assumption of nil annual stock variation in the aquifer (δ = 0), which is often verified when using time series of several years. Second, at the event scale, IGF* = IGF + δ is calculated here as defined in Sect. 2.3.3.

Le Mesnil et al. (2020) assessed multi-annual values of IGF<sub>A</sub> depth for the catchment of the present work. This IGF<sub>A</sub> value, calculated by adapting the two-stage precipitation partitioning theory of Lvovich (1979), is associated with a parameter “α”, estimating the relative impact of IGF on the rapid- and slow-flow components at an annual scale. The α values between 0 and 0.5 correspond to an annual IGF<sub>A</sub> mostly affecting the slow-flow component, whereas α values between 0.5 and 1 correspond to an annual IGF<sub>A</sub> mostly affecting the rapid-flow component; α values below 0 and over 1 indicate compensating IGF<sub>A</sub> flow, such as sur-

face loss combined with groundwater gain. To discuss the link between IGF estimates at both annual and storm-event timescales, Fig. 10 presents the storm-event IGF* depths calculated for the present work as a function of the annual IGF<sub>A</sub> depth obtained in the precedent paper for the 108 elementary catchments. Annual IGF<sub>A</sub> values were normalized by event duration for consistent comparisons.

Figure 10a shows some points with opposed annual IGF<sub>A</sub> and event IGF* signs (e.g. a losing annual water-balance-derived IGF<sub>A</sub> associated with a gaining event-derived IGF*). The point cloud does not show significant axial organization, leading to a low determination coefficient value for the linear regression. Figure 10b and c present the same relationships for catchments with α values between 0.5 and 1, i.e. annual IGF<sub>A</sub> mostly affecting the quick-flow component (Le Mesnil et al., 2020). The R<sup>2</sup> value of the linear regression has increased but is still low (0.222, not shown). Nevertheless, most catchments show consistent annual storm-event IGF* signs, falling into the upper-right or lower-left quadrants. This means that the annual estimation is good for assessing the relative impact of IGF on quick- and slow-flow components. Moreover, IGF* maps (Fig. 9) are in accordance with IGF<sub>A</sub> results by Le Mesnil et al. (2020). For example, the well-known phenomenon of Doubs river losses feeding the Loue catchment is visible at the event scale.

R<sup>2</sup> values in Fig. 10b and c show that correlations are more reliable when operating regression on groups (by geology type and study area) of K and M catchments. The two groups with a reliable relationship between annual and event-derived IGF depths are the karst catchments K and the Jura catchments J, with respective R<sup>2</sup> values of 0.664 and 0.663. The stronger relationship in K catchments can be explained by the occurrence of IGF during both recession and flood periods, whereas IGF occurred mostly during flood periods in NK catchments. This agrees with the high positive IGF* values (>50 mm per event) associated with low annual IGF<sub>A</sub>
<figure><img src="image_15.png" /><figcaption>**Figure 11.** *V*<sub>S</sub> depths versus IGF*<sup>*</sup> depths for all events and reaches, with differentiation of geology types</figcaption></figure>

values for NK. This may also explain the higher R² value of Jura catchments, as they are mostly K catchments, while the Cévennes also have NK catchments resulting in a lower R² value.

Figure 11 represents V_S versus IGF* (both expressed in depth) for all events and reaches, with different symbols according to geology type. With V_S being equal to V_S_o - V_S_i, the V_S versus IGF* relationship provides an insight into the relative proportions of total IGF* depth and quick-flow component variation. First, while IGF* values vary from -200 to 300 mm, V_S values are mostly positive, ranging from -20 to 150 mm. This indicates that for most events, when outgoing IGFs (negative IGF* values) occur, the quick-flow component still increases. This is particularly true on NK catchments, whereas some events with negative IGF* on K catchments show negative (or low positive) V_S values. In that case, a likely hypothesis is the presence of losses in karst riverbeds via sinkholes, which are not present in non-karst reaches (see such example in Charlier et al., 2019). Regarding events with positive IGF*, V_S depths are mainly inferior to IGF* ones, with values around half of the total IGF* ones. This shows that incoming IGFs during storm events feed the slow-flow component of the total streamflow in similar proportions than the quick-flow component. This is linked to the different IGF processes occurring at the catchment scale: for example, slow-flow component corresponding to aquifer drainage by the river (e.g. Bailly-Comte et al., 2012) and quick-flow component corresponding to karst spring activation (Bonacci and Bojanic, 1991; Maréchal et al., 2008).

# 6 Conclusions

We carried out a spatialized analysis of 15 easily calculable descriptors characterizing water balance, hydrograph shape and lateral exchanges for a set of 20 storm-event data at the elementary catchment scale for each of the 108 gaug-

ing stations, controlling karst and non-karst regions. The results show that karst promotes higher water infiltration, with this water being mostly retained during storm events. Karst increases characteristic flood times and limits peak-flow amplification, without affecting flood-wave celerity much. This is interpreted to be due to an interbasin groundwater flow (IGF) loss that can be high at the storm-event scale, representing around 50% of the discharge at a catchment outlet and 20% of rainfall. A spatial variability of those effects is linked to differences in karst regions: binary karst catchments mostly attenuate floods, whereas extended karst plateaus undergo alternated losses and gains. Secondary factors include climatic influence (regional variability of rainfall-event intensity) and the spatial-scale match between gain-loss processes and spacing of the gauging network. A seasonal effect has also to be considered regarding IGF magnitude and direction.

The existence of karst hydrological specificities has been known for decades but has not been quantified to a large extent, especially regarding its impacts on flood processes at the scale of the river reach. Although some research had been done on this topic, it leads to hindered modelling performance in many cases. We have quantified several important parameters for a large set of catchments, for the first time, in a spatialized study based on event-scale processes, contributing to build a common understanding at regional scale of karst behaviour during storm events and thus improving modelling and forecasting capabilities in such terrains. Although our approach is based on karst areas, it stays generic, and we hope future work will investigate other relationships between the hydrological response and physiographic characteristics of catchments, such as soil types, land use, climate, etc.
## Appendix A: Hydrograph decomposition

At each gauging station, discharge values were filtered in order to separate the quick (storm flow, $Q_S$) and slow (underground flow, $Q_U$) flow components. The quick is traditionally interpreted as the surface component, and the slow one as baseflow, which is that part of streamflow corresponding to aquifer drainage. However, in the case of karst catchments, aquifer drainage can also produce a quick-flow signal because of short transfer times through conduits. Several baseflow-separation methods exist, with most being based on graphical analysis, like the fixed-interval, sliding-interval, local-minimum, or Wallingford methods (Gustard et al., 1992; Sloto and Crouse, 1996; Rutledge, 1998; Piggott et al., 2005). Numerical approaches have also been developed (e.g. Lyne and Hollick, 1979; Eckhardt, 2005). In this study, we used an automation of the one-parameter recursive digital filter proposed by Lyne and Hollick (1979), implemented in the HydRun package (Tang and Carey, 2017). The filter equation is defined as

$$Q_S(t) = \beta Q_S(t-1) + \frac{1+\beta}{2} [Q(t) - Q(t-1)], \quad (A1)$$

with $Q_S(t)$ and $Q(t)$ being the filtered quick-flow component and total streamflow at time $t$, respectively, and $\beta$ being the filter parameter.

We chose this method as it provides consistent results, like those obtained with graphical approaches (results not shown). It can easily be automated and has only one $\beta$ parameter, fixed at 0.91 after a trial-and-error analysis on the studied catchments and considering the results of Nathan and McMahon (1990) for 186 catchments.
## Appendix B: Lateral flow simulation using the diffusive wave equation

### Diffusive wave equation

An inverse modelling approach is adopted for simulating lateral flow between two gauging stations. This approach simulates the lateral flow $Q_L$, based on measurements from two gauging stations $Q_I$ and $Q_O$.

The diffusive wave equation (DWE), accounting for lateral flow, is an approximation of the St-Venant equation that can be written as

$$ \frac{\partial Q}{\partial t} + C(Q) \left[ \frac{\partial Q}{\partial x} - q \right] - D(Q) \left[ \frac{\partial^2 Q}{\partial x^2} - \frac{\partial q}{\partial x} \right] = 0, \quad (B1) $$

where $x$ [L] is the length along the channel, $t$ [T] is the time, and celerity $C(Q)$ [LT⁻¹] and diffusivity $D(Q)$ [L²T⁻¹] are functions of the discharge $Q$ [L³T⁻¹]. The term $q(x, t)$ [L²T⁻¹] represents the lateral flow distribution. The lateral hydrograph $Q_L(t)$ is given by

$$ Q_L(t) = \int_0^t q(x, t) dx, \quad (B2) $$

with $l$ [L] being the channel length.

Moussa (1996) extended the solution of the DWE under Hayami's hypotheses (semi-infinite channel, $C(Q)$ and $D(Q)$ constant) to the case where lateral flow is uniformly distributed along the channel. Let $I(t)$ and $O(t)$ be the inlet flow minus baseflow and the outlet flow minus baseflow, respectively:

$$ O(t) = \varphi(t) + [I(t) - \varphi(t)] \cdot K(t), \quad (B3) $$

with $K(t)$ being the Hayami Kernel function defined as

$$ K(t) = \frac{l}{2(\pi D)^{1/2}} \frac{e^{\left[ \frac{Ct}{4D} \left( 2 - \frac{l}{Ct} - \frac{Ct}{T} \right) \right]}}{t^{3/2}} \quad (B4) $$

and

$$ \varphi(t) = \frac{C}{l} \int_{0}^{t} [Q_{L}(\theta) - Q_{L}(0)] d\theta. \quad (B5) $$

### The inverse problem

Under Hayami's conditions and assuming that lateral flow is uniformly distributed along the channel, Moussa (1996) proposed a solution of the inverse problem; this enables evaluation of the temporal distribution of lateral flow $Q_L(t)$ over the channel reach by knowing $I(t)$ and $O(t)$. Knowing $C$, $D$ and $l$, the lateral flow can be calculated using the following

procedure:

$$ L(t) = O(t) - I(t) \cdot K(t), \quad (B6) $$

$$ K^i(t) = K \cdot K \cdots K \text{ (i times),} \quad (B7) $$

$$ \varphi(t) = L(t) + L(t) \sum_{i=1}^{\infty} K^i(t). \quad (B8) $$

and finally the lateral flow $Q_L, C(t)$ is

$$ Q_{L, C}(t) = Q_{L}(0) + \frac{l}{C} \frac{d\varphi}{dt}. \quad (B9) $$

## Appendix C: Estimation of effective rainfall

Effective rainfall $P_{eff}$ is estimated using three different approaches, in order to provide a range for this component characterized by major uncertainties. The three approaches are based on the water-budget methods proposed by Thornthwaite (1948) and Dingman (2015), and on the GR lumped model (Edijatno et al., 1999). All three consider soil as a reservoir, used for separating the input (precipitation) into evapotranspiration and effective rainfall. The capacity of the soil reservoir $C_{max}$ is estimated with the map of available soil-water capacity from INRAE (Le Bas, 2018).

In the Thornthwaite method, water in the soil reservoir is directly available for evapotranspiration, and precipitation produces effective rainfall ($P_{eff}$) only after soil saturation. The following algorithm summarizes the method.

- If $P < E0$, the difference $E0 - P$ is subtracted from the soil-water stock, $B$, until it is empty:

- $B_t = \max(0; B_{t-1} + P_t - E0_t)$,

- $E_t = \min(E0_t; B_{t-1} + P_t)$,

- $P_{eff, t} = 0$.

- If $P > E0$, the difference $P - E0$ first feeds the soil-water stock, $B$, and then produces efficient rainfall:

- $B_t = \min(B_{max}; B_{t-1} + P_t - E0_t)$,

- $E_t = E0_t$,

- $P_{eff, t} = \max(0; B_t + P_t - E0_t - B_{max})$.

The Dingman method is similar to the previous one, with an exponential law governing water extraction for evapotranspiration from the soil reservoir.

- If $P < E0$, the difference $E0 - P$ is subtracted from the soil-water stock, $B$, following an exponential law:

- $B_t = B_{t-1} \cdot e^{\frac{-(E0_t - P_t)}{B_{max}}}$,

- $E_t = P_t + B_{t-1} - B_t$,

- $P_{eff, t} = 0$.
- If $P > E_0$, the difference $P - E_0$ first feeds the soil-water stock, $B$, and then produces efficient rainfall (as in the Thornthwaite method):

- $B_t = \min(B_{max}; B_{t-1} + P_t - E_0)$,

- $E_t = E_0$,

- $P_{eff,t} = \max(0; B_t + P_t - E_0 - B_{max})$.

The GR method is derived from the GR hydrological models (Edijatno et al., 1999) and involves a quadratic law for the water-level variation in the soil reservoir. The algorithm, summarized below, was then adapted to the BRGM "Gardenia" model (Thiéry, 2014), which has been used here.

- If $P < E_0$, the difference $E_n = E_0 - P$ is subtracted from the soil-water stock, $B$, following a quadratic law:

- $dB = ((B/B_{max})^2 - 2(B/B_{max})) \cdot dE_n$,

- $dE_t = -dB$,

- $P_{eff} = 0$.

- If $P > E_0$, the difference $P_n = P - E_0$ is partitioned into effective rainfall and soil storage following a quadratic law:

- $dB = (1 - (B/B_{max})^2) \cdot dP_n$,

- $E = E_0$,

- $dP_{eff} = (B/B_{max})^2 \cdot dP_n$.

- Integration of the differential variations provides expressions of $B_t$, $E_t$ and $P_{eff,t}$ as a function of $B_{t-1}$, $B_{max}$, and $\tanh(E_n/B_{max})$ or $\tanh(P_n/B_{max})$.

The final $P_{eff}$ value corresponds to the mean of the three estimation method results.
Appendix D

<table><caption>Table D1. List of symbols.</caption><thead><tr><th>Symbol</th><th>Unit</th><th>Formula</th><th>Description</th></tr></thead><tbody><tr><td colspan="4">Water-balance descriptors</td></tr><tr><td>V<sub>S</sub></td><td>mm</td><td>(V<sub>O,S</sub> - V<sub>I,S</sub>)/A</td><td>Storm-flow depth</td></tr><tr><td>V<sub>U</sub></td><td>mm</td><td>(V<sub>O,U</sub> - V<sub>I,U</sub>)/A</td><td>Baseflow depth</td></tr><tr><td>R<sub>C</sub></td><td>-</td><td>(V<sub>O</sub> - V<sub>I</sub>)/V<sub>P</sub></td><td>Runoff coefficient</td></tr><tr><td>R<sub>C,S</sub></td><td>-</td><td>V<sub>S</sub>/V<sub>P</sub></td><td>Storm-runoff coefficient</td></tr><tr><td>S<sub>C</sub></td><td>-</td><td>V<sub>S</sub>/V<sub>Q</sub></td><td>Storm-flow coefficient</td></tr><tr><td colspan="4">Hydrograph-shape descriptors</td></tr><tr><td>Δ/V<sub>P</sub></td><td>s<sup>-1</sup></td><td>(Q<sub>O,X</sub> - Q<sub>I,X</sub>)/V<sub>P</sub></td><td>Peak-flow variation normalized by rainfall</td></tr><tr><td>T<sub>Ri</sub></td><td>h</td><td>-</td><td>Time duration of rise from 0.5Q<sub>O,X</sub> to Q<sub>O,X</sub></td></tr><tr><td>T<sub>08</sub></td><td>h</td><td>-</td><td>Time duration at Q<sub>O</sub> > 0.8Q<sub>O,X</sub></td></tr><tr><td>T<sub>Re</sub></td><td>h</td><td>Eq. (7)</td><td>Time constant of recession (exponential model)</td></tr><tr><td>C<sub>G</sub></td><td>m s<sup>-1</sup></td><td>l/(t(G<sub>O</sub>) - t(G<sub>I</sub>))</td><td>Celerity (based on hydrograph gravity centres)</td></tr><tr><td colspan="4">Lateral-exchange descriptors</td></tr><tr><td>Δ<sub>D</sub>/V<sub>P</sub></td><td>s<sup>-1</sup></td><td>(Q<sub>I,R,X</sub> - Q<sub>I,X</sub>)/V<sub>P</sub></td><td>Diffusivity peak-flow variation normalized by rainfall</td></tr><tr><td>Δ<sub>L</sub>/V<sub>P</sub></td><td>s<sup>-1</sup></td><td>(Q<sub>O,X</sub> - Q<sub>I,R,X</sub>)/V<sub>P</sub></td><td>Exchange peak-flow variation normalized by rainfall</td></tr><tr><td>V<sub>L+</sub>/V<sub>P</sub></td><td>-</td><td>V<sub>L+</sub>/V<sub>P</sub></td><td>Lateral inflow normalized by rainfall</td></tr><tr><td>V<sub>L-</sub>/V<sub>P</sub></td><td>-</td><td>V<sub>L-</sub>/V<sub>P</sub></td><td>Lateral outflow normalized by rainfall</td></tr><tr><td>IGF<sup>*</sup>/V<sub>P</sub></td><td>mm</td><td>(V<sub>Q</sub> - V<sub>P<sub>eff</sub></sub>)/V<sub>P</sub></td><td>Interbasin groundwater flow normalized by rainfall</td></tr><tr><td colspan="4">Other variables</td></tr><tr><td>A</td><td>m<sup>2</sup></td><td>-</td><td>Catchment area</td></tr><tr><td>β</td><td>h<sup>-1</sup></td><td>-</td><td>Recession coefficient as per Maillet approximation</td></tr><tr><td>C</td><td>m s<sup>-1</sup></td><td>-</td><td>Celerity, parameter of the DWE here taken as C<sub>G</sub></td></tr><tr><td>D</td><td>-</td><td>-</td><td>Diffusivity, parameter of the DWE</td></tr><tr><td>δ</td><td>mm</td><td>-</td><td>Aquifer-storage variation, expressed as catchment depth</td></tr><tr><td>E</td><td>mm</td><td>Appendix C</td><td>Actual evapotranspiration</td></tr><tr><td>G<sub>I</sub>; G<sub>O</sub></td><td>-</td><td>-</td><td>Inlet and outlet hydrograph gravity centres</td></tr><tr><td>IGF<sup>*</sup></td><td>mm</td><td>V<sub>O</sub> - V<sub>I</sub> - V<sub>P<sub>eff</sub></sub></td><td>IGF and aquifer variation equal to streamflow excess or deficit</td></tr><tr><td>IGF<sub>A</sub></td><td>mm</td><td>-</td><td>Multi-annual IGF calculated in Le Mesnil et al. (2020)</td></tr><tr><td>l</td><td>m</td><td>-</td><td>River reach length</td></tr><tr><td>P; V<sub>P</sub></td><td>mm s<sup>-1</sup>; mm</td><td>-</td><td>Precipitation rate and depth on the elementary catchment</td></tr><tr><td>P<sub>eff</sub>; V<sub>P<sub>eff</sub></sub></td><td>mm s<sup>-1</sup>; mm</td><td>Appendix C</td><td>Effective rainfall rate and depth</td></tr><tr><td>q</td><td>m<sup>3</sup> s<sup>-1</sup></td><td>-</td><td>Streamflow at the beginning of recession limb</td></tr><tr><td>Q<sub>O</sub>; V<sub>O</sub></td><td>m<sup>3</sup> s<sup>-1</sup>; mm</td><td>-; -</td><td>Streamflow and water depth at catchment outlet</td></tr><tr><td>Q<sub>I</sub>; V<sub>I</sub></td><td>m<sup>3</sup> s<sup>-1</sup>; mm</td><td>-; -</td><td>Streamflow and water depth at catchment inlet</td></tr><tr><td>Q; V<sub>Q</sub></td><td>m<sup>3</sup> s<sup>-1</sup>; mm</td><td>Q<sub>O</sub> - Q<sub>I</sub>; V<sub>O</sub> - V<sub>I</sub></td><td>Streamflow and water depth at elementary catchment</td></tr><tr><td>Q<sub>O,U</sub>; V<sub>O,U</sub></td><td>m<sup>3</sup> s<sup>-1</sup>; mm</td><td>-; -</td><td>Slow streamflow and water depth at catchment outlet</td></tr><tr><td>Q<sub>O,S</sub>; V<sub>O,S</sub></td><td>m<sup>3</sup> s<sup>-1</sup>; mm</td><td>-; -</td><td>Quick-flow component and water depth at catchment outlet</td></tr><tr><td>Q<sub>I,U</sub>; V<sub>I,U</sub></td><td>m<sup>3</sup> s<sup>-1</sup>; mm</td><td>-; -</td><td>Slow-flow component and water depth at catchment inlet</td></tr><tr><td>Q<sub>I,S</sub>; V<sub>I,S</sub></td><td>m<sup>3</sup> s<sup>-1</sup>; mm</td><td>-; -</td><td>Quick-flow component and water depth at catchment inlet</td></tr><tr><td>Q<sub>U</sub></td><td>m<sup>3</sup> s<sup>-1</sup></td><td>Q<sub>O,U</sub> - Q<sub>I,U</sub></td><td>Slow-flow component at elementary catchment</td></tr><tr><td>Q<sub>S</sub></td><td>m<sup>3</sup> s<sup>-1</sup></td><td>Q<sub>O,S</sub> - Q<sub>I,S</sub></td><td>Quick-flow component at elementary catchment</td></tr><tr><td>Q<sub>O,X</sub></td><td>m<sup>3</sup> s<sup>-1</sup> km<sup>-2</sup></td><td>-</td><td>Peak flow at catchment outlet normalized by A</td></tr><tr><td>Q<sub>I,X</sub></td><td>m<sup>3</sup> s<sup>-1</sup> km<sup>-2</sup></td><td>-</td><td>Peak flow at catchment inlet normalized by A</td></tr><tr><td>Q<sub>L</sub></td><td>m<sup>3</sup> s<sup>-1</sup></td><td>Appendix B</td><td>Simulated lateral-exchange flow</td></tr><tr><td>Q<sub>I,R</sub></td><td>m<sup>3</sup> s<sup>-1</sup></td><td>Appendix B</td><td>Routed inlet streamflow</td></tr><tr><td>Q<sub>I,R,X</sub></td><td>m<sup>3</sup> s<sup>-1</sup></td><td>Appendix B</td><td>Routed inlet streamflow peak flow</td></tr><tr><td>Q<sub>Maillet</sub></td><td>m<sup>3</sup> s<sup>-1</sup></td><td>q · exp(-βt)</td><td>Maillet approximation of recession streamflow</td></tr><tr><td>T<sub>G</sub></td><td>h</td><td>-</td><td>Elapsed time between G<sub>I</sub> and G<sub>O</sub></td></tr></tbody></table>
*Data availability.* Hourly streamflow data were gathered from the “Banque hydro” database using a script available at https://doi.org/10.5281/zenodo.3744183 (Hakoun and Manlay, 2020).

*Author contributions.* RM and JBC were involved in conceptualization, funding acquisition and supervision. MLM gathered the data and designed the methodology with the help of JBC, RM and YC. MLM prepared the article with the help of all co-authors.

*Competing interests.* The authors declare that they have no conflict of interest.

*Acknowledgements.* The authors thank the Editor Jim Freer and three anonymous reviewers for their constructive comments.

*Financial support.* The work was funded by the French Governmental Administration for Risk Prevention (DGPR), the Service Central d'Hydrométéorologie et d'Appui à la Prévision des Inondations (SCHAPI) and the French Geological Survey (BRGM).

*Review statement.* This paper was edited by Jim Freer and reviewed by three anonymous referees.

## References

Bailly-Comte, V., Jourde, H., and Pistre, S.: Conceptualization and classification of groundwater-surface water hydrodynamic interactions in karst watersheds: Case of the karst watershed of the Coulazou River (Southern France), J. Hydrol., 376, 456–462, https://doi.org/10.1016/j.jhydrol.2009.07.053, 2009.

Bailly-Comte, V., Borrell-Estupina, V., Jourde, H., and Pistre, S.: A conceptual semidistributed model of the Coulazou River as a tool for assessing surface water-karst groundwater interactions during flood in Mediterranean ephemeral rivers: karst contribution to surface flow, Water Resour. Res., 48, W09534, 2012.

Bakalowicz, M.: Karst groundwater: a challenge for new resources, Hydrogeol. J., 13, 148–160, https://doi.org/10.1007/s10040-004-0402-9, 2005.

Basha, H. A.: Flow Recession Equations for Karst Systems, Water Resour. Res., 56, e2020WR027384, https://doi.org/10.1029/2020WR027384, 2020.

Bell, F. C. and Om Kar, S.: Characteristic response times in design flood estimation, J. Hydrol., 8, 173–196, https://doi.org/10.1016/0022-1694(69)90120-6, 1969.

Bonacci, O.: Karst hydrogeology/hydrology of dinaric chain and isles, Environ. Earth Sci., 74, 37–55, https://doi.org/10.1007/s12665-014-3677-8, 2015.

Bonacci, O. and Bojanić, D.: Rhythmic karst springs, Hydrol. Sci., J., 36, 35–47, https://doi.org/10.1080/02626669109492483, 1991.

Bonacci, O., Ljubenkov, I., and Roje-Bonacci, T.: Karst flash floods: an example from the Dinaric karst (Croatia), Nat. Hazards Earth Syst. Sci., 6, 195–203, https://doi.org/10.5194/nhess-6-195-2006, 2006.

Bouaziz, L., Weerts, A., Schellekens, J., Sprokkereef, E., Stam, J., Savenije, H., and Hrachowitz, M.: Redressing the balance: quantifying net intercatchment groundwater flows, Hydrol. Earth Syst. Sci., 22, 6415–6434, https://doi.org/10.5194/hess-22-6415-2018, 2018.

Champeaux, J. L., Laurantin, O., Mercier, B., et al.: Quantitative precipitation estimations using rain gauges and radar networks: inventory and prospects at Meteo-France, in: WMO Joint Meeting of CGS Expert Team on Surface-based Remotely-sensed Observations and CIMO Expert Team on Operational Remote Sensing, 2011.

Charlier, J.-B., Bertrand, C., and Mudry, J.: Conceptual hydrological model of flow and transport of dissolved organic carbon in a small Jura karst system, J. Hydrol., 460/461, 52–64, https://doi.org/10.1016/j.jhydrol.2012.06.043, 2012.

Charlier J.-B., Desprats J.-F., and Ladouche B.: Appui au SCHAPI 2014 – Module 1 – Rôle et contribution des eaux souterraines d'origine karstique dans les crues de la Loue à Chenecey-Buillon, BRGM/RP-63844-FR report, 109 pp., available at: https://infoterre.brgm.fr/rapports/RP-63844-FR.pdf (last access: 22 April 2020), 2014.

Charlier, J.-B., Moussa, R., Bailly-Comte, V., Danneville, L., Desprats, J.-F., Ladouche, B., and Marchandise, A.: Use of a flood-routing model to assess lateral flows in a karstic stream: implications to the hydrogeological functioning of the Grands Causses area (Tarn River, Southern France), Environ. Earth Sci., 74, 7605–7616, https://doi.org/10.1007/s12665-015-4704-0, 2015.

Charlier, J., Moussa, R., David, P., and Desprats, J.: Quantifying peakflow attenuation/amplification in a karst river using the diffusive wave model with lateral flow, Hydrol. Process., 33, 2337–2354, https://doi.org/10.1002/hyp.13472, 2019.

Cholet, C., Charlier, J.-B., Moussa, R., Steinmann, M., and Denimal, S.: Assessing lateral flows and solute transport during floods in a conduit-flow-dominated karst system using the inverse problem for the advection-diffusion equation, Hydrol. Earth Syst. Sci., 21, 3635–3653, https://doi.org/10.5194/hess-21-3635-2017, 2017.

Covino, T., McGlynn, B., and Mallard, J.: Stream-groundwater exchange and hydrologic turnover at the network scale: hydrologic turnover at the network scale, Water Resour. Res., 47, W12521, https://doi.org/10.1029/2011WR010942, 2011.

Delrieu, G., Nicol, J., Yates, E., Kirstetter, P.-E., Creutin, J.-D., Anquetin, S., Obled, C., Saulnier, G.-M., Ducrocq, V., and Gaume, E.: The catastrophic flash-flood event of 8–9 September 2002 in the Gard Region, France: a first case study for the Cévennes-Vivarais Mediterranean Hydrometeorological Observatory, J. Hydrometeorol., 6, 34–52, 2005.

De Waele, J., Martina, M. L. V., Sanna, L., Cabras, S., and Cossu, Q. A.: Flash flood hydrology in karstic terrain: Flumineddu Canyon, central-east Sardinia, Geomorphology, 120, 162–173, https://doi.org/10.1016/j.geomorph.2010.03.021, 2010.

Dingman, S. L.: Physical Hydrology, Third Edition, Waveland Press, Long Grove, Illinois, USA, 2015.
Eakin, T. E.: A regional interbasin groundwater system in the White River Area, southeastern Nevada, Water Resour. Res., 2, 251–271, https://doi.org/10.1029/WR002i002p00251, 1966.

Eckhardt, K.: How to construct recursive digital filters for baseflow separation, Hydrol. Process., 19, 507–515, https://doi.org/10.1002/hyp.5675, 2005.

Edijatno De Oliveira Nascimento, N., Yang, X., Makhlouf, Z., and Michel, C.: GR3J: a daily watershed model with three free parameters, Hydrol. Sci. J., 44, 263–277, https://doi.org/10.1080/02626669909492221, 1999.

Fan, Y.: Are catchments leaky?, WIRES-Water, 6, 6:e1386, https://doi.org/10.1002/wat2.1386, 2019.

Gárfias-Soliz, J., Llanos-Acebo, H., and Martel, R.: Time series and stochastic analyses to study the hydrodynamic characteristics of karstic aquifers, Hydrol. Process., 24, 300–316, https://doi.org/10.1002/hyp.7487, 2009.

Genereux, D. P., Jordan, M. T., and Carbonell, D.: A paired-watershed budget study to quantify interbasin groundwater flow in a lowland rain forest, Costa Rica: rainforest watershed budgets, Water Resour. Res., 41, W04011, https://doi.org/10.1029/2004WR003635, 2005.

Gustard, A., Bullock, A., and Dixon, J. M.: Low flow estimation in the United Kingdom, Institute of Hydrology, Wallingford, UK, 20–23, 1992.

Hakoun, V. and Manlay, A.: Rhydro: a set of tools to access discharge and water levels data from the French databases (Version v1), Zenodo, https://doi.org/10.5281/zenodo.3744183, 2020.

Hartmann, A., Weiler, M., Wagener, T., Lange, J., Kralik, M., Humer, F., Mizyed, N., Rimmer, A., Barberá, J. A., Andreo, B., Butscher, C., and Huggenberger, P.: Process-based karst modelling to relate hydrodynamic and hydrochemical characteristics to system properties, Hydrol. Earth Syst. Sci., 17, 3305–3321, https://doi.org/10.5194/hess-17-3305-2013, 2013.

Laaha, G. and Blöschl, G.: Seasonality indices for regionalizing low flows, Hydrol. Process., 20, 3851–3878, https://doi.org/10.1002/hyp.6161, 2006.

Le Bas, C.: "bdgsf_classe_ru.zip", Carte de la Réserve Utile en eau issue de la Base de Données Géographique des Sols de France, https://doi.org/10.15454/JPB9RB/HUY6TJ, Portail Data INRAE, V2, 2018.

Lebecherel, L., Andréassian, V., and Perrin, C.: On regionalizing the Turc-Mezentsev water balance formula: Turc-Mezentsev water balance formula, Water Resour. Res., 49, 7508–7517, https://doi.org/10.1002/2013WR013575, 2013.

Le Mesnil, M., Charlier, J.-B., Moussa, R., Caballero, Y., and Dörfliger, N.: Interbasin groundwater flow: Characterization, role of karst areas, impact on annual water balance and flood processes, J. Hydrol., 585, 124583, https://doi.org/10.1016/j.jhydrol.2020.124583, 2020.

Le Moine, N., Andréassian, V., Perrin, C., and Michel, C.: How can rainfall-runoff models handle intercatchment groundwater flows? Theoretical study based on 1040 French catchments: dealing with IGF in rainfall-runoff models, Water Resour. Res., 43, W06428, https://doi.org/10.1029/2006WR005608, 2007.

Le Moine, N., Andréassian, V., and Mathevet, T.: Confronting surface- and groundwater balances on the La Rochefoucauld-Touvre karstic system (Charente, France): closing the catchment-scale water balance – a case study, Water Resour. Res., 44, W03403, https://doi.org/10.1029/2007WR005984, 2008.

López-Chicano, M., Calvache, M. L., Martín-Rosales, W., and Gisbert, J.: Conditioning factors in flooding of karstic poljes – the case of the Zafarraya polje (South Spain), Catena, 49, 331–352, https://doi.org/10.1016/S0341-8162(02)00053-X, 2002.

Lvovich, M. I.: World water resources, present and future, GeoJournal, 3, 423–433, https://doi.org/10.1007/BF00455981, 1979.

Lyne, V. and Hollick, M.: Stochastic time-variable rainfall-runoff modelling, Institute of Engineers Australia National Conference, Vol. 1979, Institute of Engineers Australia, Barton, Australia, 1979.

Maillet, E. T.: Essais d'hydraulique souterraine & fluviale, A. Hermann, Paris, France, 1905.

Mallard, J., McGlynn, B., and Covino, T.: Lateral inflows, stream-groundwater exchange, and network geometry influence stream water composition, Water Resour. Res., 50, 4603–4623, https://doi.org/10.1002/2013WR014944, 2014.

Maréchal, J. C., Ladouche, B., and Dörfliger, N.: Karst flash flooding in a Mediterranean karst, the example of Fontaine de Nîmes, Eng. Geol., 99, 138–146, https://doi.org/10.1016/j.enggeo.2007.11.013, 2008.

Marty, R., Zin, I., and Obled, C.: Sensitivity of hydrological ensemble forecasts to different sources and temporal resolutions of probabilistic quantitative precipitation forecasts: flash flood case studies in the Cévennes-Vivarais region (Southern France): sensitivity of hydrological ensemble forecasts to PQPF, Hydrol. Process., 27, 33–44, https://doi.org/10.1002/hyp.9543, 2013.

Merz, R. and Blöschl, G.: Regionalisation of catchment model parameters, J. Hydrol., 287, 95–123, https://doi.org/10.1016/j.jhydrol.2003.09.028, 2004.

Moussa, R.: Analytical Hayami solution for the diffusive wave flood routing problem with lateral inflow, Hydrol. Process., 10, 1209–1227, https://doi.org/10.1002/(SICI)1099-1085(199609)10:9<1209::AID-HYP380>3.0.CO;2-2, 1996.

Moussa, R. and Bocquillon, C.: Criteria for the choice of floodrouting methods in natural channels, J. Hydrol., 186, 1–30, https://doi.org/10.1016/S0022-1694(96)03045-4, 1996.

Moussa, R. and Majdalani, S.: Evaluating lateral flow in an experimental channel using the diffusive wave inverse problem, Adv. Water Resour., 127, 120–133, https://doi.org/10.1016/j.adwwatres.2019.03.009, 2019.

Nathan, R. J. and McMahon, T. A.: Evaluation of automated techniques for base flow and recession analyses, Water Resour. Res., 26, 1465–1473, https://doi.org/10.1029/WR026i007p01465, 1990.

Nguyen, V. T., Dietrich, J., and Uniyal, B.: Modeling interbasin groundwater flow in karst areas: Model development, application, and calibration strategy, Environ. Modell. Softw., 124, 104606, https://doi.org/10.1016/j.envsoft.2019.104606, 2020.

Oudin, L., Andréassian, V., Perrin, C., Michel, C., and Le Moine, N.: Spatial proximity, physical similarity, regression and ungaged catchments: A comparison of regionalization approaches based on 913 French catchments, Water Resour. Res., 44, W03413, https://doi.org/10.1029/2007WR006240, 2008.

Parajka, J., Merz, R., and Blöschl, G.: A comparison of regionalisation methods for catchment model parameters, Hydrol. Earth Syst. Sci., 9, 157–171, https://doi.org/10.5194/hess-9-157-2005, 2005.

Perrin, J.-L. and Tournoud, M.-G.: Hydrological processes controlling flow generation in a small Mediterranean catch-
ment under karstic influence, Hydrol. Sci. J., 54, 1125–1140,
https://doi.org/10.1623/hysj.54.6.1125, 2009.

Piggott, A. R., Moin, S., and Southam, C.: A revised approach to the UKIH method for the calculation of baseflow, Hydrol. Sci. J., 50, 920, https://doi.org/10.1623/hysj.2005.50.5.911, 2005.

Raghunath, H. M.: Hydrology: Principles, Analysis And Design, New Age International, New Delhi, India, 2006.

Rutledge, A. T.: Computer programs for describing the recession of ground-water discharge and for estimating mean ground-water recharge and discharge from streamflow records: Update. No. 98, US Department of the Interior, US Geological Survey, 1998.

Satterthwaite, F. E.: An Approximate Distribution of Estimates of Variance Components, Biometrics Bull., 2, 110–114, https://doi.org/10.2307/3002019, 1946.

Scanlon, B. R., Mace, R. E., Barrett, M. E., and Smith, B.: Can we simulate regional groundwater flow in a karst system using equivalent porous media models? Case study, Barton Springs Edwards aquifer, USA, J. Hydrol., 276, 137–158, https://doi.org/10.1016/S0022-1694(03)00064-7, 2003.

Schaller, M. F. and Fan, Y.: River basins as groundwater exporters and importers: Implications for water cycle and climate modeling, J. Geophys. Res., 114, D04103, 2009.

Sivapalan, M., Yaeger, M. A., Harman, C. J., Xu, X., and Troch, P. A.: Functional model of water balance variability at the catchment scale: 1. Evidence of hydrologic similarity and space-time symmetry: space-time symmetry of water variability, Water Resour. Res., 47, W02522, https://doi.org/10.1029/2010WR009568, 2011.

Sloto, R. A. and Crouse, M. Y.: HYSEP: A Computer Program for Streamflow Hydrograph Separation and Analysis, ater-resources investigations report, vol. 96, p. 4040, 1996.

Tabary, P.: The New French Operational Radar Rainfall Product. Part I: Methodology, Weather Forecast., 22, 393–408, https://doi.org/10.1175/WAF1004.1, 2007.

Tang, W. and Carey, S. K.: HydRun: A Matlab toolbox for rainfall-runoff analysis, Hydrol. Process., 31, 2670–2682, https://doi.org/10.1002/hyp.11185, 2017.

Thiéry, D.: Logiciel GARDÉNIA, version 8.2, Guide d'utilisation, BRGM/RP-62797-FR report, 143 p., Orléans, available at: http://infoterre.brgm.fr/rapports/RP-62797-FR.pdf (last access: 9 March 2021), 2014.

Thornthwaite, C. W.: An Approach toward a Ratio-
nal Classification of Climate, Geogr. Rev., 38, 55–94,
https://doi.org/10.2307/210739, 1948.

Todini, E.: The ARNO rainfall – runoff model, J. Hydrol., 175, 339–382, https://doi.org/10.1016/S0022-1694(96)80016-3, 1996.

Tóth, J.: A theoretical analysis of groundwater flow in small drainage basins, J. Geophys. Res., 68, 4795–4812, https://doi.org/10.1029/JZ068i016p04795, 1963.

Vidal, J.-P., Martin, E., Franchistéguy, L., Baillon, M., and Soubey-
roux, J.-M.: A 50-year high-resolution atmospheric reanalysis
over France with the Safran system, Int. J. Climatol., 30, 1627–
1644, https://doi.org/10.1002/joc.2003, 2010.

Zanon, F., Borga, M., Zoccatelli, D., Marchi, L., Gaume, E., Bon-
nifait, L., and Delrieu, G.: Hydrological analysis of a flash
flood across a climatic and geologic gradient: The September
18, 2007 event in Western Slovenia, J. Hydrol., 394, 182–197,
https://doi.org/10.1016/j.jhydrol.2010.08.020, 2010.