River Suspended-Sand Flux Computation With
Uncertainty Estimation Using Water Samples And
High-Resolution Adcp Measurements

Jessica Marggraf, Guillaume Dramais, Jérôme Le Coz, Blaise Calmel, Benoît Camenen, David J. Topping, William Santini, Gilles Pierrefeu, François Lauters

► To cite this version:

Jessica Marggraf, Guillaume Dramais, Jérôme Le Coz, Blaise Calmel, Benoît Camenen, et al. River Suspended-Sand Flux Computation With Uncertainty Estimation Using Water Samples And High-Resolution Adcp Measurements. Earth Surface Dynamics, 2024, 12, pp.1243-1266. 10.5194/esurf-12-1243-2024. insu-04831461

HAL Id: insu-04831461

https://insu.hal.science/insu-04831461v1

Submitted on 11 Dec 2024

**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>

# River suspended-sand flux computation with uncertainty estimation using water samples and high-resolution ADCP measurements

Jessica Marggraf¹,*, Guillaume Dramais¹,*, Jérôme Le Coz¹, Blaise Calmel¹, Benoît Camenen¹, David J. Topping², William Santini³, Gilles Pierrefeu⁴, and François Lauters⁵

¹RiverLy, INRAE, 5 Rue de la Doua, Villeurbanne, 69100, France

²U.S. Geological Survey, Southwest Biological Science Center, Grand Canyon Monitoring and Research Center, 2255 N. Gemini Drive, Flagstaff, AZ 86001, USA

³IRD-GET, Institut de Recherche pour le Développement, Laboratoire GET (IRD, CNRS, UPS, CNES), Toulouse, France

⁴CACOH, CNR, 4 Rue de Chalon-sur-Saône, Lyon, 69007, France

⁵Service Etudes Eau Environnement, EDF, 134 Chemin de l'étang, Saint Martin Le Vinoux, 38950, France
★These authors contributed equally to this work.

**Correspondence:** Jessica Marggraf (jmarggra@umn.edu) and Guillaume Dramais (guillaume.dramais@inrae.fr)

Received: 18 October 2023 – Discussion started: 11 January 2024
Revised: 20 May 2024 – Accepted: 20 August 2024 – Published: 5 November 2024

**Abstract.** Measuring suspended-sand fluxes in rivers remains a scientific challenge due to their high spatial and temporal variability. To capture the vertical and lateral gradients of concentration in the cross-section, measurements with point samples are performed. However, the uncertainty related to these measurements is rarely evaluated, as few studies of the major sources of error exist. Therefore, the aim of this study is to develop a method to determine the cross-sectional sand flux and estimate its uncertainty. This SDC (for sand discharge computing) method combines suspended-sand concentrations from point samples with ADCP (acoustic Doppler current profiler) high-resolution depth and velocity measurements. The MAP (for multitransect averaged profile) method allows obtaining an average of several ADCP transects on a regular grid, including the unmeasured areas. The suspended-sand concentrations are integrated vertically by fitting a theoretical exponential suspended-sand profile to the data using Bayesian modeling. The lateral integration is based on the water depth as a proxy for the local bed shear stress to evaluate the bed concentration and sediment diffusion along the river cross-section. The estimation of uncertainty combines ISO standards and semi-empirical methods with a Bayesian approach to estimate the uncertainty due to the vertical integration. The new method is applied to data collected in four rivers under various hydro-sedimentary conditions: the Colorado, Rhône, Isère, and Amazon rivers, with computed flux uncertainties ranging between 18% and 32%. The relative difference between the suspended-sand flux in 21 cases calculated with the proposed SDC method compared to the International Organization for Standardization (ISO) 4363 standard method ranges between -40% and +23%. This method that comes with a flexible, open-source code is the first to propose an applicable uncertainty estimation that could be adapted to other flux computation methods.
# 1 Introduction

The determination of suspended-sediment load is required to estimate sediment dynamics and budgets for river restoration and monitoring, river engineering, and flood risk evaluation (Kondolf et al., 2014). Measuring and monitoring sediment loads and the associated uncertainties within a catchment constitute a major practical issue for hydrologists and river managers (Hoffmann et al., 2010). Even though suspended-sand transport is a key driver of the river evolution (Kondolf, 1997), it remains difficult to measure its concentration due to its temporal and spatial variability in the cross-section, which may account for several orders of magnitude (Armijos et al., 2017). In contrast, fine suspended sediments (< 63 µm) are relatively homogeneous throughout the cross-section, often reaching concentration variations of up to an order of magnitude (Wren et al., 2000).

The total suspended-sediment flux through a cross-section, $\Phi_{total}$ (kg s$^{-1}$), is defined as the mass of suspended sediment passing through a river cross-section per unit time:

$$ \Phi_{total} = \int_{y_{lb}}^{y_{rb}} q_{ss}(y)dy = \int_{y_{lb}}^{y_{rb}} \int_{z_a}^{h} c(y, z)v(y, z)dzdy, \quad (1) $$

where y and z are the lateral and vertical coordinates; $q_{ss}$ is the suspended-sediment discharge per vertical coordinate; $y_{lb}$ and $y_{rb}$ are the left and right boundaries of the cross-section; $z_a$ is the reference level for suspension at the top of the bedload layer, generally assumed to be the riverbed elevation; h is the water elevation; and $c(z)$ and $v(z)$ are the time-averaged suspended-sediment concentration and the water velocity perpendicular to the cross-section, respectively.

Suspended-sediment sampling and computing techniques have been developed over decades (Porterfield, 1972; Starosolsky and Rakoczi, 1981; International Organization for Standardization, 2002). Typically, these methods are based on physical water sampling to determine the suspended-sediment concentration throughout the cross-section using samples taken at different locations throughout the river cross-section. Samples may be taken following the depth-integrating method, where several nearly complete verticals at different distances from the bank are sampled, or the point sampling method, where samples are collected at different discrete water depths and distances from the bank. Different methods were proposed to estimate the suspended-sand flux through the cross-section (Lupker et al., 2011; Shah-Fairbank and Julien, 2015; Santini et al., 2019). The International Organization for Standardization (ISO) standard method (International Organization for Standardization, 2002) consists of computing the velocity-weighted cross-sectional mean concentration by combining physical samples with simultaneous velocity measurements. In this method, the cross-section is divided into $N_{seg}$ segments and for each increment $l$, the water discharge $Q_l$ and depth-

averaged velocity-weighted concentration $C_l$ are evaluated:

$$ \Phi_{total} = \sum_{l=1}^{N_{seg}} Q_l C_l. \quad (2) $$

This method is derived directly from the velocity-area method for the measurement of water discharge using current meters (International Organization for Standardization, 2009). Even though these classical discrete methods are widely accepted, they are time- and cost-consuming and sometimes difficult to deploy (Camenen et al., 2023). As they are limited to a few points or depth-averaged samples at a limited number of locations, they are characterized by a low spatial and temporal resolution. Also, the method is directly based on depth-integrated sampling with no possibility to interpolate and extrapolate results from sampled verticals to the whole cross-section. Some surrogate technologies (e.g., optical and acoustical methods) have been proposed to measure sediment properties and suspended-sand flux with a better spatial and temporal resolution (Wren et al., 2000; Gray and Gartner, 2010).

Acoustic methods using acoustic Doppler current profilers (ADCPs) have become well-established in streamflow monitoring and provide faster, safer, and more accurate acquisition of stream velocities, discharges, and depths than older current meter methods (Oberg and Mueller, 2007). In a measurement transect, data are acquired on a grid with fixed or variable cell height and many vertical ensembles. For a valid discharge measurement, several cross-sectional transects are typically acquired and processed to obtain information on discharge and velocity. Different post-processing tools have been developed such as the Velocity Mapping Toolbox (VMT) (Parsons et al., 2013) for the analysis and visualization of cross-sectional velocity data collected along multiple ADCP transects. Other examples are the discharge reviews QRev (Mueller, 2016) and QRevInt (Lennermark and Hauet, 2022) developed by an international group, which are applied to ensure discharge measurement reliability and to quantify the uncertainty in the discharge measurement (Despax et al., 2023). While these methods provide temporal averages of each ADCP ensemble, other approaches combine the ADCP measurements spatially close to each other, notably far from the instrument, where the distance between the beams is high (Vermeulen et al., 2014).

In combination with sediment sampling, ADCP measurements of flow velocity and depth can be used to compute the cross-sectional suspended-sand flux (Bouchez et al., 2011; Vauchel et al., 2017). ADCP measurements provide an increased spatial resolution throughout the cross-section compared to point velocity measurements using current meters or rating-curve estimates of the total cross-sectional discharge (Oberg and Mueller, 2007). Moreover, the acoustic backscatter measured by an ADCP, or an acoustic backscattering system (ABS), may be used to improve the spatial integration of the concentration in the cross-section. Indeed,
the acoustic backscatter can be inverted and used to measure the suspended-sand concentration (e.g., Topping and Wright, 2016; Venditti et al., 2016; Szupiany et al., 2019; Vergne et al., 2020). Several software tools have been developed to process ADCP data for estimating suspended-sand flux (Boldt, 2015; Dominguez Ruben et al., 2020) or using backscatter inversions. However, acoustic inversion techniques require many physical samples for calibration and are affected by acoustic modeling issues (Vergne et al., 2023).

Informed decisions related to sediment monitoring require reliable estimates of the uncertainty of flux measurements. However, the evaluation of the uncertainty is a difficult task because of the complexity of these measurements, mainly due to the temporal and spatial variations of the sediment concentration. Measurement uncertainty is the expression of the statistical dispersion of the values attributed to a measured quantity (JCGM, 2008). Identifying error sources and estimating uncertainty components for suspended-sand measurements have been addressed in many old reports and papers from the FISP (Federal Inter-Agency Sedimentation Project) (FISP, 1941, 1952; Colby, 1964; Guy and Norman, 1970) and more recently by the U.S. Geological Survey (USGS, Topping et al., 2011; Sabol and Topping, 2013) and others (Gitto et al., 2017). Moreover, the International Organization for Standardization (2002) standard proposes a framework to estimate the errors and uncertainty in the mean cross-sectional suspended-sand concentration determined by a point sampling method. It identifies several sources of error of random and systematic nature. These errors are notably related to the lateral integration of the concentration between the sampling verticals, thus depending on the number of verticals. Another important error source is the vertical integration of the concentration between the sampling points of a vertical, which depends on the number of sampling points along a vertical. Another error is related to the sampling time and consists of the sample's representativity of the natural fluctuations of the concentration due to turbulence. Additional error sources originate from the sampler type and the laboratory analysis. The uncertainty related to each of these error sources is estimated using a high number of samples, the average of which is considered to be the approximate true value. This value is taken as a reference and the uncertainties originating from the different error sources are estimated based on the deviation to the reference. The respective uncertainty is then determined by the difference between the measured value at a given location and the approximate true value. Even though several sources of error are addressed, the International Organization for Standardization (2002) method contains several defects in theoretical and practical respects. First, this standard is not in agreement with the framework proposed in the *Guide to the expression of uncertainty in measurement* (GUM; JCGM, 2008) defining the uncertainty propagation method, notably concerning the notations and the computation of an approximate true value. Second, the large number of additional sam-

ples required for the uncertainty analysis and taken under stable hydro-sedimentary conditions is hardly applicable in most environments due to the quick variability of hydro-sedimentary processes and technical difficulties. For example, to estimate the uncertainty due to the number of verticals (i.e., lateral integration), 15 to 20 verticals with seven point samples each are required for sections less than 100 m wide. The time needed to conduct this kind of survey is practically impossible given the temporal variability of the processes studied. Besides the International Organization for Standardization (2002) method, no other method proposes a framework addressing all commonly identified sources of error. However, some authors tried to evaluate the main sources of uncertainty. Concerning lateral integration, Colby (1964) noticed that the sand flux varies approximately (within a factor $k_1$) with the third power of the mean velocity $\bar{v}$ for a constant grain size distribution and temperature, as well as velocities ranging between 0.6 and 1.5 m s⁻¹: $\Phi_{\text{total}} = k_1 \bar{v}^3$. Based on these observations and making the required conversions, he stated that the variability of sand concentration at different sampling verticals should be closely related to the variability of $\bar{v}^2/h$, the ratio of the squared mean velocity $\bar{v}$ to the total sampled depth $h$. To ensure comparability among different sampling sections and streams, the $\bar{v}^2/h$ index, also called $\xi$, may be used:

$$ \xi = \frac{\max(\bar{v}_l^2/h_l)}{\bar{v}_z^2/h_z}, \qquad (3) $$

with $\bar{v}_l$ and $h_l$ as the depth-averaged velocity and water depth for each increment $l$ and $\bar{v}_z$ and $\bar{h}_z$ as the depth-averaged velocity and water-depth mean values for the cross-section, respectively. Based on this concept of variability, Guy and Norman (1970) prepared a nomograph that indicates the number of sampling verticals required for a desired maximum acceptable relative standard uncertainty as a function of the percentage of sand $p_s$ and $\xi$.

The issue of the appropriate sampling time and the associated time averaging has been the subject of several studies. Topping et al. (2011) analyzed the temporal variability in sediment concentration among point samples and estimated the associated uncertainties for depth-integrated measurements. In addition, Gitto et al. (2017) concluded that a 9 to 12 min sampling time was required to get a representative point sample because of the temporal variability in sediment concentration.

Common methods to estimate the mean cross-sectional suspended-sand concentration and its uncertainty are subject to various limitations such as the interpolation and extrapolation of the suspended-sand concentration towards the riverbed and riverbank or the impractical feasibility of the International Organization for Standardization (2002) uncertainty method. Therefore, the first aim of this study is to introduce a method for computing the total suspended-sand flux with a high spatial resolution. That method combines
point samples with ADCP measurements using a physically based understanding of suspended-sediment transport processes. A second aim is to provide a method to estimate the uncertainty related to this suspended-sand flux computation. Therefore, the uncertainties related to several sources of error such as discharge, lateral and vertical integration, and determined point concentrations are estimated and combined following the GUM framework.

First, the proposed sand discharge computing (SDC) method to estimate the suspended-sand concentration and the uncertainty in the sand flux is presented (Sect. 2). Then, this method is applied to four rivers across the world with different flow and sediment characteristics (Sect. 3). Finally, the methodology and presented results are discussed and further developments are suggested in Sect. 4.

We use the following notations:

- $v$ is the water velocity perpendicular to the cross-section;

- $u$ is the absolute standard uncertainty – that is, the standard deviation of the probability distribution of errors, “absolute” meaning expressed in the physical unit of the measurement (e.g., in $m^3 s^{-1}$ for discharge);

- $u'$ is the relative standard uncertainty, “relative” meaning expressed in percent of the measurement;

- $U = ku$ is the absolute expanded uncertainty, with $k$ a coverage factor taken as $k = 2$, which corresponds to a 95% probability interval if the distribution of errors is Gaussian;

- $U'$ is the relative expanded uncertainty expressed in percent of the measurement result.

# 2 Method

## 2.1 General method

The proposed sand discharge computing (SDC) method combines ADCP velocity and discharge measurements performed in multiple transects with suspended-sand concentrations obtained by point samples distributed throughout the cross-section (Fig. 1). The point sand concentrations are interpolated vertically and laterally in the cross-section using a physically based method and Bayesian modeling. In the second part of the SDC method, several sources of error are estimated using novel and literature-based approaches and combined to estimate the uncertainty in suspended-sand flux measurements.

The SDC method focuses on suspended-sand flux measurements in simple river cross-sections without tidal effects or strong secondary currents. In the case of strong secondary flow cells causing deviations in the suspended-sand

concentration from the dominant flow equilibrium, these deviations should be evaluated separately. Choosing an appropriate measurement site is essential to obtain reliable sand flux and uncertainty estimates (Edwards and Glysson, 1999).

## 2.2 Physically based method to integrate the concentration in the cross-section

### 2.2.1 Creating an ADCP multitransect averaged profile (MAP)

The multitransect averaged profile (MAP) method was developed to perform the SDC method proposed in this study (Fig. 1) and implemented in QRevInt (Lennermark and Hauet, 2022) in 2023 (see the Supplement for a detailed description). A typical ADCP discharge measurement consists of the average of individual discharge measurements from successive ADCP transects. The MAP method includes each ADCP transect to generate an averaged transect profile from the bottom to the water surface including the unmeasured areas. This method developed with Python 3 (Van Rossum and Drake, 2009) is based on the QRevInt measurement output. QRevInt provides quality analysis and quality control, which allows having clean input data. For each ADCP discharge measurement, composed of several transects, one averaged MAP profile is computed with a regular grid for the whole cross-section (including the unmeasured areas). Major differences to similar tools like to the Velocity Mapping Toolbox (VMT) (Parsons et al., 2013) are that the MAP method can compute an average profile in the absence of a GPS positioning and that it allows the extrapolation of areas unmeasured by the ADCP.

First, the MAP method defines a straight average cross-section on the selected transects of the measurement. Then, these transects are projected using an orthogonal translation on the average cross-section. At this point, each transect is interpolated on the cross-section grid (Fig. 2b). The width and height of cells can be defined by the user. Once each transect is defined on the grid, MAP overlays them to average velocity and depth in each cell of the averaged cross-section (Fig. 2c-e). Finally, velocities are extrapolated to unmeasured areas (Fig. 2f). For edge extrapolation, banks are divided into the same shape meshes as the middle streamflow. Mean primary velocity $\bar{v}_p$ (Rozovskii, 1957) on each edge's vertical is computed according to a power law:

$$ \bar{v}_p = \bar{v}_{0/l} \left( \frac{x}{L_{\text{edge}}} \right)^{\frac{1}{m_{\text{edge}}}}, \quad (4) $$

where $x$ is the distance of the vertical from the start of the bank, $L_{\text{edge}}$ is the length of the edge, $m_{\text{edge}}$ is an edge-shape exponent (2.41 for a triangular edge, 10 for a rectangular edge; edge extrapolation is not computed otherwise), and $\bar{v}_{0/l}$ is the mean primary velocity from the closest measured vertical. Then, primary velocity on each edge's mesh $v_p$ is
<figure><img src="image_3.png" /><figcaption>Figure 1. Workflow of the proposed sand discharge computing (SDC) method to estimate the suspended-sand flux distribution and the uncertainty in suspended-sand flux through a cross-section.</figcaption></figure>

computed following a power law with the QRevInt extrapolation exponent $m_{extrap}$:

$$v_p = \bar{v}_p \frac{m_{extrap} + 1}{m_{extrap}} \left( \frac{z_{cell}}{d} \right)^{\frac{1}{m_{extrap}}}, \quad (5)$$

where $z_{cell}$ is the depth to centerline of mesh and $d$ the depth on the vertical of the edge. Edge extrapolation of secondary velocity uses a linear law between the closest vertical and the edge. Vertical velocity on the edge follows the distribution of vertical velocities on the closest vertical. MAP thus generates a complete averaged profile with homogeneous cell sizes. Each cell contains information on its distance to the left bank and its depth and velocity components. Primary and secondary velocities are then transformed into stream-wise and cross-stream velocities in order to compute discharge.

## 2.2.2 Point sampling and laboratory analysis

In addition to the ADCP multitransects, point sampling distributed in the cross-section is performed (Fig. 1). Each suspended-sediment measurement follows the point sampling method and contains $m$ verticals (typically three to seven) with $N_{sam}$ samples per vertical (typically four to five). Two types of samplers, a Niskin water-trap-type sampler (instantaneous non-time-averaged sample; Filizola et al., 2009) and isokinetic samplers (US P-06; Spicer, 2019), were used in different rivers deployed from boats (Colorado, Rhône, and Amazon rivers) or cable cars (Colorado and Isère rivers). The target depth is set with a graduated tag line when deploying from the boat or using the depth information from the reel on the cable car. The sampler is equipped with a pressure sensor for post facto verification (for the Rhône, Isère, and Amazon rivers). An electrical valve allows the US P-06

samplers to collect a sample at the desired depth and for the desired sampling duration. For instantaneous samples, taken with the Niskin sampler, a traveler is sent down the rope to close the sampler. For the analysis of the suspended-sediment concentration of each sample, the standard procedure (option C) of the American Society for Testing and Materials (ASTM) (ASTM D3977, 2007) is applied, which consists of separating fine sediments and sand by wet-sieving prior to filtration of the fine sediments.

## 2.2.3 Vertical suspended-sand concentration profiles

A physically based method is applied to assign concentrations to individual cells $(i, j)$ in each sampling vertical $l$ (Fig. 1). It uses a theoretical vertical suspended-sand concentration profile estimated using a Bayesian approach to interpolate and extrapolate the sand concentrations vertically from point samples. The exponential vertical concentration profile proposed by Camenen and Larson (2008), based on a constant vertical diffusivity $\epsilon_v$ throughout the water column, is defined as

$$C(z) = C_R \exp(\alpha z), \quad (6)$$

where $C(z)$ (kg m$^{-3}$ or g L$^{-1}$) is the sediment concentration at elevation $z$ above the bed, $\alpha$ is the vertical gradient in a logarithmic scale, and $C_R$ is the bottom reference sediment concentration. Prior to the implementation of the Camenen and Larson (2008) profile, we compared empirically fitted Camenen and Larson (2008) and Rouse (1937) profiles. The differences in sand fluxes are lower than the uncertainty in sand flux of the respective measurement (see the Supplement). Moreover, the profile of Camenen and Larson (2008) is of practical interest in that it does not require arbitrarily
<figure><img src="image_4.png" /><figcaption>Figure 2. Workflow of the transect averaging procedure deployed by the multitransect averaged profile (MAP) method (based on Parsons et al., 2013).</figcaption></figure>

defining a reference level. Consequently, we used only the Camenen and Larson (2008) profile for the Bayesian modeling in the SDC method. However, other vertical concentration profiles may be used as well and included in the toolbox.

To estimate the concentration profile $C(z)$, the derived depth-averaged concentration, and its uncertainty $u'_p$ due to vertical integration (Sect. 2.3.5), the Bayesian Modeling (BaM!) method is applied (Mansanarez et al., 2019). The BaM! method is based on Bayesian inference, which allows the computation of the posterior probability of a model's parameters from their prior probability and from observations. The model can then be applied to predict the distribution of a new, unobserved data point. The posterior distribution of the parameters is computed using Bayes' theorem, and a large number (> 10 000) of realizations are sampled using an adaptive block Metropolis Markov chain Monte Carlo (MCMC) sampler (Renard et al., 2006) varying the parameters $\alpha$ and $\ln(C(z))$. A linear model is applied using logarithmic concentrations in milligrams per liter based on Eq. (6):

$$ \ln(C(z)) = \ln(C_R) - |\alpha|z. \qquad (7) $$

The first 5000 realizations are burned as a warm-up period, and then the last 5000 realizations are decimated to decrease the correlation in the results. The MaxPost profile $\ln(C_{n_0}(z))$, the best-fitting profile, is computed with the realization of parameters $n_0$ that maximizes the posterior distribution and is used for calculating the sand discharge. From this MaxPost profile, the concentration $C_{i,j}$ in each cell $(i, j)$ along vertical $l$ can be determined. The MaxPost parameters $\ln(C_{R,n_0})$ and $\alpha_{n_0}$ are retained and used for the lateral interpolation (Sect. 2.2.4).

BaM! requires the definition of the prior distribution of the equation parameters that are $\ln(C_R)$ and $\alpha$ here. Both $|\alpha|$ (to ensure increasing concentrations with depth) and $C_R$ are strictly positive, and therefore they are assumed to follow lognormal distributions with parameters $\mu$ and $\sigma$. Consequently, $\ln(C_R)$ is assumed to follow a Gaussian distribution. The parameters $\mu_\alpha$ and $\sigma_\alpha$ describing the prior distribution of $\alpha$ are the mean and standard deviation of the variable's natural logarithm, respectively. The expected values of $\alpha$ and $C_R$ are evaluated based on local hydro-sedimentary parameters (Camenen and Larson, 2008), which are determined using ADCP depth and velocity, as well as bedload measurements. The expected value of $C_R$ is calculated using the expression of Camenen and Larson (2008), which is a function of the sedimentological diameter $D_*$, the Shields parameter $\theta$, and the critical bed shear stress $\theta_{cr}$:
$$C_R = 1.5 \ 10^3 \theta \exp(-0.2 D_*) \exp\left(-4.5 \frac{\theta_{cr}}{\theta}\right). \quad (8)$$

The sedimentological diameter, or dimensionless grain size, $D_*$ is calculated as

$$D_* = \bar{D}_{50} \left( \frac{(s-1)g}{\nu^2} \right)^{1/3}, \quad (9)$$

where $\bar{D}_{50}$ is the median diameter of the sand suspension averaged over the analyzed vertical, $s = 2.65$ is the relative sediment density, $g = 9.81 \text{ m s}^{-2}$ is the acceleration due to gravity, and $\nu \approx 10^{-6} \text{ m}^2 \text{ s}^{-1}$ is the kinematic viscosity of water. The expected value of the prior distribution is then converted to $\ln(C_R)$. This reference concentration $C_R$ differs significantly from the reference concentration for a Rouse profile, where the reference concentration is sensitive to the (more or less arbitrary) choice of the reference level, adding some complexity to evaluating the priors. Similarly, the expected value of $\alpha$ can be determined as (Camenen and Larson, 2008)

$$\alpha = -\frac{6 w_s}{\sigma_t \kappa v_* h}, \quad (10)$$

where $w_s$ (m s$^{-1}$) is the settling velocity estimated following the formula of Soulsby and Whitehouse (1997), $\sigma_t$ is the turbulent Schmidt number set equal to 1 as a first approximation, $\kappa = 0.41$ is the von Kármán constant, $v_*$ is the total shear velocity, and $h$ is the water depth determined using ADCP measurements.

Defining the prior of $\alpha$ as lognormal ensures that it remains negative under all hydro-sedimentary conditions. This implies that the concentration decreases as a function of $z$ away from the bed, thereby corresponding to Rouse mechanics for suspended-sediment computing (Rouse, 1937). The concentration of the finest sizes in suspension may increase away from the bed when the concentration of suspended sediment is relatively high due to the “squeezing” effect or density stratification (Hunt, 1969; McLean, 1992), leading to possible positive $\alpha$ values. We neglect these effects since we focus on sand with relatively low concentrations. Grain size information on vertical $l$ is necessary to determine Eq. (8) and the settling velocity $w_s$ and thus to estimate both $\ln(C_R)$ and $\alpha$. In the case that they are not available, no prior parameters are defined and the model is fitted by BaM! on the observations only.

The second parameter $\sigma_\alpha$ of the lognormal distribution of $\alpha$ can be estimated by uncertainty propagation equations established following the Guide to the expression of Uncertainty in Measurement (GUM; JCGM, 2008). It is estimated based on the relative uncertainty $u'_\alpha$ of $\alpha$, supposing $\sigma_\alpha = u'_\alpha$. This approximation works well for small values ($< 0.5$) of $\sigma$ of the respective lognormal distribution. The parameter $\sigma_\alpha = u'_\alpha$ is estimated by propagation from Eq. (10):

$$u'_\alpha = \sqrt{u'_{ws}^2 + u'_{\sigma t}^2 + u'_{\kappa}^2 + u'_{v_*}^2 + u'_{h}^2}. \quad (11)$$

Only a few studies have evaluated these uncertainties. Since additional experimental measurements are beyond the scope of this article, we define, based on the literature, the uncertainty in the settling velocity $u'_{ws} = 5\%$ (Camenen, 2007), the uncertainty in the turbulent Schmidt number $u'_{\sigma t} = 20\%$ (Gualtieri et al., 2017), $u'_{\kappa} = 0$ (theoretical value with negligible variations; Smart, 2022), the uncertainty in the shear velocity $u'_{v_*} = 5\%$ (Perret et al., 2023), and the uncertainty in the elevation of the sampled point within the water column $u'_{h} = 5\%$ (Dramais, 2020). With these values, we obtain $u'_\alpha = \sqrt{0.0475} \approx 21.8\%$.

The second parameter $\sigma_{\ln(C_R)}$, the standard deviation of $\ln(C_R)$, could also be determined by an uncertainty propagation derived from the data reduction equation of $C_R$ (Eq. 8). However, it has been shown that the highest uncertainty is related to the structural uncertainty of the formula of $C_R$ itself, not to its parametric uncertainty (Camenen et al., 2014). Indeed, the dataset used to establish the semi-empirical formula of Camenen and Larson (2008) is characterized by a large scatter, with differences of about 50% between the measured and predicted concentration (other formulas also come with large structural uncertainty). Consequently, it is assumed that $\sigma_{\ln(C_R)} = u'_{C_R} = 50\%$.

### 2.2.4 Lateral interpolation

The lateral interpolation of the suspended-sand concentration to calculate $C_{i,j}$ in every cell of the MAP grid is based on a physical approach using the water depth as an index (Fig. 1). Following Camenen and Larson (2008), $C_R$ is set proportional to the local bed shear stress, which can be assumed to be proportional to the water depth $h$ if the friction slope is constant throughout the river cross-section (Khdashenas and Paquier, 1999; Camenen et al., 2011). Thus, the ratio $C_{R,j}/h_j$ of the reference concentration $C_R$ and the water depth $h$ for each column $j$ in the MAP grid is estimated through linear interpolation along the cross-section. As a first approximation, $\alpha_j$ is assumed to be independent of the local bed shear stress, since it is mostly influenced by large-scale turbulence structures (Van Rijn, 1984). $\alpha_j$ varies linearly with horizontal distance between two adjacent sampling verticals (where $\alpha$ was estimated from the concentration profiles fitted to the samples) and remains constant between the first and last sampling vertical and the edge of the cross-section.

### 2.2.5 Determination of concentration $C_{i,j}$ in each MAP cell $(i, j)$

The proposed SDC method is based on the discretization of the river cross-section by a regular grid fitted on the ADCP
data (MAP grid) composed of $N_j$ columns and $N_i$ depth cells (Fig. 1). The general idea is to assign a concentration and discharge to each cell so that a flux per cell can be obtained after multiplication. The total cross-sectional sand flux $\Phi_{total}$ is calculated by summing up the suspended-sand fluxes per cell:

$$ \Phi_{\text{total}} = \sum_{i=1}^{N_i} \sum_{j=1}^{N_j} \Phi_{i,j}, \qquad (12) $$

where $\Phi_{i,j}$ (kg s⁻¹) is the suspended-sand flux through one MAP cell $i, j$. The suspended-sand flux $\Phi_{i,j}$ can be calculated as

$$ \Phi_{i,j} = C_{i,j} Q_{i,j} = C_{i,j} u_{i,j} w_j h_{\text{cell},i,j}, \qquad (13) $$

where $C_{i,j}$ (kg m⁻³ or g L⁻¹) and $Q_{i,j}$ (m³ s⁻¹) are the suspended-sand concentration and liquid discharge through each cell ($i, j$), respectively; $u_{i,j}$ (m s⁻¹) is the normal velocity component; $w_j$ (m) is the width; and $h_{\text{cell},i,j}$ (m) is the height of the ith vertical cell in the jth column in the MAP grid. The discharge $Q_{i,j}$ through each cell is determined using the novel MAP method based on QRevInt (Lennermark and Hauet, 2022) and the suspended-sand concentration $C_{i,j}$ is determined following the novel physically based SDC method.

The parameters $\alpha$ and $C_R$ are evaluated for each MAP cell ($i, j$) applying the presented vertical and lateral integration. The suspended-sand concentration in each cell in the MAP grid is thus evaluated as

$$ C_{i,j} = \frac{z_{i,j} + h_{\text{cell},i,j/2}}{h} \int_{z_{i,j} - h_{\text{cell},i,j/2}}^{z_{i,j} + h_{\text{cell},i,j/2}} C_{R,j} \exp(\alpha_j z) dz. \qquad (14) $$

## 2.3 Estimation of the uncertainty in measurements of suspended-sand flux through a cross-section

### 2.3.1 General method

The uncertainty $U'_{\Phi}$ in measurements of the suspended-sand flux through a cross-section is based on the calculation of suspended-sand flux (i.e., Eq. 2, Fig. 1). Therefore, the flux $\Phi$ is the product of discharge $Q$ and mean cross-sectional concentration $\bar{C}$: $\Phi = Q \times \bar{C}$. Thus, $U'_{\Phi}$ can be separated into a factor related to discharge $U'_{Q}$ and one related to the concentration $U'_{C}$:

$$ U'_{\Phi} = \sqrt{U'_{Q}^{2} + U'_{C}^{2}}. \qquad (15) $$

Equation (15) is based on the hypothesis that the errors in discharge and concentration are independent; otherwise, the term has to include the associated covariances. Such an assumption that the errors are independent, however, appears to be reasonable. First, because discharge and concentration

are measured independently, the discharge is measured using ADCPs and the concentration is determined using sampling and laboratory analyses. Second, the error sources in discharge and concentration are significantly different. Velocity lateral interpolation errors are negligible due to the high spatial resolution of ADCP measurements, whereas concentration lateral interpolation errors are large. We agree that these two error components are physically correlated, but this is not a problem as the first one is negligible.

To approximate the uncertainty $U'_{\Phi}$ in the suspended-sand flux through a cross-section, uncertainties $U'_{Q}$ and $U'_{C}$ in the discharge and concentration both have to be determined (Fig. 3).

The uncertainty $U'_{Q}$ in multiple-transect ADCP discharge measurements is calculated following the OURSIN method (Despax et al., 2023) as implemented in the open-source software QRevInt (Lennermark and Hauet, 2022). $U'_{C}$ is the combination of several uncertainty components (see Fig. 3) listed in Table 1 and detailed afterwards.

### 2.3.2 Uncertainty $U'_{C}$ in the mean cross-sectional suspended-sand concentration

The uncertainty $u'_{C}$ in the mean cross-sectional suspended-sand concentration is calculated as

$$ u'_{C} = u'_{\text{sys,C}} + u'_{m} + \sum_{l=1}^{m} \frac{\Phi_{l}^{2}}{\Phi_{l}^{2}} u'_{p,l}, \qquad (16) $$

where $u'_{\text{sys,C}}$ is the uncertainty due to systematic errors in the concentration, $u'_{m}$ is the uncertainty due to the lateral integration based on the number $m$ of verticals, $u'_{p,l}$ is the total uncertainty due to the vertical integration estimated for each vertical $l$ (see Fig. 3), and $\Phi_{l}$ is the suspended-sediment flux through vertical $l$.

### 2.3.3 Uncertainty due to systematic sources of error

$$ u'_{\text{sys,C}} $$

Following the International Organization for Standardization (2002) method, the uncertainty due to systematic sources of error $u_{\text{sys,C}}$ is expressed as

$$ u'_{\text{sys,C}} = u'_{\text{sys,m}} + u'_{\text{sys,p}} + u'_{\text{sys,lab}} + u'_{\text{sys,sampler}}, \qquad (17) $$

where $u'_{\text{sys,m}}$ is the uncertainty due to the systematic error of the flux computation scheme, $u'_{\text{sys,p}}$ is the uncertainty due to the systematic error of the vertical integration, $u'_{\text{sys,lab}}$ is the uncertainty due to the systematic error of the laboratory analysis, and $u'_{\text{sys,sampler}}$ is the uncertainty due to the systematic error of the sampler type since the underlying errors are assumed to be systematic. These terms, detailed in International Organization for Standardization (2002), remain constant, independently of the increasing number of sampling points or verticals (see Table 1).
<figure><img src="image_5.png" /><figcaption>Figure 3. Flow diagram of the developed approach to estimate the uncertainty $U'_{\Phi}$ in the cross-sectional suspended-sand flux.</figcaption></figure>

<table><caption>Table 1. Uncertainty components used to compute the uncertainty $U'_{C}$ in the mean cross-sectional suspended-sand concentration.</caption><thead><tr><th>Error sources</th><th>Notation<br>(standard uncertainty)</th><th>Nature</th><th>Estimation<br>method</th></tr></thead><tbody><tr><td>Lateral interpolation</td><td>u<sub>m</sub></td><td>Systematic</td><td>Nomograph from Guy and Norman (1970)</td></tr><tr><td>Vertical interpolation</td><td>u<sub>p</sub></td><td>Random</td><td>Bayesian approach (Sect. 2.3.5): total uncertainty combining the next two components</td></tr><tr><td>Vertical interpolation</td><td>u<sub>p.param</sub></td><td>Random</td><td>Bayesian approach (Sect. 2.3.5): parametric uncertainty</td></tr><tr><td>Vertical interpolation</td><td>u<sub>p.struc</sub></td><td>Random</td><td>Bayesian approach (Sect. 2.3.5): structural uncertainty</td></tr><tr><td>Concentration measurements</td><td>u<sub>meas</sub></td><td>Random</td><td>Formula (Eq. 19) combining the following three components</td></tr><tr><td>Sampling type</td><td>u<sub>sampler</sub></td><td>Random</td><td>Fixed values (8%* and 16%)</td></tr><tr><td>Laboratory analysis</td><td>u<sub>lab</sub></td><td>Random</td><td>Formula (Eq. 20, Gordon, 2000)</td></tr><tr><td>Natural fluctuations of concentration</td><td>u<sub>nf</sub></td><td>Random</td><td>Repeated measures experiments (Sect. 2.3.8)</td></tr><tr><td>Systematic concentration errors</td><td>u<sub>sys.C</sub></td><td>Systematic</td><td>Formula (Eq. 17*) combining the four next components (3.54%)</td></tr><tr><td>Systematic lateral integration errors</td><td>u<sub>sys.m</sub></td><td>Systematic</td><td>Fixed value (1.5%*)</td></tr><tr><td>Systematic vertical integration errors</td><td>u<sub>sys.p</sub></td><td>Systematic</td><td>Fixed value (2%*)</td></tr><tr><td>Systematic laboratory analysis errors</td><td>u<sub>sys.lab</sub></td><td>Systematic</td><td>Fixed value (2%*)</td></tr><tr><td>Systematic sampler type errors</td><td>u<sub>sys.sampler</sub></td><td>Systematic</td><td>Fixed value (1.5%*)</td></tr></tbody></table>

### 2.3.4 Uncertainty $u'_m$ due to lateral integration

To facilitate the application compared to the standardized approach (International Organization for Standardization, 2002), the uncertainty $u'_m$ due to lateral integration is estimated based on Eq. (3) and the nomograph published by Guy and Norman (1970):

$$u'_{m} = 0.4 p_{s} (1.43 \xi - 1.37) m^{-0.7}, \qquad (18)$$

with $p_s$ the percentage of sand in the suspension, $\xi$ (see Eq. 3) the $\bar{v}^2/h$ index (Colby, 1964), and $m$ the number of verticals. If only the suspended-sand flux through a cross-section is measured or is of principal interest, as in our study, the percentage of sand in the equation should be assumed to be 100%, neglecting the influence of the fine sediment flux. This does not signify that no suspended fine sediments are present but allows taking into account the considerable lateral gradients of the sand suspension. It also avoids underestimating the uncertainty due to the lateral integration because the uncertainty $u'_m$ for the same sediment discharge measurement (same $m$ and $\xi$) is higher when assuming $p_s = 1$ than including fine sediments. This approach is applied in our study, although the lateral interpolation applied differs slightly, as it is based on the water depth $h$ and the parameters $C_R$ and $\alpha$ of the vertical profiles. However, this approach is assumed to be consistent with our modified lateral interpolation.

### 2.3.5 Uncertainty $u'_p$ due to vertical integration

The uncertainty $u'_p$ is determined for each vertical $l$ from the distribution of vertically integrated concentrations computed from the profiles estimated by the Bayesian approach described in Sect. 2.2.3. This uncertainty accounts for the uncertainty $u'_{meas}$ (estimated in Sect. 2.3.6) in point concentrations taken as observational data in the Bayesian inference. The integration of the previously obtained vertical concentration profiles $\ln(C_n(z))$ (Sect. 2.2.3) allows the determination of the parametric uncertainty $u'_{p.param}$. However, computing the total uncertainty $u'_p$ due to vertical integration requires the inclusion of structural errors at the elevation of the sampling points prior to the vertical integration. These structural errors are representative of the residuals between the point measurements and the exponential profiles. The structural uncertainty can be estimated from the total uncertainty $u'_p$ and the parametric uncertainty $u'_{p.param}$: $u'_{p.struc}$ as $u'_{p.struc}^2 = u'_{p}^2 - u'_{p.param}^2$.

The parametric uncertainty $u'_{p.param}$ can be determined from the distribution of concentration profiles $\ln(C_n(z))$ computed in Sect. 2.2.3 (Fig. 4a). Each of these $n$ profiles is converted to $C_n(z)$ (Fig. 4b) and linearly interpolated by applying a trapezoidal integration to determine its depth-averaged concentration $\bar{C}_n$, which is converted to $\ln(\bar{C}_n)$. Application of the entire procedure for all simulations $n$ then yields

<footnote>* Value originating from International Organization for Standardization (2002).</footnote>
a distribution of depth-averaged concentrations $ln(C_n)$. The mean value of this distribution is $ln(C)$ and the standard deviation is the uncertainty $u'_{p,\text{param}}$ based on the assumption $\sigma = u'_{p,\text{param}}$ (Fig. 4c).

To determine the structural error, the prior distribution of its standard deviation is defined as lognormally distributed with $\mu = 0$ and $\sigma = 1$ in the BaM! method. For every sampling point at the elevation $z$, a normally distributed error with mean zero and standard deviation $u'_{\text{meas}}$ is defined. An error is then drawn from this distribution and added to the estimated concentration $ln(C_n(z))$ for every simulation $n$ to obtain a modified vertical profile $ln(C_{\text{mod},n}(z))$ (Fig. 4d). In the next step, the same procedure as for the estimation of $u'_{p,\text{param}}$ is applied: conversion of $ln(C_{\text{mod},n}(z))$ to $C_{\text{mod},n}(z)$, vertical averaging to obtain $\bar{C}_{\text{mod},n}$, and conversion to $ln(\bar{C}_{\text{mod},n})$. The mean value of the resulting distribution is the mean depth-averaged concentration $ln(\bar{C}_{\text{mod}})$ and its standard deviation is the total uncertainty $u'_{p}$ due to vertical integration based on the assumption $\sigma = u'_{p}$ (Fig. 4f).

### 2.3.6 Uncertainty $u'_{\text{meas}}$ in point concentrations

As point concentration errors are accounted for in the Bayesian analysis of vertical concentration profiles, the uncertainty $u'_{\text{meas}}$ in point concentrations is already included in the uncertainty $u'_{p}$ due to vertical integration. Therefore, in contrast to the International Organization for Standardization (2002) method, $u'_{\text{meas}}$ does not explicitly appear in Eq. (16). The uncertainty $u'_{\text{meas}}$ is calculated as

$$u'_{\text{meas}} = \sqrt{u'^2_{\text{sampler}} + u'^2_{\text{lab}} + u'^2_{\text{nf}}}, \quad (19)$$

where $u'_{\text{sampler}}$ is the uncertainty due to the sampler type, $u'_{\text{lab}}$ is the uncertainty due to the laboratory analysis, and $u'_{\text{nf}}$ is the uncertainty due to natural fluctuations in sediment concentration arising from turbulence (Fig. 3).

Uncertainty $u'_{\text{sampler}}$ is due to the sampler type. Even though several comparisons have been conducted, the distribution of random errors related to a specific sampler type is difficult to assess. For example, a review of the values of $u'_{\text{sampler}}$ used in different studies is provided by Dramais (2020). In this study, the value suggested in the International Organization for Standardization (2002) standard is used for isokinetic samplers such as the US P-06: $u'_{\text{sampler}} = 8\%$. To account for the greater uncertainty arising from non-isokinetic sampling, this uncertainty is arbitrarily doubled for non-isokinetic samplers: $u'_{\text{sampler}} = 16\%$.

### 2.3.7 Uncertainty $u'_{\text{lab}}$ due to laboratory analysis

Many studies have estimated the random uncertainty related to the measurement of (fine) sediment concentration in the laboratory (e.g., by filtration). The ISO method estimates an uncertainty of 1.5% due to the random error and an uncertainty of 2% due to the systematic error (International Organization for Standardization, 2002). Based on an intercomparison study of different laboratories, Gordon et al. (2000) determined a standard uncertainty for the fine and sand fractions separately. We use the approach of Gordon et al. (2000) at the 68% confidence level and a given sand concentration $C$ (g L$^{-1}$) in the analyzed sample:

$$u'_{\text{lab}} = 1.091 C^{-0.5}. \quad (20)$$

### 2.3.8 Uncertainty $u'_{\text{nf}}$ due to natural fluctuations

To approximate the uncertainty $u'_{\text{nf}}$ due to the natural fluctuations in concentration and grain size in the point samples arising from turbulence, a simplified method, similar to the ISO method or the “at-a-point error” (APE) (Topping et al., 2011), is applied. To this end, several points are repeated at different hydro-sedimentary conditions with a time difference of less than 1 h between the first and last sample, and the suspended-sand concentration $C_i$ is calculated for each sample. One sampling point is repeated three to nine times and the mean sediment concentration $\bar{C}_{\text{rep}}$ of the respective set of measurements is determined. Based on the nomenclature of the International Organization for Standardization (2002), this mean concentration $\bar{C}_{\text{rep}}$ per set can be understood as the “approximate true value”. The relative standard deviation $u'_{\text{rep}}$ for each set of $N_{\text{rep}}$ repetitions is then calculated following International Organization for Standardization (2002):

$$u'_{\text{rep}} = \sqrt{\frac{\sum_{i=1}^{N_{\text{rep}}} \left( \frac{C_i}{\bar{C}_{\text{rep}}} - 1 \right)^2}{N_{\text{rep}} - 1}}. \quad (21)$$

Performing this calculation for all repetitions, the relative uncertainty for each set of repetitions $u'_{\text{rep}}$ can be plotted versus the mean concentration $\bar{C}_{\text{rep}}$ per set (Fig. 5).

The number of sets of repetitions and tested hydro-sedimentary conditions within this study is limited compared to the variety of sampling conditions. In the best case, these measurements should be conducted on every sampling campaign; however, in reality, this is hardly possible using the presented measurement techniques. Using the same sampling protocol, the sampling campaign with many additional samples for the uncertainty estimation would take so long that the variation in river discharge would become too great. Therefore, a constant uncertainty $u'_{\text{nf}} = 14.24\%$ is determined based on these results and applied to all point measurements, which corresponds to the median of all tested relative uncertainty $u'_{\text{rep}}$. The enlarged uncertainty of $U'_{\text{nf}} = 28.47\%$ at a 95% confidence interval roughly corresponds to the estimations made by Gitto et al. (2017) in the Fraser River (they found a 3% to 33% uncertainty range for individual 30's samples). Furthermore, it should be noted that only a small range of hydro-sedimentary conditions at a given sampling
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Workflow for the estimation of the uncertainty $u'_p$ due to vertical integration, including the estimation of the parametric uncertainty $u'_{p,\text{param}}$ in panels (a), (b), and (c) and of the total uncertainty $u'_p$ in panels (d), (e), and (f). (a, d) Vertical concentration profiles $\ln(C_n(z))$ and $\ln(C_{\text{mod},n}(z))$, respectively, sampled through Bayesian interference and including the structural error and two exemplary profiles in (d). (b, e) Vertical concentration profiles $C_n(z) = \exp(\ln(C_n(z)))$ and $C_{\text{mod},n}(z) = \exp(\ln(C_{\text{mod},n}(z)))$, respectively, with the MaxPost profile $C_{n0}(z)$ in (b). (c, f) Histograms of depth-averaged concentrations $\ln(C_n)$ and $\ln(C_{\text{mod},n})$ with the mean depth-averaged concentrations $\ln(C)$ and $\ln(C_{\text{mod}})$ as well as the standard deviations $u'_{p,\text{param}}$ and $u'_p$, respectively.</figcaption></figure>

location is sampled by this empirical approach. The uncertainties are probably higher than estimated here and may also be grain-size-dependent (Topping et al., 2011).

# 3 Application

## 3.1 Survey sites

The proposed SDC method was applied to four datasets from different rivers around the world. Each dataset includes suspended-sediment measurements following the protocol presented above and ADCP data.

### 3.1.1 Rhône River

The Rhône River is one of the major rivers of Europe, heading from the Rhône Glacier in the Alps and running through western Switzerland and southeastern France. Mostly a gravel-bed river, it is the largest silt and clay contributor to the Mediterranean sea (Delile et al., 2020). The

presented measurements were conducted near the gaging station (V3000015) at Lyon Perrache (WGS84 coordinates: 45.742344, 4.826738), France, where the Rhône River drains a catchment of about 20 300 km² with a mean annual discharge of about 600 m³ s⁻¹ (Dramais, 2020).

### 3.1.2 Isère River

The Isère River is an Alpine river and the largest tributary of the Rhône River by suspended-sediment flux (Poulier et al., 2019). At Grenoble, France (WGS84 coordinates: 45.197747, 5.768566), where the measurements were conducted (gaging station W1410010) (Némery et al., 2013), the mean annual discharge is about 180 m³ s⁻¹ with a catchment area of 5700 km².

### 3.1.3 Colorado River

The Colorado River is one of the most iconic rivers in the western United States. The measurements took place at the
<figure><img src="image_7.png" /><figcaption>**Figure 5.** The relative uncertainty $u'_{\text{rep}}$ and mean concentration $\bar{C}_{\text{rep}}$ for each set of repetitions (these repetitions include data from the Isère, Colorado, Toutle, and Cowlitz rivers; Spicer, 2019). The relative uncertainty $u'_{\text{nf}}$ and the relative expanded uncertainty $U'_{\text{nf}}$ due to natural fluctuations correspond to the average and the average multiplied by a coverage factor of $k = 2$ of all tested relative uncertainties $u'_{\text{rep}}$.</figcaption></figure>

U.S. Geological Survey (USGS) site above the Little Colorado River near the Desert View (WGS84 coordinates: 36.203484, -111.800917), Arizona, gaging station at River Mile 61. This station (number 09383100; U.S. Geological Survey, 2023) has a mean annual discharge of 306 m³ s⁻¹ and a catchment area of 296 000 km². Suspended sediments have been monitored for a long time in this area (Topping et al., 2021).

### 3.1.4 Amazon River

The Amazon River basin exceeds 6 000 000 km² in area. The Amazon River is the largest river in the world by discharge. The Manacapuru gauging station (14100000) is part of the Critical Zone Observatory HyBAm (Hydrology of the Amazon Basin) and is operated by the French National Research Institute for Sustainable Development (IRD), the Brazilian National Agency (ANA), and the Brazilian Geological Service (CPRM). This station has been used for more than 40 years by the Brazilian national hydrometric network to provide data on the Amazon (Solimões) River (WGS84 coordinates: -3.324377, -60.561183). At this station, the Amazon River watershed is approximately 2 × 10⁶ km² and average water discharge is about 103 000 m³ s⁻¹ (Filizola et al., 2009).

Those four survey sites, with various geomorphological conditions (see Table 2), were sampled according to the above-described ADCP measurement and point sampling procedures. The Isère River was sampled with an isokinetic US P-06 sampler and the other rivers with a water-trap-type sampler (Niskin).

## 3.2 Vertical suspended-sand concentration and flux profiles

Measured suspended-sand point concentrations are fitted with an exponential profile to extrapolate the concentrations to the unmeasured parts of the water column and also interpolate between points. There is substantial vertical and lateral variability in suspended-sand concentration at all study sites (Fig. 6). Indeed, different vertical gradients $\alpha$ and/or reference concentrations $C_R$ are observed among the measured verticals. The highest sand concentrations and largest gradients, with a difference of up to 3 orders of magnitude between the bottom and surface concentrations, are observed in the Amazon River (see Fig. 6d). This indicates a possible effect of the water depth on the concentration gradient through the vertical diffusion coefficient (inversely proportional to the water depth; Eq. 10) and/or due to dunes (whose height is typically proportional to the water depth and can therefore enhance the shear velocity). In contrast, the concentrations at the other sites range between 0.01 and 0.5 g L⁻¹. In the various surveys, sand concentration gradients are associated with particle size gradients, with coarser particles closer to the riverbed. The measured concentrations vary strongly at some verticals so that they do often not correspond to the fitted vertical concentration profiles, not even when taking the uncertainty $U'_{\text{meas}}$ in the point concentrations into account. This uncertainty usually varies for the presented samples between 20 % and 25 % at a 95 % confidence interval.

Vertical profiles of suspended-sand flux (Fig. 7) are determined by multiplying the suspended-sand concentration in each cell in the MAP grid with the discharge in the same cell. Similarly, the point suspended-sand fluxes are the product of the point concentration and the discharge of the surrounding cell in the MAP grid. Consequently, decreasing fluxes close to the bed, as expected based on theory, are hardly or not at all visible for most verticals sampled in the Isère (Fig. 7b) and the Colorado (Fig. 7c). Large differences between the point fluxes and the profiles notably result from poorly fitted vertical concentration profiles, e.g., vertical 34 in the Isère River and vertical 1650 in the Amazon River. In other words, when the point concentrations do not follow an exponential profile, there are large differences between point fluxes and profiles.

## 3.3 Suspended-sand flux through a cross-section

The suspended-sand concentration in each MAP cell (i, j) is calculated by applying the lateral interpolation and extrapolation of the profile coefficients $C_R$ and $\alpha$ (Sect. 2.2.3). The spatial view of the cross-sections highlights the distribution of the suspended-sand concentration (Fig. 8). Different layers in some measurements appear due to the vertical and horizontal resolution of the ADCP data, i.e., the size of the MAP cells. As the vertical integration is based on the water depth, the lateral interpolation of the profile coefficients produces
<table><caption>**Table 2.** Hydraulic and morphological characteristics of the survey sites. $Q_w$ is the total ADCP water discharge including extrapolation areas, $Q_{Meas.}$ is the ADCP-measured water discharge (excluding extrapolations areas), $\bar{v}$ is the mean cross-sectional velocity, the aspect ratio is the river width divided by mean depth, $D_{50}$ is the median diameter of suspended sand, $p_s$ the percentage of sand in the suspension, $\xi$ is an index relating stream velocity and depth, and $U'_\phi$ is the expanded uncertainty in measurements of suspended-sand flux at the 95% confidence level.</caption><thead><tr><th>River and location</th><th>Survey date</th><th>Q<sub>w</sub><br>(m<sup>3</sup> s<sup>-1</sup>)</th><th>Q<sub>Meas.</sub>/Q<sub>w</sub><br>(%)</th><th>ν&#773;<br>(m s<sup>-1</sup>)</th><th>W<br>(m)</th><th>H<br>(m)</th><th>Aspect<br>ratio (-)</th><th>Sand D<sub>50</sub><br>(µm)</th><th>p<sub>s</sub><br>(%)</th><th>ξ<br>(-)</th><th>U'<sub>Φ</sub><br>(%)</th></tr></thead><tbody><tr><td>Rhône at Lyon Perrache (France)</td><td>22 Jan 2018</td><td>2000</td><td>81</td><td>1.2</td><td>170</td><td>12</td><td>14</td><td>100–300</td><td>26.9</td><td>1.07</td><td>18.6</td></tr><tr><td>Isère at Grenoble (France)</td><td>6 Apr 2022</td><td>120</td><td>52</td><td>1.1</td><td>70</td><td>3</td><td>23</td><td>90–290</td><td>46.6</td><td>1.23</td><td>19.1</td></tr><tr><td>Colorado River Mile 61 (USA)</td><td>19 Feb 2019</td><td>370</td><td>65</td><td>1.1</td><td>100</td><td>5</td><td>20</td><td>100–130</td><td>68.2</td><td>1.6</td><td>31.1</td></tr><tr><td>Amazon at Manacapuru (Brazil)</td><td>19 Apr 2012</td><td>144 000</td><td>83</td><td>1.6</td><td>3400</td><td>43</td><td>79</td><td>180</td><td>76.6</td><td>1.33</td><td>26.6</td></tr></tbody></table>

<figure><img src="image_8.png" /><figcaption>**Figure 6.** Measured sand concentrations with uncertainty $U'_\text{meas}$ and exponential fits (represented as colored lines) for each sampling vertical using Bayesian modeling for the (a) Rhône River at Lyon Perrache, (b) Isère River at the Grenoble campus, (c) Colorado River at River Mile 61, and (d) Amazon River at Manacapuru. The colors correspond to the sampling verticals, which are indicated by their distance in meters from the left bank, called abscissa.</figcaption></figure>

high concentrations near the bed, especially when there are large water depth variations and when vertical measurements are made on the deepest parts. This is clearly observed in Fig. 8d between the central and right sampling vertical on the Amazon.

The mean cross-sectional suspended-sand concentrations $\bar{C}_{SDC}$ and fluxes $\Phi_{SDC}$ computed with the SDC method are compared to the ISO method (International Organization for Standardization, 2002) using the relative differences $\epsilon_C = (C_{SDC} - C_{ISO})/C_{ISO}$ and $\epsilon_\phi$, respectively. The results for the suspended-sand concentration are in close agreement between the two methods, with $\epsilon_C$ ranging between -2% and 3.5% for three examples (Table 3), whereas a significant concentration difference is observed between the two

methods for the Colorado computations (-15.8%). The most likely hypothesis to explain this difference is that the surface sample at the middle of the transect has a relatively low sand concentration (abscissa 50, Fig. 6c). This low concentration heavily influenced the fit of the vertical profile and reduced the flux in this part of the cross-section, which is the place of the most intense flow. This highlights one of the limitations of the method when only a few points are used for suspended-sand flux computation in the cross-section.

## 3.4 Suspended-sand uncertainty evaluation

The total uncertainty and the contribution of each error source to total variance are evaluated for the four measure-
<figure><img src="image_9.png" /></figure>

**Figure 7.** Suspended-sand fluxes interpolated in each MAP cell along the sampling verticals (lines) and measured at the sampling points for the (a) Rhône River at Lyon Perrache, (b) Isère at Grenoble, (c) Colorado River at River Mile 61, and (d) Amazon River at Manacapuru. The colors correspond to the sampling verticals, which are indicated by their distance in meters from the left bank, called abscissa.

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

**Figure 8.** Suspended-sand concentrations calculated with the SDC method in each cell of the MAP grid throughout the cross-section and location of the sampling points of the (a) Rhône River at Lyon Perrache, (b) the Isère River at the Grenoble campus, (c) the Colorado at River Mile 61, and (d) the Amazon River at Manacapuru (d).
<table><caption>**Table 3.** Mean cross-sectional suspended-sand concentrations and total fluxes for the four presented measurements using the ISO method and SDC method.</caption><thead><tr><th>Study site</th><th>C&#773;<sub>ISO</sub> (g L<sup>-1</sup>)</th><th>C&#773;<sub>SDC</sub> (g L<sup>-1</sup>)</th><th>&epsilon;<sub>C</sub> (%)</th><th>&Phi;<sub>ISO</sub> (kg s<sup>-1</sup>)</th><th>&Phi;<sub>SDC</sub> (kg s<sup>-1</sup>)</th><th>&epsilon;<sub>&Phi;</sub> (%)</th></tr></thead><tbody><tr><td>Rhône River</td><td>0.051</td><td>0.050</td><td>-2.0</td><td>102</td><td>99.3</td><td>-2.8</td></tr><tr><td>Isère River</td><td>0.085</td><td>0.088</td><td>3.5</td><td>9.2</td><td>9.5</td><td>2.9</td></tr><tr><td>Colorado River</td><td>0.120</td><td>0.101</td><td>-15.8</td><td>48.1</td><td>40.5</td><td>-16</td></tr><tr><td>Amazon River</td><td>0.304</td><td>0.298</td><td>-2.0</td><td>40284</td><td>39438</td><td>-2.1</td></tr></tbody></table>

ments (Fig. 9). The absolute uncertainty at a 95 % confidence interval with a coverage factor of $k = 2$ ranges between 19 % and 31 %. The main uncertainty components are the uncertainty $U'_p$ due to vertical integration and the uncertainty $U'_m$ due to lateral integration. The uncertainty $U'_p$ is displayed as its two components, the parametric uncertainty $u'_{p,\text{param}}$ and the structural uncertainty $u'_{p,\text{struc}}$. The parametric uncertainty $u'_{p,\text{param}}$ is determined by the information from the priors and from the sampling points used to calibrate the model, particularly their number, distribution along the vertical, and uncertainty $u'_{\text{meas}}$. Increasing the number of samples and decreasing the uncertainty $u'_{\text{meas}}$ would decrease this uncertainty $u'_{p,\text{param}}$. In contrast, the structural uncertainty $u'_{p,\text{struc}}$ is estimated from the residuals of the fit of the model to the calibration points. The farther they are from the fitted vertical concentration profile and the lower their uncertainty, the greater the uncertainty $u'_{p,\text{struc}}$.

The uncertainty $U'_m$ is estimated from index $\xi$, representing the lateral homogeneity of the cross-section in terms of depth and discharge with the percentage of sand in the suspension and the number of sampling verticals. However, as only sand concentrations are considered here, the percentage of sand was set equal to 100 %. At relatively high $\xi$ values, such as on the Colorado River, a larger number of verticals would have been required to decrease the uncertainty $U'_m$, whereas this uncertainty is relatively low in a more uniform river like the Rhône River with a low $\xi$ value. In that latter case, three verticals are sufficient to describe the lateral distribution of the concentration in the cross-section. The measurements on the Isère and Amazon rivers are both characterized by similar $\xi$. The large number of sampling verticals in the Isère River (seven) leads to a low uncertainty $U'_m$, whereas the uncertainty and small number of sampling verticals in the Amazon River (three) lead to high uncertainty $U'_m$. Nevertheless, it has to be taken into account that the lateral integration is based on the water depth, whereas the calculation of the uncertainty $U'_m$ is based on $\xi$. These two different ways to conduct the lateral integration of the suspended-sand concentration and the uncertainty estimation may affect the results. The contributions of the uncertainties $U'_Q$ in the liquid discharge and the systematic uncertainties $u'_{\text{sys}}$ to the uncertainty $U'_\phi$ in the suspended-sand flux are typically low.

<figure><img src="image_11.png" /><figcaption>**Figure 9.** Total expanded uncertainty $U'_\phi$ in suspended-sand flux and the relative contributions of uncertainty components to the total variance on the (a) Rhône River at Lyon Perrache, (b) on the Isère River at Grenoble, (c) on the Colorado at River Mile 61, and (d) on the Amazon River at Manacapuru. The five components are the expanded uncertainty $U'_Q$ in the liquid discharge, the expanded uncertainty $U'_{\text{sys}}$ due to systematic sources of error, the expanded uncertainty $U'_m$ due to the lateral integration of the concentration, the expanded parametric uncertainty $U'_{\text{param}}$ due to the vertical integration, and the expanded structural uncertainty $U'_{\text{struc}}$ due to the vertical integration.</figcaption></figure>

# 4 Discussion

## 4.1 A physically based method

The novel SDC method offers a number of advances for cross-sectional sediment flux measurement, especially the physically based integration of concentration.

For the vertical integration of concentration, the first step of the method is based on the fit of an exponential profile. The uncertainty related to the fit of this profile is estimated. As discussed in Camenen and Larson (2007), very similar results would have been obtained by fitting a Rouse profile (see
the Supplement). However, another more detailed approach (Hunt, 1969; McLean, 1992) may provide a better fit and thus modify the results. Consequently, various theoretical approaches including the effects of suspended-sediment stratification or size distribution of sediment may be integrated into the toolbox, allowing us to choose the best-fitting semi-empirical model to the sampling conditions. Currently, only the exponential profile and the Rouse profile (not presented here but included in the code) are available in the toolbox.

The goodness of the fit, and thus the structural uncertainty, depends not only on the chosen profile, but also on the number and position of the sampling points. Dramais (2020) applied the same Bayesian approach and showed that the points measured close to the bed have a great influence on the reference concentration $C_R$ value and on the slope $\alpha$ of the fitted exponential profile. In contrast, subsurface samples may in some cases bias the fit of the exponential profile. Interestingly, the best compromise was observed when sampling points are positioned close to the bed and for five samples per profile.

Depth-integrated measurements could be a solution to avoid errors due to the vertical interpolation, which occur with fitting point samples. However, with this protocol, lateral extrapolations and extrapolations in the unmeasured parts of the cross-section are not possible with a physical base. Depth-integrated samples could then be associated with a larger uncertainty.

The lateral integration of the concentration may be improved as well, even though the SDC method presents an advance by using the water depth as a proxy for the shear stress. The current method leads to artifacts of relatively high concentrations close to the riverbed within extrapolated areas (Fig. 8d). Moreover, the use of the parameters of the vertical concentration profile for the lateral integration increases the importance of a well-fitted vertical profile; otherwise, large uncertainties may result.

A major advance in the vertical integration, and particularly the lateral integration, may be made using the acoustic backscatter measured by the ADCP. One option is the use of the backscattered signal intensity as a proxy for suspended-sediment concentration, facilitating the interpolation and extrapolation between the sampling points. Moreover, several studies (Bouchez et al., 2011; Venditti et al., 2016; Szupiany et al., 2019) and some commercial software (e.g., ASET; Dominguez Ruben et al., 2020) focus on the vertical moving-boat backscatter inversion to gain information on suspended sediment. For the inversion of single-frequency applications as the ADCP, strong assumptions or calibrations are necessary to correctly estimate the concentration and grain size of silt–clay and sand-sized sediment (Vergne et al., 2023). Additionally, the issue of unmeasured areas close to the riverbed, surface, and banks persists and requires the extrapolation of the estimated concentrations, e.g., by applying theoretical suspended-sand transport formulas (Dominguez Ruben et al., 2020). However, when using the

ADCP backscatter for the calculation of the concentrations, the errors in discharge and concentration would probably not be independent anymore because the same measurement method is used and the uncertainty due to the spatial integration of the concentration, the major error source of the uncertainty in concentration, will decrease. Consequently, the formulation of the uncertainty estimation should be adapted, e.g., by including the covariances.

## 4.2 A high-resolution ADCP data-based method

Our study proposes a general method which uses high-resolution ADCP data from successive transects. Compared to existing multitransect averaging tools, the newly developed ADCP multitransect averaged profile (MAP) provides an average dataset including the unmeasured areas. MAP uses RDI or SonTek raw binary files and reduces the preprocessing error, as it uses data quality filters from QRevInt. The method may use either the bottom track or the GPS as a reference and the user can customize vertical and lateral dimensions of the resulting grid cells. The obtained regular grid then facilitates the further analysis steps.

One limitation of the SDC method is the vertical and lateral integration of concentrations, which can be limited by the ADCP data resolution. The vertical and lateral integration of concentrations is evaluated on each MAP grid cell. Although this approach allows the estimation of concentrations close to the bed and banks, the size of the grid cells is limited by the size of the ADCP cells. Consequently, if the ADCP spatial resolution is low, the resulting mean concentration may be affected. An example of this problem is provided by the data from the Colorado River (Fig. 10), where cells were set with a height of approximately 0.4 m (Fig. 10a) and 0.8 m (Fig. 10b). In some cases, it could be meaningful to adapt the ADCP cells size to increase the resolution of the measurement and consequently increase the resolution of the resulting cross-sectional estimation of the distribution of the sediment flux (e.g., Vermeulen et al., 2014).

## 4.3 A method open to various sampling protocols

Another advantage of the SDC method is its suitability to different point sampling protocols, with various numbers and locations of sampling points along the verticals and varying numbers of verticals. This flexibility is particularly useful when specific areas or depths are of major interest and require more detailed sampling or if the sampling points are not distributed in the cross-section following the ISO protocols. Also, it provides an estimation of the suspended-sand concentration close to the bed or banks in areas excluded by many methods such as the ISO method.

We evaluate the differences of the suspended-sand flux calculated using the proposed SDC method relative to the flux calculated using the ISO method using a larger dataset for Rhône River and Isère River that encompasses the de-
<figure><img src="image_12.png" /><figcaption>**Figure 10.** Influence of the vertical cell height on the cross-sectional suspended-sediment measurement in the Colorado River at River Mile 61, (a) with a cell height of 0.4 m and (b) with a cell height of 0.8 m.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 11.** Relative difference $\epsilon_{\Phi} = (\Phi_{SDC} - \Phi_{ISO})/\Phi_{ISO}$ as a function of the suspended-sand flux $\Phi_{ISO}$ determined using the ISO method for all four studied rivers. $\Phi_{SDC}$ is determined using the SDC method for all four studied rivers.</figcaption></figure>

tailed examples shown above. For these supplementary data, conditions were different but they were at the same locations. The relative difference $\epsilon_{\Phi}$ ranges for all four studied rivers between -40% and +23% with no clear relationship with the total flux (Fig. 11). These values are lower than the estimated uncertainties, underlining the importance of the errors associated with the sand flux measurement. Similar observations can be made when using an empirically fitted Rouse profile instead of the exponential profile fitted with Bayesian modeling (see the Supplement). In comparison, typical ADCP water discharge measurements are characterized by an uncertainty of 5%-12% (results from several repeated experiments performed in France; Despax et al., 2019).

The relative differences between the sand fluxes calculated following the SDC method compared to the ISO method may reach up to 40%, but no bias is visible. Extensive field and laboratory measurements under stable and known conditions are required to evaluate the performance of the two methods and determine the best among these two and other existing

methods. Divergences may originate from the difference in the sampling protocols. The ISO method is based on a sampling protocol with seven points per vertical, where the samples are taken at precise relative depths, while four samples are typically taken per vertical and at varying relative depths in the present measurements. In addition, the SDC method may be more accurate than the ISO method at high discharges or deep water depths when it is difficult to lower the sampler close to the riverbed. However, accurate flux references are lacking in rivers, so this assumption cannot be verified experimentally.

## 4.4 A first estimation of the suspended-sand flux uncertainty

This method combines existing and novel approaches to estimate the uncertainty in the suspended-sand flux taking all error sources and their contributions into account in the uncertainty budget, which represents a major advance over the existing ISO method. The method is easy and fast to apply and contains no empirical calculation, except for the uncertainty $u'_{nf}$ in natural concentration fluctuations. Compared to the ISO method, this considerably reduces the required sampling time and effort, facilitating its application even for monitoring purposes. The price of this is additional assumptions on the spatial distribution of concentrations based on a simple physical basis valid for gauging cross-sections with a simple flow. Our method allows interpolating and extrapolating concentrations with a physically based approach. Furthermore, the introduced Bayesian approach to compute $u'_{p}$ appears to be a promising way to analyze vertical concentration profiles and the related uncertainties.

This new approach has advantages, but certain limitations have already been identified. Here are some suggestions for correcting them or improving the code. Concerning the uncertainty estimation, major advances could be made by developing a more robust method for the estimation of the uncertainty $u'_{m}$ due to the lateral integration, the uncertainty $u'_{sampler}$ due to the sampler, and the uncertainty $u'_{nf}$ due to the natural fluctuations. The major issue about the estimation of the uncertainty $u'_{m}$ is its difference to the physical approach
followed for the lateral integration of the concentration. The determination of the $\bar{v}^2/h$ index, $\xi$, accounts for the water depth and the stream velocity, whereas only the water depth is taken into account for the lateral interpolation of the concentration. Developing a consistent method for both uses that is adapted to various channel geometries and lateral concentration gradients is needed, e.g., following approaches developed for index–velocity relations (e.g., Kästner et al., 2018). Similarly, robust methods to estimate the uncertainty $u'_{\text{sampler}}$ due to the sampler and the uncertainty $u'_{\text{nf}}$ due to natural fluctuations for different settings should be developed. Finally, potentially uncovered error sources such as the uncertainty in the vertical position and the uncertainty related to the total sampling duration (Topping et al., 2011; Gitto et al., 2017) should be estimated and integrated as well. Moreover, the estimation of the prior distributions of $\alpha$ and $\ln(C_R)$ to estimate the uncertainty $u'_{p}$ due to the vertical integration could be estimated as well. Other vertical concentration profiles may be included, and the best-fitting formula may be chosen, thereby limiting the uncertainty in the vertical integration.

## 4.5 An open-source method

A fully operational and open-source toolbox is available. This toolbox includes several options not presented in detail in this article, like the use of an empirically fitted Rouse profile. The code is relatively flexible and suitable for various conditions and protocols commonly applied in point sampling protocols. It allows the computation of suspended sand, but also the silt–clay concentration if available for various sampler types and deployment conditions. If silt–clay concentrations are available, the ratio of silt–clay concentration to sand concentration can be compared at each sampling point and also over the entire cross-section. The vertical position of the sampler may be determined as well as the transit time between the water surface and the final sampling depth if pressure sensor measurements are available. If available, the results of the flux measurements may be related to data from adjacent hydro-sedimentary gauging stations. If grain size data are available for several samples, they are visualized and a mean cross-sectional grain size distribution following International Organization for Standardization (2002) is calculated.

# 5 Conclusion

The new SDC method presented in this study allows meaningful determination of the suspended-sand flux through a river cross-section with uncertainty. Therefore, this method merges data from ADCP discharge measurements and point suspended-sediment samples. The SDC method includes a method for averaging several ADCP transects with discharge and velocity measurements on a regular grid in the entire cross-section is developed. Suspended-sand concentrations obtained by point sampling are then vertically interpolated by fitting a physically based exponential concentration profile and choosing the best fit using a Bayesian framework (BaM!). The lateral interpolation between the point samples and extrapolation in the unmeasured zones are performed on a physical basis. Both the vertical and lateral integrations allow the computation of the suspended-sand concentration for each ADCP grid cell and consequently the suspended-sand flux.

The toolbox presented in this article proposes a major advance in the estimation of the uncertainty in point suspended-sand sampling. It addresses several sources of error and integrates existing methods with novel approaches to propose an applicable framework. The main error sources are identified as $u_m$ due to lateral integration and $u_p$ due to vertical integration, thereby justifying the SDC method, which seeks to improve spatial integration in the whole cross-section.

The application of the methodology to several cross-sectional suspended-sand measurements conducted following different sampling protocols on four global rivers yields results that slightly differ from the ISO method (–15.9% to +2.9% suspended-sand flux difference). This approach can be easily used and is adaptable to different sampling cases; the only requirement is an ADCP discharge measurement including several transects and a point sample dataset. The data processing, analysis, and visualization toolbox are open-access and available online.

Future development may benefit from focusing on the incorporation of the acoustic backscatter measured by the ADCP to guide the vertical and lateral integration and on the development of more robust methods of estimating the uncertainties due to lateral integration, the sampler performance, and the natural fluctuations in concentration arising from turbulence.
Appendix A: Notations

Roman symbols

<table><tr><td>C</td><td>Suspended-sediment concentration (g L<sup>-1</sup>)</td></tr><tr><td>C&#773;<sub>ISO</sub></td><td>Mean cross-sectional sand concentration calculated using the ISO method (g L<sup>-1</sup>)</td></tr><tr><td>C<sub>R</sub></td><td>Bottom reference sediment concentration (g L<sup>-1</sup>)</td></tr><tr><td>C&#773;<sub>rep</sub></td><td>Mean concentration per repetition set (g L<sup>-1</sup>)</td></tr><tr><td>d</td><td>Depth of the vertical on the edge (m)</td></tr><tr><td>D<sub>*</sub></td><td>Sedimentological diameter (-)</td></tr><tr><td>C&#773;<sub>SDC</sub></td><td>Mean cross-sectional sand concentration calculated using the SDC method (g L<sup>-1</sup>)</td></tr><tr><td>D&#773;<sub>50</sub></td><td>Median diameter of the grain size distribution (m<sup>-1</sup>)</td></tr><tr><td>g</td><td>Acceleration due to gravity (m s<sup>-2</sup>)</td></tr><tr><td>i</td><td>Vertical MAP cell coordinate (-)</td></tr><tr><td>j</td><td>Lateral MAP cell coordinate (-)</td></tr><tr><td>h</td><td>Water depth (m)</td></tr><tr><td>h<sub>cell</sub></td><td>Cell height (m)</td></tr><tr><td>k</td><td>Coverage factor (-)</td></tr><tr><td>k<sub>1</sub></td><td>Coefficient relating the sand flux with the velocity (Colby, 1964) (-)</td></tr><tr><td>l</td><td>Increment, sampling vertical (-)</td></tr><tr><td>L<sub>edge</sub></td><td>Length of the edge in ADCP measurements (m)</td></tr><tr><td>ln(C<sub>n0</sub>(z))</td><td>MaxPost best-fitting vertical concentration profile (g L<sup>-1</sup>)</td></tr><tr><td>ln(C<sub>R,n0</sub>)</td><td>MaxPost parameter (ln(g L<sup>-1</sup>)</td></tr><tr><td>m<sub>edge</sub></td><td>Edge-shape exponent in ADCP measurements (-)</td></tr><tr><td>m<sub>extrap</sub></td><td>QRevInt extrapolation exponent (-)</td></tr><tr><td>n<sub>0</sub></td><td>Parameters α and ln(C<sub>R</sub>) used to calculate the ln(C<sub>n0</sub>(z)) MaxPost profile (-)</td></tr><tr><td>N<sub>rep</sub></td><td>Number of repetitions per set (-)</td></tr><tr><td>N<sub>sam</sub></td><td>Number of samples per vertical (-)</td></tr><tr><td>N<sub>seg</sub></td><td>Number of segments (-)</td></tr><tr><td>m</td><td>Number of verticals (-)</td></tr><tr><td>u</td><td>Absolute standard uncertainty</td></tr><tr><td>u'</td><td>Relative standard uncertainty (%)</td></tr><tr><td>u'<sub>α</sub></td><td>Uncertainty in the slope of the exponential profile α</td></tr><tr><td>u'<sub>C<sub>R</sub></sub></td><td>Uncertainty in the reference concentration C<sub>R</sub> (%)</td></tr><tr><td>u'<sub>ln(C<sub>R</sub>)</sub></td><td>Uncertainty in the logarithmic reference concentration ln(C<sub>R</sub>) (%)</td></tr><tr><td>u'<sub>h</sub></td><td>Uncertainty in the elevation of the sampled point within the water column (%)</td></tr><tr><td>u'<sub>lab</sub></td><td>Uncertainty due to laboratory analysis (%)</td></tr><tr><td>u'<sub>m</sub></td><td>Random uncertainty due to the lateral integration (%)</td></tr><tr><td>u'<sub>meas</sub></td><td>Uncertainty in the point concentration (%)</td></tr><tr><td>u'<sub>nf</sub></td><td>Random uncertainty due to natural fluctuations in the suspended-sediment concentration (%)</td></tr><tr><td>u'<sub>p</sub></td><td>Total uncertainty due to the vertical integration (%)</td></tr><tr><td>u'<sub>p.param</sub></td><td>Parametric uncertainty due to the vertical integration (%)</td></tr><tr><td>u'<sub>p.struc</sub></td><td>Structural uncertainty due to the vertical integration (%)</td></tr><tr><td>u'<sub>2</sub></td><td>Uncertainty in the discharge of multiple transects (%)</td></tr><tr><td>u'<sub>Q</sub></td><td>Relative standard deviation for each set of N<sub>rep</sub> repetitions (%)</td></tr><tr><td>u'<sub>rep</sub></td><td>Uncertainty due to the sampler type (%)</td></tr><tr><td>u'<sub>sampler</sub></td><td>Uncertainty due to systematic errors in the suspended-sediment concentration (%)</td></tr><tr><td>u'<sub>sys,C</sub></td><td>Uncertainty due to systematic errors in the suspended-sediment concentration (%)</td></tr><tr><td>u'<sub>sys,lab</sub></td><td>Uncertainty due to the systematic error of the laboratory analysis (%)</td></tr><tr><td>u'<sub>sys,m</sub></td><td>Uncertainty due to the systematic error of the flux computation scheme (%)</td></tr><tr><td>u'<sub>sys,p</sub></td><td>Uncertainty due to the systematic error of the vertical integration (%)</td></tr><tr><td>u'<sub>sys,sampler</sub></td><td>Uncertainty due to the systematic error of the used sampler type (%)</td></tr><tr><td>u'<sub>v*</sub></td><td>Uncertainty in the total shear velocity v* (%)</td></tr><tr><td>u'<sub>ws</sub></td><td>Uncertainty in the settling velocity w<sub>s</sub> (%)</td></tr></table>
Roman symbols

<table><tr><td>u'<sub>κ</sub></td><td>Uncertainty in the von Kármán constant κ (%)</td></tr><tr><td>u'<sub>σt</sub></td><td>Uncertainty in the turbulent Schmidt number σ (%)</td></tr><tr><td>U</td><td>Absolute expanded uncertainty</td></tr><tr><td>U'</td><td>Relative expanded uncertainty (%)</td></tr><tr><td>U'<sub>C</sub></td><td>Relative expanded uncertainty in measurements of the suspended-sediment concentration (%)</td></tr><tr><td>U'<sub>m</sub></td><td>Relative expanded uncertainty due to the lateral integration (%)</td></tr><tr><td>U'<sub>p</sub></td><td>Relative expanded uncertainty due to the vertical integration (%)</td></tr><tr><td>U'<sub>Φ</sub></td><td>Relative expanded uncertainty in measurements of the suspended-sediment flux (%)</td></tr><tr><td>U'<sub>Q</sub></td><td>Relative expanded uncertainty in measurements of the discharge (%)</td></tr><tr><td>U'<sub>Q2</sub></td><td>Relative expanded uncertainty in the discharge of multiple transects (%)</td></tr><tr><td>U'<sub>sys</sub></td><td>Relative expanded uncertainty arising from systematic sources of error (%)</td></tr><tr><td>p<sub>s</sub></td><td>Percentage of sand in the suspension (%)</td></tr><tr><td>Q</td><td>Water discharge (m<sup>3</sup> s<sup>-1</sup>)</td></tr><tr><td>q<sub>ss</sub></td><td>Suspended-sediment discharge per vertical (kg s<sup>-1</sup>)</td></tr><tr><td>s</td><td>Relative sediment density (–)</td></tr><tr><td>v</td><td>Water velocity perpendicular to the cross-section (m s<sup>-1</sup>)</td></tr><tr><td>v&#773;</td><td>Mean velocity (m s<sup>-1</sup>)</td></tr><tr><td>v&#773;<sub>p</sub></td><td>Mean primary velocity (m s<sup>-1</sup>)</td></tr><tr><td>v&#773;<sub>0/l</sub></td><td>Mean primary velocity from the closest measured vertical (m s<sup>-1</sup>)</td></tr><tr><td>v<sub>*</sub></td><td>Total shear velocity (m s<sup>-1</sup>)</td></tr><tr><td>w</td><td>Cell width (m)</td></tr><tr><td>w<sub>s</sub></td><td>Settling velocity (m s<sup>-1</sup>)</td></tr><tr><td>x</td><td>Distance of the vertical from the start of the bank (m)</td></tr><tr><td>X</td><td>East coordinate (MAP)</td></tr><tr><td>X<sub>proj</sub></td><td>East coordinate on the average cross-section (MAP)</td></tr><tr><td>Y</td><td>North coordinates (MAP)</td></tr><tr><td>Y<sub>proj</sub></td><td>North coordinate on the average cross-section (MAP)</td></tr><tr><td>y</td><td>Lateral coordinate</td></tr><tr><td>y<sub>lb</sub></td><td>Left boundary of the cross-section (m)</td></tr><tr><td>y<sub>rb</sub></td><td>Right boundary of the cross-section (m)</td></tr><tr><td>z</td><td>Vertical coordinate</td></tr><tr><td>z<sub>a</sub></td><td>Reference level for suspension at the top of the bedload layer (m)</td></tr><tr><td>z<sub>cell</sub></td><td>Depth to the centerline of the cell (m)</td></tr></table>

Greek symbols

<table><tr><td>α</td><td>Vertical gradient of the exponential vertical concentration profile (Camenen and Larson, 2008) (–)</td></tr><tr><td>α<sub>R,n<sub>0</sub></sub></td><td>MaxPost parameter (–)</td></tr><tr><td>ε<sub>v</sub></td><td>Vertical diffusivity (m<sup>2</sup> s<sup>-1</sup>)</td></tr><tr><td>ε<sub>C</sub></td><td>Relative difference in sand concentration determined by the SDC method compared to the ISO method (%)</td></tr><tr><td>ε<sub>Φ</sub></td><td>Relative difference in sand flux determined by the SDC method compared to the ISO method (%)</td></tr><tr><td>θ</td><td>Shields parameter (–)</td></tr><tr><td>θ<sub>cr</sub></td><td>Critical Shields parameter or critical bed shear stress (–)</td></tr><tr><td>κ</td><td>von Kármán constant (–)</td></tr><tr><td>μ</td><td>Parameter of the lognormal distribution</td></tr><tr><td>ν</td><td>Kinematic viscosity of water (m<sup>2</sup> s<sup>-1</sup>)</td></tr><tr><td>ξ</td><td>Ratio of the squared mean velocity to the total sampled depth (Colby, 1964)</td></tr><tr><td>σ</td><td>Parameter of lognormal distribution, standard deviation of the normal distribution</td></tr><tr><td>σ<sub>t</sub></td><td>Turbulent Schmidt number (–)</td></tr><tr><td>Φ<sub>ISO</sub></td><td>Total suspended-sand flux through a cross-section calculated using the ISO method (kg s<sup>-1</sup>)</td></tr><tr><td>Φ<sub>SDC</sub></td><td>Total suspended-sand flux through a cross-section calculated using the SDC method (kg s<sup>-1</sup>)</td></tr><tr><td>Φ<sub>total</sub></td><td>Total suspended-sand flux through a cross-section (kg s<sup>-1</sup>)</td></tr></table>
**Code availability.** The open-source toolbox is available at https://gitlab.irstea.fr/jessica.laible/analysis-solid-gauging (last access: 28 June 2023) (https://doi.org/10.5281/zenodo.13973594, Laible et al., 2023a).

**Data availability.** Data are available via https://doi.org/10.57745/NLFT7Q (Laible et al., 2023b) and https://doi.org/10.57745/YTCYSX (Marggraf et al., 2023).

**Supplement.** The supplement related to this article is available online at: https://doi.org/10.5194/esurf-12-1243-2024-supplement.

**Author contributions.** GD, JM, JLC, and BeC conceptualized this work. The methodology was developed by GD, JM, JLC, and BeC. The code and uncertainty part was written by JM and BIC. GD and JM wrote the original draft of this paper, which was edited by DT, JLC, BIC, and BeC. JLC and BeC supervised the work. All authors brought ideas into the work, participated in the field measurements, and reviewed the final draft.

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

**Disclaimer.** Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the US Government.

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

**Acknowledgements.** The authors thank Thibault Vassor and Alexis Pavaux for their help in this project and Janna Stepanian (INRAE) for her help in the sample analysis. Thanks go to Kurt Spicer (USGS) for sharing data. Thanks go to Benjamin Renard for his valuable help and review of the uncertainty part. Thanks go to the Plan Loire and Stephane Rodrigues from Tours University for lending us the US-P06.

**Financial support.** Funding for the work was provided by the Compagnie Nationale du Rhône (CNR), Electricité de France (EDF), European Regional Development Fund (ERDF), Agence de l'eau RMC, and three regional councils (Auvergne-Rhône-Alpes, PACA, and Occitanie) in the context of the Rhône Sediment Observatory (OSR, http://www.graie.org/osr, last access: 20 April 2024).

**Review statement.** This paper was edited by Kieran Dunne and reviewed by two anonymous referees.

## References

Armijos, E., Crave, A., Espinoza, R., Fraizy, P., Santos, A. D., Sampaio, F., De Oliveira, E., Santini, W., Martinez, J., 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.

ASTM D3977: Standard test method for determining sediment concentration in water samples, ASTM International, https://doi.org/10.1520/D3977-97R19, 2007.

Boldt, J. A.: From mobile ADCP to high-resolution SSC: a cross-section calibration tool, in: 3rd Joint Federal Interagency Conference on Sedimentation and Hydrologic Modeling, Reno, Nevada, 19–23 April 2015, 1258–1260, https://water.usgs.gov/osw/SALT/documents/189_Boldt.pdf (last access: 22 October 2024), 2015.

Bouchez, J., Métivier, F., Lupker, M., Maurice, L., Perez, M., Gaillardet, J., and France-Lanord, C.: Prediction of depth-integrated fluxes of suspended sediment in the Amazon River: Particle aggregation as a complicating factor, Hydrol. Process., 25, 778–794, https://doi.org/10.1002/hyp.7868, 2011.

Camenen, B.: Simple and general formula for the settling velocity of particles, J. Hydraul. Eng., 133, 229–233, 2007.

Camenen, B. and Larson, M.: A unified sediment transport formulation for coastal inlet application, Tech. rep., US Army Corps of Engineers, Engineer Research and Development Center, https://doi.org/10.21236/ada472064, 2007.

Camenen, B. and Larson, M.: A general formula for noncohesive suspended sediment transport, J. Coastal Res., 24, 615–627, 2008.

Camenen, B., Holubova, K., Lukac, M., Le Coz, J., and Paquier, A.: Assessment of Methods Used in 1D Models for Computing Bed-Load Transport in a Large River: The Danube River in Slovakia, J. Hydraul. Eng., 137, 1190–1199, https://doi.org/10.1061/(ASCE)HY.1943-7900.0000427, 2011.

Camenen, B., Le Coz, J., Dramais, G., Peteuil, C., Fretaud, T., Falgon, A., Dussouillez, P., and Moore, S. A.: A simple physically-based model for predicting sand transport dynamics in the Lower Mekong River, River Flow, Proc. River Flow conference, 3–5 September 2014, Lausanne, Switzerland, 2189–2197, ISBN 9780429069246, 2014.

Camenen, B., Dramais, G., Laible, J., Le Coz, J., Pierrefeu, G., and Lauters, F.: Quantification of continuous sand flux time-series downstream of a dam during a flushing event, Environ. Fluid Mech., https://doi.org/10.1007/s10652-023-09955-9, 2023.

Colby, B.: Discharge of sands and mean-velocity relationships in sand-bed streams, Tech. rep., U.S. Geological Service, https://doi.org/10.3133/pp462A, 1964.

Delile, H., Masson, M., Miège, C., Le Coz, J., Poulier, G., Le Bescond, C., Radakovitch, O., and Coquery, M.: Hydroclimatic drivers of land-based organic and inorganic particulate micropollutant fluxes: The regime of the largest river water inflow of the Mediterranean Sea, Water Res., 185, 116067, https://doi.org/10.1016/j.watres.2020.116067, 2020.
Despax, A., Le Coz, J., Hauet, A., Mueller, D. S., Engel, F. L., Blanquart, B., Renard, B., and Oberg, K. A.: De-composition of Uncertainty Sources in Acoustic Doppler Current Profiler Streamflow Measurements Using Repeated Measures Experiments, Water Resour. Res., 55, 7520–7540, https://doi.org/10.1029/2019WR025296, 2019.

Despax, A., Le Coz, J., Mueller, D. S., Hauet, A., Calmel, B., Pierrefeu, G., Naudet, G., Blanquart, B., and Pobanz, K.: Validation of an Uncertainty Propagation Method for Moving-Boat Acoustic Doppler Current Profiler Discharge Measurements, Water Resour. Res., 59, e2021WR031878, https://doi.org/10.1029/2021WR031878, 2023.

Dominguez Ruben, L., Szupiany, R., Latosinski, F., C., L. W., Wood, M., and Boldt, J.: Acoustic Sediment Estimation Tool-box (ASET): A software package for calibrating and processing TRDI ADCP data to compute suspended-sediment transport in sandy rivers, Comput. Geosci.-UK., 140, 104499, https://doi.org/10.1016/j.cageo.2020.104499, 2020.

Dramais, G.: Observation et modélisation des flux de sable dans les grands cours d'eau, PhD thesis, University of Lyon, https://theses.hal.science/tel-03188258 (last access: 22 October 2024), 2020.

Edwards, T. K. and Glysson, G. D.: Field methods for measurement of fluvial sediment, U.S. Geological Survey, https://pubs.usgs.gov/twri/twri3-c2/pdf/TWRI_3-C2.pdf (last access: 22 October 2024), 1999.

Filizola, N., Seyler, F., Mourão, M. H., Arruda, W., Spínola, N., and Guyot, J. L.: Study of the variability in suspended sediment discharge at Manacapuru, Amazon River, Brazil, Latin American Journal of Sedimentology and Basin Analysis, 16, 93–99, 2009.

FISP: Laboratory investigation of suspended-sediment samplers, Tech. Rep. Report No. 5, Federal Interagency Sedimentation Project, https://water.usgs.gov/fisp/docs/Report_5.pdf (last access: 22 October 2024), 1941.

FISP: The design of improved types of suspended sediment samplers, Tech. Rep. Report No. 6, Federal Interagency Sedimentation Project, https://water.usgs.gov/fisp/docs/Report_6.pdf (last access: 22 October 2024), 1952.

Gitto, A. B., Venditti, J. G., Kostaschuk, R., and Church, M.: Representative point-integrated suspended sediment sampling in rivers, Water Resour. Res., 53, 2956–2971, https://doi.org/10.1002/2016WR019187, 2017.

Gordon, J. D.: US Geological Survey Quality-assurance Project for Sediment Analysis, Tech. rep., U.S. Geological Service, 2000.

Gordon, J. D., Newland, C. A., and Gagliardi, S. T.: Laboratory Performance in the Sediment Laboratory Quality-assurance Project, 1996–98, Water-Resources Investigations Report, 99, 4184, https://qsb.usgs.gov/slqa/WRIR_99_4184.pdf (last access: 22 October 2024), 2000.

Gray, J. and Gartner, J.: Overview of selected surrogate technologies for high-temporal resolution suspended sediment monitoring, in: Proceedings of the 2nd Joint Federal Interagency Conference, Las Vegas, NV, USA, 27 June–1 July, Citeseer, https://www.usgs.gov/publications/overview-selected-surrogate-technologies-high-temporal-resolution-suspended-sediment (last access: 2 June 2022), 2010.

Gualtieri, C., Angeloudis, A., Bombardelli, F., Jha, S., and Stoesser, T.: On the values for the turbulent

Schmidt number in environmental flows, Fluids, 2, 17, https://doi.org/10.3390/fluids2020017, 2017.

Guy, H. P. and Norman, V. W.: Field methods for measurement of fluvial sediment, United State Geological Survey, Book 3, https://doi.org/10.3133/twri03C2, 1970.

Hoffmann, T., Thorndycraft, V., Brown, A., Coulthard, T., Damnati, B., Kale, V., Middelkoop, H., Notebaert, B., and Walling, D.: Human impact on fluvial regimes and sediment flux during the Holocene: Review and future research agenda, Global Planet. Change, 72, 87–98, https://doi.org/10.1016/j.gloplacha.2010.04.008, 2010.

Hunt, J.: On the turbulent transport of a heterogeneous sediment, Q. J. Mech. Appl. Math., 22, 235–246, 1969.

International Organization for Standardization: Measurement of liquid flow in open channels. Methods for measurement of characteristics of suspended sediment, ISO 4363, Geneva, Switzerland, https://www.iso.org/obp/ui/#iso:std:iso:4363:ed-3:v1:en (last access: 24 October 2024), 2002.

International Organization for Standardization: Hydrometry. Measurement of liquid flow in open channels using current-meters or floats using point velocity measurements, ISO 748, Geneva, Switzerland, 46 pp., https://www.iso.org/obp/ui/#iso:std:iso:748:ed-5:v1:en (last access: 22 October 2024), 2009.

JCGM: Evaluation of measurement data – Guide to the expression of uncertainty in measurement – JCGM, ISO, 50, 134, https://www.bipm.org/documents/20126/2071204/JCGM_100_2008_E.pdf/cb0ef43f-baa5-11cf-3f85-4dcd86f77bd6 (last access: 22 October 2024), 2008.

Kästner, K., Hoitink, A. J. F., Torfs, P. J. J. F., Vermeulen, B., Ningsih, N. S., and Pramulya, M.: Prerequisites for accurate monitoring of river discharge based on fixed-location velocity measurements, Water Resour. Res., 54, 1058–1076, https://doi.org/10.1002/2017WR020990, 2018.

Khodashenas, S. R. and Paquier, A.: A geometrical method for computing the distribution of boundary shear stress across irregular straight open channels, J. Hydraul. Res., 37, 381–388, 1999.

Kondolf, G. M.: Hungry water: effects of dams and gravel mining on river channels, Environ. Manage., 21, 533–551, 1997.

Kondolf, G. M., Gao, Y., Annandale, G. W., Morris, G. L., Jiang, E., Zhang, J., Cao, Y., Carling, P., Fu, K., Guo, Q., Hotchkiss, R., Peteuil, C., Sumi T., Wang, H.-W., Wang, Z., Wei, Z., Wu B., Wu, C., and Yang C. T.: Sustainable sediment management in reservoirs and regulated rivers: Experiences from five continents, Earths Future, 2, 256–280, 2014.

Laible, J., Calmel, B., Le Coz, J., Camenen, B., Dramais, G., and Vassor, T.: Analysis solid gauging, Zenodo [code], https://doi.org/10.5281/zenodo.13973594, 2023a.

Laible, J., Dramais, G., Camenen, B., Le Coz, J., Topping, D. J., Santini, W., and Pierrefeu, G.: Data set of solid gaugings in several rivers, Recherche Data Gouv [data set], https://doi.org/10.57745/NLFT7Q, 2023b.

Lennermark, M. and Hauet, A.: Developing a post-processing software for ADCP discharge measurement piloted by an international and inter-agency group: a unique, ambitious experience... and one that works!, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-9379, https://doi.org/10.5194/egusphere-egu22-9379, 2022.

Lupker, M., France-Lanord, C., Lavé, J., Bouchez, J., Galy, V., Mélivier, F., Gaillardet, J., Lartiges, B., and Mugnier, J.: A Rouse-
based method to integrate the chemical composition of river sediments: Application to the Ganga basin, J. Geophys. Res.-Earth, 116, 1-24, https://doi.org/10.1029/2010JF001947, 2011.

Mansanarez, V., Renard, B., Le Coz, J., Lang, M., and Darienzo, M.: Shift Happens! Adjusting Stage-Discharge Rating Curves to Morphological Changes at Known Times, Water Resour. Res., 55, 2876-2899, https://doi.org/10.1029/2018WR023389, 2019.

Marggraf, J., Camenen, B., Le Coz, J., Dramais, G., Lauters, F., Pierrefeu, G.: Suspended sediment measurements in the Isère River at Grenoble Campus, Recherche Data Gouv [data set], https://doi.org/10.57745/YTCYSX, 2023.

McLean, S.: On the calculation of suspended load for noncohesive sediments, J. Geophys. Res.-Oceans, 97, 5759-5770, 1992.

Mueller, D. S.: QRev-Software for computation and quality assurance of acoustic Doppler current profiler moving-boat streamflow measurements, U.S. Geological Survey Open-File Report, 1052, 50, Technical Manual for version 2.8, https://doi.org/10.3133/ofr20161068, 2016.

Némery, J., Mano, V., Coynel, A., Etcheber, H., Moatar, F., Meybeck, M., Belleudy, P., and Poirel, A.: Carbon and suspended sediment transport in an impounded alpine river (Isère, France), Hydrol. Process., 27, 2498-2508, https://doi.org/10.1002/hyp.9387, 2013.

Oberg, K. and Mueller, D. S.: Validation of streamflow measurements made with acoustic Doppler current profilers, J. Hydraul. Eng., 133, 1421-1432, https://doi.org/10.1061/(ASCE)0733-9429(2007)133:12(1421, 2007.

Parsons, D. R., Jackson, P., Czuba, J., Engel, F., Rhoads, B. L., Oberg, K., Best, J. L., Mueller, D., Johnson, K., and Riley, J.: Velocity Mapping Toolbox (VMT): a processing and visualization suite for moving-vessel ADCP measurements, Earth Surf. Processes, 38, 1244-1260, https://doi.org/10.1002/esp.3367, 2013.

Perret, E., Camenen, B., Berni, C., El kadi Abderrezzak, K., and Renard, B.: Uncertainties in Models Predicting Critical Bed Shear Stress of Cohesionless Particles, J. Hydraul. Eng., 149, 04023002, https://doi.org/10.1061/JHEND8.HYENG-13101, 2023.

Porterfield: Computation of fluvial-sediment discharge, Techniques of Water-Resources Investigations, https://doi.org/10.3133/twri03C3, 1972.

Poulier, G., Launay, M., Le Bescond, C., Thollet, F., Coquery, M., and Le Coz, J.: Combining flux monitoring and data reconstruction to establish annual budgets of suspended particulate matter, mercury and PCB in the Rhône River from Lake Geneva to the Mediterranean Sea, Sci. Total Environ., 658, 457-473, https://doi.org/10.1016/j.scitotenv.2018.12.075, 2019.

Renard, B., Garreta, V., and Lang, M.: An application of Bayesian analysis and Markov chain Monte Carlo methods to the estimation of a regional trend in annual maxima, Water Resour. Res. 42, W12422, https://doi.org/10.1029/2005WR004591, 2006.

Rouse, H.: Modern conceptions of the mechanics of fluid turbulence, Trans. ASCE, 102, 463-505, 1937.

Rozovskii, I. L.: Flow of water in bends of open channels, Academy of Sciences of the Ukrainian SSR: Kiev, translated from Russian by the Israel Program for Scientific Translations, Jerusalem, 1961, 1957.

Sabol, T. A. and Topping, D. J.: Evaluation of intake efficiencies and associated sediment-concentration errors in US D-77 bag-type and US D-96-type depth-integrating

suspended-sediment samplers, U.S. Geological Survey Scientific Investigations Report 2012-5208, 88 pp., 5208, 88, https://doi.org/10.3133/sir20125208, 2013.

Santini, W., Camenen, B., Le Coz, J., Vauchel, P., Guyot, J.-L., Lavado, W., Carranza, J., Paredes, M. A., Pérez Arévalo, J. J., Arévalo, N., Espinoza Villar, R., Julien, F., and Martinez, J.-M.: An index concentration method for suspended load monitoring in large rivers of the Amazonian forest, Earth Surf. Dynam., 7, 515-536, https://doi.org/10.5194/esurf-7-515-2019, 2019.

Shah-Fairbank, S. C. and Julien, P. Y.: Sediment load calculations from point measurements in sand-bed rivers, Int. J. Sediment Res., 30, 1-12, https://doi.org/10.1016/S1001-6279(15)60001-4, 2015.

Smart, G.: A base for the log law and von Karman's constant problem, J. Hydraul. Res., 60, 935-943, https://doi.org/10.1080/00221686.2022.2076164, 2022.

Soulsby, R. and Whitehouse, R.: Threshold of sediment motion in coastal environments, in: Pacific Coasts and Ports' 97: Proceedings of the 13th Australasian Coastal and Ocean Engineering Conference and the 6th Australasian Port and Harbour Conference, Christchurch, New Zealand, 7-11 September 1997, Vol. 1, Centre for Advanced Engineering, University of Canterbury, 145-150, https://search.informit.org/doi/epdf/10.3316/informit.929741720399033 (last access: 22 October 2024), 1997.

Spicer, K.: P-6 Comparison Tests. A Point-Integrating Suspended Sediment Sampler Comparison, SED-HYD Conference, Reno, NV, USA, 24-28 June 2019, http://carnettechnology.com/assets/P-6_Comparison_Tests-A-Point_Integrating_Susped_Sediment_Compa59646.pdf (last access: 22 October 2024), 2019.

Starosolsky, O. and Rakoczi, L.: Operational hydrology report (OHR), 16. Measurement of river sediments: prepared by the Rapporteur on Sediment Transport of the Commission for Hydrology, World Meteorological Organisation, https://library.wmo.int/viewer/33602/#page=1&viewer=picture&o=bookmarks&n=0&q= (last access: 22 October 2024), 1981.

Szupiany, R. N., Lopez Weibel, C., Guerrero, M., Latosinski, F., Wood, M., Dominguez Ruben, L., and Oberg, K.: Estimating sand concentrations using ADCP-based acoustic inversion in a large fluvial system characterized by bi-modal suspended-sediment distributions, Earth Surf. Processes, 44, 1295-1308, https://doi.org/10.1002/esp.4572, 2019.

Topping, D. J. and Wright, S. A.: Long-term continuous acoustical suspended-sediment measurements in rivers – Theory, application, bias, and error, U.S. Geological Survey Professional Paper 1823, 98 pp., https://doi.org/10.3133/pp1823, 2016.

Topping, D. J., Rubin, D. M., Wright, S. A., and Melis, T. S.: Field evaluation of the error arising from inadequate time averaging in the standard use of depth-integrating suspended-sediment samplers, US Geological Survey Professional Paper 1774, 95 pp., https://pubs.usgs.gov/pp/1774/pp1774.pdf (last access: 18 February 2024), 2011.

Topping, D. J., Grams, P. E., Griffiths, R. E., Dean, D. J., Wright, A. S., and Unema, J. A.: Self-limitation of sand storage in a bedrock-canyon river arising from the interaction of flow and grain size, J. Geophys. Res.-Earth, 126, e2020JF005565, https://doi.org/10.1029/2020JF005565, 2021.
U.S. Geological Survey: USGS water data for the Nation.; U.S. Geological Survey National Water Information System database [data set], https://doi.org/10.5066/F7P55KJN, 2023.

Van Rijn, L. C.: Sediment transport, part II: suspended load transport, J. Hydraul. Eng., 110, 1613–1641, 1984.

Van Rossum, G. and Drake, F. L.: Python 3 reference manual, CreateSpace, Scotts Valley, CA, 242 pp., ISBN 978-1-4414-1269-0, 2009.

Vauchel, P., Santini, W., Guyot, J. L., Moquet, J. S., Martinez, J. M., Espinoza, J. C., Baby, P., Fuertes, O., Noriega, L., Puita, O., et al.: A reassessment of the suspended sediment load in the Madeira River basin from the Andes of Peru and Bolivia to the Amazon River in Brazil, based on 10 years of data from the HYBAM monitoring programme, J. Hydrol., 553, 35–48, https://doi.org/10.1016/j.jhydrol.2017.07.018, 2017.

Venditti, J., Church, M., Attard, M., and Haught, D.: Use of ADCPs for suspended sediment transport monitoring: An empirical approach, Water Resour. Res., 52, 2715–2736, https://doi.org/10.1002/2015WR017348, 2016.

Vergne, A., Le Coz, J., Berni, C., and Pierrefeu, G.: Using a Down-Looking Multifrequency ABS for Measuring Suspended Sediments in Rivers, Water Resour. Res., 56, e2019WR024877, https://doi.org/10.1029/2019WR024877, 2020.

Vergne, A., Le Coz, J., and Berni, C.: Some Backscatter Modeling Issues Complicating the Sonar-Based Monitoring of Suspended Sediments in Rivers, Water Resour. Res., 59, e2022WR032341, https://doi.org/10.1029/2022WR032341, 2023.

Vermeulen, B., Sassi, M. G., and Hoitink, A. J. F.: Improved flow velocity estimates from moving-boat ADCP measurements, Water Resour. Res., 50, 4186–4196, 2014.

Wren, D., Barkdoll, B., Kuhnle, R., and Derrow, R.: Field techniques for suspended-sediment measurement, J. Hydraul. Eng., 126, 97–104, https://doi.org/10.1061/(ASCE)0733-9429(2000)126:2(97), 2000.