# Seabed Identification and Characterization Using Sonar

Authors: Henry M. Manik
Journal: Advances in Acoustics and Vibration
 (2012)

---

## Abstract

Application of sonar technologies to bottom acoustics study has made significant advances over recent decades. The sonar systems evolved from the simple analog single-beam and single-frequency systems to more sophisticated digital ones. In this paper, a quantified sonar system was applied to detect and quantify the bottom echoes. The increasing of mean diameter is accompanied by a higher backscattering strength. From this study, identification and characterization using sonar is possible.


---

## Body

## 1. Introduction

Sonar technologies are most effective and useful for sea-bottom exploration. They are based mainly on the measurement, process, analysis, and interpretation of the characteristics of signal reflected or scattered by the sea bottom. Sonar is also increasingly regarded as the remote-sensing tool that will provide the basis for identification, classifying, and mapping ocean resources.There are extensive literatures on the acoustic scattering from the sea bottom [1, 2]. The focus has been on low-frequency features in application such as subbottom classification [3]. Another feature of the sea-bottom scattering has been experimentally observed at a high frequency where the transmitter and receiver are not colocated [4]. This method received contributions both from the bottom surface and subbottom echoes. Most of the data were at grazing angles between 5° and 60°, but some data were collected for the interval between 1° and normal incidence (90°). They obtained results similar to those of Urick [1].One of the acoustic methods to obtain bottom scattering is to use a quantified sonar system (QSS). The QSS can measure echoes generated by reflection and scattering of sounding pulses from the bottom. The observed echo is primarily due to scattering from the water-bottom interface.


## 2. Method

### 2.1. Sonar Equation for Bottom Scattering

The bottom projection is illustrated in Figure1. The elemental backscattered power registered by the transducer is given by
(1)dPRB2=P02r-4exp(-4αr)D4SSdS,
where dPRB is elemental backscattered pressure signal from a sea bottom, P0 is source pressure level, r is range, α is absorption coefficient, D is directivity functions, and SS is bottom scattering. The elemental area dS is located at incidence angle θ, azimuthal angle Ψ, and range r, such that
(2)dS=r2tanθdθdϕ.
The echo pressure amplitude of sea bottom is obtained by integration of (1):
(3)PRB2=P0r-2exp(-4αr)ΦSS,
where Φ is equivalent beam angle for surface scattering
(4)Φ=∫02π∫θ1θ2D4tanθdθdϕ.
The length of pulse in sea water is cτ, and its leading and trailing edges make angles θ1 and θ2 as presented in Table 1.Table 1
Integration limitsθ1 and θ2 for two cases.
Scattering plane
θ
1
θ
2
Circular plane(R≤r<R+cτ/2)
0
cos
-
1
=
(
R
r
)
Circular ring(r≥R+cτ/2)
cos
-
1
=
(
R
t
-
c
τ
/
2
)
cos
-
1
=
(
R
r
)Figure 1
Principle of bottom surface scattering.The signal is amplified to give(5)ERB=PRBMGR,
where ERB is echo amplitude at preamplifier output (V), M is receiving sensitivity of transducer (V/μPa), and GR is preamplifier gain (numeric). Combining (3) and (5) we obtain
(6)ERB2=KTR2r-2exp(-4αr)ΦSS,
where KTR=PoMGR. KTR is transmitting and receiving factor and therefore
(7)SS=ERB2KTR2r-2exp(-4αr)Φ.
In decibel unit SS=10logSS. Simplified block diagram of QSS is shown in Figure 2.Figure 2
Simplified block diagram of quantified sonar system (QSS).


### 2.2. Quantified Sonar System

Quantified sonar system used in this research was PCFF80 model manufacturer by CruzPro, Ltd. (Figure3). The PCFF80 is a full-featured dual frequency (50 and 200 kHz), high-resolution personal-computer-based color fish finder that runs under windows 98, NT, 2000, XP, Vista and Win7 in both analog and DSP mode (digital signal processing). Communications Interface between transducer and PC was conducted using RS-232 serial data.Figure 3
Quantified sonar system.For data acquisition, QSS installed on the research vessel. Echo voltages were recorded on hard disc drive. The QSS was operated at a ping rate of about 40 per minute with a pulse duration of 0.4 ms and beam width of 8.5°.Calibration of QSS is a fundamental component for ensuring high-quality acoustical data. For this purpose, the acoustic system was calibrated with a 38.1 mm diameter of tungsten carbide sphere. The sphere was suspended under the boat at 0.5 m depth to obtain the transmitting and receiving factor (KTR). The target strength of the sphere at the given frequencies were calculated following Miyanohanna et al. [5] and Aoyama et al. [6].By definition, the target strength is given by the ratio between reflected sound intensity,Io, from a target and the sound intensity transmitted towards the target, Ii, referred to 1 m distance. This is regarded as identical to the ratio between the backscattering cross section, σ, for the target and the surface of a sphere with a 1 m radius [7]. The target strength can be expressed on a decibel form in the following way:
(8)TS=10log(IoIi)=10log(σ4π)(dB).
For simplicity, target strength of sphere was computed using
(9)TS=10log(a24)(dB),
where a is radius of sphere.The equivalent beam angleΨ is the solid angle which is measured in steradian at the apex of ideal conical beam which produce same echo integral. Ψ is defined mathematically as
(10)Ψ=∫0π∫ϕ=02πb4(θ,ϕ)sin(θ)dθdϕ.
In logarithmic unit, equivalent beam angle defined as EBA=10log(Ψ) which is expressed in dB relative to 1 steradian.The specifications of the QSS and calibration results are presented in Table2.Table 2
Specification of QSS and calibration data.
Parameters
Quantity
Frequency (kHz)
200
Beam width (deg)
8.5
Equivalent beam angle (dB re 1 sr)
−19.0
Band width (kHz)
4
Pulse duration (ms)
0.4
Target strength of standard sphere (dB)
−39.1
Absorption coefficient (dB/km)
45.5
Transmitting and receiving factor (dB)
50.5


### 2.3. Survey Area

An acoustic survey was conducted in conjunction with oceanographic, fisheries biology, and exploratory fishing in the Seribu Island, North Jakarta Indonesia (Figure4).Figure 4
Research location and bottom sampling point.


### 2.4. Bottom Sample Collection

Collection of bottom samples was accomplished with a system consisting of a sediment sampler (Figure5). The bottom sampler was lowered to the bottom surfaces by using a diver and entrapped the sediment. The total time of operation for one collection was about 30 minutes.Figure 5
Collection of sediment sample.Bottom samples were separated into size component using sieve separation and pipette settling procedures. Bottom material characterization was based on analysis of particles size distributions conducted during the research.


## 2.1. Sonar Equation for Bottom Scattering

The bottom projection is illustrated in Figure1. The elemental backscattered power registered by the transducer is given by
(1)dPRB2=P02r-4exp(-4αr)D4SSdS,
where dPRB is elemental backscattered pressure signal from a sea bottom, P0 is source pressure level, r is range, α is absorption coefficient, D is directivity functions, and SS is bottom scattering. The elemental area dS is located at incidence angle θ, azimuthal angle Ψ, and range r, such that
(2)dS=r2tanθdθdϕ.
The echo pressure amplitude of sea bottom is obtained by integration of (1):
(3)PRB2=P0r-2exp(-4αr)ΦSS,
where Φ is equivalent beam angle for surface scattering
(4)Φ=∫02π∫θ1θ2D4tanθdθdϕ.
The length of pulse in sea water is cτ, and its leading and trailing edges make angles θ1 and θ2 as presented in Table 1.Table 1
Integration limitsθ1 and θ2 for two cases.
Scattering plane
θ
1
θ
2
Circular plane(R≤r<R+cτ/2)
0
cos
-
1
=
(
R
r
)
Circular ring(r≥R+cτ/2)
cos
-
1
=
(
R
t
-
c
τ
/
2
)
cos
-
1
=
(
R
r
)Figure 1
Principle of bottom surface scattering.The signal is amplified to give(5)ERB=PRBMGR,
where ERB is echo amplitude at preamplifier output (V), M is receiving sensitivity of transducer (V/μPa), and GR is preamplifier gain (numeric). Combining (3) and (5) we obtain
(6)ERB2=KTR2r-2exp(-4αr)ΦSS,
where KTR=PoMGR. KTR is transmitting and receiving factor and therefore
(7)SS=ERB2KTR2r-2exp(-4αr)Φ.
In decibel unit SS=10logSS. Simplified block diagram of QSS is shown in Figure 2.Figure 2
Simplified block diagram of quantified sonar system (QSS).


## 2.2. Quantified Sonar System

Quantified sonar system used in this research was PCFF80 model manufacturer by CruzPro, Ltd. (Figure3). The PCFF80 is a full-featured dual frequency (50 and 200 kHz), high-resolution personal-computer-based color fish finder that runs under windows 98, NT, 2000, XP, Vista and Win7 in both analog and DSP mode (digital signal processing). Communications Interface between transducer and PC was conducted using RS-232 serial data.Figure 3
Quantified sonar system.For data acquisition, QSS installed on the research vessel. Echo voltages were recorded on hard disc drive. The QSS was operated at a ping rate of about 40 per minute with a pulse duration of 0.4 ms and beam width of 8.5°.Calibration of QSS is a fundamental component for ensuring high-quality acoustical data. For this purpose, the acoustic system was calibrated with a 38.1 mm diameter of tungsten carbide sphere. The sphere was suspended under the boat at 0.5 m depth to obtain the transmitting and receiving factor (KTR). The target strength of the sphere at the given frequencies were calculated following Miyanohanna et al. [5] and Aoyama et al. [6].By definition, the target strength is given by the ratio between reflected sound intensity,Io, from a target and the sound intensity transmitted towards the target, Ii, referred to 1 m distance. This is regarded as identical to the ratio between the backscattering cross section, σ, for the target and the surface of a sphere with a 1 m radius [7]. The target strength can be expressed on a decibel form in the following way:
(8)TS=10log(IoIi)=10log(σ4π)(dB).
For simplicity, target strength of sphere was computed using
(9)TS=10log(a24)(dB),
where a is radius of sphere.The equivalent beam angleΨ is the solid angle which is measured in steradian at the apex of ideal conical beam which produce same echo integral. Ψ is defined mathematically as
(10)Ψ=∫0π∫ϕ=02πb4(θ,ϕ)sin(θ)dθdϕ.
In logarithmic unit, equivalent beam angle defined as EBA=10log(Ψ) which is expressed in dB relative to 1 steradian.The specifications of the QSS and calibration results are presented in Table2.Table 2
Specification of QSS and calibration data.
Parameters
Quantity
Frequency (kHz)
200
Beam width (deg)
8.5
Equivalent beam angle (dB re 1 sr)
−19.0
Band width (kHz)
4
Pulse duration (ms)
0.4
Target strength of standard sphere (dB)
−39.1
Absorption coefficient (dB/km)
45.5
Transmitting and receiving factor (dB)
50.5


## 2.3. Survey Area

An acoustic survey was conducted in conjunction with oceanographic, fisheries biology, and exploratory fishing in the Seribu Island, North Jakarta Indonesia (Figure4).Figure 4
Research location and bottom sampling point.


## 2.4. Bottom Sample Collection

Collection of bottom samples was accomplished with a system consisting of a sediment sampler (Figure5). The bottom sampler was lowered to the bottom surfaces by using a diver and entrapped the sediment. The total time of operation for one collection was about 30 minutes.Figure 5
Collection of sediment sample.Bottom samples were separated into size component using sieve separation and pipette settling procedures. Bottom material characterization was based on analysis of particles size distributions conducted during the research.


## 3. Experimental Results

The bottom materials of sand, silt, and clay were determined using observed physical characteristics of the samples and mean diameter were calculated (Table3). Bottom images from sand and clay and echogram for each bottom type were presented in Figures 6 and 7. Figure 8 shows the bottom backscattering value for three bottom types. Figure 9 shows the example echo for sand bottom.Table 3
Classification of bottom type by particle diameter.
Bottom type
Sample point
Mean diameter
Sand (%)
Silt (%)
Clay (%)
Sand
1
289
95
2
3
2
305
96
3
1
3
292
94
3
3
Silt
4
45
7
90
3
5
52
4
92
4
6
49
3
93
4
Clay
7
10
3
90
7
8
9
5
91
4
9
11
4
90
6Underwater photography of sand (a) and clay (b).
(a)
(b)Echogram of sand (a), and clay (b).
(a)
(b)Figure 8
Bottom backscattering of sand (⚪), silt (□), and clay (△).Figure 9
Sand-bottom echo for one ping transmission.


## 4. Discussion

The quantified sonar system is useful to measure bottom backscattering (SS). We had derived bottom volume backscattering strength (SV) from SS bottom. In this study area, the increasing of mean diameter is accompanied by a higher backscattering strength. The SS of sand is higher than silt and clay by more than 10 dB. To some extent, it was possible to relate SS to mean diameter suggesting the possibility of bottom-type classification and characterization. Character of the seabed (sediment type, grain-size distribution, porosity, sediment density, sediment velocity, roughness, etc.) are embedded in the sonar echoes from the seabed.The main reason for the higher backscattering strength with larger particle size is that the porosity of sand sediment decreases as the grain size increases. As the porosity decreases, the density increases (less pore water, more mineral constituent). As the density increases, the sediment impedance increases, thus allowing more scattering from a higher impedance contrast between the overlying water and the sediment. Physically, silt and clay have a higher porosity than sand. Acoustic-bottom interaction is too complex to describe by only frequency and mean diameter. The bottom relief also determines the acoustic echo from the seabed. Because sound may penetrate into the sediments and the subbottom, the echoes can also contain information about the zone below the water-sediment interface. Increasingly, sonar technologies are being used in the future to detect, identify, characterize, and classify the sea bottom.



---

*Source: 532458-2012-11-05.xml*