Applicability of Acoustic Concentration Measurements in Suspensions of Artificial and Natural Sediments Using an Acoustic Doppler Velocimeter

For the investigation of turbulence and particles, interaction measurement systems are required, which are able to measure velocity and concentration fluctuations simultaneously. Acoustic Doppler Velocimeters (ADV) are widely used for velocity and turbulence measurements in natural and artificial flows. Based on the acoustic sonar theory, a model is presented that correlates the ADV’s Signal-to-Noise Ratio (SNR) and the suspended solids concentration of several natural (Ems Estuary, Lake Eixendorf, Lake Altmühl) and artificial sediments (Chinafill, quartz powder, bentonite, metakaolin) for the range 0.001 g/L–50 g/L. Within the presented method, the sound absorption in water and on particles is considered in a continuous approach for sampling frequencies up to 100 Hz. The widely-used log-linear relation between the SNR and the concentration, which is only valid for low concentrations, was extended for the high concentration regime. Measurement results show a similar behavior of the SNR with respect to varying suspended solid concentrations for different sediments. However, the analysis of the fit parameters shows systematic differences depending on the type of sediment. It is concluded that the proposed model is applicable as well for laboratory use as for measurements in rivers and estuaries. Finally, we discuss the reliability of the results and the methodology with regard to measurements in rivers, lakes, and estuaries.


Introduction
ADV probes provide the opportunity to measure in situ the three-dimensional flow velocity in liquids with a high temporal resolution (up to 200 Hz).This allows evaluating the velocity fluctuations in terms of the turbulent behavior in the flow.The three-dimensional velocity fluctuations u i = u i − ūi are described as the deviation of the mean value ūi .The intensity of turbulence is described by the turbulent kinetic energy k = 1  2 u i u i , which consists of the diagonal contributions of the Reynolds stress tensor, e.g., [1].
Especially in natural surface waters like rivers, channels, lakes, and estuaries, near-bottom suspended solids of different concentrations are common.For instance, in the Ems Estuary in the German Bight, Suspended Solid Concentrations (SSC) were measured up to 300 g/L [2].Cellino demonstrated that SSC influence the flow behavior within a stream, in particular the vertical structure of turbulence.It was shown that SSC can lead to turbulence damping and decreasing eddy viscosity [3].
The common method to describe the interaction between concentration and turbulence is to quantify the turbulent Schmidt-number Sc t = ν t K t with ν t being the turbulent viscosity and K t being the turbulent diffusivity.The physical behavior of both quantities is the same, as the diffusivity of particles increases with turbulence.The turbulent Schmidt-number is often used for calibrating experimental and numerical data.Depending on the type of particles, it is given as Sc t = 0.3-1.5 [4].
Applying the eddy viscosity hypothesis, ν t can be formulated by means of the extra diagonals of the Reynolds stress tensor u w and the vertical gradient of the mean velocity ∂ ū ∂z in the main stream direction: Following the same principle, K t is described by the fluctuations of the concentration c and the vertical gradient of the mean concentration.
In order to investigate the interaction of turbulence and particle concentration, it is necessary to measure the velocity and the concentration fluctuations simultaneously.Further, to detect as many turbulent eddies as possible, the measuring frequency should be as high as possible.
Salehi et al. summarized three methods for measuring the particle concentrations in currents [5]: (1) the extraction of samples at different positions with a suction device and the external analysis of the samples taken; (2) determination of the particle concentration through scattering and reflection of light (OBS); and (3) the concentration measurement by means of the reflection and scattering of an acoustic signal.The intensity of sound is attenuated by the particles within the measurement volume.Furthermore, the reflected signal is superimposed by noise, which is taken into account by the Signal-to-Noise Ratio (SNR).
As mentioned in Equation ( 2), for the investigation of the diffusivity coefficient K t , measurements of the concentration fluctuations are necessary.Applying the suction method, for instance, no information about the concentration and velocity fluctuations at the same time and location are obtained.Optical measurement devices provide a wide spectrum of recordable concentrations; however, they are not able to measure the velocity.The investigations of Fugate et al. show higher accuracy for SSC measurements in estuaries when using acoustic devices rather than optical instruments [6].Acoustic devices, like ADV probes, are used to measure the three-dimensional flow velocities with a sampling frequency f up to 200 Hz.Additionally, the SNR is recorded, which is affected by the particle concentration in the measurement volume.If it is possible to find a unique function for the SNR with respect to the concentration, then it is possible to measure velocity, turbulence, and concentration synchronously.
The foundation of the acoustic method is the sonar theory, which was first applied for small particles by Urick in 1948 [7].Urick investigated the absorption behavior of kaolin and fine silt while using a quartz crystal to send acoustic pulses.After ADV measurements became a standard technique for velocity measurements, the acoustic scattering theory again was applied for concentration measurements [8].Thorne et al. gave a summary of the advantages of acoustic methods for concentration measurements since optical devices reached their working range at this time [9].Jay et al. described difficulties when calibrating the SNR for natural rivers, since the reflected signal was strongly dependent on the particle's material, shape, and size [10].Hoitink et al. were able to calibrate the SNR for concentrations up to 1 g/L [11].However, neither study, from Jay et al. and Hoitink et al., considered the sound absorption on particles; that is why their methods are only valid up to approximately 1 g/L.
The acoustic method for concentration measurement can be applied at laboratory conditions with a homogeneous particle distribution and known sediment properties.Ha et al. and Salehi et al. presented laboratory experiments for calibrating the reflected signal intensity and the particle concentration [5,12].Nevertheless, they neglected the sound attenuation on particles, although it is a crucial parameter.Decrop et al. presented an experimental setup to measure concentrations up to 10 g/L with an ADV probe while considering sound attenuation on particles.Their measurements were conducted with a measurement frequency of 25 Hz [13].To investigate the interaction of turbulent concentration and velocity fluctuations, higher sampling frequencies are required to resolve turbulent flow behavior.
In order to investigate the interaction of turbulence and concentration, different experimental studies were realized.Best et al. conducted a laboratory experiment to measure velocities and concentrations with a Laser Doppler Anemometer (LDA).However, the application of their optical method was restricted to low concentrations [14].Cellino used a self-constructed acoustic particle flux profiler for measuring the concentration and velocity with a frequency of 16 Hz [15].A part of his research was to investigate the decay of turbulence with increasing suspended sediment concentration in a laboratory flume [3].
Moreover, Hurther et al. showed the development of an Acoustic Concentration and Velocity Profiler (ACVP), which was able to measure the two-dimensional velocity field and the concentration over a vertical profile of 25 cm [16].This device was applied by Revil-Baudard et al. for laboratory experiments on investigations of sheet-flow processes [17].They reported about maximum concentration measurement frequencies of 4.9 Hz.
The necessity of simultaneously measuring flow velocity and concentration is presented by means of the Ems Estuary (Germany, North Sea) as an example.In the Ems Estuary, simultaneous measurements of flow velocity and concentration are needed in order to describe the dynamics of fluid mud layers.Fluid mud is a high concentrated mud suspension, which has huge impacts on the maintenance of rivers and estuaries, since navigation can be impeded [18].Recently, Becker et al. showed concentration and velocity measurements of the fluid mud layers in the Ems Estuary and highlighted the unsteady and complex behavior of the mud layers with tide.Within this measurement campaign, velocity and concentration were measured by different devices.Furthermore, no turbulence measurements were presented [19].Having a measurement device that would be able to measure velocity, concentration, and turbulence simultaneously, the unsteady dynamics of such fluid mud layers could be investigated in more detail.
With developing measuring techniques, more accurate turbulence results are expected when using measurement devices with sampling frequencies as high as possible.In this case, small eddies are able to be resolved, which could not be recognized from the measurement device in former studies.Therefore, a Vectrino Profiler (Nortek, Norway) was used with a sampling frequency of 100 Hz.The Nortek Vectrino Profiler was launched in 2010 and consists of a single central transceiver and four passive receivers at an angle of 30 • towards the transceiver.A detailed description of the probe was given by Craig et al. and by Thomas et al. [20,21].
Within this paper, a method is presented to model the concentration of different suspended materials by means of the SNR with a modern ADV probe (Vectrino Profiler).The acoustic model is based on the well-known sonar theory approach, and it is applied for artificial and natural river sediments.The authors show how sound attenuation is considered in water and on sediment particles within the sonar theory-based model.In experimental studies for cohesive sediment dynamics, quartz powder, metakaolin, and Chinafill are often used [13,22,23].Therefore, the following artificial sediments were investigated: quartz powder, bentonite, metakaolin, and Chinafill.Afterwards, the results were compared with those of natural sediments from the Ems Estuary on the German Coast, Lake Altmühl, and Lake Eixendorf, both located in Bavaria, Germany.The experimental setup and the proceeding measurements are presented.Furthermore, the authors describe the analysis of the results of different sediments in order to specify the model's applicability in nature.Furthermore, it is shown that concentrations can be measured up to 50 g/L depending on the type of sediment.
This paper describes a method for simultaneously measuring velocity, turbulence, and concentration, which can be applied in laboratory investigations, as well as in rivers, estuaries, or lakes.The characteristics of the method are presented in detail, and indications are given for practical use in nature.The results of this study provide a contribution for further investigations of turbulence and diffusivity equally in the laboratory and for use in nature.

Acoustics of the ADV Probe
In the following, the theoretical aspects are described for sound and particles' interaction in water.First, the particle diameters are verified, which are applicable to the acoustic backscatter theory.Afterwards, the authors present the implementation of the general sonar theory using the backscatter intensity.The aim of this part is to derive an equation to relate the suspension concentration to an acoustic quantity, which can be measured by the ADV probe.

Rayleigh Scattering and the Intensity of the Backscatter Signal
The working principle of ADV devices is based on Rayleigh scattering, which occurs when the diameter of the particle is smaller than the wavelength of the emitted signal.In that case, the energy transfer towards the particle can be neglected.For Rayleigh scattering, the condition: has to be fulfilled, with k w = 2π λ as the wave number, λ as the wavelength, and r as the particle radius [24].Decrop et al. argued that a sufficiently high backscatter intensity only can be achieved for k w r > 0.05 [13].
In this study, we used different sediments with mean particle diameters 7 µm ≤ d 50 ≤ 20 µm.The ADV's emitter frequency was F = 10 MHz, which lead to 0.147 ≤ k w r ≤ 0.419 (Table 1).Therefore, the condition was fulfilled to apply the principles of Rayleigh scattering here.The intensity S v of the backscatter signal increases with the reflecting particle amount n b inside the measurement volume and their scattering cross-section σ bs .Assuming unidirectional reflection, Medwin and Clay presented the scattering cross-section as: The factor Υ is a material parameter, which is determined as: where ∆ρ = ρ s /ρ w is the density ratio of the particle and water and ∆E = E s /E w is the elasticity ratio.
The backscatter signal can be described in terms of n b and σ bs as the volume backscattering coefficient: which is a measure of the backscatter strength from particles.The backscatter strength increases with the number of reflecting particles and with the size of the scattering cross-section.That means that an increase of s v is equivalent to a greater backscatter intensity, which is received by the ADV probe.
The backscatter strength can be expressed in decibel units as [11,24]: Assuming spherical particles with a mass M = 4/3r 3 πρ s , the concentration: can be calculated from the number of particles inside the scattering volume.Combining Equations ( 4), ( 6) and ( 8) leads to: which describes the backscatter intensity as a function of the concentration.The backscatter intensity S v must not be mistaken with the ADV receiver intensity.S v is rather the backscatter intensity in the vicinity of the particle.On the way back to the ADV receiver, the signal intensity is exposed to further losses.

Application of Sonar Theory
The ADV transmitter emits an acoustic signal, which will be reflected on particles in the water, but it is only partly received by the probe's receiver.This process can be described by the principles of active sonar, which are illustrated in Figure 1.All quantities within the sonar theory are described in dB.The emitted signal intensity (Signal Level (SL)) is exposed to Transmission Losses (TL) in both traveling directions.Additionally, the reflected signal is influenced by the specifics of the particles (e.g., shape, size, elasticity), which are accounted for in the Target Strength (TS).The ambient sound, which is not related to the ADV probe, is called noise and is described in NL as the noise level.Within sonar theory, SE (Signal Excess) is the received backscatter intensity inclusive of noise.In order to detect the reflected signal, SE has to overcome at least the noise intensity (Noise Level (NL)).The sonar balance equation as presented in similar form in Hoitink et al. and in Salehi et al. reads [5,11]: (SE − NL) is the backscatter intensity under exclusion of noise (it is also called the Signal-to-Noise Ratio (SNR)).The target strength within the measurement volume is the backscatter intensity S v .C ADV describes the influence of additional device-specific settings.The transmission losses 2TL = 40 log 10 (R) + 2αR are accounted for in both directions, while R is the distance from the transmitter to the measurement volume [25].The effect of sound absorption is considered by: as the sum of the occurring sound absorption losses in water α w and on particles α p .The ADV device, which is used in the experiment, has an emitter frequency of F = 10 MHz.For this emitter frequency, the coefficient for absorption losses in water is given as α w = 26.9dB/m [26].Unfortunately, the particle absorption coefficient α p is not a constant.It depends on the acoustic properties of the transmitted signal and the properties of the particles.Urick presented an approach for the particle absorption coefficient in granular suspensions, which is valid for homogeneous particle distributions [7].This approach was successfully applied by Hoitink et al. and Salehi et al. when conducting ADCP (Acoustic Doppler Current Profiler) and ADV measurements [5,11].Therefore, Urick's formula for α p was used within this paper, while the constant term is called the Urick-coefficient Ur [m 2 /kg].Ur is calculated by sediment and material constants, namely the particle radius r, the solid density ρ s , the specific density ratio ∆ρ = ρ s /ρ w , the molecular viscosity ν m , the emitter frequency F, and its wave number k w : The following auxiliary variables are introduced: It is highlighted that next to sedimentological and device-specific parameters, Equation ( 12) is a function of the concentration C.

How Particles Influence the SNR
The intensity of the backscatter signal, which is received by the ADV, depends on the sediment concentration.Two effects have to be considered simultaneously.First, an increasing particle concentration will be accompanied by an increase of the overall scattering cross-section, if a homogeneous particle distribution is ensured over the volume.This leads to rising SNR values.Second, while α p is a function of the concentration, increasing sound absorption causes a decrease of the SNR.
Combining Equations ( 9) and ( 10), splitting 2αR =2α p R+2α w R, introducing the transformation parameter Π 1 , and solving for the concentration leads to: with: In order to quantify the dependency of the SNR and the concentration, we present an experiment to measure the SNR in suspensions of pre-defined concentrations.Afterwards, the measurements will be compared to the results gained by the sonar theory model.

Experimental Investigation
After deriving the relation between the concentration and the SNR, within this section, the authors present the experimental setup for calibrating the SNR to the concentration of suspension flows.After discussing the suitability of the used ADV device, the exemplary results are presented.On behalf of simplicity, the description of the experimental setup is restricted only to one sediment type, namely quartz powder.Results of the resuming sediments are shown in Section 4.  Before adding the sediment, an initial reference SNR was measured to account for the prevailing contamination in the water used.The sediment was weighed with a high precision electronic scale in order to prepare the desired concentrations.The concentration was stepwise increased from C = 0.001 g/L-50 g/L.To prevent the aggregation of particles, we initiated a strong turbulence by the stirrers and added the powder through a sieve.To control vertical homogeneous mixing conditions, measurements were conducted at four depths, equally separated.

Suitability of the ADV Probe
To determine the most unambiguous relation between the SNR and the concentration, firstly, a clear SNR signal was needed, and secondly, the deviations within the four receiver signals had to be as small as possible.Figure 3 shows the box plot evaluation of the SNR signal for three repeated measurements in a low concentrated quartz powder suspension.
Version December 12, 2018 submitted to Acoustics 8 of 18

Suitability of the ADV probe
To determine the most unambiguous relation between the SNR and the concentration, firstly a clear SNR signal was needed and secondly the deviations within the four receiver signals had to be as small as possible.Fig. 3 shows the box plot evaluation of the SNR signal for three repeated measurements in a low concentrated quartz powder suspension.The scattering of the interquartile range (IQR) was similar for all four beams and for three repetitions, within the range of 2 -3 dB.It is noticeable that beam 4 was measuring lower values than the remaining receivers.On the contrary, beam 1 was providing slightly greater values.
In the case of Run 03, this behaviour is shown in Fig. 4 where the SNR signals are presented as histograms.While fitting Gaussian curves to the histograms, the mean µ and the standard deviation Due to the narrow range of the IQR and the few detected outliers in the box plots, as well as similar mean values and low standard deviations, it was assumed that this ADV probe is satisfactory suitable for calibration measurements.

Acoustic feedback: Based on quartz powder measurements
After averaging the SNR data of the four receivers, the SNR was plotted against the quartz powder concentration in semi-logarithmic scale.While the SNR of the four receivers is already averaged, Fig. 5 illustrates the dependency of SNR on the concentration for two different measurement heights above the bottom (z 1 = 19 cm, z 2 = 15 cm).Although the measurements at different heights fluctuated in a The scattering of the Interquartile Range (IQR) was similar for all four beams and for three repetitions, within the range of 2-3 dB.It is noticeable that Beam 4 was measuring lower values than the remaining receivers.On the contrary, Beam 1 was providing slightly greater values.
In the case of Run 03, this behavior is shown in Figure 4, where the SNR signals are presented as histograms.While fitting Gaussian curves to the histograms, the mean µ and the standard deviation σ were calculated.It shows that Beams 2 and 3 provided similar mean values µ 2 = 44.96dB and µ 3 = 44.61dB, which is in between µ 1 = 46.03dB and µ 4 = 43.08 dB.Nevertheless, when averaging the data of all four receivers, the equal deviations of Beam 1 and Beam 4 were eliminated.The standard deviations of the SNR distributions were in the range of σ = 2.5-2.7 dB, which corresponded to the values obtained by the box plot evaluation.
Due to the narrow range of the IQR and the few detected outliers in the box plots, as well as similar mean values and low standard deviations, it was assumed that this ADV probe was satisfactorily suitable for calibration measurements.Histograms are given for each receiver (beam) with corresponding mean and standard deviation values.

Acoustic Feedback: Based on Quartz Powder Measurements
After averaging the SNR data of the four receivers, the SNR was plotted against the quartz powder concentration in semi-logarithmic scale.While the SNR of the four receivers was already averaged, Figure 5 illustrates the dependency of SNR on the concentration for two different measurement heights above the bottom (z 1 = 19 cm, z 2 = 15 cm).Although the measurements at different heights fluctuated in a range of about 5 dB, no systematic deviations related to the z-coordinate were noticeable.Therefore, we concluded that the prepared suspensions were well mixed, and a homogeneous concentration was apparent in the experimental box.
In semi-logarithmic scale, a characteristic linear shape can be observed up to C ≈ 1 g/L.For higher concentrations, no linear dependency can be regarded any more.On the basis of the quartz powder measurements, we propose the following subdivision of the chart in Figure 5 into three acoustically-characteristic regions.(I) The whispering zone: The whispering region is characterized by few echoes returning from the measurement volume due to low particle concentration.Within this region, SNR data reacted sensitively to small deviations of the particle amount.The log-linear behavior between the SNR and the concentration was limited in this region since the ADV probe always detected a minimum SNR value.(II) The audibility zone: Within this region, a clear echo intensity was received.The absorption losses on sediment particles were considered to be low.It is common to neglect α p for concentrations up to C ≈ 1 g/L [5].(III) The muffling zone: Due to a high concentration, the acoustic backscatter signal was intense.In this region, the log-linear shape provided no suitable data approximation.Still, the overall scattering coefficient was increasing with concentration, which led to an increase of SNR.At the same time absorption, losses on particles were becoming more dominant.This resulted in a vertex after which the SNR decreased again.When absorption losses overcome the increase of the scattering coefficient, SNR data will decrease as a matter of fact.The continuous flattening of the data points (with increasing concentration) in Figure 5 was due to the increasing influence of α p .Not knowing at which concentration α p would be negligible or not, it is proposed to consider both the absorption coefficients α w and α p irrelevant with respect to the concentration.Since α p is a function of C itself (Equation ( 12)), it will be automatically reduced in the case of low C. Therefore, a steady function SNR = f(C) has to be found.

Modeling the Signal-to-Noise Ratio
In this section, we describe the methodology of fitting the sonar equation parameters (Equation ( 13)) into the measurements.Doing so, the influence of α p is considered.Furthermore, we show how the sonar model is applied to compute the concentration by means of the SNR.In the end, we evaluate and compare the results of different kinds of sediment.

Application of the Sonar Equation
When rearranging Equation ( 13), the SNR is described as: with the three unknown parameters: Π 1 , α p = Π 2 •Ur analyt.•C and C ADV .Π 2 is interpreted as an empirical correction factor of Ur analyt. .Figure 6 describes the experimental data in comparison to the parametrized sonar function, which is based on Equation (14).The parameterization of the sonar equation was done with a non-linear least-squares method in the MATLAB R Curve-Fitting-Tool.It can be seen that the fit function and the measurements are in good agreement.The unknown parameters were found as Π 1 = 1.712,Π 2 = 4.34, and C ADV = 42.08 dB.We mention that the analytical solution of the Urick-coefficient gave Ur analyt.= 1.727 m 2 /kg, which deviated in comparison to Ur sonar = Π 2 • Ur analyt.= 7.43 m 2 /kg.Furthermore, we mention that Π 1 is a transformation coefficient of the SNR (see Equation ( 13)), which is expected to be close to 1.
When considering the region of high concentrations, the sonar model describes a flattening followed by a decrease of SNR.This behavior is explained by the absorption coefficient α p in the sonar equation.Even if the influence of α p remains small in low concentrations, the more particles are in the flow, the higher is the scattering cross-section σ bs , which increases the influence of α p .At a certain concentration, the influence of α p becomes dominant and leads to a signal damping, resulting eventually in a decrease of the SNR function.Especially, SNR measurements at higher concentrations depict this effect.A sudden decrease of SNR data was regarded for all tested sediments.The "turning point" of SNR differed with the type of sediment, thus it was expected to be material dependent.14)).This model predicts a decreasing SNR in the high concentration regime and a log-linear development in the low concentration regime.

Evaluation of the Sonar Model for Different Sediments
Next, with respect to quartz powder, we conducted measurements of three other artificial sediments, namely bentonite, metakaolin, and Chinafill.All artificial sediments were used as dry powders so that no drying or sediment preparation was necessary.In order to testify to our methodology for natural sediments, we repeated the measurements for coastal sediments from the Ems Estuary (North Sea, Germany) and for freshwater sediments from Lake Altmühl and Lake Eixendorf (Bavaria, Germany).Before using the natural sediments, they needed to be dried, crushed, and sieved (<25 mm) in order to gain the necessary powder.
Dry densities of the sediments were measured with a gas pycnometer (Pycnomatic ATC, Porotec), and mean particle diameters were calculated from grain size distributions measured with a particle sizer (partica LA-950, Horiba).A list of the resulting material constants Ur analyt., Υ, and A is given in Table 2.After repeating the calibration experiment for different sediments, the sonar model was fitted to the measurements.The results for the artificial sediments (quartz powder, bentonite, metakaolin, and Chinafill) are plotted in Figure 7.It can be seen that the sonar model was capable of reproducing the experimental datasets.
In the case of artificial sediments, it is recognizable that there was some deviation between the model and the experimental data in the log-linear regime (II: audibility zone).However, with increasing concentration, the deviations become smaller.When the concentration was sufficiently high, a decrease of the SNR was observed.The conditions to pass from increasing to decreasing SNR response (III: muffling zone) differed with the type of sediment.
Version December 12, 2018 submitted to Acoustics 12 of 18 After repeating the calibration experiment for different sediments, the sonar model was fitted to the measurements.The results for the artificial sediments (quartz powder, bentonite, metakaolin and chinafill) are plotted in Fig. 7.It can be seen that the sonar model is capable to reproduce the experimental data sets.
In case of artificial sediments it is recognizable that there is some deviation between the model and the experimental data in the log-linear regime (II: audibility zone).However with increasing concentration the deviations become smaller.When the concentration is sufficiently high, a decrease of the SNR is observed.The conditions to pass from increasing to decreasing SNR response (III: muffling zone) differs with the type of sediment.In comparison to the artificial sediments, the natural sediments from Lake Eixendorf, Lake Altmühl and the Ems Estuary are presented in Fig. 8.The resulting fit parameters are summarized in 2 Chinafill.The fit function is given for the evaluation of the sonar equation (Equation ( 14)).This model predicts a decreasing SNR in the high concentration regime and a log-linear development in the low concentration regime.
In comparison to the artificial sediments, the natural sediments from Lake Eixendorf, Lake Altmühl, and the Ems Estuary are presented in Figure 8.The resulting fit parameters are summarized in Table 3 for all measurements.Additionally, R 2 is given as the square of the correlation between the fitted and the measured values.
Even for natural sediments where the composition of sediment is far more complex compared to artificial sediments, our sonar model was able to reproduce the experimental data.Furthermore, it seems that natural sediments rather follow a log-linear approach in the low concentration regime (II: audibility zone) than artificial sediments.For natural sediments, the slope of the modeling results in the log-linear regime does not differ as much as for the artificial sediments.In the opinion of the authors, this effect is not surprising since natural sediments likely consist of similar sediment partitions.
For natural sediments, the vertices of the modeled SNR were located very close to each other, whereas the variation was greater for artificial sediments.Out of this, we conclude that there exists a sediment-and device-specific critical concentration, which separates the increasing and decreasing SNR response. .This model predicts a log-linear evolution in the low concentration regime and a decreasing SNR in the high concentration regime.
Table 3. Fit parameter and of the sonar model for different sediments.R 2 as a statistical parameter is the square of the correlation between the fitted and the measured values.14)).This model predicts a log-linear evolution in the low concentration regime and a decreasing SNR in the high concentration regime.

Discussion
We analyzed the fit parameters Π 2 and C ADV and separated the data into artificial sediments (quartz powder, bentonite, metakaolin, and Chinafill) and natural sediments (Ems Estuary, Lake Eixendorf, and Lake Altmühl).Knowing that the mentioned artificial sediments are natural materials, but artificially processed, we use this terminus to distinguish between the sediments extracted from nature.When reading the fit parameters in Table 3, differences between Π 2 and C ADV can be observed.Furthermore, it can be seen that Π 2 and C ADV can be separated into artificial and natural sediments.This is shown in Figure 9 for Π 2 and in Figure 10 for C ADV .Π 2 is described through α p, f it /α p,analyt., which is the coefficient of the fitted and the analytically-calculated particle absorption coefficient.A value close to one would imply perfect agreement between the modeled sound absorption on particles and the value based on Urick's formula in Equation (12).
Artificially-processed sediments are characterized by a very homogeneous grain size distribution, which is not or to a negligible extent contaminated by other sediments or particle diameters.The grading curve therefore is steeper than for sediments extracted from nature.The uniformity coefficient C U = d 60 /d 10 is a dimensionless number to describe the uniformity of a sediment sample.The diameters d 60 and d 10 correspond to 60% and 10% sieve passage, respectively.The sediment samples' uniformity coefficients C U are shown in Table 2.The grading curves of the artificial sediments result in lower C U -values than for natural sediments.The formula of Urick (Equation ( 12)), which is implemented in the model, is based on the assumption of homogeneous particle diameters [7].In Salehi and Strom, it is mentioned that the particle size distribution also affects α p [5].In the case of ideal sediment, where all particles have the same size, C U would be equal to one.This was not the case, no matter which sediment we used.However, in Figure 9, it is shown that α p, f it /α p,analyt.is lower when using artificial sediments with lower C U values.However, as the data for "Ems 2012" show, C U is not the only parameter that leads to a decrease of α p, f it /α p,analyt. ."Ems 2012" with C U = 3.2 results in α p, f it /α p,analyt.= 4.47, which is greater than the values of quartz powder and bentonite.
Figure 10 presents the fit parameter C ADV with respect to C U .Our results show that C ADV is not only a device-specific parameter, which can be kept as constant for every material.We conclude that it is rather a combined device-and sediment-specific parameter.The data for C ADV again can be separated into artificial and natural sediments and were in the range of C ADV = 42-58 dB for quartz powder, bentonite, Chinafill, and metakaolin.On the contrary, C ADV data for natural sediments were in the range of C ADV = 6-39 dB.To summarize, our results show that the sonar model can be applied for different sediments without losing its qualitative behavior.We also showed that the model parameters clearly indicate the variation of sediment type.It was shown that the sediment absorption coefficients α p and C ADV not only rely on the sediment diameter, but also on the grain size distribution.

Practical Use and Reliability
In order to assess the practical use, applicability, and reliability of measurements in rivers and estuaries, for instance, we highlight some main characteristics of the method.
Due to the model's dependence on the particle grain size distribution, the user in the field would need to specify the characteristic sediment that occurs at the measurement location.When Urick's sound absorption coefficient is used, the mean particle radius and the solid density are needed.
Next, with respect to the sediment-specific parameters, knowledge about the measurement device (here: ADV) is necessary.In general, the sonar model can be also applied for ADCP measurements (e.g., [11,27]), though within this study, we only discuss ADV measurements.The ADV probe (Nortek, Vectrino Profiler) that was used has a dynamic range of 60 dB.In comparison, the Vectrino+ (Nortek) has a dynamic range of 25 dB.Both probes have an emitter frequency of 10 MHz.When repeating our measurements with the Vectrino+ probe, we could not find any trend between the SNR and the concentration.Therefore, it can be concluded that a certain dynamic range is of importance to resolve the SNR response with respect to the concentration.Furthermore, it is mentioned that the presented parameters are only valid for the Vectrino Profiler with a dynamic range of 60 dB.
Following the SNR function to concentrations lower than 0.1 mg/L, one will perceive that in this case, Equation ( 14) predicts negative SNR values.In particular for C = 3 × 10 −6 g/L, the function crosses the axis of ordinate at SNR = −10 dB.Following the principles of the sonar theory, the SNR has to be positive for particle detection [25].Negative values of the SNR would imply that noise is more dominant compared to the particle backscatter signal.Thus, the region of negative SNR values should be neglected for the interpretation of the model function based on sonar theory.Additionally, it can be observed that the experimental SNR data seem to reach a minimum threshold value for concentrations close to zero.This is explained by the SNR definition, which is implemented in the software of the measuring device.Within the ADV probe, the SNR is calculated by the difference of the backscatter signal intensity, which is here called Amp b (dB), and the noise intensity, Amp n (dB): Each of the two intensities is calculated by the received counts as: with counts max = 8000 for the used ADV probe [28].This formulation leads to a minimum intensity Amp min = −78.062dB for no detected counts and to a maximum intensity Amp max = 0 dB in the case of reaching the maximum counts.When inserting Equation ( 16) into Equation ( 15), the SNR is written as: where counts b includes the counts related to particle backscattering counts p , as well as the noise counts n .
In the case of no backscatter intensity (absolutely no scattering particles in the measurement volume), the minimum SNR would be SNR min = 0 dB and fully independent of noise intensity.Indeed, this is a theoretical threshold, because there are always some scattering particles in the measurement volume.Nevertheless, there is a minimum threshold value for the SNR, which cannot be reproduced by the log-linear sonar model.Assuming well-mixed conditions in the experimental tank, depth and bottom material are not expected to influence the SNR.As the name signal-to-noise ratio says, it describes the sound intensity of the ratio of the received signal and the background noise.In practice, just before measuring, the ADV probe records the noise intensity.Later during the active measurements, this intensity is subtracted from the absolute received signal (amplitude).Of course, depth and bottom material can influence the intensity of noise.However, the SNR is not influenced by this.This makes the SNR a reliable and universal measure for different environments.However, it has to be noted that the device's technique to determine the SNR is only valid for steady noise intensity.When conducting long-term measurements with frequent changes of noise (ships, wind, etc.) during the data acquisition, the SNR will be falsified.
In our experiments, we tried to avoid flocculation of particles while the suspension was well mixed.In nature, especially in estuaries, various floc populations are observed, e.g., [29].Sound attenuation on flocs is expected to be different than on single primary particles.Flocs are a connection of solids and organic substances, which vary in size, density, and elasticity.In order to calibrate the probe for flocculated suspensions, advanced experimental setups are needed where flocs will not be destroyed during the experiment.
Summing up, the sonar model is able to reproduce the experimental data for a number of different sediments.It shows the log-linear increase of SNR with concentration, as well as the vertex behavior and the decrease of SNR for higher concentrations.This particular behavior of the SNR (increase and decrease) expects attention and special treatment when applying such a model for concentration measurements in nature.One SNR value can be assigned to two concentrations, which is a disadvantage of the model and which has to be taken into account.It reduces the measurement range of this method either to the low concentration regime (positive slope of SNR) or to the high concentration regime (negative slope of SNR).While only using the ADV probe for concentration measurements, the user will not be able to assess the appropriate concentration.Especially in nature, the usage of a second concentration measurement device (e.g., OBS, LISST, tuning fork, etc.) would be necessary to specify the concentration regime [6,11,30].Then, it is possible to utilize an ADV to measure the velocity, the turbulence, and the concentration simultaneously.

Conclusions
In the scope of the experimental investigations, a clear relation was found between the SNR and the concentration of several sediment suspensions.The method was found to be valid for maximum concentrations C max = 50 g/L, as is shown for sediments from the Ems Estuary.In higher concentrations, correlation values between two received signals dropped below 70 %, that is why those measurement results were discarded within the evaluation [31].While combining scattering and sonar theory, an advanced model was presented for computing the concentration out of the SNR for natural and artificial sediments.The derived model parameters are only valid for the ADV probe and the sediments that we used in this experiment.However, it was shown that the methodology is valid for different sediments, both natural or artificial ones, and it can be applied for further investigations of turbulent suspension flows.Since this measurement technique allows temporal high resolved (f = 100 Hz) non-invasive concentration measurements, it is meant to be a basis for further experimental investigations of turbulence and diffusivity, either in nature or at the laboratory scale.

Conflicts of Interest:
The authors declare no conflict of interest.The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; nor in the decision to publish the results.

Figure 1 .
Figure 1.Sound intensities, which are of importance to state the sonar balance equation.

Figure 2
Figure 2 illustrates the experimental setup consisting of a rectangular box, which was B 1 = 0.96 m wide and L 1 = 0.80 m long.The box was mounted on a steel bottom.The water depth in the box was H 1 = 0.30 m.Two opposite installed stirrers guaranteed homogeneously-mixed suspensions.A Nortek Vectrino Profiler ADV probe (sampling frequency f = 100 Hz) was mounted on a dynamometer wagon and adjusted in a way that the probe was box-centered and able to be translated vertically.

Figure 2 .
Figure 2. The experimental setup is shown for calibrating the ADV's SNR and the concentration.The sketch is showing the stirrer and the ADV probe in side view, and the photograph is showing the experiment in top view.

Figure 3 .
Figure 3.A box plot evaluation is shown of the SNR data for three repeated measurements in a low concentrated quartz powder suspension.Only negligible deviations of the median of each beam are recognizable.
were calculated.It shows that beams 2 and 3 provided similar mean values µ 2 = 44.96dB and µ 3 = 44.61dB, which is in between µ 1 = 46.03dB and µ 4 = 43.08 dB.Nevertheless, when averaging the data of all four receivers, the equal deviations of beam 1 and beam 4 were eliminated.The standard deviations of the SNR distributions were in the range of σ = 2.5 -2.7 dB, which corresponded to the values obtained by the box plot evaluation.

Figure 3 .
Figure 3.A box plot evaluation is shown of the SNR data for three repeated measurements in a low concentrated quartz powder suspension.Only negligible deviations of the median of each beam are recognizable.

4 7Figure 4 .
Figure 4.The evaluation of the frequency of measured SNR data is described in histograms (Run 03).Histograms are given for each receiver (beam) with corresponding mean and standard deviation values.

Figure 5 .
Figure 5. SNR measurements are shown for two different measurement heights as a function of the quartz powder concentration.Measurement results of each receiver are summarized.For better interpretation, the curve is divided into three characteristic areas (I, II, III).

Figure 6 .
Figure 6.The measured SNR data are shown as a function of quartz powder concentration.The fit function is given for the evaluation of the sonar equation (Equation (14)).This model predicts a decreasing SNR in the high concentration regime and a log-linear development in the low concentration regime.

Figure 7 .
Figure 7.The measured SNR data is shown as a function of concentration for different artificial sediments: Upper left: Quartz powder; Upper right: Bentonite; Lower left: Metakaolin; Lower Right: Chinafill.The fit function is given for the evaluation of the sonar equation (Eq.14).This model predicts a decreasing SNR in the high concentration regime and a log-linear development in the low concentration regime.

Figure 7 .
Figure 7.The measured SNR data are shown as a function of concentration for different artificial sediments: (upper left) quartz powder; (upper right) bentonite; (lower left) metakaolin; (lower right) Chinafill.The fit function is given for the evaluation of the sonar equation (Equation (14)).This model predicts a decreasing SNR in the high concentration regime and a log-linear development in the low concentration regime.

Version December 12 ,Figure 8 .
Figure 8.The measured SNR data is shown as a function of concentration for different artificial sediments: Upper left: Lake Eixendorf: Upper right: Lake Altmühl; Lower left: Ems 2012; Lower right: Ems 2015.The fit function is given for the evaluation of the sonar equation (Eq.14).This model predicts a log-linear evolution in the low concentration regime and a decreasing SNR in the high concentration regime.

Figure 8 .
Figure 8.The measured SNR data are shown as a function of concentration for different artificial sediments: upper left: Lake Eixendorf: upper right: Lake Altmühl; lower left: Ems 2012; lower right: Ems 2015.The fit function is given for the evaluation of the sonar equation (Equation (14)).This model predicts a log-linear evolution in the low concentration regime and a decreasing SNR in the high concentration regime.

Table 1 .
The different particle diameters used to investigate the validity of Rayleigh scattering.

Table 2 .
The different particle parameters used to calculate the Ur-coefficient and the material parameter A.

Table 3 .
Fit parameter of the sonar model for different sediments.R 2 as a statistical parameter is the square of the correlation between the fitted and the measured values.

Π 1 (-) Π 2 (-) C ADV (dB) R 2 (%)
Artificial sediments have lower uniformity coefficients C U than natural sediments.This results in lower values of α p, f it /α p,analyt., which means that the agreement increases towards the formula of Urick.Artificial sediments have lower uniformity coefficients C U than natural sediments.As a result, the sonar model predicts greater C ADV values.