Next Article in Journal
Toward an Automatic Quality Assessment of Voice-Based Telemedicine Consultations: A Deep Learning Approach
Previous Article in Journal
Linear Matrix Inequalities for an Iterative Solution of Robust Output Feedback Control of Systems with Bounded and Stochastic Uncertainty
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Towards Characterizing and Developing Formation and Migration Cues in Seafloor Sand Waves on Topology, Morphology, Evolution from High-Resolution Mapping via Side-Scan Sonar in Autonomous Underwater Vehicles

1
School of Information Science and Engineering, Ocean University of China, Qingdao 266000, China
2
Key Laboratory of Marine Geology and Metallogeny, Ministry of Nature Resources of People’s Republic of China, Qingdao 266061, China
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(9), 3283; https://doi.org/10.3390/s21093283
Submission received: 22 March 2021 / Revised: 1 May 2021 / Accepted: 6 May 2021 / Published: 10 May 2021
(This article belongs to the Section Remote Sensors)

Abstract

:
Sand waves constitute ubiquitous geomorphology distribution in the ocean. In this paper, we quantitatively investigate the sand wave variation of topology, morphology, and evolution from the high-resolution mapping of a side scan sonar (SSS) in an Autonomous Underwater Vehicle (AUV), in favor of online sequential Extreme Learning Machine (OS-ELM). We utilize echo intensity directly derived from SSS to help accelerate detection and localization, denote a collection of Gaussian-type morphological templates, with one integrated matching criterion for similarity assessment, discuss the envelope demodulation, zero-crossing rate (ZCR), cross-correlation statistically, and estimate the specific morphological parameters. It is demonstrated that the sand wave detection rate could reach up to 95.61% averagely, comparable to deep learning such as MobileNet, but at a much higher speed, with the average test time of 0.0018 s, which is particularly superior for sand waves at smaller scales. The calculation of morphological parameters primarily infer a wave length range and composition ratio in all types of sand waves, implying the possible dominant direction of hydrodynamics. The proposed scheme permits to delicately and adaptively explore the submarine geomorphology of sand waves with online computation strategies and symmetrically integrate evidence of its spatio-temporal responses during formation and migration.

1. Introduction

The submarine sand wave, one type of ubiquitous geomorphology distributed widely in the ocean, is highly correlative with multivariate marine environmental factors in their formation and migration process [1,2,3,4]. It is of great research value to understand sand wave dynamics in both marine science and engineering, such as the construction assessment of submarine pipelines and oil platforms [5,6]. To quantitatively identify sand wave variation in its topology, morphology, and evolution over time could symmetrically integrate evidence of its spatio-temporal responses to manifold intrinsic and external indices and help provide opportunities to explore the underlying patterns and mechanisms of their formation and migration.
Generally, submarine sand waves exhibit at a variety of scales, from small-scale megaripples to large-scale sand banks [1,7], which is actually a reflection of the variation under the influence of hydrodynamic forces [8], closely relevant to the coupling impact factors, such as the tidal characteristics, sediment grain sizes, the related residual currents, and so forth [9,10,11,12]. For decades, researchers have been working on topics with regard to the formation, detection, migration, simulation for submarine sand waves [2,5,13,14,15,16,17,18]. The existing techniques could be roughly divided into the numerical methods [5,13], mathematical modeling [19,20], field survey [21,22], remote sensing [18,23], and so forth. The classical numerical approaches could be performed by simple formulas and convenient calculation [13]. However, such empirical calculation tends to be composed of multiple parameters that must be considered, which would be sometimes very difficult to derive, and thereby possibly restrict the accuracy and applicable ranges; hence, often a rough estimation is used. The mathematical modeling could explain the generation and predict the migration of sand waves in an efficient and economical way [2,15], but no models could be universally applied. Remote sensing provides one highly efficient and dynamically observable means to conduct large-scale regional surveys of sand waves [18], but might lack showing morphological details in sand waves at small scales. The field surveys mainly employ the acoustic equipment, such as the multibeam bathymetry system, side scan sonar, sub-bottom profiler, and acoustic Doppler current profiler (ADCP), and so on to conduct sand wave surveys on site, which enables more specific investigation in the target underwater region at a higher accuracy [5,14,22]. However, usually this comes with a high cost as well as restrictions of availability by the marine environment conditions and submarine topography.
Recently, there has been increasing attention paid to the deployment of Autonomous Underwater Vehicles (AUVs) and Remote Operated Vehicles (ROV) in the field survey of submarine sand waves [24,25], e.g., Bluefin [26], Hugin [27], Remus [28], Iver2 [29], JAMSTEC [30], and so forth, which could greatly break through the restrictions of the traditional field survey driven by research vessels on a planned voyage, with more flexibility and adaptivity in applications. Underwater vehicles could be equipped with a series of sensors to complete a variety of sand wave inspection tasks and proceed at a lower speed and at a shallower depth towards the seafloor, allowing the morphological imagery at a higher resolution. AUVs are capable of carrying out a lawnmower pattern, i.e., a set of overlapping strips, which permits submarine geomorphology delicately across a certain range of underwater region rather than only one long narrow swath.
Among those categories of acoustic sensors, the side scan sonar (SSS) is one of the most widely mounted sensors on AUVs for sand wave inspections, due to its portability and imaging capacity. When the AUV proceeds along its track, the SSS emits acoustic waves to construct seafloor mapping for sand waves with the varied echo intensities. The direct use of SSS imaging could basically bear a fairly high coverage rate of seafloor mapping with time. However, acoustic echoes from SSS could be very easily interfered by various kinds of environmental noise, reverberation noise, instrument noise during their propagating process in the ocean, there will be a lot of high-frequency noises included in the waveform of sand waves, which makes it quite challenging to distinguish functional signals. Synthesizing consecutive uniform pings of echo intensities and converting these into SSS imaging might initially bring possible distortion, stretches, and compression to some extent, especially when the attitude of AUVs (roll, pitch, yaw) changes over time. The configuration of relatively poor real-time computational abilities in most AUVs may restrict online detection for sand waves, so in most cases, it still remains difficult to identify sand waves and explore the morphological and topological parameters on site from SSS imaging collected by AUVs simultaneously with high accuracy and efficiency. Hence, here we make an attempt to utilize echo intensities directly derived from SSS to help AUVs accelerate the online detection and localization process for sand waves.
So far, machine learning has been regarded as one of the most powerful models in exploring seafloor morphological and topological properties that could discriminate and approximate nonlinear function without prior assumptions. All kinds of advanced algorithms, particularly the deep learning framework, have been developed and led to progress in the understanding of SSS imaging [31,32,33]. However, the use of such deep learning algorithms could raise higher requirement to the super-computing capacities, such as high-performance GPUs, and might suffer from insufficient learning, especially when new observations arrive progressively [34,35,36]. Recently, the extreme learning machine (ELM) has attracted more and more attention by providing comparable performances with simplified architecture at a much faster speed [37,38]. The online sequential extreme learning machine (OS-ELM) inherits the merits of ELM and bears an online learning ability when observations come in succession [39], which is more suitable for few shot learning in sand wave detection and overcomes the incomplete or unstructured data acquisition by AUVs [40,41].
Meanwhile, several previous studies [42,43,44,45] have proved that the sand wave attributes reflected from the acoustic cross-section in SSS imaging could be utilized as the morphological and topological indicators of sand wave dynamics during the evolution process of the seabed load sediment transport under hydrodynamic conditions. The asymmetric shape and orientation of sand waves would reveal the ebb flow dominance in the navigation channel, and sand wave size and spacing have been analyzed, both in natural environments [46], in laboratories [47], and through morphological modeling [19]. The hydro-sedimentological conditions could be also inferred with the help of morphometric relationships between sand waves. Wynn et al. [48] used a combination of high-resolution SSS and seafloor photography to image barchan dunes and sand ripples in the Faroe-Shetland Channel and have shown that high-resolution imaging of seafloor bedforms could provide some indication of the bottom-current orientation and flow velocity, and the bottom currents are modified by dune morphology. Dijk et al. [49] determined the morphology and dynamics of sand waves from time records of SSS and multibeam imagery and found that the coastal and offshore compound sand waves were nonrelic forms that actively develop and migrate. Passchier et al. [50] investigated the influences of seasonal storm events and fair-weather periods on the stability of large-scale compound sand wave forms by SSS off the Dutch coast, and they found that multiple seasonal storm events of low intensity did not have a measurable effect on sand wave morphology or position. Crawford et al. [51] extracted seabed sand wave orientation, wavelength, and defect density parameters by a novel image processing technique based on fingerprint analysis from SSS.
In this paper, starting from the inspiration to quantitatively investigate the geomorphology in sand wave profiles along a AUV’s track, we directly extract the echo intensity of each ping from SSS imaging through time varying gain (TVG) correction, speed correction, blind zone removal, and ensemble empirical mode decomposition (EEMD) denoising to compensate the frequent motion variation in an AUV as well as the noise interference, energy attenuation, and multi-path effects over time. Both ELM learning and some classical approaches have been jointly employed and aim to potentially provide a mutually reinforcing and complementary, related solution. We establish OS-ELM architecture and feed into the echo intensity sub-sequences as inputs, respectively, from the historical survey with manually labeled annotation, as well as the online observation that newly arrives at any time, in the context of few shot learning, to accelerate online detection and localization for sand waves. Model selection is involved by assessing performances across the configuration of hyper parameters, and the developed scheme will be comprehensively evaluated by the detection accuracy, False Negative Rate (FNR), False Positive Rate (FPR), and F-Score. We propose to utilize a collection of Gaussian-type morphological templates and construct one combination of matching criterion for similarity assessment in terms of cross-correlation and RMSE through the Taylor formula. We repeatedly deployed our AUV system around Jiaozhou Bay, Qingdao, China in 2019 and 2021, within the duration of more than one year, equipped with SSS. The basic parameters describing the morphological characteristics of sand waves will be elaborately discussed and finely grained in the study region, with the help of the statistics in the envelope demodulation, zero-crossing rate (ZCR) spectrum, and the cross-correlation coefficient; and the specific morphological parameters, including wave length, wave height, asymmetric index, will be further estimated, so as to explore the topology, morphology, and evolution in sand waves at a higher resolution with online computation strategies and provide insights into the spatial-temporal evidence in the formation and migration process.
The remainder of the paper is organized as follows: Section 2, Section 3 and Section 4 respectively describe the basic principle of the side scan sonar, the basics in ELM learning, and template matching. Section 5 introduces the sand wave online detection with OS-ELM and their morphological, geometrical, and topological characterization by echo intensity. Section 6 shows the simulation experiment and result analysis. Section 7 discusses the envelope demodulation, zero-crossing rate spectrum, cross-correlation coefficients, and the specific morphological parameters including wave length, wave height, and asymmetric index, by means of echo intensity in the SSS profile.

2. Side Scan Sonar Principle

The side scan sonar (SSS) is the one of the most basic acoustic sensors employed by AUVs during marine investigation tasks. When our AUV sails forward, the SSS transmits and receives echoes at certain time intervals, then displays as the echo intensity, and then composes the acoustic imaging for the seafloor sand waves. The transducers emit an acoustic pulse with directivity. When it reaches the seabed, the reflected waves or the backscattered waves return back along the original route, and the near echo first reach the transducers. In general, the hard, rough, and protuberant seabed could have strong echoes, while the soft, flat, and concave seabed will correspond to weak echoes. The concave sand waves at the bottom obscured by the crests could receive no echoes, yielding the shadow zone. Furthermore, the echo amplitude increases at the protuberant crest. Therefore, the echo amplitude reflects the ups and downs on seafloor sand waves, as is shown in Figure 1.
The spatial relationship between the operating range and AUV height is,
R h = R s 2 H a 2
where R h is the operating range per side, and R s is the slant range. H a is the AUV height above seabed.

3. Basic ELM

So far, ELM learning has been developed to work at a much faster learning speed with a higher generalization performance in the pattern recognition. For the given N learning samples { x i , y i } i = 1 N , where x i = [ x i 1 , x i 2 , , x i n ] T and y i = [ y i 1 , y i 2 , , y i m ] T , the standard model of the ELM learning can be written as the following matrix format:
H β = Y H ( x ) = [ h 1 h 2 h L ] = [ h 1 ( x 1 ) h L ( x 1 )               h i ( x j )         h 1 ( x N ) h L ( x N ) ] = [ φ ( w 1 × x 1 + b 1 ) φ ( w L × x 1 + b L )                         φ ( w i × x j + b i )               φ ( w 1 × x N + b 1 ) φ ( w L × x N + b L ) ] N × L β = [ β 1 , β 2 , , β L ] m × L T , Y = [ y 1 , y 2 , , y N ] m × N T
where w i = [ w i 1 , w i 2 , , w i r ] T , i = 1 , 2 , L is the weight vector connecting the i th hidden neuron and the input neurons, β i = [ β i 1 , β i 2 , , β i m ] T is the weight vector connecting the i th hidden neuron and the output neurons, w i × x j denotes the inner product of w i   and   x j , j = 1 , 2 , N , and there are L hidden neurons with the activation function φ ( x ) All kinds of activation functions can be chosen here, such as the Sigmoid function, the hard-limit function, the Gaussian function, and the multiquadric function.
If the activation function φ ( x ) , w and b are all set, the only learning parameter will be β . Different from the traditional learning algorithm, ELM tends to achieve the least training error and the least norm of output weight together. According to Bartlett׳s theory [52], when the feedforward neural networks obtain smaller training errors, the norms of weights are smaller, and the generalization performance of the networks is better, β = arg min ( | | H β Y | | 2 , | | β | | ) . In order to solve the formation, both the standard optimization method and the minimal norm least square method need to be adopted. The original implementation of ELM will then be β = H + Y , where H + denotes the Moore–Penrose generalized inverse of matrix H [37]. The orthogonal projection method can be used here when H T H is nonsingular and H + = ( H T H ) 1 H T , or when H H T is nonsingular and H + = H T ( H T H ) 1 . In addition, the resulting solution tends to be more stable with a better generalization performance by adding a positive value to the diagonal of H H T or H T H [53].

4. Template Matching

The template matching process involves shifting the template over the search area and computing the similarity between the template and the window in the search area over which the template lies. These two steps are described in detail below: (1) Identify the features underlying the search area and build the parameterized template waveform. (2) Match the search area with the parameterized template waveform by setting up an optimization matching criterion. The major similarity measures are the Root Mean Squared Error (RMSE) and the cross-correlation coefficient.
Let X = [ x 1 , x 2 , x m x M ] 1 × M be the search area, and T = [ t 1 , t 2 , t m , t M ] 1 × M be the template, and the similarity measurement is by RMSE J ( X , T ) to screen out the best matches.
J ( X , T ) = m = 1 M ( x m t m ) 2 / M
Cross-correlation coefficient r ( X , T ) is taken to measure the overall similarity between waveforms and template
r ( X , T ) = C o v ( X , T ) / D ( X ) D ( T ) C ov ( X , T ) = E [ ( X E ( X ) ) ( T E ( T ) ) ] = E [ X T ] E [ X ] E [ T ]
where C ov ( X , T ) is the covariance of the search area X and template T , D ( X ) and D ( T ) stand for the variance of search area X and template T respectively, E [ ] represents the mathematical expectation.

5. Sand Wave Detection and Morphological, Geometrical, Topological Characterization with Echo Intensity

We try to develop the sand wave online detection scheme from the echo intensity in high-resolution SSS mapping via OS-ELM and quantitatively assess the morphological, topological, and geometrical characterization. A brief flow chart of our proposed strategy is shown in Figure 2, which is made up of several correlative steps, including SSS imaging acquisition and preprocessing, sand wave online detection, morphological, geometrical, topological characterization, and the analysis of morphological parameters. Firstly, time varying gain (TVG), speed correction, and ensemble empirical mode decomposition (EEMD) denoising help to compensate the frequent motion in an AUV. We then explore online detection and localization for submarine sand waves through echo intensity, with the help of the OS-ELM learning framework. Furthermore, we denote a collection of Gaussian-type morphological templates, with one integrated matching criterion for similarity assessment. Finally, we statistically discuss the envelope demodulation, zero-crossing rate (ZCR) spectrum, and cross-correlation coefficient and estimate the specific morphological parameters, including wave length, wave height, and asymmetric index.

5.1. Echo Waveform Preprocessing

5.1.1. Echo Waveform Extraction

In practice, there are actually quite a few challenges in sand wave inspections via SSS for AUVs. Besides the common problem such as noise interfere and high computational demand, sand waves with a smaller wave length, especially megaripples, tend to be very difficult to distinguish. In addition to SSS imaging, the echo intensity directly reflects the nature of the seafloor with the corresponding geoacoustical properties, and the received backscattered signals are functions of both the geometry of the rough interface in sand waves and the material properties of the ocean bottom, including the density and sound speed ratios as well as the attenuation that is due to the dissipation of acoustic energy through the interface. If we start from the analysis of the sand wave profile instead of only the entire SSS imaging, it would be more beneficial to identify the ridge lines and valley lines for sand waves, infer the direction of sand wave movement, and calculate the sand wave migrate rate, with the advantages of low calculation cost and good real-time performance. Additionally, some smaller sand waves, such as megaripples, might be more easily finely detected after each ping in much detail. We hereby try to explore an alternative scheme for sand wave detection, by utilizing the side scan lines of echo intensity instead of only the mosaicked SSS imaging to accelerate the detection and localization of sand waves.

5.1.2. Time Varying Gain Correction

Acoustic waves experience spreading and absorption losses when traveling through sea water. The further the submarine sand waves are from the SSS transducer, the weaker the echoes we would receive. The horizontal distortion of echo intensity often appears, characterized with higher values in the center region and are smaller than the expected value in the edges. The relationship between echo attenuation and propagation distance can be constructed by propagation loss and reverberation level. One kind of time varying gain (TVG) correction is hereby taken to compensate acoustic attenuation over the distance for echo intensity in each ping as follows.
x i v T V G = T V G ( v ) / T V G max × x i v
where X i = [ x i 1 , x i 2 , x i v , x i V ] is the i th ping of the echo intensity curve, consisting of multiple echo intensity values, v [ 1 , w ] is defined as the v th position of the echo intensity value along one ping, with w the total width of the echo intensity ping, x i v is the v th echo intensity value in one ping, x i v T V G is the echo intensity value after TVG correction, and the function T V G ( v ) [54,55] is denoted as,
T V G ( v ) = A × log ( C v ) + B × ( C v ) C v = R s × | w / 2 v | / ( w / 2 )
where the correction coefficient A and the absorption coefficient B are related to the beam shape as well as the emission frequency of acoustic waves, and so forth, which is adjustable according to the scenarios to be applied, C v is the actual seabed projection distance between the current AUV waypoint and the v th echo intensity location, R s is the slant length of the sonar, and w is the width of echo intensity.

5.1.3. Speed Correction

The consecutive pings of echo intensity constitute SSS imaging, which are initially mapped one by one, as if the distance between adjoining pings were uniform. Due to the pose changes, or when the AUV does not keep to a constant speed and route, AUVs possibly experience longitudinal motion instability accompanied by severely fluctuating heave and pitch. The pings sometimes are not evenly and regularly spaced, leading to significant geometrical distortions in the along-track directions. If the AUV slows down, the sand waves in smaller size might look enwidened and enlarged. On the other hand, when speeding up the AUV, it might shrink or even conceal sand waves. Furthermore, if the AUV moves at too high speeds, it can leave unscanned regions in the seafloor that are reflected by echo intensity black gaps. The speed correction is hereby taken to align to the true proportion of sand wave features and fill in the gaps by interpolation. Let h i be the voyage distance of each ping data in the side scan sonar image, which could be denoted as [54]:
h i = w ( v ¯ i + v ¯ i + 1 ) t i 4 R s = ( v ¯ i + v ¯ i + 1 ) t i 2 × w 2 R s
where ( v ¯ i + v ¯ i + 1 ) t i / 2 refers to the actual distance of AUV in the along-track directions, and w / 2 R s represents the mapping coefficient between the width of the echo intensity and the actual across-track distance in the seabed, t i is the time difference between the i th ping and the ( i + 1 ) th ping that the SSS receives echoes, v ¯ i is the instantaneous sailing speed of AUV in the i th ping, v ¯ i + 1 is the instantaneous speed of AUV in the ( i + 1 ) th ping, R s is the slant length of the sonar, w is the width of echo intensity in the i th ping.

5.1.4. EEMD Denoising

Due to the existence of the environmental noise, reverberation noise, and instrument noise in the echo intensity of SSS imaging, recovering the signal with proper filtering is always essential. Frequently, classical filtering methods, such as the Wiener filtering, Gaussian filtering, and mean filtering, have been widely used, but these are still challenging when encountering sharp signals and impulses of short duration for the nonlinear and often nonstationary type.
We therefore employ the Ensemble Empirical Mode Decomposition (EEMD) [56,57] as a self-adaptive filtering for echo intensity denoising, which could adaptively decompose the nonlinear and nonstationary waveform into the sum of components, the Intrinsic Mode Functions (IMFs) and one residual component, and then distinguish and remove IMF with noise as its main component, eliminating the mode mixing problem by adding finite white noise to the investigated signal. The brief process of EEMD denoising [58,59] is expressed as follows:
(1)
Add a white noise with the given amplitude to the original signal in echo intensity;
(2)
Perform EMD [60] to the signal in echo intensity with the added white noise to obtain P Intrinsic Mode Function (IMF) components and one residual component, whereas the definition and acquisition process of IMF are relegated to Appendix A;
(3)
Repeat with the given number of trials. In every trial, the number of IMF P , is a constant. All the trials could be written as,
X i + n s = p = 1 P I M F s p + r s
where X i is the i th original signal in echo intensity; n s is the added white noise of the s th trial; I M F s p is the p th IMF component of the s th trial; r s is the residual component of the s th trial;
(4)
Calculate the ensemble mean of all trials and the final decomposition result of the original signal X i in echo intensity by EEMD could be expressed as,
I M F p = s = 1 S I M F s p / S r ( t ) = s = 1 S r s / S X i = p = 1 P I M F p + r
where S is the number of trials, and r is the residual component of EEMD;
(5)
To eliminate the noise and maintain the integrity of useful information, the first p 1 IMFs are selected to be processed with the EEMD denoising by a universal threshold formula as follows:
σ p = m e d i a n ( | I M F p | ) / 0.6745 T p = σ p 2 ln ( N )
where σ p is the standard deviation of the noise, N is the length of I M F p .
After determining the universal threshold T p , the soft threshold function is selected to filter noise component coefficients in the high frequency IMFs
I M F t h p = { sgn ( I M F p ) ( | I M F p T p | ) , | I M F p | T p 0 ,                                                                                     | I M F p | < T p ;
(6)
Finally, EEMD denoising is given by:
X i d n = p = 1 p 1 I M F t h p + p = p 1 + 1 P I M F p + r   .

5.2. Sand Wave Online Detection

In this paper, we make an attempt to explore an online auto-detection technique for submarine sand waves through echo intensity, with the help of the OS-ELM learning framework. OS-ELM inherits the advantage of ELM, which can provide good generalization performance at an extreme learning speed, moreover, OS-ELM [39] has an online sequential learning capability that does not require retraining when new observations arrive, and especially adapts to rapid responses to AUVs dynamically. The reason why we emphasize the online computation is that we attempt to explore the data-driven mode that could perform sand wave detection during the voyage, rather than from postprocessing after sailing back. The data-driven scheme could provide the AUV system with the capacities to execute complex real-time processing and decision making, with a resulting increase in its autonomy, flexibility, and adaptability, multiplying their data collection efficiency. An AUV could hereby modify the path planning during the mission and approach regions with sand waves, which is conducive to reducing the time complexity of our survey and improving the quality and efficiency of sand wave inspection tasks. We progressively collect SSS imaging, extract the echo intensity, and then perform preprocessing such as blind zone removal, TVG correction, speed correction, EEMD denoising, and so forth and further try to establish single hidden layer feedforward neural networks (SLFNs) to derive the relevant mathematical criterion for the optimization problem of sand wave detection. The flowchart of OS-ELM learning is shown in Figure 3, with the following steps.

5.2.1. Initialization Phase

(1)
Input a small chunk of echo intensity sub-sequence 0 = { ( X j ¯ , t j ¯ ) } j ¯ = 0 N 0 . Set the number of hidden neurons as N ˜ , the training sample X j ¯ = [ x j ¯ 1 , x j ¯ 2 , , x j ¯ d ] T , the number of X j ¯ in 0 as N 0 , the target output matrix as T 0 = [ t 1 , , t j ¯ , , t N 0 ] 1 × N 0 T and class label as t j ¯ = [ t j ¯ 1 , t j ¯ 2 , t j ¯ m ] T ;
(2)
Randomly assign the input weights a k and bias b k . Set the output of the k th hidden node with respect to the input X j ¯ as activation function G ( a k , b k , X j ¯ ) = 1 / ( 1 + e ( a k X j ¯ + b k ) ) ;
(3)
Calculate the initial hidden layer output matrix H 0 :
H 0 = [ G ( a 1 , b 1 , X 1 )                         G ( a 1 , b 1 , X N 0 ) G ( a N ˜ , b N ˜ , X 1 )                       G ( a N ˜ , b N ˜ , X N 0 ) ] N 0 × N ˜ ;
(4)
Estimate the initial output weight β ( 0 ) .
Calculate the intermediate variable P 0 = ( H 0 T H 0 ) 1 and initial output weight β ( 0 ) = P 0 H 0 T T 0 by the Moore–Penrose pseudo-inverse;
(5)
Set l = 0 .

5.2.2. Sequential Learning Phase

(6)
Set the number of training data in the ( l + 1 ) th group as N l + 1 = p = 0 l + 1 N p [ ( p = 0 l N p ) + 1 ] . Input the ( l + 1 ) th echo intensity sub-sequence chunk l + 1 = { ( X j ¯ , t j ¯ ) } j ¯ = ( p = 0 l N p ) + 1 p = 0 l + 1 N p . Set the target output matrix of those N l + 1 training data as T l + 1 = [ t ( p = 0 l N p ) + 1 , , t p = 0 l + 1 N p ] T ;
(7)
Compute the partial hidden layer output matrix H l + 1
H l + 1 = [ G ( a 1 , b 1 , X ( p = 0 l N p ) + 1 )                               G ( a 1 , b 1 , X p = 0 l + 1 N p ) G ( a N ˜ , b N ˜ , X ( p = 0 l N p ) + 1 )                               G ( a N ˜ , b N ˜ , X p = 0 l + 1 N p ) ] N l + 1 × N ˜ ;
(8)
Calculate the output weight β ( l + 1 ) ,
P l + 1 = P l P l P l + 1 T ( I + H l + 1 P l + 1 H l + 1 T ) 1 H l + 1 P l
β ( l + 1 ) = β ( l ) + P l + 1 H l + 1 T ( T l + 1 H l + 1 β ( l ) )
Set l = l + 1 .
Go to Step 2.

5.2.3. Test Phase

(9)
Input test samples l * = { X ( p = 0 l 1 N p * ) + 1 * , , X j ¯ * , , X p = 0 l N p * * } . Set the number of test samples l * as N l = p = 0 l N p * [ ( p = 0 l 1 N p * ) + 1 ] ;
(10)
Predict class labels O ( l + 1 )
O ( l + 1 ) = H l + 1 * β ( l + 1 ) = [ G ( a 1 , b 1 , X ( p = 0 l N p * ) + 1 * )                               G ( a 1 , b 1 , X p = 0 l + 1 N p * * )             G ( a N ˜ , b N ˜ , X ( p = 0 l N p * ) + 1 * )                                   G ( a N ˜ , b N ˜ , X p = 0 l + 1 N p * * ) ] N l + 1 * × N ˜ β ( l + 1 )

5.3. Morphological, Geometrical, Topological Characterization

5.3.1. Morphological Template Making

Although the waveform of sand waves reflected from SSS imaging could exhibit unique and variable and sometimes rather complicated shapes, the majority forms repeatable patterns which can be grouped into a number of classes represented by typical morphological patterns, which may strongly depend on the morphological parameters such as the wave length, the wave height, and asymmetry index of submarine sand waves. We hereby denote a collection of Gaussian-type templates T m to simulate the fluctuation characteristics of sand waves:
T m = b t + exp ( - ( x μ t ) 2 / 2 σ t 2 ) / ( σ t 2 π )
where μ t is the mathematical expectation of the Gaussian function, σ t represents the standard deviation of Gaussian function, and b t determines the starting point of the Gaussian function. The example Gaussian-type templates we have established for sand wave location and morphological extraction are shown in Figure 4, where b t is defined as ( 0 . 1 + 0 . 02 × s d ) , s d [ 1 , 20 ] controls the moving stride of Gaussian-type template T m on the y-axis, μ t ranges from 600 × u d , u d [ 0 , 15 ] controls the moving stride of Gaussian-type template T m on the x-axis, and σ t is chosen as ( 1 + 0.2 × w d ) , w d [ 0 , 10 ] controls the dispersion degree of Gaussian-type templates T m .

5.3.2. Morphological Matching Criterion

We further establish one integrated template matching criterion for the morphological similarity evaluation in terms of cross-correlation r ( X , T ) and RMSE J ( X , T ) , which aims at enhancing assessment in spatial similarity.
According to Taylor’s formula, ln ( 1 + x ) = x x 2 / 2 + x 3 / 3 + + ( 1 ) n x n / n and lim x 0 ln ( 1 + x ) / x = 1 , we constrain the similarity assessment by a logarithmic function, so as to construct the following template matching criterion for the echo intensity of sand waves:
G t = r ( χ j i , T m ) + ln ( 1 + 1 / J ( χ j i , T m ) )
where χ j i = [ x j 1 i , x j 2 i , x j m i x j M i ] 1 × M is the portion of the j th sub-sequences in i th ping echo intensity used to match the template, x j m i  represents the m th echo intensity value. Additionally, the T m = [ t 1 , t 2 , t m , t M ] 1 × M is the Gaussian template, which not only considers the template matching between waveforms at the identical scales, but also the similarity measure in the context of geometric morphology. The logarithmic function is deployed to reduce the computational cost when the RMSE value is too small. The basic morphological descriptive parameters for sand waves will be hereby automatically quantified and recorded through the echo intensity during the template matching processes, including the wave height, wave length, and asymmetry index of sand waves, which would be probably helpful in further discussion on the sand wave migration rate and direction, the hydrodynamic analysis, and the stability of underwater engineering and even drive our AUV to sail along the optimal route planning to collect the specified type of sand wave more specifically.

5.3.3. Zero-Crossing Rate

We further utilized the zero-crossing rate (ZCR) to examine the ratio of the sign changes of the waveform spectrum for sand waves, i.e., the number of times from the positive to negative or a reverse direction in a given period. Regardless of the scales and amplitudes submarine sand waves may exhibit, from small-scale megaripples to large-scale sand banks, they do have the periodic morphological characteristics reflected from the waveform of echo intensity. On the contrary, the echo waveform without sand waves tends to be distributed less orderly with random noises of high-frequency. The waveform in the sand wave region could be seen as the amplitude modulation process, with the random noise of high-frequency as the carrier and the sand wave morphology of relatively low-frequency as the modulating signal. The attribute of waveform ZCR is not related to the amplitude itself, with a strong anti-interference ability. As the waveform spectrum maintains more periodicity, ZCR values calculated for a certain period approach a constant value. Intuitively, the greater number of zero crossing of echo intensity in SSS imaging is, the more rapid changes and, accordingly, a high frequency spectrum it contains. Thus, ZCR could provide indirect clues about the frequency contents for sand waves. We make an attempt to identify the topology in submarine sand waves by counting the number of times they cross zero in the echo intensity and quantify the periodic variation. An extension of ZCR could be taken as a discriminative feature descriptor of the frequency content for low sampling frequency waveforms employing the existence of underlying information related to sand waves detection and morphological characteristics. The waveform with sand waves bear a low frequency range of ZCR which can be distinguished from noise and eliminate the noise and vibration of waveforms in regions without sand waves. ZCR could be denoted as:
Z = 1 2 v = 2 V | sgn [ x v i X ¯ i ] sgn [ x ( v - 1 ) i X ¯ i ] | / ( V 1 )
where x v i and x ( v 1 ) i are the v th and ( v 1 ) th echo intensity of the sonar signal X i , X ¯ i is the mean of X i and V is the length of X i . sgn [ ] is a symbolic function as follows:
sgn [ x v i ] = { 1 , x v i > 0 - 1 , x v i < 0 .

6. Simulation Experiment and Results Analysis

In our simulation experiment, we repeatedly deployed our AUV system around Jiaozhou Bay, Qingdao, China on 6 December 2019 and 13 January 2021 for sand wave inspection, respectively starting at 36°02′38″ N, 120°16′1″ E and 36°03′27″ N, 120°15′48″ E, and ending at 36°03′43″ N, 120°16′26″ E and 36°03′0″ N, 120°16′2″ E, by the high-resolution SSS of 150 kHz frequency, 200 m operating ranges per side, 1.25 cm across track resolution, 0.6° horizontal beam width, 40° vertical beam width, as is shown in Figure 5, with the first and second AUV trajectories respectively represented in blue and red on the left, and the manually labeled sand wave region on the right. The SSS imaging log with two AUV deployments is listed in Table 1.
The diameter, length, and weight of our AUV system is respectively 324 mm, 4 m and 280 kg, as is shown in Figure 6. A high-performance Industrial Personal Computer PC104, together with the NVIDIA Jetson TX2 platform established the distributed control architecture with the super computation unit, with the sensors and measurement instrument on board, including lcocean Shark-S450D side scan sonar, u-blox NEO-M8T GPS, Linkquest Navquest micro 600 doppler velocity log (DVL), Beidou era BD-FGI920 inertial navigation system (INS), Valeport miniIPS intelligent pressure sensor, Evologics S2C acoustic communication, and so forth. With the help of the positioning and navigation, we made SSS mosaicking during postprocessing, as shown in Figure 7, where on the left, the AUV path of the first time is in blue, with the route of the stitched image in red as an example, and the corresponding mosaicking result on the right.
At the beginning, we utilized the manually labeled sand wave annotation in historic surveys by research vessel with the high resolution SSS on board for the standard references. During the process of SSS imaging acquisition by AUV, we further accumulated and refined the sand wave annotation as our ground truth. Figure 8 shows the example SSS imaging collected and the corresponding labeling, with sand waves in yellow and non sand waves in blue, respectively.
We retrieved the echo intensity in each ping and then divided these into overlapping sub-sequences with the basic length of 600. Considering the operating range of 200 m per side, we made an initial estimation to cover the seabed 25 m for each sub-sequence. In total, 327,280 echo intensity sub-sequences were chosen for OS-ELM learning and statistical analysis, with 156,306 samples labeled in the category of sand waves, and 170,974 samples in the class without sand waves. For each class, we set aside one fifth for test and verification, respectively, and the rest for training. Figure 9 lists the example echo intensity sub-sequence extraction process from raw SSS imaging, with the centered red line indicating the extraction localization of echo intensity sub-sequence from raw SSS imaging patches.
Considering the typically uneven characteristics on across-track in SSS imaging, we made TVG correction against these effects, with the empirical correction index A = 5 and empirical correction bias B = 3.1 , and the corrected example SSS imaging of Figure 8a has been shown in Figure 10. In this context, we carried out a series of TVG correction experiments to analyze the selection of the relatively optimal parameters over our study area, with the combination of typical values A = 30 , B = 20 ; A = 15 , B = 10 ; A = 5 , B = 3.1 , A = 4 , B = 3 ; A = 2 , B = 1 , and so forth. Some of the TVG correction examples are shown in Figure 11, where both overcorrection and undercorrection occurred for some pairs of parameters, resulting in either secondary gray distortion or edge texture weakening. When we chose A = 5 , B = 3.1 , the radiometric distortion could be well eliminated with TVG correction, resulting in a relatively realistic reflection of the seafloor. Instead of only generating enhanced contrast processing, TVG correction provided the acoustic energy compensation with the function of echo attenuation and propagation distance.
We verified the validity of EEMD denoising to highlight the morphology of sand wave profile waveform with higher smoothing capability as well as edge preservation, in comparison to a few classical filters, such as the median filter and Gaussian filter, as is shown in Figure 12, with the common sliding window size l f = 3  and the standard deviation of Gaussian filter σ g f = 0.8 .
We further established OS-ELM architecture and fed into training samples for sand wave detection. Regarding the model selection, the average detection accuracy has been taken to assess across the configuration of hyper parameters. The number of hidden layer nodes has been considered to avoid the learning contingency and the corresponding test time has also been recorded, as is shown in Figure 13. It has been demonstrated from the experimental results that the sand wave detection rate initially increased when the relatively low number of hidden layer nodes rose until it reached 600, with the sand wave detection rate of 95.61% and the test time of 0.0018 s on average, and then the tendency of the detection rate remained or decreased after the certain value.
We also employed the basic MobileNet architecture [61] to perform sand wave detection with SSS imaging in the context of one light-weight deep learning framework. The sand wave detection results of example SSS imaging in Figure 11a have been compared in Figure 14, where the left is from the basic MobileNet by SSS imaging, and the right is from OS-ELM by echo intensity, with sand waves in yellow and non sand waves in cyan, respectively. It can be seen from our experiment results that MobileNet aimed to the sand wave detection in a large range, while our proposed scheme could concentrate more on the details of the sand wave profile distribution in every ping, especially showing superior performances for those sand waves in smaller sizes, e.g., megaripples.
We have also evaluated the sand wave detection performances of our proposed approach comprehensively in terms of the detection Accuracy, False Negative Rate (FNR), False Positive Rate (FPR), and F1-Score, and the evaluation metrics are defined as:
F P R = F P / ( F P + T N ) F N R = F N / ( T P + F N ) P r e s s i o n = T P / ( T P + F P ) R e c a l l = T P / ( T P + F N ) A c c u r a c y = ( T P + T N ) / ( T P + F P + T N + F N ) F 1 S c o r e = 2 × [ P r e c i s s i o n × R e c a l l / ( P r e c i s s i o n + R e c a l l ) ]
where True Positive (TP) refers to the set of echo intensity sub-sequence which falls into the sand wave class and is correctly classified. False negative (FN) refers to the set of echo intensity sub-sequence which falls into the sand wave class but is misclassified as the background. True negative (TN) refers to the set of echo intensity sub-sequence which falls into the background and is correctly classified. False Positive (FP) refers to the set of echo intensity sub-sequence which falls into the background class but is misclassified as the sand waves. The sand wave detection performances of AUV inspection in December 2019 and January 2021 are listed in Table 2, Table 3 and Table 4. It is revealed in our experiment that the sand wave detection rate with echo intensity could reach 95.61% averagely, and FNR 4.36%, FPR 4.82%, F1-Score 95.32%, with the average test time 0.0018 s and the total training time 5.7168 s for sea trials, comparable to MobileNet, but at a much higher speed and low computational cost, showing its superiority particularly for sand waves at smaller scales, such as megaripples.
Generally speaking, OS-ELM architecture for sand wave detection shows its adaptivity, but most of the difficulty still tends to come from the availability of qualified echo intensity waveform from SSS profile when deploying our AUV system. We also discussed the generalization capacity of our OS-ELM learning procedure for different spots and seasons with different states of AUV. We performed the OS-ELM learning framework in the training dataset at different spots, e.g., the study area around the South China Sea, and achieved the comparable average detection rate of 90.57%. We also carried out the simulation experiment in the study area of Jiaozhou Bay with the training samples respectively collected from different seasons, i.e., in December 2019, June 2020, and January 2021, and the sand wave detection rates reached 95.42%, 94.57%, and 95.67%. Although the migration and morphology of sand waves could be related to meteorological conditions, it has been demonstrated that by now the training procedure still behaved relatively stably, with a relatively reasonable impact on sand wave detection. In practice, we usually deploy our AUV system for a sand wave detection task under the level 3 sea state and collect SSS data in the relatively stable operation stage of AUV. We gather AUV motions of surge, sway, heave position, and angular rotations of roll, pitch, and yaw angle and also help with the compensation and correction of the position offset of the echo intensity in the along- and across-track directions to the navigation trajectory. Therefore, the learning procedure have been designed to perform relatively adaptively for some general situation, in sand wave detection applications.
We not only examined the existence of sand waves in the study area, but also focused on investigating what type of morphological properties the sand waves would probably possess. Therefore, we also established a variety of deformable slide Gaussian templates to match the sand wave profile waveform with the similarity metric. Some example Gaussian template matching process has been shown in Figure 15, with an average calculation speed of 0.53 s for each echo intensity sub-sequence in locating all the matching sand waves.
When the AUV system sails through the sand wave region, it is possibly not aligned perfectly with ridge lines. If there exists an angle between the direction vector of AUV with the front wave axis, on the basis of the edge detection in both horizontal and vertical directions, we could try to first estimate the horizontal tilt angle α and vertical tilt angle β by Radon transform, and then make use of shear transformation for the geometric correction of the entire SSS image before morphological template matching as follows,
[ x ˜ y ˜ ] = [ 1 sin ( α × π / 180 ) sin ( β × π / 180 ) 1 ] × [ x ˜ y ˜ ]
where x ˜ and y ˜ are the abscissa and ordinate of the image pixels in the original SSS imaging, x ˜ and y ˜ are the abscissa and ordinate of transformed image pixels in the newly generated SSS imaging. Figure 16 shows one rotation transformation example to SSS imaging with sand wave regions, before applying morphological Gaussian templates.

7. Discussion

From the perspective of echo intensity measurement in SSS imaging, sand waves appear highly irregular in size, shape, and spacing. In order to investigate the variability of geometrical properties in sand waves, it is of great importance to identify the locations and main directions of the crests and troughs for sand waves in their individual SSS profiles during the hydrodynamic and sediment dynamic processes.
We were hereby inspired to make the upper and lower envelopes for sand wave profile to compose the reconstructed boundaries within which the maxima and minima signals in echo intensity are contained in the following steps: (1) identify all the local extrema; (2) interpolate between extrema by a cubic spline. Since the geomorphology of sand waves would produce vibration signals of low frequency in echo intensity, the envelope demodulation process reflects the amplitude changes and direction of the sand wave signals with the extreme of the waveform, detect the vibration with low signal-to-noise ratio from massive noises by filtering, and further help determine sand waves types.
We further examined the zero-crossing rate spectrum for the example envelope demodulation, as is shown in Figure 17, where the crests and troughs could also be recognized by mapping between zero upcrossings and zero downcrossings in the bedform SSS profile with stability. We also made a statistical comparison of ZCR extracted from the echo intensity with or without sand waves in our study area when we deployed our AUV to Jiaozhou Bay for SSS imaging collection during the sea trial in January 2021, as is shown in Figure 18, and found that ZCR remained between 0.3–0.4 in the regions with sand waves and ranged from 0.4–0.6 for the regions without sand waves, with an average calculation speed of 0.56 ms for each echo intensity sub-sequence. It has been estimated from the experimental result that, compared to the regions with sand waves, those regions without sand waves might most probably bear the relatively large zero-crossing rate percentage and the denser zero-crossing spectrum in higher frequency, implying that the echo intensity with sand waves is of stronger regularity in a relatively low frequency, while the ZCR without sand waves tends to be distributed more disorderly and varies more frequently, so that we could further exclude the influence of noise interference in SSS imaging. To employ ZCR might possibly provide one of the distinguishing attributes in analyzing morphological properties for sand waves, which could also help reduce dimension, remove redundancy, improve the computational cost.
Furthermore, we established cross-correlation coefficients to explore a broad class of relationships among upper and lower envelope demodulation and evaluated the strength via the covariance matrix, either for a sample of echo intensity or a population within the whole research area. Figure 19 refers to the statistics of cross-correlation coefficients between the upper and lower envelopes in SSS profiles on both sand wave and non sand wave regions for our study area, with an average calculation speed of 0.12 s for each echo intensity sub-sequence. It has been demonstrated that for the local regions with the sand waves, the upper and lower envelopes co-moved with each other synchronously, and the cross-correlation coefficient was relatively large, implying that the amplitude rises and falls in sand waves could roughly reflect the fluctuation of submarine topography through the nearly periodic variation of echo intensity in SSS profiles. It has been calculated that when the cross-correlation coefficient was greater than 0.75, the possibility of the existence of sand waves approached 90% in terms of the detection precision. In the regions without sand waves, the cross-correlation coefficients remained relatively small, exhibiting the disaccordance and inconsistence in their distribution, which indicates the randomness and independence with more disorderly variation.
We further discussed how to refine the specific morphological parameters by means of echo intensity in SSS profile, including wave length, wave height, asymmetric index. The hydrodynamic instability generates coherent flow structures and contribute to the formation of sand waves. The wave length and wave height are the most intuitive reflections of sand wave morphology. Let the span of a wave be the wave length S L , which could be denoted as:
S L = S L R H × R H
where S L is the length of echo intensity sub-sequence in SSS imaging, R H represents the maximum number of echo intensity per side in SSS imaging, R H stands for the real distance that SSS projects acoustic waves onto the seafloor on each side, which could be derived by the spatial relationship between the operating range per side and AUV height above the seabed.
We calculated the statistics for the wave length estimation as above. Figure 20 shows the wave length calculation for the example SSS imaging in Figure 10a, and we inferred that the wave length around that region primarily ranged between 10 m and 25 m, responsible for 78% of all types of sand waves, which shows a significant consistency with the historic surveys executed by research vessels with a high resolution SSS on board. More precisely, in review of the whole region after sand wave detection, the sand waves at the wave length of 10–15 m in rose color accounted for 46% of all sand waves, the sand waves at the wave length of 15–25 m in red color almost reached to 32%, while the sand waves at the wave length less than 10 m in green color constituted 22%. However, the proposed scheme is more applicable and beneficial to those morphological parameter calculations for sand waves when the wave ridge lines in survey is parallel to AUV path planning.
We also analyzed the statistical morphological relationships between the wave length estimation in every single waveform and the corresponding local maxima echo intensity value among the local minima for sand waves, as is shown in Figure 21. It could be seen from the scatter plot that the sand waves in that region could be further roughly classified into two types according to their relationship distribution, i.e., the shorter wave lengths within 10 m with smaller maxima echo intensity less than 0.025, and the longer wave lengths between 10–30 m with the larger maxima echo intensity around 0.032, which could be characterized and approximated respectively by two linear functions with the slope 3.35, 1.14 and the intercept 0.64, 10.71, regarding 1000 times the echo intensity versus the wave length.
Although the wave height of sand waves is unable to be retrieved from SSS imaging accurately, it has been reported in previous studies [3,62,63] that the wave height of sand waves exhibits the tendency to increase with the increasing wave length, which indicated that the two morphological descriptive indices of sand waves potentially possessed significant positive correlation with each other under certain conditions.
Flemming [63] statistically investigated 1491 sand waves all over the world and proposed the an equilibrium geometrical relationship between a wave height S H and the wave length S L as follows,
S H = 0.0677 ( S L ) 0.8098 .
Flemming [63] also put forward an upper limit S H max for the wave height by:
S H max = 0.16 ( S L ) 0.84 .
In a different equilibrium geometrical relationship between the wave height S H and the wave length S L proposed by Franzetti et al. [62], about 500 sand waves in the continental shelf of western Brittany were identified and measured as follows,
S H = 0.0139 ( S L ) 1.164 .
Notwithstanding that a large number of sand waves in various formation environments have been considered, there are still examples not well described by the above Flemming’s empirical relationships; e.g., Van Landeghem et al. [64,65] have observed a very large dune with a length of 435 m and a height of about 35.5 m at a water depth of 90 m in the Irish Sea which had a height larger than Flemming’s global mean regression would predict. Generally, Flemming’s wave height prediction could be regarded as one of the currently most authoritative and commonly used formulas for the two morphological descriptive indices of sand waves globally. According to the wave length estimation results in our study region, we applied the empirical relationship in Flemming’s formula and deduced the possible wave height estimation around 0.77 m and the upper limit about 1.98 m for SSS imaging in the study region. To be specific, the wave height and upper limit within the wave length 10 m were estimated to be 0.32 m and 0.82 m, while for the wave length between 10–30 m, 0.86 m and 2.23 m, which showed good agreement with the previous research on the morphological characteristics of submarine sand waves around our study region in Jiaozhou Bay [66,67].
The internal structure of the sand waves apparently depends on the strength and degree of asymmetry of the governing currents, and the relatively asymmetrical currents would shape sand waves with one side too steep. The bigger the difference between the upstream and downstream slope projection in the sand wave is, the bigger the flow velocity difference will be reflected, which is much closer to the one-directional flow, and the activity of sand waves would rise. Conversely, the activity of sand waves becomes weak with better stability. We hereby tried to calculate the values of the upstream and downstream slope projection in echo intensity and explore the asymmetric index of sand waves, as is shown in Figure 22. The ratio between the upstream slope projection S u and the downstream slope projection S d be one asymmetric index R ,
R = S u / S d
The asymmetrical index estimation for sand waves in the above study region is shown in Figure 23. It has been demonstrated from our experimental results that those sand waves with the asymmetrical index R < 0.5 in red color were calculated to be more likely to incline westward, the sand waves with the asymmetrical index 0.5 < R < 1.5 in rose color look to be relatively stable and balanced between the upstream and downstream slope projection from the echo intensity with the relatively weak activity in movement, and the sand waves with the asymmetrical index R > 1.5 in green color were inferred to possibly lean eastward, respectively accounting for 72%, 22%, and 6% in the study region. As a whole, from the perspectives of morphology, the study region of sand waves would most probably be influenced by the westward tidal current compared to the eastward ones, implying the direction of hydrodynamic forces that represent the dominant ones in driving the sand wave movements, with a slightly higher strength compared to the other directions.
It is generally accepted that the lee-stoss asymmetry of sand waves is a sign for sand wave migration in the direction of the steeper slope [68]. We have therefore deployed our AUV twice around the same study region in December 2019 and January 2021 and collected SSS imaging within the duration of more than one year to verify the estimated direction of sand wave migration. Regarding the intersection of AUV trajectories of both sea trials in the yellow star, the comparison of echo intensity in SSS imaging is illustrated in Figure 24, with the left collected in December 2019 and the right collected in January 2021, and the echo intensity curves refer to the ping in the red line located on SSS imaging. It has been demonstrated from SSS profiles that compared to the echo intensity curves before and after one year, there is indeed exhibited a westward migration tendency of the sand waves, with the asymmetry index less than 0.5, and the steep slopes were assumed to shape the westward. It has been reported that the tidal current in Jiaozhou Bay could be characterized as a regular semidiurnal one, towards the westward direction at the flood tide and to the eastward direction at the ebb tide, where the current at the flood tide has been supposed to be generally stronger than the one at the ebb tide [69]. Li et al. stated that there is a similar type of hydrodynamic environment around Jiaozhou Bay and demonstrated the impact in the tidal flow and tidal flux with the cross-bay bridge built [69]. Bian et al. [67] studied the distribution, morphology, and migration of submarine sand waves around Jiaozhou Bay by means of the multibeam bathymetric data, SSS data, and seismic profiles and found that the sand waves, mainly exhibiting straight crest type and crescentic type, moved westward under the control of the tidal current. Our findings through morphological parameters with echo intensity in SSS imaging agree with the previous field surveys and proved to some extent the effectiveness of the proposed scheme.

8. Conclusions

Submarine sand waves exhibit ubiquitous geomorphology at a variety of scales, including sand banks, sand waves, and megaripples, highly reflecting the coupling variation of multivariate marine environmental factors, such as the tidal characteristics, sediment grain sizes, related residual currents, and so forth. From the morphological and topological perspectives, submarine sand waves declare a uniform wavelength at equilibrium and develop defects while adjusting to variation. In this paper, we quantitatively explore the formation and migration cues in sand waves on topology, morphology, and evolution from high-resolution mapping via a side scan sonar (SSS) mounted in an Autonomous Underwater Vehicle (AUV), in favor of online sequential Extreme Learning Machine (OS-ELM). We first make an attempt to utilize echo intensities directly derived from SSS to help accelerate online detection and localization for sand waves, with the help of time varying gain (TVG) correction, speed correction, ensemble empirical mode decomposition (EEMD) denoising, and so on to compensate for the frequent motion variation on the AUV in synthesizing consecutive uniform pings of SSS imaging over time. We then denote a collection of Gaussian-type morphological templates to capture the fluctuation characteristics of sand waves and establish one integrated matching criterion for similarity assessment in terms of cross-correlation and RMSE. We statistically investigate the envelope demodulation, zero-crossing rate (ZCR) spectrum, and cross-correlation coefficients to examine the high variability of sand waves regarding their size, shape, and spacing and further potentially estimate the specific morphological parameters, including wave length, wave height, and asymmetric index, by means of echo intensity in the SSS profile. In our simulation experiment, we repeatedly deployed our AUV system around Jiaozhou Bay, Qingdao, China in 2019 and 2021, within the duration of more than one year, equipped with SSS of 150 kHz frequency, 200 m operating ranges per side, 1.25 cm across track resolution, 0.6° horizontal beam width, 40° vertical beam width, as well as a high-performance PC and GPU for control and online computation, DVL, INS, digital compass, and depth gauge for navigation. We established OS-ELM architecture, retrieved 327280 overlapping echo intensity sub-sequences, fed into 196368 training samples, and comprehensively evaluated the detection accuracy, False Negative Rate (FNR), False Positive Rate (FPR), and F1-Score with the sand wave annotation in historic surveys by research vessel and the increased manually labeling for the standard references. The model selection was elaborated by assessing the detection performances across the configuration of hyper parameters to determine the number of hidden layer nodes. It is revealed in our experiment that the sand wave detection rate with echo intensity could reach 95.61% averagely, and FNR 4.36%, FPR 4.82%, F1-Score 95.32%, with the average test time of 0.0018 s for sea trials, which is comparable to the image-based detection performances in the light-weight deep learning framework, e.g., MobileNet, but at a much higher speed and low computational cost, showing its superiority particularly for sand waves at smaller scales, such as megaripples. The statistics of ZCR and the cross-correlation coefficient between envelope demodulation illustrated differences in sand wave and non sand wave regions, implying the regularity, periodicity, and synchronousness in the low-frequency waveform of sand waves. We inferred that the wave length around that region primarily ranged between 10 m and 25 m, responsible for 78% of all types of sand waves, characterized and approximated respectively by linear functions with the slope 3.35, 1.14 and the intercept 0.64, 10.71, and the sand waves in the study would most probably exhibit the influences of the westward currents compared to the eastward ones, implying the possible dominant direction of hydrodynamics in driving the sand wave migration, which shows a strong consistency with previous studies. The proposed scheme permits to delicately explore the submarine geomorphology of sand wave variation, allowing the morphological imagery at a higher resolution and online computation strategies, greatly breaking through the restrictions of field survey with more flexibility and adaptivity, and symmetrically integrating the geometrical, topological, and morphological evidences yielded from echo intensity in SSS imaging, which would help provide opportunities for identifying the underlying dominantly correlative factors in nature and potentially suggest the leading roles in the responses of sand waves during their formation and migration.

Author Contributions

Conceptualization, R.N. and L.Z.; methodology, R.N., L.Z., F.Y. and S.R.; validation, R.N., L.Z., X.G., X.L. and B.H.; formal analysis, R.N., L.Z., X.G., X.L. and B.H.; investigation, F.Y., X.L., B.H., L.Z. and S.R.; resources, B.H. and F.Y.; writing—original draft preparation, R.N., L.Z. and X.G.; writing—review and editing, R.N. and L.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key R&D Program (2019YFC1408304), the National Program of International S&T Cooperation (2015DFG32180), the National High-Tech R&D 863 Program (2014AA093410), the Natural Science Foundation of China (31202036), the National Science & Technology Pillar Program (2012BAD28B05), the National Key R&D Program (2016YFC0301400) and the Natural Science Foundation of China (41376140).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available from the author upon reasonable request. The data are not publicly available right now due to project privacy.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

  • Each Intrinsic Mode Function (IMF) of the ith ping echo intensity X i is defined as a function that must satisfy the following conditions,
    (1)
    In the whole data set, the number of extrema and the number of zero-crossings must be equal or differ at most by one;
    (2)
    At any point, the mean value of the envelope defined by the local maxima and the envelope defined by the local minima is zero.
  • The EMD decomposition of the i th ping echo intensity  X i is according to the following steps:
    (1)
    Initialization: r s 0 = X i , p = 1 ;
    (2)
    Computation of the p th IMF, I M F s p :
    (a)
    Initialization: h 0 = r s ( p 1 ) , j ˜ = 1
    (b)
    Find the local maxima and minima of h j ˜ 1
    (c)
    Interpolate the local minima (resp. maxima) to get e 1 (resp e 2 ) by cubic spline
    (d)
    Compute the mean of these envelops: E = ( e 1 + e 2 ) / 2
    (e)
    Obtain the j ˜ th component h j ˜ by h j ˜ = h j ˜ 1 - E
    (f)
    If the above two conditions of IMF are satisfied, then I M F s p = h j ˜ ; else, j ˜ = j ˜ + 1 ;
    (3)
    Calculate residual component r s p = r s ( p 1 ) - I M F s p ;
    (4)
    If r s p is not monotonic, go to step (2); otherwise, the decomposition is complete.
When the decomposition is complete, we can write the i th ping echo intensity  X i as X i = p = 1 P I M F s p + r s P , P .

References

  1. Bellec, V.K.; Be, R.; Bjarnadóttir, L.R.; Albretsen, J.; Elvenes, S. Sandbanks, sandwaves and megaripples on Spitsbergenbanken, Barents Sea. Mar. Geol. 2019, 416, 105998. [Google Scholar] [CrossRef]
  2. Wang, Z.; Liang, B.; Wu, G.; Borsje, B.W. Modeling the formation and migration of sand waves: The role of tidal forcing, sediment size and bed slope effects. Cont. Shelf Res. 2019, 190, 103986. [Google Scholar] [CrossRef]
  3. Blondeaux, P.; Vittori, G. A model to predict the migration of sand waves in shallow tidal seas. Cont. Shelf Res. 2016, 112, 31–45. [Google Scholar] [CrossRef]
  4. Zang, Z.; Liang, C.; Gao, F. Application of ROMS for Simulating Evolution and Migration of Tidal Sand Waves. In Proceedings of the Sixth International Conference on Asian and Pacific Coasts, Hong Kong, China, 14–16 December 2011. [Google Scholar]
  5. Knaapen, M.A.F.; Henegouw, C.N.V.B.; Hu, Y.Y. Quantifying bedform migration using multi-beam sonar. Geo Mar. Lett. 2005, 25, 306–314. [Google Scholar] [CrossRef] [Green Version]
  6. Giri, S.; Shimizu, Y. Numerical computation of sand dune migration with free surface flow. Water Resour. Res. 2006, 42. [Google Scholar] [CrossRef] [Green Version]
  7. Ye, Y. Submarine Sand Waves and Sand Ridges. Mar. Geo-Hazards China 2017, 523–554. [Google Scholar] [CrossRef]
  8. Fredsøe, J. On the development of dunes in erodible channels. J. Fluid Mech. 1974, 64, 1–16. [Google Scholar] [CrossRef]
  9. King, E.; Conley, D.; Masselink, G.; Leonardi, N.; McCarroll, R.; Scott, T. The impact of waves and tides on residual sand transport on a sediment-poor, energetic, and macrotidal continental shelf. J. Geophys. Res. Ocean. 2019, 124, 4974–5002. [Google Scholar] [CrossRef]
  10. Nienhuis, J.H.; Perron, J.T.; Kao, J.C.; Myrow, P.M. Wavelength selection and symmetry breaking in orbital wave ripples. J. Geophys. Res. Earth Surf. 2014, 119, 2239–2257. [Google Scholar] [CrossRef] [Green Version]
  11. Cappietti, L.; Chopard, B. A Lattice Boltzmann Study of the 2d Boundary Layer Created by AN Oscillating Plate. Int. J. Mod. Phys. C 2006, 17, 39–52. [Google Scholar] [CrossRef]
  12. Bagnold, R.A. Motion of waves in shallow water. Interaction between waves and sand bottoms. Proc. R. Soc. Lond. Series A Math. Phys. Sci. 1946, 187, 1–18. [Google Scholar]
  13. Besio, G.; Blondeaux, P.; Frisina, P. A note on tidally generated sand waves. J. Fluid Mech. 2003, 485, 171–190. [Google Scholar] [CrossRef] [Green Version]
  14. Santen, R.B.V.; Swart, H.E.D.; Dijk, T.A.G.P.V. Sensitivity of tidal sand wavelength to environmental parameters: A combined data analysis and modelling approach. Cont. Shelf Res. 2011, 31, 966–978. [Google Scholar] [CrossRef]
  15. Borsje, B.W.; Roos, P.C.; Kranenburg, W.M.; Hulscher, S.J.M.H. Modeling tidal sand wave formation in a numerical shallow water model: The role of turbulence formulation. Cont. Shelf Res. 2013, 60, 17–27. [Google Scholar] [CrossRef]
  16. Venditti, J.G.; Church, M.; Bennett, S.J. Morphodynamics of small-scale superimposed sand waves over migrating dune bed forms. Water Resour. Res. 2005, 41, W10423–W10437. [Google Scholar] [CrossRef] [Green Version]
  17. Grasmeijer, B.; Kleinhans, M. Observed and predicted bed forms and their effect on suspended sand concentrations. Coastal Eng. 2004, 51, 351–371. [Google Scholar] [CrossRef]
  18. Li, Y.; Ma, L.-M.; Yang, J.-S.; Shi, A. Study on stability of sand waves by satellite sensing. In Proceedings of the First Asian and Pacific Coastal Engineering Conference, APACE2001, Dalian, China, 18–21 October 2001; pp. 850–856. [Google Scholar]
  19. Németh, A.A.; Hulscher, S.J.; de Vriend, H.J. Modelling sand wave migration in shallow shelf seas. Cont. Shelf Res. 2002, 22, 2795–2806. [Google Scholar] [CrossRef]
  20. Borsje, B.W.; Kranenburg, W.M.; Roos, P.C.; Matthieu, J.; Hulscher, S.J.M.H. The role of suspended load transport in the occurrence of tidal sand waves. J. Geophys. Res. Earth Surf. 2014, 119, 701–716. [Google Scholar] [CrossRef] [Green Version]
  21. Houthuys, R.; Trentesaux, A.; Wolf, P.D. Storm influences on a tidal sandbank’s surface (Middelkerke Bank, southern North Sea). Mar. Geol. 1994, 121, 23–41. [Google Scholar] [CrossRef]
  22. Minor Salvatierra, M.; Aliotta, S.; Susana Ginsberg, S. Morphology and dynamics of large subtidal dunes in Bahia Blanca estuary, Argentina. Geomorphology 2015, 246, 168–177. [Google Scholar] [CrossRef]
  23. Wei, S.; Wang, M.; Li, X.; Pichel, W.G. Ocean sand ridge signatures in the Bohai Sea observed by satellite ocean color and synthetic aperture radar measurements. Remote Sens. Environ. 2011, 115, 1926–1934. [Google Scholar]
  24. Durán, R.; Guillén, J. Continental Shelf Landforms. In Submarine Geomorphology; Springer: Cham, Switzerland, 2018; pp. 185–206. [Google Scholar]
  25. Guerrero Sánchez, Q. Bedforms and Associated Sediment Dynamics on the Inner Shelves at Different Spatio-Temporal Scales. Ph.D. Thesis, Universitat Politècnica de Catalunya, Barcelona, Spain, 2019. [Google Scholar]
  26. Williams, D.P.; Baralli, F.; Micheli, M.; Vasoli, S. Adaptive underwater sonar surveys in the presence of strong currents. In Proceedings of the 2016 IEEE International Conference on Robotics and Automation (ICRA), Stockholm, Sweden, 16–21 May 2016; pp. 2604–2611. [Google Scholar]
  27. Hagen, O.K.; Anonsen, K.B. Terrain referenced navigation using sand dunes and sand ripples in presence of tides. In Proceedings of the OCEANS 2018 MTS/IEEE Charleston, Charleston, SC, USA, 22–25 October 2018; pp. 1–6. [Google Scholar]
  28. Kennish, M.J.; Haag, S.M.; Sakowicz, G.P.; Tidd, R.A. Side-Scan Sonar Imaging of Subtidal Benthic Habitats in the Mullica River–Great Bay Estuarine System. J. Coast. Res. 2004, 45, 227–240. [Google Scholar] [CrossRef]
  29. Mhammdi, M.N.N.; Geawhari, M.A.G.M.A.; Trakadas, A.T.A.; Huff, L.H.L. AUV for Shallow Water Hydrographic Surveys: Detecting Bedforms in Oued Loukkos Estuary (Morocco); European Association of Geoscientists & Engineers: Houten, The Netherlands, 2013. [Google Scholar]
  30. Baba, T.; Matsumoto, H.; Kashiwase, K.; Hyakudome, T.; Kaneda, Y.; Sano, M. Micro-bathymetric Evidence for the Effect of Submarine Mass Movement on Tsunami Generation During the 2009 Suruga Bay Earthquake, Japan. In Submarine Mass Movements and Their Consequences. Advances in Natural and Technological Hazards Research; Springer: Dordrecht, The Netherlands, 2012. [Google Scholar]
  31. Wang, X.; Jiao, J.; Yin, J.; Zhao, W.; Han, X.; Sun, B. Underwater sonar image classification using adaptive weights convolutional neural network. Appl. Acoust. 2019, 146, 145–154. [Google Scholar] [CrossRef]
  32. Einsidler, D.; Dhanak, M.; Beaujean, P.P. A Deep Learning Approach to Target Recognition in Side-Scan Sonar Imagery. In Proceedings of the OCEANS 2018 MTS/IEEE Charleston, Charleston, SC, USA, 22–25 October 2018; pp. 1–4. [Google Scholar]
  33. Wu, M.; Wang, Q.; Rigall, E.; Li, K.; Zhu, W.; He, B.; Yan, T. ECNet: Efficient convolutional networks for side scan sonar image segmentation. Sensors 2019, 19, 2009. [Google Scholar] [CrossRef] [Green Version]
  34. Sung, F.; Yang, Y.; Zhang, L.; Xiang, T.; Torr, P.H.S.; Hospedales, T.M. Learning to Compare: Relation Network for Few-Shot Learning. In Proceedings of the 2018 IEEE/CVF Conference on Computer Vision and Pattern Recognition, Salt Lake City, UT, USA, 18–23 June 2018. [Google Scholar]
  35. Najafabadi, M.M.; Vilanustre, F.; Khoshgoftaar, T.T.; Seliya, N.; Wald, R.; Muharemagic, E. Deep learning applications and challenges in big data analytics. J. Big Data 2015, 2, 1. [Google Scholar] [CrossRef] [Green Version]
  36. Chen, X.-W.; Lin, X. Big data deep learning: Challenges and perspectives. IEEE Access 2014, 2, 514–525. [Google Scholar] [CrossRef]
  37. Huang, G.B.; Zhu, Q.Y.; Siew, C.K. Extreme learning machine: Theory and applications. Neurocomputing 2006, 70, 489–501. [Google Scholar] [CrossRef]
  38. Huang, G.B.; Zhou, H.; Ding, X.; Rui, Z. Extreme Learning Machine for Regression and Multiclass Classification. IEEE Trans. Syst. Man Cybern. Part B 2012, 42, 513–529. [Google Scholar] [CrossRef] [Green Version]
  39. Liang, N.Y.; Huang, G.B.; Saratchandran, P.; Sundararajan, N. A Fast and Accurate Online Sequential Learning Algorithm for Feedforward Networks. IEEE Trans. Neural Netw. 2006, 17, 1411–1423. [Google Scholar] [CrossRef]
  40. Rong, H.J.; Huang, G.B.; Sundararajan, N.; Saratchandran, P. Online Sequential Fuzzy Extreme Learning Machine for Function Approximation and Classification Problems. IEEE Trans. Syst. Man Cybern. Part B 2009, 39, 1067–1072. [Google Scholar] [CrossRef]
  41. Sun, Y.; Yuan, Y.; Wang, G. An OS-ELM based distributed ensemble classification framework in P2P networks. Neurocomputing 2011, 74, 2438–2443. [Google Scholar] [CrossRef]
  42. Twichell, D.C. Bedform distribution and inferred sand transport on Georges Bank, United States Atlantic Continental Shelf. Sedimentology 1983, 30, 695–710. [Google Scholar] [CrossRef]
  43. Collins, M.B.; Shimwell, S.J.; Gao, S.; Powell, H.; Taylor, J.A. Water and sediment movement in the vicinity of linear sandbanks: The Norfolk Banks, southern North Sea. Mar. Geol. 1995, 123, 125–142. [Google Scholar] [CrossRef]
  44. Reynaud, J.v.; Tessier, B.; Proust, J.o.; Dalrymple, R.; Lericolais, G. Eustatic and hydrodynamic controls on the architecture of a deep shelf sand bank (Celtic Sea). Sedimentology 2010, 46, 703–721. [Google Scholar] [CrossRef]
  45. Hennings, I.; Lurin, B.; Vernemmen, C.; Vanhessche, U. On the behaviour of tidal current directions due to the presence of submarine sand waves. Mar. Geol. 2000, 169, 57–68. [Google Scholar] [CrossRef]
  46. van Dijk, T.A.G.P.; Lindenbergh, R.C.; Egberts, P.J.P. Separating bathymetric data representing multiscale rhythmic bed forms: A geostatistical and spectral method compared. J. Geophys. Res. 2008. [Google Scholar] [CrossRef]
  47. Venditti, J.G.; Church, M.; Bennett, S.J. On the transition between 2D and 3D dunes. Sedimentology 2005, 52, 1343–1359. [Google Scholar] [CrossRef]
  48. Wynn, R.B.; Masson, D.G.; Bett, B.J. Hydrodynamic significance of variable ripple morphology across deep-water barchan dunes in the Faroe–Shetland Channel. Mar. Geol. 2002, 192, 309–319. [Google Scholar] [CrossRef]
  49. Dijk, T.A.G.P.V.; Kleinhans, M.G. Processes controlling the dynamics of compound sand waves in the North Sea, Netherlands. J. Geophys. Res. Earth Surf. 2005, 110. [Google Scholar] [CrossRef]
  50. Passchier, S.; Kleinhans, M.G. Observations of sand waves, megaripples, and hummocks in the Dutch coastal area and their relation to currents and combined flow conditions. J. Geophys. Res. Earth Surf. 2005. [Google Scholar] [CrossRef]
  51. Crawford, A.; Skarke, A. Automatic detection of sand ripple features in sidescan sonar imagery. In Proceedings of the 2014 Oceans, St. John’s, NL, Canada, 14–19 September 2014; pp. 1–5. [Google Scholar]
  52. Bartlett, P.L. The sample complexity of pattern classification with neural networks. IEEE Trans. Inf. Theory 1998, 44, 525–536. [Google Scholar] [CrossRef] [Green Version]
  53. He, B.; Xu, D.; Nian, R.; van Heeswijk, M.; Yu, Q.; Miche, Y.; Lendasse, A. Fast Face Recognition Via Sparse Coding and Extreme Learning Machine. Cogn. Comput. 2014, 6, 264–277. [Google Scholar] [CrossRef] [Green Version]
  54. Han, C.H.; Zhang, J.M.; Liang, J.F.; Yin, R.G. Design and implementation of side-scan sonar data management system. Mar. Sci. Bull. 2011, 30, 188–193. [Google Scholar]
  55. Jianhu, Z.; Jun, Y.; Hongmei, Z.; Junxia, M. A New Radiometric Correction Method for Side-Scan Sonar Images in Consideration of Seabed Sediment Variation. Remote Sens. 2017, 9, 575. [Google Scholar]
  56. Jia, Y.; Li, G.; Dong, X.; He, K. A novel denoising method for vibration signal of hob spindle based on EEMD and grey theory. Measurement 2021, 169, 108490. [Google Scholar] [CrossRef]
  57. Wu, Z.; Huang, N.E. Ensemble empirical mode decomposition: A noise-assisted data analysis method. Adv. Adapt. Data Anal. 2009, 1, 1–41. [Google Scholar] [CrossRef]
  58. Mao, Q.; Fang, X.; Hu, Y.; Li, G. Chiller sensor fault detection based on empirical mode decomposition threshold denoising and principal component analysis. Appl. Therm. Eng. 2018, 144, 21–30. [Google Scholar] [CrossRef]
  59. Jin, T.; Li, Q.; Mohamed, M.A. A novel adaptive EEMD method for switchgear partial discharge signal denoising. IEEE Access 2019, 7, 58139–58147. [Google Scholar] [CrossRef]
  60. Huang, N.E. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. Lond. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  61. Howard, A.G.; Zhu, M.; Chen, B.; Kalenichenko, D.; Wang, W.; Weyand, T.; Andreetto, M.; Adam, H. MobileNets: Efficient Convolutional Neural Networks for Mobile Vision Applications. arXiv 2017, arXiv:1704.04861. [Google Scholar]
  62. Franzetti, M.; Roy, P.L.; Delacourt, C.; Garlan, T.; CancouëT, R.; Sukhovich, A.; Deschamps, A. Giant dune morphologies and dynamics in a deep continental shelf environment: Example of the banc du four (Western Brittany, France). Marine Geol. 2013, 346, 17–30. [Google Scholar] [CrossRef] [Green Version]
  63. Flemming, B.W. Zur klassifikation subaquatischer, stromungstransversaler Transportkorper. Boch.geol.u.geothchn.arb 1988, 29, 44–47. [Google Scholar]
  64. Landeghem, K.J.J.V.; Wheeler, A.J.; Mitchell, N.C.; Sutton, G. Variations in sediment wave dimensions across the tidally dominated Irish Sea, NW Europe. Mar. Geol. 2009, 263, 108–119. [Google Scholar] [CrossRef]
  65. Landeghem, K.V.; Baas, J.H.; Mitchell, N.C.; Wilcockson, D.; Wheeler, A.J. Reversed sediment wave migration in the Irish Sea, NW Europe: A reappraisal of the validity of geometry-based predictive modelling and assumptions. Mar. Geol. 2012, 295–298, 95–112. [Google Scholar] [CrossRef]
  66. Bian, S.-H.; Xia, D.-X.; Chen, Y.; Zhao, Y.-X. The classification characteristics and developing factors of the sand waves at the mouth of Jiaozhou Bay. Period. Ocean Univ. China 2006, 36, 327–330. [Google Scholar]
  67. Zhao, Y.; Liu, B.; Li, X.; Wu, J.; Su, T. Topography Feature and Migration of Submarine Sand Waves in Jiaozhoy Bay Mouth. Oceanol. Limnol. Sin. 2006, 37, 464. [Google Scholar]
  68. Knaapen, M.A.F. Sandwave migration predictor based on shape information. J. Geophys. Res. Earth Surf. 2005, 110. [Google Scholar] [CrossRef] [Green Version]
  69. Li, P.; Li, G.; Qiao, L.; Chen, X.; Shi, J.; Gao, F.; Wang, N.; Yue, S. Modeling the tidal dynamic changes induced by the bridge in Jiaozhou Bay, Qingdao, China. Cont. Shelf Res. 2014, 84, 43–53. [Google Scholar] [CrossRef]
Figure 1. SSS characterization in sand wave inspections.
Figure 1. SSS characterization in sand wave inspections.
Sensors 21 03283 g001
Figure 2. The flow chart of the proposed scheme.
Figure 2. The flow chart of the proposed scheme.
Sensors 21 03283 g002
Figure 3. The flowchart of OS-ELM learning.
Figure 3. The flowchart of OS-ELM learning.
Sensors 21 03283 g003
Figure 4. Example Gaussian-type templates.
Figure 4. Example Gaussian-type templates.
Sensors 21 03283 g004
Figure 5. AUV trajectories and the manually labeled sand wave region. (a) AUV trajectories. (b) The manually labeled sand wave region.
Figure 5. AUV trajectories and the manually labeled sand wave region. (a) AUV trajectories. (b) The manually labeled sand wave region.
Sensors 21 03283 g005
Figure 6. AUV deployment.
Figure 6. AUV deployment.
Sensors 21 03283 g006
Figure 7. AUV path and SSS mosaicking. (a) AUV path. (b) SSS mosaicking.
Figure 7. AUV path and SSS mosaicking. (a) AUV path. (b) SSS mosaicking.
Sensors 21 03283 g007
Figure 8. Example SSS annotation. (a) Example SSS imaging. (b) Labeling.
Figure 8. Example SSS annotation. (a) Example SSS imaging. (b) Labeling.
Sensors 21 03283 g008
Figure 9. Example echo intensity sub-sequence extraction from raw SSS imaging.
Figure 9. Example echo intensity sub-sequence extraction from raw SSS imaging.
Sensors 21 03283 g009
Figure 10. TVG correction for example SSS imaging. (a) The corrected example SSS imaging 1. (b) The corrected example SSS imaging 2.
Figure 10. TVG correction for example SSS imaging. (a) The corrected example SSS imaging 1. (b) The corrected example SSS imaging 2.
Sensors 21 03283 g010
Figure 11. TVG correction performance comparison. (a) Raw data. (b) A = 30 , B = 20 . (c) A = 15 , B = 10 . (d) A = 5 , B = 3.1 . (e) A = 4 , B = 3 . (f) A = 2 , B = 1 .
Figure 11. TVG correction performance comparison. (a) Raw data. (b) A = 30 , B = 20 . (c) A = 15 , B = 10 . (d) A = 5 , B = 3.1 . (e) A = 4 , B = 3 . (f) A = 2 , B = 1 .
Sensors 21 03283 g011
Figure 12. Example sand wave profile waveform and denoising. (a) Raw echo intensity sub-sequence. (b) Median filtering. (c) Gaussian filtering. (d) EEMD denoising.
Figure 12. Example sand wave profile waveform and denoising. (a) Raw echo intensity sub-sequence. (b) Median filtering. (c) Gaussian filtering. (d) EEMD denoising.
Sensors 21 03283 g012
Figure 13. Sand wave detection rate and time regarding to the number of hidden layer nodes in OS-ELM. (a) Test recognition rate. (b) Test time.
Figure 13. Sand wave detection rate and time regarding to the number of hidden layer nodes in OS-ELM. (a) Test recognition rate. (b) Test time.
Sensors 21 03283 g013
Figure 14. Sand wave detection for example SSS imaging. (a) Basic MobileNet with SSS imaging. (b) OS-ELM with echo intensity.
Figure 14. Sand wave detection for example SSS imaging. (a) Basic MobileNet with SSS imaging. (b) OS-ELM with echo intensity.
Sensors 21 03283 g014
Figure 15. Example sand wave profile waveform matching.
Figure 15. Example sand wave profile waveform matching.
Sensors 21 03283 g015aSensors 21 03283 g015b
Figure 16. Geometric correction example. (a) Raw SSS imaging. (b) Raw echo intensity values. (c) SSS imaging after geometric correction. (d) Echo intensity values after geometric correction.
Figure 16. Geometric correction example. (a) Raw SSS imaging. (b) Raw echo intensity values. (c) SSS imaging after geometric correction. (d) Echo intensity values after geometric correction.
Sensors 21 03283 g016
Figure 17. Zero-crossing spectrum. (a) Without sand wave. (bd) With sand waves.
Figure 17. Zero-crossing spectrum. (a) Without sand wave. (bd) With sand waves.
Sensors 21 03283 g017
Figure 18. ZCR statistics. (a) With sand waves. (b) Without sand waves.
Figure 18. ZCR statistics. (a) With sand waves. (b) Without sand waves.
Sensors 21 03283 g018
Figure 19. Cross-correlation statistics. (a) With sand wave. (b)Without sand wave.
Figure 19. Cross-correlation statistics. (a) With sand wave. (b)Without sand wave.
Sensors 21 03283 g019
Figure 20. Wave length estimation distribution for sand waves.
Figure 20. Wave length estimation distribution for sand waves.
Sensors 21 03283 g020
Figure 21. Morphological relationship of wave length estimation versus local maxima echo intensity for sand waves.
Figure 21. Morphological relationship of wave length estimation versus local maxima echo intensity for sand waves.
Sensors 21 03283 g021
Figure 22. Asymmetric index definition.
Figure 22. Asymmetric index definition.
Sensors 21 03283 g022
Figure 23. Asymmetrical index estimation distribution for sand waves.
Figure 23. Asymmetrical index estimation distribution for sand waves.
Sensors 21 03283 g023
Figure 24. Example echo intensity in SSS imaging at the intersection of two AUV trajectories within the duration of more than one year. (a) SSS imaging in December 2019. (b) SSS imaging in January 2021. (c) Echo intensity curve regarding the ping in red on SSS imaging (a). (d) Echo intensity curve regarding the ping in red on SSS imaging (b).
Figure 24. Example echo intensity in SSS imaging at the intersection of two AUV trajectories within the duration of more than one year. (a) SSS imaging in December 2019. (b) SSS imaging in January 2021. (c) Echo intensity curve regarding the ping in red on SSS imaging (a). (d) Echo intensity curve regarding the ping in red on SSS imaging (b).
Sensors 21 03283 g024
Table 1. SSS imaging with AUV deployment in Jiaozhou Bay.
Table 1. SSS imaging with AUV deployment in Jiaozhou Bay.
YearNameStart and End TimeDurationPing Number
2019Jiaozhou Bay 110:31:55–12:44:002:12:0521868
2021Jiaozhou Bay 214:00:23–16:32:522:32:2925181
Table 2. Sand wave detection performance on Jiaozhou Bay 1.
Table 2. Sand wave detection performance on Jiaozhou Bay 1.
Jiaozhou Bay 1FNR
(%)
FPR (%)Accuracy (%)F1-Score
(%)
Training Time (s)Testing Time (s)
MobileNet V34.936.3794.6994.079013.12470.0201
Ours4.415.0795.4295.292.47650.0013
Table 3. Sand wave detection performance on Jiaozhou Bay 2.
Table 3. Sand wave detection performance on Jiaozhou Bay 2.
Jiaozhou Bay 2FNR (%)FPR (%)Accuracy (%)F1-Score
(%)
Training Time (s)Testing Time (s)
MobileNet V34.655.4995.3295.2710802.67450.0247
Ours4.284.7395.6795.383.21790.0021
Table 4. Sand wave detection performance on Jiaozhou Bay 1&2.
Table 4. Sand wave detection performance on Jiaozhou Bay 1&2.
Jiaozhou Bay 1 & 2FNR
(%)
FPR (%)Accuracy (%)F1-Score (%)Training Time (s)Testing Time (s)
MobileNet V34.736.2594.8794.4718704.46570.0239
Ours4.364.8295.6195.325.71680.0018
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nian, R.; Zang, L.; Geng, X.; Yu, F.; Ren, S.; He, B.; Li, X. Towards Characterizing and Developing Formation and Migration Cues in Seafloor Sand Waves on Topology, Morphology, Evolution from High-Resolution Mapping via Side-Scan Sonar in Autonomous Underwater Vehicles. Sensors 2021, 21, 3283. https://doi.org/10.3390/s21093283

AMA Style

Nian R, Zang L, Geng X, Yu F, Ren S, He B, Li X. Towards Characterizing and Developing Formation and Migration Cues in Seafloor Sand Waves on Topology, Morphology, Evolution from High-Resolution Mapping via Side-Scan Sonar in Autonomous Underwater Vehicles. Sensors. 2021; 21(9):3283. https://doi.org/10.3390/s21093283

Chicago/Turabian Style

Nian, Rui, Lina Zang, Xue Geng, Fei Yu, Shidong Ren, Bo He, and Xishuang Li. 2021. "Towards Characterizing and Developing Formation and Migration Cues in Seafloor Sand Waves on Topology, Morphology, Evolution from High-Resolution Mapping via Side-Scan Sonar in Autonomous Underwater Vehicles" Sensors 21, no. 9: 3283. https://doi.org/10.3390/s21093283

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop