Stream Hydrochemical Response to Flood Events in a Multi-Lithological Karstic Catchment from the Pyrenees Mountains (SW France)

: Hydrological variations hold a signiﬁcant inﬂuence over the water chemistry in the karstic critical zone. In this context, the karstic Baget Catchment (BC) has been monitored at a high resolution over two years at the outlet in order to set up a typology of the ﬂood events. The objective was to assess the multiple streamwater physico-chemical patterns in response to hydrological variations, streamﬂow component (quick-response, subsurface, and baseﬂow) and lithological contributions, and biogeochemical processes. The karstic catchment exhibited an impulsive response to ﬂood events in relation to the typical structural and morphological characteristics of the karst. In addition, this response was constrained by the magnitude of the rainfall and the preceding hydroclimatic conditions. The variability of the dissolved load in streamﬂow was closely associated with the characteristics of the weathered rocks and the hydrological conditions throughout the year. Two simple indicators allow to characterize the concentration–discharge relationships with different hysteresis patterns on a set of ﬂoods with various intensities and shapes of the hydrograph and under different hydrological conditions before the ﬂood. Almost all elements exhibited either clockwise loops or more complex behaviors, suggesting a higher overall concentration when the major water contribution comes from the quick-response ﬂow (karst and surface runoff fraction). Besides, the epikarst ﬂushing under dry conditions led counterclockwise hysteresis patterns for calcium (Ca 2+ ) and bicarbonate (HCO 3 − ) which revealed an overall chemostatic behavior as a result of carbonate dissolution in the karst. On the contrary, sulfate (SO 42 − ) exhibited the widest relative variation during ﬂooding and showed a signiﬁcant sensitivity to the dilution process with increasing discharge. For medium ﬂood episodes (Q max < 4.4 m 3 · s − 1 ), an overall concentration increase or chemostatic behavior could be observed during the rising limb of the hydrograph. On the opposite, under extreme ﬂood episodes (Q max > 8.3 m 3 · s − 1 ) occurring after several rain events, a dilution pattern was noted for all elements originating from rock weathering. Finally, high-frequency sampling during storm events improved the understanding of the factors controlling the hydrochemical dynamic in karstic catchments.

sic in the northern upper section bordered also by silicate rocks (Cretaceous flysch of Arbas-other lithologies-and Paleozoic Lachein schists) ( Figure 1C, Appendix A). Sulfide minerals (meanly pyrite, FeS) are associated as a minor phase to limestones and flyschs [14,36], while sulfate minerals (mainly gypsum, CaSO4·2H2O) from the Upper Triassic represent only 0.2% of the BC surface area. The experimental BC is part of the French Karst Observatory Network (SNO Karst) [40] and belongs to the observatories of the French Research Infrastructure called OZCAR (Observatoire de la Zone Critique: Application et Recherche) [41], included in the European Research Infrastructure eLTER (Long Term Ecosystem Research in Europe) [42,43].
The hydrochemical survey was carried out at the BC outlet (Baget outlet, Figure 1B): (i) a weekly sampling between 1978 and 2006 by the CNRS laboratory of Moulis [44]; (ii) every 6 months by BRGM Occitanie [45] from 2007 to 2014; and (iii) a biweekly sampling that started in 2014 and is currently ongoing by the Laboratoire écologie fonctionnelle et environnement [17]. In addition, high-frequency sampling was conducted between February 2019 and December 2020 during the flood events. Rainfall and rainwater samples were also collected monthly from September 2019 to December 2020 at the meteorological station of Balagué (42 • 58 06.5 N 1 • 00 39 E, elevation 658 m, Figure 1B). Rainfall between January 2019 and August 2019 was calculated using the daily rainfall registered at the closest Meteo-France meteorological station of St. Girons (elevation 414 m), situated 8.3 km downstream from the outlet of BC. The very significant relationship (R 2 = 0.80; p < 0.001; n = 365, with a slope of 1.39) between the rainfall at Balagué and at St. Girons stations was used for this calculation.
The physico-chemical parameters, such as water temperature (Water T • ), specific conductivity (Cond), dissolved oxygen (DO), and pH, were measured at the time of sampling using a portable multi-parameter (WTW probe, Xylem Analytics Germany Sales GmbH & Co., Weilheim, Germany) at the Baget outlet. This site is equipped with a mechanical limnigraph (OTT 20 1/5; Loveland, CO, USA) and a float-type water-level sensor (OTT Thalimedes; Loveland, CO, USA). Several gauges served to plot the rating curve and thus estimate the discharge from the water height. Discharge data at the outlet since 1968 are freely accessible [46]. In addition, since April 2014, water level and different physico-chemical parameters (Water T • , pH, Cond, Turbidity, DO, nitrate and chloride) were measured continuously at a high frequency (10-minute interval) in the mid-depth of the water column (at the outlet) using a multi-parameter probe (YSI 6920V2-01, YSI Incorporated, Yellow Springs, OH, USA) equipped with multiple sensors and an atmospheric pressure correction. Moreover, during the study period (2019-2020), water samples were collected every two hours using an automatic sampler and analyzed in the laboratory.
Water samples were collected in a clean and dry 1 L HDPE bottle, triple rinsed with streamwater before sampling, and stored in darkness at 5 • C until analyses. All samples were filtered using Millipore 0.22 µm cellulose nitrate membranes. Alkalinity was determined using a standard acid titration method with HCl 0.02N and a Metrohm titrant (716 DMS Titrino, Metrohm, Riverview, FL, USA). Major anion concentrations (SO 4 2− , NO 3 − , and Cl − ) were analyzed by ion chromatography (Dionex apparatus ICS 5000+, Thermo Fisher Scientific, Waltham, MA, USA). Major cation (Ca 2+ , Mg 2+ , Na + , and K + ) and silica concentrations were measured using an acidified sample with 16N HNO 3 , later analyzed by an inductively coupled plasma optical emission spectrometer (ICP-OES; Iris Intrepid II XLD, Thermo Electron, Thermo Fisher Scientific, Waltham, MA, USA). The total dissolved solids content (TDS) was calculated by adding up all the dissolved species mentioned above. Dissolved organic carbon (DOC) was measured on a Shimadzu TOC 5000 analyzer (Shimadzu Corporation, Kyoto, Japan), coupled to the CO 2 detection by an infrared gas analyzer (NDIR).
Regarding the quality assurance (QA), the protocols for sampling, filtration, pretreatment, and storage of the stream-and rainwater before analysis were always controlled. Concerning the quality control (QC), all the methods of analysis were controlled following French (NF), European (EN), and international (ISO) standards. In addition, the physicochemical analysis laboratory of the Laboratoire écologie fonctionnelle et environment takes part every year to the ICP-Water Intercomparison with the Norwegian Institute for Water Research (NIVA), which follows a QA/QC procedure.

Data Treatment 2.3.1. Hydrograph Separation Method
The streamflows were split into the quick-response flow (karst and surface runoff fraction), the subsurface flow (epikarst and infiltration fraction), and the baseflow (groundwater fraction) ( Figure 2).
The onset of a storm period was established when the discharge increased suddenly, exhibiting a marked inflection (Point A, Figure 2A). This sudden change was assessed when the discharge change rate (DCR, Equation (1)) was greater than 100%.
where Q i and Q i+1 are the daily discharge at time-step i and i + 1. The end of the storm period was determined when the DCR was greater than −20% and the base discharge was lower than 300 L·s −1 in BC (Point D, Figure 2A). Moreover, small floods leading to discharge variation only within ±1 m 3 ·s −1 were not considered.
into straight lines by using the logarithm of the discharge. Then, the descending limb corresponding to the recession period is divided into three straight lines with different slopes through the determination of two inflection points (points C and D, Figure 2A corresponding to the end of the quick-flow and subsurface-flow contribution, respectively [24]. Then, the separation line between the subsurface runoff with the quick-response runoff and with the baseflow was drawn by connecting points A-C and A-D, respectively.  (2)).
The embedded small graph in Figure 2B, exhibits the relationship between the discharge contribution of each flow component (%) and the total stream discharge during the recession period. The black points on the Figure 2B correspond to water sampling.

Analyses of C-Q Relationships
The temporal variations in both discharge and concentrations can be analyzed also using the relationships between concentrations and discharges, which generally follow a power-law function (Equation (3) where Qi and Ci are the instantaneous discharges and dissolved element concentrations, respectively, at time-step i. The exponent b is a useful indicator of the chemodynamic features of the dissolved element behavior. If b < 0, the dilution effect produced by Q is the  (2)).
The embedded small graph in Figure 2B, exhibits the relationship between the discharge contribution of each flow component (%) and the total stream discharge during the recession period. The black points on the Figure 2B correspond to water sampling.
The quantitative analysis of the recession curves is derived from the research of Maillet [47], who indicated that the recession curve of the hydrograph during a storm period can be characterized by an exponential decrease of the discharge (Equation (2)), suggesting a linear relationship between the piezometric head and flowrate [48]. (2) where Q t is the instantaneous discharge at time-step t (in days in our study), Q 0 is the initial discharge at time zero, and α is the recession coefficient. This coefficient depends on the characteristics of each flow associated with the reservoir it drains, allowing to differentiate the different streamflow components on the basis of their "α" coefficient. The hydrograph separation method proposed by Probst [24,49,50] was adapted for karst systems, assuming that the quick-response flow (karst and surface runoff fraction), the subsurface flow (epikarst and infiltration fraction), and the baseflow (groundwater fraction) occur simultaneously. The first step is to transform the curves of the hydrograph into straight lines by using the logarithm of the discharge. Then, the descending limb corresponding to the recession period is divided into three straight lines with different slopes through the determination of two inflection points (points C and D, Figure 2A corresponding to the end of the quick-flow and subsurface-flow contribution, respectively [24]. Then, the separation line between the subsurface runoff with the quick-response runoff and with the baseflow was drawn by connecting points A-C and A-D, respectively.

Analyses of C-Q Relationships
The temporal variations in both discharge and concentrations can be analyzed also using the relationships between concentrations and discharges, which generally follow a power-law function (Equation (3)) [51]: where Q i and C i are the instantaneous discharges and dissolved element concentrations, respectively, at time-step i. The exponent b is a useful indicator of the chemodynamic features of the dissolved element behavior. If b < 0, the dilution effect produced by Q is the major regulator on C; if b = 0, C exhibits a chemostatic behavior, without influence of Q variation; if b > 0, an enrichment takes place on the solutes when the discharge increases.

Hysteresis Analyses
During a flood, the C-Q relationship can take the shape of a straight line or a loop, called hysteresis [52]. This hysteresis reflects the different concentration responses to the same discharge during the rising and falling limbs of the hydrograph.
To characterize the hysteresis patterns, we used two simple descriptors that describe the relative changes in dissolved concentrations (∆C), the direction of rotation, and the relative amplitude of the loop pattern (∆I).
where C b and C p are the instantaneous dissolved element concentrations at the base flow and during the peak flow of the flood event, respectively. Positive and negative values of ∆C indicate a concentration increase or dilution, respectively [53]. ∆C values ranging between −10 and 10 can be considered as corresponding to a neutral response discharge increase.
The normalized data were also used to compare the hysteresis behavior of each parameter within and between different flood events (Equations (5) and (6)).
where Q i and C i are the instantaneous discharge and dissolved element concentrations at time-step i, respectively. The Q min , Q max and C min , C max correspond to the minimum and maximum values of the discharges and concentrations during a complete flood event, respectively. The relative amplitude of the loop pattern (∆I) can be calculated as follows: where C i−RL_norm and C i−FL_norm refer to the normalized concentrations on the rising limb (RL) and on the corresponding falling limb (FL) of the hydrograph, respectively, for each 25% increase in discharge [54,55]. Positive and negative values of ∆I indicate clockwise or counterclockwise loops, respectively. For ∆I values between −10 and 10, one can consider that there is no loop or ambiguous rotational response to increasing discharge.

Method of Flux Calculation
The total flux (F total ) was estimated by adding all fluxes corresponding to each flood event, estimated between two sampling intervals (i to i + 1) using the discharge-weighted average concentration (C j ) and the water volume (V j ), according to the following equations: with C i : concentration at time i; Q i : discharge at time i; and V j : volume of water flow between i and I + 1.

Hydrochemical Survey during the Period 2019-2020
The summary of the hydrochemical and physicochemical characteristics of the streamwater are reported in Table 1. All parameters reacted to different extents to the hydrological variations, as indicated by their coefficient of variation (CV = σ/X). The discharge, water, and air temperatures ( Figure 3) exhibited an obvious seasonality during the two years of survey, but their CV is very different: the highest variation was for the discharge (CV = 1.03) and the lowest one was for the water temperature (CV = 0.02). The study period is 20% wetter (Q = 0.53 m 3 ·s −1 ) compared to mean of the last 40 years (Q = 0.44 m 3 ·s −1 ) [17]. The total rainfall amount was 1543 mm and 1632 mm for the years 2019 and 2020, respectively. The mean volume of rainfall was 21.0 × 10 6 m 3 per year. Open field precipitation is circumneutral, and very low concentrated in the major dissolved elements (<25 µmol·L −1 ). Na + and Cl − are the dominant ions whereas the Mg 2+ , K + , SO 4 2− , and NO 3 − concentrations are very low (mean concentration <6 µmol·L −1 , Appendix B). The streamwater in BC is mildly alkaline and pH showed low fluctuations (7.9 ± 0.2). The water temperature is highly stable (10.1 ± 0.2 • C) throughout the study period. Conductivity is related to total dissolved solids (R 2 = 0.4, p < 0.01). The mean concentration of the different elements in the streamwater exhibited the following abundance order during the study period: for cations, calcium (Ca 2+ , 88.8% of total cationic concentration in µeq·L −1 ) >> magnesium (Mg 2+ , 9.7%) >> sodium (Na + , 1.1%) > potassium (K + , 0.4%); and for anions, bicarbonate (HCO 3 − , 93.5% of total anionic concentration in µeq·L −1 ) >> sulphate (SO 4 2− , 4.4%) > chloride (Cl − , 1.3%) > nitrate (NO 3 − , 0.8%). The predominance of Ca 2+ , Mg 2+ , and HCO 3 − , which represent 96% of the ionic charge, is related to the dominant carbonate lithology of the BC [14,17]. The discharge-weighted mean concentration of sulphate has significantly decreased during the study period (0.15 µeq·L −1 ) compared to the period 1978-2018 (0.35 µeq·L −1 ) [17]. In addition, Ca 2+ and HCO 3 − were the most stable ions, while SO 4 2− , K + , NO 3 − , and DOC exhibited a higher CV (Table 1). The HCO 3 − is the dominant species of DIC according to the water pH and its concentration is much higher than the world river average concentration of 850 µmol·L −1 [56].  (a) This mean discharge was calculated using the sampling interval method (see Equation (9)) for the entire sampling period. (b) TDS = Ca 2+ + Mg 2+ + Na + + K + + HCO3 − + SO4 2− + Cl − + NO3 − + H4SiO4.

Figure 3.
High frequency temporal series of the stream discharge (Q), rainfall (P), and physicochemical parameters (air temperature (air T°), water temperature (water T°), pH, conductivity (Cond), and TDS (total dissolved solids)) during the two years 2019 and 2020. Black solid lines represent the 5-day daily moving average for air and water temperature. Dots represent punctual streamwater sampling during the study period (N = 297). The numbers circled refer to the major floods sampled during the study period. Flood events 2, 4, 5, 9, 10, and 11 (light blue circles) were considered further as typical examples. High frequency temporal series of the stream discharge (Q), rainfall (P), and physicochemical parameters (air temperature (air T • ), water temperature (water T • ), pH, conductivity (Cond), and TDS (total dissolved solids)) during the two years 2019 and 2020. Black solid lines represent the 5-day daily moving average for air and water temperature. Dots represent punctual streamwater sampling during the study period (N = 297). The numbers circled refer to the major floods sampled during the study period. Flood events 2, 4, 5, 9, 10, and 11 (light blue circles) were considered further as typical examples.
The concentrations of the major elements, represented by TDS, show significant temporal variations in the BC ( Figure 3). The average TDS (285.5 mg·L −1 ) is around 3 times the world river TDS average (97 mg·L −1 ) [56]. Compared to other rivers draining carbonates in the world, the mean TDS in BC is higher than those of the Wujiang River in the Yunnan-Guizhou plateau (265 mg·L −1 ) [35], the Ganges and Indus rivers in the Himalayas (164 mg·L −1 ) [57,58], or the Mackenzie River in the Rocky Mountains (160 mg·L −1 ) [59]. However, compared to rivers in temperate karst terrains, the TDS values measured in BC are similar to those measured in rivers of the Alps Mountains (253 mg·L −1 during spring and 378 mg·L −1 during winter) [60] and lower than in the Jura Mountains (378 mg·L −1 ) [61].
The Net Inorganic Charge Balances (NICB % = 200 × (TZ + − TZ − )/ (TZ + + TZ − )) between the sum of anions (TZ − ) and the sum of cations (TZ + ) (both in meq·L −1 ) are on average 4.0 ± 3.5%, with 70% of the samples having an NICB of less than 5%. This indicates that the contribution of unanalyzed ions, such as organic anions, play a negligible role in the water chemistry.

Flood Hydrological Characteristics
During the study period, eleven hydrological events were observed ( Figure 3, Table 2). The duration of these events was variable but closely within the range of influence on discharge, between 10 and 20 days, as shown by Mangin [36] for the BC. Table 2. Summary of the mean hydrological characteristics for each flood event. Q and Qs refer to discharge (in m 3 ·s −1 ) and specific discharge (in mm). "Time before" refers to the period between the last flood event and the beginning of the current flood event analyzed, "Q mean-24h " to the average discharge the day before the start of the flood event, "Time start Q" to the date and hour of the beginning of the flood, "Start Q" to the discharge value at the beginning of the flood, "Peak Q" to the discharge value at the peak flow, "End Q" to the discharge value at the end of the flood, "Duration" to the total event, and "Total Qs" to the average specific discharge during the flood. For time and duration, the unit is in hours.

Flood
Time The contrasted flood events, 2, 9, 10, and 11 ( Figure 3, Table 2), were analyzed in more detail. These flood events were chosen based on (i) the period of the hydrological year: beginning of the hydrological cycle (floods 9 and 10), first high-water level period (flood 11), and the last flood of the hydrological cycle (flood 2); (ii) the intensity/magnitude of the flood event reflected by the differentiated discharge peaks (2.13, 2.29, 4.42, and 8.25 m 3 ·s −1 , respectively) and specific discharges (44.6, 63.7, 92.5, and 165.3 mm, respectively); (iii) similar discharge values (Q mean-24h ) preceding the flood event; and (iv) the potential influence of past and/or successive floods.
During the last three studied flood events (9, 10, and 11), the pH of the open field precipitation was between 6 and 7, and Ca 2+ is a major element in the rainwater together with Na + and Cl − ions (<25 µmol·L −1 ). Unlike flood events 10 and 11, flood event 2 followed a relatively dry period straight after a major flood, such as for floods 1 (discharge peak of 7.0 m 3 ·s −1 ) and 4 (discharge peak of 9.4 m 3 ·s −1 ). Flood 9 was preceded by the smallest flood identified (discharge peak of 1.9 m 3 ·s −1 with a total duration of 5 days) whose influence on the water chemistry was minor.

C-Q Relationships
The control of discharge over the dissolved element concentrations can be observed by comparing the temporal variations inof TDS and of discharge (Figure 3), and by the relationships between the concentrations and discharges (C-Q) ( Figure 4).
The C-Q relationships indicated consistent patterns between some elements, such aslike Ca 2+ and HCO 3 − inon one way, and Cl − and K + inon the other way. A similar pattern was also detected between Mg 2+ , Si, and Na + , buteven slightly less evident. Nevertheless, the global relationships were not significant for Ca 2+ , HCO 3 − , Cl − , K + , + and NO 3 − , showing globally a chemostatic behavior during the whole period. On the contrary, the C-Q relationship for SO 4 2− is significant and follows a dilution curve. The slope (b = −0.45) of the power-law function (C = aQ b ) for SO 4 2− tends to a pure dilution curve (b = −1). The Mg 2+ (b = −0.07), Si (b = −0.07), and Na + (b = −0.05) concentrations underwent a significant dilution. The DOC was the only parameter with a clear enrichment with increasing discharge, exhibiting a positive significant relationship (b = 0.19). Among those global trends, specific patterns can however be observed during the flood events ( Figure 4). For instance, flood events 10 (October) and 11 (December) exhibited concentrations of Ca 2+ and HCO 3 − much higher ( Figure 4) than during flood events 5 (February) and 2 (May). Flood event 5 displayed the lowest concentration of Na + and DOC but the highest concentration of Cl − , K + , and Si ( Figure 4).  (3)). The studied floods events are mentioned, as well as the background data discriminated according to four periods.

C-Q Patterns for Representative Floods
Considering all the data together, no clear relationships could be observed between C and Q ( Figure 4), but single flood events exhibited more obvious C/Q patterns ( Figure  5). The overall patterns can be illustrated for flood events 2 and 11 ( Figure 5), which in essence represented similar initial hydrological conditions observed for flood 9 and for floods 4, 5, and 10.
Loop trajectories of Ca 2+ and HCO3 − exhibited clockwise (floods 4, 5, 10, and 11) but also counterclockwise responses (floods 2 and 9). For the other dissolved elements, loop trajectories displayed clockwise responses or more complex behaviors, such as a figureof-eight shape ( Figure 5).  (3)). The studied floods events are mentioned, as well as the background data discriminated according to four periods.

C-Q Patterns for Representative Floods
Considering all the data together, no clear relationships could be observed between C and Q ( Figure 4), but single flood events exhibited more obvious C/Q patterns ( Figure 5). The overall patterns can be illustrated for flood events 2 and 11 ( Figure 5), which in essence represented similar initial hydrological conditions observed for flood 9 and for floods 4, 5, and 10.
ing concentrations up to about 2.5 m ·s and then increasing concentrations ( Figure 5B and 4C). This behavior was also observed for flood 4, unlike floods 5, 9, and 10 whose peak discharges were lower than 2.5 m 3 ·s −1 .
An overall dilution was observed for the SO4 2− concentration ( Figure 5D), except during the first peak of the minor flood event 11 as well as in flood 5 and 10 (Qpeak < 1.5 m 3 ·s −1 ). The DOC concentration showed a clear increase with discharge for all flood events ( Figure  5F), except for flood 10.

Typology of Flood Events
Based on the hydrological events sampled in BC, the main hydrochemical patterns were characterized during flood events, depending on the initial hydrological conditions before the flood, the intensity, and the shape of the flood hydrograph (Table 3). Loop trajectories of Ca 2+ and HCO 3 − exhibited clockwise (floods 4, 5, 10, and 11) but also counterclockwise responses (floods 2 and 9). For the other dissolved elements, loop trajectories displayed clockwise responses or more complex behaviors, such as a figure-of-eight shape ( Figure 5).
A decrease in Ca 2+ and HCO 3 − concentrations with increasing discharges was observed during the major peak flow (8.3 m 3 ·s −1 ) of flood 11. Similar behavior was observed during the major peak flow (9.4 m 3 ·s −1 ) of flood 4. For the other floods, such as flood 2, a slight variation was observed in the Ca 2+ and HCO 3 − concentrations ( Figure 5A). During flood events 2 (peak Q = 4.4 m 3 ·s −1 ) and 11 (peak Q = 8.3 m 3 ·s −1 ), the Na + , Si, Mg 2+ , and K + concentrations showed a similar trend during the recession period: decreasing concentrations up to about 2.5 m 3 ·s −1 and then increasing concentrations ( Figures 4C and 5B). This behavior was also observed for flood 4, unlike floods 5, 9, and 10 whose peak discharges were lower than 2.5 m 3 ·s −1 .
An overall dilution was observed for the SO 4 2− concentration ( Figure 5D), except during the first peak of the minor flood event 11 as well as in flood 5 and 10 (Q peak < 1.5 m 3 ·s −1 ). The DOC concentration showed a clear increase with discharge for all flood events ( Figure 5F), except for flood 10.

Typology of Flood Events
Based on the hydrological events sampled in BC, the main hydrochemical patterns were characterized during flood events, depending on the initial hydrological conditions before the flood, the intensity, and the shape of the flood hydrograph (Table 3). Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak.

Initial dry Hydrological Conditions
The duration of the previous recession period was higher than 95 days and lower than 35 days for the dry and wet hydrological conditions, respectively.
Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and rela-tively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak.

Initial dry Hydrological Conditions
The duration of the previous recession period was higher than 95 days and lower than 35 days for the dry and wet hydrological conditions, respectively.
Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and rela-tively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak.

Initial dry Hydrological Conditions
The duration of the previous recession period was higher than 95 days and lower than 35 days for the dry and wet hydrological conditions, respectively.
Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
Hysteresis loops The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak.

Initial dry Hydrological Conditions
The duration of the previous recession period was higher than 95 days and lower than 35 days for the dry and wet hydrological conditions, respectively.
Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak.

Initial dry Hydrological Conditions
The duration of the previous recession period was higher than 95 days and lower than 35 days for the dry and wet hydrological conditions, respectively.
Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak.

Initial dry Hydrological Conditions
The duration of the previous recession period was higher than 95 days and lower than 35 days for the dry and wet hydrological conditions, respectively.
Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak.

Initial dry Hydrological Conditions
The duration of the previous recession period was higher than 95 days and lower than 35 days for the dry and wet hydrological conditions, respectively.
Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak.

Initial dry Hydrological Conditions
The duration of the previous recession period was higher than 95 days and lower than 35 days for the dry and wet hydrological conditions, respectively.
Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
Na + , Si, Mg 2+ The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak.

Initial dry Hydrological Conditions
The duration of the previous recession period was higher than 95 days and lower than 35 days for the dry and wet hydrological conditions, respectively.
Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak.

Initial dry Hydrological Conditions
The duration of the previous recession period was higher than 95 days and lower than 35 days for the dry and wet hydrological conditions, respectively.
Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak.

Initial dry Hydrological Conditions
The duration of the previous recession period was higher than 95 days and lower than 35 days for the dry and wet hydrological conditions, respectively.
Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

K +
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

SO 4 2−
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
Cl − , NO 3 − The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

DOC
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO3 − concentrations and clockwise variations for Na+, Si, and Mg 2+ . The SO4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with nonsignificant concentration variation. Meanwhile, SO4 2− , K + , Cl − , and NO3 − showed a concentration increase response, unlike the DOC concentration. Table 3. Summary of the main C-Q responses in BC according to distinct hydrological conditions. The schematic representation of the hysteresis loop corresponds to the C-Q behavior during the major peak. Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
(a) The duration of the previous recession period was higher than 95 days and lower than 35 days for the dry and wet hydrological conditions, respectively.
The initial dry conditions correspond to long periods of very low rainfall and relatively high evapotranspiration. The floods occurring after these dry conditions are typically of short duration and intensity. Their chemical response to increasing discharge is characterized by slight anti-clockwise variations for the Ca 2+ and HCO 3 − concentrations and clockwise variations for Na + , Si, and Mg 2+ . The SO 4 2− concentration is systematically diluted while K + and DOC are concentrated. Cl − and NO 3 − exhibit different clockwise behaviors, depending on the intensity of the flooding.
During the initial wet conditions, the C-Q relationships displayed an overall dilution (except for DOC) if the intensity of the flood event is extremely high associated with a long period of intermittent rainfalls before the major peak. Under low flooding intensity, Ca 2+ , HCO 3 − , Na + , Si, and Mg 2+ concentrations displayed a clockwise hysteresis with non-significant concentration variation. Meanwhile, SO 4 2− , K + , Cl − , and NO 3 − showed a concentration increase response, unlike the DOC concentration.
Based on the C/Q patterns ( Figure 5) and on the typology of flood events (Table 3), we selected the representative flood events 2 and 11 for more in-depth analysis. These two floods represent the hydrological and hydrochemical responses under dry and wet hydrological conditions, respectively. Besides, the shapesshape of these hydrographs are less complex and allows better separation of the different streamflow components.

Hydrochemical Patterns and Streamflow Components during Floods 2 and 11
HydrologicalThe hydrological event 2 ( Figure 6A) started on 24 May and ended 10 days later, with most rainfall on the first day. The major peak flow was reached 25 hours after the rainfall peak in hourly precipitation (5.2 mm on 24 May, 06:00). Several small rainfall events perturbed the discharge during the recession limb of the flood event.

Sources and Major Processes Controlling Streamwater Chemistry
The streamwater chemical signature is derived from atmospheric, biological, and anthropogenic inputs and rock weathering [62,63], with the carbonate dissolution as the dominant process controlling BC streamwater chemistry [17]. The C-Q ( Figure 4) and flux- The hydrograph separation of flood event 2 (Figures 2A and 6A) exhibited three distinct recession coefficients (α), respectively, for the quick-response flow (α q : 0.64 to 0.66), the subsurface flow (α s : 0.19), and the baseflow (α b : 0.11 to 0.12). The analysis of the other flood events in BC (Figure 3) showed that the different α values (with t in days) ranged between 0.40 and 0.80 for α q , from 0.10 to 0.20 for α s , and between 0.04 and 0.12 for α b .
During the first part of the recession, the quick flow-response represented the main contribution (>50%) to streamflow during one-third (3.5 days) of the flooding episode up to a stream discharge higher than 1.5 m 3 ·s −1 . Its contribution decreased progressively while the subsurface flow became dominant for a stream discharge close to 1 m 3 ·s −1 ( Figure 2B, Figure 6A). Likewise, the baseflow dominance corresponds to a stream discharge lower than 0.65 m 3 ·s −1 .
Water temperature and turbidity increased (up to 0.4 • C and 180 NTU, respectively) with increasing discharge and then returned to the initial values. The increasing discharge caused a sharp decrease in conductivity (up to a minimum of 297 µS·cm −1 ) followed by a return to previous conditions, except after the last rainy episode when it showed a positive synchronism as discharge increases (29 May). The pH presented a small increase with the main discharge peak and remained stable until 29 May, and then it gradually decreased. During the first part of the falling limb, the Ca 2+ and HCO 3 − concentrations and the molar ratio Ca/Na increased, unlike the Mg 2+ , SO 4 2− , and Si concentrations and the molar ratio SO 4 /HCO 3 . The concentrations of the first group of elements decreased or stabilized till the last rainfall episode of 29 May, after which the Ca 2+ and HCO 3 − concentrations decreased sharply. For the second group of elements, they increased from 27 May onwards, when the contribution of the quick-response flow decreased markedly ( Figure 6A). A peak in the DOC concentration is observed with the peak flow, followed by a slight and regular decrease.
Hydrological event 11 started on 6 December till 20 December 2020 ( Figure 6B). During this event, two minor flow peaks preceded the major flow peak. The first flow peak was reached 6 hours after the first highest hourly rainfall (7.8 mm on 10 December, 17:00), while the major flow peak was reached 8 hours after the second highest hourly rainfall (7.8 mm on 12 December, 10:00). Note that a long period of light rains (42% of total rainfall event during 4.5 days) preceded these rainfall and flow peaks. Contrary to event 2, the recession limb of this flood event was not perturbed by new rainfall events. The quick-response flow was the main compartment (>50% of total stream discharge) during more than one-third (5.9 days) of the flooding episode up to a stream discharge around 1.2 m 3 ·s −1 . The baseflow dominance corresponds to a stream discharge lower than 0.75 m 3 ·s −1 ( Figure 6B).
The turbidity increased up to 136 NTU as the water level rose and then it decreased back to the initial conditions. Water temperature decreased with the first discharge peak and showed a clear intermittent variation until 17 December. From the beginning of the storm to the higher peak flow, conductivity had an intermittently decreasing trend (after an initial sudden increase), but with a more significant decrease after the first peak flow (11 December). The overall trend of the main dissolved element concentrations is downward. Silica and Na + increased again when the baseflow contribution increased. The DOC concentration increased smoothly with a maximum at the main peak flow and then returned intermittently to the initial values ( Figure 6B).

Sources and Major Processes Controlling Streamwater Chemistry
The streamwater chemical signature is derived from atmospheric, biological, and anthropogenic inputs and rock weathering [62,63], with the carbonate dissolution as the dominant process controlling BC streamwater chemistry [17]. The C-Q ( Figure 4) and flux-Q (Figure 7) relationships can be used to identify changes in sources and to assess the sensitivity of the solute concentrations and fluxes to runoff during changing hydrological conditions [64]. A general concentration increase or chemostatic pattern (except for SO 4 2− ) could be observed in the rising limb of the flood episodes, such as in floods 2, 4-the first flood peak-9, or 10 (Q max < 4.4). Whereas a dilution behavior was observed for all dissolved elements coming from the rock weathering in extreme flood episodes, such as floods 4-second flood peak-or 11 (Q max > 8.3) ( Figure 5). These patterns suggest that the dilution effect due to the discharge increase can be affected by other processes, supplying additional and variable inflow contributions of the solute elements with changes in discharge [65]. Among the power-law functions that can be calculated in BC, the lowest b-value corresponds to the C-Q relationship for SO 4 2− , which potentially comes from the gypsum dissolution and to a lesser extent from pyrite oxidation [17]. Although gypsum, probably the most soluble mineral present in BC, represents 66% of the SO 4 2− in streamwater [39], its limited point distribution in the catchment promotes a strong dilution effect. Other elements exhibited a slight dilution, such as Mg 2+ , Na + , and Si, which come from the dolomite dissolution (Mg 2+ ) and the silicate weathering (Mg 2+ , Na + , and Si). The dilution is more noticeable when the discharge increases for Mg 2+ and Si as opposed to Na + (Figure 4). The explanation lies in that rainfall (mean Na + rain = 21.8 µmol·L −1 ) can be also considered as an additional important source for Na + in streamwater (mean Na + streamwater = 40.4 µmol·L −1 ), unlike the other two elements. Another ion coming from precipitation is Cl − , whose origin is mainly rainfall, since no evaporitic deposit of halite (NaCl) have been identified in BC [38].
The main products of carbonate dissolution, Ca 2+ and HCO 3 − , showed an overall chemostatic behavior as the discharge increases ( Figure 4). Such behavior could initially be based on a counterbalance between their low concentration in rainwater input and the rapid kinetic of carbonate dissolution with biogenic soil CO 2 [29], indicating control by a process-limited regime. The high concentration of Ca 2+ and HCO 3 − at the beginning of the hydrological year (flood events 10 and 11) suggested a piston effect of more concentrated water, probably stored in the epikarst before the flood event ( Figure 4).
DOC is an important component of the terrestrial and aquatic ecosystem dynamics and is affected by the temperature, the soil leaching, and the streamwater flow pathways associated with the discharge changes [35]. In the BC streamwater, DOC exhibited the highest b-value in its C-Q relationship ( Figure 4) and a clear upward trend in the DOC concentration ( Figure 5) and flux (Figure 7) as the discharge increases. This suggests multiple sources mobilized during the flood events and a regulation by a transport-limited regime. Ladouche et al. [20] observed that the streamwater DOC response corresponds to the leaching of the upper soil layer or saturated areas containing organic matter compounds in forested catchments. In addition, the enrichment in K + as discharge increases could be due to cation exchange in soils [66] and to its high mobility/availability within the biomass [67].
Fluxes of Ca 2+ , HCO 3 − , Na + , and Si showed a strong correlation with discharge, suggesting that hydrological conditions dominate rock weathering fluxes (Figure 7). The sensitivity of the dissolved element concentrations to discharge changes is illustrated by the power law exponents "b" of the F-Q relationships, ranging from −1 (dilution) to 0 ("chemostasis"). Ca 2+ and HCO 3− fluxes lie just above the power law line 0, indicating higher carbonate dissolution rates as discharge increases for all flood events. The lower temperatures in winter and spring during the periods of the highest flows have increased the carbonate dissolution kinetics, leading to a greater amount of Ca 2+ and HCO 3 − released into the stream [26]. In addition, during flooding, physical erosion of rocks and soils released fine carbonate particles that could be dissolved in surface runoff during their transfer to the streamwater, and participate in the Ca 2+ and HCO 3 − concentration increase [12,64]. The Na + and Si values lie mainly below the power law line 0, indicating a slight dilution effect, which depends on the discharge but also on the preceding hydrological conditions. For instance, Si fluxes displayed clear dilution during flood events with high peak discharges (flood events 2 and 11). Conversely, Na + fluxes exhibit chemostatic behavior for flood events (10 and 11) occurring at the beginning of the hydrological year, in contrast to the dilution effect observed for those occurring in the middle (flood event 5) or at the end of the hydrological year (flood event 2). The dilution effect also observed for Mg 2+ (flood events 2 and 11) and more obviously for SO 4 2− fluxes (all flood events) is due to their continuous and rapid mobilization, which contributed to deplete their source reservoirs [29]. The higher concentrations of K + and Cl − during flood 5 ( Figure 4) lead to higher fluxes of these dissolved elements with slopes greater than C b (baseflow concentration) or even C m (mean concentration) (Figure 7). This suggests an origin of these elements from forest recycling [68], coupled with a release from some punctual livestock activities in the BC [39] that accumulated in the soil. to their continuous and rapid mobilization, which contributed to deplete their source reservoirs [29]. The higher concentrations of K + and Cl − during flood 5 ( Figure 4) lead to higher fluxes of these dissolved elements with slopes greater than Cb (baseflow concentration) or even Cm (mean concentration) (Figure 7). This suggests an origin of these elements from forest recycling [68], coupled with a release from some punctual livestock activities in the BC [39] that accumulated in the soil. Cb is the concentration when the discharge and flux were minimal during the study period (2019-2020). The dashed brown line is the trend line for the streamwaterstream water samples, whose slope is the mean concentration of each dissolved element (Cm). When Cm > >Cb, the cluster of points is above the line 1:1 (DOC, K + , Cl − , Ca 2+ , + and HCO3 − ); − ), on the contrary, when Cm<Cb, the cluster of points is below the line (SO4 2− , Na + , Mg 2+ , + and Si).

Contribution of the Different Reservoirs
The karstic catchments show an impulsive response to flood events [34]. The graphical hydrograph separation that waswhich has been performed allowed to quantify the contributions of the different streamflow components (quick-response, subsurface, and The dashed brown line is the trend line for the streamwaterstream water samples, whose slope is the mean concentration of each dissolved element (C m ). When Cm > >Cb, the cluster of points is above the line 1:1 (DOC, K + , Cl − , Ca 2+ , + and HCO 3 − ); − ), on the contrary, when Cm<Cb, the cluster of points is below the line (SO 4 2− , Na + , Mg 2+ , + and Si).

Contribution of the Different Reservoirs
The karstic catchments show an impulsive response to flood events [34]. The graphical hydrograph separation that waswhich has been performed allowed to quantify the contributions of the different streamflow components (quick-response, subsurface, and baseflow) to the total discharge (Appendix A). The contribution of the quick-response flow (karst flow and surface runoff fraction) represents 63% and 75% for floods 2 and 11, respectively. The magnitude of the rainfall and preceding hydroclimatic conditions play a pivotal role in the karst systems (Table 3) [34]. The first rainfalls of the hydrological year are mainly stored in the soil (or the epikarst) due to the preceding dry summer conditions (flood events 3, 8, and 9). For instance, the rainfall amount was measured at 212.0 and 105.2 mm during floods 3 (October 2019) and 8 (September 2020), while the total streamflow was estimated at 76.3 and 19.1 mm, i.e., only 36% and 18% of the total rainfall, respectively. In addition, Mangin [36] reported that the flood events around March reconstitutes the water reserves in this karst system. This replenishment of water reserves in the aquifer can be observed during the second period of the high waterflow level, such as flood events 5 and 6, whose total streamflow represented 67% and 73% of the total rainfall.
The typology of the flood events initially proposed in Table 3 was improved by using two simple indicators. They allow to quantify the shape of the C-Q relationships, based on relative changes in dissolved concentrations (∆C) and on the rotation direction and the relative amplitude of the loop pattern (∆I) (Figure 8). For most elements and flood events, the main loop rotation direction is clockwise (∆I > 10), due to a higher concentration during the rising limb of the hydrograph (when the major water contribution comes from the quick-response flow) than during the falling limb. However, for Ca 2+ and HCO 3 − , the C-Q relationship exhibited an anticlockwise loop pattern (∆I < −10) under dry hydrological conditions, indicating higher concentrations during the falling limb of the hydrograph (when the major water contribution comes from the subsurface flow). The contribution of the epikarst to this subsurface flow explains this hydrochemical response. Indeed, water can be stored and concentrated into the epikarst and then flushed with the first rain events [4,[68][69][70][71][72].
During flood event 11, a general dilution (∆C < −10) is observed for all elements, except for DOC. This dilution can be attributed to high discharges with intermittent rainfall inputs ( Figure 6B) combined with a continuous concentration decrease in the different hydrological reservoirs for more than 4 days before the major discharge peak (pre-flood event in Figure 6B). During the other floods, Na + and Si originating both from silicate weathering and Mg 2+ from dolomite dissolution exhibited a similar response pattern ( Figure 8) due to the main location of their lithological sources in the downstream part of the catchment, which is the most hydrologically active, to the north (Lachein tributary to the north of the La Hille losing stream) and south (Cayssau and Lasquert creek) of the Lachein stream ( Figure 1). Furthermore, the increase of K + concentration (Figure 8) may be related to its biogenic source [67,68] but also to the same order of magnitude of its concentration in rainwater (mean K + rain = 5.5 µmol·L −1 ) and streamwater (mean K + streamwater = 13.6 µmol·L −1 ) (Appendix B).
The in-depth investigation of the two major flood events (2 and 11) with different initial hydrological conditions, selected on the basis of the typology of flood events (Table 3, Figure 8), allowed to identify the contribution of the different compartments to the streamflow; i.e., surface runoff, karst, epikarst, infiltration, and baseflow. Their respective contributions during each flood event were evaluated regarding four representative hydrological periods (pre-flood event, rising period, recession period, and post-flood event) by coupling a detailed hydrograph with discrete sampling chemographs ( Figure 6). several Devon catchments and by other authors on several calcareous basins [17,29,34,74], a mobilization of solute elements stored during the summer in the karst system and released by the flushing effect explained the variation in the conductivity and concentration of these elements (coming from the rock weathering) as discharge slightly increases. During the rising period (RiP), up to 94% of the total streamflow is supplied by the quick-response flow ( Figure 6). The sudden increase in discharge was accompanied by an increase in turbidity ( Figure 6) and DOC (Table 3, Figure 8), whose peaks were almost simultaneous. The water entering the system, i.e.,. rainfall, can modify the streamwater temperature, increasing (flood 2, spring) or decreasing this signal (flood 11, autumn). This was illustrated by an upward water temperature signal during the flood event 2 and downward values during the flood event 11. Conductivity shows two antagonist patterns as discharge increases. On the one hand, a dilution effect by weakly concentrated rainwater in a period of low CO2 production and therefore of reduced rock weathering [14,17] explain the reduction in streamwater conductivity (flood 2). On the other hand, a pistonflow behavior on the last reserves of concentrated water mobilized during the rising limb of the first major peak (10 December, 22:30) accounts for the slight increase in streamwater conductivity (flood 11). However, these reserves are then depleted, and a general dilution occurs as the discharge continues to increase. This is demonstrated by a global drop in conductivity up to the highest peak flow (−9% regarding the initial value) and a decrease in concentrations of dissolved elements originating from rock weathering: Ca 2+ (ΔC = −13%), Mg 2+ (−45%), HCO3 − (−9%), SO4 2− (−85%), and Si (−14%). The dilution by continuous heavy precipitation is more effective on inputs supplied by silicate weathering and gypsum dissolution than on inputs originating from carbonate dissolution, as evidenced by a 14% increase in the Ca/Na ratio and 80% decrease in the SO4/HCO3 ratio, compared to the initial conditions. At the pre-flood event (pFE), two scenarios are clearly differentiated ( Figure 6). In flood event 2, the total streamflow is supplied by the baseflow as attested by the undisturbed signal (turbidity, water temperature, conductivity, and pH). The opposite scenario took place during flood event 11, where the hydrochemical and physicochemical signals were disturbed by a long period of light rainfalls, which increased the stream discharge and decreased the relative contribution of the baseflow fraction to the total streamflow (around of 50%). This can be confirmed by the increase in conductivity, Ca 2+ , Mg 2+ , HCO 3 − , and Si before the minor discharge peaks, and their subsequent gradual decrease, particularly visible during the first minor peak. As was evidenced by Walling and Foster [73] on several Devon catchments and by other authors on several calcareous basins [17,29,34,74], a mobilization of solute elements stored during the summer in the karst system and released by the flushing effect explained the variation in the conductivity and concentration of these elements (coming from the rock weathering) as discharge slightly increases.
During the rising period (RiP), up to 94% of the total streamflow is supplied by the quick-response flow ( Figure 6). The sudden increase in discharge was accompanied by an increase in turbidity ( Figure 6) and DOC (Table 3, Figure 8), whose peaks were almost simultaneous. The water entering the system, i.e., rainfall, can modify the streamwater temperature, increasing (flood 2, spring) or decreasing this signal (flood 11, autumn). This was illustrated by an upward water temperature signal during the flood event 2 and downward values during the flood event 11. Conductivity shows two antagonist patterns as discharge increases. On the one hand, a dilution effect by weakly concentrated rainwater in a period of low CO 2 production and therefore of reduced rock weathering [14,17] explain the reduction in streamwater conductivity (flood 2). On the other hand, a pistonflow behavior on the last reserves of concentrated water mobilized during the rising limb of the first major peak (10 December,22:30) accounts for the slight increase in streamwater conductivity (flood 11). However, these reserves are then depleted, and a general dilution occurs as the discharge continues to increase. This is demonstrated by a global drop in conductivity up to the highest peak flow (−9% regarding the initial value) and a decrease in concentrations of dissolved elements originating from rock weathering: Ca 2+ (∆C = −13%), Mg 2+ (−45%), HCO 3 − (−9%), SO 4 2− (−85%), and Si (−14%). The dilution by continuous heavy precipitation is more effective on inputs supplied by silicate weathering and gypsum dissolution than on inputs originating from carbonate dissolution, as evidenced by a 14% increase in the Ca/Na ratio and 80% decrease in the SO 4 /HCO 3 ratio, compared to the initial conditions.
The recession period (ReP) was separated into two sub-periods: the first one controlled by the quick-response flow (ReP-Karst), and the second one by the subsurface flow together with the baseflow (ReP-nKarst) ( Figure 6). The sub-period ReP-Karst was characterized by a stable (flood 2) or a decreasing (flood 11) conductivity signal. A dilution effect observed in flood event 11 was evidenced by a global decrease in the concentrations of all the dissolved elements. Conversely, flood event 2 showed an increase in Ca 2+ and HCO 3 − concentrations but a decrease in Mg 2+ and Si concentrations. The rapid mobilization of dissolved elements and the hydrological conditions cause a decrease in concentration in the reservoirs and hinder an immediate return to the initial conditions in the stream, in particular for source-limited elements, such as SO 4 2− (Figure 6) [29]. The sub-period ReP-nKarst was characterized by (1) an accelerated decrease in the contribution of the quick-response flow; and (2) an unchanged evolution of the concentration or a slight trend to return to their initial values. The weathering contribution of the different rocks to the streamwater chemistry also changed. A higher contribution from silicate weathering is demonstrated by a pronounced and constant decrease of the Ca/Na molar ratio (from 46 to 35) coupled to an increase in Si concentrations. A higher relative contribution from gypsum dissolution was underlined by stronger SO 4 2− concentrations and a slight increase in the SO 4 /HCO 3 molar ratio (from 0.01 to 0.03) ( Figure 6). The change in the trend of these ratios at the beginning of the sub-period ReP-nKarst thus marked the beginning of the increase in the relative contribution of lithological formations containing silicate and gypsum. Therefore, for both flood events (Figure 6), the first inputs of water and dissolved elements reaching the catchment outlet come mainly from the draining of the karst carbonate area (quick-response flow). The inputs from the other lithologies reach the outlet later on, becoming relatively more important as the quick-response flow decreases. Finally, the end of the flood event (eFE) was marked by a general trend towards a return to initial conditions ( Figure 6).

Conclusions
This paper revealed hydrological and hydrochemical behaviors in response to stream discharge variations in a multi-lithological karstic catchment located in the French Pyrenees Mountain. These changes are complex and showed different patterns regulated by the complexity of the fractured geological substratum, the magnitude of the rainfall, and the hydroclimatic conditions preceding the flood event. A high solute river load during a flood event at the beginning of the hydrological year suggest an initial piston effect of more concentrated waters in the different compartments of the catchment, which is then diluted with the increase in discharge. The fast mobilization and the antecedent hydrological conditions deplete progressively the karst and epikarst reserves of the dissolved elements and delay a quick return to the initial conditions. During the flood events, the rapid kinetics of the carbonate dissolution with biogenic soil CO 2 offsets the dilution effect by rainwater, suggesting control by a process-limited regime, as illustrated by an overall chemostatic behavior of Ca 2+ and HCO 3 − . Besides, the slower outputs and/or reduced amount sources of SO 4 2− led to a strong rainwater dilution as the discharge increases. DOC was mobilized with increasing discharge, most likely from the upper soil horizons (litter fall and organic layers), suggesting regulation by a transport-limited regime. Extreme flood events (Q max > 8.3 m 3 ·s −1 ) associated with prolonged rainfall events contributed to a dilution pattern in rock weathering products during the rising limb, even though an overall concentration increase or chemostatic response was noticed for smaller flood events (Q max < 4.4 m 3 ·s −1 ). The contribution from non-carbonate lithological compartments was shown during the second part of the recession period in both floods. The hysteresis analysis, based on two simple indicators, revealed an overall clockwise hydrochemical response, indicating higher concentrations when the major water contribution comes from the karst and the surface runoff. However, carbonate dissolution products showed an anti-clockwise response under dry hydrological conditions, suggesting a flushing of more concentrated water into the epikarst during the subsurface flow. This research points out the importance of high-frequency hydrochemical monitoring coupled with a wellintegrated discrete sampling, and of the use of a combination of hydro-physico-chemical tracers. This combined approach is a pivotal tool to better understand the dynamics of water-solute element coupling mechanisms in the critical zone of such karstic multilithological catchments.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript:     Figure 1C.
Appendix B