# Higher Approximations to Study Statistical Characteristics of Waves in Multiscale Inhomogeneous Media.

1. IntroductionThe presence of random inhomogeneities in a radio signal path can have a profound effect on functioning of radiotechnical systems. Fluctuations (scintillations) in a signal in the randomly inhomogeneous ionosphere, in particular, reduce resolving power of the global Earth satellite monitoring and the rate of data transfer in communication systems. On the other hand, observations of scintillation of radio signals from such extraterrestrial sources as stars and artificial Earth satellites contribute significantly to our knowledge of the inhomogeneous structure ofthe ionosphere and troposphere.

The methods for describing scintillation developed in the 60-80s of the last century used the geometrical optics (GO) and Born and Rytov approximations [1-9], as well as the phase-screen method [2, 3, 7, 10-12]. The GO approximation model can only be adopted to describe the effects of inhomogeneities with transverse scales exceeding the Fresnel radius. The effects of inhomogeneities with scales smaller than the Fresnel radius (Fresnel diffraction effects) cannot be correctly described within the GO model.

It is possible to take account of the Fresnel diffraction effects using the perturbations theory. The Born approximation (single scattering theory) uses the perturbation theory to describe the field itself, and in the Rytov approximation, the perturbation theory is employed to derive the complex phase associated with the field logarithm. However, the Rytov and Born approximations take the diffraction Fresnel effects into account only for weak fluctuations; these approximations are not applicable in multipath emergence and in strong focusing on random caustics. Most of the approaches to simulation of the latter effects are based on integral field representations. When the inhomogeneities are in a restricted area away from the source and the observer and the intensity of inhomogeneities is rather weak, it is possible to use GO [10, 11] or Rytov [12] approximations near the inhomogeneous area. To get the field integral representation in this situation, accounting for the diffraction effects and caustics arising away from inhomogeneities, often the Kirchhoff method [2, 5-7, 13] is applied, which goes to the phase-screen method when the inhomogeneous layer thickness is reduced [2, 3, 7].

This approach is extended by the path integral method [14, 15] and the multiple phase-screen method [11, 16, 17]. However, the path integral method represents a solution in the form of a two-dimensional integral continuum. And the second method involves numerical integration, finite but still large enough (up to 20 for optical radiation in turbulent atmosphere [18]) number of phase screens, which makes it difficult to use these methods for diagnostic purposes.

In the Markov approximation, equations for statistical moments [2] are obtained from the wave equation in the small-angle approximation. Unfortunately, using analytics, we can find solutions only to equations for the mean field and for the spatial correlation function [2]. Equations for other moments can be solved with either numerical [19] or asymptotic methods [20].

Another approach involves the Maslov and interference integral methods [21, 22] and relies on the expanding of the solution of the wave equation in plane waves. This approach requires knowledge of a ray field structure. In a number of problems such as wave propagation in random media and when diagnosing inhomogeneous media, this structure is not always known.

Combinations of the Maslov and interference integral methods were used in solving equations for the two-point random function (averaging of which gives the coherence function) [23] and Markov equations for statistical moments [24]. The same approach enabled solving the equation for the product of the desired field and the known solution (termed as "joker wave" in [25]) in extended 3D space for 2D problem [26] and in the extended 5D space for 3D problem [26-29]. The resulting integral representation in the form of a double weighted Fourier transform (DWFT) describes the wave field, considering for simultaneous effects of inhomogeneities of different scales [26-29].

We will use here the DWFT method [26-30], which in the small-angle approximation describes both strong fluctuations associated with the random caustics that form in a large-scale inhomogeneous medium, and the Fresnel type diffraction effects typical for wave propagation in a medium with small-scale inhomogeneities. Unlike the Maslov and interference integral methods, where the Fourier transform is performed in the receiver or transmitter coordinates, respectively, the DWFT method uses the Fourier transform with respect to both receiver and transmitter coordinates. As a result, in the DWFT integral representation, the wave field is expanded in partial waves that do not have amplitude fluctuations in the first approximation. This not only makes it easier to find the statistical moments of the field, but also shows a greater validity range of the DWFT method in statistical problems compared to the Maslov and interference integral methods.

Let us remark that, in the DWFT method, as in other GO generalizations to generate a wave solution, we use "spanning a wave-field cloth on a ray framework" [21]. In the first DWFT approximation [26-30], the contribution of low-amplitude inhomogeneities to the partial wave phase takes the integral from permittivity perturbation along the unperturbed (direct) ray. For such an effect that inhomogeneities of different scales impose on the wave field we shall call the first-order effect. In addition, in a multiscale medium, large-scale inhomogeneities distort the field structure of a wave incident on a small-scale inhomogeneity. At the same time, keeping in mind the several order difference of scales and amplitudes of various inhomogeneities, the distortion of an incident wave on the scales of small-scale inhomogeneities can be significant. The effects which require trajectory variations to be taken into account when it passes through an inhomogeneity will be called the second-order effects.

Here we use the DWFT method when analyzing statistical characteristics of a wave propagating in a multiscale randomly inhomogeneous medium, with and without regard to the second-order effects. The next, second, section describes the DWFT method. The third section employs the perturbation theory in DWFT and shows transition of the second DWFT approximation to the second GO approximation. The fourth section focuses on the use of the first DWFT approximation to derive expressions for the mean field and the two-frequency mutual coherence function. In the fifth section, the second DWFT approximation considers the second-order effects. Also, the expression is derived for the intensity correlation function, and we compare it asymptotically over long distances, with regard to the second-order effects, with the known results of the path integral method. The conclusion summarizes the main results.

2. Method of Double Weighted Fourier Transform (DWFT)

Let us consider the field of a harmonic (with time dependence of the form exp(-i[omega]t)) point source located at the point [r.sub.0] = {[x.sub.0], [y.sub.0], [z.sub.0]} = {[[rho].sub.0], [z.sub.0]}. In a scalar approximation, this field in an inhomogeneous medium is a solution of equation

[DELTA]E (r, [r.sub.0]) + [k.sup.2][epsilon] (r) E (r, [r.sub.0]) = [A.sub.0][delta] (r - [r.sub.0]), (1)

where permittivity [epsilon](r) = [bar.[epsilon]](r) + [??](r) includes the background (deterministic) [bar.[epsilon]](r) and random [??](r) components, respectively; [A.sub.0] is the amplitude of the incident spherical wave; k = [omega]/c = 2[pi]/[lambda] is the wave number; c and [lambda] are the light speed and the wavelength in free space, respectively. Let [bar.[epsilon]](r) = 1; that is, we will consider the effects of random inhomogeneities in the absence of a background inhomogeneity. Quite often, the size [l.sub.[epsilon]] of inhomogeneities exceeds the wavelength such that [kl.sub.[[epsilon]] [much greater than] 1. In this case, a wave is scattered along its propagation path, for example, along the axis z, and usually [2-4] by substituting

E (r, [r.sub.0]) = U (r, exp [ik (z - [z.sub.0])] (2)

and from elliptical equation (1) we pass on to the parabolic equation for the "complex amplitude" U(r, [r.sub.0]):

2ik [partial derivative]/[partial derivative]z U ([rho], z) + [[partial derivative].sup.2]/[partial derivative][[rho].sup.2] U ([rho], z) + [k.sup.2][??](r) U ([rho], z) = 0 (3)

with the initial condition

U ([rho], [z.sub.0]) = -i[A.sub.0][delta] (r - [r.sub.0])/(2k). (4)

In the absence of inhomogeneity with [??](r) = 0 and [A.sub.0] = 1, the solution of problem (3)-(4) takes the form of

U([rho], z) = [U.sub.0] ([rho], z) = -[A.sub.0] [(4[pi](z - [z.sub.0])).sup.-1] exp {ik0.5 [([rho] - [[rho].sub.0]).sup.2]/(z - [z.sub.0])}. (5)

Let us introduce a five-variable function as the product of the unknown function U([[rho].sub.1]; z) and the known function [U.sub.2] * ([[rho].sub.2], z).

V ([[rho].sub.1], [[rho].sub.2], z) = U([[rho].sub.1], z) [U.sub.2] * ([[rho].sub.2], z). (6)

For function [U.sub.2] * ([[rho].sub.2], z), sometimes called a "joker wave" [25], we will take [Usub.0.sup.*]( [[rho].sub.2], z), which satisfies the equation

2ik[partial derivative][U.sub.2.sup.*]/[partial derivative]z - [[partial derivative].sup.2][U.sub.2.sup.*]/[partial derivative][[rho].sup.2.sub.2] = 0, [U.sub.2.sup.*] ([[rho].sub.a], [z.sub.0]) = i[delta] ([[rho].sub.2] - [[rho].sub.20])/(2k). (7)

Taking into account (3)-(4) and (7), we can write an equation for V([[rho].sub.1], [[rho].sub.2], z):

[mathematical expression not reproducible] (8)

In (8), we introduce new variables:

[[rho].sub.+] = ([[rho].sub.1] + [[rho].sub.2])/2,

[[rho].sub.-] = [[rho].sub.1] - [[rho].sub.2],

[[rho].sub.0+] = ([[rho].sub.01] + [[rho].sub.02]/2,

[[rho].sub.0-] = [[rho].sub.01] - [[rho].sub.02]. (9)

Let us turn from the function V([[rho].sub.1], [[rho].sub.2], z) to its Fourier image with respect to variables [[rho].sub.0+] and [[rho].sub.-]:

[mathematical expression not reproducible] (10)

[mathematical expression not reproducible] (11)

From (8), we get the equation for [??](s, p, z):

[mathematical expression not reproducible] (12)

It is easy to find a solution to (12) in the absence of inhomogeneities, that is, if [??] = 0,

[mathematical expression not reproducible] (13)

In the GO method and its generalizations, solution of the wave equation with large fc is presented as a Debye series in inverse powers of fc. Here, we will also find the solution of problem (12) in the form of

[??] = [[infinity].summation over (0)] [A.sup.(m)]/[(ik).sup.m] exp (ik[phi]). (14)

Substituting (14) into (12) and equating coefficients of equal power in fc, we derive a system of recurrent equations for [phi] and [A.sup.(m)]. If we limit (14) to the first term, we get

[??] [approximately equal to] [A.sup.(0)] exp (ik[phi]) = A exp (ik[phi]), (15)

where the phase path (eikonal) [phi] and the amplitude A satisfy equations

[mathematical expression not reproducible] (16)

[mathematical expression not reproducible] (17)

Initial conditions for eikonal (16) and transfer (17) equations can be obtained from (13):

[mathematical expression not reproducible] (18)

Solving equations in partial first-order derivatives (16)-(17) can be reduced to solving equations of characteristics [1], and thus from (2), (6), (10), we obtain the following expression:

[mathematical expression not reproducible] (19)

where

[mathematical expression not reproducible] (20)

In (20), integration is performed along the characteristics (rays) that satisfy ray equations

[mathematical expression not reproducible] (21)

with initial conditions

[[rho].sub.+] (z) = [[rho].sub.-] ([z.sub.0]) = 0

[P.sub.-] ([z.sub.0]) = s,

[p.sub.+] (z) = p. (22)

Our solution (19) contains "excess" variables [[rho].sub.02] and [[rho].sub.2]. They can be expressed through coordinates [[rho].sub.01] and [[rho].sub.1]. This corresponds to the transition to a certain surface in the five-dimensional space {[[rho].sub.1], [[rho].sub.2], z}. In particular, just like in [26], we can choose this connection from the condition that there are no transverse coordinates [[rho].sub.+], [[rho].sub.0-] in the arguments of the function [??]([rho], z):

[[rho].sub.+] + [[rho].sub.0-]/2 = ([[rho].sub.1] + [[rho].sub.2] + [[rho].sub.01] - [[rho].sub.02])/2 = 0 (23)

In this case, having set [[rho].sub.1] = [rho], [[rho].sub.01] = [[rho].sub.0] and taking (5) into account, (19)-(21) can be written as follows:

[mathematical expression not reproducible] (24)

where

[mathematical expression not reproducible] (25)

[mathematical expression not reproducible] (26)

[mathematical expression not reproducible] (27)

[mathematical expression not reproducible] (28)

is a field in the absence of inhomogeneities.

3. The Perturbation Theory in DWFT Solution

To use solution (24), ray equations (27) with the initial conditions (22) should be solved. Since the exact solution of ray equations for an arbitrary random field [??](r) is unknown, we have to use approximate methods.

Supposing the permittivity fluctuations [??](r) typical of dispersion [[sigma].sup.2.sub.[epsilon]] = <[[??}.sup.2]> are small, that is, considering

[[sigma].sup.2.sub.[epsilon]] [much less than] 1, (29)

let us use the perturbation theory here.

Assume

[[rho].sub.+,-] (z') = [[bar.[rho]].sub.+,-] (z') + [[??].sub.+,-] (z'), (30)

where

[absolute value of ([[??].sub.+,-] (z'))] = O([[sigma].sub.[epsilon]). (31)

[[bar.[rho]].sub.+,-](z') are solutions of ray equations (27) and (22) in the absence of inhomogeneities ([??] = 0), and the solutions have the form

[[bar.[rho]].sub.+] (z') = p (z' - z),

[[bar.[rho]].sub.-] (z') = s (z' - [z.sub.0]),

[[bar.[rho]] (p, s, z') = [[bar.[rho]].sub.+] (z') + [[bar.[rho]].sub.-] (z')

= p (z' - z) + s (z' - [z.sub.0]). (32)

For the trajectory variation [[bar.[??]].sub.+,-] (z'), there are equations

[mathematical expression not reproducible] (33)

with zero initial conditions

[[??].sub.+] (z) = [[??].sub.-] ([z.sub.0]) = [[??].sub.+] (z) = [[??].sub.-] ([z.sub.0]) = 0. (34)

Use iterations to find approximate solutions (33)-(34):

[mathematical expression not reproducible] (35)

Inserting (30), (32), and (35) in (24)-(27) and considering that [gamma] = O([[??].sup.2]), we obtain

[mathematical expression not reproducible] (36)

where

[bar.[phi]] (s, p) = -2ps (z - [z.sub.0]) - 2s[[rho].sub.0] + 2p[rho], (37)

[mathematical expression not reproducible] (38)

[mathematical expression not reproducible] (39)

If trajectory variations are less than the minimum size of inhomogeneities

[??] (z') < [l.sub.[epsilon]] min (40)

in (38), we can expand the function [mathematical expression not reproducible] in powers [??](z') and get

[mathematical expression not reproducible] (41)

Taking (41) into account, expression (36) gives the second DWFT approximation, in the integrand of which the terms O([[absolute value of ([??])].sup.3]) are discarded in phase and O([[absolute value of ([??])].sup.2]) in the amplitude, respectively.

Compare this approximation to the second GO approximation. To do this, we calculate (36), using the stationary phase method, assuming for simplicity that [rho] = [[rho].sub.0] = 0:

[mathematical expression not reproducible] (42)

where

[mathematical expression not reproducible] (43)

[mathematical expression not reproducible] (44)

Exactly the same expression was obtained in the GO approximation in [31]. The first and second terms on the right-hand side of (43) describe the phase path variations caused by a phase velocity change and by a trajectory variation, respectively. Typically, in statistical problems the first approximation is used, and the next approximation is calculated mostly to determine the validity range of the first approximation [4]. However, in some problems, estimation of the second approximation is of great importance. For example, in global satellite navigation systems, the first approximation describes a first-order error that is eliminated in dual-frequency measurements, and the error in these measurements is determined by the next approximation [31].

From derivation of second GO approximation (42), we can see that, for its applicability, the condition of stationary phase method applicability should be met. This condition goes to the condition of the excess of sizes of inhomogeneities over the Fresnel radius

[r.sub.fr] = (z - z')(z' - [z.sub.0])/[k (z - [z.sub.0])] < [l.sub.[epsilon]]. (45)

In addition, condition (40) that sizes of inhomogeneities exceed the trajectory variations should be fulfilled (see also [1]).

4. Statistical Characteristics of the Wave Field in the First DWFT Approximation

In the first approximation, when condition (40) is met, the DWFT solution is

[mathematical expression not reproducible] (46)

where

[mathematical expression not reproducible] (47)

Let us consider [??]([rho], z) a real normal quasi-homogeneous random field with a zero mean and correlation function

[mathematical expression not reproducible] (48)

where [DELTA][rho] = [[rho].sub.1] - [[rho].sub.2], [[rho].sub.[eta]] = ([[rho].sub.1] + [[rho].sub.2])/2, [xi] = [z.sub.1] - [z.sub.2], [eta] = ([z.sub.1] + [z.sub.2])/2.

With respect to the [??]([rho], z) field quasi-homogeneity, the first two variables of the function [[PSI].sub.[epsilon]]([DELTA][rho], [xi], [[rho].sub.[eta]], [eta]) change much faster than the others.

4.1. Mean Field. After averaging (46), we get

[mathematical expression not reproducible] (49)

where

[mathematical expression not reproducible] (50)

[mathematical expression not reproducible] (51)

Here we take into account the fact that the size of an inhomogeneous area exceeds the correlation radius greatly and use the delta-correlation [??]([rho], z) (see [2], for example), which is related to the fact that in the small-angle approximation when integrating (51) we can neglect the [rho]([xi], [eta]) dependence on [xi].

Considering the absence of the [[sigma].sup.2.sub.[phi]] dependence on s, p, from (49), we obtain

<E([rho], z)> = [E.sub.0] ([rho], z) exp {-[[sigma].sup.2.sub.[phi]]/2}. (52)

This result for the mean field was derived earlier in the DWFT approximation for the 2D-inhomogeneous medium in [26]. It can also be obtained in the Markov and first phase (i.e., when amplitude fluctuations are neglected) GO approximations [2].

4.2. Two-Frequency Mutual Coherence Function. The two-frequency mutual coherence function for the two spherical waves is derived as follows:

[[GAMMA].sub.E] ([[rho].sub.1], [[rho].sub.2], [[omega].sub.1], [[omega].sub.2])

= <E([[rho].sub.1], [[omega].sub.1], z, [z.sub.0]) [E.sup.*] ([[rho].sub.2], [[omega].sub.2], [[omega].sub.2], z, [z.sub.0])>. (53)

By substituting (46) in (53) and averaging, we get

[[GAMMA].sub.E] ([[rho].sub.1], [[rho].sub.2], [[omega].sub.1], [[omega].sub.2]) = [E.sub.0] ([[rho].sub.1], [[omega].sub.1]) [E.sub.0.sup.*] ([[rho].sub.2], [[omega].sub.2]) [[??].sub.2] ([[rho].sub.1], [[rho].sub.2], [[omega].sub.1], [[omega].sub.2]), (54)

where

[mathematical expression not reproducible] (55)

[mathematical expression not reproducible] (56)

After averaging (56), we obtain

[mathematical expression not reproducible] (57)

[mathematical expression not reproducible] (58)

Substituting integration variables in (57) [s.sub.1,2] = [bar.s]s [+ or -] s/2, [p.sub.1,2] = [bar.p] [+ or -] p/2 and taking integrals over [bar.s] and [bar.p], we get

[mathematical expression not reproducible] (59)

Here

[mathematical expression not reproducible] (60)

and the notations are introduced (see [19], for example):

[k.sub.0] = ([k.sub.1] + [k.sub.2])/2,

X = ([k.sub.2] - [k.sub.1])/2[k.sub.0],

[[??].sub.0] = [k.sub.0] ([X.sup.-1] - X),

[DELTA][[rho].sub.0] = [[rho].sub.01] - [[rho].sub.02]. (61)

A spatial coherence function can be obtained from two-frequency mutual coherence function (53), (59)-(60) by setting the frequency separation X to zero:

[mathematical expression not reproducible] (62)

Expression (62), like formula (52) for the mean field, coincides with the results of the Markov approximation [2], and it can also be obtained from the phase GO approximation.

For small but finite X, when analyzing the narrowband signals propagation, expression (60) can be calculated asymptotically, and from (53) we obtain

[mathematical expression not reproducible] (63)

where

[mathematical expression not reproducible] (64)

The condition of applicability of asymptotic expression (63) is

[[??].sub.fr] < [l.sub.[epsilon]], (65)

where

[mathematical expression not reproducible] (66)

For narrow-band signals, when X [much less than] 1, condition (66) is much weaker than the condition of GO approximation (45).

5. Second-Order Effects in an Inhomogeneous Medium

Here we account for the fact that, in inhomogeneous media, with distance from the source, the trajectory fluctuations grow; at a distance great enough, these fluctuations become larger than sizes of some inhomogeneities, and condition (40) that was used above is no more applicable, at least for small-scale inhomogeneities. So here, in (38) we leave the trajectory variations in the argument. On the other hand, using (36), we neglect phase addition [[??].sub.1](s, p). Thus, we use an expression similar to (47)

[mathematical expression not reproducible] (67)

but with another phase fluctuation:

[mathematical expression not reproducible] (68)

where [bar.[rho]](z') and [??](z') are described by expressions (32) and (35), respectively.

As we can see from (68), averaging of expression (67), required to calculate statistical characteristics, is complicated by the nonlinear dependence of (68) on [??]([rho], z). To solve this problem, we use a hybrid approach [5, 9]. We assume that the random permittivity field [??]([rho], z) can be represented as the sum of two independent random fields:

[??]([rho], z) = [[??].sub.1] ([rho], z) + [[??].sub.2] ([rho], z), (69)

where [[??].sub.1]([rho], z) and [[??].sub.2]([rho], z) correspond to large-scale and small-scale inhomogeneities, respectively. With that, we account for the trajectory variations only when estimating the contribution of small-scale inhomogeneities. In this case,

[[??].sub.d] (s, p) = [[??].sub.1] (s, p) + [[??].sub.2] (s, p), (70)

[mathematical expression not reproducible] (71)

[mathematical expression not reproducible] (72)

Furthermore, we consider the fact that, due to a rapid drop in the inhomogeneities spectrum, the contribution of small-scale inhomogeneities to trajectory variations in (35) can be neglected:

[mathematical expression not reproducible] (73)

After substituting (70)-(73) in (67), we get

[mathematical expression not reproducible] (74)

5.1. Mean Field. Since [[??].sub.2] is a linear part of exponent argument in (74), being independent of [[??].sub.1], averaging (71) only over realizations of [[??].sub.2] yields

[mathematical expression not reproducible] (75)

where, with the delta-correlation of [[??].sub.2], just as in derivation of (50), we have

[mathematical expression not reproducible] (76)

After averaging (75) over realizations [[??].sub.1], we get

<E ([rho], z)> = [E.sub.0] ([rho], z) exp {[[sigma].sup.2.sub.[phi]]/2}, (77)

where

[mathematical expression not reproducible] (78)

From the comparison between (77) and (52), it is apparent that the consideration of the incident wave distortions does not affect the mean field.

5.2. Two-Frequency Mutual Coherence Function. Substituting (74) in (55) gives

[mathematical expression not reproducible] (79)

By averaging the resulting expressions over realizations [[??].sub.2], we get

[mathematical expression not reproducible] (80)

where

[mathematical expression not reproducible] (81)

We take into account the fact that for small-scale inhomogeneities

[mathematical expression not reproducible] (82)

and rewrite (80):

<[W.sub.2] ([[rho].sub.1], [[rho].sub.2], [[omega].sub.1], [[omega].sub.2]) | [[??].sub.1]> [approximately equal to] [W.sub.2d1] + [W.sub.2d2], (83)

where

[mathematical expression not reproducible] (84)

[mathematical expression not reproducible] (85)

In (84), [W.sub.2] is a function (56) used above to investigate the two-frequency mutual coherence function, without regard to second-order effects for [??] = [[??].sub.1].

When calculating (85), we have to take into account the relationship of the spectrum and the correlation function

[mathematical expression not reproducible] (86)

Besides, considering the smoothness of large-scale inhomogeneities, we assume

[mathematical expression not reproducible] (87)

As a result, by substituting (86)-(87) in (85), we obtain

[mathematical expression not reproducible] (88)

We average (88) over large-scale fluctuations [[??].sub.1] and turn to variables [s.sub.1,2] = [bar.s] [+ or -] s/2, [p.sub.1,2] = [bar.p] [+ or -] p/2 and evaluate integrals over [bar.s] and [bar.p]. After summing up with [W.sub.2d1], we get

[mathematical expression not reproducible] (89)

where

[mathematical expression not reproducible] (90)

If we assume

[kappa] [absolute value of ([z.sub.+] - [eta])]/([[??].sub.0][l.sub.[epsilon]]) [much less than] 1, (91)

that is, neglect the dependence of (90) on [kappa], from (89), we get previously obtained expression (57). This means that when condition (91) is met, the second-order effects, which are related to the front curvature of a wave incident on an inhomogeneity, are not significant in the two-frequency mutual coherence function.

Note that, with the bandwidth decreasing, [[??].sub.0] grows, and for [DELTA][omega] = [[omega].sub.2] - [[omega].sub.1] = 0, that is, for a monochrome signal, conditions (91) are always met. It follows that the second-order effects do not affect the behavior of the spatial correlation function.

5.3. Intensity Spatial Correlation. As the previous paragraph shows, the second-order effects are not significant in the behavior of the field spatial correlation function. Now let us consider the spatial correlation function for intensity.

As the intensity equals I(r) = [[absolute value of (E(r))].sup.2], the correlation function for intensity equals

[mathematical expression not reproducible] (92)

From (62), it follows that <I([rho], z)> = <[[absolute value of (U ([rho], z))].sup.2]) = [[GAMMA].sub.E](0, 0) = [[absolute value of ([E.sub.0])].sup.2]. To determine the function

[[GAMMA].sub.I] ([[rho].sub.1], [[rho].sub.2], z) = <[W.sub.4] ([[rho].sub.1], [[rho].sub.2])> = <I ([[rho].sub.1], z) I ([[rho].sub.2], z)), (93)

we will use DWFT solution (74) with respect to the second-order effects

[mathematical expression not reproducible] (94)

where

[mathematical expression not reproducible] (95)

[mathematical expression not reproducible] (96)

By averaging (95) in small-scale inhomogeneities and taking their small value into account, we get

[mathematical expression not reproducible] (97)

where

[mathematical expression not reproducible] (98)

As above, in the analysis of the two-frequency mutual coherence function, we take into account relationship (86) of the correlation function with the spectrum and the smoothness of large-scale inhomogeneities (87):

[mathematical expression not reproducible] (99)

After substituting (99) in (97) and averaging (97) over [[??].sub.1] in the resulting integral representation for [[??].sub.I]([[rho].sub.1], [[rho].sub.2], z) = ([[??].sub.4]([[rho].sub.1], [[rho].sub.2])>,we turn to new integration variables

[s.sub.1,2] = [q.sub.1] + [v.sub.1]/2 [+ or -] {[u.sub.1] + [[tau].sub.1]/2),

[s.sub.3,4] = [q.sub.1] [+ or -] [v.sub.1]/2 [+ or -] ([u.sub.1] - [[tau].sub.1]/2),

[p.sub.1,2] = [q.sub.2] + [v.sub.2]/2 [+ or -] ([u.sub.2] + [[tau].sub.2]/2),

[p.sub.3,4] = [q.sub.2] - [v.sub.2]/2 [+ or -] ([u.sub.2] - [[tau].sub.2]/2) (100)

and evaluate integrals over [q.sub.1,2] and [u.sub.1,2]. As a result, we obtain

[mathematical expression not reproducible] (101)

where

[mathematical expression not reproducible] (102)

[mathematical expression not reproducible] (103)

[mathematical expression not reproducible] (104)

[mathematical expression not reproducible] (105)

[mathematical expression not reproducible] (106)

[mathematical expression not reproducible] (107)

[bar.[rho]] (p, s, z') = p (z' - z) + s (z' - [z.sub.0]). (108)

At small distances, when condition (91) is met and the incident wave distortion effects can be neglected, (101)-(107) yield

[mathematical expression not reproducible] (109)

From expression (109), we can obtain the known results of the phase-screen method, Rytov and geometrical optics approximations for intensity correlation and scintillation index.

With increasing layer thickness, when (91) is not met, we use the approach of [14] for analyzing (101)-(107). We have to bear in mind that the maximum contribution at great distance belongs to the neighborhoods of the areas [v.sub.1] = [v.sub.2] = 0 and [[tau].sub.1] = [[tau].sub.2] = 0. Resulting from the asymptotic evaluation of integrals in (102)-(106) with [absolute value of (z - [z.sub.0])] [right arrow] [infinity] for normalized correlation function [K.sub.I]([DELTA][rho]) = [[PSI].sub.I]([DELTA][rho])/[<I>.sup.2] = [[GAMMA].sub.I]([DELTA][rho], z)/[<I>.sup.2] - 1, we get

[K.sub.I] ([DELTA][rho]) [approximately equal to] [B.sup.(1).sub.I] ([DELTA][rho]) + [B.sup.(2).sub.I] ([DELTA][rho]) + [B.sup.(2).sub.I] ([DELTA][rho]), (110)

where

[mathematical expression not reproducible] (111)

[mathematical expression not reproducible] (112)

[mathematical expression not reproducible] (113)

[mathematical expression not reproducible] (114)

Hence, for the scintillation index [[beta].sup.2] of a wave propagating in a multiscale randomly inhomogeneous medium, we get

[mathematical expression not reproducible] (115)

With respect to the filtering effect of the coefficient {1 - cos[[[kappa].sup.2] ([z.sub.+] - z)([z.sub.+] - [z.sub.0])/(k(z - [z.sub.0]))]} in (112), (113), and (115), we can substitute the spectrum of small-scale inhomogeneities [mathematical expression not reproducible] with the full spectrum [mathematical expression not reproducible]. Since the contribution of small-scale inhomogeneities is small, we can substitute [[bar.[PSI]].sub.[epsilon]1] (0, [eta]) [approximately equal to] [[bar.[PSI]].sub.[epsilon]] (0, [eta]) = [[bar.[PSI]].sub.[epsilon]1](0, [eta]) + [[bar.[PSI]].sub.[epsilon]2] (0, [eta]) in (112), (113), and (115). As a result, we get expressions for the correlation function [K.sub.}([DELTA][rho]) and scintillation index [[beta].sup.2], for a plane wave such that [z.sub.0] [right arrow] -[infinity] go to results [14] of asymptotic computation of the strong expressions obtained by the path integral method.

6. Conclusion

The use of the DWFT method and hybrid approach together enabled us to derive expressions for statistical characteristics of the wave field in a multiscale randomly inhomogeneous medium. Upon that we took into account the second-order effects associated with the ray trajectory variations within sizes of small-scale inhomogeneities. It is shown that, similar to the Markov approximation, the mean field and the spatial coherence function coincide with their phase geometrical optics approximations. Within the DWFT, we were able to find a solution for the frequency coherence function and to show that, for narrow-band signals, the second-order effects do not significantly affect the frequency correlation. Moreover, it is shown that, with decreasing bandwidth, phase geometrical optics approximation can be used to calculate frequency correlation.

A more significant role in the intensity correlation belongs to the second-order effects. Analysis of the obtained expressions showed that while the second-order effects can be neglected for the small thickness of the inhomogeneous layer, these effects become noticeable as the layer thickness increases. It is interesting to note that the resulting intensity correlation function at large distances goes to the results obtained earlier by the path integral method. This shows that the intensity fluctuations of a wave propagating in a multiscale randomly inhomogeneous medium at great distance are significantly affected by the second-order effects, that is, the incident wave distortion effects.

https://doi.org/10.1155/2018/1570407

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The author is very grateful to S. V. Filippova for her assistance in preparing the English version of the manuscript. This work was supported by a project part of the Government Assignment for Scientific Research from the Ministry of Education and Science, Russia (no. 3.903.2017/4.6).

References

[1] Y. A. Kravtsov and Y. I. Orlov, Geometrical Optics of Inhomogeneous Media, Springer, Berlin, Germany, 1990.

[2] S. M. Rytov, Y. A. Kravtsov, and V. I. Tatarskii, Introduction to Statistical Radiophysics. Vol. 4, Wave Propagation Through Random Media, Springer, New York, NY, USA, 1989.

[3] A. Ishimaru, Wave Propagation and Scattering in Random Media. Vol 2. Multiple Scattering, Turbulence, Rough Surfaces and Remote Sensing, Academic Press, New York, NY, USA, 1978.

[4] V. I. Tatarskii, "The Effect of a Turbulent Atmosphere on Wave Propagation," Tech. Rep., National Technical Information Service, Springfield, Va, USA, 1971.

[5] Y. A. Kravtsov, S. I. Feyzulin, and A. G. Vinogradov, Radio Wave Propagation through The Earth's Atmosphere, M.: Radio i Svyaz, 1983.

[6] V. E. Kunitsyn and E. D. Tereshchenko, Ionospheric Tomography, Springer, Berlin, Germany, 2003.

[7] Y. N. Barabanenkov, Y. A. Kravtsov, S. M. Rytov, and V. I. Tamarskii, "Status of the theory of propagation of waves in a randomly inhomogeneous medium," Soviet Physics Uspekhi, vol. 13, no. 5, pp. 551-575, 1971.

[8] B. N. Gershman, L. M. Erukhimov, and Y. Y. Yashin, Wave Phenomena in The Ionosphere and Space Plasma, M.: FIZMATLIT, 1984.

[9] Y. A. Kravtsov, "Propagation of electromagnetic waves through a turbulent atmosphere," Reports on Progress in Physics, vol. 55, no. 1, article 002, pp. 39-112, 1992.

[10] E. N. Bramley, "The diffraction of waves by an irregular refracting medium," Proceedings of the Royal Society A Mathematical, Physical and Engineering Sciences, vol. 225, pp. 515-518, 1954.

[11] L. M. Erukhimov and V. P. Uryadov, "Frequency correlation of radio-wave fluctuations behind n random screens," Radiophysics and Quantum Electronics, vol. 11, no. 12, pp. 1050-1058, 1972.

[12] V. E. Gherm, N. N. Zernov, and H. J. Strangeways, "Propagation model for transionospheric fluctuating paths of propagation: simulator of the transionospheric channel," Radio Science, vol. 40, no. 1, pp. RS1003-RS10003, 2005.

[13] R. L. Fante, "Two-position, two-frequency mutual-coherence function in turbulence," Journal of the Optical Society of America, vol. 71, no. 12, pp. 1446-1451, 1981.

[14] V. I. Klyatskin, Stochastic Equations and Waves in Random Inhomogeneous Media, 1980, [Translation: Ondes et Equations Stochastique dans les Milieux Aleatoirement Non-homogenes, Editions de Physiques, Besancon, 1985 (In French)].

[15] M. I. Charnotskii, J. Gozani, V. I. Tatarskii, and V. U. Zavorotny, "IV wave propagation theories in random media based on the path-integral approach," Progress in Optics, vol. 32, pp. 203-266, 1993.

[16] V. P. Kandidov, "Monte Carlo method in nonlinear statistical optics," Physics-Uspekhi, vol. 39, no. 12, pp. 1243-1272, 1996.

[17] D. L. Knepp, "Multiple phase-screen calculation of the temporal behavior of stochastic waves," Proceedings of the IEEE, vol. 71, no. 6, pp. 722-737, 1983.

[18] J. M. Martin and S. M. Flatte, "Intensity images and statistics from numerical simulation of wave propagation in 3-D random media," Applied Optics, vol. 27, no. 11, pp. 2111-2126, 1988.

[19] C. H. Liu and K. C. Yeh, "Frequency and spatial correlation functions in a fading communication channel through the ionosphere," Radio Science, vol. 10, no. 12, pp. 1055-1061, 1975.

[20] A. A. Bitjukov, V E. Gherm, and N. N. Zernov, "On the solution of Markov's parabolic equation for the second-order spaced frequency and position coherence function," Radio Science, vol. 37, no. 4, pp. 19-1-19-9, 2002.

[21] Y. A. Kravtsov and Y. I. Orlov, Caustics, Catastrophes and Wave Fields, Springer, Berlin, Germany, 1999.

[22] M. V. Tinin, "Quasi-geometrical optics methods in the theory of wave propagation through a randomly inhomogeneous layer," Waves in Random and Complex Media, vol. 8, no. 3, pp. 329-338, 1998.

[23] R. Mazar, "High-frequency propagators for diffraction and backscattering in random media," Journal ofthe Optical Society of America A: Optics and Image Science, and Vision, vol. 7, no. 1, pp. 34-46, 1990.

[24] S. Frankenthal, "The mutual coherence function in a scattering channel--a two scale solution," The Journal of the Acoustical Society of America, vol. 85, no. 1, pp. 104-113, 1989.

[25] G. Samelsohn, "High-frequency directed wave propagators: a path integral derivation," IEEE Transactions on Antennas and Propagation, vol. 61, no. 11, pp. 5637-5648, 2013.

[26] Y. A. Kravtsov and M. V. Tinin, "Representation of a wave field in a randomly inhomogeneous medium in the form of the double-weighted Fourier transform," Radio Science, vol. 35, no. 6, pp. 1315-1322, 2000.

[27] M. V. Tinin, S. N. Kolesnik, and Y. A. Kravtsov, "An integral representation for the field of the wave propagating in a medium with random inhomogeneities of different scales," Proceedings of SPIE--The International Society for Optical Engineering, vol. 4678, pp. 605-609, 2002.

[28] Y. A. Kravtsov and M. V. Tinin, "Ray based diffraction tomography of the ionosphere and laboratory inhomogeneous plasma," Cosmic Research, vol. 41, no. 4, pp. 357-358, 2003.

[29] Y. A. Kravtsov, M. V. Tinin, and S. I. Knizhnin, "Diffraction tomography of inhomogeneous medium in the presence of strong phase variations," Journal of Communications Technology and Electronics, vol. 56, no. 7, pp. 831-837, 2011.

[30] M. V. Tinin and S. I. Knizhin, "Eliminating the effects of multipath signal propagation in a smoothly inhomogeneous medium," Radiophysics and Quantum Electronics, vol. 56, no. 7, pp. 413-421, 2013.

[31] B. C. Kim and M. V. Tinin, "Contribution of ionospheric irregularities to the error of dual-frequency GNSS positioning," Journal of Geodesy, vol. 81, no. 3, pp. 189-199, 2007.

M. V. Tinin (iD)

Irkutsk State University, 20 Gagarin Blvd, Irkutsk 664003, Russia

Correspondence should be addressed to M. V Tinin; mtinin@api.isu.ru

Received 29 April 2017; Accepted 4 January 2018; Published 7 February 2018

Academic Editor: Xavier Leoncini

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Tinin, M.V. |

Publication: | Advances in Mathematical Physics |

Article Type: | Report |

Geographic Code: | 1USA |

Date: | Jan 1, 2018 |

Words: | 6026 |

Previous Article: | A Note on Finsler Version of Calabi-Yau Theorem. |

Next Article: | From the Kinematics of Precession Motion to Generalized Rabi Cycles. |

Topics: |