The Potential of Isotopic Tracers for Precise and Environmentally Clean Stream Discharge Measurements

: Accurate discharge measurement is mandatory for any hydrological study. While the “ve-locity” measurement method is adapted to laminar ﬂows, the “dilution” method is more appropriate for turbulent streams. As most low-gradient streams worldwide are neither laminar nor turbulent, a methodological gap appears. In this study, we demonstrate that the application of the “dilution” method to a low-gradient small stream gives very satisfactory results in addition to revealing surface/subsurface processes. A variety of chemical and isotopic tracers were injected into the stream (anions, ﬂuorescent dyes, and chloride and hydrogen isotopes). We report the ﬁrst use of 37 Cl for stream discharge measurement and show that 37 Cl and 2 H can be reliably used as quantitative tracers. Discharge uncertainty calculations show that deuterium is the most accurate tracer method used. We also compare the differences in the tailing part of the restitution curves of tracers and investigate the role of transient surface and hyporheic zones in solute transport in light of a simple transport modelling approach. We conclude that isotopic tracers can be used as “environmentally friendly” tracers for discrete stream discharge measurements and that the application of multi-tracers tests in rivers opens the path to a better understanding of surface–subsurface interaction processes.


Introduction
Surface waters discharge measurement is a fundamental aspect of any hydrological and hydrogeological study.Indeed, accurate and precise stream discharge measurement is critical for water resource monitoring; water budget computation; hydraulic infrastructure design, such as weirs, culverts and artificial channels; flood monitoring and forecasting; and groundwater inflows calculations [1][2][3][4].Different field approaches can be used to determine stream discharge [2,5].One of the most used techniques is based on the use of current meters.The discharge is deducted by integrating the flow velocity distribution over a river cross-section area.Discharge measurement uncertainties using current meters largely depends on the velocity measure uncertainties and is usually between 2% and 20%, depending on measurement conditions [2,6,7].This technique is best suited for open channels with near-laminar flow.Indeed, the representativeness of velocity measurements is very limited in turbulent streams [8] or in the presence of dead storage zones and backwaters [9].The second technique is the so-called "gauging by tracer dilution" method.This technique consists of injecting a known mass of tracer upstream and measuring its complete restitution over time downstream [5].Several key hypotheses must be verified: (i) complete dissolution and conservation of the tracer, (ii) homogeneity of tracer concentration across the river section at the measurement point downstream (ensured by transverse and vertical dispersion), and (iii) the discharge must be under a permanent regime during the entire experiment duration [8].The use of the tracer dilution technique is especially appropriate in shallow, turbulent streams where hydrodynamic dispersion ensures a fast tracer homogenization, i.e., in contexts not suitable for traditional current meters.
The main limitation of the "dilution" technique is the presence of transient storage zones, both at the surface and in the subsurface.Indeed, it is now widely accepted that streams exchange intimately with surrounding groundwater bodies in a variety of climates [10][11][12][13][14], especially in temperate, wet climates [14].Groundwater and surface waters act as a continuum, and exchange in the hyporheic zone, i.e., the portion of sediments surrounding the stream that is permeated with stream water, is significant [15].Transport in the hyporheic zone is thus much slower than in the main stream channel, allowing biological and geochemical processes to occur [16].In addition, the transport in the hyporheic zone, as well as in surface transient storage zones, results in a tailing effect in restitution curves of artificial tracer injection, which might limit the "dilution" method [1].Exchanges with the hyporheic zone are mainly controlled by streambed porosity and permeability.They are known to be mainly advective but can be controlled by diffusion in the presence of fine streambed sediments, such as clays, in calm streams [15].Hyporheic exchanges are typically modelled using the classical advection dispersion equation (ADE) to which a sink term representing the subsurface exchanges with the stream is added [15,[17][18][19][20].Then, hyporheic flow represents the portion of stream water that is transported through the hyporheic zone before returning to the stream.In hydrological modelling, model parameters are calibrated against field data of tracer tests, often thanks to fluorescent dyes and salts tracers.Most low-gradient small streams worldwide consist of a succession of small pools and riffles, with very contrasted flow velocities, important streambed macro-topography, and the presence of transient storage zones.The flow in these rivers is not exclusively laminar or exclusively turbulent, which indicates a gap in streamflow measurement methodology.
Globally, anthropogenic activities have various negative impacts on stream quality, mainly resulting in biodiversity loss and eutrophication issues [21,22].In addition, it is well known that groundwater pumping can decrease water levels in rivers [23].In this context, the application of the "dilution" method through injection of salts (e.g., NaCl) can be problematic for aquatic ecosystems [24] and in rivers exploited for drinking water because their use will temporarily damage the stream quality.It is thus relevant to search for more environmentally friendly tracers for stream discharge measurements.For a given river discharge, the quantity of tracer has to be minimized but will depend on the initial stream concentration and the accuracy of the measurement.Chemical elements that are naturally found in small quantities and that can be analysed with a very high accuracy can therefore be injected in sufficiently small quantities to minimize their impact on the environment.While less used, "heavy" natural stable isotopes that are naturally present in (very) small amounts can be measured with a high precision using mass spectrometry.Hence, deuterium ( 2 H) is an ideal tracer of the water molecule since it is naturally bound to it and is fully conservative.It has already been used in many artificial tracer injections, especially in various groundwaters studies [25] but not yet for river discharge measurements.Rarely used, chloride isotopes are interesting tracers because of the high solubility of chlorine and the fact that it is conservative in aquatic systems [26].Chlorine has two common stable isotopes, 35 Cl and 37 Cl, and high-accuracy mass spectrometry measurements are very promising for environmental studies.The application of these "cleaner" isotopic tracers opens new perspectives to the understanding of environmental processes.
The objectives of this study are multiple: (i) to compare both the velocity and dilution methods for discharge measurements of a small stream, (ii) to assess the potential of isotopic tracers for more environmentally friendly discharge measurements, and (iii) to document surface-subsurface processes in the light of a multi-tracer approach.To complete these objectives, a diverse array of tracers, including electrolytes (NaBr, 37 Cl-enriched salt) conservative anions (bromide [Br − ] and chloride [Cl − ]), stable isotopes (deuterium [ 2 H] and chlorine [ 37 Cl]) and fluorescent dyes (uranine, tinopal) were injected into a small stream (Ruisseau Rousse, Québec, Canada) and analysed.We compare two approaches for stream discharge determination and hypothesize that the "dilution" method can be employed with satisfactory results for low-gradient small streams.Finally, we compare the tailing effect of tinopal, 2 H, 37 Cl and Br − in restitution curves to assess the exchange processes between the stream main channel and transient storage zones.The simultaneous injection of tracers having different chemical and physical (size) properties enabled us to gain a better understanding of these interactions/processes.measurements are very promising for environmental studies.The application of these "cleaner" isotopic tracers opens new perspectives to the understanding of environmental processes.

Description of the Field Site
The objectives of this study are multiple: (i) to compare both the velocity and dilution methods for discharge measurements of a small stream, (ii) to assess the potential of isotopic tracers for more environmentally friendly discharge measurements, and (iii) to document surface-subsurface processes in the light of a multi-tracer approach.To complete these objectives, a diverse array of tracers, including electrolytes (NaBr, 37 Cl-enriched salt) conservative anions (bromide [Br − ] and chloride [Cl − ]), stable isotopes (deuterium [ 2 H] and chlorine [ 37 Cl]) and fluorescent dyes (uranine, tinopal) were injected into a small stream (Ruisseau Rousse, Québec, Canada) and analysed.We compare two approaches for stream discharge determination and hypothesize that the "dilution" method can be employed with satisfactory results for low-gradient small streams.Finally, we compare the tailing effect of tinopal, 2 H, 37 Cl and Br − in restitution curves to assess the exchange processes between the stream main channel and transient storage zones.The simultaneous injection of tracers having different chemical and physical (size) properties enabled us to gain a better understanding of these interactions/processes.

Description of the Field Site
Field work occurred in Oka, Québec, Canada in August 2022 in a small stream under baseflow conditions (coordinates: 45°30′15.8″N/74°02′45.9″W-see Figure 1 below) on a sunny day without any precipitation.In addition, the last precipitation event occurred 8 days before the test, roughly 4 times the recession time of this small river.The basement geology of the watershed consists of intrusive plutonic rocks (charnockite, mangerite) from the Precambrian age and Mesozoic carbonatite.It is covered by quaternary glaciomarine deposits (silts to gravel sands) that shelter the alluvial aquifer sustaining the river [27].The local climate is continental humid, characterized by long, cold, and snowy winters and hot and wet summers.The watershed is mainly cultivated, but its southern part is occupied by the Parc National d'Oka, a protected natural area.The stream, called Ruisseau Rousse, drains the watershed and flows to the Lac des Deux Montagnes before ultimately joining the St. Lawrence River.The studied river transect extends 105 metres from the injection site (IS) to the restitution site (RS) (Figure 2).Similar to numerous streams in intermediate gradient plains, it The studied river transect extends 105 metres from the injection site (IS) to the restitution site (RS) (Figure 2).Similar to numerous streams in intermediate gradient plains, it is characterized by a succession of pools and small riffles with very contrasted flow velocities, river widths, and river depths.The average river width was estimated to be 2 m (2.6 m at the studied cross-section), with the maximal cross-sectional depth varying between 10 and 60 cm.
is characterized by a succession of pools and small riffles with very contrasted flow velocities, river widths, and river depths.The average river width was estimated to be 2 m (2.6 m at the studied cross-section), with the maximal cross-sectional depth varying between 10 and 60 cm.

Field Work
The solution used for the tracing experiment contained approximately 20 L of river water to which a variety of tracers was added: 535.13 ± 0.01 g of Cl − with a δ 37 Cl value of + 2.18‰ ± 0.03 vs. SMOC (supplier: BioBasic, Markham, ON, Canada), 788.60 ± 0.01 g of Br − (supplier: Fischer Scientific, Ottawa, ON, Canada), 403.09 ± 0.20 g of 2 H as deuterated water 2 H2O (supplier: Isowater, Collingwood, ON, Canada), 2.42 ± 0.01 g of uranine, and 16.50 ± 0.01 g of tinopal (supplier: ThermoFischer, Waltham, MA, USA).The sample was placed in a sealed bucket until the salts dissolved completely, which took a few minutes.The total dissolved solids (TDS) of the mixture was 126.4 g/L.It was then injected instantaneously at IS.We ensured tracer homogeneity across the river width by manually fanning the river at IS shortly after injection (Figure 2).The injected volume has no effect regarding the total stream discharge.
Stream water velocity was measured using a SonTek FlowTracker at RS in 55 points across the river section to assess its distribution with a high spatial resolution, minutes before the injection occurred.At RS, a total of five field multimeters (WTW) and 2 Solinst ® LTC loggers were used to monitor the distribution of salt tracers (using electrical conductivity) across the section at a time step of 5 s.All of the equipment was calibrated against a 1413 µS/cm standard in the morning of the field experiment.The mean uncertainty in electrical conductivity measurements was reported to be 1.1% of the measured value.Water fluorescence was monitored at RS in the centre of the stream (green dot in Figure 3), at a time step of 15 s using a field fluorometer (Tetraedre GGUN FL-30).Samples for anions (Cl − , Br − ) and isotopic compositions ( 2 H, 37 Cl) were taken in the centre of stream at two points (Fluo.and WTW3) at 10 and 30 cm of depth (distance = 120 cm from the left bank, see Figure 3) in 250 mL HDPE bottles.The first discrete samples were taken every 2 min, and the time interval of sampling was progressively reduced up to 30 s when approaching to the arrival of the peak of tracer.It was then increased up to 5 min during the tailing part.

Field Work
The solution used for the tracing experiment contained approximately 20 L of river water to which a variety of tracers was added: 535.13 ± 0.01 g of Cl − with a δ 37 Cl value of +2.18‰ ± 0.03 vs. SMOC (supplier: BioBasic, Markham, ON, Canada), 788.60 ± 0.01 g of Br − (supplier: Fischer Scientific, Ottawa, ON, Canada), 403.09 ± 0.20 g of 2 H as deuterated water 2 H 2 O (supplier: Isowater, Collingwood, ON, Canada), 2.42 ± 0.01 g of uranine, and 16.50 ± 0.01 g of tinopal (supplier: ThermoFischer, Waltham, MA, USA).The sample was placed in a sealed bucket until the salts dissolved completely, which took a few minutes.The total dissolved solids (TDS) of the mixture was 126.4 g/L.It was then injected instantaneously at IS.We ensured tracer homogeneity across the river width by manually fanning the river at IS shortly after injection (Figure 2).The injected volume has no effect regarding the total stream discharge.
Stream water velocity was measured using a SonTek FlowTracker at RS in 55 points across the river section to assess its distribution with a high spatial resolution, minutes before the injection occurred.At RS, a total of five field multimeters (WTW) and 2 Solinst ® LTC loggers were used to monitor the distribution of salt tracers (using electrical conductivity) across the section at a time step of 5 s.All of the equipment was calibrated against a 1413 µS/cm standard in the morning of the field experiment.The mean uncertainty in electrical conductivity measurements was reported to be 1.1% of the measured value.Water fluorescence was monitored at RS in the centre of the stream (green dot in Figure 3), at a time step of 15 s using a field fluorometer (Tetraedre GGUN FL-30).Samples for anions (Cl − , Br − ) and isotopic compositions ( 2 H, 37 Cl) were taken in the centre of stream at two points (Fluo.and WTW3) at 10 and 30 cm of depth (distance = 120 cm from the left bank, see Figure 3) in 250 mL HDPE bottles.The first discrete samples were taken every 2 min, and the time interval of sampling was progressively reduced up to 30 s when approaching to the arrival of the peak of tracer.It was then increased up to 5 min during the tailing part.

Analytical Methods
Since discharge computation requires concentration data, the electrical conductivity signal was converted into an added salts concentration time series.To do so, a sample of the river water was obtained prior to the injection (1L in HDPE bottle).Salts were carefully added in the same proportion as the injection mix (i.e., 40.4% of NaBr and 59.6% of 37 Cl-enriched salt) and progressively dissolved, and the relationship between electrical conductivity and added dissolved salts was assessed for each field probe.The field fluorometer was calibrated at the laboratory on the day before the tracing experiment using distilled water to which uranine and tinopal were added in known quantities (7 points ranging from 0 to 160 ppb for uranine and 10 points ranging from 0 to 1300 ppb for tinopal).Anion (chloride, bromide) analyses were performed using liquid ion chromatography at the Université du Québecà Montréal, QC, Canada (Geotop-UQAM lab).

Analytical Methods
Since discharge computation requires concentration data, the electrical conductivity signal was converted into an added salts concentration time series.To do so, a sample of the river water was obtained prior to the injection (1L in HDPE bottle).Salts were carefully added in the same proportion as the injection mix (i.e., 40.4% of NaBr and 59.6% of 37 Clenriched salt) and progressively dissolved, and the relationship between electrical conductivity and added dissolved salts was assessed for each field probe.The field fluorometer was calibrated at the laboratory on the day before the tracing experiment using distilled water to which uranine and tinopal were added in known quantities (7 points ranging from 0 to 160 ppb for uranine and 10 points ranging from 0 to 1300 ppb for tinopal).Anion (chloride, bromide) analyses were performed using liquid ion chromatography at the Université du Québecà Montréal, QC, Canada (Geotop-UQAM lab).
The raw field data for all the tracers can be found as Supplementary Material (Excel spreadsheet S1).

Data Interpretation
Traditionally, when using a current meter, the discharge is deducted by integrating the flow velocity distribution over the whole river cross-section area: where Q is the stream discharge [L 3 T −1 ], and v the water velocity field [LT −1 ] that is a function of width (x-[L]) and depth (y-[L]).The variables x and y are used to define the river cross-section area, A.
Using the "dilution" method requires verification of the following hypotheses: (i) complete dissolution and conservation of the tracer, (ii) homogeneity of tracer concentration across the river section (ensured by transverse and vertical dispersion), and (iii) stream discharge under a permanent regime during the entire experiment duration [8].When such conditions are met, the discharge is computed according to the following formula [1,25]: where M tracer is the mass of tracer that is injected [M], and C(t) is the added tracer concentration at time t [ML −3 ].The time t = 0 corresponds to the time of injection, while T is the total duration of the experiment.Discharge measurements using the "dilution" method can only be overestimated if one or more of the three assumptions are not met.Usually, the denominator is calculated using a simple numerical integration scheme, such as that noted in the following: where i refers to the i th concentration measurement of the tracer (at time t i ), and N indicates the total number of measurements.Practically, the first sample (C 1 ) is taken at time 0, and the N th sample (C N ) is taken at time T. Typically, when using field probes, the measurement timestep (∆T) is constant, and Equation ( 5) is represented as follows: In addition, the uncertainty propagation calculation was performed based on Equations ( 4) and ( 5): where ∆ refers to the measure uncertainty (1σ) in any given measurement.This assumes that there is no uncertainty in the timing of the measurement.Of note, in the case of a constant sampling timestep ∆T (e.g., electrical conductivity monitoring), Equation (7) simplifies to the following: Tracer experiments allow computation of stream discharge, but they also allow estimation of flow parameters, such as dispersion, velocity, and the influence of slower velocity zones.Here, flow parameters were determined by modelling the tracers' restitution curves using a simple binary 1-D dispersion mixing model.A similar lumped parameter model approach was used to simulate the transport of an artificial tracer release in the Rhine river and is described in Leibundgut et al. [25].Considering a conservative tracer, the classical 1-D lumped dispersion model can be expressed as follows [31]: Here, Here, C is the tracer concentration over time [ML −3 ], D is the longitudinal dispersion parameter [L 2 T −1 ], x the distance between the injection and measurement points [L], v the mean stream velocity [LT −1 ], and C inj is the injection function [ML −3 ].Here, it is a Dirac-type function with a value of C 0 at t = 0.
In this binary mixing model approach, the measured tracer concentration was simulated by adding a second flow component representing the slower fraction: The model was distributed over time and implemented in Excel (time step of 10 s).It was fitted with tracer restitution curves by minimising the normalized root mean square deviation defined as (where C meas is the measured value, C sim is the simulated value).The parameters of the model were fitted using the restitution curves of 2 H, 37 Cl, and Br − .The Excel spreadsheet containing the model as well as its documentation can be downloaded as Supplementary Materials (Excel spreadsheet S2 and Text file S1 respectively).

River Section Water Velocity Distribution
The survey performed using the FlowTracker allowed visualization of the water velocity distribution at the cross-section scale.The FlowTracker method yielded a discharge of 50.2 ± 3.6 L/s (Table 1).Zones of apparent stagnant water are identified near the banks (Figure 4), while most of the flow is in the centre of the stream at depths between 20 and 40 cm.Such discrepancies in the river velocities are expected to play a significant role in solute transport and on our tracing experiment.Their role on the interpretation of the tracer injection restitutions is discussed below.

Using Salt Tracers to Build a Representative Restitution Curve for the Whole Section
A total of 7 restitution curves for electrical conductivity were obtained in the field.Restitution curves-expressed in mg/L of added salt tracers-are shown in Figure 5.The peak time is defined for each restitution curve as the time elapsed between the injection time and when the maximum concentration was measured.Tracer homogeneity across the river width was ensured at the injection point by laterally mixing the river water seconds after the injection.Except for WTW5, all restitution curves are very similar (Figure 5) in terms of peak time.This indicates that the transverse dispersion occurring along the reach is sufficient to ensure a good tracer homogeneity.The probe WTW5 was installed near the left bank in a zone of very low velocity with backwaters (Figure 4).Its peak height is lower than others and delayed with an increased peak width, which indicates tracer exchange between the main channel and transient storage zones [32,33].

Using Salt Tracers to Build a Representative Restitution Curve for the Whole Section
A total of 7 restitution curves for electrical conductivity were obtained in the field.Restitution curves-expressed in mg/L of added salt tracers-are shown in Figure 5.The peak time is defined for each restitution curve as the time elapsed between the injection time and when the maximum concentration was measured.Tracer homogeneity across the river width was ensured at the injection point by laterally mixing the river water seconds after the injection.Except for WTW5, all restitution curves are very similar (Figure 5) in terms of peak time.This indicates that the transverse dispersion occurring along the reach is sufficient to ensure a good tracer homogeneity.The probe WTW5 was installed near the left bank in a zone of very low velocity with backwaters (Figure 4).Its peak height is lower than others and delayed with an increased peak width, which indicates tracer exchange between the main channel and transient storage zones [32,33].Equations ( 4) and ( 5) were used to compute discharge values for each curve.In addition, the uncertainty propagation calculation was performed using the formula given in Equation (7).Discharge values are between 49.8 ± 1.1 and 56.9 ± 1.3 L/s.Most of them agree with the value obtained with the FlowTracker.However, the values obtained at measuring points WTW1, WTW2, and WTW3 are slightly above the value measured with the FlowTracker (Table 1).We suspect that there could be either a slight underestimation of the flow velocity measurements using the FlowTracker or an overestimation of the quantity of tracer.However, the latter is unlikely as it would shift all the obtained values and the quantities of tracer used were also thoroughly measured.
In all strictness, the second hypothesis of the use of Equation (4) (i.e., homogeneity of tracer concentration across the river section) is not entirely verified.We show below that Equation (4) can still be applied, with minimal effect on discharge calculations.To do so, we calculated a tracer output function averaged for the whole cross-section.This new restitution function (called "composite" in the following) includes the whole mass of tracer and is calculated as follows.Each restitution curve is supposed to be representative of the local tracer output.The river section was cut into seven sub-areas following Thiessen tillage (Figure 4), each of which were represented by one conductivity probe.The interpolated water velocity distribution presented in Figure 4 was then used to compute the discharge related to each sub-area.The latter is computed by multiplying the interpolated mean flow velocity by the surface of the sub-area.Typically, sub-areas near the banks have Equations ( 4) and ( 5) were used to compute discharge values for each curve.In addition, the uncertainty propagation calculation was performed using the formula given in Equation (7).Discharge values are between 49.8 ± 1.1 and 56.9 ± 1.3 L/s.Most of them agree with the value obtained with the FlowTracker.However, the values obtained at measuring points WTW1, WTW2, and WTW3 are slightly above the value measured with the FlowTracker (Table 1).We suspect that there could be either a slight underestimation of the flow velocity measurements using the FlowTracker or an overestimation of the quantity of tracer.However, the latter is unlikely as it would shift all the obtained values and the quantities of tracer used were also thoroughly measured.
In all strictness, the second hypothesis of the use of Equation (4) (i.e., homogeneity of tracer concentration across the river section) is not entirely verified.We show below that Equation (4) can still be applied, with minimal effect on discharge calculations.To do so, we calculated a tracer output function averaged for the whole cross-section.This new restitution function (called "composite" in the following) includes the whole mass of tracer and is calculated as follows.Each restitution curve is supposed to be representative of the local tracer output.The river section was cut into seven sub-areas following Thiessen tillage (Figure 4), each of which were represented by one conductivity probe.The interpolated water velocity distribution presented in Figure 4 was then used to compute the discharge related to each sub-area.The latter is computed by multiplying the interpolated mean flow velocity by the surface of the sub-area.Typically, sub-areas near the banks have a lower discharge than the sub-areas located in the middle of the flow.The composite restitution function (see Figure 5) was then calculated as follows: where i refers to the sub-area related to each probe (7 being the total number), C is the salt concentration [ML −3 ] and f i is the fraction of total discharge passing by the sub-area i (adimensional).This cross-sectional flow-averaged function is very close to restitution functions obtained from the main flow channel since the weight of the dead zones is very limited.Indeed, the discharge and uncertainty propagation calculations yielded Q composite = 53.3± 6.5 L/s, which is very close to discharges calculated over the crosssection.This strongly suggests that the "dilution" method is adapted to such small lowgradient streams with contrasted velocities.Furthermore, the calculation proposed in this work (Equation ( 12)) is adapted to any stream, especially if restitution functions are much different from the banks to the middle and requires (i) a precise survey of streamflow velocities across the whole section and (ii) the monitoring of tracers at various spots across the river section.The advantage of this function is that it verifies all the hypotheses of Equation ( 3) since it simulates the section-integrated homogenised tracer concentration dynamics.However, its main limitation is obviously the uncertainty propagation.While individual restitutions yield discharges values with uncertainties between 2.2% and 3.2%, the composite restitution gives uncertainties of 12.2%.

Isotopic Tracers as New Tracers for Discrete and Clean Discharge Measurements-Comparison with Other Tracers
Salt tracers are very useful for stream discharge measurements because they can easily be found, are inexpensive and are easy to monitor in the field thanks to electrical conductivity probes.However, their use can be very problematic in watersheds sheltering endangered ecosystems and species or used for surface drinking water production.In such contexts, it is worth to search for more adapted artificial tracers that can be injected in sufficiently low quantities that do not harm the environment and do not draw public attention like fluorescent dyes.
Stable isotopes tracers are now widely used in hydrology and hydrogeology as natural tracers [25,26,34,35].Their use as artificial tracers is also reported in the literature, especially in the unsaturated zone and for tracing groundwater flows over short distances (see [25] and references therein).However, their potential for discrete, clean, and precise stream discharge measurements has not been addressed to date.In this work, heavy stable isotopes 2 H (deuterated water, 2 H 2 O) and 37 Cl (isotopically enriched salt) were used as conservative tracers for quantitative interpretation of river discharge and comparison with other tracers.
The injected mass of each isotopic tracer was determined using the known mass of injected deuterated water and salt and their isotopic contents.Indeed, the deuterated water has a 2 H abundancy of 99.90 ± 0.05%, and the 37 Cl-enriched salt has a δ 37 Cl value of +2.18 ± 0.03‰ vs. SMOC, yielding totals of 403.09 ± 0.20 g of 2 H and 131.17 ± 0.09 g of 37 Cl.Stream discharge was computed using Equations ( 4) and ( 5), and uncertainty propagation was performed using Equation (7) (details for isotopes are also presented in Annex).Results are presented in Table 1 and the restitution curves of tracers can be found in Figure 6.
The highest discharge is obtained for uranine (60.6 ± 0.7 L/s).Only 2.42 g of uranine was injected, and a degradation of roughly 0.3 g could explain the shift obtained in discharge values.It is unlikely to have been considerably affected by photodegradation as the experiment time is less than 2 h, but it could be explained by an underestimation of the injected mass linked to the quality of the tracer and the fact that the calibration of the fluorometer was performed prior the field test using distilled water.Tinopal yielded a discharge value of 53.0 ± 0.5 L/s, which is in the same range as the FlowTracker measurement and electrical conductivity monitoring values.The discharge value obtained with deuterium is 51.1 ± 0.4 L/s, while 37 Cl data yielded a discharge of 39.4 ± 2.3 L/s.One must be aware that δ 37 Cl measures uncertainties that are usually of the order of ±0.1‰, which is greater than the ±0.03‰ used in our calculations.These higher uncertainties would yield a total uncertainty of ±9.0 L/s (±22.9% of the discharge value).However, even though the discharge value obtained with 37 Cl is consistent with the discharge obtained using Cl − , it is very different from other tracers.The chloride and 37 Cl signals seem more concentrated relatively to other tracers' signals.The highest discharge is obtained for uranine (60.6 ± 0.7 L/s).Only 2.42 g of uranine was injected, and a degradation of roughly 0.3 g could explain the shift obtained in discharge values.It is unlikely to have been considerably affected by photodegradation as the experiment time is less than 2 h, but it could be explained by an underestimation of the injected mass linked to the quality of the tracer and the fact that the calibration of the fluorometer was performed prior the field test using distilled water.Tinopal yielded a discharge value of 53.0 ± 0.5 L/s, which is in the same range as the FlowTracker measurement and electrical conductivity monitoring values.The discharge value obtained with deuterium is 51.1 ± 0.4 L/s, while 37 Cl data yielded a discharge of 39.4 ± 2.3 L/s.One must be aware that δ 37 Cl measures uncertainties that are usually of the order of ± 0.1‰, which is greater than the ±0.03‰ used in our calculations.These higher uncertainties would yield a total uncertainty of ± 9.0 L/s (±22.9% of the discharge value).However, even though the discharge value obtained with 37 Cl is consistent with the discharge obtained using Cl − , it is very different from other tracers.The chloride and 37 Cl signals seem more concentrated relatively to other tracers' signals.
Strictly speaking, for the same mass of a tracer, the use of 2 H as an artificial tracer allows for a more precise discharge measure than the use of 37 Cl (Figure 7).The calculations used for Figure 7 are presented in the Appendix.However, even though deuterium analyses are widely affordable, the purchase of deuterated water is subject to governmental limitations, so one must be aware of potential difficulties in finding such a tracer in addition to the need for in-house produced enriched standards.On the other hand, while Strictly speaking, for the same mass of a tracer, the use of 2 H as an artificial tracer allows for a more precise discharge measure than the use of 37 Cl (Figure 7).The calculations used for Figure 7 are presented in the Appendix.However, even though deuterium analyses are widely affordable, the purchase of deuterated water is subject to governmental limitations, so one must be aware of potential difficulties in finding such a tracer in addition to the need for in-house produced enriched standards.On the other hand, while chloride is simple to find and cheap, 37 Cl-enriched salts are much harder to find.In addition, one must be aware that the background noise of 2 H in rivers is commonly of the order of 140 to 156 ppm (for a δ 2 H ranging from −100‰ to 0‰ vs. VSMOW) and that the 37 Cl background noise of rivers with high chloride concentrations might be high because of the relatively high natural abundance of 37 Cl ( 37 R SMOC = 0.324).In the quantities that were injected, deuterated water can be considered as the best tracer to measure stream discharge because it does not modify or alter the qualitative and ecological state of the stream while being a robust and extremely precise tracer.
discharge because it does not modify or alter the qualitative and ecological state of stream while being a robust and extremely precise tracer.  3Cl) fo studied river (the blue dots represent this experiment).It is important to note that the x-axis r to the actual isotopic tracer's mass alone and not 37 Cl-enriched salt or deuterated water.Detai these calculations can be found in Appendix A.

Transient Storage Zone Exchanges Inferred from Multi-Tracer Data
When a tracer is injected into the stream, advective-dispersive transport will resu solute tracer homogenization over the river width due to transverse dispersion and a c sical Gaussian form of the restitution function due to longitudinal dispersion.Howe the tracer will also exchange with transient storage zones, such as dead zones near banks and/or the hyporheic zone at the bottom of the stream, and will be slowly relea to the main flow channel [17,18,20].Transport in the hyporheic zone is known for b mainly advective [15].Since tracer residence times in the transient storage zones are t cally several orders of magnitude higher than those of the tracer in the stream alone formation on hyporheic exchanges can be found on the tailing part of the breakthro curves.
Here, we took advantage of having multiple tracer restitutions to study the inte tions between the stream and the hyporheic zone in more detail.Normalized br through curves (C/Cmax) of tinopal, 2 H, 37 Cl, Br − , and Cl − show different behaviours du the tailing part of the experiment (Figure 8).If we define persistence as the affinity of tracer to stay in the stream for long times, bromide is the least persistent tracer, w chloride (and especially 37 Cl) is the most persistent. 2 H lies in between, and being pa the water molecule, is the ideal tracer of the latter.It seems that the surface solute trans is affected by one (or more) hydrological unit(s) that acts as a filter, filtering the entra of the tracer and thus controlling the residence time of tracer within it.Such a hydrolog unit could be a surface transient storage zone and/or hyporheic zone.  3Cl) for the studied river (the blue dots represent this experiment).It is important to note that the x-axis refers to the actual isotopic tracer's mass alone and not 37 Cl-enriched salt or deuterated water.Details on these calculations can be found in Appendix A.

Transient Storage Zone Exchanges Inferred from Multi-Tracer Data
When a tracer is injected into the stream, advective-dispersive transport will result in solute tracer homogenization over the river width due to transverse dispersion and a classical Gaussian form of the restitution function due to longitudinal dispersion.However, the tracer will also exchange with transient storage zones, such as dead zones near the banks and/or the hyporheic zone at the bottom of the stream, and will be slowly released to the main flow channel [17,18,20].Transport in the hyporheic zone is known for being mainly advective [15].Since tracer residence times in the transient storage zones are typically several orders of magnitude higher than those of the tracer in the stream alone, information on hyporheic exchanges can be found on the tailing part of the breakthrough curves.
Here, we took advantage of having multiple tracer restitutions to study the interactions between the stream and the hyporheic zone in more detail.Normalized breakthrough curves (C/C max ) of tinopal, 2 H, 37 Cl, Br − , and Cl − show different behaviours during the tailing part of the experiment (Figure 8).If we define persistence as the affinity of the tracer to stay in the stream for long times, bromide is the least persistent tracer, while chloride (and especially 37 Cl) is the most persistent. 2 H lies in between, and being part of the water molecule, is the ideal tracer of the latter.It seems that the surface solute transport is affected by one (or more) hydrological unit(s) that acts as a filter, filtering the entrance of the tracer and thus controlling the residence time of tracer within it.Such a hydrological unit could be a surface transient storage zone and/or hyporheic zone.
Since the tracers used in this study have different chemical and physical properties, we can expect them to behave differently regarding hyporheic and surface dead zone exchanges.Exchanges between the stream and the subsurface (e.g., the hyporheic zone) are controlled by the streambed permeability and porosity [15] but might also be controlled by the effective accessible pore size of tracers.Bromide restitution shows a lower tailing effect, meaning that it is less affected by such exchanges than 2 H and 37 Cl.The breakthrough curves of the latter two indicates that the release of chloride is slower than that of deuterium.This makes sense, as chloride transport in a porous media might be slowed due to its size compared to 2 H. Tinopal was expected to fall below the Br − curve since it is the largest molecule used (chemical formula: C 40 H 40 N 12 O 8 S 2 ).Because of its high susceptibility to be absorbed to suspended particles, such as clays and organic matter, and its low effective solubility [36,37], tinopal is often not considered a totally conservative tracer.This results in a prolonged tailing effect, that might explain why it is more "persistent" than bromide.Since the tracers used in this study have different chemical and physical properties, we can expect them to behave differently regarding hyporheic and surface dead zone exchanges.Exchanges between the stream and the subsurface (e.g., the hyporheic zone) are controlled by the streambed permeability and porosity [15] but might also be controlled by the effective accessible pore size of tracers.Bromide restitution shows a lower tailing effect, meaning that it is less affected by such exchanges than 2 H and 37 Cl.The breakthrough curves of the latter two indicates that the release of chloride is slower than that of deuterium.This makes sense, as chloride transport in a porous media might be slowed due to its size compared to 2 H. Tinopal was expected to fall below the Br − curve since it is the largest molecule used (chemical formula: C40H40N12O8S2).Because of its high susceptibility to be absorbed to suspended particles, such as clays and organic matter, and its low effective solubility [36,37], tinopal is often not considered a totally conservative tracer.This results in a prolonged tailing effect, that might explain why it is more "persistent" than bromide.
Bromide, 37 Cl, and 2 H restitution curves were fitted to the binary mixing model described in Equations ( 9)-(11) (Figure 9).The f1 parameter that represents the fraction of the fast component as well as its velocity, v1, were set as the same for all tracers because this part of the flow is controlled by advective transport and should not vary from one tracer to the other.However, parameters for the slow component (dispersion and velocity) were set to vary to model the differences in tailing seen in Figure 8 and account for the different behaviours of tracers.The set of parameters that provided the best fit is given in the Table 2.The transient zone storage fraction, f2 (f2 = 1 − f1), is 0.42.This means that it exerts a very strong influence in the transport characteristics of the dissolved solutes.It is thought to be the fraction of the flow that is affected by both surface transient storage and transport within the hyporheic zone.Bromide and deuterium breakthrough curves are well modelled, especially at longer times.However, 37 Cl is affected by another phenomenon.Indeed, we saw earlier that the discharge measured using both chloride and its heavy isotope was somewhat less than that noted with other tracers.In addition, the model predicts a higher concentration than what was measured in the field.It is still unclear why chloride behaves in such a way.Bromide, 37 Cl, and 2 H restitution curves were fitted to the binary mixing model described in Equations ( 9)-(11) (Figure 9).The f 1 parameter that represents the fraction of the fast component as well as its velocity, v 1 , were set as the same for all tracers because this part of the flow is controlled by advective transport and should not vary from one tracer to the other.However, parameters for the slow component (dispersion and velocity) were set to vary to model the differences in tailing seen in Figure 8 and account for the different behaviours of tracers.The set of parameters that provided the best fit is given in the Table 2.The transient zone storage fraction, f 2 (f 2 = 1 − f 1 ), is 0.42.This means that it exerts a very strong influence in the transport characteristics of the dissolved solutes.It is thought to be the fraction of the flow that is affected by both surface transient storage and transport within the hyporheic zone.Bromide and deuterium breakthrough curves are well modelled, especially at longer times.However, 37 Cl is affected by another phenomenon.Indeed, we saw earlier that the discharge measured using both chloride and its heavy isotope was somewhat less than that noted with other tracers.In addition, the model predicts a higher concentration than what was measured in the field.It is still unclear why chloride behaves in such a way.

Conclusions
Solute transport in a small low-gradient stream was studied under baseflo tions.The stream is characterized by contrasting channel characteristics with a s of small riffles and pools.The cross-sectional dynamics of salt tracers (monito the electrical conductivity signal) revealed a very good homogenization of the tr cept in a large dead zone near one of the banks.The computation of a compos restitution, integrating the cross-sectional dynamics, showed that the "dilution is still applicable in such small, low-gradient, heterogeneous streams.Stable iso and 37 Cl are reliable tracers for discharge measurements in small streams.In

Conclusions
Solute transport in a small low-gradient stream was studied under baseflow conditions.The stream is characterized by contrasting channel characteristics with a succession of small riffles and pools.The cross-sectional dynamics of salt tracers (monitored using the electrical conductivity signal) revealed a very good homogenization of the tracers, except in a large dead zone near one of the banks.The computation of a composite tracer restitution, integrating the cross-sectional dynamics, showed that the "dilution" method is still applicable in such small, low-gradient, heterogeneous streams.Stable isotopes 2 H and 37 Cl are reliable tracers for discharge measurements in small streams.In addition,

Field
work occurred in Oka, Québec, Canada in August 2022 in a small stream under baseflow conditions (coordinates: 45 • 30 15.8 N/74 • 02 45.9 W-see Figure 1 below) on a sunny day without any precipitation.In addition, the last precipitation event occurred 8 days before the test, roughly 4 times the recession time of this small river.The basement geology of the watershed consists of intrusive plutonic rocks (charnockite, mangerite) from the Precambrian age and Mesozoic carbonatite.It is covered by quaternary glacio-marine deposits (silts to gravel sands) that shelter the alluvial aquifer sustaining the river [27].The local climate is continental humid, characterized by long, cold, and snowy winters and hot and wet summers.The watershed is mainly cultivated, but its southern part is occupied by the Parc National d'Oka, a protected natural area.The stream, called Ruisseau Rousse, drains the watershed and flows to the Lac des Deux Montagnes before ultimately joining the St. Lawrence River.

Figure 1 .
Figure 1.Map of Canada (on the left) and of the Ruisseau Rousse watershed (on the right).The red dots in both maps symbolise where the field work occurred.

Figure 1 .
Figure 1.Map of Canada (on the left) and of the Ruisseau Rousse watershed (on the right).The red dots in both maps symbolise where the field work occurred.

Figure 2 .
Figure 2. (On the left): picture of the stream at the injection site (under a small culvert).We ensured tracer homogeneity by manually fanning the water under the culvert.(On the right): picture of the restitution site with the equipment used to monitor solute concentration.

Figure 2 .
Figure 2. (On the left): picture of the stream at the injection site (under a small culvert).We ensured tracer homogeneity by manually fanning the water under the culvert.(On the right): picture of the restitution site with the equipment used to monitor solute concentration.

Figure 3 .
Figure 3. River working section (looking upstream) with equipment location.

Figure 4 .
Figure 4. Distribution of interpolated river water velocity (n = 342).A value of zero was forced at the streambed and above the water surface for the spatial interpolation.Velocities are expressed in m.s −1 .

Figure 4 .
Figure 4. Distribution of interpolated river water velocity (n = 342).A value of zero was forced at the streambed and above the water surface for the spatial interpolation.Velocities are expressed in m.s −1 .

Hydrology 2024 , 19 Figure 5 .
Figure 5. Upper graph: the seven salts concentration restitution functions and the composite restitution function.Lower graph: zoom-in around the peak time showing slight differences in peak heights.Background initial concentrations were removed, which explains the zero concentration at the beginning of the experiment.

Figure 5 .
Figure 5. Upper graph: the seven salts concentration restitution functions and the composite restitution function.Lower graph: zoom-in around the peak time showing slight differences in peak heights.Background initial concentrations were removed, which explains the zero concentration at the beginning of the experiment.

Figure 7 .
Figure 7. Discharge measure uncertainty as a function of mass of injected tracer ( 2 H,37 Cl) fo studied river (the blue dots represent this experiment).It is important to note that the x-axis r to the actual isotopic tracer's mass alone and not37 Cl-enriched salt or deuterated water.Detai these calculations can be found in Appendix A.

Figure 7 .
Figure 7. Discharge measure uncertainty as a function of mass of injected tracer ( 2 H,37 Cl) for the studied river (the blue dots represent this experiment).It is important to note that the x-axis refers to the actual isotopic tracer's mass alone and not37 Cl-enriched salt or deuterated water.Details on these calculations can be found in Appendix A.

Table 1 .
Discharge values measured and associated uncertainties.Note that anion (Cl − , Br − ) uncertainties were not available.

Table 2 .
Best-fit set of parameters for the modelling of solute transport affected by diffusive transfer in the hyporheic zone.