Discover the dynamics of UHF propagation in evaporation ducts in this detailed article. We present key measurement results, theoretical comparisons, and insights into perturbation theory for normal waves in a stratified troposphere. Topics include problem formulation, linear and smooth distortion, height functions, and the effects of linear-logarithmic profiles near the sea surface.
The most important characteristic of the troposphere in terms of radio wave propagation is a refractivity profile averaged over the coordinates tangential to the sea surface, i. e., an M-profile in a “stratified” troposphere:
where z is the height above the sea surface, the gradient 0,157 N-units m–1 takes into account the curvature of the earth in the case of normal refraction. According to radio-meteorological and refractometer world-wide data, the gradient of the modified refractivity gm = dM/dz is less than the critical gradient gc = –0,157 N-units m–1 for more than 50 % of the time of observation over the sea surface. In those conditions the surface M-inversion forms an evaporation duct, which signifi- cantly affects the process of radio wave propagation by trapping the waves radiated in a selective frequency band.
The simplest characteristic of the ducting properties of the surface based M-inversion is the “critical” wavelength λc defined as:
and:
Parameter ς is introduced in Ref. as ς = z/Zs thus providing a meaningful definition for:
with q(0) = 4, q(1) = 0.
Frequently observed surface M-inversions of height up to 15 m are capable of ducting radiowaves in a cm band with wavelength λ ≤ λc ~ 5 cm. With increasing wavelength the impact of the evaporation duct is weakened, though the attenuation of the waves with wavelength λ > λc may still be significantly less than in the case of normal refraction. Fock showed that the attenuation of the field in the shadow region follows the exponential law exp (-γx), where x is the distance from the horizon, and the magnitude of the attenuation exponent γ depends on the curvature of the M-profile at the minimum point at z = Zs:
where Θ is a coefficient determined by the imaginary part of the propagation con- stanttn.
The above estimates are qualitative in nature. To calculate the field we need to solve the boundary problem (26), (27).
In a stratified troposphere, calculation of the electromagnetic field beyond the horizon is truncated to the determination of the complex propagation constants En and the height-gain functions χn(z) of the normal waves. The total field is then composed of a superposition of the normal waves, which converge in a shadow region. The analytical solution to the problem is known only for a few etalon problems with a limited number of selected M-profiles. In the general case the modal formalism will require a numerical solution of the characteristic equation for propagation constants and, in many cases, numerical integration of the differential equations for the height-gain functions.
This article is arranged as follows: First in Section “Some Results of Propagation Measurements and Comparison with Theory”, we discuss some results of the propagation measurements and comparison with the prediction based on existing propagation models in order to highlight the current status of the theory and the remaining problems in modelling the propagation phenomena. Then we introduce the perturbation theory for the normal waves and propagation constants in Section “Perturbation Theory for the Spectrum of Normal Waves in a Stratified Troposphere”. The perturbation theory provides a means of analytic study of the spectrum of the propagation constants and the height-gain functions with small variations of the M-profile from the etalon profile for which the solution is known.
Section 4.3*** is dedicated to the determination of the spectrum of normal waves in an evaporation duct. The solution is limited to the case of a stratified troposphere that is commonly used in existing radio coverage prediction systems.
In Section 4.4*** we will study the impact of random fluctuations in refractive index on propagation inside an evaporation duct. The results suggest that the impact may be significant in the frequency range 10 GHz and above.
Section 4.5*** deals with the height gain structure of the field inside and outside the evaporation duct for the case of scattering on the turbulent fluctuations in the refractive index. It is shown that scattering on random inhomegeneities of the refractive index leads to a smoother height dependence of the field beyond the horizon compared with duct theory. The results seem to be closer to observations.
Some Results of Propagation Measurements and Comparison with Theory
Systematic studies of radiowave propagation over the sea surface started in the 1940s, driven by the development of radar, and are still underway due both to the importance of the problem and advances in technology and applications. Several comprehensive programs of radio-meteorological measurements performed in the 1940s and 1950s still provide a benchmark reference for further studies due to the exceptionally high quality of the obtained results.
Here we briefly review the major results. The propagation measurements were performed at two frequencies 3 GHz and 10 GHz over a sea path extending to several hundred kilometres and were supported by shipboard meteorological measurements of the temperature, water vapor and air pressure. The restored M-profile was then averaged over several measurements along the propagation path and that averaged profile was approximated by a linear-exponential M-profile (Peceris’s model) in order to perform a theoretical calculation of the propagation loss and to estimate the attenuation rate of the signal beyond the horizon. A sample of the comparison of the experimental attenuation rates with Peceris’s duct model is presented in Figure 1 for 3 and 10 GHz measurements.

As observed from Figure 1, the measured attenuation rates of the 3 GHz signal tend persistently to be less than predicted from the linear-exponential model of an evaporation duct. In contrast, the measured attenuation rate significantly exceeded the theoretical estimates. The conclusion drawn in 1949 was that there was a significant quantitative discrepancy between the theory and the measurements.
A second revision of the reference data was performed in the 80s. A new advanced computer-based model of an evaporation duct was implemented to calculate the propagation constants and height gain functions. For comparison of the theory and measurements it was decided to limit to the distances to less than 150 km. The analysis of the measured data led to the conclusion that the evaporation ducts restored at the time of measurement may explain the behavior of the signal in that sub-range, while at larger distances the observed signal levels might be explained by either single-scattering theory in the case of low level signals or the presence of the elevated M-inversion. The last was difficult to analyse since there were no adequate radiosound measurements performed at the time of the radio measurements.
Figure 2 shows the result of a comparison of the same measured data as for Figure 1 at 3 GHz with a computer-based calculation for bilinear approximation of the averaged M-profile. As observed, the measured attenuation rates are still less than those predicted from the evaporation duct model. Major discrepancies are observed for lower duct heights, in conditions of unstable stratification.

One of the measured samples is shown in Figure 3, where the restored M-profile reveals an evaporation duct with parameters: Zs = 6 m, and ΔM = 2 N-units. Such discrepancy tends to be persistent in other measurements and cannot just be explained by unavoidable errors in measurements of the meteorological data.

In principle, the bilinear model tends to underestimate the attenuation rate compared with the more realistic linear-logarithmic model of the evaporation duct. This leads to the conclusion that the other models of propagation may provide a somewhat intermediate situation when the evaporation duct is not sufficiently strong to support the trapping mechanism by itself, however, negative gradients of the refractivity provide significant enhancement of the propagation beyond the horizon. One of the unusual but possible mechanisms is studied in Chapter 6***.
It should be also noted that, in the presence of a strong evaporation duct with heights exceeding 10 m, the waveguide mechanism is more pronounced at a frequency of 3 GHz, and comparison of the evaporation duct model with measurements is rather satisfactory. A typical example is shown in Figure 4 for an evaporation duct with height Zs = 14 m, and ΔM = 7 N-units. The heights of the transmitting and receiving antennas are 18 m and 6 m, respectively.

Figure 5 shows the comparison of the revised attenuation rates for frequency 10 GHz. As observed, the spread of the data points relative to the bisector is wider than in Figure 2 for the 3 GHz data.

In most cases the measured attenuation rates are higher than predicted and sometimes the actual propagation mechanism can be questioned, for example in Figure 6. In general, the evaporation ducts with a height of over 10 m are supposed to provide a waveguide mechanism at 10 GHz with attenuation rate of several hundredth dB km–1 (chiefly due to a finite impedance of the sea surface). On the other hand, the attenuation rates measured during the experiment were never less than 0,2 dB km–1 at 10 GHz, even in the presence of a very strong and stable with distance evaporation duct.

Figure 7 shows the results obtained when the height of the receiving antennas was changed during the experiment. The evaporation duct structure was stable over the distance and the assumption of the uniformity of the duct in a horizontal plane was quite reasonable during the experiment. The evaporation duct with parameters Zs = 14 m, and ΔM = 7 N-units forms two trapped modes at frequency 10 GHz.

The interference between these modes is responsible for the ripples in the theoretical curve for the received signal at a height of 12 m. At a height of 6 m the interference is not pronounced since the second mode has a minimum at this height. A receiver placed above the evaporation duct at a height of 18 m is supposed to receive the em field by penetration through a potential barrier, or, in other words, by a “leak” of the trapped modes. Theoretically predicted signal levels at 18 and 6 m should then differ by 20 dB, as shown in Figure 7. On the other hand, the measured data do not reveal any significant difference with the height of the receiving antenna. This situation is not unique and is observed in other experiments. One possible explanation of the above phenomena is provided in Section 4.5***.
We may conclude that significant development in modelling the evaporation duct structure in the lower part of the marine boundary layer resulted in the availability of computer-based tools capable of reasonably good prediction of the shipboard radar coverage in nearly real time over a wide range of frequency bands from VHF to EHF.
The statistical assessment of the evaporation duct contribution also demonstrates unexpectedly high signal levels beyond the horizon in the case of weak ducting conditions at frequencies 0,6 GHz and 3 GHz, where the evaporation ducts with heights less than 10 m are too weak to ensure a single trapped mode. Nonetheless, the received signal tends to be higher than predicted, as shown in Figure 8.

It was suggested that the observed high level signals might be caused by propagation mechanisms other than an evaporation duct. One of the conclusions emphasised is that the height dependence of the received field is not as significant as may follow from the evaporation duct theory and, basically, high altitude antennas are preferable. In terms of the optimal frequency the results of the measurements suggest the 10–20 GHz band is an optimal band where the evaporation duct makes a strong impact without being counteracted by gaseous absorption at higher frequencies.
A clear demonstration of the advances in radar coverage prediction and range detection enhancement due to an evaporation duct for 18 GHz radar measurements. Similar results for propagation experiments in a range from 3 to 94 GHz. The predicted values of the path loss are depicted as a curve in Figure 9.

The ripples for the higher duct’s height are caused by the presence of several trapped modes and corresponding restructuring of the height gain function of the total trapped field in the waveguide. While overall agreement between prediction and measurements is very good the prediction tends to underestimate the path loss for lower duct heights. It might be noted that the scattering on the rough surface is modelled in a rather conservative way using the Kirchhoff approximation, which is applicable to a coherent field component. Nonetheless, an additional attenuation factor, like scattering at random fluctuations in the refractive index, might be taken into account for possible improvement in the prediction model. It may also be pointed out that the spread of the data points in Figure 9 tends to be more close and symmetrical relative to the prediction curve, which may again suggest the importance of the full treatment of the scattering mechanism. We recall these arguments later in Section 4.3*** where the scattering mechanism is described.
Another experiment in the millimetre wave range for 94 GHz propagation. The results obtained demonstrated the validity of the evaporation duct mechanism in this frequency band and gave good agreement between the measurements and the theory of the evaporation duct, with some underestimation of predicted path loss with an average value of error of 10 dB. It might be noted that while the prediction model has made use of the Kirchhoff theory of the scattering at a rough sea surface, the model did not take into account scattering on turbulent fluctuations in the refractive index.
Perturbation Theory for the Spectrum of Normal Waves in a Stratified Troposphere
The presence of stratified inhomogeneities of the refractive index results in a considerable change in the field structure in the region of geometric shadow. A complete solution to the problem of wave propagation in a stratified troposphere, in which the refractive index depends only on the height h above the earth’s surface, was obtained by Fock. He showed that for large values of the parameter m = (ka/2)1/3 ≫ 1, here k is the wave number and a is the earth’s radius, the problem can be simplified by transition from a spherically stratified to a plane-stratified medium. Such a simplification is provided through the use of the Understanding Parabolic Approximation in Wave Propagation: Analytical Methods and Applicationsparabolic approximation and introduction of the modified permittivity:
Even in this simplified formulation an analytical solution can be obtained only for a few standard problems with a limited set of permittivity profiles εm(h). In the general case the analysis of the propagation conditions comes down to a numerical solution for an eigenvalue problem with complex propagation constants. In several cases the function εm(h) differs little from the “etalon” function, as for instance in the case of weak refraction. Therefore, some need arises to construct a perturbation theory for “open” system for which the spectrum of eigenvalues (propagation constants is complex and eigenfunctions (height factors for the Wave Field Fluctuations in Random Media over a Boundary Interfacewave field χn(h)) grow exponentially with height, χn(h) → ∞, h → ∞. A similar problem arises in quantum mechanics with the study of the decay of quasi-stationary states. Perturbation theory for that case was developed by Zeldovitch, although the essential assumption was made on the finiteness of the potential (the analogue of εm(h)) at (h) → ∞. Such a condition is not satisfied in the problem of wave propagation in a stratified troposphere:
Here we describe a generalisation of the perturbation theory to the case of potentials unlimited at (h) → ∞.
Problem Formulation
The field attenuation factor V for a point source in the spherically stratified medium can be represented by superposition of normal waves:
where we introduce dimensionless coordinates:
- k0 – is the wave number in a vacuum, and;
- D – is the distance along the earth’s surface from the source of radiation.
The transmitter (point source) is situated at height h0 and the receiver at height h, respectively. The height functions χ(γ, tn) are governed by the equation:
and satisfy the following boundary conditions:
Equations (4) and (5) determine a discrete spectrum of propagation constants tn. We represent the modified refractive index U(γ) = m2(ε(γ)-1)+γ, as shown in Figure 10, in the form U(γ) = U0(γ)+δU(γ), where U0(γ) is the unperturbed index of refraction for which the solution to the boundary value problem Eqs. (4) and (5) is known. We shall treat δU(γ) as the perturbations to the etalon profile U0(γ).

We introduce the logarithmic derivative:
through which the height function is expressed as:
For zn(γ) we can obtain the equation:
We seek zn(γ) and tn in the form:
The solution to Eq. (10) has the form:
Substituting Eq. (12) into Eq. (11), we obtain the equation for correction to the propagation constant:
where:
Equation (13) differs from the usual equation of perturbation theory in two aspects:
- first, there is aunder the integration instead of the square of the absolute value |χn(γ)|2 and;
- second, instead of the norm, which does not exist in the present case, there is the finite term N equal to:
In the case of the troposphere the variation of the modified refractive index becomes linear at large enough heights and Eq. (14) takes the form:
In the case of normal refraction:
Equation (14) is then reduced to:
Calculating the integral appearing in Eq. (15):
and using the boundary condition (5), we obtain:
Linear Distortion
Let us consider a linear distortion:
where 1+μ3 is the gradient of the modified refractive index in the height range 0 < γ ≤ H as shown in Figure 10. The equation for the spectrum tn of the normal modes was Basics of Fock’s Theory – the Hartree-Fock Approach to Evaporation Ductsderived by Fock, its solution cannot be obtained analytically and requires application of numerical methods. Application of perturbation theory in form (17) allows investigation of the tn-spectrum in a considerably simple way.
When the stronger inequality H2 τn ≪ 1 applies we can obtain from Eq. (20):
As observed from Eq. (21), increase in the height H or depth (1+μ3)H of inversion in the refractive index leads to a decrease in the real and imaginary components of the propagation constant tn, that corresponds to an increase in the phase velocity of the normal wave and a decrease in its attenuation.
The applicability of the perturbation theory is limited by a small variation to the propagation constant:
Inequality (22) in conjunction with H2τn ≪ 1 yields the limits of applicability of Eq. (21):
and, for Eq. (20):
Let us consider another limiting case, when the height of distortion is large H ≫ τn. We represent the integral in Eq. (17) as the sum of integrals over the arc of radius H in the sector of angles (0, π/3) and over the ray arg γ = π/3:
In the integral I1 we make a change in the variable γ = τejπ/3 and use the property of the Airy function w1(xejπ/3) = 2ejπ/3 ν(-x):
The upper limit of integration in Eq. (25) can be extended to infinity since the Airy function ν(τ-τn) decreases exponentially for τ > τn. Then, using Eq. (16) we obtain:
where:
The main contribution to the integral (27) comes from the region of small φ where the only linear terms can be retained in the expansion of the integrand. Then we obtain:
Substituting Eqs. (28) and (26) into Eq. (17) we obtain:
As observed from Eqs. (29) and (21), the correction factor dtn grows exponentially with number n since the imaginary part of the propagation constant is Im tn~n3/2, and hence the inequality (22) ceases to be satisfied, starting with certain n. It can be shown, however, that this result is a consequence of the “non-physical” nature of the broken linear profile. For smooth functions U(γ) that actually describe realistic behavior of the refraction index, the dtn-dependence on a wave-number n has a fundamentally different character (see, for instance, Eq. (40) below).
In fact, performing integration by parts in Eq. (17), we can write:
where:
If δU(γ) is a function continuous along y altogether with its derivatives and which decreases sufficiently at infinity, then integration in Eq. (30) will result in values of δtn which decrease for large n.
In the case where δU(γ) is abruptly reduced to zero at γ = H (a discontinuity of the first kind), a singularity of a delta-function type appears in the integrand in Eq. (30) at this point, and at large tn it make the main contribution to the integral:
This result can be generalized to the case of discontinuity of the kth derivative:
Smooth Distortion
Let the function δU(γ) be such that F(s) has no singularities in the right half plane s or on the imaginary axis. Then the integration contour over ds can be shifted into the left half plane, while the contour in the integral over dγ can be deformed into a ray with arg(γ) = π/3. Substituting the asymptotic form (19) into Eq. (33) we obtain:
In line with the above assumptions, Re s < 0, so the second term in Eq. (32) vanishes. Using an inverse Laplace transformation, we obtain:
It should be noted that Eq. (35) could be used when the convergence sector of the function δU(γ) is wider than π/3.
Let us consider the application of Eq. (35) to the example of perturbation given by δU(γ) = U0 exp (-γ/H). We set τn ≫ M1 and H ≪ 1. Substituting the asymptotic form ν(γ-τn), we obtain:
When H2 τn ≪ 1 Eq. (37) converges to a simple form:
In the opposite case, H2 τn ≫ 1, we obtain:
Height Function
Linear-Logarithmic Profile at Heights Close to the Sea Surface
Consider another specific kind of perturbation to the refractivity profile, which can be applied to the estimation of the impact from large gradients of humidity in the immediate vicinity of the sea surface. Refractometer measurements or measurements of the meteorological parameters are difficult to perform at these heights because of the roughness of the sea surface. On the other hand, this layer can be characterised by extremely rapid changes in refractivity due to the gradients of humidity. As follows from hydrodynamic theory of the evaporation duct, the modified refractivity profile in the layer close to the sea surface can be modelled by a linear-logarithmic function:
where Ys = kZs/m and Zs is the evaporation duct’s height, γr = kzr/m, and zr is the height above the sea surface below which it is not feasible to make the measurements. Normally zr is associated with a sea roughness parameter. Now we consider the impact from the gradient of humidity in a layer below γr as a perturbation δU to a linear profile:
In fact, we may approximate the initial profile (41) with a linear segment in the interval 0 < γ < γr as discussed later in Section 4.3***, and regard the combined profile as the “true” profile while considering the logarithmic part as the perturbation.
We may now evaluate the reference height zr below which the exact shape of the profile does not make a difference to the calculation of the field attenuation beyond the horizon. The impact perturbation is negligible if the variation in attenuation of the first mode due to perturbation (42) is small, i.e:
In the range of distances x ~ 100 km for frequencies of the order of 10 GHz, we may see that the above inequality is satisfied with zr ≤ 2 m. Therefore, we may state that the refractivity measurements can be performed from heights above zr = 2 m and the M-profile close to the surface at heights below zr can be approximated by a linear function or any other function approaching a finite value at z = 0.