.
Our site needs your help!
Site categories

Excitation of Waves in a Continuous Spectrum and Evaporation Duct with Two Trapped Modes

Join Our Telegram (Seaman Community)

Explore the phenomenon of wave excitation in evaporation ducts with a continuous spectrum, including the case of two trapped modes, and gain insights into the role of statistical inhomogeneity in shaping wave propagation patterns.

In the absence of fluctuations in refractive index, δε, the signal strength in the shadow region above the evaporation duct is exponentially small because of the smallness of the “sub-barrier leakage” of the trapped modes.

Excitation of Waves in a Continuous Spectrum in a Statistically Inhomogeneous Evaporation Duct

The contribution of the waves of the continuous spectrum, initially excited by the transmitter, can be neglected because of the diffraction attenuation in the shadow region. The scattering of the waveguide trapped modes by random non-uniformities of the refractive index gives rise to incoherent exchange of energy between trapped modes and excitation of waves in the continuous spectrum. In addition, the initially excited waves of the continuous spectrum may reach beyond the horizon due to a scattering in the upper layers of the troposphere (Booker–Gordon’s single scattering mechanism). In this section, we study the contribution of the waves of the continuous spectrum, whose source is a waveguide mode, to the intensity of the field, i. e. the mechanism of the excitation of waves of the continuous spectrum by trapped modes of the evaporation duct scattered on the imhomogeneities in the refractive index.

Recalling the equations for the coherence function from the previous section, we concentrate on investigation of the spectral amplitude of the component of the coherence function related to a continuous spectrum:

g~cx, E1, E2, p, γ=12π0dz10dz2dY·Γx, Y, γ, z1, z2 ΨE1*z1 ΨE2 z1 exp jpY              Eq. 1

As is well known, in the presence of random fluctuations of the refractive index, δε, the coherent component of the waves of the continuous spectrum attenuates with distance due to incoherent multiple scattering. The decrement of attenuation in the coherent component is given by:

γc=πk22d2κΦε0, κ              Eq. 2
where
Φε0, κ
is the spectrum density of the fluctuations in δε. As observed in Understanding Parabolic Approximation in Wave Propagation: Analytical Methods and Applications“Green Function for a Parabolic Equation in a Stratified Medium”, in the case of free-space propagation with random irregularities of δε, the attenuation of the coherent component in the field intensity is completely compensated by the inflow of the incoherently scattered field from other directions, and the total intensity of the wave field does not reveal exponential attenuation. Since the propagation above the evaporation duct is free-space-like we may use the above argument in support of neglecting multiple scattering of the wave of the continuous spectrum in a space above the evaporation duct.
However, the effect of multiple scattering of the wave of the continuous spectrum can obviously be neglected if we demand a smallness of the attenuation of the coherent component of the CS spectrum wave, i. e., γc, xa ≪ 1, over the path xa of the wave from the point where scattering of the wave guide mode into the CS wave occurred up to the point of observation. Evaluating the maximum value of the distancexa as
xa max=2az,
where z is the height of the observation point, we obtain the limitation
z2aγc21.

In the case of Kolmogorov’s turbulence:

γc=2,66k2Cn2L05/3.
Assuming
Cn2=1015 cm2/3
, L0 = 104 cm, k = 2 cm–1, we obtain the result that the maximum heights under consideration should not exceed the value zmax ~ 103 m. Thus, the approximation in which the estimate is obtained below for the waves of the continuous spectrum is applicable for altitudes of the receiving point z satisfying the inequality z ≤ zmax.

Let us study the intensity of the field of the waves of the continuous spectrum:

Icx, z=Γx, ρ, ρ, ρ=0, z,

setting the y-coordinates of the receiving and transmitting antennas to zero. Then, integrating Coherence Function in a Random and Non-uniform AtmosphereEquation and employing Coherence Function in a Random and Non-uniform Atmospherethe expression for the component of the coherence function associated with waves of the discrete spectrum, Γd, we obtain:

Icx, z=ka216xϕdz02 exp γdx·0xdxdqQ<q/2dQBq, Q ΨQ+q/2z ΨQq/2*z exp jqx2k+γdx                Eq. 3

where we have introduced q = E1-E2, Q = (1/2)(E1+E2):

Bq, Q=dκγdκzΦε0, κγ, κz·0dz10dz2 ΨQ+q/2*z1 ΨQq/2z2ϕdz1ϕd*z2 exp jκzz1z2.                Eq. 4

The integral B(q, Q) determines the coupling between the waves of the discrete and continuous spectra due to the scattering on the inhomogeneities of δε. The major contribution to B(q, Q) comes from the spectral components of δε with the wavenumbers:

κz>μ=dU/dz1/3

with z < Zs and the region of the altitudes z1,2 in the neighbourhood of the turning point zd of the mode of the discrete spectrum, where U(zd) = Ed. We use the WKB approximation for the waves of the continuous spectrum and the Airy function representation for the waveguide mode in the vicinity of the turning point:

ΨEz=2πUzE1/4 exp j3π4 sin 0zdzUzE,                Eq. 5
ϕdz=1Ndνμzzd, Nd=0dz ϕdz2.

The spatial spectrum of the fluctuations in the refractive index we define in the form (Atmospheric Boundary Layer and Basics of the Propagation Mechanismsin this equation) for the case of turbulent fluctuations δε isotropic in the (x, y)-plane:

Φεκ=0,033C||2ακ||2+α2κz211/6            Eq. 6

we substitute into Eqs. (4) and (5) the integral representation of the Airy function:

νt=12πdξ exp jξ33+jξt

and require the inequalities:

κz>μ and 2Edmk1              Eq. 7

to hold, which is necessary for the entire scheme of the solution obtained in Spectrum of Normal Waves in an Evaporation Duct“Detailed Wave Spectrum Analysis in Evaporation Ducts”. The meaning of inequality (7) is that we do not take into account the scattering of the leaky modes, assuming that the waveguide field is composed only of trapped (non-attenuated) modes. The B(q, Q) can then be truncated to:

Bq, Q=Cε2α5/3·0,089sin qΛdQqNd2EdQ4/3,   κκd10,053cos qΛdQμ3Nd2EdQ5/6,   κκd1              Eq. 8

where:

Λd=120zddzUzQ,     κd=qdQ,     κ=U0QEdQ          Eq. 9

Parameter κd determines the Bragg scattering angle θs from the waveguide mode into the wave of the continuous spectrum with energy Q at the height of the turning point z = zd, i. e., 2 sin θs/2 = κd/k. The parameter Δk is a difference between the magnitude of the scattering angle at the surface z = 0 and at the altitude of the turn-ing point z = zd. Thus the inequality Δκ/κd ≪ 1 requires that the change in the scattering angle owing to refraction in the scattering volume, i. e. in the region of the mode’s localization, be much smaller than the scattering angle itself at the altitude zd.

The relationship between the characteristic “energyQ and the sliding angle θ of the wave front of the wave of the continuous spectrum relative to the surface z = Zs is given by the equation Q = -k2/m2 tan2θ. The parameter Λd(Q) determines the distance along x traversed by the wave of the continuous spectrum with “energyQ as it propagates from z = 0 to z = zd. Thus, B(q, Q) defines the angular distribution of the scattered field, i. e., the scattering function in the “directionsE1, E2.

To calculate the intensity of the field lc(x, z) we employ the WKB approximation (4) and (7) in Eq. (2). We leave aside the constant terms for later and concentrate on the integration in Eq. (2). The integration over the “energy” difference q leads to the appearance of terms that place the limits on the region of the distances x from which the scattered field is collected in the receiving antenna. With
U0qd/qd1
we obtain:
Icx, z~0dQqdQ4/3UzQ1/20xdx exp γdx·12θΛd+Λxk12θΛdΛxkθΛdxk cos 20zdzUzQ       Eq. 10
and for
U0qdqd1 :
Icx, z~0dQqdQ5/6UzQ1/20xdx exp γdx·12δxk+ΛdΛ+12δxkΛd+Λ+12δxkΛdΛδxkΛd cos 20zdzUzQ                Eq. 11

where:

θx=1,     x>01/2,   x=00,     x<0

and:

Λ=ΛQ, z=120zdzUzQ.

We introduce the dimensionless variables τ = -Qm2/k2, τd = Edm2/k2, h = kz/m, β = m2γs/k and the modified profile U(h) = m2/k2U(z). We then require that the inequality:

kγdΛd01              Eq. 12

holds. Since Λd(0) ≥ Λd(Ed), where d(Ed) the length of the cycle of a waveguide mode and the contributing values of Q are limited by the inequality |Q| ≤ Ed, the condition (12) demands the smallness of the attenuation of the waveguide mode at the distance equal to its cycle. It may be noted that in the opposite case there is no waveguide mode structure to sustain scattering, the mode will be destroyed during one cycle.

Finally we obtain the intensity of the wave of the continuous spectrum normalized at the intensity of the free-space field:

Jcx, z=Icx, zx2a2=0,078k5/3 xCε||2 α5/3m11/3ϕdz02Nd2·Kh exp γdx             Eq. 13
and the function K(h), which commands that the height distribution of the scattered field is given by two equations obtained from Eqs. (10) and (11). With
U0τd/τd1 :
Kh=0dττd+τ4/3Uh+τ1/2Λdτ·exp βΛτ, hcos20hdhUh+τ,              Eq. 14
and for
U0τd/τd1 :
Kh=πdUdh2/3 h=hd0dττd+τ5/6 Uh+τ1/2·exp βΛτ, hcos20hdhUh+τ.              Eq. 15

We shall now study the total field in the shadow region as being the result of composition of the waves of the discrete and continuous spectra excited due to the random scattering of the waveguide modes. For the intensity of the waveguide mode the relation below follows:

Jd=x8πkNd4 ϕdz02 ϕdz2 exp γdx.              Eq. 16

We define the total intensity as:

Jtot=Jd+Jc=x8πkNd2 ϕdz02 exp γdx Sz             Eq. 17

where:

Sz=ϕdz2Nd2+0,139k2/3 Cε2α5/3 m11/3 Kz.            Eq. 18

The function S(z) determines the height distribution of the total intensity of the field. We examine the behavior of the function S(z) in the case of a bilinear dependence of U(z). We assume the following values for the parameters of the problem:

  • k = 2 cm–1;
  • a = 8 500 km;
  • Zs = 11 m;
  • ΔεM = εM(0)-εM(Zs) = 5,8·106.

The modified profile U(h) can be defined as:

Uh=μ13Hsh,   hHshHs,        h>Hs           Eq. 19

where:

  • Hs = kZs/m = 2,42;
  • μ13=m2εM/Hs=2;
  • τd=U0μ12τ1=1;
  • U(0) = m2ΔεM = 4,34;
  • τ1 = 2,338.

Consider the height dependence of the waveguide mode, i. e., the first term in S(z). Thus far in the analysis of the waveguide field we have neglected the leakage of the trapped waves of the discrete spectrum through the potential barrier. The attenuation caused by the effect of the leakage is accounted for in the imaginary part of Ed, which can be written as:

δd=Im Ed=μ2μ124τ1 exp 43τd3/21+μ13μ13.             Eq. 20

Sub-barrier leakage has virtually no effect on the height structure of the trapped mode inside the waveguide channel for h < Hs, where:

ϕdh=1Ndντdμ12μHsh.             Eq. 21

Outside the waveguide channel the height structure of ϕd(h) corresponds to an outgoing wave with amplitude proportional to the leakage factor δd:

ϕdh=μ11/3μNdτ11/4 δd1/2 w1 τd+jδdμ2hHs.             Eq. 22

Examining the second term in S(z) given by Eq. (18), we can express the coefficient in front of K(z) in terms of the non-dimensional attenuation coefficient b using Coherence Function in a Random and Non-uniform AtmosphereEquation for γd from the previous section:

β=m2γd/k=0,264k2/3 Cε||2α5/3 m11/3 μ15/9.

Then for S(z) we obtain:

Sz=ϕdz2Nd2+0,492μ15/9 βKz.             Eq. 23

In principle, in order to determine the coefficient β it is necessary to know the magnitude of the structure constant Cε and the anisotropy parameter of the irregularities of the refractive index in the volume of the waveguide channel. In some cases when the radio signal strength is measured, the measured data may be deployed to estimate the height dependence of the received signal. In particular, when the magnitude of the field attenuation per unit length γx = dI/dx is known from the measured signal, then γd is related to γx via γx = 4,34γd, and given knowledge of the average M-profile we can construct a theoretical height structure of the field intensity S(z).

The γx magnitude for centimeter waves lies in the range γx = 0,2–0,5 dB km–1. The corresponding limits for β are β = 0,207–0,52 for k = 2 cm–1 and a = 8 500 km (normal refraction above the evaporation duct). Thus under real conditions the coefficient in front of K(z) in Coherence Function in a Random and Non-uniform AtmosphereEquation is of the order of unity and, consequently, the contribution of the waves of the continuous spectrum to S(z) can be significant.

Figure 1 shows the result of the calculation of the resulting height distribution 10 log S(h) based on Eq. (23) for β = 0,207 and β = 0,52 with U(h) given by Eq. (19) and its parameters defined above.

Height structure
Fig. 1 Height structure of the wave field in the presence of an evaporation duct: f = 10 GHz, Curve 1, no scattering on inhomogeneities of the refractive index, curves 2 and 3, scattering is taken into account

The figure also shows the height dependence of the first term in Eq. (23), i. e., the contribution of the waves of the discrete spectrum alone. As observed from the figure, the intensity of the field above the turning height hd (in this case hd = 1,67) is a contribution from the waves of the continuous spectrum, i. e. second term in Eq. (23). Apparently, scattering on random fluctuations in the refractive index leads to a significant change in the height dependence of the signal strength; in the presence of random fluctuations de in the refractive index there is no sharp exponential decay in the average signal strength outside the duct and for sufficiently strong fluctuations the δε field outside the duct may increase to the order of magnitude of the field inside, thus revealing a rather smooth unpronounced height dependence compared to the case of deterministic duct propagation.

At relatively large altitudes, h = Hs ≫ τd, the function S(h) has the following asymptotic form:

Sh~14hHsτd exp 43τd3/2μ13+1μ13+2δdμ2hHsτd + 1,58μ15/9βτd1/3hHs exp βhHs.              Eq. 24
Beginning with the altitudes
h1=Hs+τd+μ42δd2
and
h2=Hs+14β2,
the exponential growth predominates over the cylindrical divergence in either the first or second terms, respectively.

The exponential growth of those terms in Eq. (24) is due to the arrival of the waves from shorter distances x < x, at which the field of the waveguide mode was exponentially large, compared with its value at distance x. The second term makes the main contribution in Eq. (24), since the exponential factor β exceeds the sub-barrier leakage factor δ by at least an order of magnitude.

As mentioned before, we ignored the scattering of the waves in the space above the waveguide. While, in principle, it should be taken into account for the sake of consistency, the qualitative character of the total field will not be expected to change drastically. We may assume however that additional scattering of the waves of the continuous spectrum in the space above the evaporation duct will smooth the exponential dependence of the scattered field in Eq. (24) for large altitudes. At the same time additional divergence of the waves due to a broadening of the angular spectrum can be compensated by the arrival of the scattered waves from the shadow region, i. e. from distances
x<xa/mhHs.

Evaporation Duct with Two Trapped Modes

Finally, we provide a closed form solution for the field intensity in the case when the evaporation duct can trap two modes. In this case the intensity in the two-mode waveguide is given by:

Jx, z, z0=4πxI1·S1+I2·S2               Eq. 25

where:

I1x=A1eλ1x+A2eλ2x,
I2x=A3eλ1x+A4eλ2x,
λ1,2=Γ1+Γ22±12Γ2Γ124W12,
A1=I10λ2+Γ2W12I20λ1λ2,
A2=I10λ1+Γ1W12I20λ1λ2,
A3=I20λ2+Γ2W12I10λ2λ1,
A4=I20λ1+Γ1W12I10λ1λ2,
Γ1=1,603·As, Γ2=3,98·As, W12=0,316·As, and As=k2/3 Cε2α5/3m11/3.

Coefficients In0 = |χn(h0)|2 are the coefficients of the excitation of the mode with number n and function Sn(h) is a height distribution of the total intensity of the nth mode:

Snh=χnh2+0,492Γnμ15/3Pnh.             Eq. 26

The height gain function χn is given by:

χnh=μ1τnνtnUhμ12,       h<Hs             Eq. 27
χnh=δ1/2w1tnUhμ22,       hHs

where δ is the coefficient of penetration through the potential barrier into the space above the duct:

δ=μ14τnexp43μ23+μ13μ13μ23tn3/2,           Eq. 28
τn=32πn142/3.

The function Pn(h) is given by the integral:

Pnh=π/μ120dstn+s5/6Us+s1/2eΛs, hcos2Qs, h           Eq. 29

where:

Qs, h=23μ13U0+s3/2Uh+s3/2,     h<Hs           Eq. 30
Qs, h=231μ13U0+s3/2μ13+μ23μ13μ23+Uh+s3/2,     hHs
Λs, h=1μ13U0+s1/21μ13Uh+s1/2,   h<Hs          Eq. 31
Λs, h=1μ13U0+s1/2μ13+μ23μ13μ23+Uh+s1/2,   hHs.

It should be noted that the above formulas actually represent the incoherent sum of two trapped modes. Under very moderate assumptions for the intensity of the atmospheric turbulent fluctuations (as in Figures 2 and 3 below) the interference term in the coherence function may be neglected since phase fluctuations in this case are strong enough to make the above approximation valid.

Atmospheric turbulent fluctuation
Fig. 2 Range dependence of the signal strength at 10 GHz in the presence of an evaporation duct with two trapped modes

On Figure 2 the transmitter antenna is mounted at 10 m above the sea surface. Curves (1) and (2) correspond to the range dependence of the received field strength for the receiving antennas at the heights inside (1) and outside (2) the duct for the case of a uniform duct in the absence of random fluctuations in the refractive index. Curves (3) and (4) correspond to the same antenna elevations for a non-uniform duct with random and anisotropic fluctuations of the refractive index:

Cε2=1014 cm2/3, α=0,1.
Signal strength
Fig. 3 Height-gain dependence of the signal strength at 10 GHz in the presence of an evaporation duct with two trapped modes at 100 km from the transmitter

On Figure 3 The transmitter antenna is mounted at 10 m above the sea surface. Curve (1) shows the received field versus the height in the absence of random fluctuations in the refractive index. Curve (2) is a height-gain dependence for a non-uniform duct with random and anisotropic fluctuations of the refractive index:

Cε2=1014 cm2/3, α=0,1.

To illustrate the effect of scattering on random fluctuations in a two-mode evaporation duct we calculated the signal strength at frequency 10 GHz in the evaporation duct with parameters Zs = 15 m and ΔM = 5 N-units. Figure 2 shows a range dependence of the signal strength in the absence and presence of the fluctuations in the refractive index. As observed from the figure, at a height 10 m inside the duct the composition of the two trapped modes produces an interference pattern, creating fades. The second curve shows the range dependence at the height 150 m where the second trapped mode is dominant and therefore no interference pattern is observed. The height gain pattern in this case is typical for one that follows from the model of the evaporation duct, namely, revealing the maximum of the signal strength inside the duct, Figure 3. In the presence of fluctuations we observe significant attenuation of the field with distance as well as considerable changes in the height gain distribution, namely, the height structure becomes largely uniform.

In the range of frequencies from 3 GHz to 20 GHz the evaporation duct in most cases can support only a few trapped modes. The calculation of the field strength for a finite number of modes can be performed in a way similar to that described in this section. The surface-based and elevated duct created by an advection mechanism and normally associated with strong elevated inversion of temperature may trap a hundred modes and the above approach is inefficient. In that case other methods based on the application of diffusion theory for a parabolic equation can be employed to solve the problem of scattering in a multimode waveguide.

Author
Author photo - Olga Nesvetailova
Freelancer
Literature
  1. Jensen, N.O., Lenshow, D.H. An observational investigation on penetrative convection. J.Atmos.Sci., 1978, 35(10), 1924–1933.;
  2. Monin, A.S., Yaglom, A.M. Statistical Fluid Mechanics, Vol.1, MIT Press, Cambridge, MA, 1971, p. 769.;
  3. Gossard, E.E. Clear weather meteorological effects on propagations at frequencies above 1 GHz, Radio Sci., 1981,16(5), 589–608.;
  4. Tatarskii, V.I. The Effects of Turbulent Atmo- sphere on Wave Propagation, IPST, Jerusalem, 1971.;
  5. Gavrilov, A.S., Ponomareva, S.M. Turbulence Structure in the Ground Level Layer of the Atmosphere. Collected Data, Meteorology Series, No.1, Research Institute for Meteorological Information, Obninsk (in Russian), 1984.;
  6. Kukushkin, A.V., Freilikher, V.D. and Fuks, I.M. Over-the-horizon propagation of UHF radio waves above the sea , Radiophys. Quantum Electron. (translated from Russian), Consultant Bureau, New York , RPQEAC 30 (7), 1987, 597–620.;
Footnotes
Sea-Man

Did you find mistake? Highlight and press CTRL+Enter

Март, 21, 2025 73 0
Add a comment


Notes
Text copied