China’s First Demonstration of Cobalt-rich Manganese Crust Thickness Measurement in the Western Pacific with a Parametric Acoustic Probe

Cobalt-rich manganese crusts (CRCs) are important as a potential mineral source that could occur throughout the Pacific on seamounts, ridges, and plateaus. We built a prototype parametric acoustic probe to complete the task of in-situ thickness measurements to estimate the volumetric distribution of deep-sea mineral. The prototype is designed with dual-channels for receiving the primary and secondary signal, which lays a foundation for improving the thickness extraction algorithm. Considering that the signal quality is degraded by the system interference and ambient noise, some improvements to the algorithm are proposed by including the wavelet-based envelope extraction method and the adaptive estimation strategy based on the dual-channel information. Additionally, wavelet regression is applied to reduce the measuring noise assuming that the CRCs have local thickness invariability. The algorithm is suitable for the CRCs with the structure of the multilayers at the top surface and one single layer at the bottom surface. A laboratory experiment is performed to validate the effectiveness of the algorithm. The experiments carried out on the China Ocean 51th voyage in the Western Pacific Ocean on Aug 30, 2018, are described and the data obtained by using the sit-on-bottom stationary measurement are processed to validate the design of the prototype.


Introduction
Cobalt-rich manganese crusts (CRCs) are important as a mineral resource on the seafloor for copper, cobalt, nickel, platinum, manganese, thallium, tellurium, and other metals [1]. CRCs could occur throughout the Pacific on seamounts, ridges, and plateaus, where currents have kept the rocks, swept clean of sediments, at least intermittently, for millions of years [1]. Crusts form at water depths of about 400 to 4000 m, with the thickest and most Co-rich crusts occurring at depths of about 800 to 2500 m, which may vary on a regional scale. Gravity processes, sediment cover, submerged and emergent reefs, and currents control the distribution and thickness of crusts on seamounts [1][2][3][4].
Cobalt has many special properties and plays an important role in industry. For example, its properties allow for increasing chip performance, lead it to have a major place in the semiconductor manufacturing industry. As cobalt works better at higher densities, it may be called to also play a decisive role in the new and innovative artificial intelligence field. It is estimated that the resources of cobalt in the CRCs may be comparable with the land-based resources. Rising demand for minerals and metals, in tandem with the depletion of land-based resources, has led to a surge of interest in marine It is very common to measure the thickness of the target by using the reflected mechanical waves. Such methods have been widely used in multiple applications, such as geological exploration [13,14] and medical acoustics [15][16][17]. In this study, we emphasize on the work for improving the performance of the thickness measurements, considering that PPPAAP17 is a parametric array acoustic probe designed for the application of estimating the resources of CRCs. In practice, the necessary parameters, i.e., the sound velocity, the arrival time of the top surface, and the arrival time of the bottom surface, are calculated by processing the received signal that contains complicated context with the top-bottom peak detection. The reliability of the calculation is affected by noise in the recorded signal arising due to scattering, multi-path reflection, and seafloor features such as local inclusions inside the crust [5]. Hence, considering that the signal quality is degraded by the system interference and ambient noise, the dual-channel signal information and the inner connection between the large amounts of the recorded data are exploited. Some improvements are introduced to the processing algorithm, including the adaptive candidate selection strategy and the wavelet techniques [18,19]. For clarification, the key ideas are listed here: • The envelope is extracted by using wavelet transformation based on the dual-channel information which gives better stability compared to the traditional Hilbert envelope extraction method.

•
The determination of the arrival time of the top surface not only relies on an adaptive candidate selection strategy suitable for real-time processing logic but the auxiliary information of secondary channel, while that of the bottom surface relies on the adaptive candidate selection strategy, the geometric constraints, and the auxiliary information of the primary channel.

•
The wavelet regression is applied to further reduce ambient noise. The reason behind this is the assumption that the CRCs have the local thickness invariability. The utilization of neighborhood information of the large amounts of recorded data can further improve the robustness.
In this paper, we pay attention to multiple aspects, i.e., system design, processing algorithm, the laboratory experiment, and the real experiment. For the remaining sections, the problem formation is presented firstly. Afterward, the system description which includes the hardware configuration and the array design is briefly presented. Next, the processing algorithm is elucidated. Finally, the processing and analysis of the experiments are demonstrated, and the comparisons are given to verify the effectiveness of the system and the algorithm.

System Description
The parametric acoustic probe is designed to perform accurate remote measurements of CRCs with thicknesses greater than 30 mm and less than 350 mm. To satisfy both the requirements of along-track resolution and the measuring range, the proper trade-off should be made. We have known that finer along-track resolution needs the shorter wavelength of the signal, whereas the stronger penetrating capability needs the longer wavelength of the signal considering the attenuation. The probe transmits a 1 MHz, high-frequency amplitude modulated signal to generate a narrow 100 kHz beam that penetrates the target, which could achieve a narrow beam and increase the range of operation of the system. An acoustic beam of frequency 100 kHz on its target has been found to provide a suitable balance between resolution and signal attenuation when measuring CRCs with thicknesses in the range of interest [20].
As shown in Figure 1, PPPAAP17 mainly consists of two receiving transducers, two filtering and compensation boards, the transmitting transducer array, the transmit control board, the data collection and the controlling board, and the host computer. The real photo of the electronic processing unit is presented in Figure 2a, with the transmit control board, the filtering and compensation board, and the data collection and controlling board inside. When the system works, the host computer sends the transmitting parameter and command to the data collection and controlling board and the gain compensation parameters to the filtering and compensation board with the Ethernet protocol, respectively. Afterward, the data collection and controlling board will forward the information to the transmit control board to periodically drive the transmitting transducer array to fire. In the meantime, the sampling, storage, and controlling board is synchronized to receive the reflected signal from two separate receiving transducers, with one fed into the primary channel and the other fed into the secondary channel. Considering that the probe, PPPAAP17, is designed for ROV currently, the size of the array should satisfy the requirements of both mass and volume. Note that the array has an impact on the across-track resolution since the interference caused by the scattering waves depends on the beamwidth. In other words, the cross-track resolution is finer if the beamwidth is narrower. We have performed many acoustic field simulations for multiple arrays. As shown in Figure 2b, the transmitting array design is optimal at the aspect of the penetrating capability and the across-track resolution.
At the aspect of receiving, the primary signal of high-frequency and the secondary parametric signal of low-frequency is obtained by feeding the signal received by the transducer centered at the array to two designed filters, which is marked in red. The two designed filters are used to improve the Signal to Noise Ratio (SNR) of the signal. The gain compensation part can be adjusted by the data collection and controlling board to obtain suitable signal amplitude. The filter and sampling parameters are listed in Table 1. Note that the secondary channel has two different modes with different Band Pass Filter (BPF) and sampling parameters. The sampled data are transferred and stored in a real-time way and processed by using the on-line controlling and processing software to obtain thickness information with an effective algorithm.  Considering that the probe, PPPAAP17, is designed for ROV currently, the size of the array should satisfy the requirements of both mass and volume. Note that the array has an impact on the across-track resolution since the interference caused by the scattering waves depends on the beamwidth. In other words, the cross-track resolution is finer if the beamwidth is narrower. We have performed many acoustic field simulations for multiple arrays. As shown in Figure 2b, the transmitting array design is optimal at the aspect of the penetrating capability and the across-track resolution.
At the aspect of receiving, the primary signal of high-frequency and the secondary parametric signal of low-frequency is obtained by feeding the signal received by the transducer centered at the array to two designed filters, which is marked in red. The two designed filters are used to improve the Signal to Noise Ratio (SNR) of the signal. The gain compensation part can be adjusted by the data collection and controlling board to obtain suitable signal amplitude. The filter and sampling parameters are listed in Table 1. Note that the secondary channel has two different modes with different Band Pass Filter (BPF) and sampling parameters. The sampled data are transferred and stored in a real-time way and processed by using the on-line controlling and processing software to obtain thickness information with an effective algorithm.

Processing Algorithm
The aim of the algorithm is to automatically calculate the thickness of the CRCs with dualchannel information, assuming that the CRCs are with the structure of the multilayers at the top surface and only one single layer at the bottom surface. For this algorithm, it is demonstrated to handle the case of two layers at the top surface, while the case of three can be deduced in the same way. The CRCs with more than three layers at the top surface are rare. The processing method, as illustrated in Figure 3, is developed to extract the thickness of the CRCs. The algorithm mainly contains five steps based on the integrated utilization of the dual-channel information, i.e., primary signal Step 2. Determine the arrival time of the top surface with the auxiliary information of the secondary channel. The central idea of this step is that the top arrival time is calculated by using the primary signal with the help of the secondary signal. Such a method has some advantages. Firstly, compared with the secondary signal, the reflected primary signal of the top surface has a relatively higher SNR since the received power of the primary signal is much larger than the secondary signal. Moreover, the resolution of the primary signal is also finer than the secondary signal. The reason that we could use the secondary signal for reference is that the reflection will happen simultaneously on dual-channel. Another important consideration is based on the modeling of the structure of the CRCs. Multiple reflections will occur within the top surface, thus leading to multiple peaks appearing in the primary envelope. The energy of the peaks varies even in the same position which is impacted by the environment, which contributes to inaccuracy in the arrival time estimation. We give a simple way to handle this problem with robust performance.
Specifically, as marked yellow in Figure 3, this step contains several sub-steps. Firstly, the algorithm checks whether the SNR of the primary signal SNRp satisfies the reliable condition to perform further processing:

Processing Algorithm
The aim of the algorithm is to automatically calculate the thickness of the CRCs with dual-channel information, assuming that the CRCs are with the structure of the multilayers at the top surface and only one single layer at the bottom surface. For this algorithm, it is demonstrated to handle the case of two layers at the top surface, while the case of three can be deduced in the same way. The CRCs with more than three layers at the top surface are rare. The processing method, as illustrated in Figure 3, is developed to extract the thickness of the CRCs. The algorithm mainly contains five steps based on the integrated utilization of the dual-channel information, i.e., primary signal S p (t) and the secondary signal S d (t). The main steps are listed as follows: Step 1. Extract the envelope of S p (t) and that of S d (t) with the wavelet method to obtain W S p (t) and W(S d (t)), respectively.
Step 2. Determine the arrival time of the top surface with the auxiliary information of the secondary channel. The central idea of this step is that the top arrival time is calculated by using the primary signal with the help of the secondary signal. Such a method has some advantages. Firstly, compared with the secondary signal, the reflected primary signal of the top surface has a relatively higher SNR since the received power of the primary signal is much larger than the secondary signal. Moreover, the resolution of the primary signal is also finer than the secondary signal. The reason that we could use the secondary signal for reference is that the reflection will happen simultaneously on dual-channel. Another important consideration is based on the modeling of the structure of the CRCs. Multiple reflections will occur within the top surface, thus leading to multiple peaks appearing in the primary envelope. The energy of the peaks varies even in the same position which is impacted by the environment, which contributes to inaccuracy in the arrival time estimation. We give a simple way to handle this problem with robust performance.
Specifically, as marked yellow in Figure 3, this step contains several sub-steps. Firstly, the algorithm checks whether the SNR of the primary signal SNR p satisfies the reliable condition to perform further processing: where max[·] denotes calculating the peak value of the signal, E[·] denotes the average operator, and the primary signal threshold, ksnr p , is a constant larger than one. Secondly, the primary adaptive threshold can be set to the maximum value according to: where k p is a constant less than one Next, the maxima vector search for the top surface, which contains all the candidates of arrival time value, needs to be done. The window is used to partition the signal into some parts, whose size can be empirically set as 200 µs. By using the thresholds (2), the maxima vector is masked by removing the points that are lower than the threshold to reduce the number of candidates. If the masked vector with the length of n, N p = [N 1 , N 2 , . . . , N n ] is not empty, we can select the candidates based on some criterion. One typical criterion is using the descending selection, i.e., choosing the current maximum from the masked maxima vector first, and the new maxima vector can be given as The next step is to sequentially fetch two adjacent candidates from the candidate vector: They can be called candidate pair if they satisfy the envelope and the arriving time constraints: where TH N denotes the resolution parameter and can be set according to the resolution. Regardless of the result of the candidate pair or one candidate, the next step is to judge whether the selected point Ndes i p ∈ Ndes p satisfies the dual-channel constraints (DC) by using: where N j d denotes the arbitrary element of the filtered maxima vector N d for secondary channel, f sL and f sH denotes the sampling frequency of the primary channel and secondary channel, respectively, and TH ∆t denotes the threshold for arrival time secondary for dual-channel. In other words, Equation (7) exploits the principle of almost simultaneous reflection of the primary signal and secondary signal at the top surface. The division by the sampling frequency converts the discrete number to the arrival time. If the selected point satisfies the constraints, the arrival time of the top surface t top will be determined. Otherwise, the algorithm removes it and continues to go through all the elements of the new maxima vector.
If the estimation satisfies the dual-channel condition, we can determine the arrival time by different cases as: The reason behind Equation (8) is that we need to solve the ambiguity of the reflection of the multiple layers that occurred at the surface of the CRCs. If we know that two reflections with similar energy happen at the top surface, we will take the arrival time of the latest one as the estimation. The reason behind Equation (8) is that we need to solve the ambiguity of the reflection of the multiple layers that occurred at the surface of the CRCs. If we know that two reflections with similar energy happen at the top surface, we will take the arrival time of the latest one as the estimation. If only one reflection exists or two reflections with the very different energy, we will only take the estimation of the only candidate's arrival time of the top surface.
Step 3. Determine the arrival time of the bottom surface with the auxiliary information of the secondary channel. Firstly, similar to the processing of the primary channel, the algorithm checks whether the SNR of the secondary signal SNR d satisfies the reliable condition to perform further processing with the expression: where ksnr d is a constant larger than one. Secondly, the adaptive threshold can be set to the maximum value according to: where k d is a constant less than one. Next, the following steps, i.e., maxima vector search for bottom surface, maxima vector filtering, and selection of the arrival time candidate, are performed to determine the arrival time of the bottom surface. They are very similar to the processing of the primary channel. Note that the maxima vector filtering for secondary signal also utilizes the arrival time of the top surface as a reference. After such steps, we can get the bottom arrival candidates Ndes d with the criterion of descending.
It is worth noting that in such steps we exploit the geometric condition to reduce the interferences. The geometric condition is based on the fact that the CRCs of the minimum thickness of 30 mm, d min , can be of most interest as a potential mineral resource, and those of the maximum thickness of 350 mm, d max , can be practically found. In other words, the arrival time candidate of the bottom surface: where c represents the average sound velocity within the sampling CRC and the factor of two means the two-way traveling time. Finally, we can get the estimation of the arrival time of the bottom surface, t bottom .
Step 4. Obtain the thickness of CRCs Y i by using the following equation: By combining every scalar d i measured at time i, we can get the thickness vector Y by: where M denotes the total record number, F w denotes the working frequency, and T denotes the measuring duration.
Step 5. Remove the spikes that exist in the thickness vector by using the wavelet regression method. In this study, the details of the key issues will be discussed in the following sections.

Velocity Measurement
Since the sound velocity of CRCs is not a constant, a feasible way is to measure the values of several samples that are growing in a specific area and take the averaged values as the reference velocity. Here, we adopt the insert substitution method to measure the arrival time difference to calculate the sound velocity. The average sound velocity can be calculated by, where c 0 denotes the underwater sound velocity, L denotes the thickness of the test sample, and ∆t denotes the time difference, i.e., the arrival time without inserting the sample minus that after inserting the sample. The sound velocity measurement can be implemented for multiple times and the average velocity is taken as the reference value.

Dual-Channel Information
Compared with other algorithms that only use one-channel information for processing, the proposed algorithm utilizes dual-channel information which can provide additional information to extract the thickness information more precisely and robust. When the system works, the secondary signal often contains the interferences, making it difficult to estimate the reflection time. Compared with the secondary signal, the primary signal with a centroid frequency of about 1 MHz has a better resolution and less interference when it encounters the top surface. The reflection time of the primary signal not only can be used for the final calculation of the thickness but provides important prior information to determine the reflection on the bottom surface. The secondary signal has a great penetrating capability, which is used for determining the reflection at the bottom surface.

Wavelet-Based Envelope Extraction
The envelope extraction is a key step to obtain the top and bottom reflection time. Hilbert transformation is widely used in many applications related to location or distance measurement. However, extracting envelope with Hilbert transformation is limited by several factors. Firstly, it could be only applied to the narrow-band signal. For the thickness extraction application, the filter for secondary signal has configurations of different bandwidth. Secondly, it has the disadvantage of extracting the high-frequency component as a part of the envelope, which further causes spines. Such a phenomenon will largely impact the estimation of the top or bottom arrival time.
Hence, we use the complex continuous wavelet transform for analysis. The complex wavelet families contain Complex Gaussian Wavelets, Complex Morlet Wavelets, Complex Frequency B-Spline Wavelets, and Complex Shannon Wavelets. To obtain the envelope of the signal f (t), we will adopt the Complex Shannon Wavelets with scale and shift parameter, i.e., where f b is the bandwidth parameter, f c is a wavelet center frequency, a is the scale parameter, and b is the shift parameter. Generally, the value of ( f b , f c ) varies from multiple combinations such as (1, 1.5), (1, 1), (1, 0.5) and so on. In the thickness application, we empirically choose (1, 1.5) as the parameters. Considering the computational time, a can be set to be 32 typically. We will apply the complex continuous wavelet transform to the signal, and we will get the following equation: where W r f (a, b) and W i f (a, b) denote the real and imaginary parts, respectively. Therefore, the envelope W f can be explicitly represented as:

Wavelet Regression
When we estimate the thickness after obtaining the top surface and bottom surface arrival time, another important step is to remove the spikes that exist in the thickness signal. General wavelet thresholding is not capable of handling such signals which may have a small number of discontinuities. The method of wavelet regression which penalizes lack of smoothness in the reconstruction is very simple but can maintain sharp changes in the function and suppressing noise spikes.
Assume that the thickness vector which is represented as, Y i = THK i + n i , i = 1, 2, . . . , N, where THK i denotes the extracted thickness, n i denotes the measuring noise, and N is the total number of measurements. After wavelet regression, we can getTHK i = a i . In vector notation, we can get a = W T TWY, where T is the diagonal matrix with diagonal elements given by: and W is the discrete wavelet transform matrix. The purpose of the regression is to establish the zeros of the diagonal of T. Then we can build and minimize the proper function to determine the form of T [21]. For example, we can penalize the residual sum of squares and the second derivative of the estimator. After wavelet regression, we can get a good thickness estimation that exhibits good properties in terms of both fit and smoothness, while the sharp changes still retain well.

Laboratory Experiment with the Phantom
To systematically evaluate the effectiveness of the method, a laboratory experiment for measuring a phantom of known thickness was implemented. The experiment configuration is presented in Figure 4. The phantom of the POM material with the thickness of 100.00 mm is placed 1.30 m distance away from the transmitting transducer array.
In this experiment, the filter parameter is set to be mode 2. The primary data and secondary data contain 160 records, which means the total measuring time is 3.20 s. Taking the 100th signal of the record for example, the primary and secondary channels, as shown in Figure 5a,b, are of good quality with the maximum amplitude value of 2.53 V and 2.93 V, respectively. Note that the comparisons of the wavelet method and Hilbert method for the dual-channel are given in Figure 5c,d, respectively. Considering that the primary channel has a relatively narrow bandwidth, the difference is not obvious in the waveform although the former has a better quantitative performance than the latter. Figure 5d shows that the envelope extraction of the secondary signal is not robust by the Hilbert method since there are many small high-frequency spikes, while that of the wavelet method is very stable.
(a) Note that the comparisons of the wavelet method and Hilbert method for the dual-channel are given in Figure 5c,d, respectively. Considering that the primary channel has a relatively narrow bandwidth, the difference is not obvious in the waveform although the former has a better quantitative performance than the latter. Figure 5d shows that the envelope extraction of the secondary signal is not robust by the Hilbert method since there are many small high-frequency spikes, while that of the wavelet method is very stable. Note that the comparisons of the wavelet method and Hilbert method for the dual-channel are given in Figure 5c,d, respectively. Considering that the primary channel has a relatively narrow bandwidth, the difference is not obvious in the waveform although the former has a better quantitative performance than the latter. Figure 5d shows that the envelope extraction of the secondary signal is not robust by the Hilbert method since there are many small high-frequency spikes, while that of the wavelet method is very stable. The processing parameters are listed in Table 2. The SNR of the primary signal is 50.46 dB, and it is much higher than ksnr p . After performing the adaptive threshold masking and the priority calculation with the criterion of amplitude descending, the new maxima vector of the length of five can be obtained as: Ndes p = [8740, 17453, 6824, 8837, 15546]. Taking the first candidate for example, it is processed as one candidate since it does not satisfy the condition of the candidate pair. The next step is to judge whether the selected point Ndes i p ∈ Ndes p satisfies the dual-channel constraints by Equation (7). By comparing the real arrival time difference of 9 µs with the threshold TH ∆t , we can verify that the first candidate satisfies the constraints and the value of t top is obtained as 1748.0000 ms. Table 2. Parameters for the laboratory measurement with the phantom.

Parameters Symbols Value
The primary/secondary signal threshold ksnr p /ksnr d 6 dB/6 dB The primary/secondary adaptive threshold k p /k d 0.02/0.10 The threshold for the arrival time secondary TH ∆t 20 µs The threshold for the candidate pair TH N 150 The SNR of the secondary signal is 44.13 dB, and it is much higher than ksnr n . After performing the adaptive threshold masking, maxima vector filtering with the reference of the top surface, and the priority calculation with the criterion of amplitude descending, the new maxima vector can be obtained as: Ndes d = [5868, 3682, 3680, 3513, 3521, 6998, 5840, 7174, 6035, 7201, 5821, 6657].
Taking the first candidate for example, it is removed since the calculated thickness of 1280.80 mm exceeds the constraints by using Equation (7). By looping the algorithm, we can obtain the second candidate satisfies the constraints and the value of the bottom t bottom is obtained as 1841.0000 ms. The processing parameters are listed in  (7). By comparing the real arrival time difference of 9 μs with the threshold t TH Δ , we can verify that the first candidate satisfies the constraints and the value of top t is obtained as 1748.0000 ms.
Taking the first candidate for example, it is removed since the calculated thickness of 1280.80 mm exceeds the constraints by using Equation (7). By looping the algorithm, we can obtain the second candidate satisfies the constraints and the value of the bottom bottom t is obtained as 1841.0000 ms.

(a)
The number of the measurements The results of the arrival time of the top surface and bottom surface, and the estimated thickness are presented in Figure 6a-c, respectively. The quantitative comparison of the results with the different method is given in Table 3. The sound velocity of the phantom is measured as 2160.0 m/s in the laboratory. It can be seen that the estimated average thickness by wavelet method is closer to the ground-truth thickness of 100.0 mm and the standard deviation is smaller. Next, we use the Haar wavelet of level 6 for wavelet regression. The regression estimate result is 100.55 mm for the Wavelet method, while that is 95.18 mm for the Hilbert method. Since the laboratory experimental environment is good, the improvement of the regressed result is not very obvious. The final obtained results based on the 160 repeated measurements verify the effectiveness of the algorithm.  Figure 6a-c, respectively. The quantitative comparison of the results with the different method is given in Table 3. The sound velocity of the phantom is measured as 2160.0 m/s in the laboratory. It can be seen that the estimated average thickness by wavelet method is closer to the ground-truth thickness of 100.0 mm and the standard deviation is smaller. Next, we use the Haar wavelet of level 6 for wavelet regression. The regression estimate result is 100.55 mm for the Wavelet method, while that is 95.18 mm for the Hilbert method. Since the laboratory experimental environment is good, the improvement of the regressed result is not very obvious. The final obtained results based on the 160 repeated measurements verify the effectiveness of the algorithm.

The Real Experiment in the Western Pacific Ocean
The experiments were carried out on the China Ocean 51th voyage in the Western Pacific Ocean on Aug 30, 2018. The sea-trials were performed in the Western Pacific Ocean with the sit-on-bottom stationary measurements aiming to investigate the performance of the acoustic probe and the algorithm. As shown in Table 4, during the survey several measurements were taken at station 10 in the CRCs mountain regions of the interest. As shown in Table 4, we take the data of S1 for example to demonstrate the data process and analysis. A brief explanation of the measuring way, i.e., the sit-on-bottom stationary measurement, is given here. When performing the sit-on-bottom stationary measurement, HAIMA ROV lands fixed at the shoulder of the CRC mountain. As shown in Figure 7a, the electronic processing unit of PPPAAP17 marked with a red dashed box lies at the bottom of the HAIMA Remotely Operated Vehicle (ROV). The transmitter transducer array and receiving transducer array of PPPAAP17 marked with a blue dashed box is mounted on the hydraulic system of HAIMA ROV, as shown in Figure 7b. It can be seen that the vertical distance from the array to the ground could be adjusted. For the stationary measurement, the transmitter transducer array and receiving transducers mounted on HAIMA ROV are about 300 mm away from the seabed in the vertical direction.

The Real Experiment in the Western Pacific Ocean
The experiments were carried out on the China Ocean 51th voyage in the Western Pacific Ocean on Aug 30, 2018. The sea-trials were performed in the Western Pacific Ocean with the sit-on-bottom stationary measurements aiming to investigate the performance of the acoustic probe and the algorithm. As shown in Table 4, during the survey several measurements were taken at station 10 in the CRCs mountain regions of the interest. As shown in Table 4, we take the data of S1 for example to demonstrate the data process and analysis. Table 4. The data description of the experiment in the Western Pacific.

Station Name
No. Start Time (UTC+8) END TIME 10 Sit-on-bottom 1(S1) 2018-8-30-13-24-22 2018-8-30-13-41-31 A brief explanation of the measuring way, i.e., the sit-on-bottom stationary measurement, is given here. When performing the sit-on-bottom stationary measurement, HAIMA ROV lands fixed at the shoulder of the CRC mountain. As shown in Figure 7a, the electronic processing unit of PPPAAP17 marked with a red dashed box lies at the bottom of the HAIMA Remotely Operated Vehicle (ROV). The transmitter transducer array and receiving transducer array of PPPAAP17 marked with a blue dashed box is mounted on the hydraulic system of HAIMA ROV, as shown in Figure 7b. It can be seen that the vertical distance from the array to the ground could be adjusted. For the stationary measurement, the transmitter transducer array and receiving transducers mounted on HAIMA ROV are about 300 mm away from the seabed in the vertical direction. The parametric acoustic system, PPPAAP17, transmits the signal for a fixed spot every 20 ms and it samples and keeps the first 2 ms data for each pulse. The data of the primary channel and the The parametric acoustic system, PPPAAP17, transmits the signal for a fixed spot every 20 ms and it samples and keeps the first 2 ms data for each pulse. The data of the primary channel and the secondary channel were separately received and stored. It should be mentioned that the transducer of the secondary channel outstretches 20 mm compared with the primary one, which makes the offset for the reflection instant of the dual-channel. Such a feature causes a fixed time delay difference for the dual-channel when receiving the acoustic echo signal, which could be assumed as the system error. It could be measured in advance and then compensated in the algorithm.
The feasibility of the proposed method and the comparison with the reference method [4] will be given by processing the data acquired by the one stationary measurements, i.e., Sit-on-bottom 1(S1). The real scene captured by the underwater camera is given in Figure 8. secondary channel were separately received and stored. It should be mentioned that the transducer of the secondary channel outstretches 20 mm compared with the primary one, which makes the offset for the reflection instant of the dual-channel. Such a feature causes a fixed time delay difference for the dual-channel when receiving the acoustic echo signal, which could be assumed as the system error. It could be measured in advance and then compensated in the algorithm. The feasibility of the proposed method and the comparison with the reference method [4] will be given by processing the data acquired by the one stationary measurements, i.e., Sit-on-bottom 1(S1). The real scene captured by the underwater camera is given in Figure 8. We select 500 records of the primary data and secondary data, whose duration time is 10.0 s. To better illustrate the features of the reflected signal, we take the 500th primary signal for depiction. As shown in Figure 9a, the primary signal has a good SNR value of 24.88 dB. The signal captured with the red dashed box contains the reflection instant of the top surface. As shown in Figure 9b, the normalized primary signal is zoomed with the extracted envelopes by different methods. Although not very obvious, it still can be seen that the envelope extracted with the wavelet method is better than the Hilbert method at the aspect of the robustness when the high-frequency distortion exists from the areas marked with the dotted box. In addition, we can see that two peaks of the envelope appear at the time of 0.4192 ms and 0.4266 ms, respectively, which means that below the top surface exists another layer. The phenomenon corresponds to the multilayer on the top surface. Figure 9c reveals that the secondary signal also has a good SNR value of 31.13 dB, while the signal has obvious interference which occurs at a very short time after transmission. The interference is marked with the dashed box (A) and the two major reflections are surrounded by the dashed box (B), as shown in Figure 9d. We select 500 records of the primary data and secondary data, whose duration time is 10.0 s. To better illustrate the features of the reflected signal, we take the 500th primary signal for depiction. As shown in Figure 9a, the primary signal has a good SNR value of 24.88 dB. The signal captured with the red dashed box contains the reflection instant of the top surface. As shown in Figure 9b, the normalized primary signal is zoomed with the extracted envelopes by different methods. Although not very obvious, it still can be seen that the envelope extracted with the wavelet method is better than the Hilbert method at the aspect of the robustness when the high-frequency distortion exists from the areas marked with the dotted box. In addition, we can see that two peaks of the envelope appear at the time of 0.4192 ms and 0.4266 ms, respectively, which means that below the top surface exists another layer. The phenomenon corresponds to the multilayer on the top surface. Figure 9c reveals that the secondary signal also has a good SNR value of 31.13 dB, while the signal has obvious interference which occurs at a very short time after transmission. The interference is marked with the dashed box (A) and the two major reflections are surrounded by the dashed box (B), as shown in Figure 9d The first peak (b1) marked with an arrow corresponds to the top reflection and the second peak (b2) is the corresponding tail that is composed of many scattered waves within the footprint when the transmitting wave arrives at the surface. The third peak (b3) and the third peak (b4) correspond to the bottom reflection and scattering, respectively. It is worth noting that the time difference between the first reflection of the dual-channel, which corresponds to the fixed offset and will be compensated at the final step for thickness estimation. Compared with the primary channel, the secondary channel only has one peak at the top surface since the low-frequency envelope has a stronger penetrating capability.
To compare the performance of different methods, we present the results of the arrival time of top and bottom surface and the estimated thickness in Figure 10a,b, respectively. As shown in Figure  10a, the two lines marked with A and B denotes the results obtained by the reference method. Due to the fact that the signal is severely interfered, the simple peak search algorithm improperly takes the first two peaks of relatively large amplitude. The proposed method is more robust even when the interference occurs. Besides, as marked by the black dashed circle, it is very obvious that the proposed method has more stable results than others with no multilayer logic.
The comparisons of estimated results and the visible thickness of the CRCs are drawn in Figure  10b. It is worth noting that the estimation of the methods with no multilayer logic has biases than the others since they always extract the first peak that occurred at the top surface. By calculation, we have known that in the presence of noise the reference method could only obtain the wrong estimation which exceeds 500mm. Before we could perform the comparison, the interference is first removed. As shown in Figure 10b, we could obtain the estimation of the reference method is about 220 mm. The first peak (b1) marked with an arrow corresponds to the top reflection and the second peak (b2) is the corresponding tail that is composed of many scattered waves within the footprint when the transmitting wave arrives at the surface. The third peak (b3) and the third peak (b4) correspond to the bottom reflection and scattering, respectively. It is worth noting that the time difference between the first reflection of the dual-channel, which corresponds to the fixed offset and will be compensated at the final step for thickness estimation. Compared with the primary channel, the secondary channel only has one peak at the top surface since the low-frequency envelope has a stronger penetrating capability.
To compare the performance of different methods, we present the results of the arrival time of top and bottom surface and the estimated thickness in Figure 10a,b, respectively. As shown in Figure 10a, the two lines marked with A and B denotes the results obtained by the reference method. Due to the fact that the signal is severely interfered, the simple peak search algorithm improperly takes the first two peaks of relatively large amplitude. The proposed method is more robust even when the interference occurs. Besides, as marked by the black dashed circle, it is very obvious that the proposed method has more stable results than others with no multilayer logic.
The comparisons of estimated results and the visible thickness of the CRCs are drawn in Figure 10b. It is worth noting that the estimation of the methods with no multilayer logic has biases than the others since they always extract the first peak that occurred at the top surface. By calculation, we have known that in the presence of noise the reference method could only obtain the wrong estimation which exceeds 500mm. Before we could perform the comparison, the interference is first removed. As shown in Figure 10b, we could obtain the estimation of the reference method is about 220 mm.
To quantitatively assess the results, we list the statistics in Table 5. It can be seen that the simple method with Hilbert/No Multilayer/has a large bias of 11.40 mm and a standard deviation of 17.30 mm given that the visible thickness (#) of 255.00 mm. In contrast, using the simple method with Hilbert/No Multilayer/R we could optimize the stability and obtain the estimation with the standard deviation of 1.23mm. We have known that wavelet regression can optimize the results of all methods. As for the reference method, the estimation has the largest bias and also a large standard deviation, suggesting that it fails to resolve the multilayer case and only achieves the poor precision based on only one channel. Finally, we could see that the proposed method with Wavelet/No Multilayer/R obtains the most accurate results with the mean value of 256.89 mm and the standard deviation of 1.20 mm. To quantitatively assess the results, we list the statistics in Table 5. It can be seen that the simple method with Hilbert/No Multilayer/has a large bias of 11.40 mm and a standard deviation of 17.30 mm given that the visible thickness (#) of 255.00 mm. In contrast, using the simple method with Hilbert/No Multilayer/R we could optimize the stability and obtain the estimation with the standard deviation of 1.23mm. We have known that wavelet regression can optimize the results of all methods. As for the reference method, the estimation has the largest bias and also a large standard deviation, suggesting that it fails to resolve the multilayer case and only achieves the poor precision based on only one channel. Finally, we could see that the proposed method with Wavelet/No Multilayer/R obtains the most accurate results with the mean value of 256.89 mm and the standard deviation of 1.20 mm.
After analyzing the results, we can conclude that the proposed method has several advantages over the reference method. Firstly, the envelope extracted with wavelet transform tends to be more robust even when the SNR decreases. Secondly, based on the multilayer model and local continuity assumption, the proposed method could recognize the correct peak of the top surface reflection even though the value of the peaks varies with the underwater circumstance. In contrast, the reference method could only obtain the first refection peak since it exceeds the threshold, thus leading to  After analyzing the results, we can conclude that the proposed method has several advantages over the reference method. Firstly, the envelope extracted with wavelet transform tends to be more robust even when the SNR decreases. Secondly, based on the multilayer model and local continuity assumption, the proposed method could recognize the correct peak of the top surface reflection even though the value of the peaks varies with the underwater circumstance. In contrast, the reference method could only obtain the first refection peak since it exceeds the threshold, thus leading to inaccurate estimation. In addition to this, the proposed method has a better resolution considering that the sampling frequency of the primary channel is much larger than the secondary one. Finally, using the dual-channel information the proposed method could simply remove the interferences while the reference method could not. From multiple aspects, the results have proven the feasibility of the proposed method.

Conclusions
This paper described the importance of the CRCs as a potential mineral necessity source and the necessity of mapping and estimating the volumetric distribution of deep-sea mineral. The prototype acoustic probe, i.e., PPPAAP17, was described in detail, aiming to collect and process the data both for the thickness measurement and for the material recognition in the future. As one of the highlights, the prototype was designed with the dual-channel for receiving the primary and secondary signal. Considering that the signal quality is degraded by the system interference and ambient noise, the thickness extraction algorithm was proposed with some improvements, such as the wavelet-based envelope extraction method, the adaptive strategy of the determination of the bottom and top surface arrival time based on the dual-channel information, and the wavelet regression to reduce measuring noise. The algorithm is suitable for the assumption that the CRCs are with the structure of the multilayers at the top surface and one single layer at the bottom surface. Such an assumption is reasonable since it is suitable for some of the real scenes of the CRCs. The algorithm considers the connection of the dual-channel signal, i.e., the features of the CRCs will be shared. The wavelet regression considers the fact of the local invariability of a large amount of data. The laboratory experiment has been demonstrated in the tank and has validated the effectiveness of the processing method. Besides this, the real experiment on the China Ocean 51th voyage in the Western Pacific Ocean on Aug 30, 2018, was carried out by using the sit-on-bottom stationary measurement and the data were obtained. The processing results have shown that the proposed algorithm has very good performance as a standard variation of 1.20 mm (S1) with wavelet regression. As for further research, we will consider applying the Empirical Mode Decomposition (EMD) method to the problem mentioned above. The EMD method decomposes a signal into modes having physical meaning, and it allows for envelope extraction in each mode separately. Besides, the work could be carried out at the aspect of building the more complex thickness model describing the multilayer feature at the bottom surface.