MISCELLANEOUS TECHNICAL ARTICLES BY Dr A R COLLINS

Underwater Sound Propagation

Temperature and sound velocity profiles

The speed of sound in the ocean is a function of temperature, pressure and salinity. Eqn (1) shows the dependence of sound speed C, on these factors: $$ C(T, P, S) = 1449.2 + 4.6\ T + 0.055\ T^2 + 1.39\ (S - 35) + 0.016\ D \tag{1} $$

where: \(C\) = sound speed in m/sec,
\(T\) = water temperature in ° Celsius,
\(D\) = depth below the surface in metres,
\(S\) = salinity in parts per thousand.

Typical ocean salinity, \(S\), is fairly constant at 35 ppt and so the sound speed formula reduces to a function of temperature and depth.

From Eqn 1 it can be seen that the sound speed is strongly dependant on the water temperature \(T\). Solar radiation warms the surface of the ocean, wave action tends to mix the warm surface water generating a surface layer of fairly uniform temperature. In temperate and tropical waters this isothermal surface layer may extend from 20 m to 200 m deep. The water temperature below this layer rapidly drops away. Such a region of falling temperature is called a negative thermocline.

The ocean depths, where little or no sunlight penetrates is of uniform temperature, varying only 1 or 2°C over thounsands of meters of depth. In these isothermal regions the sound speed, \(C\), is largly detirmined by the pressure which increases uniformly with depth. The rate of increase is 0.016 m/sec per metre of depth.

Fig 1a shows temperature as a function of water depth, a bathythermograph, for the first 1000m below the surface. This bathy's mixed layer and thermocline are clearly defined and are typical of warmer ocean regions.

Fig1

Figure 1a. Typical Eastern Pacific winter bathythermograph

Figure 1b. Sound velocity profile corresponding to the bathy shown in Fig 1a.

In Fig 1 the isothermal layer extends down about 75m. The temperature in this layer is 17°C.

Fig 1b shows the corresponding sound speed profile. In the isothermal layer the sound speed is about 1512 m/sec while below 100m the sound speed has fallen to below 1500 m/sec due to the decreased temperature. The sound speed then rapidly falls with the lower temperature, until the deeper layers where the pressure dominates and the sound speed increases again.

In higher latitude, colder waters less solar radiation results in a small isothermal surface layer or none at all. A typical cold water profile is shown in Fig 2a. The temperature continuously declines with depth.

The corresponding sound profile, shown in Fig 2b is typical of high latitudes. The marked differences in acoustic conditions in different locations indicate that sonar conditions in one area may not be applicable to another area or season.

Knowledge of the temperature profile and hence the sound speed profile, enables the prediction of sound transmission paths. Such information serves as input to the tactical environment in which submarine and surface navies operate.

Fig1

Figure 2a. Bathy with no isothermal surface layer.

Figure 2b. Sound speed profile corresponding to the bathythermograph of Fig 2a.

Ships employ eXpendable BathyThermographs (XBT) to measure the water temperature as a function of depth. These probes sink at a known rate. The temperature measurements are sent back to the ship's receiver via a fine copper wire unreeling from the probe and are logged or displayed (on an appropriate scale) as sound speed versus depth.

Submarines use a similar device, the SSXBT, which is ejected from the dived submarine, floats to the surface and then sinks at a known rate.

Sound ray tracing

The propagation paths of sound in the ocean may predicted using much the same techniques employed to predict the reflection and refraction of light rays.

In crossing the interface between the different media, such as from air to glass, the light ray is refracted or bent. The angle of refraction depends on the ratio of the speed of light in the two media, this is Snells law. A similar process occurs with sound transmission in the ocean.

For calculating the ray paths, the ocean may be considered as a set of thin layers of water of different sound speeds.

Fig 3 illustrates Snells Law for refraction of sound rays.

Fig3

Figure 3. Diagram of the refraction of sound rays.

Referring to Fig 3, Snells law for sound in water may be stated as:

$$ \frac{\cos(B_1)}{\cos(B_2)} = \frac{C_1}{C_2} \tag{2} $$

Hence, if the angle of a ray to the horizontal is \( B_1 \) in a region of sound speed \( C_1 \), then the refracted ray angle, \( B_2 \) as it passes into region with sound speed \( C_2 \) can be calculated.

From Eqn 2, if \( C_1 > C_2 \) then \( cos(B_1) > cos(B_2) \) and hence \( B_1 < B_2 \), i.e. for a ray passing into a slower sound speed layer, the ray angle is larger, it bends away from the horizontal. Conversely, if \( C_1 < C_2 \) then \( B_1 > B_2 \) and the ray is shallower, nearer to horizontal. When the sound speed increases the ray angle tends to decrease.

Snells law may be restated for acoustic rays: "sound rays passing from one medium to another bend toward the medium with the slower wave speed"

If a bathy has been taken, the sound speed C as a function of temperature and depth is known and the angles through which sound rays will bend as they pass through different depths can be calculated and the ray paths plotted.

Fig 4 shows some ray paths which may occur for the bathy in Fig 1.

Fig4

Figure 4. Possible ray paths for the bathy in Fig 1.

For a sound source located in the isothermal layer near the surface, the sound speed will increase with depth hence there a ray starting downwards will be bent back toward the surface. There will be a ray launch angle which just becomes horizontal at the bottom of the mixed layer. This ray and all smaller angle rays are trapped in the isothermal layer as they don't cross into the negative thermocline. The isothermal layer forms a sonar duct for these rays.

The launch angle of this critical ray is given by putting \( B_1 = 0 \) and hence \( cos(B_1) = 1 \) in Eqn (2), and substituting the appropriate values for the sound speed.

$$B_2 = \cos^{-1}\left(\frac{C_2}{C_1}\right)$$

For a source depth of 5 m, a typical depth for a surface ship active sonar, the source angle for rays to be trapped in the duct is given by:

$$B_s = \cos^{-1}\left(\frac{1512 + 0.017 \times 5}{1512 + 0.017 \times 60}\right) = 2^\circ$$

The same ray angle applies for rays that start in an upwards direction as these are merely reflected at the sea surface and then pass through the source depth at 2° declination. Hence rays that leave the sonar within 2° of horizontal are trapped within the duct (see Fig 6).

It follows that rays which start at steeper angles can penetrate below the layer. A ray that starts with a steeper angle escapes the duct and is immediately refracted downward towards the colder water, slower sound speed region. It follows that there is a region through which rays don't pass, called the SHADOW ZONE. This is the region where a submarine may be relatively invisible to a surface ship active sonar.

Deep ocean propagation

Water temperature deep in the oceans is very nearly constant and in such regions the sound speed is dominated by its dependence on the water pressure which increases with depth. Rays which escape the duct are initially refracted downward toward the colder water which has a lower sound speed. Depending on the depth of water the sound may hit the sea floor and be reflected back up, but if the water is sufficiently deep, the sound penetrates to depths where the sound speed has increased significantly. The rays begin to bend back toward the slower sound speed at shallower depths without hitting the bottom.

For this REFRACTED transmission mode to occur, the deep water sound speed must at least equal the sound speed at the source. The sound rays starting horizontal near the surface will become horizontal again at this depth and then bend upward toward the slower sound speed layers above. Rays starting out at steeper angles require deeper water with even higher sound speeds to become horizontal and then turn back toward the surface.

Once the sound has been refracted upwards it travels to the source depth where it again encounters higher sound speeds and bends back to the horizontal. The same process can then be repeated.

For sound rays transmitted near the surface to loop back to the surface without encountering the bottom the deep water sound speed must equal the surface sound speed. The higher sea surface temperature, the deeper the water must be for the deep sound speed to equal the surface speed. So in temperate zones this type of transmission can only occur in deep ocean basins. For example, a surface temperature of 23°C which indicates a sound speed of about 1530 m/sec and the equivalent water depth required is a little over 4000 m. So the phenomenon can occur at reasonable high water temperatures in the ocean basins which are typically 5000 m deep, shallower water requires cooler surface temperatures for the effect to develop.

Fig 5 shows a bathy east of Hawaii in very deep water.

Fig1

Figure 5a. Bathy east of Hawaii showing a typical January temperature profile in deep water.

Figure 5b. Sound speed profile showing the formation of a low sound speed layer in deep water.

Convergence Zones

When refracted transmission occurs the rays returning to the surface tend to do so at a more or less constant range from the source resulting in a higher level sound level being observed in an annulus about the source. The sound reflecting off the surface then repeats the cycle looping back to the surface at regular distances of about 50 to 70 km, depending on the surface water temperature. The regions of higher sound strength are termed CONVERGENCE ZONES.

Fig 6 shows a plot of convergence zone formation with the bathy from Fig 5.

Fig6

Figure 6. Deep water sound ray paths showing convergence zones.

Bottom Bounce

In a region where the water depth is too shallow for refracted transmission paths, sound rays will intersect the sea floor and be reflected. More energy is reflected if the sound speed in the bottom material is hard with much faster sound speed than water. A rocky bottom will reflect sound more efficiently than a muddy bottom.

Thus reflection depends on the characteristics of the bottom and is generally frequency sensitive, with the result that energy in higher frequency bands is reflected better than energy at lower frequencies. Nevertheless BOTTOM BOUNCE can be utilized to advantage.

Fig 7 shows such an application. Energy from a near surface source, a hull mounted active sonar is initially refracted down and intersects the sea floor, where it is reflected. Depending on the ray angle, this reflected energy may be able to reach a region that could not be ensonified by direct transmission paths. Active sonars can use this mode to ensonify the shadow zone.

Fig7

Figure 7. Ray paths using bottom bounce.

SOFAR propagation

If the sound source is below the layer, then a variation of the refracted propagation paths may be sustained. This is illustrated by Fig 7. Because the sound speed increases both above and below the source depth, rays become trapped in this deep sound channel without suffering a lossy reflection at bottom or surface. This is known as SOFAR propagation. The acronym being derived from Sound Fixing And Ranging but indicating that sound in this low loss transmission mode travels so far.

Fig8

Figure 8. SOFAR propagation for deep launched rays.

Spreading laws and propagation loss

Intensity, \(I\), of sound is a measure of energy flowing across a unit area perpendicular to the direction of propagation. As sound energy travels through the ocean and spreads over a larger volume, the intensity of the sound decreases.

The term TRANSMISSION LOSS describes the drop in sound energy level as it propagates from one point in the ocean to another. The transmission loss is defined as the ratio of the acoustic intensity at the particular point, to the level at some reference position. If the reference is the source level of a transmitter, then the reference position is taken at a distance of 1m from the source (1 yard in the Imperial system). Hence, the transmission or propagation loss at a point P is given by,

$$\frac{I_p}{I_{ref}}$$

where \( I_p \) and \( I_{ref} \) are the respective acoustic intensities at point P and at the reference position.

As sound energy travels from the source it spreads out over the surface of a sphere which increases in radius \(R\). For a given amount of energy emitted the intensity of the sound will be determined by the size of the sphere over which it has spread. At a distance \(R\) from the source the Intensity, \(I\) is inversely proportional to the surface area of sphere. $$I = \frac{k}{R^2} $$ where \(k\) is a constant.

If \(R\) is the distance to the reference point, 1 metre from the source, then the transmission loss for spherical spreading, expressed in decibels (dB) is given by: $$TL = 20 \log (R)$$

The source level of sonar transmitters are quoted as the level 1 m from the transmitter, for a high power anti-submarine sonar this is typically about 220dB re 1µPa. This formula for transmission loss due to spherical spreading implies that every time the range doubles, the sound intensity falls by 6dB. So 2m from the source the level is only 214dB and so on. Doubling the range 10 times gives a sound level 60dB down at 1024m, ie. approximately 1km from the source the level has dropped to 160db. At 16km range, the level would be down to 136dB.

Another factor which affects propagation in the ocean is ABSORPTION. Sound energy is dissipated by a lossy interaction with some of the salts in the water. The process is non-linear with the rate of absorption varying roughly as the square of the frequency, as well as being a function of salinity and temperature. A plot of the absorption of sound in sea water according to temperature and frequency are shown in Fig 8.

The effect of the extra absorption loss is to modify the transmission loss equation so that it now becomes: $$TL = 20 \log (R) + aR $$ where \(a\) is the absorption in dB per metre.

Fig9

Figure 9. Absorption of sound in sea water (20°C).

The absorption of sound energy causes a fall in sound pressure level for every kilometre the sound travels in addition to the fall due to spherical spreading. This additional absorption loss is small to negligible for frequencies less than 1kHz. But for higher frequency sonars the loss is very significant. For a sonar operating at 15 kHz the loss is 2dB per km. So at 16km there will be a 32dB loss in addition to the 84dB spreading loss.