.
Our site needs your help!
Site categories

Atmospheric Boundary Layer and Basics of the Propagation Mechanisms

Join Our Telegram (Seaman Community)

Explore the Atmospheric Boundary Layer and its impact on signal propagation in this comprehensive article. Discover the Standard Model of the Troposphere, non-standard mechanisms such as Evaporation Ducts and Elevated M-inversions, and the influence of dielectric permittivity fluctuations. Gain valuable insights into how these factors affect communication and environmental conditions.

The troposphere is the lowest region of the atmosphere, about 6 km high at the poles and about 18 km high at the equator. In this book we study radio wave propagation along the ocean and can reasonably assume that all processes of propagation occur in a lower region of the earth’s troposphere. That lower region and the atmospheric conditions are of most importance for the subject under study here.

From the perspective of radio communications/propagation we limit our objective to an investigation of the impact of the atmospheric structure on the characteristics of the radio signal propagating through the atmospheric turbulence. All known methods of solutions to a similar problem are based on the separation of the spacetime scales of the variations in both the refractive index n and electromagnetic field
E, H
in two domains, described in terms of deterministic and stochastic methods.

Intuition suggests that the spectrum of turbulent variations in the refractive index n will have the energy of its fluctuations confined to a limited space–time domain or, at least, have a clear minimum and, desirably, a gap spread over a significant interval in the time–space domain. It is apparent that the horizontal scales of variations in refractive index larger than the length of the radio propagation path have no immediate effect on the characteristics of the radio signal and rather affect its long term variations over the permanent path. This comes down to an upper boundary of the spatial variations of refractive index in a horizontal plane of about 100 km. The vertical irregularities are of most importance since they are responsible for the refraction and scattering of the radio waves in troposphere. However, there are some natural limitations on the region of the troposphere which might be of interest in its impact on radio propagation. The troposphere is naturally divided into two regions: the lower part of the troposphere, commonly called an atmospheric boundary layer, and the area above, called clear atmosphere.

The electric properties of the troposphere can be characterised by the dielectric permittivity ε or the refractive index
n=ε.
The numerical value of the non-dimensional parameter n is pretty close to unity, however even a relatively small deviation of the refractive index from unity may have significant impact on radio wave propagation. Therefore, common practice is to use another definition of the refractive index N = (n–1) × 106 instead of n, measurable in so-called N-units.

The refractive index N, also called the refractivity, has the following relationships with atmospheric pressure p, temperature T and humidity, the mass-fraction of the water vapor, q, in the air:

N=AN·pT1+BNqT           Eq. 1

where AN = 77,6 N-units × K hPa–1, BN = 7733 K–1. The components p, T and q are random functions of the coordinates and time. The stochastic behavior of the meteorological parameters p, T, q and, therefore, the refractive index N is caused by atmospheric turbulence.

There are several reasons to separate the region of the first 1–2 km of atmosphere over the earth’s surface, called the atmospheric boundary layer, ABL. The upper boundary of the ABL is seen as the height at which the atmospheric wind changes direction due to a combined effect of the friction and Carioles force.

Among those reasons are:

  • The interaction of the earth’s surface and atmosphere is especially pronounced in this region.
  • The meteorological parameters such as temperature, humidity and wind speed experience daily variations in this region due to apparent cyclic variations in the sun’s radiation due to the earth’s rotation.
  • The ABL can be regarded as an area constantly filled with atmospheric turbulence. This is quite opposite to the atmospheric layer above the ABL, the so-called region of clear atmosphere, where turbulence is present only in isolated spots.
  • The border between the clear atmosphere and the ABL is clearly pronounced with sharp variations in all meteorological parameters, as illustrated in Figure 1.

The spectrum of turbulent fluctuations in the atmospheric boundary layer is extremely wide: the linear scales of the variations range from a few millimetres to the size of the earth’s equator, the time scales from tens of milliseconds to one year.

Atmospheric pressure
Fig. 1 Meteorological parameters in the atmospheric boundary layer as function of height (atmospheric pressure, p):
Humidity (a), temperature (b) and wind speed (c). All parameters experience sharp variations at the upper boundary of the atmospheric boundary layer

Studies of the energy spectrum of the turbulent fluctuations of the meteorological parameters (temperature, humidity, pressure and wind speed) have shown that the energy spectrum reveals three distinct regions:

  • Large scale quasi-two-dimensional fluctuations in a range of frequencies from 10–6–10–4Hz, the meso-meteorological minimum with low intensity of the fluctuations in the range 10–3–10–4Hz and;
  • A small scale three-dimensional fluctuation region with frequencies above 10–3Hz.

Figure 2 shows the energy spectrum of fluctuations of the horizontal component of the wind speed in the atmosphere, where the ordinate corresponds to the product of the spectrum density S(ω) and the cyclic frequency ω of the variations in one of the meteorological components, and the abscissa corresponds to the frequency ω = 2π/T, T being a period of variations.

Energy spectrum
Fig. 2 Energy spectrum of the fluctuations in the wind speed in the atmospheric boundary layer

As observed in Figure 2, there are two major extremes of the function ω · S(ω) the high frequency maximum corresponds to a linear scale of turbulence of the order of tens and hundreds of meters, the low frequency maximum has a time scale of 5–10 days which is caused by synoptic variations (cyclones and anti-cyclones), the respective horizontal scale is thousands of kilometres and the vertical scale is of the order of 10 km. There is also an extended minimum in the spectrum ω · S(ω) that corresponds to the fluctuations with respective horizontal scales from 1 to 500 km and is called the meso-pause. The region of low frequency variation is called the macrorange while the region to the left of the meso-pause (high frequency variations) is called the micro-range.

The nature of the atmospheric turbulence is different in these two regions: in the macro-range the synoptic processes can be regarded as two-dimensional variations, while in the micro-range, with scales up to a hundred meters, the turbulence is three-dimensional and locally uniform. The mezo-pause is a transition region where a combined mechanism is observed. It is important to notice that by describing small-scale three-dimensional fluctuations in the micro-range region one can use Taylor’s hypothesis of “frozen turbulence” which allows a transformation from time to space-fluctuation scales by means of:

L = 2pv/f,

where:

  • L is the spatial scale of the irregularities;
  • f is the frequency of time variations in the refractive index, and;
  • v is the mean speed of the incident flow.
The basic conclusion that follows from the above observations is that, to some extent, the refractive index N and the dielectric permittivity ε can be presented as a sum of a slow varying component
ε0r
regarded as a quasi-deterministic function of the coordinates
r={x, y, z}
and the random component
δε(r).
As observed from Figure 2, the quasi-deterministic component
ε0(r)
still varies in the horizontal plane and the energy of variations in the meso-pause minimum is not negligible. However, these variations have less impact on radio wave propagation than either variations of
δε(r).
in the micro-range or over-the-height variations in the “deterministic” component which may be responsible for a ducting in the troposphere.
Mathematically, the problem of radio wave propagation in a randomly inhomogeneous medium comes down to solving a stochastic wave equation with dielectric permittivity
ε(r, t)
which is a random function of coordinates and time. In many cases the process of propagation of a monochromatic wave in the troposphere can be considered in a quasi-steady state approximation, i. e. “frozen” in time.

It is then convenient to represent the dielectric permittivity in the form:

ε0(r)ε0(z)+δε(r),

where:

ε0(z)=ε0(r).
The angular brackets denote averaging over the ensemble of the realisations of
ε(r).
In fact, the mean characteristic of the tropospheric dielectric permittivity
ε(r)
is commonly understood as a large-scale structure homogeneous in the horizontal plane and practically non-varying over the time over which the signal measurements have been performed and then, as a result of mathematical idealisation,
ε(r)=ε0(z).

In radio-meteorology mean characteristics of the meteoparameters are usually understood to be the values obtained by averaging over a 30 min interval, i. e. averaging is performed over a frequency interval the lower limit of which is positioned within the limits of the meso-meteorological minimum.

The average characteristic obtained in this way is usually a function of the height z above the surface (sea, ground) and varies slowly with the horizontal coordinates and time. Assuming ergodicity and Taylor’s hypothesis, such averaging over a time interval corresponds to the averaging over the ensemble of the realisations of
εr, t.
As a compromise in analytical studies of radio wave propagation, a common approach is to neglect the residual variations in
εr
over the horizontal coordinates, i. e. to regard the
ε0rε0(z).
This assumption results in the introduction of the traditional model of a stratified atmosphere, in which the average values of refractive index N vary over the vertical coordinate z, the height above the ground. This traditional model provides some basis for a classification of the radio wave propagation mechanisms, in particular, a separation of the propagation into two classes: standard and non-standard. The following section below provide a brief analysis of the standard and non-standard models, while other section deals with a statistical model for a random component of the refractive index.

Standard Model of the Troposphere

The standard mechanism of radio wave propagation is classified under the condition where the average vertical gradients of the refractive index
γN=dN/dz
are close to the value
γNst=39 Nunits km1.

Such conditions of refractivity constitute a model of standard linear atmosphere defined as:

N=28939·z               Eq. 2

and are applicable for heights less than 2 km. Let us consider this model in detail involving a geometrical optic presentation for wave propagation. Let us define “ray” as a normal to a wavefront propagating through the medium with varying refractivity n(z). As is known, the ray bends in such a medium and the bending is defined by Snell’s law.

Shell law
Fig. 3 Illustration of Snell’s law

Introducing the horizontally stratified medium in terms of the set of thin layers with value of refractivity ni, i = 0, 1, …, such as illustrated in Figure 3, Snell’s law can be written as:

ni cos ϕi=const           Eq. 3

where ϕi is a sliding angle. Introducing the differentials of the ray direction dz, dS in the ith layer and differentiating both sides of Eq. (3) with respect to S:

cos ϕidnidSni sin ϕidϕidS=0.             Eq. 4

Substituting dS = -dz/sin ɸi, we obtain:

dϕidS=cos ϕi(dn/dz)ni.                  Eq. 5

The radius of curvature at any point, Ri = dS/dϕi, and using Eq. (5) it results that:

Ri=nicos ϕi1(dn/dz).               Eq. 6

For the standard atmosphere with dn/dz = 39×10–6 km–1, the radius of curvature is given by:

Ri=25 000nicos ϕi               Eq. 7

If the launch angle ɸi is close to the horizontal, the ratio ni/cos ɸi ≈ 1 and a propagation path can be described as a circle of radius R = 25 000 km, Figure 4(a). By comparison, the radius of the earth’s curvature is α = 6 370 km, and 1/α = 157×10–6km–1. When the curvatures of both the propagation path and the earth are reduced by 39×10–6, as in Figure 4(b), the propagation path has an effective curvature of zero (which is a straight line) and that hypothetical earth has an effective curvature of (157-39)×10–6km–1 = 118×10–6km–1.

Effective earth radius
Fig. 4 Introduction of the “effective” radius of the earth:
a) Ray refraction in a “normal” atmosphere with “true” earth radius α = 6 370 km. b) Effective ray refraction in case when the difference in curvatures of both ray and earth surface in (a) is compensated by introduction of the modified earth radius αe = 8 500 km

The equivalent radius of the sphere αe can therefore be defined as:

1αe=1α1R=1αγN·106=1118·106 km=8 500 km,               Eq. 8

approximately, αe = 4/3α.

The modified refractive index can be defined as follows:

nm(z)=n(z)+1/α.                Eq. 9

It is apparent that the second term in Eq. (9) is a compensation for the earth’s curvature. The modified refractivity M(z) is then given by:

Mz=nmz1·106=Nz+0,157z.              Eq. 10

As observed from Eq. (10) for the case of standard refraction, when γN = dN/dz = -39 N-units km-1, γM = dM/dz = 118 N-units km-1. As explained in Section 2.1***, modified refractivity comes logically from the parabolic approximation to the wave equations and effectively results in the introduction of the flat earth and additional curvature of the “rays” associated with radio waves propagating at low-angles along the earth’s surface, as shown schematically in Figure 5.

Refractivity on the ray
Figure 5 Impact from modified refractivity on the ray trajectories from a point source

This traditional linear model of the troposphere provides some means of classification of the propagation mechanisms based on the value of the gradient of the refractivity in the lower troposphere, both conventional γN = dN/dz and modified refractivity γM = dM/dz. A very broad classification is concerned with the separation of the propagation into two classes: standard and non-standard. Figure 6 represents schematically the ray traces for three basic sub-classes of standard mechanisms of propagation in the troposphere:

  • Standard refraction, when the vertical gradient of refractivity is very close to
    γNs=39 km1, γMs=118 km1.
  • Sub-refraction, when γN > –39 N-units km–1, γM > 118 N-units km–1.
  • Super-refraction, when γN < –39 N-units km–1, γM < 118 N-units km–1.

As seen from Figure 6 the sub-refraction and super-refraction are associated with bending the waves outwards and inwards to the earth’s surface respectively.

Different refractions
Fig. 6 Classification of the propagation mechanisms
The standard propagation conditions can be characterised by the presence of two regions: the so-called line-of-sight (LOS) region, where the radio signal in the receiver has contributions from the direct wave and the wave reflected from the earth’s surface; and the shadow region where the radio signal arrives as a result of diffraction over the earth surface. These two regions are separated at the range of horizon which can be defined as
2αe·z
, where the parameter z is the height of the observation point above the earth’s surface. It should be noted that for a standard refraction
2αe=4,12·103 m1/2
compared with the value for optics of 3,83·103m1/2. For frequencies above 20 GHz the standard mechanisms include attenuation of the radio waves due to absorption by atmospheric gases (water vapor and oxygen).

Non-standard mechanisms of radio wave propagation are related to some anomalies in the vertical distribution of the refractive index, namely the gradients cM become negative. These anomalies are in turn caused by some extremes in the vertical profiles of either the humidity or temperature.

Under such conditions the waves can be trapped inside a localised region associated with such negative gradients of modified refractivity, this effect is called “ducting” and the radio wave propagation mechanism is somewhat similar to propagation in a waveguide.

Non-standard Mechanisms of Propagation

The term M-inversion is used to characterize the area of negative gradients of the M-profile, where dM/dz < 0. When M-inversions exist over some height interval, a radio-waveguide can be formed. The upper boundary of the waveguide is the height H at which dM/dz = 0, while the lower boundary zmin is given by M(zmin) = M(H), see Figure 9 later. The waveguide is elevated if zmin > 0. If zmin = 0 or the equation M(zmin) = M(H) has no solution, then the waveguide is classified as a surface wave-guide, Figure 7.

Graph of evaporation
Fig. 7 Evaporation duct, M-profile

Most often M-inversions are formed in the lower and upper regions of the atmospheric boundary layer because of the local interaction between the turbulent air mass and either the surface or the free atmosphere above the ABL.

Evaporation Duct

One of the most important types of radio waveguide is the evaporation duct above the water surface which has a nearly permanent presence in the subtropical regions of the Pacific Ocean. It is formed as a result of the air in contact with the sea surface becoming saturated by water vapor at the surface temperature.

The water vapor content decreases logarithmically with height and the resulting M-profile has the form:

M(z)=N(zr)Hα lnzzr+αz              Eq. 11

where:

  • zr is the roughness parameter of the sea surface and;
  • α = 106/α = 0,157 m–1 and α = 6 370 km, the earth’s radius.

The M-profile (11) is can be constructed from the data of a standard hydrometeorological measurement. It describes quite well the height dependence of the M-profile under stable and neutral conditions of stratification in turbulence. The studies performed show that the thickness of the duct (inversion height H) is a function of the stability parameter of the turbulence at the water level layer and is functionally related to the inversion depth ΔM = M(zr) – M(H), Figure 7.

It should be noted that no analytical solution to the wave equation exists for the M-profile in form (11). While a numerical solution can be adopted with some modifications to Eq. (11), in some applications it might be preferable to replace Eq. (11) by a profile for which the analytical solution to a wave equation is known.

In particular, a common approach is to use a piecewise linear M-profile given by:

M(z)=M(H)+α(zH),                   z>HM(H)+M(Hz)/(Hzr), zH            Eq. 12

Comparison of the results obtained with both the piecewise linear approximation and the continuous M-profile (12) is presented in Section 4.1***. It has shown a satisfactory agreement of the results of calculations for the lower order modes. In fact, the discrepancies observed for the high order modes in terms of their propagation constants and therefore the height functions are not really important because those modes are significant in the transition area between the shadow and LOS regions, where, in turn, the mode presentation is hard to utilise because of the very slow convergence of the modal series. On the other hand, it was observed that the errors introduced by replacement of the continuous profile by a piecewise approximation for the lower modes do not exceed the errors produced by uncertainties in determining the values of the refractive index from data of standard hydro-meteorological measurements. The overall conclusion therefore is that the approximation (12) is a valid practical approach to analytically evaluate the propagation mechanism in the evaporation duct. In computer based prediction systems, where input radio-meteorological information can be collected in nearly real-time with substantial precision, the hydro-dynamical model of the evaporation duct Eq. (11) can be easily implemented. The common problem of calculating the propagation constants of the trapped (or significant) modes can be reduced to pre-calculated look-up tables and interpolation.

Trajectory of the rays
Fig. 8 Ray trajectories in evaporation duct

Figure 8 shows a schematic ray bending in the evaporation duct, which allows over-the-horizon radio wave propagation.

Elevated M-inversion

The most frequent cause of abnormally high signal level far beyond the horizon (200–500 km) is associated with the presence of elevated M-inversions. Those are formed in the subtropical region above the ocean as a result of powerful motion downwards in the center of an anticyclone and are usually called depression layers. The elevated M-inversions, also called elevated waveguides or inversion layers, are characterized by lower height of inversion Hi, the thickness of the layer with negative M-profile (M-inversion) ΔH = Hu-Hi and the depth of inversion, the so-called M-deficit ΔMi = M(Hi)-M(Hu), as shown in Figure 9.

M-inversion
Fig. 9 Elevated M-inversion:
a) Elevated tropospheric duct in the height interval Hu > z > Hb; b) Surface based tropospheric duct in the height interval Hu > z > 0

Figure 10 shows schematic ray trajectories in an elevated duct when the source is placed inside its boundaries.

Traces of the rays
Fig. 10 Ray traces in the elevated duct

The observations indicate that characteristic values for the depression layers are as follows: Hi ≈ 200-1 500 m, ΔH ≈ 100-300 m, and ΔMi ≈ 10-30 units. The depression layers are normally observed as stable atmospheric formations existing over a wide area of the ocean surface. The height of inversion usually coincides with the height of the atmospheric boundary layer. The elevated M-inversion normally separates a convective boundary layer from the free atmosphere above, and in the direct vicinity of M-inversion turbulence alternates with laminar flow. Due to a penetration of convective elements from the lower atmospheric boundary layer the local temperature and gradients of the wind speed increase and as a result Kelvin–Helm-holtz waves develop. When these waves collapse they generate an intense turbulent structure.

Characteristic scales of the turbulent inhomogeneities in the depression layer can be estimated as follows: horizontal scale L||~500–1 000 m, vertical scale Lz ~100 m.

There is another important type of elevated M-inversion which is often formed in coastal regions due to the formation of a surface duct as a result of heating of the sea surface at sunrise and elevation of the surface M-inversion at night time. The elevated M-inversion is then unstable due to a developed convection.

Random Component of Dielectric Permittivity

A comprehensive theory of random fields and methods of their study are provided in Refs. Below we summarize the basic assumptions to be used further in this category to describe the random field of dielectric permittivity
ε(r, t).
When the random field
δε(r, t)
is stationary and uniform in time and space respectively, a two-point statistical characteristic of this field can be described via a space–time correlation function:
B(ρ, τ)=δεr+ρ, t+τδεr, t.                Eq. 13
It is assumed that < δε > = 0, and, as a consequence of uniformity and stationarity,
δε(r, t)
is invariant to an offset in the space and time domain, and its correlation function (13) depends only on the difference in space
ρ
and time τ between two points of observation.
The space–time correlation function
Bρ, τ
is related to the spectrum density G(ω, k) of the field
δεr, t
via a Fourier transform:
Bρ, τ=dω  dkGω, k·ejkρωτ.                        Eq. 14
The function G(ω, k) is also called a
ωk
density, where
k
is the wave vector of fluctuation in δε.

The reverse transform is given by:

Gεω, k=12π3d3ρdτBρ, τejkρωτ.                  Eq. 15
Both the spatial spectrum Φε and the frequency (energy) S(ω) spectrum of the fluctuations δε can be expressed via
ωk
density:
Φεk=d3ρBρ, 0·ejkρ,                    Eq. 16
Sεω=dkGω, k=dτB0, τ·ejwτ.                    Eq. 17
Consider the case when all time variations in
δεr, t
are caused by transfer of the spatial inhomogeneties
δε~r
with constant velocity
ν
, then
δεr, t=δε~rνt.
If a motionless uniform random field
δε~r
is described by a spatial correlation function:
bε~ρ=δε~r+ρ

then:

Bτ, ρ=bε~ρντ.                Eq. 18
Substituting Eq. (18) into Eq. (15) we can find a relationship between
ωk
the density of the “frozen” field
δεr, t
and the spatial spectrum density
Φε~
of the still field
δε~r
:
Gω, k=Φε~kδωkν.                  Eq. 19
Using Eqs. (16) and (19) we find that both the spatial spectrum of the frozen field
δεr, t
and that of the motionless field
δεr
are the same.
Φεk=Φε~k.                    Eq. 20
This is an important relationship since all further studies in this book are based on the hypothesis of the “frozen” turbulent fluctuations of dielectric permittivity
δεr, t
and all calculations are performed with a motionless field of
δεr
instead of
δεr, t
and the symbol ~ can be omitted for
δεr
.

Locally Uniform Fluctuations

The random field of dielectric permittivity
δεr
can be regarded as uniform only as a very rough approximation. For example, the statistical characteristics of
δεr
in the atmospheric boundary layer depend on the height above the ground and, therefore, another approach, based on the introduction of the structure function, should rather be used to describe the turbulent fluctuation of the
δεr
.
The structure function is an average of the square of the module of the increment of the fluctuations
δεr
and is defined as:
Dr1, r2=δεr1δεr22=δεr1δεr2)2.               Eq. 21
The module sign is removed in Eq. (21), given the assumption that
δε~r
has a real value. The hypothesis of the local uniformity is based on the assumption that the increment of the fluctuations of the random field de at two points
r1
and
r2
is caused chiefly by inhomogeneities of δε with a sizel less than the distance between two points of observations:
lr1r2.
While the inhomeogeneities of the large scale may exist, they produce nearly the same perturbation in the field δε in points
r1
,
r2
and, therefore, a negligible difference in the perturbation in δε. When this hypothesis holds, the field de is called “locally uniform” or, using an analogy with a random process with stationary first increment, we can call δε a random field with uniform increment. As a consequence, the structure function
Dr1, r2
depends only on the spatial increment
ρ=r1r2:
Dr1, r2=Dρ.                  Eq. 22

The important advantage of using a structure function is that it makes sense even when the correlation function does not exist, Bε(0) → ∞. However, if it does exist then the following relationship is valid:

Dερ=2Bε0Bερ.                    Eq. 23

When Bε(∞) → 0, Dε(∞) → 2Bε(0) and:

Bερ=12Dε()Dερ.                   Eq. 24

The structure function Dε is related to a spatial spectrum Φε via the relationships:

Dερ=2dkΦεk1cos kρ,                     Eq. 25
Φεk=116π3k2d3ρ·kDερ sin kρ.                  Eq. 26
When the random field is statistically isotropic, the increment of the fluctuations in δε depends only on the distance between two points
r1
and
r2
, then
Dερ=Dερ
and
Φεk=Φεk
, where
ρ=r1r2, k=k
are the modules of
ρ and k
respectively. For locally uniform and isotropic fluctuations of δε Eqs. (25) and (26) can be transformed into the following ones:
Dερ=8π01sin kρkρΦε(k)k2dk,                   Eq. 27
Φεk=14π2k30sin kρkρ·cos kρdDερdρdρ.                  Eq. 28
The Eqs. (19) and (20) are valid for a “frozen” and locally uniform random field
δεr, t.
The equation for the time–space correlation function (18) can be substituted by a similar equation for a space–time structure function:
Dερ, τ=Dερντ.                 Eq. 29
The turbulent fluctuations of dielectric permittivity
δεr
in a free atmosphere can be modelled as a locally uniform random field. The structure function Dε(ρ) of fluctuations de obeys a law of 2/3 in the Kolmogorov–Obukhov model for large enough ρ:
DερCε2ρ2/3, ρl0                   Eq. 30
where
Cε2
is the structure constant, a measure of fluctuation intensity and l0 is the internal scale of the turbulence (the inhomeogeneties with the scale less than l0 dissipate rapidly). For small ρ the structure function follows the square law:
DερCε2ρ2,    ρl0.                      Eq. 31

As observed from Eq. (30) the structure function Dε(ρ) is unlimited by a scale of fluctuations and allows the existence of turbulent bursts of any scale, which is contradictory to the applicability of the Kolmogorov–Obukhov model which is itself limited to an inertial interval of the locally isotropic turbulence.

In order to limit the external scales of the fluctuations to an inertial interval of the turbulence and to make the model consistent one can introduce a maximum scale L0 of the turbulent fluctuations and therefore a “saturation” feature in the behavior of the structure function:

Dερ=Cε2L02/3,      ρL0.                    Eq. 32
In fact, a saturation means a limited value of Dε(∞) and, therefore, a limited value of variance
σε2=Bε0
of the fluctuations δε. In this case those values are related as follows:
Cε2L02/3=2σε2.

A combined structure function for locally uniform and isotropic fluctuation in δε can be written as follows:

Dερ=Cε2l04/3ρ2, ρl0Cε2ρ2/3, l0ρL0Cε2L02/3, ρL0              Eq. 33
Within an inertial interval of locally isotropic turbulence the spatial spectrum
Φεk
is described by the Kolmogorov–Obukhov model:
Φεk=0,033 Cε2 k11/2 expk2km2.                 Eq. 34
Here
k
is the wave vector of fluctuation in δε,
k=k, Cε2
is the structure constant, a measure of the fluctuation intensity, km = 5,92/l0 and l0 is the internal scale of the turbulence.

The applicability of Kolmogorov’s model (34) for the atmospheric boundary layer over the ocean surface has been studied in Ref. The basic observations are as follows:

  • A one-dimensional spectrum of the turbulent fluctuations in the horizontal plane follows Kolmogorov’s model up to scales: L|| ~10-100 m. The spectrum interval over which the postulates of a locally isotropic turbulence are valid depends on the stability of the atmosphere and decreases greatly when the stability intensifies. In the case of very stable stratification associated with the presence of M-inversions the turbulent structure reveals significant anisotropy of the large-scale fluctuations.
  • At heights above the surface comparable with the heights of the swells, the state of the sea surface reveals a significant impact on a form of the spectrum of turbulent fluctuations. This manifest itself in the form of flares in spectral density at the frequencies corresponding to the energy-bearing frequencies of the sea waves. The amplitude of the wave perturbations decreases with height and has no practical effect on the form of the spectrum at heights of 15-20 m above the surface.
  • The intensity of the high-frequency turbulent fluctuations depends parametrically on the height above the surface. In a low-frequency domain (the buoyancy and energy intervals) one must include additional parameters in the model of spectrum: the thickness of the atmospheric boundary layer and some parameters characterizing the roughness of the sea surface.

It can be stated that despite numerous studies, a physically justifiable model for a three-dimensional spectrum of fluctuations in the index of refraction within the atmospheric boundary layer does not yet exist. The anisotropy of the spatial fluctuations in the refractive index is commonly introduced via the anisotropy parameter α = Lz/L|| in the following model of the spectrum:

Φεk=0,063 σε2 Lz L||21+k||2 L||2+kz2 Lz211/6                  Eq. 35
where the
σε2
is a variance of the fluctuations. It can be noted that in the case of turbulence Lz, L|| → ∞,
σε2
→ ∞. The generalisation of Kolmogorov’s model to anisotropic turbulence can be given as follows:
Φεk=0,033 Cε2αk2+α2kz211/6=0,033 Cεz2α5/3k2+α2kz211/6                   Eq. 36

where Cε⟂ and Cεz are structure constants in the horizontal plane and over height respectively, as observed they are not independent and are related via:

Cε=α1/3 Cεz.               Eq. 37
In practical measurements the time of averaging is always finite and so variance of the fluctuations
σε2
and scales L, Lz is observed. In fact model (35) can be trans formed into model (36) with the substitution
σε2=1,9 Cε2L2/3
, assuming that the external scale of the fluctuations in dielectric permittivity L is confined within the micro-range region of the spectrum of turbulence.

The observations suggest that the value of the anisotropy parameter α itself is a function of the external scale of the turbulent fluctuations and can be estimated to be of the order of α ~0,1-0,2 for Lz ~1 m at the heights close to the water level. At the upper boundary of ABL the parameter a may vary from a ~0,8 for L|| ~15 m to α ~0,005 for L|| ~10 000 m.

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

Февраль, 26, 2025 68 0
Add a comment


Notes
Text copied