The waveguide invariant summarizes the pattern of constructive and destructive interference between acoustic modes propagating in the ocean waveguide. For many sonar signal-processing schemes, it is essential to know the correct numerical value for the waveguide invariant. While conventional beamforming can estimate the ratio between the waveguide invariant and the range to the source, it cannot unambiguously separate the two terms. In the present work, striation-based beamforming is developed. It is shown that the striation-based beamformer can be used to produce an estimate for the waveguide invariant that is independent of the range. Simulation results are presented.

Consider an underwater passive sonar scenario where a horizontal array measures the acoustic field produced by a distant source. Mapped as a function of frequency ω and range r, the observed acoustic intensity I often exhibits striations. These striations, alternating bands of high and low intensity, result from constructive and destructive interference between the propagating acoustic modes supported by the ocean waveguide. Chuprov1 developed a simple formula for the slope dω/dr of the striations:

(1)

where β is the so-called waveguide invariant. Typically, one would know the frequency and measure the slope allowing the β/r ratio to be calculated.

While the ratio between the waveguide invariant and the range may be readily determined, isolating the range from the numerical value for β is more problematic. In practice, one might simply assume that β is known and then calculate r from the β/r ratio. The canonical value in shallow water is β = 1, but this choice may prove far too coarse. If both the source and array are below the thermocline, both simulations2 and experiments3,4 have shown that β = 1.3 or 1.4 is more appropriate. Erroneously assuming β = 1 when it is actually β = 1.3 translates directly into a 30% error in the estimated value for the range.

In the present work, a method is proposed for isolating the waveguide invariant and producing an estimate for its numerical value that is independent of the range. The method relies on twice beamforming the acoustic pressure measured on a horizontal array due to a distant source. The pressure field is first beamformed conventionally at a single frequency to achieve a focus at the correct bearing of the source. The field is then beamformed again but with each point on the array evaluated at a slightly different frequency. The frequency offset is selected based on the observed slope of the intensity striations. It is shown that the focus will shift from what is observed with conventional beamforming and that the extent of the shift gives a measure of β that is independent of r. The method also gives an indirect coarse estimate of the source’s depth.

To begin the development, assume a point source at r0=(r0cosφ0,r0sinφ0,z0) in the ocean waveguide. The resulting field is measured on a horizontal line array of length L oriented along the y-axis, -L/2<y<L/2, at depth za. Assuming a range-independent environment and an exp(-iωt) convention, the measured complex field is expanded in normal modes as

(2)

where each mode Ψm has corresponding horizontal wavenumber ξm and certain unimportant scaling factors are neglected. It is assumed that the source is coherent over the narrow frequency band of interest. The range r to a point on the array is

(3)

to the term second order in y that expresses wavefront curvature.

Since I=|p|2, the slope of the observed intensity follows from Eqs. (1)–(3):

(4)

where the distant source is assumed not too near broadside, φ0=0, so the wavefront curvature term can be neglected. Solving the separable differential equation yields

(5)

that may be approximated for frequencies near the center frequency ω0 as

(6)

Equation (6) says that the level curves of constant intensity as observed along the array will map onto straight lines of slope S. Consequently, incoherent processing using only the intensity is sufficient to estimate S. The numerical value for S might be calculated by spectral methods, as originally suggested by Chuprov,1 or perhaps by using Radon transforms.5 Regardless of the specific method used to do the calculation, the remainder of this report assumes that the value for S is known.

Having first processed the measured field incoherently to deduce S, it is now processed coherently by beamforming. Define the beamformer output in look direction φ as

(7)

For a conventional Bartlett beamformer, the weights are

(8)

where c0 is the reference sound speed. Substituting Eqs. (2), (7), and (8) and rearranging terms yields

(9)

The terms slowly varying in y and ω are lumped into Cm. The rapidly varying phase terms are included in Dm:

(10)

where the frequency dependence of the horizontal wavenumber ξm has been made explicit.

In conventional beamforming, the frequency is fixed at some constant value, say, ω=ω0. In the present work, a striation-based beamformer is considered where different points along the array are evaluated at different frequencies. The frequency varies linearly across the array based on the observed slope S of the intensity striations. To maintain generality, let the frequency across the array be

(11)

where α is an adjustable parameter. For conventional beamforming, α = 0 and the frequency is constant across the array. For beamforming along the striations, α = 1; see Eq. (6). The frequency-dependent modal wavenumber across the array can be approximated by a Taylor’s series expansion. To first order,

(12)

where vm and um are the phase and group velocities, respectively, for mode m at ω=ω0. The expansion in Eq. (12) should retain the quadratic term to be strictly consistent with the expansion in Eq. (3). The quadratic term in expansions like Eq. (12) is usually neglected; see the discussion in Ref. 6 for the limitation associated with this approximation.

Combining the above equations yields for the striation-based beamformer

(13)

where terms have been grouped in like order of y with coefficients

(14)
(15)
(16)

It is instructive to consider the linear coefficient θ1 in greater detail. In transitioning between the two forms of Eq. (15), the definition of the striation slope S, Eq. (6), has been applied. Since S is proportional to the ratio β/r0, the product Sr0 appearing in Eq. (15) is independent of r0. The key point is that θ1 therefore depends on the waveguide invariant but is independent of range; β has been isolated from r0 by striation-based beamforming.

Equation (13) can be evaluated exactly in terms of Fresnel integrals.7 To gain greater insight into the result, it is useful to assume that the source is in the far field of the receiving array. In the far field, the quadratic term in Eq. (13) is neglected and the integral reduces to

(17)

where sinc(x)sin(x)/x. The magnitude is a maximum when θ1 = 0 implying focusing when the beamformer look direction φ satisfies

(18)

where the approximate form ignores the small differences between c0 and the phase and group velocities to yield a result that is independent of the mode index.

Equations (6) and (18) suggest a three step processing strategy for estimating the waveguide invariant. First, process the field incoherently to estimate the slope S of the intensity striations. Second, beamform conventionally with α = 0. The beamformer will focus at the true bearing of the source, φ=φ0. Finally, perform striation-based beamforming with α0. The focus will shift to some new φφ0 but, since the true φ0 has already been estimated, the location of the new focus can be used in Eq. (18) to solve for β. The value for the waveguide invariant is thereby estimated without knowing the range or depth of the source.

The numerical value for the adjustable parameter α should be selected with some care. The choice α = 1 is attractive because it means a strong signal across the array without fades. This should improve performance in a noisy environment. When the source is at a bearing too far from array broadside, however, the choice α = 1 will fail to produce a focus. For example, no beamformer look direction φ will satisfy Eq. (18) if α = β = 1.0 and |φ0|>30. If conventional beamforming (α = 0) reveals a source far removed from broadside, the striation-based beamformer should use a value for α that ensures that Eq. (18) will have a solution. In addition, α must be selected so that the total frequency shift across the array, αSL, is within the coherence bandwidth of the source.

A specific numerical simulation serves to illustrate the striation-based beamforming concept. Consider an environment typical of the East China Sea as observed during a 2001 experiment.8 The 105 m deep water column has a 40 m deep surface mixed layer above a sharp thermocline. The sound speed contrast between the surface mixed layer and bottom of the water column is 12 m/s. Assume a 128 element horizontal receiving array with 0.5 m spacing (L = 63.5 m) positioned at depth za = 60 m. The source is positioned at range r0 = 5 km and bearing relative to array broadside φ0  = 20°. The initial source depth is z0 = 60 m so that both the source and the receiving array are positioned below the part of the thermocline with a sharp gradient. With both the source and the receiving array below the main thermocline, the low-order acoustic modes that have an upper turning depth below the sea surface will contribute strongly to the acoustic field. In such a scenario, previous results2–4 show that β > 1.0 might be expected.

Figure 1 shows the simulated intensity as calculated along the array over a 50 Hz band. The calculation was made using a normal mode code and the plot has 20 dB dynamic range. Striations are plainly evident. A plausible estimate for the slope of the striations is S = 0.089 Hz/m implying a SL = 5.7 Hz frequency shift over the length of the array. Superimposed on the image are lines corresponding to α = 0 and α = 1 using reference frequency ω0/2π = 1 kHz in Eq. (11).

A discrete Bartlett beamformer was implemented and Fig. 2 shows the results. Beamforming conventionally (α = 0), the figure shows a focus at look direction φ = φ0 = 20° as expected. The plot is normalized so the conventional B(φ0) = 1.0. The field was then beamformed along the striation (α = 1) with the same normalization applied and the result plotted. As expected, the output at the focus is higher for the striation-based beamformer. More importantly, the location of the focus has shifted to φ = 52°. Using φ = 52° and φ0 = 20° in Eq. (18) yields a value for the waveguide invariant β = 1.3, a value consistent with experimental results3,4 when both source and receiver are below the main thermocline.

The simulation was repeated with all the same parameters except the source depth was changed to z0 = 10 m. With the source in the surface mixed layer, the low-order modes that cause β >1.0 are not excited. Striations are again evident, but with a shallower slope than is observed in Fig. 1; S = 0.069 Hz/m is a reasonable value. The conventional and striation-based beamformers focus at φ = φ0 = 20° and φ = 43°, respectively. Using these angles in Eq. (18) yields β = 1.0 for the shallow source.

Because the striation-based beamformer is able to quantify β, it also provides a coarse estimate of the source depth, at least to the extent that it can distinguish between shallow and deep sources. A caveat is that the present analysis has assumed a range-independent environment. Range dependence from such environmental factors as internal waves will cause the acoustic modes to couple7 and make source depth estimation more problematic.

To summarize, while conventional beamforming can determine a source’s bearing and the β/r0 ratio, it is unable to separate the waveguide invariant β from the range r0. Striation-based beamforming breaks the β/r0 ambiguity and provides an independent estimate for the waveguide invariant.

This work was supported by the Office of Naval Research.

1.
S. D.
Chuprov
, “
Interference structure of a sound field in a layered ocean
,” in
Ocean Acoustics: Current State
, edited by
L. M.
Brekhovskikh
and
I. B.
Andreevoi
(
Nauka
,
Moscow
,
1982
), pp.
71
91
.
2.
D.
Rouseff
and
R. C.
Spindel
, “
Modeling the Waveguide Invariant as a Distribution
,” in
Ocean Acoustic Interference Phenomena and Signal Processing
, edited by
W. A.
Kuperman
and
G. L.
D’Spain
(
AIP
,
New York
,
2002
), pp.
137
148
.
3.
W. S.
Hodgkiss
,
H. C.
Song
,
W. A.
Kuperman
,
T.
Akal
,
C.
Ferla
and
D. R.
Jackson
, “
A long-range and variable focus phase-conjugation experiment in shallow water
,”
J. Acoust. Soc. Am.
105
,
1597
1604
(
1999
).
4.
J. E.
Quijano
,
L. M.
Zurk
, and
D.
Rouseff
, “
Demonstration of the invariance principle for monstatic active sonar
,”
J. Acoust. Soc. Am.
123
,
1329
1337
(
2008
).
5.
M.
Born
and
E.
Wolf
,
Principles of Optics
, 7th ed. (
Cambridge University Press
,
Cambridge
,
1999
), pp.
219
227
.
6.
J. A.
Stratton
,
Electromagnetic Theory
(
McGraw-Hill
,
New York
,
1941
), pp.
330
340
.
7.
Handbook of Mathematical Functions
, edited by
M.
Abramowtiz
and
I. A.
Stegun
(
U. S. Govt. Printing Office
,
Washington, D.C.
,
1964
), p.
304
.
8.
D.
Rouseff
and
D.
Tang
, “
Internal wave effects on the ambient noise notch in the East China Sea: model/data comparison
,”
J. Acoust. Soc. Am.
120
,
1284
1294
(
2006
).