1. Introduction
It is well known that the parameters of natural waveguides (ocean, atmosphere), as well as optical, electrodynamic, plasma, dielectric, and other waveguides, to one degree or another, are subject to random variations [
1,
2]. For example, during the propagation of an acoustic signal in the sea, one of the general sources of fluctuations of the main parameter, that is the speed of sound, is the passage of internal gravity waves (IGW) [
3,
4,
5]. This circumstance is especially significant in the coastal (shelf) zones of the ocean, where a shallow water waveguide with a thermocline in the water column and a penetrable, absorbing bottom occurs for a low-frequency sound signal. As a rule, weak random perturbations of the speed of sound in such a waveguide cause a field of diffuse (background) internal waves that form on the “tails” of numerous collapses of nonlinear internal waves (solitons) when the latter passes into the shallow part of the sea shelf from the deepwater region [
6]. Traditionally, on the basis of the central limit theorem, it is assumed that background internal waves lead to Gaussian perturbations of the sound speed in the sea. So, the influence of such fluctuations of the sound speed on the propagation of an acoustic signal in an oceanic waveguide has been studied for the past 40 years in numerous studies on this subject [
7,
8,
9,
10,
11,
12,
13,
14,
15,
16,
17,
18,
19,
20], both theoretically and experimentally. Modern ideas about the effect of random inhomogeneities of the marine environment on the sound fields are based on the results of an approximate theory described in [
7,
8,
9], as applied to the description of weak perturbations in the sound speed in the deep ocean when the internal waves travel with a special type of spectrum (Garrett-Munk spectrum). Subsequent attempts to transfer this theory, called the diffusion approximation, to sound propagation in the shelf zones of the ocean [
11,
12,
13,
14,
18,
19,
20,
21,
22,
23,
24,
25,
26] showed that the theory poorly describes the situation if there are not fluctuations that are too weak, not small losses in the medium, but actual and also at typical distances for the shallow sea, not exceeding 50–100 km. There is nothing surprising in the fact that for the conditions of sound propagation in the shallow sea, fundamentally different than in the deep ocean, it became necessary to make adjustments to the previous theory in order to find out new statistical patterns. The first steps in this direction were made in [
11,
21], where for weak fluctuations of the sound speed, in particular, strong fluctuations of the sound field were established at relatively small distances from the source and unusual changes were found in the law of decay of the field average intensity. We emphasize that new results are inherent in shallow water areas, where the bottom effect is significant. In the deep ocean, where acoustic waves propagate within the underwater sound channel, such patterns are not registered [
9,
10,
15,
16,
17]. This paper briefly presents the previously obtained results of solving a two-dimensionally inhomogeneous statistical problem of sound propagation in a fluctuating shallow water waveguide with horizontal boundaries [
21,
22,
23,
24,
25,
26]. Further, on the basis of the developed approach, examples of non-Gaussian probability distributions of sound velocity fluctuations are considered, and for these cases the results of a statistical simulation of the average intensity of a low-frequency acoustic signal are presented. A comparative analysis showed that the conclusions obtained for Gaussian fluctuations of the speed of sound are also valid for the considered examples with non-Gaussian fluctuations. Moreover, non-Gaussian fluctuations may even somewhat enhance the previously established effects of attenuation of the decay of the average signal intensity in the waveguide.
2. Formulation of the Statistical Problem and Some Analytical Results
In a cylindrical coordinate system (
r,
z), we consider a marine waveguide with horizontal boundaries consisting of a water layer and a bottom layer of liquid sediments. Suppose that in the water column with a constant density ρ
0, the average (regular) sound velocity
c0(
z) experiences weak two-dimensional fluctuations δ
c(
r,
z) =
c(
r,
z) −
c0(
z), |δ
c/c0| << 1. The layer of liquid sediments is homogeneous with constant values of density, sound velocity, and absorption: ρ
1,
c1, β
1. The random sound velocity field ε(
r,
z) = −2δ
c(
r,
z)/
c0 is characterized by an average zero value 〈ε(
r,
z)〉 = 0, and is described by some anisotropic correlation function
Bε(
r1−
r2,
z1−
z2) ≡ 〈ε(
r1,
z1)ε(
r2,
z2)〉 (hereinafter, the angle brackets mean averaging over the ensemble of realizations). This function has amplitude
Bε(0,0) = σ
ε2 (dimensionless intensity of fluctuations) and is characterized by scales of spatial correlation in the horizontal direction
Lr and along the depth
Lz. The specific form of the correlation function does not limit the generality of the research findings, but is usually dictated by the convenience of analytical evaluations [
23,
24,
25].
The pressure field
p(
r,
z) and components of velocity {
u,w} of a point source with a frequency ω, located within a waveguide at
r = 0,
z =
z0, satisfy linear equations of acoustics with random coefficients:
Variations in the density of the medium ρ involved in the Equations (1) for acoustic frequencies above 1 Hz do not affect the sound propagation [
4,
5], therefore in underwater acoustics it is usually assumed in the water layer that ρ = ρ
0, and in the bottom we will also assume ρ = ρ
1.
The pressure field
p(
r,
z) in the wave zone can be represented by decomposition into local modes of an irregular waveguide (density in the water ρ = ρ
0 = 1 kg/m
3):
Eigenfunctions φ
m at the surface (
z =
H) and at the bottom (
z = 0) of an ocean satisfy the following boundary conditions: φ
m(
r,
H) = 0, φ
m(
r,0) +
gm(
r)φ′
m(
r,0) = 0, where
gm(
r) characterizes the impedance of the penetrable bottom, and the square of the wave number
k2(
r,
z) =
k02(
z)[1+ε(
r,
z)] is a random function due to fluctuations in the speed of sound ε(
r,
z),
k0=ω/
c0(
z). In [
27,
28,
29], it was shown that in an irregular waveguide ignoring the backscattered field (that is, in the one-way propagation or forward scattering approximation), the amplitudes of modes
Gm(
r) satisfy the following quadrature representation:
where κ
m r >> 1,
am = φ
m(0,
z0)/2,
Am =
iam[2π
iκ
m(
r)
r]
−1/2,
m = 1,2…. In this representation
is the element of an skew-symmetric matrix,
Vmn (
r) = −
Vnm (
r),
Vnn = 0, describing mode coupling due to horizontal variations caused by fluctuations in the speed of sound. The neglect of backscattering in the problem under consideration is justified due to the smallness of ε(
r,
z) [
27,
28]. To perform analytical estimates, we use an equation that obeys the local eigenvalues κ
m(
r), as well as an expression for the mode coupling matrix, written via derivatives of the waveguide parameters [
23,
24,
25,
30]:
These relations follow from the original equations and they are valid in the case of horizontal changes in the medium due to fluctuations in ε(
r,
z). For actual observed values of ε(
r,
z), random variations of eigenfunctions and eigenvalues of modes are very small, so in the first order of the small perturbation method with respect to ε we can assume φ
m(
r,
z) ≈ φ
0m(
z), where
is the eigenfunction of the deterministic (unperturbed) problem, when ε(
r,
z) = 0. In this approximation, for eigenvalues we obtain a linear functional dependence on fluctuations:
If we assume
, where κ
0m = κ
m (0), ε(0,
z) = 0, then taking into account (4)
. Due to the smallness of random variations of δκ
m2(
r) the following relations are also valid:
Similarly, from (3) follows an approximate expression for the mode coupling matrix:
For the mode amplitudes
Gm(
r) solution (2), taking into account (5) and (6), is written (
an/
am ≈ 1):
The second term in the exponent (7) contains the integral of
r from the inhomogeneities and, therefore, characterizes the effect of fluctuations of the speed of sound on the acoustic field accumulating as the distance increases. In the third term related to the coupling of modes, there is no such integral. This implies that the coupling of modes is some variable-sign addition to the solution, weakly dependent on distance. This suggests that in this situation the adiabatic approximation to the solution for mode amplitudes
describes the main statistical effects, which are just accumulating with distance. Calculations in many cases confirm this conclusion. We note that in the framework of the theory of [
7,
8,
9,
13,
14,
15,
16,
17,
20] instead of (1′), the solution is sought through the eigenvalues and eigenfunctions of the unperturbed waveguide for ε(
r,
z) = 0) κ
0m, φ
0m. In addition, for the mode amplitudes
Gm(
r) and mode coupling matrix
Vmn (
r) approximate equations are used:
In the first order of the perturbation method with respect to ε, the diagonal term of the matrix
V in (8) gives a solution that is close to adiabatic (the first term of the exponent in (2), based on (4), (5)
. However, as shown in [
31], for actual parameters of a shallow sea waveguide with a not too rigid bottom, the difference between solution (8) and more accurate (2) can be significant and grows with distance.
The second statistical moment of the pressure field (average intensity), represented by incoherent and coherent sums, has the form:
The simplest analysis of the structure of sums in (9) shows that for small fluctuations ε(
r,
z), substantial changes (with respect to the deterministic problem) in both sums are possible only due to accumulating statistical effects in exponential terms
GnGm∗(
r). Due to sound absorption in marine sediments (β
1) and penetration into the bottom, the horizontal wave number of modes is always complex
. We emphasize that it is precisely the sequential consideration of this circumstance, as compared with all known works, that allows analytical interpretation, and in general cases, with the help of numerical modeling, to establish new features of the behavior of average intensity and moments of higher order. It is known that the first, incoherent sum [
21,
22] makes the largest contribution to the mean intensity value. The oscillating coherent sum in (9), which is responsible for the intermode beats in some field realizations, after averaging, does not contribute to the solution when moving away from the source. Consider the asymptotic dependence on
r of the terms of the incoherent series (9), assuming that the fluctuations ε are Gaussian. Averaging over the ensemble of realizations of the function
GnGn∗(
r), with (2′) taken into account, allows us to obtain the following expression [
25]:
where notations are entered:
Let us consider two limiting cases of the horizontal scales of fluctuations:
r << Lr and
r >> Lr. In the first case of small distances with relatively large horizontal scales of fluctuations, it can be approximately assumed that the correlation function of ε depends only on
z:
Bε(
r1−
r2,
z1−
z2)
≈ B1(
z1−
z2). In the second case of large distances, the correlation function
Bε(
r1−
r2,
z1−
z2) is “sharp”, so it can be replaced with an effective
Bε(
r1−
r2,
z1−
z2)
≈ δ(
r1−
r2)
B2(
z1−
z2) [
2,
25]. Here,
B1(
z1−
z2) and
B2(
z1−
z2) are functions determined from the condition of conservation of the normalization of new effective correlation functions:
B1(
z1−
z2) =
Bε(0,
z1−
z2), and the function
B2(
z1−
z2) = C
Bε(0,
z1−
z2), where the constant C is determined from the following identity:
. Thus, replacing it in the first of formulas (11) the initial correlation function with its effective analogues, we obtain the following asymptotics:
To illustrate the formulas (12), we choose the anisotropic correlation function in the form
. For example, a similar function was used in [
14] to describe the effect of internal waves on sound propagation in the Florida Strait. In this case, one can accurately integrates over
r the expression (11) for
Λnn(
r):
Obviously, the asymptotic expressions (12) follow from the formula (13). Thus, on the basis of (12) it is clear that the function Λ
nn(
r) changes from linear, with
r/Lr >> 1, to quadratic, with
r/Lr << 1. So, the curve describing the dependence of the average field intensity 〈|
p|
2〉 on the distance is higher than the curve |
p|
2, corresponding to the solution of the deterministic problem. This rise in average intensity is determined by the exponent in formula (10). At short distances from the source
r Lr, this rise is faster than at long distances. Obviously, in the layered problem (
Lr → ∞), the effect is observed along the entire path of sound propagation [
21,
22]. In this case, as shown above (estimates (12)), the asymptotics considered do not depend on the specific type of correlation function
Bε(
r1−
r2,
z1−
z2) and, accordingly, on the form of the spectral density of inhomogeneities. Therefore, reference to the Garrett-Munk spectrum [
3,
4,
5,
6,
7,
8,
9], adopted when studying the passage of sound through random inhomogeneities in the form of background internal waves, is not mandatory, but from an analysis point of view it is undesirable because of the complex empirical appearance of this spectrum. Unfortunately, it is possible to obtain transparent analytical expressions (10)–(13) only under the assumptions made above about the Gaussian fluctuations and the small contribution of mode coupling to the acoustic field. Thus, in the general case, numerical simulation is of particular interest.
3. Results of Numerical Simulation
Below are examples of statistical simulation for Gaussian and non-Gaussian fluctuations in the speed of sound. As a model for calculations, a shallow water three-layer waveguide with a regular thermocline and two-dimensional fluctuations of the speed of sound (
Figure 1) was taken.
Fluctuations are described by the exponential correlation function
Bε(
r1−
r2, z1−
z2), given above. The waveguide characteristics are as follows: Depth
H = 50 m, surface layer 15 m thick, where
c0(
z) = 〈
c(
z)〉 = 1525 m/s, and intensity of fluctuations 〈(δ
c/
c0)
2〉 = 10
−6, bottom layer 10 m thick, where
c0(
z) = 1500 m/s and 〈(δ
c/
c0)
2〉 = 10
−6, a layer of linear thermocline 25 m thick with 1500 m/s ≤
c0(
z) ≤ 1525 m/s, and 〈(δ
c/
c0)
2〉 = 10
−5. ρ
1,
c1, β
1 were the parameters of a homogeneous liquid bottom were used different for the representativeness of the research. These hydrological conditions are typical (with the exception of constant depth) for the autumn observation period in the shelf areas of the Sea of Japan, when a pronounced thermocline region is formed in the water layer. Active hydrodynamic processes occur at these horizons, randomly disturbing the average sound velocity profile. The considered formulation of the problem corresponds to the passage of a sound signal along the propagation path of internal waves (major mode), or other hydrodynamic disturbances. For calculations, a sound frequency of 500 Hz was chosen as the reference. In this situation, 4–11 propagating (trapped) modes are formed in the waveguide, depending on the bottom penetrability, or only leaky modes are present. In the presence of a thermocline, the behavior of horizontal wave numbers κ
m(
r) (eigenvalues of modes) on the complex plane (
Re(κ
m),
Im(κ
m)), is characterized by substantial non-monotonicity, unlike a homogeneous water layer (Pekeris waveguide model), for which the consistent growth of κ
m on the complex plane is typical. The presence of fluctuations enhances non-monotony: in many individual realizations, the modes with higher numbers (2nd–5th) can propagate even better than the 1st one. This feature leads to the rapid development of strong fluctuations of the field in the waveguide, established in [
11,
14,
22,
32], which is manifested in a significant dispersion of the levels of individual intensity realizations. Average sound field intensity 〈
I〉 = 〈|
p|
2〉 was calculated by averaging over an ensemble of 1000–2000 random realizations. The following three examples show the effect of Gaussian fluctuations in the speed of sound. In
Figure 2 there is a graph of transmission losses for a waveguide with bottom parameters
c1 = 1600 m/s, (ρ
1/ρ
0) = 2, β
1 = 0.01. In such a waveguide, there are 11 propagating modes that fluctuate and couple with each other. This example is illustrative, since the model bottom chosen is sufficiently rigid, with the result that mode attenuation with a distance becomes noticeable at relatively large distances
r > 10–15 km. The curves in the graphs are given starting from a distance of 200 m from the source. It is clearly seen that, compared with the deterministic model, there is a noticeable slowdown in the decay of the average intensity of the sound field with a distance.
As a result, the level of average intensity gradually rises above the level of averaged (over spatial oscillations) intensity corresponding to deterministic propagation conditions (dashed curve). Moreover, the stronger this rise is expressed, the larger the horizontal scale
Lr of the correlation of inhomogeneities, in full accordance with the above analytical estimates for Gaussian fluctuations. So, for the scale
Lr = 50 km, the level rise reaches 23 dB to a distance of 50 km, for the scale
Lr = 1 km, the level rise at a distance of 50 km is 3−7 dB depending on the observation horizon. We emphasize that the transmission losses for the scale
Lr = 50 km are presented purely for demonstration purposes, to confirm the analytical estimates. Typical horizontal scales for modes of internal wave perturbations are within the range of 1–10 km [
8,
9,
14,
18]. Based on the estimates made (9)–(12), it is also obvious that the effect is stronger, the greater the intensity of fluctuations, and the greater the modal absorption coefficient of sound by the bottom, since in this case the imaginary parts of the eigenvalues fluctuate more strongly (see (4) and (5)). The latter circumstance takes place, in particular, when the radiation frequency increases. In contrast, lowering the frequency of the sound reduces all statistical effects. Another important parameter on which the magnitude of the effects under consideration depends is the degree of penetrability of the waveguide bottom for sound waves, which is determined by the value of the refractive index at the water-bottom interface (
c0(0)/
c1)(1 +
iβ
1) [
26].
Above, a fairly rigid bottom boundary was considered. We now turn to waveguides with greater bottom penetrability, and we will call the bottom boundary “rigid” if
c0(0) <
c1, and “soft” in the contrary case
c0(0) >
c1.
Figure 3 demonstrates the transmission losses in the waveguide with the parameters:
c1 = 1530 m/s, (ρ
1/ρ
0) = 1.5, β
1 = 0.01. In this case, four propagating modes are excited in the unperturbed water layer, and four leaky modes were additionally taken into account in the calculations (the reference was carried out to the Pekeris cut on the complex plane of κ). The imaginary parts of the horizontal wave numbers κ
m of modes for a given waveguide are noticeably superior to those of the previous waveguide, so the sound field decays with distance much faster. The source location corresponds to 8 m from the bottom, where the maximum of the first eigenfunction of the deterministic waveguide is located; two observation horizons are considered: 8 m and 24 m from the bottom.
Figure 3 shows that the slowdown in the decay of the average intensity is manifested already at fairly close distances of 8–10 km from the source, and to a distance of 30 km it reaches 20–23 dB.
In
Figure 4 there is a graph of transmission losses in the waveguide with the “soft” bottom:
c0(0) >
c1 = 1200 m/s, (ρ
1/ρ
0) = 1.5, β
1 = 0.01. For this waveguide, no propagating mode is excited within the water layer, the calculations took into account eight leaky modes, which form the acoustic field in the deterministic waveguide at distances from the source of interest for statistical analysis. In this case, the sound field decays even faster than in the previous example (dashed curve in
Figure 4), and the absorption value in the bottom β
1 does not play a practical role and can be set equal to zero. As can be seen from
Figure 4, the effect of slowing the decay of the average intensity is noticeable already at distances of 5–7 km from the source, and to a distance of 30 km it can exceed 60 dB. The significant effect in this waveguide compared with the previous ones is due to stronger fluctuations of the modal wave numbers κ
m. The consequence of this is the rapid development of strong fluctuations of the intensity of the sound field, which on average reduces the transmission losses.
Let us now compare the simulation results for waveguides with Gaussian statistics of fluctuations ε(
r,
z), presented above, with results for similar waveguides, but having non-Gaussian fluctuations ε. Consider random fields of the form:
where ε
1(
z) is Gaussian random process and ε
2(
r) is non-Gaussian one. These processes are statistically independent and have the same correlation function as before:
. One of the well-known non-Gaussian random processes is the telegraph process, which is widely used in various fields of economics, mathematics, physics, and optics [
33]. For the description of random perturbations in underwater acoustics, the telegraph process was not widespread, but it can serve as a model of weak nonlinear cnoidal waves having a discontinuous amplitude, which are often present in the composition of background internal waves on the sea shelf. This process is defined as follows: ε
2(
r) = σ
2(−1)
n(0,r), where n(0,
r) is the Poissonian random process with a probability distribution P
n(r1,r2)=n = 〈n(
r1,
r2)〉
n exp[-〈n(
r1,
r2)〉]/n!, where 〈n(
r1,
r2)〉 is the average number of points on the interval (
r1,
r2), 〈n(
r1,
r2)〉 = ν|
r1-
r2|. If at
r = 0, ε
2(0) = ±σ
2 is equiprobable, then ε
2(
r) is a stationary process, having the average 〈ε
2(
r)〉 = 0, and the correlation function
,
Lr = (2ν)
−1.
In
Figure 5 for two horizons, the transmission losses are presented in a waveguide with Gaussian fluctuations, as in
Figure 3, and with fluctuations of the form (14) using the telegraph random process ε
2(
r) described above. It is clearly seen that the telegraph process noticeably underlines the weakening of the decay of the average intensity in the waveguide. The corresponding curves begin to diverge at
r > 8–10 km, and the loss attenuation to 30 km, caused by the difference in the fluctuation field from the Gaussian one, is 10 dB or more. In this case, the curves become somewhat more cut due to the discontinuous nature of the telegraph process ε
2(
r). For a waveguide with a “soft bottom” (
Figure 6), the picture is similar, only the discrepancy with Gaussian average intensity curves increases towards attenuation of propagation losses to a distance of 30 km, additional attenuation is 20 dB or more for different observation horizons.
The following example demonstrates the situation with a log-normal probability distribution of the fluctuations ε
2(
r). Log-normal distributions are quite often encountered in statistical wave problems [
2,
33]. Main feature of the such distribution in comparison with Gaussian one is the asymmetry and the presence of slowly falling tails, which leads to the appearance in the ensemble of random realizations of rare but strong emissions in intensity. The probability density of the generated process ε
2 is shown in
Figure 7. As before, ε
2 has zero mean value and is described by an exponential correlation function. The median of the distribution is non-zero, but the distribution parameters are chosen so that the fluctuations are adequate in intensity to those observed in the experiment and correspond to the values described above. As applied to the propagation of low-frequency sound on the sea shelf, log-normal fluctuations in the speed of sound can be viewed as the result of the influence of transient hydrodynamic perturbations leading to the formation of background internal waves (formed on the tails of numerous collapses of nonlinear internal waves), which are commonly described by Gaussian processes.
From the transmission loss curves given below, it can be seen that the effect of reducing the decay of the average intensity is preserved also with log-normal fluctuations in the speed of sound. Compared with Gaussian fluctuations (thin curves in
Figure 8), at distances of 10–20 km the effect can be 2–4 dB, depending on the observation horizon.