Steady state charge neutral models of the magnetopause

Plane models of the magnetopause are investigated under the assumption that ionospheric electrons are able to short-circuit electric fields (exact charge neutrality). Using the Vlasov theory a general method is presented for constructing distribution functions that lead to given magnetic field and tangential bulk velocity profiles. As an example we describe the magnetic field transition in terms of error functions and obtain particle distributions in explicit form, including bulk velocities.It is thus shown that bulk velocities in the direction of the magnetic field do not necessarily lead to a non-equilibrium magnetopause which investigations by Parker and Lerche seem to suggest.Of the European Space Research Organisation (ESRO).     

Force-free equilibrium at magnetic neutral points

It is shown that, at neutral points of force-free magnetic fields, the electric current density must vanish. This property is independent of whether the neutral points are isolated or (e.g.) fill lines or surfaces. One implication is the fact that in a cold pressureless plasma the formation of neutral current sheets cannot be adiabatically slow. The field-line topology in the neighbourhood of neutral points is discussed. At neutral points of force-free magnetic fields in general three constant- surfaces, defined by the equation ×B=B, with the same value of intersect orthogonally. If, during a time-development, the magnetic field gradient matrix B _i/x _j becomes singular at a neutral point, the field topology can change qualitatively — in general connected with the merger of two or more neutral points into one and/or the splitting up of one neutral point into several others. This can be interpreted as implying the transition from a quasi-static evolution to a dynamical state in which magnetic energy is released.     

Stability of the magnetopause

The Magnetosphere-Solar Wind boundary is treated as a tangential discontinuity between collisionless plasma described by C.G.L. equations. An estimate is made, analytically, of the minimum shear speed required to render the interface unstable for propagation in any arbitrary direction. It is found that even systems, which are, for all shear speeds, stable towards propagation parallel to the field, are rendered unstable by propagation away from the field, in some velocity domain. Numerical evaluations for parameters characteristic of the Solar wind-Magnetosphere boundary show that the interface could be unstable even under relatively quiet conditions.     

The structure of the magnetopause

The solar wind is a magnetized flowing plasma that intersects the Earth's magnetosphere at a velocity much greater than that of the compressional fast mode wave that is required to deflect that flow. A bow shock forms that alters the properties of the plasma and slows the flow, enabling continued evolution of the properties of the flow on route to its intersection with the magnetopause. Thus the plasma conditions at the magnetopause can be quite unlike those in the solar wind. The boundary between this “magnetosheath” plasma and the magnetospheric plasma is many gyroradii thick and is surrounded by several boundary layers. A very important process occurring at the magnetopause is reconnection whereby there is a topological change in magnetic flux lines so that field lines can connect the solar wind plasma to the terrestrial plasma, enabling the two to mix. This connection has important consequences for momentum transfer from the solar wind to the magnetosphere. The initiation of reconnection appears to be at locations where the magnetic fields on either side of the magnetopause are antiparallel. This condition is equivalent to there being no guide field in the reconnection region, so at the reconnection point there is truly a magnetic neutral or null point. Lastly reconnection can be spatially and temporally varying, causing the region of the magnetopause to be quite dynamic.     

Wave-coupling at the magnetopause

The paper considers wave coupling for an arbitrary direction of propagation on the basis of single fluid hydromagnetic equations appropriate for a rarefied plasma. The analysis is used to study the transfer of solar wind momenta into the magnetosphere. It is found that wave refraction is significant only during disturbed conditions for waves travelling with the wind. Enhanced reflection of waves might be important even under quiet conditions in the flanks of the magnetosphere.     

On equilibrium rotation of the central object of an accretion disc

The hydrodynamic interaction of an accretion disc with its central object is reanalysed within the framework of the slim-disc approximation. Arguments are presented against an interpretation of the total angular momentum flux as an eigenvalue of the system. A simple intuitive consideration is provided, which shows that the central object may be in a state of stationary rotation even if the disc imposes the constraint of a finite angular momentum flux into it. It is argued that equilibrium rotation is characterized by vanishing viscous torque rather than by zero total angular momentum flux. As a consequence, the central object can be in a state of stationary rotation below the break-up limit, although its angular momentum increases. Despite accretion, even for positive total angular momentum flux and subcritical rotation, central objects are spun down within a considerable range of their parameters. The results are illustrated by application to FU Orionis systems.     

On some problems concerning the calculation of the form of the magnetopause and the electric field on the magnetopause in the framework of the continuum medium model

In order to understand the reason of the existence of the electric field in the magnetosphere, and for the theoretical evaluation of its value, it is necessary to find the solution of the problem of determination of the magnetosphere boundary form in the frameworks of the continuum medium model which takes into account part of the magnetospheric plasma movement in supporting the magnetospheric boundary equilibrium. A number of problems for finding the distribution of the pressure, the density, the magnetic field and the electric field on the particular tangential discontinuity is considered in the case when the form of discontinuity is set (the direct problem) and a number of problems for finding the form of the discontinuity and the distribution of the above-mentioned physical quantities on the discontinuity is considered when the law of the change of the external pressure along the boundary is set (for example, with the help of the approximate Newton equation). The problem which is considered here, which deals with the calculation of the boundary form and with the calculation of the distribution of the corresponding physical quantities on the discontinuity of the 1st kind for the compressible fluid with the magnetic field with field lines which are perpendicular to the plane of the flow in question, concerns the last sort of problems. The comparison of the results of the calculation with the data in the equatorial cross-section of the magnetosphere demonstrates that the calculated form of the boundary, the value of the velocity of the return flow and the value of the electric field on the magnetopause, agree satisfactorily with the observational data.     

The three-dimensional shape of the magnetopause

Nearly 1000 magnetopause crossings from HEOS-2, HEOS-1, OGO-5 and 5 IMP space-craft covering most of the northern and part of the southern dayside and near-Earth tail magnetopause (X >?15 R_E) have been used to perform a detailed study of the three-dimensional shape and location of the magnetopause. The long-term influence of the solar wind conditions on the average magnetopause geometry has been reduced by normalising the radial distances of the observed magnetopause crossings to an average dynamical solar wind pressure. Best-fit ellipsoids have been obtained to represent the average magnetopause surface in geocentric solar ecliptic (GSE) and (as a function of tilt angle) in solar magnetic (SM) coordinates. Average geocentric distances to the magnetopause for the 1972–1973 solar wind conditions (density 9.4 cm^?3, velocity 450 km s^?1) are 8.8 R_E in the sunward direction, 14.7 R_E in the dusk direction, 13.4 R_E in the dawn direction and 13.7 R_E in the direction normal to the ecliptic plane. The magnetopause surface is tilted by 6.6° ± 2° in a direction consistent with that expected from the aberration effect of the radial solar wind. Our data suggest that the solar wind plasma density and the interplanetary magnetic field (IMF) orientation affect the distance to the polar magnetopause, larger distances corresponding to higher plasma density and southward fields. Our best-fit magnetopause surface shows larger geocentric distances than predicted by the model of Choe et al. [Planet Space Sci. 21, 485 (1973).] normalised to the same solar wind pressure.     

The stability of the pristine magnetopause

A MHD theory of combined Kelvin-Helmholtz (KH) and Rayleigh-Taylor (RT) instabilities for a transition layer with two different scale lengths (Δ and δ for the variation of velocity/magnetic fields and density, respectively) is presented. The study is motivated by reports of magnetopauses with no low latitude boundary layer, in which a sharp density drop over a distance δ?Δ is observed (“pristine” magnetopauses (J. Geophys. Res. 101 (1996) 49). The theory ignores compressibility effects and applies to subsonic regions of the dayside magnetopause. The RT effect is included to account for temporary periods of acceleration of the magnetopause, caused by sudden changes of the solar wind dynamic pressure. For small wavelengths λ, such that δ?λ?Δ, a WKB solution shows that the velocity gradient operates, together with magnetic tensions, to attenuate or even stabilize the Rayleigh-Taylor instability within a certain wavelength range. An exact dispersion relation for flute modes, valid for all λ, in the form of a fourth order polynomial for the complex frequency ω, is obtained from a model with a constant velocity gradient, dV/dy within Δ, and with δ→0. Flute modes are possible because of the existence of bands of very small magnetic shear on the dayside magnetopause (J. Geophys. Res. 103 (1998) 6703). The exact solution allows for a study of the change of the action of the velocity gradient with λ from the long-λ range where dV/dy is KH destabilizing to the short-λ range where dV/dy produces a stabilizing effect. Both, the WKB approximation and the well known tangential discontinuity model (Δ→0) are recovered as limiting cases of the exact solution. Properties of the KH and RT instabilities, for different density ratios on either side of the magnetopause, are described. For flute modes, at very small λ the RT instability grows faster and becomes the dominant effect. However, it is shown that the growth rate remains bounded at a finite value as λ→0, when a theory with a finite δ model is considered. To study configurations with finite, arbitrary, δ/Δ ratios, the MHD perturbation equations are solved numerically, using hyperbolic tangent functions for both the density and velocity transitions across the magnetopause. To examine the influence of different δ/Δ ratios on the growth rates of KH and RT, calculations are performed for different δ/Δ, with and without acceleration, and for two different density ratios. It is found that the general features exhibited by the constant dV/dy model, are confirmed by these numerical solutions. The stability of pristine magnetopauses, and the possibility of observing some theoretical predictions during magnetopause crossings in ongoing missions, are discussed.     

Two-dimensional equilibrium of a neutral sheet in a gravitational field

We solve the two-dimensional MHD-equations to find the equilibrium structure of a neutral sheet having its axis parallel to the gravitational field in an atmosphere. The evolution of an initially plane sheet which is a self-consistent, non-equilibrium solution of the MHD-equations, is followed until static equilibrium of forces is obtained. The effect of field line anchoring in a cool dense layer at the bottom of the sheet is taken into account.Presently at The Auroral Observatory, Boks 953, 9001 Tromsø, Norway.     

The Kelvin-Helmholtz instability on the magnetopause

Conditions for the development of Kelvin-Helmholtz (K-H) waves on the magnetopause have been known for more than 15 years; more recently, spacecraft observations have stimulated further examination of the properties of K-H waves. For amagnetopause with no boundary layer, two different modes of surface waves have been identified and their properties have been investigated for various assumed orientations of magnetic field and flow velocity vectors. The power radiated into the magnetosphere from the velocity shear at the boundary has been estimated. Other calculations have focused on the consequences of finite thickness boundary layers, both uniform and non-uniform. The boundary layer is found to modify the wave modes present at the magnetopause and to yield a criterion for the wavelength of the fastest growing surface waves. The paper concludes by questioning the extent to which the inferences from boundary layer models are model dependent and identifies areas where further work is needed or anticipated.     

On the MHD instability of the Earth's magnetopause and its geophysical effects

This paper is concerned with the Kelvin-Helmholtz magnetohydrodynamic (MHD) instability near the low-latitude boundary layer. It is argued that the instability may be responsible both for viscous interaction of solar wind with the magnetosphere and for generation of surface waves over the range of geomagnetic pulsation frequencies. The influence of the inclination angle of the IMF vector with respect to the Earth-Sun line upon the instability growth-rate value is studied numerically. On the basis of the results the morning maximum of surface waves and of geomagnetic pulsations as well as the dependence of the latter on IMF orientation are discussed.     

Tearing modes in the magnetopause current sheet

A simple method is proposed to investigate the stability of a charge neutral magnetopause current sheet with respect to the tearing-mode instability. This method may serve as a useful tool in understanding the processes of local opening of the closed magnetosphere.     

On the equilibrium figures of an ideal rotating liquid in the post-Newtonian approximation of general relativity

Maclaurin's P-ellipsoid, which is an equilibrium figure in the post-Newtonian approximation of general relativity, is constructed in the neighbourhood of Maclaurin's classical ellipsoid. Its shape and rotation velocity are investigated. It is shown that in the case of a P-ellipsoid with the mass and the angular momentum of a basic ellipsoid the effects of general relativity reduce the rotation velocity and decrease its volume.     

On the equilibrium figures of an ideal rotating fluid in the post-newtonian approximation of general relativity

A stability criterion is given for the equilibrium form of an ideal rotating fluid in the post-Newtonian approximation. This generalizes the known Lyapunov criterion in classical dynamics. The sphere stability is also investigated and it is shown that it is stable only whenR>22.2R _g (R is the relativistic sphere radius,R _g the Schwarzschild radius).     

On neutral sheets in the solar wind

The structure and dynamics of neutral sheets in the solar wind is examined. The internal magnetic topology of the sheet is argued to be that of thin magnetic tongues greatly distended outward by the expansion inside the sheet. Due to finite conductivity effects, outward flow takes place across field lines but is retarded relative to the ambient solar wind by the reverse J×B force. The sheet thickness as well as the internal transverse magnetic field are found to be proportional to the electrical conductivity to the inverse one third power. Estimating a conductivity appropriate for a current carried largely by the ions perpendicular to the magnetic field, we find sheet dimensions of the order of 500km representative for the inner solar corona. For a radial field of strength 1/2G at 2R , the transverse field there is about 2 × 10^–3G and decreases outward rapidly.The energy release in the form of Joulean dissipation inside the sheet is estimated. It is concluded that ohmic heating in current sheets is not a significant source of energy for the overall solar wind expansion, mainly because these structures occupy only a small percentage of the total coronal volume. However, the local energy release through this mechanism is found to be large - in fact, over 7 times that expected to be supplied by thermal conduction. Therefore, ohmic heating is probably a dominant energy source for the dynamical conditions within the sheet itself.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Interaction of interplanetary MHD shock waves with the magnetopause

The interaction between a shock-wave and the magnetopause is formulated on the basis of one-dimensional magnetohydrodynamics. The magnetopause is assumed to be a tangential discontinuity, and the magnetic field is limited to the case of perpendicularity. Both the forward and reverse shocks' impact on the magnetopause are considered and analyzed separately. The forward shock-magnetopause interaction results in a transmitted shock, a tangential discontinuity, and a simple rarefaction wave. The reverse shock-magnetopause interaction creates a transmitted shock, a tangential discontinuity, and a reflected wave. The propagation of an SSC signal which is related to an interplanetary shock-induced geomagnetic storm's onset-time on Earth is discussed in general terms. It was found in earlier work (Shen and Dryer, 1972) that the propagation velocity of an inter-planetary shock is decreased by about 1015% following its impact with the earth's bow shock; the present study shows that its velocity is then suddenly increased by a factor of two to three after impact with the magnetopause. The fast propagating shock-wave inside the magnetosphere degenerates into a hydromagnetic wave as it advances into an increasing intensity of the distorted dipole geomagnetic field.     

Reflection and refraction of hydromagnetk waves at the magnetopause

Reflection and transmission coefficients of MHD waves are obtained at a stable, plane interface which separates two compressible, perfectly conducting media in relative motion to each other. The coefficients are evaluated for representative conditions of the quiettime, near-Earth magnetopause. The transmission coefficient averaged over a hemispherical distribution of incident waves is found to be 1–2 per cent. Yet the magnitude of the energy flux deposited into the magnetosphere in a day averaged over a hemispherical distribution of waves having amplitudes of say 2–3 gamma, is estimated to be of the order 10²² erg. Therefore the energy input of MHD waves must contribute significantly to the energy budget of the magnetosphere. The assumption that the boundary surface is a tangential discontinuity with no curvature limits the present theory to hydromagnetic frequencies higher than about 10⁻¹ Hz. The ion gyrofrequencies for the models assumed here lie above 2 × 10⁻¹ Hz. Therefore the present treatment applies to MHD waves near 10⁻¹ Hz.