Phenomenological model of suspended sediment
transport in a small catchment

Amande Roque-Bernard, Antoine Lucas, Eric Gayer, Pascal Allemand, Céline

Dessert, Eric Lajeunesse

► To cite this version:

Amande Roque-Bernard, Antoine Lucas, Eric Gayer, Pascal Allemand, Céline Dessert, et al. Phenomenological model of suspended sediment transport in a small catchment. Earth Surface Dynamics, 2023, 11 (3), pp.363-381. 10.5194/esurf-11-363-2023. hal-04089800

HAL Id: hal-04089800

https://hal.science/hal-04089800

Submitted on 5 May 2023

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

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

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

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

# Phenomenological model of suspended sediment transport in a small catchment

Amande Roque-Bernard¹, Antoine Lucas¹, Eric Gayer¹, Pascal Allemand², Céline Dessert¹, and
Eric Lajeunesse¹

¹Institut de physique du globe de Paris, Université Paris Cité, CNRS, 75005, Paris, France

²Laboratoire de Géologie de Lyon, Terre Planètes Environnement, Université de Lyon,
Université Lyon 1 & ENS Lyon & CNRS, UMR 5276, 69100 Villeurbanne, France

**Correspondence:** Amande Roque-Bernard (a.roqber@gmail.com) and Antoine Lucas (lucas@ipgp.fr)

Received: 7 September 2022 – Discussion started: 22 September 2022

Revised: 16 March 2023 – Accepted: 7 April 2023 – Published: 5 May 2023

**Abstract.** We develop a phenomenological model of suspended sediment transport on the basis of data acquired in the Capesterre river, which drains a small tropical catchment in Guadeloupe. The model correctly represents the concentration of suspended sediment during floods, provided that the relation between concentration and water level forms a counterclockwise loop. In the model, the properties of the sediment and of the river are all lumped into four parameters: a settling velocity related to the size of the suspended sediment, a threshold water level which acts as a proxy for the threshold shear stress, a characteristic entrainment rate, and a dimensionless exponent. The value of the parameters changes from one flood to the next, probably reflecting changes in the characteristics of the river and the fine sediment. Finally, a test of the model against data acquired in a small catchment in the French Alps suggests that the model is versatile enough to be used in diverse hydrological settings.

## 1 Introduction

Rivers transport sediments according to their particle size. Coarse particles, like gravel or pebbles, bounce and roll above the sediment bed (Bagnold, 1973). Finer sediments, like silt, clay, or even sand, behave differently: their falling velocity is comparable to the fluctuations of the flow velocity induced by turbulence, and they remain suspended in the water column, trapped by turbulent eddies (Van Rijn, 1984). They thus advance with the flow, until they eventually settle back on the riverbed (Phillips et al., 2019). This suspended load is often the main contribution to the sediment discharge that a river carries out of its watershed (Turowski et al., 2010; Liu et al., 2011). It is therefore a major component of the erosion of continental surfaces (Summerfield and Hulton, 1994; Syvitski et al., 2003). Yet the matter suspended in a river does not entirely consist of sediment; it also includes nutrients, particulate organic matter, carbon, micro-plastics, colloidal particles, and various contaminants (e.g., D'Avignon et al., 2022). A suspended load thus affects the quality of water and

riverine ecological habitats (Suttle et al., 2004; Battin et al., 2008; Lloret et al., 2013; Koiter et al., 2013). It is therefore critical to assess this suspended load.

The simplest way to estimate the suspended load carried by a river is to filter a sample of water collected in the stream at a point where water is well mixed. The weight of the filter, once it is dried, gives a measure of the concentration of sediment, which, combined with the flow discharge, yields the rate of suspended sediment transport (Bierman and Montgomery, 2014). In practice, however, this long and tedious procedure is inappropriate for high-frequency measurements. Instead, in situ monitoring of the suspended load often relies on the use of a turbidimeter, an instrument capable of measuring the turbidity of the river at a high frequency and over long periods of time (Turowski et al., 2010; Esteves et al., 2019). Turbidity measures the amount of light scattered by the suspended particles in the water column. It is thus a convenient proxy for the suspended load concentration. Turbidity, however, also depends on the size distribution of the sus-
pended particles, on their shape, and on their chemical com-
position. Its conversion into a sediment concentration thus
requires an on-site calibration, which ultimately relies on the
filtering of water samples (Minella et al., 2008).

Field measurements show that the concentration of sus-
pended sediment fluctuates with the river discharge and cul-
minates during floods. Based on this observation, it is tempt-
ing to use discharge as a proxy for the suspended load (e.g.,
Ahn and Steinschneider, 2018). Yet, the relation between dis-
charge and concentration is not univocal: when observed at
the scale of a single flood event, it often exhibits a hysteretic
loop (e.g., Williams, 1989). These loops are observed under
various geological and climatic settings, independently of the
size of the catchment that the river drains (Langlois et al.,
2005; Bača, 2008; Eder et al., 2010; Ziegler et al., 2014).

Many factors combine to shape this hysteretic behavior. First of all, the suspended load concentration adjusts to the local shear stress over a characteristic time that depends on the flow depth and the particle settling velocity (Claudin et al., 2011). The resulting delay between discharge and sediment concentration induces a counterclockwise loop. In gravel-bed rivers, fine particles are often trapped below a layer of coarse sediments such as pebbles, a phenomenon known as armoring (Frey and Church, 2009; Ferdowsi et al., 2017a). The riverbed then acts as a sediment buffer that stores and releases fine particles according to its own dynamics, altering the shape of the concentration-discharge relationship (Orwin and Smart, 2004; Turowski et al., 2010; Park and Hunt, 2017; Guillon et al., 2018; Misset et al., 2019b). Fine particles, however, do not always originate from the riverbed. During storms, hillslope runoff, landsliding, and bank erosion may also feed the river with a significant quantity of fine particles (Hovius et al., 2000; Allemand et al., 2014). In this case, the spatial distribution of rainfall in the catchment area and the distance between the sediment sources and the sampling point influence the shape of the concentration-discharge relationship in a complex way (e.g., Asselman, 1999; Smith and Dragovich, 2009; Misset et al., 2019a). Finally, if the velocity of the flow that carries the suspended sediment is smaller than the celerity of the flood wave, the resulting delay between the discharge and the sediment peak induces a counterclockwise loop. This effect may dominate the discharge-concentration relationship if the distance traveled by the suspended particles is long enough (Klein, 1984; Nistor and Church, 2005).

The variety of processes involved in the transport of sus-
pended sediment has led to the development of different
types of models (Gao, 2008; Vercruysse and Grabowski,
2019). Following the pioneering work of Rouse (1939), a
large body of work focused on the formulation of mathe-
matical models that explain how the balance between turbu-
lent diffusion and sedimentation determines the vertical pro-
file of the suspended sediment concentration in a turbulent
flow (Van Rijn, 1984; Garcia and Parker, 1991; Wright and
Parker, 2004; Claudin et al., 2011). Combined with equations

for the flow, these so-called “hydromorphodynamic models”
provide a satisfactory physical picture of the transport of sus-
pended sediment in rivers (Van Rijn, 2007; Bouchez et al.,
2010; Armijos et al., 2017). However, their use to predict
sediment transport requires the calibration of a large number
of hydrological parameters as well as precise knowledge of
the river topography and discharge (Lepesqueur et al., 2019).

Hydrological models do not incorporate the input of fine sediment from the hillslopes surrounding the river. Yet, in a small catchment, this contribution may represent a significant fraction of the total sediment yield during rainfall (Gao, 2008). The need to account for hillslope processes motivated the development of “distributed models”. Distributed models break down the catchment into several regions. In each of them, the susceptibility to rainfall-induced soil erosion is parameterized by a series of indices that describe local properties such as soil characteristics, land cover, land use, and topography (Wischmeier and Smith, 1978; Renard et al., 1991, 2017). Based on the value of these indices, the model solves a system of equations that account for various hillslope and hydrological processes, and it predicts the sediment yield at the catchment outlet (De Aragão et al., 2005). These models explicitly account for the location of the sediment sources inside the catchment. However, they involve a large number of adjustable parameters, at the risk of overfitting the data (White, 2006).

Recently, the increase in the volume of environmental data
has led to the development of methods based on data min-
ing and artificial intelligence (Vercruysse and Grabowski,
2019). The principle is to use multivariate data analysis to
identify the variables that primarily control the concentration
of suspended sediment and reconstruct the function that re-
lates them together. This approach – which relies on methods
ranging from quantile regression forests to artificial neural
networks – successfully predicts the concentration of sus-
pended sediment in various contexts (Cobaner et al., 2009;
Bilotta et al., 2012; Zimmermann et al., 2012). In particular,
it allows us to identify correlations between complex hys-
teresis patterns and antecedent hydro-meteorological condi-
tions (Perks et al., 2015). This method, however, is essen-
tially a “black box” approach, and prediction often comes at
the cost of understanding (Vercruysse and Grabowski, 2019).
Besides, it requires the acquisition of large datasets over mul-
tiple timescales, which restricts its application to a few well-
instrumented catchments.

Despite this impressive body of work, modeling and pre-
dicting the entrainment, propagation, and deposition of sus-
pended particles are still challenging tasks. Here, we use data
acquired in the Capesterre river, a gravel-bed river that drains
a small watershed on the southeast coast of Basse-Terre Is-
land (Guadeloupe archipelago), to develop a phenomenolog-
ical model of suspended sediment transport. We begin with
a description of the Capesterre catchment, which is equipped
with a gauging station that records high-frequency measure-
ments of the river discharge and turbidity (Sect. 2). Based
<figure><img src="image_3.png" /></figure>

**Figure 1.** Capesterre catchment on Basse-Terre Island, Guadeloupe, French West Indies. (a) Topographic map of the catchment. Inset: map of Guadeloupe (IGN RGE ALTI® 10 m). Red dots locate the catchment outlet on both maps (16°04′19″N, 61°36′33″W). (b, c) Images of the river at low and high stages, respectively. (d) Concentration of suspended sediment versus turbidity in the Capesterre river. (e) Water level and (f) concentration of suspended sediment from 1 February to 1 April 2021.

on the data collected there, we develop a phenomenologi-
cal model of suspended sediment transport, which explic-
itly accounts for the exchange of small particles between
the riverbed and the water column (Sect. 3). We then test
this model against the data collected in Capesterre, discuss
the physical meaning of the model parameters, and demon-
strate its versatility by applying it to a second catchment,
the Draix–Laval catchment in the French Alps (Sect. 4). Fi-
nally, we discuss the strengths and limitations of the model
(Sect. 5) before concluding.

2 The Capesterre catchment

2.1 Field site and measurements

We begin our investigation in the Capesterre river, located on Basse-Terre Island. This volcanic island of the Guadeloupe archipelago belongs to the subduction arc of the Lesser Antilles (Feuillet et al., 2002) (Fig. 1a). The Basse-Terre climate is tropical: daily temperatures range between 24 and 28 °C, and the annual rainfall rate is about 5200 mm yr⁻¹. The combination of this tropical climate with a steep volcanic relief produces high erosion rates, which range between 800 and 4000 t km⁻² yr⁻¹ (Rad et al., 2006; Dessert et al., 2015). These values place Basse-Terre Island among

the fastest eroding spots on Earth (Summerfield and Hulton,
1994).

On Basse-Terre, rainfall is intermittent and occur mainly
as short, high-magnitude events. As a result, the discharge
of rivers varies abruptly, with frequent flash floods trigger-
gered by tropical rainfall and hurricanes (Fig. 1b and c).
These extreme climatic events, particularly frequent during
the rainy season from June to January, trigger landslides, rock
avalanches, and debris flows, which are the main drivers of
erosion (Allemand et al., 2014).

The Capesterre river drains a steep watershed on the wind-
ward side of the active Soufrière volcano (Fig. 1a). The mean
annual rainfall rate is 5700 mm yr⁻¹. However, topography
induces a strong orographic gradient, and the intensity of
rainfall is highly heterogeneous within the catchment.

The Capesterre catchment is underlain by an andesitic
bedrock, aged from 400 to 600 kyr (Samper et al., 2007).
Soils mostly consist of thin andosols of thickness typically
less than 1 m (Colmet-Daage and Bernard, 1979; Lloret et al.,
2016). The Capesterre river flows over 19.7 km, from its
headwater on the flanks of the active volcano, at an alti-
tude of 1390 m a.s.l., down to the Capesterre village, where
it discharges into the Atlantic Ocean. Its channel is made
of bedrock, partly covered by a thin layer of alluvial sedi-
ment. 3 km from the sea, the river suddenly turns alluvial as
its slope gradually decreases.
The Observatoire de l'Eau et de l'Érosion aux Antilles (ObsERA) operates a gauging station at the site of La Digue, at an altitude of 189 m a.s.l., a few hundred meters upstream of the point where the river turns alluvial (Fig. 1b and c). At the station, the river drains a catchment with an area of 16.4 km², almost entirely located within the boundaries of the National Park of Guadeloupe. There, a thick rainforest limits the input of sediment from hillslopes, and anthropogenic forcing is weak. Based on a 4-year water-sampling campaign, Lloret et al. (2013) estimated that the mean flux of suspended matter is about 153 t km⁻² yr⁻¹. About 10% of this flux consists of particulate organic carbon.

La Digue's station is equipped with a pressure sensor (CS451, Campbell Scientific Inc.) which measures the river stage relative to the fixed datum defined by a staff gauge installed by the Direction de l'Environnement, de l'Aménagement et du Logement (DEAL-Guadeloupe). In addition, a turbidimeter (OBS3+, Campbell Scientific Inc.) measures the turbidity of the water. Both pressure and turbidity sensors are connected to a data logger (Campbell CR800), which has been recording their respective measurements every 5 min since 2013 (Fig. 1e and f). The station is also equipped with an automatic water sampler (6712 Full-Size Portable Sampler) triggered by a pressure probe (ISCO 720 Submerged Probe module). This device collects water samples when the river stage exceeds a threshold set by the operator (about 50 cm). Filtration of these samples allows us to measure the concentration of suspended material C (expressed in mg L⁻¹) and to calibrate the relation between the latter and the turbidity T (NTU) (Fig. 1d). A fit of a power law through our data yields $C = \alpha T^n$, with $\alpha = 0.26 \pm 0.04$ mg L⁻¹ and $n = 1.23 \pm 0.04$.

The data acquired in the Capesterre river show that the water stage fluctuates as floods follow each other (Fig. 1e). Each flood starts with an abrupt increase of the water stage – the latter may rise from 0.1 to 1 m in less than 1 h – followed by a slow recession which lasts for 12 to 20 h. In between two floods, the water stage remains low (i.e., 10 to 20 cm). The concentration of suspended material, calculated from turbidity, is highly intermittent (Fig. 1f): it is virtually zero 98% of the time and rises only during floods, during which it may reach up to 1430 mg L⁻¹.

## 2.2 Concentration of suspended material

Sediment transport in the Capesterre river occurs only during floods. We therefore extract four flood events from our dataset and observe the evolution of the water level and of the concentration of suspended sediment in more detail (Figs. 2 and 3). All four floods exhibit the same simple shape: a steep rise followed by a slow decline. Visual inspection of the data reveals that this pattern is representative of the entire dataset, with multi-peak events representing less than a few percent of floods.

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

**Figure 2.** Three floods in the Capesterre river for which the concentration versus water-level relation forms a counterclockwise loop: (a) from 19:00 LT on 5 September to 10:00 LT on 6 September 2019, (b) from 22:00 LT on 22 July to 10:00 LT on 23 July 2019, (c) from 10:00 LT on 3 November to 00:00 LT 4 November 2021. Left panels: time series of the water level (blue line) and the concentration of suspended sediment (black line) measured at the gauging station. The dashed vertical lines indicate the moments when the concentration reaches its peak value. Right panels: relation between concentration and water level. Gray arrows indicate the direction of the hysteresis loops. On each panel, the red line is the concentration predicted from the best-fit model.

During a flood, the concentration of suspended sediment follows the fluctuations of the water stage: it is equal to zero before the flood, rises when the water stage exceeds a threshold of about 20 cm, increases with the water stage until the flood peak, and finally decreases back to zero during the recession limb of the flood (Fig. 2, left panels).

A closer look at the data, however, reveals a time lag between the flood peak and the concentration peak. For instance, for the three floods shown in Fig. 2, the concentration peak occurs about 30 min after the flood peak. A plot of the concentration as a function of the water stage better highlights this time lag between concentration and water level: the concentration versus water-level relation forms a counterclockwise hysteretic loop (Fig. 2, right panels).
<figure><img src="image_5.png" /><figcaption>Figure 3. Flood recorded in the Capesterre river from 07:00 LT on 15 January to 00:00 LT on 16 January 2021. During this event, the concentration versus water-level relation formed a clockwise loop. (a) Time series of the water level (blue line) and (b) the concentration of suspended sediment (black line) measured at the gauging station. (c) Relation between concentration and water level. Gray arrows indicate the direction of the hysteresis loops. The dashed vertical line indicates the moment when the concentration reaches its peak value. On each panel, the red line is the concentration predicted from the best-fit model.</figcaption></figure>

Because we deduce the concentration of suspended sediment from the river turbidity, interpreting the origin of concentration loops requires some caution. Indeed, concentration loops may sometimes be an artifact caused by a change of the turbidity-concentration relationship in the course of a flood (Landers and Sturm, 2013). This happens when the size distribution of the suspended load evolves significantly during a flood. The ObsERA observatory used a LISST-StreamSide (Laser In Situ Scattering and Transmissometry, Sequoia Scientific Inc.) to measure the size distribution of the particles in suspension in the Capesterre river during a flood in 2010. Their data reveal a bimodal distribution with two peaks (see Sect. 5.1.1). The first one, around 5 µm, corresponds to the wash load and remains remarkably stable during the flood. The position of the second peak varies between 23 and 32 µm. Assuming that this flood is representative, we expect no significant change of the turbidity-concentration relationship during a flood. Moreover, direct measurements of the concentration of suspended sediment, based on manual water sampling, confirm that the concentration versus water-level relation forms hysteretic loops (Lloret, 2010; Lloret et al., 2013). Based on these observations, we therefore attribute the hysteresis of the turbidity of the Capesterre river to a hysteresis of its concentration of suspended sediment.

In the Capesterre river, the flood peak is not systematically ahead of the concentration peak. Sometimes, it is the concentration peak that precedes the flood peak, and the concentration versus water-level relation then forms a clockwise hysteretic loop (Fig. 3).

To determine the proportion of clockwise versus counterclockwise loops in the Capesterre river, we first extract individual flood events from our dataset. To do so, we define a flood as a period during which the water stage is higher than a

threshold set to 20 cm. This value roughly corresponds to the threshold above which the flow is strong enough to carry suspended sediment. Based on this definition, we isolate 217 individual floods between May 2019 and December 2021.

We must now determine the direction of the concentration versus water-level relation for each of these floods. Inspired by Langlois et al. (2005) and Misset et al. (2019a), we define the hysteresis index for each flood as

$$I_H = \int \tilde{C} d\tilde{h}, \quad (1)$$

where we introduce the normalized water level, $\tilde{h} = (h - h_{min})/(h_{max} - h_{min})$, and the normalized concentration, $\tilde{C} = (C - C_{min})/(C_{max} - C_{min})$. $h_{min}$ and $h_{max}$ are the minimum and maximum values of the water level $h$ during the flood. Similarly, $C_{min}$ and $C_{max}$ are the minimum and maximum values of the concentration $C$ during the flood. With these definitions, both the normalized water level and the normalized concentration vary between 0 and 1. Accordingly, the hysteresis index $I_H$ ranges between -1 and 1. Negative values correspond to counterclockwise hysteresis, while positive ones indicate clockwise loops.

To characterize the proportion of clockwise versus counterclockwise loops, we just need to compute the hysteresis index $I_H$ for each of the 217 floods of our catalog. We find that this index ranges between -0.719 and 0.452, with a mean value of -0.030, a median of -0.021, and a standard deviation of 0.178. A plot of the cumulative distribution of the hysteresis index reveals that 61 % of the floods detected in Capesterre are characterized by a counterclockwise loop (Fig. 4a).

To detect a possible influence of the season, we categorize the flood events by wet (177 events) and dry (44 events) seasons. We then compute the distributions of the hysteresis index for each season and find that clockwise hysteresis and counterclockwise hysteresis are equally present in both seasons, with median values very close to zero: -0.007 for the dry season and -0.030 for the wet one (Fig. 4b).

To summarize, about 61 % of the floods of the Capesterre river exhibit a counterclockwise hysteresis loop, with no influence of the season. In the next section, we therefore focus on the counterclockwise hysteretic loops, and we formulate a simple model that accounts for their shape.

## 3 Phenomenological model of suspended sediment transport

In this section, we develop a model to account for the evolution of the concentration of suspended sediment in the Capesterre river. Given the complexity of the problem, our objective is not to establish a comprehensive physical theory of suspended sediment transport but rather to derive a phenomenological equation that reproduces our field measurements. We therefore start with a series of simplifying assumptions, driven by field observations.
<figure><img src="image_6.png" /><figcaption>Figure 4. (a) Cumulative distribution function (CDF) of the hysteresis index $I_H$ in Capesterre. Insets: cartoons illustrating the physical meaning of the hysteresis index $I_H$. The magnitude of $I_H$ is equal to the colored area, while its sign depends on the direction of the hysteresis loop: it is positive for a clockwise loop (blue area) and negative for a counterclockwise one (red area). (b) Box plot of the hysteresis indices according to the season. The box extends from $Q_1$ to $Q_3$ of the data. Its central line represents the median ($Q_2$). The whiskers extend from the box by $1.5 \times (Q_3 - Q_1)$. Fliers are outliers.</figcaption></figure>

We first note that the Capesterre catchment is densely vegetated. Accordingly, we assume that, during a flood, the quantity of fine sediment that hillslopes deliver to the river is negligible compared to that originating from the riverbed. We therefore treat the riverbed as the sole source of sediment. We furthermore assume that the quantity of fine sediment stored in the riverbed is so large that the entrainment of sediment is not limited by supply, a condition often referred to as “transport limited” (Dietrich et al., 2003). Of course, these two assumptions only hold at the scale of a flood event. Over the longer term, hillslope processes gradually replenish the river with sediments, compensating for the erosion of the bed during floods. The timescale of this process is uncertain, but it is likely much longer than a few weeks. During extreme rainfall events, landslides or rock avalanches as well as debris flows may also feed the river with large amounts of sediment. When this happens, the sediment bed cannot be

considered the only source of suspended sediment, and our model might fail to represent its concentration.

In the Capesterre river, the water level rarely exceeds 1 m, a value comparable to the median grain size of its sediment bed, on the order of 10 cm. As suspended particles mainly consist of silt and fine sand of size less than 100 µm (see Sect. 5.1.1), we expect that the turbulence induced by the bed roughness is high enough to homogenize the concentration of suspended sediment in the river, an assumption supported by images of the river during floods (Fig. 1c). Accordingly, we neglect any vertical or lateral gradient of the concentration. We discuss in more detail the limits of this assumption in Sect. 5.1.1.

Within the previous set of assumptions, conservation of the mass of suspended sediment reads

$$ \frac{\partial}{\partial t}(WhC) + \frac{\partial}{\partial x}(uWhC) = W(E - D), \quad (2) $$

where C is the concentration of suspended sediment (mass per unit volume), h is the water level in the river, u is the average flow velocity, W is the river width, x is the streamwise coordinate along the river course, and t is the time. E is the entrainment rate, i.e., the mass of sediment entrained from the riverbed per unit time and area. Conversely, D denotes the deposition rate, i.e., the mass of sediment deposited on the riverbed per unit time and area.

To compute the concentration of suspended sediment, we need to supplement Eq. (2) with expressions for the entrainment and deposition rates. Particles in suspension in the river settle under the action of gravity. In the Capesterre river, the concentration of particles never exceeds $0.9 \text{ g L}^{-1}$, a value small enough to neglect interactions between particles. In this dilute regime, we expect that the settling velocity of a particle is independent of the concentration. It depends, however, on the particle size. The suspended load usually consists of grains of different sizes, which thus settle with different velocities. However, following Van Rijn (1986), we further simplify our model and assume that the decline in total sediment concentration can be represented with a single representative settling velocity $V_s$. We shall discuss the limits of this assumption in more detail in Sect. 5.1.1. For the moment, we use it to derive the deposition rate which reads

$$ D = C V_s. \quad (3) $$

Let us now turn our attention to the entrainment rate. The river entrains sediment when the shear stress it exerts on its bed exceeds a threshold value (Shields, 1936). Calculating this shear stress requires additional equations, which also means additional assumptions and parameters. To avoid this inconvenience and to keep the model as simple as possible, we note that an increase of water level often implies an increase of bottom shear stress. In the following, we shall therefore use water level as a proxy for shear stress. Indeed, field data show that the Capesterre river carries sediment
only when the water level exceeds a threshold of approximately $h_t \approx 20$ cm. Accordingly, we look for an expression of the entrainment rate in the form $E = f(h - h_t)$, where $f$ is some unknown function. For lack of additional constraint, we choose the simplest possible form and propose the following:

$$E = \epsilon \left( \frac{h}{h_t} - 1 \right)^n H \left( \frac{h}{h_t} - 1 \right), \quad (4)$$

where $\epsilon$ is a characteristic entrainment rate (mass per unit time and bed area), $n$ is a dimensionless exponent, and $H$ is the Heaviside function. Keeping in mind that the water level $h$ acts as a proxy for the shear stress, Eq. (4) is similar to expressions of the entrainment rate commonly used in the literature (Bagnold, 1956; Van Rijn, 1986, 2007; Claudin et al., 2011).

Combining the expression of the deposition (Eq. 3) and erosion (Eq. 4) rates with the mass balance (Eq. 2) yields the following equation:

$$\frac{\partial(hC)}{\partial t} = \left[ \epsilon \left( \frac{h}{h_t} - 1 \right)^n H \left( \frac{h}{h_t} - 1 \right) - C V_s \right] - \frac{1}{W} \frac{\partial(WuhC)}{\partial x}, \quad (5)$$

where we assume that the river width does not change significantly over the duration of a flood. The evolution of the concentration of suspended sediment (Eq. 5, left term) is thus controlled by the balance between two terms. The first one (Eq. 5, right member, first term) is a local term which accounts for the difference between particle entrainment and deposition. The second one (Eq. 5, right member, second term) is an advective term, which accounts for the streamwise variations of the flux of suspended sediment $F = WuhC$.

Streamwise variations of the suspended sediment flux can result from an inhomogeneous distribution of sediment sources along the river course. We have, however, explicitly excluded this possibility from the model assumptions. The other possible cause of streamwise variations is the formation of a flood wave: during rainfall, groundwater and overland flow discharge water into the river and generate a surge that propagates downstream (Alsdorf et al., 2005; Guérin et al., 2019). This flood wave induces variations of the flow depth, velocity, and therefore of the sediment flux in the streamwise direction, whose effects are described by the spatial derivative of the advective term (Lepesqueur et al., 2019).

Before trying to solve Eq. (5), we perform a scaling analysis to estimate the magnitude of each of its terms (Barenblatt and Isaakovich, 1996). First, we note that Eq. (5) involves three parameters: the threshold water level $h_t$, the particle settling velocity $V_s$, and the characteristic entrainment rate $\epsilon$. Accordingly, we introduce the rescaled water level $\tilde{h} = h/h_t$, the rescaled concentration $\tilde{C} = (V_s/\epsilon)C$, and the rescaled time $\tilde{t} = t/\tau_s$, where we define the characteristic settling time $\tau_s = h_t/V_s$. The latter corresponds to the time nec-

essary for a particle to settle at the representative settling velocity $V_s$ over a height equal to the threshold height $h_t$. To turn Eq. (5) into a dimensionless equation, we still need a characteristic scale for the streamwise coordinate $x$. The $x$ derivative in Eq. (5) corresponds to the streamwise gradient of flux induced by the formation of a flood wave. We thus expect $x$ to scale like the characteristic length of the flood wave $\ell \sim U\tau_f$, where $\tau_f$ is the duration of the flood and $U$ is the characteristic flow velocity. We therefore define the rescaled streamwise coordinate as $\tilde{x} = x/(U\tau_f)$ and the rescaled flow velocity as $\tilde{u} = u/U$. With these new variables, we write the conservation of mass (Eq. 5) in its dimensionless form:

$$\frac{\partial(\tilde{h}\tilde{C})}{\partial \tilde{t}} = (\tilde{h} - 1)^n H (\tilde{h} - 1) - \tilde{C} - \frac{\tau_s}{\tau_f} \left[ \frac{1}{W} \frac{\partial}{\partial \tilde{x}} (W \tilde{u} \tilde{h} \tilde{C}) \right]. \quad (6)$$

Scaling analysis thus reveals that the relative magnitude of the advective term depends on the ratio between the characteristic particle settling time and the duration of the flood. When the flood lasts much longer than the settling time, the advective term can be neglected: the flood wave is so large that the streamwise gradient of the sediment flux becomes negligible. The evolution of the concentration of suspended material is then controlled by the local balance between the entrainment and the deposition rates. In the Capesterre river, field measurements indicate that (1) the threshold for grain entrainment is approximately $h_t \approx 20$ cm, and (2) the median grain size of the suspended load is about 40 µm (see Sect. 5.1.1). Assuming that the density of the particles is $\rho_s = 2700$ kg m$^{-3}$, this corresponds to a settling velocity of the order of $V_s \sim 1.5 \times 10^{-3}$ m s$^{-1}$. With these values, the characteristic settling time, $\tau_s \approx 130$ s, is very short compared to the typical duration of a flood, $\tau_f \gtrsim 12$ h, so that $\tau_s/\tau_f \lesssim 3 \times 10^{-3}$. Therefore, we neglect the advective term in the following.

With this approximation, and coming back to dimensional variables, the conservation of mass becomes

$$\frac{d\phi}{dt} = \epsilon \left( \frac{h}{h_t} - 1 \right)^n H \left( \frac{h}{h_t} - 1 \right) - V_s \frac{\phi}{h}, \quad (7)$$

where we introduce the mass of suspended sediment per unit area $\phi = hC$. Neglecting the advective term simplifies the model considerably: it reduces it to the ordinary differential equation in Eq. (7). The latter describes how the concentration of suspended sediment evolves in response to the water level $h(t)$, which acts as a forcing function. To compute the sediment concentration, we use the ODE function from SciPy (Virtanen et al., 2020) to solve Eq. (7) numerically. This procedure yields the mass of suspended sediment per unit of bed area, $\phi$, from which we deduce the concentration $C = \phi/h$.

Our model reduces the dynamics of suspended sediment transport to an exchange of particles between the bed and the river, driven by the water level. As the latter must rise before the river can entrain any sediment, our theory can only
produce counterclockwise loops of the concentration versus water-level relation. Finally, we did not establish Eq. (7) on a rigorous physical ground but derived it from phenomenological considerations. Equation (7) should thus be considered an ansatz, whose validity will depend on its ability to represent field data.

Equation (7) involves four parameters: the characteristic entrainment rate $\epsilon$, the threshold water level $h_t$, the settling velocity $V_s$, and the exponent $n$. To determine the ability of our model to reproduce field measurements, we must find the values of these four parameters that best fit the data. To this aim, we use the optimization procedure described in the next section.

# 4 Application of the model to the Capesterre river

## 4.1 Model adjustment procedure

To test the model against field data, we select a time period and extract the corresponding data. The latter form a collection of $N$ discrete values of concentration $C_i^d$ and water level $h_i^d$ measured at times $t_i$, with $i = 1, ..., N$. We then solve Eq. (7) numerically using the water-level data, $h_i^d$, as the forcing function. The resulting numerical solution depends on the values of the four parameters: $\epsilon$, $h_t$, $V_s$, and $n$. To compare it to the field data, we compute the theoretical concentration, $C_i^m$, for each of the $N$ discrete times $t_i$ and estimate the distance between the model and the data from the chi-square function:

$$ \chi^2 = \frac{1}{N - n_p} \sum_i (C_i^d - C_i^m)^2, \qquad (8) $$

where $n_p = 4$ is the number of parameters in our model.

To achieve the best fit between the model and the data, we need to determine the values of the parameters that minimize the chi-square function. For this purpose, we use the trust-region reflective algorithm method, implemented in the scipyoptimize.curve_fit function, and used in the model fitting wrapper of LMFIT (Non-Linear Least-Squares Minimization And Curve Fitting for Python) Python library (Newville et al., 2016).

This optimization procedure is sensitive to the duration of the time period over which we apply it. In the next section, we discuss the calibration of our model on field data over a short period of time corresponding to a single flood event. The case of longer time periods, which encompass several flood events, will be discussed in Sect. 4.3.

## 4.2 Calibration of the model on a single flood event

Our model can only produce counterclockwise loops of the concentration versus water-level relation. We therefore test it first on the three floods of Fig. 2, which exhibit such loops. To do so, we set the initial values of the parameters to $\epsilon = 3 \text{ mg m}^{-2} \text{ s}^{-1}$ for the entrainment rate, $V_s = 0.001 \text{ m s}^{-1}$

for the settling velocity, $h_t = 15 \text{ cm}$ for the threshold water level, and $n = 1$ for the exponent of the entrainment law. We then apply the optimization procedure described in the previous section and determine the parameters that best fit the data (Table 1). Despite its simplicity, the model reproduces the evolution of the concentration of suspended sediment measured in the field surprisingly well (Fig. 2).

A plot of the theoretical concentration as a function of the water level reveals that the model also accounts reasonably well for the hysteretic loop of the concentration versus water-level relation (Fig. 2, right). For the three flood events that we analyze, the model better represents the recession than the rise of the concentration. This may result from a bias of the optimization procedure which favors the recession limb, as the latter contains a greater number of data points than the flood rise.

The best-fit parameters vary from one flood to the other (Table 1). The threshold water level thus ranges from $h_t = 15$ to 23 cm, in agreement with the estimate based on a direct visualization of the data (Fig. 2). The settling velocity $V_s$ varies between $V_s = 6.2 \times 10^{-4}$ and $1.39 \times 10^{-3} \text{ m s}^{-1}$. As discussed later, these values are compatible with the settling velocity of quartz grains of size between 10 and 100 µm. The entrainment rate $\epsilon$ ranges between 0.641 and $21.1 \text{ mg m}^{-2} \text{ s}^{-1}$. We are not aware of any direct measurement of the rate of entrainment of fine particles from a sediment bed and therefore cannot assess whether these values are realistic or not. Finally, the exponent $n$ varies between 1.1 and 2.3.

Out of curiosity, we also tested our model against the flood of Fig. 3, for which the concentration versus water-level relation forms a clockwise hysteresis. As expected, the model utterly fails to represent the concentration of suspended sediment during this flood: it underestimates the concentration peak by about 45% and the concentration versus water-level relation is a mere line, instead of a hysteresis.

To summarize, Eq. (7) correctly models the transport of suspended sediment during floods, provided that the latter presents a counterclockwise loop of the relation between concentration and water level. Yet, the best-fit parameter values vary from one flood to the other. This raises the question of the ability of the model to represent the evolution of the concentration over a long time with a unique set of parameters. We address this problem in the next section.

## 4.3 Calibration of the model on a series of successive floods

To test the ability of Eq. (7) to represent sediment transport over a time longer than the duration of a single flood, we select a period of 4 d, marked by the occurrence of four successive floods, all of them characterized by a counterclockwise loop (Fig. 5). We then apply the optimization procedure and determine the parameters that best fit the data (Table 1). This unique set of parameters reproduces the evolution of the con-
<table><caption>Table 1. Values of the best-fit parameters, $\epsilon$, $h_t$, $n$, $V_s$, and hysteresis index $I_H$ of the floods displayed in Figs. 2, 3, 5, and 7. The grain size $d_s$ is inferred from the settling velocity $V_s$ based on Eq. (12). $\emptyset$ indicates irrelevant values ($I_H > 0$).</caption><thead><tr><th>Figure</th><th>Catchment</th><th>Period</th><th>I<sub>H</sub></th><th>ε [mg m<sup>-2</sup> s<sup>-1</sup>]</th><th>h<sub>t</sub> [cm]</th><th>n</th><th>V<sub>s</sub> [m s<sup>-1</sup>]</th><th>d<sub>s</sub> [µm]</th></tr></thead><tbody><tr><td>2a</td><td>Capesterre</td><td>19:00 LT on 5 Sep–10:00 LT on 6 Sep 2019</td><td>-0.194</td><td>21.1 ± 0.1</td><td>17.3 ± 0.1</td><td>1.580 ± 0.008</td><td>1.385 ± 0.006 × 10<sup>-3</sup></td><td>40.50</td></tr><tr><td>2b</td><td>Capesterre</td><td>22:00 LT on 22 Jul–22:00 LT on 23 Jul 2019</td><td>-0.236</td><td>16.12 ± 0.08</td><td>22.5 ± 0.1</td><td>1.122 ± 0.005</td><td>7.54 ± 0.08 × 10<sup>-4</sup></td><td>29.39</td></tr><tr><td>2c</td><td>Capesterre</td><td>10:00 LT on 10 Mar–00:00 LT on 11 Mar 2021</td><td>-0.278</td><td>0.641 ± 0.001</td><td>15.02 ± 0.03</td><td>2.300 ± 0.003</td><td>6.22 ± 0.01 × 10<sup>-4</sup></td><td>26.59</td></tr><tr><td>3</td><td>Capesterre</td><td>07:00 LT on 15 Jan–00:00 LT on 16 Jan 2021</td><td>+0.312</td><td>10<sup>-12</sup></td><td>0.34 ± 0.03</td><td>5.23 ± 0.08</td><td>0.26 ± 0.02</td><td>Ø</td></tr><tr><td>5</td><td>Capesterre</td><td>18:00 LT on 6 Feb–23:00 LT on 10 Feb 2021</td><td>-0.168</td><td>23.2 ± 0.09</td><td>20.0 ± 0.1</td><td>2.031 ± 0.007</td><td>9.71 ± 0.06 × 10<sup>-3</sup></td><td>125.43</td></tr><tr><td>7a</td><td>Draix–Laval</td><td>06:00 LT on 29 May–11:00 LT on 29 May 2016</td><td>-0.044</td><td>40.6 ± 0.1</td><td>7.22 ± 0.03</td><td>0.876 ± 0.004</td><td>3.12 ± 0.01 × 10<sup>-3</sup></td><td>63.29</td></tr><tr><td>7b</td><td>Draix–Laval</td><td>15:00 LT on 29 May–20:00 LT on 29 May 2016</td><td>-0.249</td><td>15.3 ± 0.1</td><td>10.5 ± 0.01</td><td>0.500 ± 0.004</td><td>5.08 ± 0.07 × 10<sup>-4</sup></td><td>23.94</td></tr></tbody></table>

concentration during the four successive floods reasonably well. The agreement is, however, not perfect: the model, for example, underestimates the amplitude of the second concentration peak by about 30% and overestimates the amplitude of the last one by 34% (Fig. 5).

To better assess the quality of the fit, we compute the mass of suspended sediment exported from the catchment by integrating the sediment flux over the entire duration of our time series, $M = \int Cqdt$, where $q$ is the flow discharge. By integrating the data, we obtained a mass of $M = 1.239 \times 10^3$ kg. The same calculation, conducted with the model, yields $M = 1.177 \times 10^3$ kg, which is about 5% smaller than the data. The model therefore provides a reasonable estimate of the mass exported out of the catchment.

Despite this encouraging result, we note that the model better predicts the data when its parameters are optimized on a single flood event, rather than on a series of floods. In short, the best-fit parameter values change from one flood to the next. This raises the question of their physical meaning, a topic we discuss in the next section.

# 5 Discussion

## 5.1 Physical meaning of the model parameters

The transport of suspended sediment depends on the properties of the sediment, on the specifics of the flow, and on the configuration of the catchment and the river. In our model, these properties are all lumped into the four parameters of Eq. (7). A change in the value of these parameters from one flood to the next therefore reflects a change of these properties.

The characteristic entrainment rate $\epsilon$, the exponent $n$, and the threshold of water level, $h_t$, for example, parameterize the expression of the rate at which particles are entrained from the bed (Eq. 4). Any change of these parameters thus reflects a modification of the conditions in which fine sediment is entrained from the bed. Physics suggests that the entrainment rate depends on the threshold of grain entrainment, on the sediment availability, and on the structure of the sediment bed. In our model, the threshold of water level, $h_t$, accounts for the threshold of grain entrainment. We therefore hypothesize that the characteristic entrainment rate $\epsilon$ and the exponent $n$ might be related to the sediment availability and the organization of the sediment bed. Explicitly formalizing this relation is, however, a difficult problem. Instead, we now turn our attention to the settling velocity $V_s$ and the threshold water level $h_t$.

### 5.1.1 Settling velocity

Our model represents the decline in sediment concentration after the flood peak with a single settling velocity $V_s$. In the Capesterre river, the concentration of suspended sediment is always small enough to neglect interactions between parti-
<figure><img src="image_7.png" /><figcaption>Figure 5. (a) Water level (blue line) and (b) concentration of suspended sediment in the Capesterre river between 18:00 LT on 6 February and 23:00 LT on 10 February 2021. (c) Concentration of suspended sediment versus water level. Black line: concentration measured from the turbidity. Orange line: concentration calculated from the flow depth based on an empirical rating curve. Red line: concentration predicted from the best-fit model.</figcaption></figure>

cles, and the settling velocity primarily depends on the grain size (Sect. 3). It is therefore tempting to try to infer the size of the suspended sediments from the value of $V_s$ deduced from the adjustment of the model to field data. To do this, we consider the suspended load to be a mixture of homogeneous particles of identical size $d_s$ and density $\rho_s$, settling in water of density $\rho = 1000 \text{ kg m}^{-3}$ and kinematic viscosity $\nu = 10^{-6} \text{ m}^2 \text{ s}^{-1}$.

Following Abraham (1970), we write the drag force exerted on a particle as

$$F_d = \frac{\pi}{8} C_d \rho d_s^2 V_s^2, \qquad (9)$$

where $C_d$ is a drag coefficient. The latter depends on the settling velocity through the Reynolds number $Re = V_s d_s / \nu$:

$$C_d = \left[ C_{\infty}^{1/2} + \left( \frac{24}{Re} \right)^{1/2} \right]^2, \qquad (10)$$

where $C_\infty$ is a constant that is about 1 in the case of natural grains (Andreotti et al., 2013). Equations (9) and (10) provide a general expression of the drag coefficient that represents the drag force from $Re = 0$ to $Re \approx 5000$ reasonably well (Abraham, 1970). In particular, at small Reynolds numbers, the drag coefficient simplifies into $C_d = 24/Re$, and one recovers the classical Stokes formula, $F_d = 3\pi \mu d_s V_s$. Conversely, for high Reynolds numbers, the drag force reduces to $F_d = \frac{\pi}{8} C_\infty \rho d_s^2 V_s^2$, as expected in a turbulent flow.

Balancing the drag force with the reduced weight of the particle, $F_g = (\pi/6)(\rho_s - \rho)gd_s^3$, we obtain the following expression for the settling velocity:

$$V_s = 6 \frac{\nu}{d_s} \left[ \left( \sqrt{\frac{(\rho_s - \rho)gd_s^3}{27\rho\nu^2}} + 1 \right)^{1/2} - 1 \right]^2, \qquad (11)$$

where $g$ is the acceleration of gravity. Inverting this equation yields the grain size as a function of the settling velocity:

$$d_s = \frac{24\nu}{\left[ \left( \sqrt{\frac{32(\rho_s - \rho)g\nu}{\rho V_s}} + \frac{V_s}{4} \right)^{1/2} - \left( \frac{V_s}{4} \right)^{1/2} \right]^2}. \qquad (12)$$

The best we can hope for is that Eq. (12) provides a value of the grain size that is representative of the suspended load. Indeed, our model assimilates the suspended load to a mixture of homogeneous particles of identical size. In reality, suspended load often consists of particles of different sizes and densities, which settle at different velocities. The conversion of settling velocity into a grain size only works as long as the size distribution of the grains in suspension is narrow.

Moreover, our model assumes that turbulence is high enough to homogenize the concentration of suspended sediment in the water column. In fact, laboratory experiments show the existence of a vertical gradient of suspended sediment concentration, even when the bed roughness is of the order of the flow depth (Grams and Wilcock, 2007). If such a concentration gradient forms in our river, the model prediction will likely overestimate the settling velocity by a factor equal to the ratio of the near-bed concentration to the averaged one. In this case, Eq. (12) underestimates the grain size.

Keeping these limitations in mind, we use Eq. (12) to estimate the grain size in the Capesterre river. The adjustment of our model to the field data yields the settling velocities reported in Table 1. Setting the density of sediment to $\rho_s = 2700 \text{ kg m}^{-3}$, we use Eq. (12) to turn these values into grain sizes. We find that the latter falls between 26 and 125 µm (Table 1), in agreement with the range of suspended sediment sizes reported in the literature (Sheldon et al., 1972; Wilcock et al., 2009).

For lack of simultaneous measurement of turbidity, grain size, and water level in the Capesterre river, we cannot compare the grain size calculated from the model with that measured in the field. In 2010, however, the observatory ObsERA
installed a LISST-StreamSide in the river and kept it running for a few months. The LISST-StreamSide (Laser In Situ Scattering and Transmissometry, Sequoia Scientific Inc.) is a laser particle sizer which measures the concentration of suspended particles in 32 logarithmically spaced size classes from 2 to 381 µm (Agrawal and Pottsmith, 2000). Assuming that the sediment size has not drastically changed since 2010, we compare the grain sizes measured by the LISST to that obtained from a fit of our model to the three 2019 floods of Fig. 2. Figure 6 shows the evolution of the grain size distribution during a flood, on 30 October 2010. We observe a smooth transition between two different behaviors, depending on the grain size. Below a diameter of about 10 µm, the concentration weakly depends on the water level: it is small, but not zero, before the flood, increases during the flood rise, and relaxes after the flood peak but does not return to zero over the duration of the recording. In this size range, particles settle very slowly on the bed, if at all. Based on these observations, we assimilate the particles of size less than 10 µm to the wash load, although recent investigations suggest that this conceptualization of wash load is oversimplified (Ren and Packman, 2007; Dallmann et al., 2020).

Above about 10 µm, the concentration of particles follows the evolution of the water level: it increases during the flood rise and relaxes to zero after the flood peak (Fig. 6b and c). We therefore interpret the corresponding population of particles as the suspended load. During the flood, their size ranges from 10 to 168 µm, a range that is consistent with the values calculated from the model. Encouraged by this result, we now turn our attention to the threshold water level.

### 5.1.2 Threshold water level

Our results show that the flow entrains fine sediment from the bed when the water level exceeds a threshold of approximately $h_t \approx 20$ cm. Although this value is computed for only four floods, it is, by far, the most robust result of our paper: direct observation of the data confirms that the concentration of suspended sediment rises only if the water level exceeds this threshold of about 20 cm, independently of the magnitude of the flood or the time of the year (Fig. 1e and f).

In our model, the water level is a proxy for the shear stress $\tau$ that the river exerts on its bed (Sect. 3). For sediment to be entrained by the flow, this stress must exceed a threshold $\tau_t$, a condition usually expressed in terms of the dimensionless Shields stress:

$$ \frac{\tau}{(\rho_s - \rho)gd_s}\theta_t, \qquad (13) $$

where $\theta_t$ is the threshold Shields stress, a dimensionless number that varies with the grain size and the flow regime (Shields, 1936; Andreotti et al., 2013). The fourth parameter of our model, the threshold water level, should therefore be linked to the threshold Shields stress.

To determine this link, we first note that the width of the Capesterre river is constant over a few hundred meters up-

<figure><img src="image_8.png" /><figcaption>Figure 6. Data acquired with a LISST-StreamSide in the Capesterre river on 30 October 2010. (a) Water level and (b) concentration of particles of median size 34.8 µm (class sizes between 31.9 and 37.7 µm) as a function of time. (c) Evolution of the concentration (color code) as a function of time (horizontal axis) for each of the 32 classes of grain sizes (vertical axis). Dashed line: transition between wash and suspended load. Dotted line: concentration profile displayed in (b).</figcaption></figure>

stream of the gauging station. There, assuming that the flow is uniform, we approximate the threshold shear stress by

$$ \tau_t = \rho g h_t S, \qquad (14) $$

where S is the slope of the river in the downstream direction. The latter, estimated from a digital elevation model (spatial resolution: 5 m), is about $S = 6.3 \times 10^{-2}$ at the La Digue station. Combining Eqs. (13) and (14) yields the relation between the threshold water level and the threshold Shields stress:

$$ h_t = \frac{(\rho_s - \rho)d_s\theta_t}{\rho S}. \qquad (15) $$

In Sect. 5.1.1, we found that the size of the suspended particles varies between 26 and 125 µm. The corresponding threshold Shields stresses range from $\theta_t \approx 0.1$ for the largest grains to $\theta_t \approx 0.2$ for the smallest ones, which are sensitive to cohesion forces (Shields, 1936; Van Rijn, 1984; Andreotti et al., 2013; Dunne et al., 2022). The threshold water level,
deduced from Eq. (15), thus varies between $h_t \approx 10^{-3}$ and $h_t \approx 10^{-2}$ cm. These values are inconsistent with the threshold $h_t \approx 20$ cm estimated from the model.

As the bed of the Capesterre river is covered with centimeter-size sediments (Fig. 1b), we suspect that this inconsistency might be the signature of bed armoring, an effect commonly observed in gravel-bed rivers as well as in aeolian systems (Ferdowsi et al., 2017b; Gao et al., 2016). The riverbed is said to be armored when a layer of coarse sediment overlies finer material, preventing it from being entrained in the flow, unless the armoring particles move first (Misset et al., 2021). If this scenario holds, the threshold of suspended sediment transport should coincide with that of the coarse particles. Measurements of the grain-size distribution of the bed of the Capesterre river, at the La Digue station, indicate that the riverbed is predominantly made of gravel and pebbles, with a median grain size of $d_{50} \approx 10$ cm. The corresponding threshold Shields stress is $\theta_t \approx 3 \times 10^{-3}$, for which Eq. (15) predicts a threshold water level $h_t \approx 10$ cm, comparable to the value deduced from our model. We therefore conclude that, in the Capesterre river, the threshold of suspended sediment transport is set by that of the coarse particles.

This result does not imply that fine sediments are stored in the subsurface, under a layer of coarse sediment. In fact, in the Capesterre river, silt and sand particles form tiny patches at the surface of the bed. These patches are trapped in the narrow space between neighboring pebbles. As the bed roughness is high (coarse sediment has a median size is about 10 cm), silts and sands are effectively screened from the flow and remain trapped between coarse particles until the latter start moving.

## 5.2 Application of the model to a small alpine catchment

Despite its simplicity, our model accounts for the transport of suspended sediment during floods, provided that the relation between concentration and water level forms a counterclockwise loop. We also found that two parameters of the model, the settling velocity and the threshold water level, can be used to estimate the size of the suspended sediment and to detect possible armoring of the bed. So far, however, we have only tested the model on data acquired in the Capesterre river. To test its versatility, we now apply it to the description of suspended sediment transport in a completely different geological context, that of the Laval catchment.

The Laval catchment, located in the southern part of the French Alps, is a small catchment of area $0.86 \text{ km}^2$, drained by the Laval stream (Fig. 7). Its elevation ranges from 850 to 1250 m a.s.l., and the annual rainfall rate is about $900 \text{ mm yr}^{-1}$; it experiences heavy rainfall events during spring and summer and less intense but longer rainfall in autumn. Besides, this catchment can also experience snowfall during winter (Ariagno et al., 2022). Unlike Capesterre, vegetation covers only 32% of the Laval catchment (Car-

riere et al., 2020) (Fig. 7a). The latter is underlain by easily erodable Middle Jurassic black marls, leading to the formation of steep-slope badlands within the catchment (Ariagno et al., 2022). A gauging station operated by the Draix-Bléone Critical Zone Observatory monitors the concentration of suspended sediment and the water level in the stream. The corresponding data are available on the Draix-Bléone Critical Zone Observatory website.

Figure 7 displays the water level and the concentration of suspended sediment during two floods, extracted from the catalog of the Draix-Bléone Critical Zone Observatory. The water level ranges from 5 to 25 cm, and the concentration reaches up to $35 \text{ g L}^{-1}$, which is about 1000 times higher than the maximum concentration measured in Capesterre. During these two floods, the concentration versus water-level relation forms a counterclockwise hysteretic loop, which makes them, in principle, compatible with our model (Fig. 7b and c, right panels).

Encouraged by this observation, we apply the optimization procedure introduced in Sect. 4.1 and determine the parameters that best fit the data (Table 1). As in Capesterre, the model represents rather well the evolution of the concentration of suspended sediment measured in the field (Fig. 7b and c) with similar $\chi^2$ values ($< 0.04$). Setting aside the exponent $n$, we note that the best-fit parameter values fall within the same ranges as those obtained in Capesterre (Table 1).

A comprehensive evaluation of the model would require a systematic test against data measured not only in Draix-Laval but also in several others catchments. Such a work is beyond the scope of the present paper. For the time being, we simply note that the model seems able to reproduce sediment transport in a context different than that of Capesterre.

## 5.3 Strengths and limitations of the model

In this paper, we adopt a systemic approach: instead of focusing on details, we treat the river as a uniform reservoir of homogenous sediment, and we reduce the dynamics of suspended sediment to an exchange of particles between the bed and the flow, driven by the water level. The resulting model describes the transport of suspended sediments by means of an ordinary differential equation, while the characteristics of the sediment and the river – such as the grain size distribution, the availability of sediment, or the threshold shear stress necessary to set sediment in motion – are all lumped into four parameters: a representative settling velocity, a threshold water level, a characteristic entrainment rate, and a dimensionless exponent.

Despite its simplicity, the model represents reasonably well the transport of suspended sediment in two small catchments, Capesterre and Draix-Laval, provided that the relation between concentration and water level forms a counterclockwise loop. The simple exchange of particles between the bed and the flow may therefore suffice to form these
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Floods in the Draix-Laval catchment. (a) Aerial view and boundaries of the catchment (© Google map); the red dot locates the outlet where the Laval station is (44°08′26.7″ N 6°21′39.4″ E). (b, c) Two floods recorded on 29 June 2016, from 06:00 to 11:00 LT and from 15:00 to 20:00 LT, respectively. Left panels: time series of the water level (blue line) and the concentration of suspended sediment (black line) measured at the gauging station. Right panels: relation between concentration and water level. Gray arrows indicate the direction of the hysteresis loops. On each panel, the red line is the concentration predicted from the best-fit model.</figcaption></figure>

loops, which are not necessarily the signature of a particu-
lar configuration of sediment sources or of a difference in
speed between the flow and the flood wave.

In many catchments, the sediment yield of the river is
calculated from the flow stage – or, equivalently, the flow
discharge – based on a sediment rating curve (Bierman and

Montgomery, 2014). In contrast to our model, this approach
assumes that the sediment concentration is a one-to-one func-
tion of the water stage and cannot represent hysteretic loops
of the relation between concentration and water level. Yet,
out of curiosity, we compare the prediction of our model with
the concentration estimated from the sediment rating curve
deduced from direct measurements of the concentration of suspended sediment in the Capesterre river. For the series of floods displayed in Fig. 5, the rating curve systematically overestimates the concentration measured by the turbidity sensor. As a result, the sediment yield calculated from the rating curve exceeds the one measured with the turbidimeter by about 54%. By comparison, the prediction of the model falls within 5% of the data. These observations do not necessarily disqualify the use of the sediment rating curve. Indeed, the latter is constructed from samples collected by an automatic water sampler. The latter collects 24 samples, one every 15 to 60 min, as soon as the river stage exceeds a threshold set by the operator. In the dataset at our disposal, this threshold varies between 30 and 80 cm, which is way above the threshold of entrainment in the Capesterre river. As a result, most concentration data belong to the falling rather than the rising part of the hysteresis. This bias likely explains why, in our case, the rating curve overestimates the concentration of suspended sediment in the river. Unfortunately, we do not have enough measurements of the concentration to correct this bias.

By construction, our model can only underestimate the sediment load during floods that develop a clockwise relation of concentration versus water level. Yet, even during such an event, the sediment yield calculated from the rating curve exceeds the one measured with the turbidimeter by about 32%, while the prediction of the model falls within 3% of the data (Fig. 3).

Obviously, our model is far too simple to incorporate all the processes involved in the transport of suspended sediment. Yet, simplicity has its advantages. First, our model relies on only four parameters. As a result, its computational cost is very small, compared to the one required by more complex approaches, such as hydromorphodynamic or distributed models (Vercruysse and Grabowski, 2019; Lepesqueur et al., 2019). Indeed, the calibration of our model against a single flood, such as the ones of Fig. 2, only takes a few minutes of computational time. This small number of parameters also reduces the risk of data overfitting, which, sometimes, hinders the use of more complex models (Gao, 2008).

Secondly, each parameter of our model admits a physical interpretation: the representative settling velocity is related to the size of the suspended sediment; the threshold water level acts as a proxy for the threshold shear stress; the characteristic entrainment rate and the dimensionless exponent parameterize the ability of the flow to entrain sediment out of the bed. One may therefore use these parameters to infer information about the properties of the catchment, however imperfect this information might be. In the Capesterre river, for example, the values of the parameters calculated from a fit of the data to our model change from one flood to the next, probably reflecting changes of the river and sediment properties. Changes of the settling velocity $V_s$, for example, reflect changes of the size distribution of the suspended sediment.

Last, but not least, the calibration of our model requires a single field station. This makes it particularly suitable for the analysis of data acquired in small catchments, as the latter are often monitored with a single station at the basin outlet (Gaillardet et al., 2018).

Simplicity comes at a cost, and our model suffers from many limitations. First, it is unable to represent clockwise loops of the concentration versus water-level relationship. In the Capesterre river, such events amount to about 40% of the floods and are therefore far from anecdotal. In practice, this inability to account for clockwise loops limits the use of the model to short time intervals, over which the relationship of concentration to water-level remains counterclockwise.

Based on field observations, we suspect that, in the Capesterre river, clockwise loops form when the stock of fine particles stored in the sediment bed is too small. In this supply-limited configuration, the entire stock of fine sediment is rapidly entrained by the flow, and the concentration reaches its maximum prior to the hydrograph peak (Williams, 1989). We could incorporate this effect into our model by introducing a second equation that accounts for the quantity of fine sediment stored in the bed. In this improved model, the rate at which the flow entrains sediment from the bed would depend not only on the water level but also on the quantity of fine sediment available on the bed. Such an improvement is a work in progress.

Equation (7) assumes that the contribution of the advective term is negligible. This assumption, which considerably simplifies the model, holds as long as the representative settling time of the particles is short compared to the duration of the flood, $h_t/V_s \ll \tau_f$ (Sect. 3). This condition depends on the size of the suspended sediment and on the river depth. It is satisfied in the Capesterre river. It might not be the case in other rivers.

As all theories, our model is susceptible to predicting the correct result for the wrong reasons. In lowland rivers, for example, the slope is shallow and flood waves may travel faster than the particles suspended in the flow. The resulting concentration versus water-stage relation then forms a counterclockwise loop. We suspect that a fit of our model to such an event would correctly represent the data. In this case, however, our model, which does not incorporate the propagation of flood waves, would be right for the wrong reasons, and the parameter values would be meaningless. Such a configuration requires more sophisticated equations that explicitly account for the flow (Lepesqueur et al., 2019). Our model should thus be used with caution and its parameters be interpreted with care.

Finally, the model that we propose is far too simple to incorporate all the mechanisms that may influence the transport of suspended sediment. In particular, it assimilates the suspended load to a suspension of grains of uniform and constant size and density. Yet, recent investigations suggest that flocculation and complex interaction with the bed may change the size distribution of fine sediment over the
course of a flood (Ren and Packman, 2007; Dallmann et al., 2020; Lamb et al., 2020). As a matter of fact, the chemical composition and the size of the particles suspended in the Capesterre river change with time: at the beginning of a flood, the suspended load is dominated by litter debris and Allophanes; sand particles appear later, when the flow discharge is sufficiently high (Céline Dessert, personal communication, 2022). Modeling such a complex behavior is a difficult problem, far beyond the abilities of the model presented in this work.

## 6 Conclusions

Based on data acquired in the Capesterre river, a small tropical river in Guadeloupe, we develop a phenomenological model that describes the transport of suspended sediment. Instead of focusing on details, our model adopts a systemic approach, and it treats the river as a uniform reservoir of homogeneous sediment entrained by the flow. It accounts for neither the propagation of a flood wave nor the consequences of an inhomogeneous sediment source. In addition, the model assumes that the entrainment of sediment is not limited by supply, a condition often referred to as “transport limited”. In short, we reduce the transport of suspended sediment to a uniform exchange of particles between the bed and the flow, driven by the water level.

Despite these simplifications, the model represents reasonably well the transport of suspended sediment in two small catchments, Capesterre and Draix-Laval, provided that the relation between concentration and water level forms a counterclockwise loop. The model, however, cannot represent clockwise loops of the concentration versus water-level relationship. This restricts its use to short time intervals, during which the concentration versus water-level relationship remains counterclockwise.

Our model describes the transport of suspended sediments by means of an ordinary differential equation, while the characteristics of the sediment and the river – such as the grain size distribution, the availability of sediment, or the threshold shear stress necessary to set sediment in motion – are all lumped into four parameters: a settling velocity related to the size of the suspended sediment, a threshold water level which acts as a proxy for the threshold shear stress, a characteristic entrainment rate, and a dimensionless exponent, which are related to the sediment availability. This simplicity gives the model a low computational cost. In addition, calibration of the model requires a single field station. This makes it particularly suitable for the analysis of data acquired in small catchments, as the latter are often monitored with a single station at the basin outlet (Gaillardet et al., 2018).

In the Capesterre river, the value of the parameters deduced from a fit of the model to the field data changes from one flood to the next, probably reflecting changes of the river and sediment properties. Changes of the settling velocity $V_s$,

for example, reflect changes of the size distribution of the suspended sediment. Similarly, the threshold water level, $h_t$, is a proxy for the threshold stress necessary to entrain sediment from the bed. As for the entrainment rate $\epsilon$ and the exponent $n$, we suspect that they are linked to the availability of fine sediment in the riverbed. If these conjectures prove to be correct, our model might offer the opportunity to detect river-wide changes, like modifications of the sediment size and availability induced by bank incision or landsliding.

**Code availability.** The code is available at https://doi.org/10.5281/zenodo.7893930 (Roque-Bernard et al., 2023).

**Data availability.** Obsera data are available at https://in-situ.theia-land.fr/ (OZCAR, 2023) and DRAIX data are available at https://bdoh.irstea.fr/DRAIX/ (INRAE, 2023).

**Author contributions.** ARB drove the science, wrote the model, and produced the figures. AL and EL initiated the research. EL provided the original concept of the model. AL provided his expertise in inverse problems and modeling. EG provided the background on erosion in tropical environments. PA and CD provided data and expertise that support this research. All authors contributed to the paper writing and sharing ideas.

**Competing interests.** The contact author has declared that none of the authors has any competing interests.

**Disclaimer.** Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Acknowledgements.** The data used in this paper were collected by two observatories of the OZCAR research infrastructure: the Observatoire de l'Eau et de l'Érosion aux Antilles (ObsERA, INSU-CNRS, http://webobsera.ipgp.fr/, last access: 10 May 2022), located in Guadeloupe, and the Draix-Bléone Critical Zone Observatory (https://bdoh.irstea.fr/DRAIX/, last access: 17 May 2022), located in the French Alps. These data are freely available from the website of each observatory. The authors thank the anonymous referees and associate editor for their insightful review that helped to improve the paper. The authors thank the members of ObsERA, especially Thierry Kitou, Vincent Robert, Olivier Crispi, Jolaine Ajax, and Brenda Lamaille. We are also grateful to Caroline Le Bouteiller for sharing her knowledge and expertise on the Laval catchment.

**Financial support.** This research has been supported by the Labex UnivEarthS (grant nos. ANR-10-LABX-0023 and ANR-18-IDEX-0001).
**Review statement.** This paper was edited by Paola Passalacqua and reviewed by two anonymous referees.

## References

Abraham, F. F.: Functional dependence of drag coefficient of a sphere on Reynolds number, Phys. Fluids, 13, 2194–2195, 1970.

Agrawal, Y. C. and Pottsmith, H. C.: Instruments for particle size and settling velocity observations in sediment transport, Mar. Geol., 168, 89–114, 2000.

Ahn, K.-H. and Steinschneider, S.: Time-varying suspended sediment-discharge rating curves to estimate climate impacts on fluvial sediment transport, Hydrol. Process., 32, 102–117, 2018.

Allemand, P., Delacourt, C., Lajeunesse, E., Devauchelle, O., and Beauducel, F.: Erosive effects of the storm Helena (1963) on Basse Terre Island (Guadeloupe – Lesser Antilles Arc), Geomorphology, 206, 79–86, https://doi.org/10.1016/j.geomorph.2013.09.020, 2014.

Alsdorf, D., Dunne, T., Melack, J., Smith, L., and Hess, L.: Diffusion modeling of recessional flow on central Amazonian floodplains, Geophys. Res. Lett., 32, L21405, https://doi.org/10.1029/2005GL024412, 2005.

Andreotti, B., Forterre, Y., and Pouliquen, O.: Granular media: between fluid and solid, Cambridge University Press, ISBN 9781107034792, 2013.

Ariagno, C., Le Bouteiller, C., van der Beek, P., and Klotz, S.: Sediment export in marly badland catchments modulated by frost-cracking intensity, Draix–Bléone Critical Zone Observatory, SE France, Earth Surf. Dynam., 10, 81–96, https://doi.org/10.5194/esurf-10-81-2022, 2022.

Armijos, E., Crave, A., Espinoza, R., Fraizy, P., Santos, A. D., Sampaio, F., De Oliveira, E., Santini, W., Martinez, J. M., Autin, P., Pantoja, N., Oliveira, M., and Filizola, N.: Measuring and modeling vertical gradients in suspended sediments in the Solimões/Amazon River, Hydrol. Process., 31, 654–667, https://doi.org/10.1002/hyp.11059, 2017.

Asselman, N. E.: Suspended sediment dynamics in a large drainage basin: the River Rhine, Hydrol. Process., 13, 1437–1450, 1999.

Bača, P.: Hysteresis effect in suspended sediment concentration in the Rybárik basin, Slovakia, Hydrolog. Sci. J., 53, 224–235, https://doi.org/10.1623/hysj.53.1.224, 2008.

Bagnold, R.: The nature of saltation and of Bedload transport in water, P. Roy. Soc. Lond. A, 332, 473–504, 1973.

Bagnold, R. A.: The flow of cohesionless grains in fluids, Philos. T. Roy. Soc. Lond. A, 249, 235–297, 1956.

Barenblatt, G. I. and Isaakovich, B. G.: Scaling, self-similarity, and intermediate asymptotics: dimensional analysis and intermediate asymptotics (Cambridge Texts in Applied Mathematics), Cambridge University Press, https://doi.org/10.1017/CBO9781107050242, 1996.

Battin, T. J., Kaplan, L. A., Findlay, S., Hopkinson, C. S., Marti, E., Packman, A. I., Newbold, J. D., and Sabater, F.: Biophysical controls on organic carbon fluxes in fluvial networks, Nat. Geosci., 1, 95–100, 2008.

Bierman, P. R. and Montgomery, D. R.: Key concepts in geomorphology, W. H. Freeman and Company Publishers, A Macmillan Higher Education Company, ISBN 978-1319059804, 2014.

Bilotta, G. S., Burnside, N. G., Cheek, L., Dunbar, M. J., Grove, M. K., Harrison, C., Joyce, C., Peacock, C., and Davy-Bowker,

J.: Developing environment-specific water quality guidelines for suspended particulate matter, Water Res., 46, 2324–2332, 2012.

Bouchez, J., Lajeunesse, E., Gaillardet, J., France-Lanord, C., Dutra-Maia, P., and Maurice, L.: Turbulent mixing in the Amazon River: The isotopic memory of confluences, Earth Planet. Sc. Lett., 290, 37–43, 2010.

Carriere, A., Le Bouteiller, C., Tucker, G. E., Klotz, S., and Naaim, M.: Impact of vegetation on erosion: Insights from the calibration and test of a landscape evolution model in alpine badland catchments, Earth Surf. Proc. Land., 45, 1085–1099, 2020.

Claudin, P., Charru, F., and Andreotti, B.: Transport relaxation time and length scales in turbulent suspensions, J. Fluid Mech., 671, 491–506, https://doi.org/10.1017/S0022112010005823, 2011.

Cobaner, M., Unal, B., and Kisi, O.: Suspended sediment concentration estimation by an adaptive neuro-fuzzy and neural network approaches using hydro-meteorological data, J. Hydrol., 367, 52–61, 2009.

Colmet-Daage, F. and Bernard, Z.: Contribution à l'Atlas des départements d'Outre-mer: Guadeloupe, Carte des sols de la Guadeloupe, Grande-Terre, Marie-Galante. Carte des pentes et du modelé de la Guadeloupe, Grande-Terre, Marie-Galante, ORSTOM, Antilles, https://horizon/documentation.ird.fr/exl-doc/pleins_textes/divers16-01/29909.pdf (last access: May 2022), 1979.

Dallmann, J., Phillips, C. B., Teitelbaum, Y., Sund, N., Schumer, R., Arnon, S., and Packman, A. I.: Impacts of suspended clay particle deposition on sand-bed morphodynamics, Water Resour. Res., 56, e2019WR027010, https://doi.org/10.1029/2019WR027010, 2020.

D'Avignon, G., Gregory-Eaves, I., and Ricciardi, A.: Microplastics in lakes and rivers: an issue of emerging significance to limnology, Environ. Rev., 30, 228–244, https://doi.org/10.1139/er-2021-0048, 2022.

De Aragão, R., Srinivasan, V. S., Suzuki, K., Kadota, A., Oguro, M., and Sakata, Y.: Evaluation of a physically-based model to simulate the runoff and erosion processes in a semiarid region of Brazil, Sediment Budgets, 2, 85–93, 2005.

Dessert, C., Lajeunesse, E., Lloret, E., Clergue, C., Crispi, O., Gorge, C., and Quidelleur, X.: Controls on chemical weathering on a mountainous volcanic tropical island: Guadeloupe (French West Indies), Geochim. Cosmochim. Ac., 171, 216–237, 2015.

Dietrich, W. E., Bellugi, D. G., Sklar, L. S., Stock, J. D., Heimsath, A. M., and Roering, J. J.: Geomorphic transport laws for predicting landscape form and dynamics, Geophys. Monogr. Ser., 135, 103–132, https://doi.org/10.1029/135GM09, 2003.

Dunne, K. B. J., Arratia, P. E., and Jerolmack, D. J.: A New Method for In Situ Measurement of the Erosion Threshold of River Channels, Water Resour. Res., 58, e2022WR032407, https://doi.org/10.1029/2022WR032407, 2022.

Eder, A., Strauss, P., Krueger, T., and Quinton, J. N.: Comparative calculation of suspended sediment loads with respect to hysteresis effects (in the Petzenkirchen catchment, Austria), J. Hydrol., 389, 168–176, https://doi.org/10.1016/j.jhydrol.2010.05.043, 2010.

Esteves, M., Legout, C., Navratil, O., and Evrard, O.: Medium term high frequency observation of discharges and suspended sediment in a Mediterranean mountainous catchment, J. Hydrol., 568, 562–574, 2019.
Ferdowsi, B., Ortiz, C. P., Houssais, M., and Jerolmack, D. J.: Riverbed armouring as a granular segregation phenomenon, Nat. Commun., 8, 1–10, 2017a.

Ferdowsi, B., Ortiz, C. P., Houssais, M., and Jerolmack, D. J.: Riverbed armouring as a granular segregation phenomenon, Nature Commun., 8, 1363, https://doi.org/10.1038/s41467-017-01681-3, 2017b.

Feuillet, N., Manighetti, I., Tapponnier, P., and Jacques, E.: Arc parallel extension and localization of volcanic complexes in Guadeloupe, Lesser Antilles, J. Geophys. Res., 107, 2331, https://doi.org/10.1029/2001JB000308, 2002.

Frey, P. and Church, M.: How river beds move, Science, 325, 1509–1510, 2009.

Gaillardet, J., Braud, I., Hankard, F., Anquetin, S., Bour, O., Dorfliger, N., De Dreuzy, J.-R., Galle, S., Galy, C., Gogo, S., Gourcy, L., Habets, F., Laggoun, F., Longuevergne, L., Le Borgne, T., Naaim-Bouvet, F., Nord, G., Simonneaux, V., Six, D., Tallec, T., Valentin, C., Abril, G., Allemand, P., Arènes, A., Arfib, B., Arnaud, L., Arnaud, N., Arnaud, P., Audry, S., Bailly, Bailly Comte, V., Batiot, C., Battais, A., Bellot, H., Bernard, E., Bertrand, C., Bessière, H., Binet, S., Bodin, J., Bodin, X., Boithias, L., Bouchez, J., Boudevillain, B., Bouzou Moussa, I., Branger, F., Braun, J. J., Brunet, P., Caceres, B., Calmels, D., Cappelaere, B., Celle-Jeanton, H., Chabaux, F., Chalikakis, K., Champollion, C., Copard, Y., Cotel, C., Davy, P., Deline, P., Delrieu, G., Demarty, J., Dessert, C., Dumont, M., Emblanch, C., Ezzahar, J., Estèves, M., Favier, V., Faucheux, M., Filizola, N., Flammarion, P., Floury, P., Fovet, O., Fournier, M., Francez, A. J., Gandois, L., Gascuel, C., Gayer, E., Genthon, C., Gérard, M. F., Gilbert, D., Gouttevin, I., Grippa, M., Gruau, G., Jardani, A., Jeanneau, L., Join, J. L., Jourde, H., Karbou, F., Labat, D., Lagadeuc, Y., Lajeunesse, E., Lastennet, R., Lavado, W., Lawin, E., Lebel, T., Le Bouteiller, C., Legout, C., Lejeune, Y., Le Meur, E., Le Moigne, N., Lions, J., Lucas, A., Malet, J. P., Marais-Sicre, C., Maréchal, J. C., Marlin, C., Martin, P., Martins, J., Martinez, J. M., Massei, N., Mauclerc, A., Mazzilli, N., Molénat, J., Moreira-Turcq, P., Mougin, E., Morin, S., Ndam Ngoupayou, J., Panthou, G., Peugeot, C., Picard, G., Pierret, M. C., Porel, G., Probst, A., Probst, J. L., Rivière, A., Robain, H., Ruiz, L., Sanchez-Perez, J. M., Santini, W., Sauvage, S., Schoeneich, P., Seidel, J. L., Sekhar, M., Sengtaheuanghoung, O., Silvera, N., Steinmann, M., Soruco, A., Tallec, G., Thibert, E., Valdes Lao, D., Vincent, C., Viville, D., Wagnon, P., and Zitouna, R.: OZCAR: The French network of critical zone observatories, Vadose Zone J., 17, 1–24, 2018.

Gao, P.: Understanding watershed suspended sediment transport, Prog. Phys. Geogr., 32, 243–263, 2008.

Gao, X., Narteau, C., and Rozier, O.: Controls on and effects of armoring and vertical sorting in aeolian dune fields: A numerical simulation study, Geophys. Res. Lett., 43, 2614–2622, 2016.

Garcia, M. and Parker, G.: Entrainment of bed sediment into suspension, J. Hydraul. Eng., 117, 414–435, 1991.

Grams, P. E. and Wilcock, P. R.: Equilibrium entrainment of fine sediment over a coarse immobile bed, Water Resour. Res., 43, W10420, https://doi.org/10.1029/2006WR005129, 2007.

Guérin, A., Devauchelle, O., Robert, V., Kitou, T., Dessert, C., Quiquerez, A., Allemand, P., and Lajeunesse, E.: Stream-discharge surges generated by groundwater flow, Geophys. Res. Lett., 46, 7447–7455, 2019.

Guillon, H., Mugnier, J.-L., and Buoncristiani, J.-F.: Proglacial sediment dynamics from daily to seasonal scales in a glaciated Alpine catchment (Bossons glacier, Mont Blanc massif, France), Earth Surf. Proc. Land., 43, 1478–1495, 2018.

Hovius, N., Stark, C. P., Hao-Tsu, C., and Jiun-Chuan, L.: Supply and removal of sediment in a landslide-dominated mountain belt: Central Range, Taiwan, J. Geol., 108, 73–89, 2000.

INRAE: Observatoire Draix-Bleone, https://bdoh.irstea.fr/DRAIX/, last access: 4 May 2023.

Klein, M.: Anti clockwise hysteresis in suspended sediment concentration during individual storms: Holbeck Catchment, Yorkshire, England, Catena, 11, 251–257, 1984.

Koiter, A. J., Lobb, D. A., Owens, P. N., Petticrew, E. L., Tiessen, K. H., and Li, S.: Investigating the role of connectivity and scale in assessing the sources of sediment in an agricultural watershed in the Canadian prairies using sediment source fingerprinting, J. Soils Sediments, 13, 1676–1691, 2013.

Lamb, M. P., de Leeuw, J., Fischer, W. W., Moodie, A. J., Venditti, J. G., Nittrouer, J. A., Haught, D., and Parker, G.: Mud in rivers transported as flocculated and suspended bed material, Nat. Geosci., 13, 566–570, https://doi.org/10.1038/s41561-020-0602-5, 2020.

Landers, M. N. and Sturm, T. W.: Hysteresis in suspended sediment to turbidity relations due to changing particle size distributions, Water Resour. Res., 49, 5487–5500, 2013.

Langlois, J. L., Johnson, D. W., and Mehuys, G. R.: Suspended sediment dynamics associated with snowmelt runoff in a small mountain stream of Lake Tahoe (Nevada), Hydrol. Process., 19, 3569–3580, https://doi.org/10.1002/hyp.5844, 2005.

Lepesqueur, J., Hostache, R., Martínez-Carreras, N., Montargès-Pelletier, E., and Hissler, C.: Sediment transport modelling in riverine environments: on the importance of grain-size distribution, sediment density, and suspended sediment concentrations at the upstream boundary, Hydrol. Earth Syst. Sci., 23, 3901–3915, https://doi.org/10.5194/hess-23-3901-2019, 2019.

Liu, Y., Métivier, F., Gaillardet, J., Ye, B., Meunier, P., Narteau, C., Lajeunesse, E., Han, T., and Malverti, L.: Erosion rates deduced from Seasonal mass balance along an active braided river in Tian Shan, Solid Earth Discuss., 3, 541–589, 2011.

Lloret, E.: Dynamique du carbone organique dans des petits bassins versants tropicaux: exemple de la Guadeloupe, PhD thesis, Paris 7, https://www.theses.fr/2010PA077068 (last access: May 2022), 2010.

Lloret, E., Dessert, C., Pastor, L., Lajeunesse, E., Crispi, O., Gaillardet, J., and Benedetti, M. F.: Dynamic of particulate and dissolved organic carbon in small volcanic mountainous tropical watersheds, Chem. Geol., 351, 229–244, https://doi.org/10.1016/j.chemgeo.2013.05.023, 2013.

Lloret, E., Dessert, C., Buss, H. L., Chaduteau, C., Huon, S., Alberic, P., and Benedetti, M. F.: Sources of dissolved organic carbon in small volcanic mountainous tropical rivers, examples from Guadeloupe (French West Indies), Geoderma, 282, 129–138, 2016.

Minella, J. P., Merten, G. H., Reichert, J. M., and Clarke, R. T.: Estimating suspended sediment concentrations from turbidity measurements and the calibration problem, Hydrol. Process., 22, 1819–1830, 2008.

Misset, C., Recking, A., Legout, C., Poirel, A., Cazilhac, M., Esteves, M., and Bertrand, M.: An attempt to link suspended load
hysteresis patterns and sediment sources configuration in alpine catchments, J. Hydrol., 576, 72–84, 2019a.

Misset, C., Recking, A., Navratil, O., Legout, C., Poirel, A., Cazil-hac, M., Briguet, V., and Esteves, M.: Quantifying bed-related suspended load in gravel bed rivers through an analysis of the bedload-suspended load relationship, Earth Surf. Proc. Land., 44, 1722–1733, 2019b.

Misset, C., Recking, A., Legout, C., Bakker, M., Gimbert, F., Geay, T., and Zanker, S.: Using continuous turbidity and seismic measurements to unravel sediment provenance and interaction between suspended and bedload transport in an Alpine catchment, Geophys. Res. Lett., 48, e2020GL090696, https://doi.org/10.1029/2020GL090696, 2021.

Newville, M., Stensitzki, T., Allen, D. B., Rawlik, M., Ingargi-ola, A., and Nelson, A.: LMFIT: Non-linear least-square mini-mization and curve-fitting for Python, Astrophysics Source Code Library, ascl-1606, https://lmfit.github.io/lmfit-py/ (last access: November 2020), 2016.

Nistor, C. J. and Church, M.: Suspended sediment transport regime in a debris-flow gully on Vancouver Island, British Columbia, Hydrol. Process., 19, 861–885, 2005.

Orwin, J. F. and Smart, C.: Short-term spatial and temporal patterns of suspended sediment transfer in proglacial channels, Small River Glacier, Canada, Hydrol. Process., 18, 1521–1542, 2004.

OZCAR: https://in-situ.theia-land.fr/ (last access: 4 May 2023), 2023.

Park, J. and Hunt, J. R.: Coupling fine particle and bedload transport in gravel-bedded streams, J. Hydrol., 552, 532–543, 2017.

Perks, M., Owen, G., Benskin, C. M. H., Jonczyk, J., Deasy, C., Burke, S., Reaney, S., and Haygarth, P. M.: Dominant mechanisms for the delivery of fine sediment and phosphorus to fluvial networks draining grassland dominated headwater catchments, Sci. Total Environ., 523, 178–190, 2015.

Phillips, C. B., Dallmann, J. D., Jerolmack, D. J., and Packman, A. I.: Fine-Particle Deposition, Retention, and Resuspension Within a Sand-Bedded Stream Are Determined by Streambed Morphodynamics, Water Resour. Res., 55, 10303–10318, 2019.

Rad, S., Louvat, P., Gorge, C., Gaillardet, J., and Allègre, C. J.: River dissolved and solid loads in the Lesser Antilles: new insight into basalt weathering processes, J. Geochem. Explor., 88, 308–312, 2006.

Ren, J. and Packman, A. I.: Changes in fine sed-iment size distributions due to interactions with streambed sediments, Sediment. Geol., 202, 529–537, https://doi.org/10.1016/j.sedgeo.2007.03.021, 2007.

Renard, K. G., Foster, R., Weesies, G., and Porter, G.: C “RUSLE: Revised universal soil loss equation”, J. Soil Water Conserv., 46, 30–33, 1991.

Renard, K. G., Laflen, J., Foster, G., and McCool, D.: The revised universal soil loss equation, in: Soil erosion research methods, Routledge, 105–126, ISBN 9780203739358, 2017.

Roque-Bernard, A., Lucas, A., and Lajeunesse, E.: Model for Trans-port of fine particles in a steep relief river (Version 1), Zenodo [code], https://doi.org/10.5281/zenodo.7893930, 2023.

Rouse, H.: Experiments on the Mechanics of Sediment Suspension, in: Proceedings of the Fifth International Congress for Applied Mechanics, John Wiley & Sons, New York, 550–554, 1939.

Samper, A., Quidelleur, X., Lahitte, P., and Mollex, D.: Timing of effusive volcanism and collapse events within an oceanic

arc island: Basse-Terre, Guadeloupe archipelago (Lesser Antilles Arc), Earth Planet. Sc. Lett., 258, 175–191, 2007.

Sheldon, R., Prakash, A., and Sutcliffe Jr., W.: The size distribution of particles in the ocean, Limnol. Oceanogr., 17, 327–340, 1972.

Shields, A.: Application of similarity principles and turbulence research to bed-load movement, California Institute of Tech-nology, Pasadena, CA, https://resolver.caltech.edu/CaltechKHR: HydroLabpub167 (last access: 4 May 2023), 1936.

Smith, H. G. and Dragovich, D.: Interpreting sediment delivery pro-cesses using suspended sediment-discharge hysteresis patterns from nested upland catchments, south-eastern Australia, Hydrol. Process., 23, 2415–2426, 2009.

Summerfield, M. and Hulton, N.: Natural controls of fluvial denuda-tion rates in major world drainage basins, J. Geophys. Res.-Solid, 99, 13871–13883, 1994.

Suttle, K. B., Power, M. E., Levine, J. M., and McNeely, C.: How fine sediment in riverbeds impairs growth and survival of juvenile salmonids, Ecol. Appl., 14, 969–974, 2004.

Svivitski, J. P., Peckham, S. D., Hilberman, R., and Mulder, T.: Pre-dicting the terrestrial flux of sediment to the global ocean: a plan-etary perspective, Sediment. Geol., 162, 5–24, 2003.

Turowski, J. M., Rickenmann, D., and Dadson, S. J.: The partition-ing of the total sediment load of a river into suspended load and bedload: a review of empirical data, Sedimentology, 57, 1126–1146, 2010.

Van Rijn, L. C.: Sediment Transport, Part II: Suspended Load Transport, J. Hydraul. Eng., 110, 1613–1641, https://doi.org/10.1061/(ASCE)0733-9429(1984)110:11(1613), 1984.

Van Rijn, L. C.: Mathematical modeling of suspended sediment in nonuniform flows, J. Hydraul. Eng., 112, 433–455, 1986.

Van Rijn, L. C.: Unified view of sediment transport by currents and waves. II: Suspended transport, J. Hydraul. Eng., 133, 668–689, 2007.

Vercruysse, K. and Grabowski, R. C.: Temporal variation in sus-pended sediment transport: linking sediment sources and hydro-meteorological drivers, Earth Surf. Proc. Land., 44, 2587–2599, https://doi.org/10.1002/esp.4682, 2019.

Virtanen, P., Gommers, R., Oliphant, T. E., Haberland, M., Reddy, T., Cournapeau, D., Burovski, E., Peterson, P., Weckesser, W., Bright, J., van der Walt, S. J., Brett, M., Wilson, J., Millman, K. J., Mayorov, N., Nelson, A. R. J., Jones, E., Kern, R., Larson, E., Carey, C. J., Polat, I., Feng, Y., Moore, E. W., VanderPlas, J., Laxalde, D., Perktold, J., Cimrman, R., Henriksen, I., Quintero, E. A., Harris, C. R., Archibald, A. M., Ribeiro, A. H., Pedregosa, F., van Mulbregt, P., and SciPy 1.0 Contributors: SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python, Nat. Meth., 17, 261–272, https://doi.org/10.1038/s41592-019-0686-2, 2020.

White, S.: Sediment yield prediction and mod-eling, Encyclopedia of hydrological sciences, https://doi.org/10.1002/0470848944.hsa089, 2006.

Wilcock, P., Pitlick, J., and Cui, Y.: Sediment transport primer: es-timating bed-material transport in gravel-bed rivers, Gen. Tech. Rep. RMRS-GTR-226, US Department of Agriculture, Forest Service, Rocky Mountain Research Station, Fort Collins, CO, p. 78, https://doi.org/10.2737/RMRS-GTR-226, 2009.
Williams, G. P.: Sediment concentration versus water discharge during single hydrologic events in rivers, J. Hydrol., 111, 89–106, https://doi.org/10.1016/0022-1694(89)90254-0, 1989.

Wischmeier, W. H. and Smith, D. D.: Predicting rainfall erosion losses: a guide to conservation planning, Department of Agriculture, Science and Education Administration, 537 pp., https://naldc.nal.usda.gov/download/CAT79706928/PDF (last access: January 2022), 1978.

Wright, S. and Parker, G.: Flow resistance and suspended load in sand-bed rivers: simplified stratification model, J. Hydraul. Eng., 130, 796–805, 2004.

Ziegler, A. D., Benner, S. G., Tantasirin, C., Wood, S. H., Sutherland, R. A., Sidle, R. C., Jachowski, N., Nullet, M. A., Xi, L. X., Snidvongs, A., Giambelluca, T. W., and Fox, J. M.: Turbidity-based sediment monitoring in northern Thailand: Hysteresis, variability, and uncertainty, J. Hydrol., 519, 2020–2039, https://doi.org/10.1016/j.jhydrol.2014.09.010, 2014.

Zimmermann, A., Francke, T., and Elsenbeer, H.: Forests and erosion: Insights from a study of suspended-sediment dynamics in an overland flow-prone rainforest catchment, J. Hydrol., 428, 170–181, 2012.