Transformation of the Spatial Spectrum of Scattered Radio Waves in the Conductive Equatorial Ionosphere

: The oblique radio wave incidence on a turbulent equatorial conductive collision plasma layer is considered. The “Compensation Effect” has been discovered by us. A complex refractive index of the equatorial terrestrial ionosphere has been derived for the ﬁrst time. Second-order statistical moments of the spatial power spectrum (SPS) of scattered radio waves are obtained for the ﬁrst time using the WKB method, taking into account the asymmetry of the problem: the inclined incidence of the wave on a plasma boundary and the asymmetry of the magneto-ionic parameters. It was established for the ﬁrst time that a certain direction exists along which the inclined incidence radio wave on a plasma layer and the anisotropy parameters of a magnetoplasma compensate each other. This result will have great practical application in communication. In this case, the SPS of scattered radio waves neither widens nor is its maximum displaced. The behavior of this spectrum versus distance propagated by radio waves in the conductive equatorial ionosphere is analyzed numerically for different penetration angles and anisotropy factors of asymmetric anisotropy electron density irregularities. It was shown that the anisotropy factor of elongated plasmonic structures has a substantial inﬂuence on the “Compensation Effect” of scattered ordinary and extraordinary waves penetrating in the conductive collision ionospheric plasma the slab. Numerical calculations are carried out for the anisotropic Gaussian correlation function applying IRI experimental data.


Introduction
Randomly varying magneto-ionic plasma parameters containing electron concentration inhomogeneities have a significant impact on communication accuracy and navigation in a wide frequency band.Randomly varying plasma parameters have a substantial influence on the radio waves.Conductivity, collision between the plasma particles, complex refractive index, and its imaginary part cause significant impact on the statistical characteristics of scattered radio waves propagating in the conductive collision ionospheric magnetized (COCOIMA) plasma, particularly in the equatorial region of the terrestrial atmosphere.Asymmetry of the task can also lead to the increase in turbulence in magnetoplasma.The direction of both the geomagnetic field and an inclined incident of radio wave on a conductive turbulent collision magnetoplasma is a reason for the asymmetry of the problem.Inclined incident plane radio waves on the boundary of absorptive plasma layer penetrating into the turbulent plasma become inhomogeneous.
In many papers [1][2][3][4], second-order statistical characteristics of the multiple scattered radio waves in the terrestrial ionosphere have been investigated analytically and numerically using experimentally observed data.It was shown [5] that at small angles of incidence wave on a surface between vacuum and isotropic collisional turbulent plasma, the spatial (angular) power spectrum (SPS) of a multiple scattered radio wave monotonically broadens with distance approaching a certain asymptotic value [5].New phenomena arise increasing at increasing asymmetry of the problem: the incline incidence of radio waves on the anisotropic collision plasma slab.Different spectral components of the SPS attenuate variously; the width of the SPS anomalously broadens and becomes asymmetric with respect to its maximum.The peak of the SPS maximum is non-monotonically shifted to the normal relative to the interface between vacuum and collisional magnetoplasma.This is a consequence of the asymmetry of the task.Similar phenomena arise not only at the inclined incident wave on a surface; such anisotropy can also be an internal property of a turbulent plasma itself.
Turbulent conductive collisional magnetoplasma is a randomly inhomogeneous absorptive medium.The inclined incidence of a plane radio wave on a semi-infinite collision magnetoplasma layer with electron concentration fluctuations was considered [6][7][8] in a small-angle scattering approximation.Spectrum width varies nonmonotonically with distance from a plasma boundary.The widening of the SPS of scattered radio waves and shift of its maximum for both anisotropic Gaussian and power-law correlation spectral functions of electron concentration fluctuations in the turbulent collisional magnetoplasma were analyzed [7][8][9], applying the Wentzel-Kramers-Brilluen method.In this case, the width of the spectrum of a scattered electromagnetic wave substantially exceeds the width of the collisionless plasma.The "Compensation Effect" in the polar terrestrial ionosphere was considered in [10].
Physical processes in the polar and equatorial ionospheres are very different.Statistical moments of scattered radio waves propagating in the equatorial terrestrial ionosphere are investigated for the first time.They include the following anisotropy parameters: anisotropy of the ionosphere, taking into account Hall's, Pedersen, and longitudinal conductivity fluctuations, the direction of an external magnetic field, anisotropy of elongated plasmonic structures due to the diffusion processes in the filed aligned and field perpendicular directions containing anisotropy factor and the inclination angle of anisotropic ionospheric irregularities with respect to the geomagnetic lines of forces, incline incidence of wave on a turbulent plasma layer.A Complex refractive index of the equatorial terrestrial ionosphere has been obtained for the first time.Evaluation of the SPS of scattered electromagnetic waves propagating in the COCOIMA plasma is investigated in this paper using the complex ray-(optics) approximation.Compensation conditions are established for the experimentally observed power-law spectral function of electron density irregularities applying the statistical ionospheric model IRI.

Radio Wave Propagation in a Homogeneous Conductive Magnetized Plasma
The electric field satisfies the wave equation: where: k 0 = ω/c is the wavenumber of an incident wave having frequency ω; ∆ is the Laplacian, ) are the second rank permittivity and conductivity tensors of the COCOIMA plasma, respectively.Let a plane radio wave be incident from vacuum on a semi-infinite homogeneous COCOIMA magnetoplasma layer; Y axis coincides with an external homogeneous magnetic field H 0 .Conductive homogeneous collision magnetoplasma have permittivity tensor components: where Here, ϕ is the polar angle between of the wave vector k on the XOZ plane.
If the directions of radio wave propagation and the homogeneous imposed magnetic field are perpendicular (θ = 90 0 ) and ϕ = 0 0 , we obtain two roots only the Hall's and Pedersen conductivities give the contribution in Equation ( 3), E y = 0 .Substitution (4) into (3) gives (E z /E x ) = ∓ i.Consequently, we have a right-hand circularly polarized (RHCP) wave (upper sign) and a left-hand circularly polarized (LHCP) wave (lower sign) propagating along the homogeneous magnetic field.Analyses show that k I I < k I , the rotation is clockwise and, hence, the Faraday rotation angle is For an incident 3 MHz frequency wave, θ F = 0.01.At quasi-longitudinal propagation of wave (θ = 0 0 ), we have the following roots: Contrary to the previous case, the second root contains the longitudinal conductivity.For nonconductive The second root is the same.Ionospheric plasma is an important subject of research in the field of radio physics for radio wave propagation.If θ is an angle between H 0 ||Y and k 0 ⊂ YOZ vectors, s << ε i j , σ i j , complex refractive index of the conductive collision equatorial ionospheric plasma was obtained in [11].
After some transformations from Equation (1), we obtain where: 5) has two zero and two poles, wherein each feature corresponds to one of four major waves (modes) that may exist in the plasma.These modes in the conductive equatorial plasma can be classified as follows: (1) Propagation of wave perpendicular to the external magnetic field (2) Wave propagation along the external magnetic field (θ = 0) these are weakly dumping radio waves propagating in the equatorial terrestrial ionosphere.
In the absence of an external magnetic field and in a non-conductive plasma nonconductivities (σ ij = 0), all waves are reduced to the ordinary wave.At intermediate angles between 0 and π/2, propagating waves will represent some combinations of these major waves.Group velocities v gr of both ordinary and extraordinary (O− and E−) waves generally do not coincide with the phase velocity v ph = ω k/k 2 .The angle Θ between v gr and k in a homogeneous medium, as is known, is determined by the relation [1][2][3] At small angles θ (quasi-longitudinal propagation) N 1,2 has a simple form Upper sign corresponds to the ordinary wave (O-wave), lower sign-to the extraordinary wave (E-wave).It is easy to make sure that the ratio tg Θ 1 / tg Θ 2 determined by the last ratio less than zero, i.e., the group velocities of the O−wave and E−wave deviate from the phase in opposite directions.
It can be concluded that at vertical sounding of the ionosphere, the areas of reflection of O−wave and E−wave are spaced apart in the horizontal direction.Calculations show that at vertical sounding of the ionosphere in both the southern and northern hemispheres, the trajectory of O−wave deviate from the vertical to the equator, and the E−wave in the opposite direction.
If radio waves propagate along the external magnetic field in the equatorial ionosphere, from Equations ( 8) and ( 9) we obtain the RHCP waves are reflected in the ionosphere, while LHCP waves penetrating the ionospheric layers will propagate in the upper ionosphere; Ω e = e H 0 /m e c is the electron gyrofrequency.

Statistical Moments of the Phase Fluctuations
The main theoretical tool describing short wavelength electromagnetic waves propagation in the terrestrial atmosphere is the ray theory.It is one of the most important methods among other asymptotic methods for a number of reasons: one of them is its simplicity and the possibility of obtaining an analytical solution for a wide range of problems that cannot be investigated by accurate or other asymptotic methods.
Let us introduce the Cartesian system of coordinates; a plane layer of a turbulent COCOIMA plasma is radiated by a plane radio wave, refracting at an angle θ 1 relative to the normal to this layer (it is related to the angle of incidence θ i by Snell's law).This is the asymmetry of the problem.The thickness of a plasma layer is L. At the equatorial latitudes, the geomagnetic field is nearly horizontal; therefore, we will assume that the homogeneous magnetic field vector H 0 is directed along the y-axis and the wave vector of an incident wave k 0 ; both are located in the YOZ plane (the main plane).Rectangular components of a complex wave vector on the Y and Z coordinate axes in vacuum are purely real values: Refracting at the interface vacuum-plasma, this wave becomes inhomogeneous, as it obeys the boundary conditions, the Z-projection of its wave vector remains real, but the We assume that the characteristic spatial scale of electron concentration irregularities exceeds the wavelength of an incident wave l >> λ and L >> l.This allows the use the ray-(optics) approximation for the investigation of statistical moments of a scattered radio wave field.
Electron concentration fluctuations in a plasma layer cause fluctuation of a scattered field in the observation point.At small-angle scattering in the ray-(optics) approximation, second-order statistical characteristics are determined by fluctuation of a complex phase ϕ(r) [3] satisfying the eikonal equation k 2 = k 2 0 N 2 , k(r) = − ∇ ϕ complex wave vector is a function of position.Complex refractive index of the COCOIMA plasma contains components of the wave vector N 2 (r) = N 2 (n(r), k x , k z , ω).The eikonal equation can be rewritten as [8] (k where, k ⊥ = k ⊥ (k x , k z ).We can use the following series expansions for the wave vector and the following phase: • leading to the relationship for the phase fluctuation considering only first order small terms of electron density fluctuations n 1 /n 0 .After some algebraic transformations, we obtain integrating Equation ( 15) along the complex characteristics, taking into account that Here, θ is an angle between H 0 and the wave vector of a refracted wave.
Irregularities in the equatorial ionosphere are highly field aligned, and therefore the spread of F turbulence is 2D.At some scale length the spectrum is anisotropic, i.e., the ionosphere is much more narrowly bounded vertically than horizontally.
Taking into account the boundary condition, ϕ 1 (z = 0) = 0 and expand the function ϕ 1 in a two-dimensional Fourier integral; for the phase fluctuations, we obtain where, ; for brevity, the index '0 will be omitted everywhere.Phase fluctuations allow for the investigation of second order statistical characteristics of scattered radio waves.A transverse correlation function of the phase fluctuation may be easily obtained where, ρ z and ρ x are distances between observation points spaced apart at small distances in the main and perpendicular planes, V n (k) is the arbitrary correlation function of electron concentration fluctuations.
At strong fluctuations of the phase < ϕ 1 ϕ * 1 > >> 1, it can be assumed, as usual [3], that they are normally distributed.If the regular phase difference is negligible, the transverse correlation function of the complex field is described by the formula all derivatives of the correlation function of the phase are taken at ρ x = ρ z = 0. Using the 2D Fourier transformation, we obtain the APS having great practical importance.This statistical characteristic is the same as the ray intensity (brightness) in the radiation transport equation [3] Second-order statistical characteristics ∆ z ≡ ∆ k z describe the displacement of maximum of the SPS due to the random variation of electron density, x > determines the broadening of the SPS in the YOZ and XOY planes, respectively.
Substituting (18) into Equation ( 22), we yield the double integrals depend only on the spatial spectrum of electron concentration fluctuations but not on the strength of the fluctuations.As it was mentioned above, the asymmetry of the problem includes all anisotropic parameters of the turbulent conductive magneto-plasma and oblique incidence of the wave.Particularly (a) at inclined incidence wave on the isotropic absorptive plasma ∂k y2 /∂k 0z = −k 0z /k y2 , Υ 2 = 0; (b) at oblique incidence of the wave on a magnetoplasma slab, we can apply Formula (19).From these formulas, follow that anomalous broadening of the APS and shift of its maximum are clearly manifested at Υ 2 k z L > 1. Taking into account that ∂ N 2 /∂ k z ∼ ∂ N 2 /∂ θ, we can directly apply formula (5) differentiating by angle θ.From the Equations ( 18) and ( 24), follow that at Υ 2 = 0 exists the compensation direction, along which the SPS neither broadens nor does its maximum shift.At different incident angles on a plasma layer, the SPS will anomalously broaden, but the mean value of the SPS will displace to the compensation direction, increasing distance travelling by radio wave in the equatorial region of the COCOIMA plasma.

Numerical Calculations
Artificial ionospheric inhomogeneities generating across the geomagnetic field have scales ranging from a meter up to more kilometers.The statistical ionospheric model IRI recommended by URSI as a basic model in the study and prediction of ionospheric propagation of radio waves [12,13] is the most proven technology for describing the spatial distribution of electron density in the equatorial ionosphere.
Observing the characteristic linear scale of electron density irregularities in F region of the equatorial ionosphere is in the range from 100 m to 10 km.Satellite experiments consistently observe a power law spectrum of plasmonic structures, with the spectral index in the interval −1 to −3.
Experimental observation of the equatorial ionospheric irregularities distribution was analyzed by GPS signals registering by ground-based radar stations and by Ionospheric Scintillation Monitors (ISM).Irregularities in the equatorial ionosphere have a bubbles form, created near the magnetic equator on the bottom side of the F2 layer, rising along the geomagnetic lines of forces.These bubbles are generated due to a Rayleigh-Taylor-like instability.Experimentally observing the phase spectral index p describes the strength of electron density irregularities.According to the ISM observations, index p is in the interval 2.5 ≤ p ≤ 3.5.Experimental observations of the anisotropic electron density irregularities in the F-region of the ionosphere show that they have power-law spectrum with a power index p [14,15].The 3D spectral function can be written as follows: where, C 2 n is the mean-square deviation of electron density, K ν (x) is McDonald function, r 0 is the inner scale of turbulence, L 0 = 2π/k 0 is the outer scale.In the interval k 0 r 0 << k r 0 << 1 correlation function is described by the formula [14,15]: where Γ(x) is the gamma function.
In analytical calculations, we apply the power-law spectral function: where, is the anisotropy factorthe ratio of longitudinal and transverse characteristic linear sizes of plasma irregularities.The shape of electron density irregularities has a spheroidal form due to diffusion processes in the field align and field perpendicular directions.Small-scale plasmonic structures of the ionospheric F region are mainly field-aligned [16][17][18][19].
In numerical calculations, we used experimental data on ionospheric parameters for an altitude of 300 km [12,13].
Figure 1 illustrates the deviation of ray trajectories of both O-and E-waves (∆h in meters) relative to the direction connecting the source and the receiver versus the distance between observation points (in km-s) (see formula (12).This figure illustrates the phenomenon of the "fountain effect" of the O-and E-waves, which was revealed in the approximation of geometric optics.In the area of the geomagnetic equator, the geomagnetic field is almost parallel to the Earth's surface.A global empirical ionospheric model IRI (International Reference Ionosphere) was used in numerical calculations of the equatorial latitudes.The behavior of O-and E-waves in the equatorial region is different.In this region of the ionosphere, the plasma rises and gradually turns north in the northern hemisphere and south in the southern hemisphere, which is caused by an increase in the inclination of geomagnetic lines of forces on both sides of the geomagnetic equator.As a result, maxima (or crest) of electron concentration are formed on both sides of the geomagnetic equator, i.e., north and south crests of equatorial anomaly or crests corresponding to ordinary and unusual waves.This effect significantly affects the operation of radio communication, radio navigation, location, etc.The phase spectral index p depends on the state of the ionosphere.In a calm ionosphere, parameter p varies within interval 1.7 ≤ p ≤ 2, while in a highly perturbed ionosphere, the p value can increase to 4 and even to 6.At spectral index p = 2, the variance of the phase fluctuations of the received navigation radio signal is significantly less than at spectral index p ≥ 3, and their difference increases in proportion to the RMS deviation of electronic concentration fluctuations in small-scale irregularities.
Figure 3 shows the 3D correlation function of the phase fluctuations as a function of distances between the observation points orthogonal (η x ) and parallel (η y ) to the magnetic field lines.

Conclusions
The statistical characteristics of the O-and E-radio waves are calculated, applying complex geometrical optics approximation.Anisotropy of the task contains inclined incidence of a radio wave on a plasma boundary, direction of the geomagnetic field, anisotropy of ionospheric conductivities: Hall's, Pedersen, and longitudinal conductivities, and anisotropy of ionospheric plasmonic structures.
A complex refractive index of the COCOIMA plasma of the equatorial region has been obtained for the first time.At radio waves propagation in a homogeneous conductive equatorial ionosphere, it is shown that there exist two weakly dumping RHCP and LHCP spiral-type waves containing Pedersen and Hall's conductivities.Calculating the Faraday angle, it was established that the rotation of the polarization plane is clockwise.The "fountain effect" of the O-and E-waves has been revealed.Crests corresponding to the Oand E-waves on both sides of the geomagnetic equator have been revealed.
Numerical analyzes were carried out for the experimentally observing anisotropic power-law spectral function of electron density fluctuations with the spectral index p = 3. Knowledge of the parameter p allow to determine anisotropy factor of elongated ionospheric irregularities measuring the backscattering cross-sections of radio waves with vertical and inclined soundings of the same ionospheric region (σ 1 /σ 2 ) ∼ l 1−p || .Asymmetry of the task: both the anisotropy factor of stretched electron density irregularities and ionospheric conductivities in the equatorial ionosphere play an important role in the "Compensation Effect"; particularly, they have significant influence on the Ewave.Radar sounding of the same area in the equatorial COCOIMA plasma at different frequencies and incidence angles makes it possible to define ionospheric parameters and the characteristic linear scales of elongated plasmonic structures.It was shown that at fixed anisotropy factor of elongated plasmonic structures and different penetration angles of Oand E-radio waves (asymmetry of the problem) in the COCOIMA plasma, the "Compensation Effect" takes place at small distances from a plasma boundary.At a fixed penetration angle in the plasma slab, the varying anisotropy factor "Compensation Effect" appears at big distances from a plasma boundary in both main and perpendicular planes.
Investigating the features of electromagnetic waves propagation in the equatorial ionosphere allows to solve the problems of radio communication, radio navigation, as well as the problems of diagnosing of anisotropic ionospheric plasmonic structures.

Figure 1 .
Figure 1.Relation of O-and E-waves ray paths deviation to observation points.

Figure 2
Figure2demonstrates the dependence of root mean square (RMS) deviation of the phase random variation versus the RMS deviation of electron concentration fluctuations.The phase spectral index p depends on the state of the ionosphere.In a calm ionosphere, parameter p varies within interval 1.7 ≤ p ≤ 2, while in a highly perturbed ionosphere, the p value can increase to 4 and even to 6.At spectral index p = 2, the variance of the phase fluctuations of the received navigation radio signal is significantly less than at spectral index p ≥ 3, and their difference increases in proportion to the RMS deviation of electronic concentration fluctuations in small-scale irregularities.Figure3shows the 3D correlation function of the phase fluctuations as a function of distances between the observation points orthogonal (η x ) and parallel (η y ) to the magnetic field lines.Figure 4 illustrates the displacement of the maximum of the SPS of the O-(blue curves) and E-waves (red curves) in the conductive collision magnetized plasma at different penetration angle θ and anisotropy factor χ = 2. SPS for the O-wave reaches maximums at L/l || = 49, θ = 5 0 ; and L/l || = 71, θ = 15 0 .In this case, O-waves tend to the compensation direction at L/l || ≈ 260.SPS for E-wave has maxima at L/l || = 60, θ = 10 0 ; and L/l || = 51, θ = 23 0 .The curves corresponding to E-wave reach the compensation direction at L/l || ≈ 210.The "Compensation Effect" of the displacement of the SPS for both waves takes place at a distance L/l || ≈ 800.

Figure 4
Figure2demonstrates the dependence of root mean square (RMS) deviation of the phase random variation versus the RMS deviation of electron concentration fluctuations.The phase spectral index p depends on the state of the ionosphere.In a calm ionosphere, parameter p varies within interval 1.7 ≤ p ≤ 2, while in a highly perturbed ionosphere, the p value can increase to 4 and even to 6.At spectral index p = 2, the variance of the phase fluctuations of the received navigation radio signal is significantly less than at spectral index p ≥ 3, and their difference increases in proportion to the RMS deviation of electronic concentration fluctuations in small-scale irregularities.Figure3shows the 3D correlation function of the phase fluctuations as a function of distances between the observation points orthogonal (η x ) and parallel (η y ) to the magnetic field lines.Figure 4 illustrates the displacement of the maximum of the SPS of the O-(blue curves) and E-waves (red curves) in the conductive collision magnetized plasma at different penetration angle θ and anisotropy factor χ = 2. SPS for the O-wave reaches maximums at L/l || = 49, θ = 5 0 ; and L/l || = 71, θ = 15 0 .In this case, O-waves tend to the compensation direction at L/l || ≈ 260.SPS for E-wave has maxima at L/l || = 60, θ = 10 0 ; and L/l || = 51, θ = 23 0 .The curves corresponding to E-wave reach the compensation direction at L/l || ≈ 210.The "Compensation Effect" of the displacement of the SPS for both waves takes place at a distance L/l || ≈ 800.

Figure 2 .
Figure 2. RMS deviation of σ ϕs versus RMS deviation of σ n for different phase spectral index p.

Figure 3 .
Figure 3. Three-dimensional correlation function of the phase fluctuations.

Figure 4 .
Figure 4. Shift of maximum of the SPS at non-dimensional space parameter.

Figure 5
Figure 5 shows the "Compensation Effect" for the broadening of both waves in the conductive equatorial terrestrial ionosphere at θ = 5 • .The broadening of the SPS for Owave in the main plane has a maximum at χ = 6, propagating distance L/l || ≈ 3.4 • 10 3 ; the broadening tends to the compensation direction at a distance L/l || ≈ 8 • 10 3 .For E-wave, the broadening of the SPS at the anisotropy factor χ = 6 has a maximum at L/l || ≈ 2.2 • 10 3 and tends to the compensation direction at L/l || ≈ 10 4 .The "Compensation Effect" of the broadening of the SPS for E-waves is revealed earlier than for the O-waves.Figure 6 depicts the "Compensation Effect" for the broadening of both waves in the conductive equatorial region of the terrestrial ionosphere at θ = 20 • .Broadening of the SPS corresponding to the O-wave reaches to the compensation direction at small anisotropy factor of electron density irregularities.The "Compensation Effect" for Σ x is revealed at distance L/l || ≈ 5 • 10 3 for O-wave.A varying anisotropy factor in the interval χ = 22 ÷ 25, the "Compensation Effect" is observed at distances L/l || ≈ 2 • 10 4 .
Figure 5 shows the "Compensation Effect" for the broadening of both waves in the conductive equatorial terrestrial ionosphere at θ = 5 • .The broadening of the SPS for Owave in the main plane has a maximum at χ = 6, propagating distance L/l || ≈ 3.4 • 10 3 ; the broadening tends to the compensation direction at a distance L/l || ≈ 8 • 10 3 .For E-wave, the broadening of the SPS at the anisotropy factor χ = 6 has a maximum at L/l || ≈ 2.2 • 10 3 and tends to the compensation direction at L/l || ≈ 10 4 .The "Compensation Effect" of the broadening of the SPS for E-waves is revealed earlier than for the O-waves.Figure 6 depicts the "Compensation Effect" for the broadening of both waves in the conductive equatorial region of the terrestrial ionosphere at θ = 20 • .Broadening of the SPS corresponding to the O-wave reaches to the compensation direction at small anisotropy factor of electron density irregularities.The "Compensation Effect" for Σ x is revealed at distance L/l || ≈ 5 • 10 3 for O-wave.A varying anisotropy factor in the interval χ = 22 ÷ 25, the "Compensation Effect" is observed at distances L/l || ≈ 2 • 10 4 .

Figure 5 .
Figure 5.The broadening of the SPS versus non-dimensional spatial parameter in the main plane.

Figure 6 .
Figure 6.The broadening of the spectrum perpendicular to the main plane.
2, ω p (r) = 4 π n e (r) e 2 /m e 1/2 is the plasma frequency, v(r) = ω 2 p (r)/ω 2 , ν e f f is the ef-fective collision frequency of electrons with other plasma particles, n e (r) = n 0 + n 1 (r) is the electron density, n 0 is a homogeneous background term, n 1 (r) is a random function of position, n 1 << n 0 ; e and m e are charge and mass of an electron, respectively.⊥ , σ H , and σ || are the Pedersen, Hall's, and longitudinal conductivities, respectively, ν e , i is the electron or ion collision frequency with the neutral molecules, ω e and ω i are an electronic and an ionic angular gyrofrequencies, respectively; m e and m i are the mass of an electron and an ion, respectively.Equation (1) yields the set of equations for the COCOIMA plasma