icsinPpr|Dsuso ae icsinPpr|Dsuso ae | Paper Discussion | Paper Discussion | Paper Discussion | Paper Discussion Nonlin. Processes Geophys. Discuss., 1, 1431–1464, 2014 www.nonlin-processes-geophys-discuss.net/1/1431/2014/ doi:10.5194/npgd-1-1431-2014 NPGD © Author(s) 2014. CC Attribution 3.0 License. 1, 1431–1464, 2014
This discussion paper is/has been under review for the journal Nonlinear Processes Self-organization of in Geophysics (NPG). Please refer to the corresponding final paper in NPG if available. ULF electromagnetic wave structures
Self-organization of ULF electromagnetic G. Aburjania et al. wave structures in the shear flow driven dissipative ionosphere Title Page Abstract Introduction G. Aburjania1,2, K. Chargazia1,2, O. Kharshiladze2, and G. Zimbardo3 Conclusions References 1 I. Vekua Institute of Applied Mathematics, Tbilisi State University, 2 University str., Tables Figures 0143 Tbilisi, Georgia 2M. Nodia Institute of Geophysics, Tbilisi State University, 1 Aleqsidze str., 0193 Tbilisi, Georgia J I 3 Physics Department, University of Calabria, Ponte P. Bucci, Cubo 31 C, 87036 Rende, Italy J I
Received: 18 July 2014 – Accepted: 28 July 2014 – Published: 26 August 2014 Back Close
Correspondence to: K. Chargazia ([email protected]) Full Screen / Esc and O. Kharshiladze ([email protected])
Published by Copernicus Publications on behalf of the European Geosciences Union & the Printer-friendly Version American Geophysical Union. Interactive Discussion
1431 icsinPpr|Dsuso ae icsinPpr|Dsuso ae | Paper Discussion | Paper Discussion | Paper Discussion | Paper Discussion Abstract NPGD This work is devoted to investigation of nonlinear dynamics of planetary electromag- netic (EM) ultra-low-frequency wave (ULFW) structures in the rotating dissipative iono- 1, 1431–1464, 2014 sphere in the presence of inhomogeneous zonal wind (shear flow). Planetary EM 5 ULFW appears as a result of interaction of the ionospheric medium with the spatially in- Self-organization of homogeneous geomagnetic field. The shear flow driven wave perturbations effectively ULF electromagnetic extract energy of the shear flow increasing own amplitude and energy. These pertur- wave structures bations undergo self organization in the form of the nonlinear solitary vortex struc- tures due to nonlinear twisting of the perturbation’s front. Depending on the features G. Aburjania et al. 10 of the velocity profiles of the shear flows the nonlinear vortex structures can be either monopole vortices, or dipole vortex, or vortex streets and vortex chains. From ana- lytical calculation and plots we note that the formation of stationary nonlinear vortex Title Page structure requires some threshold value of translation velocity for both non-dissipation Abstract Introduction and dissipation complex ionospheric plasma. The space and time attenuation specifi- Conclusions References 15 cation of the vortices is studied. The characteristic time of vortex longevity in dissipative ionosphere is estimated. The long-lived vortices transfer the trapped medium particles, Tables Figures energy and heat. Thus they represent structural elements of turbulence in the iono-
sphere. J I
J I 1 Introduction Back Close
20 The interaction between the solar wind (SW) and terrestrial magnetosphere is the Full Screen / Esc primary driver of many processes and phenomena occurring in the magnetosphere
and consequently, in the ionosphere. However, many of the energy transfer processes Printer-friendly Version have a sporadic/bursty character, and observations have highlighted that the vortex- like plasma flow structures are common in the Earth’s magnetosphere. They prevail Interactive Discussion 25 in the nightside plasma sheet (Hones, 1978, 1981, 1983; Keiling et al., 2009) and at the flank magnetospheric low-latitude boundary layers (LLBLs) (Fairfield et al., 2000; 1432 icsinPpr|Dsuso ae icsinPpr|Dsuso ae | Paper Discussion | Paper Discussion | Paper Discussion | Paper Discussion Otto and Fairfield, 2000; Hasegawa et al., 2004). Plasma vortex-like flows have also been observed on the middle to high-latitude boundary of the outer radiation belt by NPGD the Cluster spacecraft fleet (Zong et al., 2009). The plasma flow vortex found in the 1, 1431–1464, 2014 magnetotail is characterized by pronounced vortical motion in the plane that is approx- 5 imately parallel to the ecliptic plane. Vortex-like plasma flows in the plasma sheet are thought to be important in transportation of the kinetic energy from fast flow or bursty Self-organization of bulk flows (BBFs) in the magnetotail to the near Earth region (Snekvik et al., 2007; Keil- ULF electromagnetic ing et al., 2009). A sudden braking and/or azimuthal deflection of BBFs may generate wave structures the plasma flow vortices at the boundary between the magnetotail plasma sheet and G. Aburjania et al. 10 the inner magnetosphere as suggested by Hasegawa (1979) and Vasyliunas (1984). Keika et al. (2009) have shown that plasma vortices are formed near the region where the earthward flows slow down and turn in azimuthal directions. The theoretical rela- Title Page tion between the field-aligned current (FAC) and plasma vorticity showed that FACs are generated to transport transverse momentum along magnetic field lines (Pritchett Abstract Introduction 15 and Coroniti, 2000). Keiling et al. (2009) presented a scenario in which the plasma flow Conclusions References vortices in the plasma sheet generated the FAC of the substorm current wedge (SCW) at the beginning of the substorm expansion phase and coupled to the ionosphere, Tables Figures causing the ionospheric vortices. In the present paper, we continue investigation of a special type of internal waves, J I 20 which appear in the ionosphere under the influence of the spatially inhomogeneous geomagnetic field and the Earth’s rotation velocity (Aburjania et al., 2002, 2003, 2004, J I 2007). In the previous works generation mechanism and self-organization into the vor- Back Close tex structures is studied. Now, we are interested in large-scale (planetary) ultra-low- Full Screen / Esc frequency (ULF) electromagnetic (EM) vortex structures in the ionospheric medium 25 (consisting of electrons, ions and neutral particles), which have a horizontal linear scale 3 Printer-friendly Version Lh of order 10 km and higher, a vertical scale Lv of altitude scale order H(Lv ≈ H). Observations (Gershman, 1974; Gossard and Hooke, 1975; Kamide and Chian, Interactive Discussion 2007) show also, that spatially inhomogeneous zonal winds (shear flows), produced by nonuniform heating of the atmospheric layers by solar radiation, permanently exist
1433 icsinPpr|Dsuso ae icsinPpr|Dsuso ae | Paper Discussion | Paper Discussion | Paper Discussion | Paper Discussion in the atmosphere and ionosphere layers. Herewith, investigation of the problem of gen- eration and evolution of ionospheric EM ULF electromagnetic wave structures taking NPGD into account the inhomogeneous zonal wind (shear flow) becomes important. 1, 1431–1464, 2014
2 The governing equations Self-organization of ULF electromagnetic 5 We choose our model as a two-dimensional β-plane with sheared flow. Since the 3 wave structures length of planetary waves (λ ≥ 10 km) is comparable with the Earth’s radius R, we investigate such notions in approximation of the β-plane, which was specially devel- G. Aburjania et al. oped for analysis of large-scale processes (Pedlosky, 1978), in the “standard” coor- dinate system. In this system, the x axis is directed along the parallel to the east, 10 the y axis along the meridian to the north and the z axis – vertically upwards (the Title Page local Cartesian system). For simplicity, the equilibrium velocity V 0, geomagnetic field Abstract Introduction H0, perturbed magnetic field h and frequency of Earth’s rotation Ω0 are given by for- mulas V 0 = V 0(y)ex, H0(0,0,−Hp cosθ), h(0,0,hz), Ω0(0,0,Ω0 cosθ). Here and else- Conclusions References e e e e H −5T where ( x, y , z) denotes a unit vector, p = 5 × 10 is the value of geomagnetic Tables Figures 15 field strength in the pole and we suppose that geomagnetic colatitude θ coincides with a geographical colatitude θ0. In the ionosphere the large-scale motions are quasi- J I horizontal (two-dimensional) (Aburjania et al., 2002, 2003, 2006) and hydrodynamic velocity of the particles V = (Vx,,Vy ,0). The fluid is assumed to be incompressible and J I ψ V ψ e therefore a stream function can be defined trough = [∇ , z]. Medium motion is Back Close 20 considered near the latitude ϕ0 = π/2 − θ0. Not considering any more detail in the new under review branches of planetary waves Full Screen / Esc (see Aburjania et al., 2002, 2003, 2004, 2007, 2011) we would like to note that be- ginning with the altitude of 80 km and higher, the upper atmosphere of the Earth is Printer-friendly Version
a strongly dissipative medium. Often when modelling large-scale processes for this re- Interactive Discussion 25 gion of the upper atmosphere, effective coefficient of Rayleigh friction between the iono- spheric layers is introduced. The role of the ion friction rapidly increases at the altitudes above 120 km (Kelley, 1989; Kamide and Chian, 2007) and its analytical expression 1434 icsinPpr|Dsuso ae icsinPpr|Dsuso ae | Paper Discussion | Paper Discussion | Paper Discussion | Paper Discussion coincides with the Rayleigh friction formula (Aburjania and Chargazia, 2007). There- fore, often during a study of large-scale (103–104) km, ULF (10–10−6) s−1 wavy struc- NPGD tures in the ionosphere, we will apply the well-known Rayleigh formula to dissipative 1, 1431–1464, 2014 force F = −ΛV , assuming the altitudes above (80–130) km Λ ≈ 10−5 s−1 (Dickinson, 5 1969; Gosard and Hooke, 1975), and the altitudes above 130 km Λ = Nνin/Nn, where N and Nn denote concentrations of the charged particles and neutral particles, νin is fre- Self-organization of quency of collision of ions with molecules (Gershman, 1974; Kelley, 1989; Al’perovich ULF electromagnetic and Fedorov, 2007). wave structures The governing equations of the considered problem are the closed system of magne- G. Aburjania et al. 10 tohydrodynamic equations of the electrically conducting ionosphere (Gershman, 1974; Kelley, 1989; Aburjania et al., 2004, 2007; Al’perovich and Fedorov, 2007). The solution of the temporal evolution of inhomogeneously sheared flow reduces to solution of the Title Page set of nonlinear partial differential equations for ψ and magnetic field perturbation, hz (see Aburjania et al., 2002): Abstract Introduction ∂ ∂ 00 ∂ψ ∂h Conclusions References 15 + V (y) ∆ψ + β − V + C + Λ∆ψ = J(ψ,∆ψ), (1) ∂t 0 ∂x 0 ∂x H ∂x Tables Figures ∂ ∂ ∂ψ ∂h + V (y) h − β + δ · C = J(ψ,h). (2) ∂t 0 ∂x H ∂x H ∂x J I Here J I ∂2Ω0 1 ∂ 2Ω0 sinθ0 β = = − (2Ω ) = , Back Close ∂y R ∂θ 0 R eH 2 Full Screen / Esc eN ∂H0z N p 00 d V0(y) 20 β = = − sinθ < 0, V (y) = , H ρc ∂y N MRc 0 0 d 2y n Printer-friendly Version eN c ∂H0z cHp h = hz, CH = = − sinθ0 < 0, Interactive Discussion NnMc 4πeN ∂y 4πeNR ∂2 ∂2 ∂a ∂b ∂a ∂b ∆ = + , J(a,b) = · − · . (3) ∂x2 ∂y2 ∂x ∂y ∂y ∂x 1435 icsinPpr|Dsuso ae icsinPpr|Dsuso ae | Paper Discussion | Paper Discussion | Paper Discussion | Paper Discussion
ρ = NnM is density of neutral particles; m and M are masses of electrons and ions NPGD (molecules); e is the magnitude of the electron charge; c is the light speed. Further 1, 1431–1464, 2014 we consider a motion in neighborhood of fixed latitude (θ = θ0). The dimensionless pa- 5 rameter δ is introduced here for convenience. In the ionospheric E region (80–150) km, where the Hall effect plays an important role, this parameter is equal to unity (δ = 1). Self-organization of In the F region (200–600) km, where the Hall effect is absent, δ turns to zero (δ = 0). ULF electromagnetic The system of Eqs. (1) and (2), at corresponding initial and boundary conditions, wave structures describes nonlinear evolution of the spatial two-dimensional large-scale ULF electro- G. Aburjania et al. 10 magnetic perturbations in sheared incompressible ionospheric E and F regions. From the Eqs. (1) and (2) we determine the temporal evolution of the energy of wavy structures, E(x,y,t) Title Page Z Z ∂E ∂ψ ∂ψ 2 = V 0(y) dxdy − Λ |∇ψ| dxdy. (4) Abstract Introduction ∂t 0 ∂x ∂y Conclusions References 15 We note that in the absence of zonal flow (V0 = 0) and Rayleigh friction (Λ = 0) the wavy structure energy is conserved. Tables Figures Therefore, the existence of sheared zonal flow can be considered as the presence J I of an external energy source. One can see, that it (term with V0(y) in Eq. 4) feeds the medium with external source of energy for generation of the wave structures (devel- J I 20 opment of the shear flow instability). The shear flows can become unstable transiently until the condition of the strong relationship between the shear flows and wave per- Back Close turbations is satisfied (Chagelishvili et al., 1996; Aburjania et al., 2006), e. i. the per- Full Screen / Esc turbation falls into amplification region in the wave number space. Leaving this region, e. i. when the perturbation passes to the damping region in the wave vector space, Printer-friendly Version 25 it returns an energy to the shear flow and so on (if the nonlinear processes and self- organization of the vortex structure will not develop before) (Aburjania et al., 2006). Interactive Discussion The experimental and observation data shows the same (Gossard and Hooke, 1975; Pedlosky, 1987; Gill, 1982). 1436 icsinPpr|Dsuso ae icsinPpr|Dsuso ae | Paper Discussion | Paper Discussion | Paper Discussion | Paper Discussion 3 Criterion of instability of ULF PEMW wave structures in the ionopashere with the shear flow NPGD
The character of the shear flow much defines the evolution of the wave perturbation in 1, 1431–1464, 2014 the medium. Therefore, the shear flow in the hydrodynamics and in magnetohydrody- 5 namics are often unstable (Gossard and Hooke, 1975; Mikhailovskii, 1974; Timofeev, Self-organization of 00 2 2 2000). Existing of the term proportional to V0 = d V0/dy in Eq. (1) is related with the ULF electromagnetic criterion (condition) of instability of the shear flow. In the linear approximation for the wave structures small scale perturbations of the form ψ(x,y,t),h(x,y,t) = (ψ1(y),h1(y))exp(ikxx − iω) from the Eqs. (1) and (2) follows the equation of Ohr–Zomerfeld G. Aburjania et al.
2 ! 2 ! k β − V 00 Λ d 2 d 2 x 0 10 i − k ψ + − k ψ − ψ 2 x 1 2 x 1 1 Title Page ω − kxV0 dy dy ω − kxV0 k2C β Abstract Introduction + x H H ψ = 0. (5) 1 Conclusions References (ω − kxV0)[ω − kx(CH + V0)] Neglecting the dissipation terms (Λ → 0) from Eq. (5) we get Tables Figures