The assessment of an acoustic Doppler velocimetry profiler from a user’s perspective

Acoustic Doppler velocimetry profilers (ADVPs) are widely used in both experimental and field studies because of their robustness in velocity measurements. The acquired measurements do not only offer estimates of the local and instantaneous flow velocity at the interrogated measurement volume, but can also be further processed for the estimation of the bed surface shear stresses, thus they are finding a wide range of applications ranging from water engineering to geomorphology and eco-hydraulics. This study aims to evaluate the performance of an ADVP in obtaining hydrodynamics measurements under fixed flow conditions, with various probe configurations. To this goal, a robust search is conducted where ADVP probe settings are sequentially altered. A number of assessment criteria are used including qualitative observations, such as checking the shape of the velocity profile, as well as quantitative error metrics, including signal-to-noise ratio, correlations and number of spikes. Further, estimation of the bed shear stresses computed by means of using the log Law of the Wall and turbulent kinetic energy, allow obtaining a better understanding of the uncertainties involved and the importance of making a better informed choice with respect to the probe configuration settings. Thus, the methodology and performance metrics provided herein, although presented for a given flow, can generally be applied from practitioners and researchers alike.


Introduction
Acoustic Doppler velocimetry (ADV) is widely used as one of the most versatile and robust flow diagnostics tools for both laboratory and field studies across a range of research and applied themes spanning engineering eco-hydraulics and geomorphology. A range of specific ADV probes with varying specifications are readily available for use by professionals and researchers. ADV measures the sampling volume of a water body in the water flow via the Doppler shift caused by the reflection of the pulse transmitted from the centre transducer (the centre beam) and received by the four receiving beams (Nortek 2012;NortekAS 2015;Kraus et al. 1994). An acoustic Doppler velocimetry profiler (ADVP), i.e. Vectrino II, has the same working principle, but enables profile measurements, which may involve interrogating more than one measurement points, at the same time. Thus, ADVPs may typically have additional parameters that the user may need to define in the probe settings, which even though desirable, may leads to increased uncertainties if there is no generally acceptable guidance to enable their optimal use.
Many research works about the noise associated with acoustic signal of ADV have been done, such as Khorsandi et al. (2012), García et al. (2005), Anderson and Lohrmann (1995). It is suggested that Doppler noise might increase the value of flow hydrodynamics computed by time series of instantaneous velocities, such as turbulence intensities and turbulent kinetic energy (Anderson and Lohrmann 1995). Signal-to-noise ratio (SNR), as an important indicator of data quality, is defined as SNR = 20log 10 Amplitude signal Amplitude noise , and it is given the unit of dB (Nortek 2012). SNR sweet spot, where the measurement has the highest SNR is located 5 cm below the transmitter (Nortek 2012). At the same time, the correlations between the successive pings are suggested to be above 90% (NortekUSA 2011).

3
A wide range of post-processing methods have been developed for velocity time series to improve the performance of ADV. Some of these methods are numerical analyses of the time series, and some of them involve indicators of data quality like SNR and correlations (Goring and Nikora 2002;Parsheh et al. 2010;Hurther and Lemmin 2001;Jesson et al. 2013Jesson et al. , 2015. Most of these methods introduce removing or replacing the spikes. In fact, any post-processing algorithm may influence the integrity of the velocity representation. Compared to post-processing the velocity time series, enhancing the data quality during recording is also important to enhance the data quality (Voulgaris and Trowbridge 1998). To achieve this, some studies have identified the relation between Doppler acoustic seeding material and Doppler noise, as an inherent factor of this technique (Khorsandi et al. 2012).
Certain hydrodynamics measurements taken by ADVP (Vectrino II) are compared with other flow velocity measuring equipment, such as PIV and ADV, like Vectrino I, (Zedel and Hay 2011;Craig et al. 2010). Near bed measurements taken by ADVP may be influenced by the bottom echo, and may be associated with very high SNR, which can hardly reflect the flow velocity (Rusello 2012). Vectrino I has been used in recording flow velocities over a few decades. Vectrino II, as an updated version of Vectrino I, is more customised, and it enables the highest frequency of 100 Hz, which can better represent the flow hydrodynamics. However, in practice using ADVP under certain default configurations can easily result in obtaining flow diagnostics that are nonrepresentative of the real flow conditions. This appears to be true for most probe settings but even more so for those with which higher temporal resolution can be achieved. Vectrino II has a few indexes in the probe configuration, and different probe settings can influence the velocity result to some extents. In order to identify the best probe configuration, a number of velocity profiles are taken at the same location of a flume, under the same flow rate, but with different probe configurations.

Methods
Preliminary observations in using this probe revealed that there is a varying level of dependency on a number of parameters in the probe configuration, which are not detailed in the user manual, and lack a definite guide for the user. Subsequently users of this equipment may end up underutilizing or using it in a manner that returns inaccurate results. To this goal a series of laboratory experiments have been conducted, under the same open channel flow conditions, using a profiler (ADVP Vectrino II from Nortek®) aiming to cover the full range of probe configuration combinations that can be used in practice. Flow velocities have been measured via the software Nortek Multi-Instrument Data Acquisition System (Vectrino Profiler). The operation of the probe is configured via the configuration panel at the probe's proprietary software and is governed by a range of parameters ( Fig. 1). For this study a range of values for the configuration parameters are chosen to be tested, so that their impacts on flow diagnostics can be identified.
The experiment has been conducted in a recirculating glass walled flume with sufficient length (> 10 m) in the Water Engineering Laboratory at the University of Glasgow. A series of acrylic panels are placed on the side of the flume with an angle of 17° to simulate riverbank (Fig. 2). It is also shown in Fig. 2 that a vegetation array consist of some acrylic rods (with a uniform diameter of 6 mm, spaced 160 mm each other and configured linearly) are placed on the top of the riverbank via a series of support panels, starting from 4.5 m upstream of the measurement location. The width of the main channel is 1 m, and the flume bed consists of a few centimetres high layer of sand with nominal diameter from 0.5 to 2.36 mm (d 50 = 1 mm) covering the whole flume width (Fig. 2). The sand bed is designed to be fixed during all experiment runs. The velocity profiles have been taken at the same location (6 m downstream of the inlet, and 345 mm away from the sidewall; the measurement location is denoted by a red bar in Fig. 2) under the exactly same flow conditions, but different probe configurations. The flow rate is about 0.0105 m 3 /s, which is calculated from the average measurements of electromagnetic flow meters attached to the inlet pipes via which water is circulated in the flume. The flow depth is constantly equal to 142 mm along the flume length, which results in a mean flow velocity of about 0.060 m/s. The Reynolds number is about 7000, and the flow is turbulent and the turbulent boundary layer is fully developed at the measurement cross section.
The probe is joined with a vertical gauge, which enables the probe to move vertically. The measurement height can be measured by using the bottom check function of Vectrino II. For each experiment, single or profile (measuring multiple points at one time) measurements have been taken to reconstruct a velocity profile. While taking the velocity profiles, the probe has been adjusted to different heights so that the measurement locations correspond to the specific flow depth as shown in the velocity profile plots. Given that the flow assessed herein is statistically the same along the longitudinal direction, and given that the probe's geometry remains fixed, we expect there is no effect on the flow measurements. An effort is applied in keeping the orientation of the probe as consistent as possible and aligning vertically to the main flow direction. Some preliminary tests show that the velocity time series longer than 1 min has negligible errors in both mean velocity and standard deviation. So, depending on various frequencies, with the minimum measurement time is 1 min, about 6000 samples were constantly exerted in the velocity time series for all configurations. For instance, the probe records velocity for 1 min at a sampling rate of 100 Hz, and records velocity for 4 min at a frequency of 25 Hz. In addition, a velocity profile taken by ADV (Vectrino I) is also presented in this study, for reference.
The effect of tested parameters, including profile length (is shown in Fig. 3, adjusted by Range to first cell and Range to last cell), Range to fist cell, Sampling rate, Ping algorithm, Transmit pulse size and Cell size, on the sensitivity and the accuracy of the obtained results is assessed. The range of values that is chosen for the test indexes are given in Table 1. This initial assumption of the values of the parameters is made according to the user manual and after preliminary tests. Only one index is changed each time, which allows for checking the effect of each of those indexes on the velocimetry results. A naming system is devised by a combination of six digits, which denotes the value of the aforementioned parameters (Table 1). For instance, 213,211 means that a 10 mm long profile, which consists of 10 cells at one measurement, with both Transmit pulse size and Cell size of 1 mm measuring the velocity at 100 Hz, and it has Max interval Ping algorithm, and starts at 40 mm below the probe centre beam. Those with 'Single' as Profile length indicates that the probe measures only one point at one measurement, which is defined as point measurements, and its actual length is dependent on the length of the Cell size (cell size and pulse length are usually the same for two pulse-coherent method, but the larger cell sizes can be obtained by averaging the pulse length returns at different depths). The elevation of the middle point of the sampling volume denotes the height of the velocity measurement.
The experiment has been conducted in two phases. In the preliminary phase, the maximum and the minimum value of each index are chosen, respectively, to consist of 48 profiles. Because the restrain due to the water depth, probe configurations with profile length of 30 mm and Range to first cell of 50 mm can only take two profiles, and these probe configurations are not experimented. After acquiring the velocity profiles in this phase, a visual check has been done (details can be seen in the next section). The result clearly showed that Adaptive Ping algorithm always lead to random results (extremely high errors), where the detail can be seen from Table 2. Similar findings also can be found in (Rusello 2012). So the Ping algorithm should be fixed to be Max interval in the refined phase of the experiment. With additional considerations of the results from the preliminary phase, Range to first cell can be fixed to 40 mm, because very small variations have been found in this phase, and Transmit pulse size and Cell size are changed correspondingly. Then, in the refined phase, other values of the index are tested in a similar way, which consisted of another 12 probe configurations.

Results
All results are processed to velocity profile analysis, including visual check and fitness check to a standard logarithmic function. The results of velocity profiles with qualified logarithmic shape are then processed to shear stress and turbulent kinetic energy analyses. Some spikes occur in the time series due to the inherent problems of the measuring technology or insufficient seeding material for specific probe configurations. So all velocity time series obtained in this study are despiked as it is described in (Goring and Nikora 2002), and the identified spikes are removed. The impact of spikes on the result is discussed separately.

Velocity profiles
By visually checking the velocity profile plots, the results obtained from the preliminary phase are identified as 'smooth' (Fig. 4a, b), 'curved' (Fig. 4c) and 'random' (Fig. 4d). It can be seen that all measurements with Adaptive Ping algorithm are identified as random. The influence of echoes from boundaries can affect the result up to approximately 10 mm above the sand bed (up to the 5th point from the bed in the profile). Different from smooth and curved profiles, at the near bed region of random profiles, the velocity varies a lot (111144 and 113141 in Fig. 4d). Even from observation the magnitude of velocity for the smooth profiles varies up to 25% (Fig. 4a, b). Figure 3c shows that many profile measurements result in curved velocity profile shape in each profile measurement. For example, experiment 411241 and 411214 ( Fig. 4c) only show an abnormally small velocity at the near bottom portion of each profile, but experiment 411211, 413211 and 413241 ( Fig. 4c) show the highest velocity at the sweet spot (50 mm below the centre beam, Fig. 3) and abnormally smaller velocity at distal to the sweet spot ( Fig. 4c).
Some additional quantitative analyses of the velocity profiles by assessing their fitness to logarithmic functions are also taken to assess the quality of the velocity profile. Velocity profiles are plotted in log scale to see how they fit a line of u = p 1 x log (z) + p 2 where z denotes the elevation of velocity measurement, u denotes the time-averaged velocity, p 1 and p 2 identify the slope of the function and its intersection to the axis, respectively. The first five points from the bed are not taken into account for analysing, because the result are influenced by the strong echo from the boundary as it has been mentioned before. Some indicators are involved in the function fitness check stage, including sum of squared errors (SSE, defined by Eq. 1), coefficient of determination ( R 2 , defined by Eqs. 2 and Table 1 Naming system used in this experiment *'Profile length' is basically the value of Range to last cell minus the value of Range to first cell, 'Single' indicates point measurements. The actual length may differ with the number shown here for profile measurements, and it will be the nearest integer number, which is smaller than the displayed value  where K denotes the total number of velocity points in the velocity profile; u i denotes the time-averaged velocity; û i is calculated from the best-fit logarithmic function by using the corresponding elevation. R 2 is defined by Eqs. 2 and 3: where U denotes the profile mean velocity, n is number of samples in each velocimetry record, and K is the total number of velocity measurements in the velocity profile. R-squared closer to 1 means that the velocity profile has a better regression fit to its best-fit logarithmic function. RMSE is defined as: RMSE describes the error in addition to SSE with considerations of the number of velocity points in the velocity profile. Some errors are extremely high because of those random results ( Table 2). The median value of each error indicator is chosen as the critical value. In other words, the first 50% among all profiles from low error to high error for the specific indicator are processed for the next step. Those profiles meet the critical value of all these three error indicators are identified as good-fit results. Profile names marked with '*' indicated the good-fit profiles, which are plotted in Fig. 4a.
It can be seen from Table 2, that all good results from logarithmic function check are parts of 'smooth' profiles identified by visual check. The results of the preliminary phase show that the average SSE, R 2 and RMSE are 1.75, 0.15 and 0.26 for all Adaptive Ping algorithm, and are 1.12, 0.62 and 0.15 for all Max interval Ping algorithm. R 2 as an indicator of the regression of the profile shape to a standard Random velocity profiles are also observed in Thomas et al. (2017). Even among the good-fit results, the values of y 1 and y 2 vary a lot, which means the logarithmic shape that they follow is different, which can significantly influence the near bed shear stress estimation, for example done under the log Law of the Wall assumption. However, both the magnitude and the shape of the stream-wise velocity profiles with different Range to first cell are generally the same (i.e. 111211 and 121211; 113241 and 123241 have similar p 1 and p 2 ). In the refined phase is fixed to 40 mm below the centre beam, which reduces restraints due to the flume geometry, and enables more measurements. According to the instruction by the probe manufacturer, the Transmit pulse size and the Cell size should normally correspondingly change, and no clear impact of the correlations between these two parameters is observed in the preliminary experiment phase, so in the refined phase, these two indexes are changed correspondingly to the same value.
Most velocity profiles obtained in the refined phase are identified as smooth by visual check, except experiment 212222 and 311222, where clear gaps between the two distal ends of each profile measurement can be observed (Fig. 5). The strong echo from the boundary causes near zero result until about 0.01 m above the bed for most measurements. Larger cell size may lead to longer boundary affected region (e.g. profile 311244 and 312244 in Fig. 5), as the midpoint of the cell is considered as its elevation in this plot. The velocity profile obtained by Vectrino I is also smooth, but it clearly follows a different logarithmic shape comparing with those of most good results taken by Vectrino II. Similar function checks are applied to refined phase as well, and the results are shown in Table 3. The overall quality of results for the refined phase is much higher than it for the preliminary phase, so the same reference values of SSE, R 2 and RMSE used in the preliminary phase are applied, and the results are given in Table 3.
Most good-fit profiles by function check are smooth, unless experiment 211222 and 311222. Though these profile are good-fit of the logarithmic function, it can hardly use in research studies, because clear biases exist between each profile measurement. From a view of all good-fit results among the whole test, p 1 and p 2 vary from approximately 19-75, and from approximately − 6 to − 2.6, respectively. These indicate that both the logarithmic function and its intersection to the axis may have quite big variations, which can be reflected in the magnitude of the shear stress estimated by the log Law of the Wall and the computed   (NortekUSA 2011). The value of SNR for the Vectrino I is constantly above 18 dB, which is above the threshold of it suggests by many research studies (Khorsandi et al. 2012;McLelland and Nicholas 2000). However, such guidance for point measurements taken by Vectrino II has not been made available yet. The result also shows that SNR and correlations are not independent of probe configurations. In this experiment, the flow condition remains exactly the same, but the mean SNR varies from 6 to 30 dB, and the mean correlations vary from 58 to 96% (Tables 2 and 3). SNR and correlations indicate the noise level of the measurement, lower SNR and correlations may result in an increase number of spikes of the velocity time series. However, good SNR and correlation do not have direct link to good quality of results (i.e. experiment 111111 and 121114).

Near bed shear stress
Many research studies have used velocity profiles obtained by either ADV or ADVP to estimate the near bed shear stress by using the log Law of the Wall ( log ). For the same bed roughness and same flow condition, as in this experiment, log should be same for all experiment runs. The impact of various probe settings, based on the profiles that are good-fit to logarithmic functions and log are presented, which can also reflects the shape of the velocity gradient of the profile. Equation 5 is the method introduced in Sumer (2007) which involves Nikuradse's equivalent sand roughness ( k s ): where is the Von Karman constant ( ≈ 0.4 ), z denotes the height of the velocity measurement; assuming z � ∕k s = 0.2 -this value is chosen within the range because it can cover more measuring points in the logarithmic function (Sumer 2007), z ′ denotes the distance between the tops of the elements comprising the bed surface and the location where the bed shear stress is zero; u * denotes the shear velocity. Equation 5 can be rearranged to Eq. 6: which can be used to acquire u * and k s . Similarly, the 5 points close to the bed are not taken into account. Thus, both k s and u * can be worked out, and then log can be estimated as: u * = √ log , where is the water density. According to this method, both k s and log for all 'smooth' profiles that are good-fit to a logarithmic function are estimated. k s varies from sub millimetres to 8.7 mm (Table 4). Even the velocity profile is smooth, different probe settings can still influence k s and log estimations a lot. In this method, the values of u * and k s may change correspondingly, and both of these two results are going to be assessed. k s has been empirically suggested a range of its magnitude, which is well recognised from 1 to 5 times of the median size of the bed material (Sumer 2007;Camenen et al. 2009). The ratio of k s and d 50 is presented in Table 4. Generally, k s obtained by profilebased measurements is higher than it computed by pointbased measurement data. Turbulent kinetic energy ( k t ) can be related to near bed shear stress as well. k t is defined in Eq. 7: where u ′ , v ′ and w ′ denote the fluctuation component of velocity time series of streamwise, transverse and vertical directions, respectively, , u is the time-averaged velocity, u i is the instantaneous velocity. k t estimation involves u ′ , v ′ and w ′ , which may have higher requirement of the data quality (Voulgaris and Trowbridge 1998). k t is estimated by using the time series of the three components of the local flow velocity at z ≈ 0.01 m ( z∕h ≈ 0.07 ), where boundary effects are avoided, and k t reaches its maximum value (Biron et al. 2004). It can be approximated as k t = C b u 2 , where C b is a coefficient relevant to bed roughness, and U is the depth mean velocity (Yang et al. 2016). C b has been estimated by Yang et al. (2016) through both experiment and theory, the obtained values are 0.02 and 0.025, respectively, for the sand bed material with sizes between 0.6 and 2 mm, which is similar to this experiment setup. Among all good-fit velocity profiles, the mean flow velocity estimated by the velocity profile have negligible variations, so U is constantly equal to 0.063 m/s in this study. The bed shear stresses can then be estimated as (Biron et al. 2004;Yang et al. 2016). Turbulent kinetic energy estimated by C b has been calculated from both of these two values, 0.025 u 2 (= 0.019 Pa) and 0.020 u 2 (= 0.015 Pa), and they highly match log estimated for good logarithmic fit profiles (Table 4). Probe configurations that meet both k s and k t check are marked with * in Table 4.

Physical context and interpretation of the effects of the probe settings
It is generally important to note that as a first consideration, one should always keep in mind the spatiotemporal scales of the relevant flows and the range of fluctuations of the assessed variables, towards making a selection for a velocimetry probe. The research theme of geophysical flows measurement and instrumentation thrives with the implementation of a range of technologies from optical (such as laser Doppler velocimetry or particle image velocimetry) to acoustics and ultrasonics. Consideration of the flow scales and the physical application to being examined will help assess the best technology for each case. As an example, for the low mixing, wall bounded steady flow tried herein ( U ≪ 1.2 m/s), the Vectrino II ADVP has been used. It is noted that for applications with relatively higher flows (e.g. U ≫ 1 m/s), different velocimetry probes may be preferable. Even though the production of this ADVP has been discontinued, the similar experiments can be conducted with more modern probes (such as Ubertone's ADVP). It is also important to note that bi-static Doppler Velocimetry profilers using a narrow-band pulse-coherent technology are limited by the maximum velocity profile range. Broadband incoherent technologies overcome this issue (as testing of more recent acoustic probes from Ubertone has demonstrated, e.g. Conevski et al. 2020). Despite these apparent limitations, the results and observations presented herein are useful both in terms of the general search framework applied towards validating produced ADVP data using sound first principles, as well as a best practice guide for current users of the Vectrino II ADVP probe.
The spatiotemporal resolution and accuracy of point or profile measurements depend on the shape and size of the measurement volume, which in turn is defined by the orientation of the acoustic emitting beams (which are fixed for a given probe) and the probe settings being selected. Vectrino II has a projected beam that becomes wider with the distance from the face of the emitter. This implies that profile measurements at cells more distant from the transducer will be contaminated (due to the receiver probing velocity from adjacent cells), resulting in relatively more erroneous measurements. This is avoided from recently produced devices that have a narrower beam and wider receiving beams, resulting in more focused velocity profiles, and minimising the velocity contamination (see for example the Ubertone's ADVP, Conevski et al. 2020). The above consideration may be used in an attempt to explain the "banana"-like curved shape of the velocity profiles in Fig. 4c. Considering that for the flow examined herein the sweet spot is at about 50 mm from the face of the probe, Fig. 4c demonstrates capturing velocity profile measurements that deviate more from the logarithmic profile the more distant from the sweet spot they are. Generally, good results are found for sampling rates 25-50 Hz (e.g. Figure 4a), and even single point 1 mm velocimetry measurements at the sweet spot are found to be much better quality than the profiler results (Fig. 4c). Seeding of the flow should be sufficient so that the signal-to-noise ratio (SNR) will be around the range of 30-40 (especially around the sweet spot location). However, it needs be noted that the SNR was found to be highly dependent on the probe configuration. Thus, as a probe configuration is selected as being the optimal for the flow the seeding needs be checked to ensure sufficiently high SNR. In the experiments studied herein the flow is steady and boundary is fixed. However, for transient processes, such as seen in moving boundaries, e.g. over a wave induced flow field (Gaeta et al. 2020), high SNR and correlation may be indicative of a reflective boundary, implying that other means need be sought after (such as to analyse particles scattering and flow velocity spatiotemporal patterns in parallel with videography), and that high correlation is not an imperative goal (especially in cases where flow turbulence may be such where the acoustic pulses get to become uncorrelated).
Regarding the choice for the ping setting, it is thought that Adaptive ping may be preferable for cases of long-term deployment where the flow conditions may be changing or where there is no control or knowledge of the boundary of the flow. This works by sending relatively short space pings to estimate flow velocity and then trying to reduce the ping interval to match the velocity range, while considering the distance to the solid boundary. Minimal interval is a much preferable setting, where many more pings are sent over the same time period. Max interval is matching the interval of the pings to the velocity range of the measurements, using very long gap between the pings that gives the narrower velocity range that can be measured.

Dependance of the flow metrics on probe settings
As it is shown in Fig. 4 and Table 1, some experiments with high SNR and correlations but result in random velocity profile, and some results with relatively lower SNR and correlations provide more reasonable velocity profiles. For example, experiment 113114 has high averaged SNR (about 29 dB) and high averaged correlations (about 95%), however even the time series looks smooth and few spikes (Fig. 6a) and made the velocity profile random (Fig. 4). From those good-fit results, though the mean correlation remains above 85%, the mean SNR varies from 12 to 16 dB for some experiments (see Tables 2 and 4). The mean SNR values calculated with inclusion of the measurements at the near bed region, in which the echo from the boundary can increase the value of SNR, which means the averaged SNR for some points is even smaller. Therefore, SNR and correlations cannot fully reflect the data quality. Table 4 also shows the mean percentage of spikes of all measurements in the profile with the specific probe configuration. Low SNR and correlations can result in the high percentage of spikes (Table 4), but they have negligible influence in the mean flow velocity calculated from the despiked velocity time series. All velocity data used in this study are after running the despiking code, which has little impact on the mean velocity calculation and the relevant analysis, such as velocity profile assessment and estimation of log . For instance, experiment 113241 has the maximum percentage of data points defined as spikes, but the averaged difference of the mean velocity is about 2 × 10 −7 m/s (about 0.0003% of the mean velocity) before and after despiking (Fig. 6b), and for others among results that good-fit to a logarithmic function (e.g. Fig. 6c, d), this difference is much smaller.
Sampling rate is another key factor that affected the number of spikes, and the result shows that higher frequency of sampling rate can cause relatively more spikes in the velocity time series (Table 4). However, sampling rate does not have direct link to SNR and correlations, and this study also shows that the percentage of spikes is more sensitive to sampling rate. For example, experiment 111241 and 113241 contain 4.1% and 7.2% spikes, respectively, in which SNR and correlations of them are both small (see Table 2). In general, the number of spikes can be largely reduced by adding more seeding material (Khorsandi et al. 2012;. In this experiment, the flow condition remains exactly same, so the percentage of spikes reflects the required minimum seeding concentration in the flow. Meanwhile, the measured k t is independent of the percentage of spikes. Both experiment 113241 and experiment 123241 have 7.2% spikes and similar SNR and correlations conditions, but more negative results can be observed in experiment 113241 (Fig. 6b, c), which do not significantly affect the mean velocity and velocity profile gradient ( log for both are similar), but TKE for 113241 is about three times of the TKE for experiment 123241. The processed results for both experiment 123241 and 312244 have reasonable distribution of velocity (Fig. 6c, d), but the magnitude of k t varies a lot, and the former is closer to t . This indicates that metrics of SNR, Correlations (CC), and percentage of spikes while being indicative of the signal quality, they may not entirely be sufficient in representing how well the flow field is being sampled to produce reliable flow readings.
Both log and TKE vary a lot, even they are only computed by the velocity profiles that are good-fit to a logarithmic function. For specific configurations of Vectrino II (like experiment 111241), TKE is a few times higher than other probe configurations (such as experiment 111244 and 312244). The highest log is more than three times of the lowest. Comparing to the influence of spikes on hydrodynamic 1 3 variables, the impact of probe configuration is more significant. Thus, it is even more important to record velocity data with a good probe configuration than developing post-processing algorithms, because probe configurations can also significantly influence the mean velocity profile in addition to other computed results by using instantaneous velocities.
The time series of the flow velocity at a certain location could be obtained using acoustic Doppler velocimetry based on the acoustic principle. However, the instantaneous flow velocity readings are highly dependent on different parameters from the seeding to the probe configurations, depending on the flow conditions, as it has been shown in this work, which results in having uncertainties in the local flow velocities. The field of acoustic Doppler velocimetry has certain advantages over other velocimetry tools and one of its applications is the estimation of bed shear stress for the assessment of the destabilization of bed surface. The bed shear stresses estimated using the flow velocity measurements is used to estimate the Shields stresses and hence assess the conditions that can result in sediment entrainment using the Shields criterion, or the Shields graph. However, there is a huge uncertainty in the assessment of riverbed destabilization using the Shields criterion due to the certainties discussed previously, the fact that the estimated bed shear stresses are dependent on the methodology used for their estimations (Yager et al. 2018), and the fact that the Shields criterion for sediment entrainment involves up to two orders of magnitude variability. Therefore, the use of other novel technologies, such as inertial sensors and instrumented pebbles (Al-Obaidi et al. 2020;Al-Obaidi and Valyrakis 2021a, b), allows for more robust assessment of bed surface destabilization directly and more efficiently than using Doppler-based flow velocimetry in addition to assessing the distribution of energetic flow structures sufficient Fig. 6 Distribution of unprocessed and processed instantaneous velocity time series at about 1/3 of the flow depth for the experiment with: a high SNR and correlations; b low SNR and correlations; c more spikes identified among the results well follow logarithmic functions; d less spikes identified among the results well follow logarithmic functions to mobilize sediments (Al-Obaidi and Valyrakis 2021a, b;Valyrakis et al. 2013).

Optimal probe configurations
No matter what changes for other indexes, all results show that Max Ping interval always made the SNR smaller than it with Adaptive Ping interval ( Table 2). All profiles that fitted logarithmic functions follow the trend that for point-based measurements have larger values of p 1 and smaller p 2 comparing to those profile-based measurements. That results in both smaller log and smaller k s . Cell size and Transmit pulse size should be treated carefully. If take the height of the midpoint of Cell size as the elevation of the measurement, the bed boundary affected region can become longer, however this is not clearly delivered in the user manual.
The four best probe configurations of Vectrino II are selected after assessing k s , log and TKE , and the velocity profiles are plotted in Fig. 7, which consists of two pointbased measurements (experiment 111222 and 112244) and two profile-based measurements (experiment 411244 and 211222). The shape of the velocity profiles obtained by ADVP with these four probe configurations does not have big variation in p 1 and p 2 shown in Table 3. By using any of these four probe configurations, the recorded and computed flow hydrodynamic variables, such as mean velocity, log and k t well agree to each other, which also shows a level of accuracy and consistency of Vectrino II. This study shows a high level of dependency between probe configurations and velocity results, so depending on the type of analyses by the measured velocity, a test should be done to identify a proper probe configuration for the given flow condition before using Vectrino II. According to the aforementioned findings, Ping algorithm should be Max interval, and Range to first cell can be fixed either 40 mm or 50 mm, because it will not significantly affect the results. It is always preferred to choose high Sampling rate in many experimental works, however, a despiking algorithm may need to be run for configurations with high Sampling rate, which may cause more spikes along the unprocessed velocity time series. In addition, a note should be made that increasing either Cell size or Transmit pulse size may increase the length of boundary affected zone.

Framework for assessing optimal probe configurations
To identify an optimal probe configuration for near bed shear stress estimation under the given flow, the bed boundary influenced region should be detected first. Measurements within this region always have near zero velocity measurements in all directions and very high SNR with small fluctuations along the whole time series. At this region, the values of SNR for the four receiving beams do not match each other, and one or some of them are clearly higher. Measurements within this region do not reflect the actual flow velocity. Then, the quality of velocity time series needs to be evaluated. The time mean velocity after running an appropriate despiking algorithm should be almost the same as the raw mean velocity. Alternatively, the difference of the root-mean-square error between it for velocity time series processed once and twice should be negligible. The quality of the velocity profile can be assessed by comparing it with standard logarithmic functions. The smaller error indicates the higher quality of the profile. log and TKE can be estimated separately, and the difference between log and TKE is empirically suggested to be constantly smaller than 0.5 log + TKE . Benefits can also be obtained from comparing the computed result with the theoretical estimation of the near bed shear stress, if this is available.
Velocity profiles can be taken at the same location with different probe settings (choose a range of the values of profile length, Sampling rate, Cell size and Transmit pulse size depending on the aim and the requirement of the experiment, and test them). Visual check and check of fitness to a logarithmic function should be done. Analyses by using time mean velocities, such as depth mean velocity and shear stress estimated by the log Law of the Wall, can be achieved then. For analyses involving instantaneous velocities in time series, an additional analysis of turbulent kinetic energy (like comparing to the values of log and TKE ) is needed, and can follow the process demonstrated in this study. The percentage of spikes in the velocity time series does not significantly affect the time mean velocity, however it may influence the instantaneous velocity analyses in some extents, which can also be reflected by TKE . So, shear stresses estimated by Fig. 7 Plots of velocity profiles for the best probe configurations assessed in this study. Because of the boundary effects, the profile starts from the sixth point from the channel bed different methods with proper probe configurations should not vary a lot, and these values should well agree to them estimated from the theory.

Conclusions
Velocity profiles are taken at the same location and under the same flow condition by Vectrino II with different probe configurations. Although the production of the specific tested ADVP Vectrino II probe has been discontinued, these are still in use from a plethora of practitioners and researchers alike. Thus, it is expected that the practical findings regarding the specific probe configurations presented herein may still be of value, while caution is needed in using the proper velocimetry probe and settings for different flow dynamics and flow measurement circumstances. Analysis of the experimentally obtained results shows that the time-averaged flow velocity, log and TKE vary significantly (the measured time-averaged velocity result can vary more than 50% of the mean flow velocity estimated from the flowrate, see Fig. 4; and the variation between the highest and the lowest shear stress can be up to 8 times, see Table 4), which shows that the sensitivity of the flow diagnostics to the probe configuration parameters needs be taken into account, before any measurements or estimations are attempted. The assessable metrics of data quality, SNR and correlations, suggested by the probe manufacture can only reflect the quality of velocity time series in certain level, in terms of number of spikes, but may not result in consistent measurements of the flow hydrodynamics. SNR and correlations are also sensitive to probe configurations. Given the inability of these metrics to offer a comprehensive understanding for the suitability of the probe configuration under a range of various settings, the use of a set of checks, flow metrics and estimated parameters, is demonstrated in assessing the optimal configuration. Specifically, comparing the shape of the obtained velocity profile to the expected (e.g. logarithmic, for wall bounded flows) can be performed by visual check. In addition to this, goodness of fit metrics to the expected theoretical velocity profile function can be also used. Computed values of TKE from the instantaneous velocities, as well as values of log and k s estimated from the velocity profiles, can be used to demonstrate the importance of carefully considering the optimal probe parameters. Some good probe configurations (for example those with sampling rates from 25 to 50 Hz, and single point measurement vs profiler measurement) are selected for the flow conditions assessed herein experimentally and the framework and metrics used to define these is presented. Certain trends can be observed from the analysis of the experiment results, consistent to the recent literature.
Adaptive Ping algorithm always leads to random velocity profiles (see Fig. 4d) and not preferable for normal wallbounded flows with fixed boundaries. Min interval ping is preferable for flows up to 1.2 m/s. Point-based measurements always associate with lower SNR, and they result in smaller log and smaller k s . Larger Cell size can make the boundary affected region longer. Overall, the probe's configuration parameters plays an important role on the accuracy of data obtained from Vectrino II, and this study offers a framework and metrics for assessing their optimal selection, following a set of focused pilot tests before the acquisition of experimental data for a given flow, as demonstrated in this study for specific flow conditions. Funding This research has been supported by the Royal Society (Research Grant RG2015 R1 68793/1), the Royal Society of Edinburg (Crucible Award), and the Carnegie Trust for the Universities of Scotland (project 7066215).

Data availability
The data that support the findings of this study are available on request from the corresponding author.

Conflict of interest
The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.