Learn how hybrid representation techniques are applied using Fock’s contour integral to analyze the attenuation factor W(ξ, h, h0).
In this article we use a Fock’s contour integral representation for the attenuation factor W(ξ, h, h0). As shown in Understanding Parabolic Approximation in Wave Propagation: Analytical Methods and Applications“Parabolic Approximation to the Wave Propagation”, both the contour integral representation and the expansion over the set of eigenfunctions of the continuous spectrum are equivalent in representing the attenuation factor W(ξ, h, h0) in the absence of random inhomogeneities of the refractive index. From practical point of view, in a deterministic problem the manipulations with a contour integral are a bit less cumbersome than eigenfunctions of the continuous spectrum, since the latter contains terms, such as waves χ–(t, h) coming from infinity that produce a negligible contribution to the integral for positive ξ in the absence of random fluctuations.
We consider a slightly more general case of the M-profile for the hybrid ray-mode representation, as shown in Figure 1.

The contour integral for the attenuation factor W(ξ, h, h0) can be written as:
The function K(h, h0, t) satisfies the equation:
and the conditions:
where q = jmZ and Z is a surface impedance of the interface h = 0. The integration contour embraces singularities of the integrand lying in the first quadrant. The precise form of K depends on the shape of the profile U(h) and the relative values h, h0, Hs and Hi. If both the transmitter and the receiver are within the lower inversion layer, i. e. h0 < h < Hs, then K takes the form:
where f1 and f2 are the independent solutions of uniform equation (2) specified within 0 < h < Hs, and:
The values Rg and R are defined as:
The functions φ1 and φ2 are independent solutions of the uniform equation (2) for Hs < h < Hi, and ψ1 is the solution for h > Hi which satisfies the condition:
The derivative sign f ′ in Eqs. (6)–(8) means the derivative over the h-variable. Apparently, Rg, R and Ri can be regarded as the reflection coefficients at h = 0, h = Hs and h = Hi, respectively.
Now we can get to the hybrid representation of the received field, as manifested in the beginning of this article. As a first step toward implementing this program we will express the factor (1-RiRg)-1 in terms of a partial geometric series:
and substitute Eq. (10) into the integral (1):
Here:
and N is an arbitrary integer whose choice is dictated by the length of the propaga- tion path. It represents the maximum number of reflections from the upper “boundary” h = Hi to be taken into account.
If there is no elevated M-inversion (cf. the dashed profile in Figure 1), the attenuation function W(ξ, h, h0) would take the form:
where Rs denotes the reflection coefficient at the “upper wall” h = Hs of the surface duct, i.e:
Of the two independent functions above Hs, i. e. φ1 and φ2, only φ1 will remain because of the “new” condition:
To isolate the set of eigenwaves of a solitary surface duct from the solution (11), we multiply and divide the first term of Eq. (11) by the “resonant denominator” of Eq. (13), i. e., (1-RsRg), and represent R as:
with:
Now consider the function g(h, t) defined as the solution of:
with the conditions:
where U1(h) corresponds to the dashed M-profile of Figure 1. We can obtain that:
and hence Eq. (19) is equivalent to the characteristic equation 1-RsRg = 0 specifying the modes of the surface duct. Further, we bring the integrand of Eq. (11) to the form:
Representing (1-RsRg)-1 in the second and third term as another partial geometric series, gen- erally speaking with a number of terms different from Eq. (10), we assemble similar terms and arrive at:
where term W01 represents the sum of the normal modes of the surface duct:
and:
where:
The term T2 can be interpreted as the field arriving at the observation point due to a single reflection from the elevated inversion while T3 is a set of waves multiply reflected between the boundaries h = Hi, h = Hs, and h = 0. T4 is a remainder term whose smallness will be secured by the proper choice of N for every specific separa- tion from the source. It should be emphasized that in bringing Eq. (11) to the form of Eq. (23) we made no approximation, and hence Eq. (23) is an exact representation of the field in the structure analysed (to be more precise, it is exact in the same sense as the initial formula Eq. (1), which itself corresponds to a parabolic equation approximation).
Read also: Basics of Fock’s Theory – the Hartree-Fock Approach to Evaporation Ducts
The set of normal modes associated with the surface duct, the entire infinite spectrum, has been separated from Eq. (1) solely by transforming the integrand (4). The fields given by T2 and T3 so far cannot be interpreted in the ray optical form, since rays appear only at the stage of asymptotic evaluation of the integrals.
As for the numerical truncation of the modal series, the question seems almost trivial, since in the absence of a reflecting upper wall the higher-order modes are evanescent and their respective eigenvalues have the ordering:
As observed from Eqs. (25) and (26), the integrands T2 and T3 contain no singularities. Depending on the number of significant saddle points, each of these integrals splits into several terms representing waves which are radiated either up or down from the transmitting antenna. Further along the path, they are reflected from the boundaries h = 0, h = Hs and h = Hi to arrive at an observation point along some of the ray trajectories of Figure 2.
For example, T2 splits into two terms, each corresponding to a wave that is reflected once from the elevated layer. The first leaves the source in an upward direction, reaches a reflection point in the upper layer, and comes to the observation point along the ray 1, Figure 5.7. Another emerges downwards from the source downward, undergoes reflections from the earth’s surface and the elevated layer, and arrives at the observation point along trajectory 2, Figure 2.

The different components of the angular spectrum of waves excited by the transmitter dipole are subject to a kind of spatial “filtration” in the nonuniform troposphere. We can define the critical angle θ0 as the angle at which the ray QO′ emerging from the source turns at the height h = Hs, as shown in Figure 2. The rays with departure angles θ < θ0 are trapped by the surface duct; those with θ > θ0 leave it freely.
The ray with θ = θ0+δ, δ → 0 comes back to earth, upon reflection from the elevated layer, at the maximum range attainable to a single reflection from the elevated layer, we define this range as ξ1m. Real rays arriving at greater distances than ξ1m can do so only through a higher number of reflections from the “boundaries” h = Hi, h = 0 and h = Hs. Generally, each group of rays of a given reflection multiplicity, l, has its own “horizon”, i. e. the limiting distance still corresponds to real stationary points in the integrand of T2 (for l = 1) and/or T3 (l > 1). It is in the vicinity of the lth horizon that the contribution of l times the reflected wave is the highest.
Indeed, the grazing angle with respect to the elevated layer assumes the lowest possible value when the observation point is near the horizon. As a result, one can distinguish, along the path, characteristic zones of the “first hop” (i. e. near n = ξ1m where is the dominant contribution to T2+T3+T4, for the N terms in T3, only those with l ≤ N “work“). Moving beyond the lth horizon, an observer finds himself in an umbral zone where the field of the l+1 hop has not yet formed while that of the lth hop can no longer propagate except by diffraction. The magnitude of this diffractional field is determined by the small contribution of complex rays corresponding to the complex stationary points of T2 and/or T3, but mainly by the value of the remainder integral T4. For evaluating this latter it seems convenient to represent it as:
with:
As can be seen, T41 is the contribution of waves N+1 times reflected from the upper layer, and hence it should be small near the horizon of l ≤ N-reflected waves.
The number of reflections, N, which are to be taken into account at a specific distance from the transmitter can be estimated as follows: Let us assume for simplicity that the reflection coefficient Ri depends on the grazing angle θ as:
The saddle point of Eq. (33) is:
Recalling the relation between the grazing angle and t, i. e.:
and the fact that real saddle points of real rays satisfy the inequality:
we can obtain from Eq. (35) an estimate at t → 0 for the grazing angle θ of N times the reflected waves:
To have Ri(θmin) ≠ 0, it is necessary that:
The order of magnitude of θ0 is 10–2, and that of (2Hi/a)1/2 ≅ 5·10-3 to 10-2, whence the number of terms to be retained in Eq. (23) near the horizon l = 1 is N = 1–3.
We leave evaluation of the magnitude of term T42 in Eq. (31) until the next section, noting here that its contribution cannot be described in pure ray optical terms. In fact this term describes a peculiar physical effect of secondary excitation of the surface (evaporation) duct.
A better illustration of the relative importance of the different terms in Eq. (23) can be achieved through numerical analysis. Figure 3 represents the calculated attenuation function for f = 10 GHz in a two-channel structure of the refractivity. The surface M-inversion is characterised by Zs = 9 m, ΔM = M(0)-M(Zs) = 2 M-units. The elevated M-inversion lies at the height Zi = 500 m and has been approximated with a discontinuity ΔMi = 10 M-units of no thickness (i. e., Zk = Zi).

An approximation like this obviously results in an overestimated value of the reflection coefficient Ri(t). The heights of the transmitting and receiving antennas are z0 = 8,5 m and z = 5 m, respectively, apparently both inside the evaporation duct. The normalised impedance q of the earth’s surface is assumed to have a value corresponding to that of sea water: Req = 6, Im q = 94. The calculated imaginary part of the first-mode eigenvalue is found to be Im t1 = 0,4. For a given combination of the antenna heights the optical horizon would be at ≈ 20 km. As observed from Figure 3, there is a region beyond that horizon (extending to about 60 or 70 km) where the field level is determined by the lower-order mode of the surface duct alone. Further along the path, the contribution of waves reflected from the elevated layer becomes substantial. In the zone of the first hop (i. e. 60 km ≤ x ≤ 200 km) the maximum field strength becomes close to the free space level. The remainder term T4 reaches a –40 dB level near distance 200 km with N = 1.
The mechanism of the radio wave propagation in the above two-channel system can be summarised as follows. At relatively smaller distances from the transmitter the dominant contribution to the received signal strength comes from evaporation duct propagation alone. Further along the path, the waves, departing with angles exceeding the critical angle θ0 and reflected by the elevated M-inversion, contribute to the received field at an appropriate distance of the single hop. These waves interfere with each other as well as with the waves propagating through the evaporation duct producing the deep fades caused by large phase difference. For example, in Figure 3, the deepest fades are observed in the range 60–100 km where the wave propagating through the evaporation duct and superposition of the waves reflected from the elevated M-inversion are comparable in amplitude. Then, along the distance, the field reflected from the elevated M-inversion takes over and provides the major contribution to the received field strength.
Secondary Excitation of the Evaporation Duct by the Waves Reflected from an Elevated Refractive Layer
Referring to the discussion in the previous section and Figure 2, we state that the rays with departure angles from the source less than critical angle θ0 are trapped in the surface duct, those with θ > θ0 leave it freely. We also note that rays with θ > θ0 reach the observation point at the distances ξ < ξ1max. At distances greater than the horizon of first hop (ξ1max, single reflection) there are no single reflected rays launched at the angles of departure θ > θ0. Hence, between the first and second hop there is a range of distances where the receiving field is formed by another mechanism represented by term T42 in Eq. (31) and evaluated in this section.
Consider the range of distances ξ ≥ ξ1max concentrated on trapped waves/rays. The trapped waves with departure angles θ < θ0 turn inside the surface duct, as shown in Figure 3. Since the potential barrier is finite the trapped waves can leak to the space above the surface duct as shown schematically in Figure 4, their respective phases are complex due to the leaking mechanism distributed along the distance. The field created by these waves in the space above the surface duct can be imagined as a cluster of rays distributed along the distance and sliding along the tangents to circles of radius r = a+h2, where the height h2 is given by the relation 2m2 10–6 M(h2) = t for h2 > Zs, where t is a dimensionless “energy” corresponding to the angle θ.

In evaluating the maximum magnitude of T42 and therefore the impact from sec- ondary excitation at ξ ≥ ξ1max we shall set N = 1. Taking into account that the observation point lies within the surface layer, we will represent the integral (31) as a sum of residues at the zeros of 1–RsRg = 0.
Using Eq. (16) we can isolate the coefficient of reflection from the upper boundary of the surface duct Rs in R and present the denominator of Eq. (31) in the form:
The term v may be interpreted as a coupling factor between the lower and upper ducts. For low attenuating modes (Im tn ≪ 1) it can be estimated as:
where μ = |dU/dh| is a gradient of profile U(h) for 0 < h < Hs. The greater Re tn the lower is the turning height h1, the lesser the amount of leaking from the surface duct and, in accordance with a reciprocity theorem, such waves are harder to generate by a source distributed across the reflective layer but outside the surface duct.
To finish evaluation of T42, we write a partial geometric series for (1-RsRg)-1 (1-RRg)-1 with the help of formula (39), namely:
and substitute it in Eq. (31).
Now we restrain the evaluation by the first term of the expansion (41) using the same arguments of angular filtration as in the previous section. Then we assume that Hi-Hs ≫ 1 and a/Hi·10-6 ΔM ≪ 1, where ΔM = M(z = 0)-M(Zs) is the M-deficit of the surface duct. Therefore, we assume that the waves reflected from the elevated refractive layer can be regarded as a plain waves on their arrival at the level of the surface duct, and the horizons of departure of the leaked waves h2 in Figure 4 are not to too different from Hs compared with Hi, i. e., max(h2)/Hi ≪ 1. This allows us to use the same magnitude of the reflection coefficient |Ri(t)| ≈ |Ri(0)| for all rays from the cluster of leaking waves, expanding the amplitude and phase of the reflection coefficient into a series over powers of t. Taking into account that the observation point lies within the surface layer, we will represent the integral (31) as a sum of residues at the zeros of (1-RsRg)2 = 0.
Here S is the number of modes trapped in the surface duct.
Equation (42) provides an estimate of the secondary excitation of the surface duct by a limited number of the modes of the surface duct, “trapped”, from the point of view of geometrical optics but, in fact, “leaking” from the duct along the distance of propagation due to the finite thickness of the potential barrier in U(h). These waves, then being reflected from the elevated M-inversion, reach the surface duct from above and “leak” back into the wave guide reproducing itself in the surface duct. Of course, the symmetry in the secondary excitation of the surface duct relies on the uniformity of the refractive index in the horizontal plane along the distance. Analysing Eq. (42) we may note that for the shorter wavelength λ of the electromagnetic radiation, the trapping condition of the surface duct becomes stronger and leaking of the trapped modes is reduced, as is the attenuation of the ducted field with distance. This leads to an even smaller contribution from the secondary excited field described by term T42, Eq. (42).
Increase in the wavelength λ leads, on the one hand, to a greater attenuation of the ducted field and stronger leaking of the trapped modes from the duct and, on the other hand, to better conditions for reflection from the elevated layer and stronger secondary excitation of the surface duct. The major contribution to the secondary guided field is from the modes with |Re tn| ≪ 1. These modes have relatively high attenuation compared with deeply trapped modes, however they still provide a significant contribution at smaller distances from the source and, being reflected from the elevated layer, at distances ξ ≥ ξ1max. Evaluating this result numerically, we have found that the secondary guided field given by Eq. (42) may reach a maximum level of –50 dB relative to a free space condition with ideal reflection from the elevated layer, |Ri| = 1.
The rather exotic though observable situation is shown in Figure 5 where we have extremely strong elevated M-inversion and surface-based evaporation duct simultaneously. Here Zs = 18 m, ΔM = M(0) – M(Zs) = 6 M-units, Zs = 500 m, Zk = 600 m, ΔMs = 60 M-units. The calculation was performed for radio frequency f = 3 GHz and the following combination of antenna heights:
- z0 = 16 m;
- z = 10 m and 50 m, respectivly.
Figure 5 shows that the field reflected from elevated M-inversion dominates at distances of the order of the reflection cycle length Λi despite there being lots of trapped modes of the elevated duct.

At shorter distances the total field is a combination of the modes of the evaporation and elevated ducts with predominance of the evaporation duct. It also might be observed that at distances exceeding the first hop distance but still shorter than the second hop distance (in the above case with x > 200 km) the secondary excited field of the evaporation duct may provide a substantial contribution revealing a distinct change in the attenuation rate with distance and somewhat repeating the field at short range.
Figure 2 shows the ray trajectories for this situation. As observed, in the vicinity of the first hop the waves reflected from the elevated M-inversion provide a dominant contribution to the received field.