Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Network Optimisation and Performance Analysis of a Multistatic Acoustic Navigation Sensor
Next Article in Special Issue
Development of Scientific Fishery Biomass Estimator: System Design and Prototyping
Previous Article in Journal
A Review of Microelectronic Systems and Circuit Techniques for Electrical Neural Recording Aimed at Closed-Loop Epilepsy Control
Previous Article in Special Issue
Chladni Figures in Modal Analysis of a Double-Panel Structure
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mapping Small-Scale Horizontal Velocity Field in Panzhinan Waterway by Coastal Acoustic Tomography †

1
Ocean College, Zhejiang University, Zhoushan 316021, China
2
Laboratory for Marine Geology, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266061, China
*
Author to whom correspondence should be addressed.
This paper is an extension version of the conference paper: Guo, Y.; Huang, H.; Wu, J.; Bo Xu Zhu, C.; Wang, Z.; Feng, R. Short-range velocity field observation by coastal acoustic tomography. In Proceedings of the Oceans 2019 mts/IEEE Seattle IEEE, Seatlle, WA, USA, 27–31 October 2019.
Sensors 2020, 20(19), 5717; https://doi.org/10.3390/s20195717
Submission received: 20 August 2020 / Revised: 20 September 2020 / Accepted: 5 October 2020 / Published: 8 October 2020
(This article belongs to the Special Issue Intelligent Sound Measurement Sensor and System)

Abstract

:
Mapping small-scale high-precision velocity fields is of great significance to oceanic environment research. Coastal acoustic tomography (CAT) is a frontier technology used to observe large-scale velocity field in the horizontal slice. Nonetheless, it is difficult to observe the velocity field using the CAT in small-scale areas, specifically where the flow field is complex such as ocean ranch and artificial upwelling areas. This paper conducted a sound transmission experiment using four 50 kHz CAT systems in the Panzhinan waterway. Notably, sound transmission based on the round-robin method was recommended for small-scale CAT observation. The travel time between stations, obtained by correlation of raw data, was applied to reconstruct the horizontal velocity fields using Tapered Least Square inversion. The minimum net volume transport was 8.7 m3/s at 12:32, 1.63% of the total inflow volume transport indicating that the observational errors were acceptable. The relative errors of the range-average velocity calculated by differential travel time were 1.54% (path 2) and 0.92% (path 6), respectively. Moreover, the inversion velocity root-mean-square errors (RMSEs) were 0.5163, 0.1494, 0.2103, 0.2804 and 0.2817 m/s for paths 1, 2, 3, 4 and 6, respectively. The feasibility and acceptable accuracy of the CAT method in the small-scale velocity profiling measurement were validated. Furthermore, a three-dimensional (3-D) velocity field mapping should be performed with combined analysis in horizontal and vertical slices.

1. Introduction

The velocity of seawater is one of the most fundamental basic data used in ocean science research and engineering [1]. Ocean acoustic tomography (OAT) is an advanced oceanographic technology first proposed by Munk and Wunsch in 1979 [2]. This technique can make a simultaneous mapping (snapshot) of time-varying subsurface structures of current velocity and sound speed using an underwater sound channel [3,4]. Coastal acoustic tomography (CAT) is the further advancement of OAT that has been developed in coastal areas. The CAT system was developed by the Hiroshima University Group as a new observation tool utilized in oceanographic observation from the shore since 1995 [5]. Consequently, a growing number of studies on this technology have been conducted. Initially, more comprehensive, multi-station measurements with CAT systems (more than five stations) were conducted by the Hiroshima University Group initially [6,7,8,9]. Of note, CAT has been highly applied in measuring various coastal–sea phenomena, including the tidal current [10,11], residual current [12,13,14], internal solitary waves [15], internal tides [16], tidal bores [17] and the coastal upwelling [18,19] et al. This was attributed to its advantages, including low cost, compact system, simple instrument operation and easy execution on board. Additionally, it is worth noting that CAT reconstructs velocity structures with few sensors covering a wide range of survey zone [11], the traditional observation range (station-to-station distance), which might be less than 50 km in the coastal seas shallower than 100 m [2].
Many CAT applications in large scale areas have achieved satisfactory outcomes. For instance, a 2009 study Zhu et al. [20] performed the CAT experiment for mapping the tidal currents in Zhitouyang Bay with seven CAT stations. Here, the volume transport for the entire tomography domain, considered as 22 × 22 km, was calculated to be nearly equal to zero, implying that the observational errors were quite small. In 2015, an experiment using 11 CAT stations ranged between 2.4 km and 16.1 km was performed by Zhang et al. [11] in the Dalian Bay where the detailed horizontal distributions of tidal and residual currents were well reconstructed with the standard inversion method. Nevertheless, the data assimilation was not to be implemented since the tomography field was not well surrounded by shorelines. In 2018, Chen M et al. [21] provided CAT with the latest feature to realize the real-time data transmission from an offshore subsurface station to a shore station with mirror-transponder functionality. Furthermore, a 2019 investigation by Syamsudin et al. observed the subsurface structures of internal solitary waves in the Lombok Strait with two bottom-moored CAT systems over a path length of 18.286 km [15]. In this study, the four-layer subsurface temperature structures by the inverse analysis were reconstructed where the phase relations of internal solitary waves to tides were initially observed. However, the overestimation of travel times occurred since the range-dependent distribution of the sound speed was not considered in the simulation. Additionally, data assimilation acted as an effective partner of CAT in predicting coastal–sea environmental variations [2]. Zhu et al. [22] presented the application of an unstructured triangular grid to the Finite-Volume Community Ocean Model, to assimilate CAT data, showing that missing data exhibited lesser impact on assimilation than on inversion. These large-scale measurements using CAT in bays, ports, semi-enclosed bays, and inland seas were significantly considered and achieved satisfactory progress [23,24,25]. On the other hand, the applications of CAT in small-scale areas are largely understudied.
Since 2010, a group of Kawanisi proposed that the fluvial acoustic tomography (FAT) system extends the applications of CAT to even shallower waters ranging from mountainous rivers to the mouth of estuaries in coastal regions (maximum 10 m deep) [26,27]. As a consequence, the application of CAT in shallow waters such as the marine ranch, artificial upwelling, and Autonomous Underwater Vehicle (AUV) observation has been widely investigated, and particularly the short-range velocity measurements have gradually been an important focus of research. The observation area of those specific targets is usually within 1 km (much smaller than the traditional CAT observation range) [28,29,30]. Luo et al. [31] conducted a field experiment in Delaware Bay with two CAT stations at a distance of 387 m to observe the average velocity along the propagation path. Elsewhere, Li et al. [32] experimented with a laboratory pool using two CAT systems at a distance of 22.73 m and observed a unidirectional flow of 0.8 m/s, which was in line with current findings from simulation studies. Satisfactory results were obtained in laboratory tests, nonetheless, it is still necessary to conduct tests on seas where the environment is more complex. Although the CAT has been used in the observation of the flow field in the small-scale areas, the available evidence remains scanty.
Complex hydrological environment, high noise level, and complicated flow field are among the observation environments that markedly influence sound propagation in the small-scale area compared to the large-scale CAT observation. At present, there is no systematic study on CAT in small-scale areas. Moreover, there are few existing CAT observations involving the snapshots of two-dimensional (2-D) velocity fields in the field experiment, which is vital to intelligent monitoring and management for small-scale areas.
Due to the large environmental noise, M-sequence with an efficient autocorrelation property was generally adopted. The observation results would be directly influenced by the number of cycles used by an M-sequence bit, the order of M-sequence, the carrier frequency, and the signal repetition times [26,33]. Therefore, a reasonable signal design is vital in short-range CAT observation. In the previous study, we designed a research outline and obtained part of the short-range velocity field by CAT, which validated the feasibility of the preliminary assumption [34]. This study conducted a quantitative analysis to acquire more comprehensive information on the small-scale velocity field. Moreover, quantitative error analysis and systematic discussion were performed to confirm the reliability of the results. A small-scale experiment with four CAT systems for reconstructing 2-D horizontal velocity fields was conducted on January 20, 2019, in the Panzhinan waterway. The distances between CAT systems in the survey zone were significantly short with the shortest distance being 85.1 m. To obtain high resolution in short-range observation, a high carrier frequency of 50 kHz was selected, which was the first short-range sea test conducted with a 50 kHz CAT system. A sound transmission based on the round-robin method was proposed to adopt the short-range CAT with the 12th order M sequence, which can be beneficial to ensure the normal signal transmission and obtain high signal-to-noise-ratio (SNR) for small-scale CAT observation in the high-noise environment. The feasibility and considerable accuracy of the CAT method in the small-scale velocity profiling measurement were verified by several methods.
The rest of the paper is structured as follows. In Section 2, the inversion method used to calculate the velocity fields in a horizontal slice is introduced. Experimental settings and ray simulation are also discussed. Section 3 includes the correlation of the data and the results of horizontal velocity fields. The net volume transport and comparisons of the CAT results are performed to examine the accuracy. The concluding remarks are provided in Section 4.

2. Method

2.1. The Forward and Inverse Problems

When the sound travels between two subsurface points, it is influenced by the current velocity and sound speed along a transmission path, resulting in different travel times in different directions. A pair of stations in the CAT experiment was composed of two acoustic stations (S1 and S2) (Figure 1). The reciprocal sound transmission between acoustic stations S1 and S2 was sketched in Figure 1. The transmission sound line between S1-S2 is projected onto the horizontal plane, and u is the projection of the velocity v in the direction of the transmission sound line. The travel time t can be formulated as follows:
t = d s C 0 + Δ C + v n
where C0 represents the reference sound speed, v is the velocity vector and n is the unit vector along the transmission path. ΔC is the deviation from C0 caused by temperature variations. Δ C C 0 and u = v n C 0 are usually satisfied in the sea. The Taylor expansion is adopted to the denominator of Equation (1):
t = d s C 0 ( 1 Δ C ± v n C 0 + ( Δ C ± v n ) 2 C 0 2 )
Taking the second Taylor expansion:
t = d s C 0 ( 1 Δ C ± v n C 0 ) + Δ ( d s ) C 0 ( 1 Δ C ± v n C 0 ) = ( Δ C ± v n )   d s C 0 2 + 1 C 0 ( 1 Δ C ± v n C 0 ) Δ ( d s )
The second or higher-order terms can be ignored. Thus:
t ± = d s C 0 ( 1 Δ C ± v n C 0 )
By taking the difference of t+ and t, the differential travel time Δ τ i becomes:
Δ τ i = t + t 2 v i C 0 2 d s
The 2-D velocity fields on a horizontal slice can be represented by the stream function, which expanded into a Fourier function series [22]:
ψ ( x , y ) = a x + b y + k = 0 N x l = 0 N y { A k , l c o s 2 π ( k x L x + l y L y ) + B k , l s i n 2 π ( k x L x + l y L y ) } = j = 1 ( N x + 1 ) ( N y + 1 ) D j Q j ( x , y )
where D = { D j } = { a , b , A 00 , , B 00 , A 01 , B 01 , , A N x N y , B N x N y } and
Q ( x , y ) = { Q j } = { x , y , 1 , 0 , c o s 2 π y L y , s i n 2 π y L y , , c o s 2 π ( N x x L x + N y y L y ) , s i n 2 π ( N x x L x + N y y L y ) }
.
Equation (6) reduces in a matrix form:
y = E x + e
where y = { Δ τ i } the travel time vector is obtained in the experiment, x = { D j } is the unknown coefficient vector to be solved, e is the error vector, and E = { 2 C 0 2 j = 1 ( N x + 1 ) ( N y + 1 ) 0 L i Q j c o s θ 2 i d x } is the coefficient matrix. There is an inevitable error vector e in practice at the right end of Equation (7), which is caused by the observation error and the inaccurate model description error.
Equation (7) is generally ill-posed often manifested in the instability of the solution, which cannot be solved using the usual linear algebraic methods. Therefore, the Tapered Least Square method accompanied by the L-curve method is required to solve the inversion problem [17]. The cost function J, which comprised data misfit and smoothness measure of solution vector x, is given by Equation (8).
J = n T n + α 2 x T x = ( y E x ) T ( y E x ) + α 2 x T x
The coefficient α represents the damping factor updated during each sound transmission process to make solutions much flexible to trace the dynamic underwater environments.
By singular value decomposition (SVD), transform matrix E reduces:
E = i = 1 N s λ i u i v i T = U Λ V T
where u and v represent the singular vector of E and λ is the corresponding eigenvalue. U = [ u 1 u 2 u N s ] , V = [ v 1 v 2 v N s ] , and
Λ = ( λ 1 0 0 λ N s )
where Ns depicts the number of nonzero singular values. By taking the x-derivative of J and setting it to zero, an expected solution can be obtained to minimize J through a singular value decomposition, which is expressed by Equation (10).
x ^ = i = 1 N s λ i ( u i T y ) v i λ i 2 + α 2
The expected error vector e ^ is
e ^ = y E x ^ = { I E ( E T E + α 2 I ) 1 E T } y
In this paper, the Tapered Least Square method was also used to damp or filter out noise in the high-frequency component thereby stabilizing the obtained solution.

2.2. The Influence of Position Accuracy

As shown in Figure 1, the reciprocal travel times along a transmission path are formulated as:
t 1 = L C m + V m
t 2 = L C m V m
where t1 and t2 denote the travel times from S1 to S2 and from S2 to S1, respectively. Cm and Vm denote the range-average sound speed and range-average velocity along the path, respectively. L is the horizontal distance between the S1 and S2. Solving the coupled Equations (12) and (13), we obtain:
V m = C m 2 2 L Δ t
By calculating the range-average velocity Vm of each path, a matrix can be formed to solve the velocity of any point in the observation domain. In addition, we can prove the correctness of CAT inversion results by comparing Vm and the projection of the 2-D flow velocity field in this path, where Δt = t2t1. t1 and t2 are nearly equal to tm. Taking the total derivatives of Equation (14), we obtain:
δ V m V m = δ L L + δ ( Δ t ) Δ t
When the global positioning system (GPS) clock is used, the clock accuracy is as small as 0.6 μs and the second term on the right-hand side of Equation (15) is negligibly small. For L = 100 m and δL = 2 m [coming from GPS positioning errors], δ L L = 0.02 and δV = 0.02Vm, showing that positioning error is not an important factor of velocity measurement. This means that the effect of positioning errors on velocity errors is negligibly small in the small-scale CAT observation. In contrast, temperature field measurement requires position correction using the method proposed by Zhang [11].

2.3. Experiment Settings

A CAT field experiment lasting 4 hours was conducted in the Panzhinan waterway near the Zhairuoshan Island, Zhoushan. The map of Zhairuoshan Island with adjacent regions and the CAT station array is shown in Figure 2. The northern part of the Zhairuoshan Island was magnified in the right panel, which showed the 4 CAT stations (S1, S2, S3, and S4) and one conductivity–temperature–depth (CTD) cast location in the black and red solid circles, respectively. The yellow solid line was the transmission path of the sound line. The upper left of the right panel contained the shipboard acoustic Doppler current profiler (ADCP) and CTD, respectively. The water depths varied from 20 m at the northeastern part to 30 m at the southwestern part within the observation area. The ADCP was fixed on the ship with a fixator for collecting the velocity and bottom terrain data in a bottom-tracking mode. Additionally, limited information on the flow field was collected by the shipboard ADCP survey along several ship tracks crossing the observation region. The CTD cast was performed to measure the sound velocity and temperature profile in the observation region.
The CAT systems and batteries were arranged on the shore base (Figure 3). The acoustic transceivers were hoisted by a rope under the water around 1-2 m together with a counterweight to prevent large position drift of acoustic transceiver due to the influence of the current. Whereas, the major components of the system, such as electronic housing, batteries, and the Global Positioning System (ANN-MS-0-005 GPS) antennas were placed on the shore. All stations were equipped with a GPS of high clock accuracy (less than 0.6 μs error) to synchronize the clock timing during transmission, reception, and AD conversion to correctly reconstruct the velocity fields. In addition, GPS provided clocked signals and provided station positioning data. The distances between the stations are listed in Table 1.
In the short-range experiment, the signal design was vital and needed to be comprehensively considered. In the survey zone, it was inevitably interfered by various noises such as seabed topography and ship navigation, which required high precision and resolution of CAT. To obtain high resolution and high precision in small-scale area, a transceiver transducer with a central frequency of 50 kHz was adopted by the CAT systems. The acoustic transceiver used in the experiment configured signal transceiver modules for both transmitting and receiving acoustic signals, which implied more convenience for system integration of CAT.
Time resolution tr can be represented as:
t r = Q f
where Q represents the number of cycles used by an M-sequence, and f represents the carrier frequency (i.e., 50 kHz).
The transmission loss was increased with the increase of the frequency. For sake of the correct reception of signals and high SNR, the signal design of the short-range CAT test was indispensable. The length of the Mn is (2n − 1), and the time length of the signal ts can be expressed as:
t s = ( 2 n 1 ) Q f
where n represents the order of M sequence. Equation (17) indicated that the signal length increases with the increase of the order of M sequence.
S N R = 20 log 2 n 1
In Equation (18), the value of SNR is positively associated with the order of M sequence. Based on the preparatory experiment, different signals were tried to extract and match the multi-peak value. In addition, considering the significant environmental noises generated by ADCP-ship navigation and marine power generation platform, strong current and turbid water quality, one period of the 12th order M sequence was used to modulate with the 50 kHz carrier signal. During the experiment, the signal was transmitted at intervals of every 4 minutes from each acoustic transceiver and the signal repetition times (N) was 8.
According to Equation (17), when the 12th order M sequence was used, Q = 3, N = 8, ts = 1.9656 s. C0 = 1486.5 m/s, the transmission distance should be 2921.86 m as synchronous transmission used. Since the acoustic transceiver cannot perform the function of transmitting and receiving signals concurrently, the distance between CAT stations should be larger than the transmission distance to ensure the normal operation of the synchronous transmission. Therefore, in the case of short-range tests, the ordinary synchronous transmission cannot be obtained using a high order of M sequence.
As shown in Table 2, a sound transmission based on a round-robin method was proposed to adopt the short-range test, which can benefit the normal reception of signals and obtain high SNR. Moreover, it solves the problem that each acoustic transceiver cannot perform the function of transmitting and receiving signals concurrently due to the short distance of the CAT station.
Arrival signals transmitted from each acoustic station were identified from the multiple stations based on M-sequences, which was allocated before the experiment. The sampling time window was enlarged to increase the data sampling number and decrease the influence of position drifting. Through the complex demodulation, the received signals were stored in the internal micro SD card. By taking a multiperiod transmission through the ensemble average of the successively received data, the SNR of the received data can be further improved.

2.4. Ray Simulation

Figure 4a shows the ray patterns of S1-S3 in the range-independent ray simulation, which used the range-averaged sound–speed profile obtained from the CTD data during the CAT experiment. The three ray-patterns include: (1) the DR (direct) rays; (2) the SBR (Surface-bottom reflected) rays; and (3) the SR (Surface reflected) rays. The rays were almost straight lines in short-distance transmission, which can be applied to the approximate simulated station-to-station distances using the ray length (Figure 4a). According to the ray simulation, it was found that the suitable rays in which reflection times were 1 or less could not be received when the depth was more than 25 m.
The first group of rays passing near the surface constructed the first arrival peaks, the characteristic of which map the 2-D horizontal velocity fields. The distances between the stations were calculated by the GPS data from the CAT systems. Despite the inevitable errors due to a few factors, including the GPS positioning that might cause distance error, the impact of position drifting on the horizontal flow field observation was negligibly small. Thus, even though extremely high-precision position information was not attained, the 2-D horizontal velocity fields could be reconstructed.
Based on pre-investigations, the velocity of seawater in the Panzhinan waterway was fast causing the seawater at different depths to be well mixed, there was no apparent temperature stratification in the survey zone and the temperature of the survey zone changed little with the depth. Therefore, the information obtained from a 1-point CTD cast should be sufficient to infer the small-scale tomographic domain in the present experiment. The CTD profile data including temperature and sound speed profiles are shown in Figure 4b and c, respectively. The temperature increased with depth due to surface cooling in winter followed by the sound speed. Notably, water in the observation region was well homogenized horizontally because of the strong tidal currents. Hence, temperature stratification was significantly weak and assumed to be constant. Though the variation of salinity has a certain influence on the sound velocity, except for the river mouth, salinity variation is nearly equal to zero [2].

3. Result and Discussion

3.1. The Signal Correlation

Travel times were determined by arrival peaks in the correlation pattern obtained from the cross-correlation of the received data with one period of the M sequence used in the transmission. The largest arrival peak method was used to identify and trace the significant arrival peaks in correlated patterns of received data (i.e., identify the largest arrival peak over the SNR threshold). Figure 5 illustrates the stack diagram of the correlation patterns obtained sequentially for each transmission path under the SNR threshold of 20. The number of significant arrival peak was only one, and the first arrival peak also was the largest arrival peak. The first-arrival peaks were identified although the peak heights gradually decreased with battery power exhaustion (Figure 5).
The SNR was the ratio of the peak height to the mean height of the received data in this work [2]. From this definition, SNR was slightly underestimated because signal peaks were included in an estimate of the noise level, and the travel time can be accurately identified. This definition of SNR was convenient for comparison and peak extraction. The travel time was obtained by identifying the largest arrival peaks in the cross-correlation pattern of the received data with the transmitted signal. As shown in Figure 5, the identification of the cross-correlation peaks from S2-S4 (path 5) was significantly unsuccessful. The red circles (arrival peaks) were randomly distributed in time windows. This could be due to the renewable energy platform in the ocean between S2 and S4, which might have negatively influenced signal transmission. Thus, path 5 was not used in the later inversion.

3.2. Horizontal Velocity Fields

Figure 6 shows the time plots of the differential travel time for five station pairs. The differential travel time ranged from 0.02 to 0.6 ms, corresponding to variation in range-average velocity with a dependence on the distance between each station pair. The differential travel time was used to reconstruct horizontal velocity distributions and calculate the range-average velocity as comparative data.
The differential travel time data were used in the inverse analysis. The stream function for displaying the 2-D current field was expanded into the Fourier function series, and the unknown expansion coefficients were determined by the Tapered Least Square method accompanied by the L-curve method. A satisfactory inversion outcome was obtained when L-curve was folded at the maximum curvature point. The inversion results provided snapshots of horizontal velocity distributions at different times. The velocity structures between 12:21:00 and 13:04:00 on January 20, 2009, are shown in Figure 7. The current was strong with a maximum velocity of 1.52 m/s in the observed region and gradually changed with time. First, the velocity of the flow field decreased followed by an increase from 12:21 to 13:04. The velocity in the whole flow field was lowest (0.25 m/s) at 12:43. The northward current in the central region was separated into NW and NE currents in the left-side and right-side regions, respectively. Overall, the flow field distribution in the observation region was evident and it rapidly changed.
The ADCP velocity data with vector plots along the tracks are diagramed in Figure 8. In Figure 8b, the thin arrow in the lower-left is the 1 m/s velocity fields, and the inversion is consistent with the ship-mounted ADCP velocity data, which corresponds to the time when the ADCP ship was on each track (Table 3).
The sound signal was transmitted at intervals of every 4 mins from each acoustic transceiver in the CAT system, while the sampling frequency of ship ADCP was set to 1 s generally. Since the sampling frequencies between CAT and ADCP were different, the directly instantaneous comparison between CAT and ADCP results was not significant enough. It was reasonable to calculate the ADCP average data for comparison. The comparisons of the ADCP average data and the range-average velocity Vm are shown in Figure 9, taking S1–S3 and S3–S4 as examples, which had high coincidence with the ADCP measurement path. The ADCP measurement data along S1–S3 (path 2) and S3–S4 (path 6) are shown in Figure 9a,c, respectively. After averaging the ADCP measurement data, it was found that the ADCP average data and the corresponding range-average velocity Vm exhibited satisfactory consistency in Figure 9b,d, respectively. The relative error of S1–S3 (path 2) and S3–S4 (path 6) were 1.54% and 0.92%, respectively.
A comparison between ADCP and the projection of the 2-D flow field obtained by the CAT inversion on the path was also necessary to validate the results. The range-average velocity Vm calculated by differential travel time (Δt) was compared to the CAT inversion results. Figure 10 illustrates a comparison of the velocities of five different paths at nine moments, which corresponded to the horizontal velocity fields (Figure 7). Each graph in Figure 10 shows the average velocity of five different paths at a certain time, where the blue hollow and the red triangle solid lines represent the inversion results and the range-average velocity Vm, respectively. The omission of path 5 was due to the influence of the oceanic renewable energy platform (Figure 10). The platform affected signal transmission by either stopping a direct sound path or generating significant noise. From the results, the SNRs of the cross-correlation peaks were weak along the S2–S4 paths (Figure 5).
The root-mean-square difference (RMSD) between the CAT inversion and the range-average velocity Vm were 0.5163, 0.1494, 0.2103, 0.2804 and 0.2817 m/s for path 1, 2, 3, 4 and 6, respectively. The RMSD of the S1-S2 (path 1) was relatively large, perhaps because there was an oceanic renewable energy platform between S2 and S4. The sound signals of S2-S4 (path 5) were directly blocked by the oceanic renewable energy platform so that the effective time information was not obtained, while the data quality of S2-S3 (path 4) was acceptable (the relative error was 0.2804). Therefore, it was not only the influence of oceanic renewable energy platform but also the layout of stations. The layout of the stations was not precise enough since the sound path was too dense between three stations (S1, S2, and S4) and the signal quality was poor due to the influence of the platform. In general, the relative errors of the inversion results were approximately 10% to 20% which were acceptable, verifying the feasibility of small-scale CAT observation (Figure 10).

3.3. Cross-Section Volume Transport

It is reasonably assumed that vertical velocity is negligibly small in most regions of the sea. Though there are fluctuations for a long time, the sea level could be considered as uniform in the instantaneous state. Furthermore, the 3D structure of the current field was vertical homogeneous since the vertical fluctuation was actually negligibly small. The accuracy of the CAT data was validated by the equivalence of volume transportation.
The volume transport across the entire section was estimated using the Vp and the vertical section areas:
Q = V p A
where Vp represents the velocity perpendicular to the transmission path, and A represents the cross-sectional area of the vertical section along a sound transmission path.
To examine the volume balance inside the closed domain, the net volume transport across the 2-m perpendicular upper part of the transects at the periphery was calculated (Table 4). The inflow was set to be positive whereas the outflow was negative. Flow rates calculated corresponding to six different moments are summarized as net volume transport in Table 4. The inflow occurred through the S3–S4 (path 6) while the outflow occurred through S1–S2 (path 1) and S1–S3 (path 2). Only S2–S4 (path 5) fluctuated from outflow (minus) to inflow (plus). The minimum net volume transport estimated for the entire topographic domain reached 8.7 m3/s at 12:32, which was 1.63% of the total inflow volume transport. The average net volume transport of −12.6 m3/s was insignificant compared to the maximum inflow volume transport of 619.94 m3/s. The total net volume transport was nearly equal to zero, showing a balance of transport between the inflow and outflow. Generally, our findings indicate sufficient accuracy of the present small-scale CAT experiment.
Notably, the CAT results were consistent with ADCP results despite a few deviations. Several possible reasons might account for the observed error, including: (1) The lower SNR at the later period of the experiment might be due to the exhaustion of battery power. (2) The layout of the stations was not precise enough since the sound path was too dense in three stations (S1, S2, and S4) and the signal quality was poor due to the influence of the platform. (3) The inevitable position drift of the acoustic transceiver with the current showed an adverse effect on the precision. (4) The amount of data was too small to be processed by the moving average.

4. Conclusions

In conclusion, we performed the CAT experiment for mapping small-scale velocity structures using four CAT stations. Sound transmission based on the round-robin method was proposed for small-scale CAT observation. Moreover, the reciprocal sound transmission was performed along six transmission lines spanning the four CAT stations, while CTD casts were conducted to obtain temperature profiling as well as sound velocity data. Range-independent ray simulation based on the CTD data was performed to determine ray paths. The tomographic inversion based on the Tapered Least Square method was adopted to reconstruct horizontal velocity distributions with a differential travel time, and obtain continuous velocity variations mapping in a horizontal slice. The net volume of transport was calculated to examine the volume balance. We found that the estimated minimum net volume transport was 8.7 m3/s at 12:32, which was 1.63% of the total inflow volume transport. The average net volume transport was quite small compared to the maximum inflow volume transport (−12.6 m3/s vs 619.94 m3/s), indicating acceptable observational errors. The relative errors of the range-average velocity calculated by (S1-S3 and S3-S4) were 1.54% and 0.92%, respectively. The CAT results were in line with the range-average velocity calculated by differential travel time, which RMSE were 0.5163, 0.1494, 0.2103, 0.2804, and 0.2817 m/s for path 1, 2, 3, 4, and 6, respectively. Although the accuracy of results needed to be improved compared to the short-range laboratory experiment, the results from this first application of high-frequency 50 kHz CAT system were promising and verified that CAT might be a solution for small-scale observation. Additionally, we provided a reference for the signal design idea in small-scale CAT observation.
More emphasis is essential on the station layout, position drifting, arrival peak identification, the presence or absence of simultaneous transmission among other aspects to improve the accuracy of small-scale CAT observation. Moreover, further research should place much premium on the mapping of 3-D velocity fields based on the horizontal and vertical slices analyses.

Author Contributions

The work described in this article is the collaborative development of all authors. H.H. contributed to the idea of experiments setup and designed the methods. Y.G. carried out the experiments. X.X. and Y.G. made contributions to data collection, processing and analysis. X.X. and H.W. participated in the drafting and revising of the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Wenhai Program of QNLM under Grant number 2017WHZZB0302.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Timothy, F.; Lin, Y.; Newhall, A.; Helfrich, K.; Lynch, J.; Zhang, W.; Lermusiaux, P.; Wilkin, J. Multiscale multiphysics data-informed modeling for three-dimensional ocean acoustic simulation and prediction. J. Acoust. Soc. Am. 2019, 146, 1996–2015. [Google Scholar]
  2. Kaneko, A.; Zhu, X.; Lin, J. Coastal Acoustic Tomography; Academic Press: Cambridge, MA, USA; Elsevier: Amsterdam, The Netherlands, 2020. [Google Scholar]
  3. The Ocean Tomography Group. A demonstration of ocean acoustic tomography. Nature 1982, 299, 121–125. [Google Scholar] [CrossRef]
  4. Munk, W.; Worceser, P.F.; Wunsch, C. Ocean Acoustic Tomography; Cambridge University Press: New York, NY, USA, 1995; p. 433. [Google Scholar]
  5. Zheng, H.; Yamaoka, H.; Gohda, N.; Noguchi, H.; Kaneko, A. Design of the acoustic tomography system for velocity measurement with an application to the coastal sea. J. Acoust. Soc. Jpn. 1998, 19, 199–210. [Google Scholar] [CrossRef] [Green Version]
  6. Park, J.-H.; Kaneko, A. Computer simulation of coastal acoustic tomography by a two-dimensional vortex model. J. Oceanogr. 2001, 57, 593–602. [Google Scholar] [CrossRef]
  7. Yamaguchi, K.; Lin, J.; Kaneko, A.; Yamamoto, T.; Gohda, N.; Nguyen, H.-Q.; Zheng, H.A. Continuous Mapping of Tidal Current Structures in the Kanmon Strait. J. Ocean. 2005, 61, 283–294. [Google Scholar] [CrossRef]
  8. Yamaoka, H.; Kaneko, A.; Park, J.-H.; Zheng, H.; Gohda, N.; Takano, T.; Zhu, X.-H.; Takasugi, Y. Coastal acoustic tomography system and its field application. IEEE J. Ocean. Eng. 2002, 27, 283–295. [Google Scholar] [CrossRef]
  9. Kaneko, A.; Yamaguchi, K.; Yamamoto, T.; Gohda, N.; Zheng, H.; Fadli, S.; Lin, J.; Nguyen, H.-Q.; Matsuyama, M.; Hachiya, H.; et al. A coastal acoustic tomography experiment in Tokyo Bay. Acta Oceanol. Sin. 2005, 24, 86–94. [Google Scholar]
  10. Chen, M.; Zhu, Z.; Zhang, C.; Zhu, X.; Wang, M.; Fan, X.; Zhao, R.; Lin, J.; Kaneko, A. Mapping Current Fields in a Bay Using a Coast-Fitting Tomographic Inversion. Sensors 2020, 20, 558. [Google Scholar] [CrossRef] [Green Version]
  11. Zhang, C.; Zhu, X.; Zhu, Z.; Liu, W.; Zhang, Z.; Fan, X.; Zhang, R.; Dong, M.; Wang, M. High-precision measurement of tidal current structures using coastal acoustic tomography. Estuar. Coast. Shelf Sci. 2017, 93, 12–24. [Google Scholar] [CrossRef]
  12. Zhu, X.; Zhu, Z.; Guo, X.; Ma, Y.; Fan, X.; Dong, M.; Zhang, C. Measurement of tidal and residual currents and volume transport through the Qiongzhou Strait using coastal acoustic tomography. Cont. Shelf Res. 2015, 108, 65–75. [Google Scholar] [CrossRef]
  13. Zhang, C.; Kaneko, A.; Zhu, X.H.; Howe, B.M.; Gohda, N. Acoustic measurement of the net transport through the Seto Inland Sea. Acoust. Sci. Technol. 2016, 37, 10–20. [Google Scholar] [CrossRef]
  14. Lin, J.; Kaneko, A.; Gohda, N.; Yamaguchi, K. Accurate imaging and prediction of Kanmon Strait tidal current structures by the coastal acoustic tomography data. Geophys. Res. Lett. 2005, 32, 14607. [Google Scholar] [CrossRef]
  15. Syamsudin, F.; Taniguchi, N.; Zhang, C.; Hanifa, A.; Li, G.; Chen, M.; Mutsuda, H.; Zhu, Z.; Zhu, X.; Nagai, T.; et al. Observing Internal Solitary Waves in the Lombok Strait by Coastal Acoustic Tomography. Geophys. Res. Lett. 2019, 46, 10475–10483. [Google Scholar] [CrossRef] [Green Version]
  16. Syamsudin, F.; Chen, M.; Kaneko, A.; Adityawarman, Y.; Zheng, H.; Mutsuda, H.; Hanifa, A.D.; Zhang, C.Z.; Auger, G.; Wells, J.C.; et al. Profiling measurement of internal tides in Bali Strait by reciprocal sound transmission. Acoust. Sci. Technol. 2017, 38, 246–253. [Google Scholar] [CrossRef]
  17. Zhu, X.H.; Zhang, C.; Wu, Q.; Kaneko, A.; Fan, X.; Li, B. Measuring discharge in a river with tidal bores by use of the coastal acoustic tomography system. Estuarine. Coast. Shelf Sci. 2012, 104, 54–65. [Google Scholar] [CrossRef]
  18. Zhang, C.; Kaneko, A.; Zhu, X.H.; Gohda, N. Tomographic mapping of a coastal upwelling and the associated diurnal internal tides in Hiroshima Bay, Japan. J. Geophys. Res. Ocean. 2015, 120, 4288–4305. [Google Scholar] [CrossRef]
  19. Chen, M.; Kaneko, A.; Lin, J.; Zhang, C. Mapping of a Typhoon-Driven Coastal Upwelling by Assimilating Coastal Acoustic Tomography Data. J. Geophys. Res. Ocean. 2017, 122, 7822–7837. [Google Scholar] [CrossRef]
  20. Zhu, X.H.; Kaneko, A.; Wu, Q.; Zhang, C.; Taniguchi, N.; Gohda, N. Mapping tidal current structures in Zhitouyang Bay, China, using coastal acoustic tomography. IEEE J. Ocean. Eng. 2013, 38, 285–296. [Google Scholar] [CrossRef]
  21. Chen, M.; Fadli, S.; Kaneko, A.; Gohda, N.; Howe, B.M.; Mutsuda, H.; Aruni, D.; Zheng, H.; Huang, C.-F.; Taniguchi, N.; et al. Realtime offshore coastal acoustic tomography enabled with mirror transpond functionality. IEEE J. Ocean. Eng. 2020, 45, 645–655. [Google Scholar] [CrossRef]
  22. Zhu, Z.N.; Zhu, X.H.; Guo, X.; Fan, X.; Zhang, C. Assimilation of coastal acoustic tomography data using an unstructured triangular grid ocean model for water with complex coastlines and islands. J. Geophys. Res. Ocean. 2017, 122, 7013–7030. [Google Scholar]
  23. Taniguchi, N.; Li, Y.W.; Huang, C.F. Moving ship tomography experiments in a coastal sea. J. Acoust. Soc. Am. 2016, 140, 3075. [Google Scholar] [CrossRef]
  24. Hanifa, A.D.; Syamsudin, F.; Zhang, C.Z.; Mutsuda, H.; Chen, M.; Zhu, X.H.; Kaneko, A.; Taniguchi, N.; Li, G.M.; Zhu, Z.N.; et al. Tomographic measurement of tidal current and associated 3-h oscillation in Bali Strait. Estuar. Coast. Shelf Sci. 2020, 236, 106655. [Google Scholar] [CrossRef]
  25. Huang, C.F.; Li, Y.W.; Taniguchi, N. Mapping of ocean currents in shallow water using moving ship acoustic tomography. J. Acoust. Soc. Am. 2019, 145, 858–868. [Google Scholar] [CrossRef] [PubMed]
  26. Razaz, M.; Kawanisi, K.; Kaneko, A.; Ioan, N. Application of acoustic tomography to reconstruct the horizontal flow velocity field in a shallow river. Water Resour. Res. 2015, 51, 9665–9678. [Google Scholar] [CrossRef] [Green Version]
  27. Kawanisi, K.; Razaz, M.; Kaneko, A.; Watanabe, S. Long-term measurement of stream flow and salinity in a tidal river by the use of the fluvial acoustic tomography system. J. Hydrol. 2010, 380, 74–81. [Google Scholar] [CrossRef]
  28. Jiang, X.; Gao, Z.; Shang, W.; Liu, C.; Sun, Z. Implementation of data management and analysis system for marine ranching. Remote Sensing and Modeling of Ecosystems for Sustainability XIII. In Proceedings of the SPIE Optical Engineering + Applications, San Diego, CA, USA, 28 August–1 September 2016. [Google Scholar]
  29. Zhao, Z.; Zhang, Y.; Yang, W.; Chen, D. High Frequency Ocean Acoustic Tomography Observation at Coastal Estuary Areas. Adv. Ocean Acoust. 2012, 1495, 360–367. [Google Scholar]
  30. Huang, H.; Guo, Y.; Wang, Z.; Shen, Y.; Wei, J. Water Temperature Observation by Coastal Acoustic Tomography in Artificial Upwelling Area. Sensors 2019, 19, 2655. [Google Scholar] [CrossRef] [Green Version]
  31. Luo, J.; Karjadi, E.; Badiey, M. High-frequency broadband acoustic current tomography in shallow water. J. Acoust. Soc. Am. 2008, 123, 3914. [Google Scholar] [CrossRef]
  32. Li, G.; Ingram, D.; Kaneko, A.; Chen, M. Vertical underwater acoustic tomography in an experimental basin. J. Acoust. Soc. Am. 2017, 141, 3656. [Google Scholar] [CrossRef]
  33. Thorne, R.E. Review: Experiences with shallow water acoustics. Fish. Res. 1998, 35, 137–141. [Google Scholar] [CrossRef]
  34. Guo, Y.; Huang, H.; Wu, J.; Bo Xu Zhu, C.; Wang, Z.; Feng, R. Short-range velocity field observation by coastal acoustic tomography. In Proceedings of the Oceans 2019 mts/IEEE Seattle IEEE, Seatlle, WA, USA, 27–31 October 2019. [Google Scholar]
Figure 1. Sound transmission.
Figure 1. Sound transmission.
Sensors 20 05717 g001
Figure 2. Map of the Panzhinan waterway near the Zhairuoshan Island with adjacent regions (left panel) and the CAT station array at a magnified scale (right panel).
Figure 2. Map of the Panzhinan waterway near the Zhairuoshan Island with adjacent regions (left panel) and the CAT station array at a magnified scale (right panel).
Sensors 20 05717 g002
Figure 3. Experimental settings between the two stations; the acoustic transceivers were suspended down to 1–2 m under the surface with an additional weight by a fixed mount, coastal acoustic tomography (CAT) systems were synchronized with the GPS.
Figure 3. Experimental settings between the two stations; the acoustic transceivers were suspended down to 1–2 m under the surface with an additional weight by a fixed mount, coastal acoustic tomography (CAT) systems were synchronized with the GPS.
Sensors 20 05717 g003
Figure 4. (a) Ray patterns of S1-S3 in the range-independent ray simulation; (b) Sound–speed profile data measured using conductivity–temperature–depth (CTD); (c) Temperature profile data measured using CTD.
Figure 4. (a) Ray patterns of S1-S3 in the range-independent ray simulation; (b) Sound–speed profile data measured using conductivity–temperature–depth (CTD); (c) Temperature profile data measured using CTD.
Sensors 20 05717 g004
Figure 5. Stack diagram of the correlation patterns with first-arrival peaks marked with red circles.
Figure 5. Stack diagram of the correlation patterns with first-arrival peaks marked with red circles.
Sensors 20 05717 g005
Figure 6. Time plots of the differential travel time for five station pairs.
Figure 6. Time plots of the differential travel time for five station pairs.
Sensors 20 05717 g006
Figure 7. Horizontal velocity fields mapping.
Figure 7. Horizontal velocity fields mapping.
Sensors 20 05717 g007
Figure 8. ADCP velocity data with vector plots along the tracks. (a) five ADCP ship tracks. (b) Vector plots for ADCP velocity.
Figure 8. ADCP velocity data with vector plots along the tracks. (a) five ADCP ship tracks. (b) Vector plots for ADCP velocity.
Sensors 20 05717 g008
Figure 9. (a,c) are the ADCP measurement data along S1–S3 and S3–S4, respectively; (b,d) are the comparison of the ADCP average data and the range-average velocity Vm. The solid blue lines represent the Vm and the red dots represent the ADCP average value calculated from ADCP measurement data, respectively.
Figure 9. (a,c) are the ADCP measurement data along S1–S3 and S3–S4, respectively; (b,d) are the comparison of the ADCP average data and the range-average velocity Vm. The solid blue lines represent the Vm and the red dots represent the ADCP average value calculated from ADCP measurement data, respectively.
Sensors 20 05717 g009
Figure 10. Comparison between the CAT inversion and the range-average velocity Vm at nine moments.
Figure 10. Comparison between the CAT inversion and the range-average velocity Vm at nine moments.
Sensors 20 05717 g010
Table 1. Station-to-Station distances.
Table 1. Station-to-Station distances.
Path123456
StationS1-S2S1-S3S1-S4S2-S3S2-S4S3-S4
Distance156.0 m276.6 m239.7 m300.6 m85.1 m328.0 m
Table 2. Sound Transmission based on the round-robin method.
Table 2. Sound Transmission based on the round-robin method.
StationS1S2S3S4
Time
1 minSend
2 min Send
3 min Send
4 min Send
5 minSend
Table 3. Acoustic Doppler current profiler (ADCP) path and its corresponding time.
Table 3. Acoustic Doppler current profiler (ADCP) path and its corresponding time.
Time1st2nd3rd4th5th
Start12:1312:2612:3812:4612:56
End12:2412:3712:4412:5513:01
Table 4. Cross-section volume Transport.
Table 4. Cross-section volume Transport.
PathS1-S2S1-S3S2-S4S3-S4Net Volume TransportProportion of Total Inflow
Time
12:21−273.57−185.87−0.83439.60−20.674.70%
12:32−277.58−247.8252.87481.238.71.63%
12:39−86.56−200.7165.04238.7016.465.42%
12:49−127.42−67.17−20.81200.21−15.17.54%
12:54−380.40−243.31−12.49609.71−26.494.34%
13:04−406.35−252.118.92611.02−38.526.21%
(Unit: m3/s).

Share and Cite

MDPI and ACS Style

Huang, H.; Xie, X.; Guo, Y.; Wang, H. Mapping Small-Scale Horizontal Velocity Field in Panzhinan Waterway by Coastal Acoustic Tomography. Sensors 2020, 20, 5717. https://doi.org/10.3390/s20195717

AMA Style

Huang H, Xie X, Guo Y, Wang H. Mapping Small-Scale Horizontal Velocity Field in Panzhinan Waterway by Coastal Acoustic Tomography. Sensors. 2020; 20(19):5717. https://doi.org/10.3390/s20195717

Chicago/Turabian Style

Huang, Haocai, Xinyi Xie, Yong Guo, and Hangzhou Wang. 2020. "Mapping Small-Scale Horizontal Velocity Field in Panzhinan Waterway by Coastal Acoustic Tomography" Sensors 20, no. 19: 5717. https://doi.org/10.3390/s20195717

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