An electrodynamic model of the solar wind interaction with the ionospheres of Mars and Venus

The electrodynamic model for the solar wind interaction with non-magnetic planets. (Cloutier and Daniell, Planet. Space Sci.21, 463, 1973; Daniell and Cloutier, Planet. Space Sci.25, 621, 1977) is modified to include the effects of non-ohmic currents in the upper ionosphere. The model is then used to calculate convection patterns induced by the solar wind in the ionospheres of Mars and Venus. For Mars the observations of the neutral mass spectrometer or Vikings 1 and 2 provided the neutral atmosphere. Model calculations reproduced the retarding potential analyzer data and indicate that the ionosphere above about 200 km is probably controlled by convection rather than chemistry or diffusion. For Venus a model atmosphere based on Dickenson and Ridley, J. Atmos. Sci.32, 1219 (1975) and Mayr et al., J. geophys. Res.83, 4411 (1978) was used. The resulting model calculations were compared to radio occultation data from Mariners 5 and 10 and Venera 9 which represent extremes in the variability of the upper Cytherean ionosphere. The model calculations are shown to fall within this variation. These results represent the state of the theory immediately prior to the Pioneer-Venus encounter.     

An empirical model for solar cycle changes in mesospheric structure at longitudes 44–77°E

Rocketsonde temperature changes during the solar cycle 1971–1982 for the height range 50–80 km at four stations, which have previously been shown by Mohanakumar [(1985) Planet. Space Sci.33, 795] to correlate with the Zürich sunspot number, are further analysed to obtain a height- and latitude-dependent formula for winter and summer seasons. The analysis confirms the strong correlation previously obtained and the magnitude of the effect but provides no evidence of a N-S hemispherical difference (as previously reported) other than that associated with season. Formulae are derived for the corresponding changes in In(pressure) and geostrophic W-E wind and are evaluated at a height above the main heating region, at 90 or so km. Changes in pressure for a change in sunspot number ΔR_n = 100 amount to 18–35% according to latitude. The corresponding changes in geostrophic W-E wind are calculated and discussed in comparison with observed zonal wind changes with solar cycle. The model presented here for solar cycle temperature changes differs markedly from that obtained by numerical modelling by Garcia et al. [(1984) Planet. Space Sci.32, 411].     

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.     

Large density fluctuations in the martian ionosphere as observed by the Mars Express radar sounder

The MARSIS (Mars Advanced Radar for Subsurface and Ionosphere Sounding) instrument on the Mars Express spacecraft provides both local and remote measurements of electron densities and measurements of magnetic fields in the martian ionosphere. The density measurements show a persistent level of large fluctuations, sometimes as much as a factor of three or more at high altitudes. Large magnetic field fluctuations are also observed in the same region. The power spectrums of both the density and magnetic field fluctuations have slopes on a log-log plot that are consistent with the Kolmogorov spectrum for isotropic fluid turbulence. The fractional density fluctuation, Δn_e/n_e, of the turbulence increases with altitude, and reaches saturation, Δn_e/n_e ∼ 1, at an altitude of about 400 km, near the nominal boundary between the ionosphere and the magnetosheath. The fluctuations are usually so large that a well-defined ionopause-like boundary between the ionosphere and the solar wind is seldom observed. Of mechanisms that could be generating this turbulence, we believe that the most likely are (1) solar wind pressure perturbations, (2) an instability in the magnetosheath plasma, such as the mirror-mode instability, or (3) the Kelvin-Helmholtz instability driven by velocity shear between the rapidly flowing magnetosheath and the ionosphere.     

Joule heating of Io's ionosphere by unipolar induction currents

The unexpectedly large scale height of Io's ionosphere (Kliore, A., et al., 1975, Icarus24, 407–410) together with the relatively large molecular weight of the likely principal constituent, SO₂ (Pearl, J., et al., 1979, Nature280, 755–758), suggest a high ionospheric temperature. Electrical induction in Io's ionosphere due to the corotating plasma bound to the Jovian magnetosphere is one possible source for attainment of such high temperatures. Accordingly, unipolar induction models were constructed to calculate ionospheric joule heating numerically. Heating rates produced by highly simplified models lie in the range 10^?9 to 10^?8 W/m³. These heating rates are lower than those determined from uv photodissociative heating models (Kumar, S., 1980, Geophys. Res. Lett.7, 9–12) at low levels in the ionosphere but are comparable in the upper ionosphere. The low electrical heating rate throughout most of the ionosphere is due to the power limitation imposed by the Alfvén wings which complete the electrical circuit (Neubauer, F.M., 1980, J. Geophys. Res.85, 1171–1178). Contrary to the pre-Voyager calculations of Cloutier, P. A., et al. (1978, Astrophys. Space Sci.55, 93–112), our numerical results show that the J × B force density due to unipolar induction currents in the ionosphere is much less than the gravitational force density when the combined mass of the neutral species is included. The binding and coupling of the ionosphere is principally due to the relatively dense (possibly localized) neutral SO₂ atmosphere. In regions where the ions and neutrals are collisionally coupled the ionosphere will not be stripped off by the J × B forces. However at a level above that (to which the ions move by diffusion only) the charged species would be removed. Thus there appears to be no need to postulate the existence of an intrinsic Ionian magnetic field as suggested by Kivelson, M. G., et al. (79, Science 205, 491–493) and Southwood, S. J., et al. (1980, J. Geophys. Res., in press) in order to retain the observed ionosphere.     

Low energy ions in co-rotating interaction regions at 1 AU: Evidence for statistical ion acceleration

The possibility of the statistical acceleration of solar wind ions to energies above 10 keV in the vicinity of co-rotating high speed solar wind streams by scattering from hydromagnetic waves is considered. We find that this process may occur only in the compressed fast stream plasma within the interaction region between the stream interface and the trailing edge, and may account for the energetic ion enhancements observed in this region by Richardson and Zwickl (Planet. Space Sci. 32, 1179, 1984). When statistical acceleration occurs in the outer heliosphere, the accelerated ions may provide a source population for acceleration at the co-rotating reverse shock.     

Dependence of the magnetopause position on the southward interplanetary magnetic field

The distance to the dayside magnetopause is statistically analyzed in order to detect the possible dependence of the dayside magnetic flux on the polarity of the interplanetary magnetic field. The effect of changing solar wind pressure is eliminated by normalizing the observed magnetopause distances by the simultaneous solar wind pressure data. It is confirmed that the normalized size of the dayside magnetosphere at the time of southward interplanetary magnetic field is smaller than that at the time of northward interplanetary magnetic field. The difference in the magnetopause position between the two interplanetary field polarity conditions ranges from 0 to 2R_E. Statistics of the relation between the magnetopause distance and the magnetic field intensity just inside the magnetopause testifies that the difference in the magnetopause position is not due to a difference in the magnetosheath plasma pressure. The effect of the southward interplanetary magnetic field is seen for all longitudes and latitudes investigated (|λ_GM|? 45°, |φ_SM|? 90°). These results strongly suggest that a part of the dayside magnetic flux is removed from the dayside at the time of southward interplanetary magnetic field.     

Solar Modulation of Cosmic Rays during the Declining and Minimum Phases of Solar Cycle 23: Comparison with Past Three Solar Cycles

We study solar modulation of galactic cosmic rays (GCRs) during the deep solar minimum, including the declining phase, of solar cycle 23 and compare the results of this unusual period with the results obtained during similar phases of the previous solar cycles 20, 21, and 22. These periods consist of two epochs each of negative and positive polarities of the heliospheric magnetic field from the north polar region of the Sun. In addition to cosmic-ray data, we utilize simultaneous solar and interplanetary plasma/field data including the tilt angle of the heliospheric current sheet. We study the relation between simultaneous variations in cosmic ray intensity and solar/interplanetary parameters during the declining and the minimum phases of cycle 23. We compare these relations with those obtained for the same phases in the three previous solar cycles. We observe certain peculiar features in cosmic ray modulation during the minimum of solar cycle 23 including the record high GCR intensity. We find, during this unusual minimum, that the correlation of GCR intensity is poor with sunspot number (correlation coefficient R=?0.41), better with interplanetary magnetic field (R=?0.66), still better with solar wind velocity (R=?0.80) and much better with the tilt angle of the heliospheric current sheet (R=?0.92). In our view, it is not the diffusion or the drift alone, but the solar wind convection that is the most likely additional effect responsible for the record high GCR intensity observed during the deep minimum of solar cycle 23.     

Solar Energetic Particle Events in the 23rd Solar Cycle: Interplanetary Magnetic Field Configuration and Statistical Relationship with Flares and CMEs

We study the influence of the large-scale interplanetary magnetic field configuration on the solar energetic particles (SEPs) as detected at different satellites near Earth and on the correlation of their peak intensities with the parent solar activity. We selected SEP events associated with X- and M-class flares at western longitudes, in order to ensure good magnetic connection to Earth. These events were classified into two categories according to the global interplanetary magnetic field (IMF) configuration present during the SEP propagation to 1 AU: standard solar wind or interplanetary coronal mass ejections (ICMEs). Our analysis shows that around 20 % of all particle events are detected when the spacecraft is immersed in an ICME. The correlation of the peak particle intensity with the projected speed of the SEP-associated coronal mass ejection is similar in the two IMF categories of proton and electron events, ≈?0.6. The SEP events within ICMEs show stronger correlation between the peak proton intensity and the soft X-ray flux of the associated solar flare, with correlation coefficient r=0.67±0.13, compared to the SEP events propagating in the standard solar wind, r=0.36±0.13. The difference is more pronounced for near-relativistic electrons. The main reason for the different correlation behavior seems to be the larger spread of the flare longitude in the SEP sample detected in the solar wind as compared to SEP events within ICMEs. We discuss to what extent observational bias, different physical processes (particle injection, transport, etc.), and the IMF configuration can influence the relationship between SEPs and coronal activity.     

Allowance for the geometry of solar wind interaction with the Earth’s magnetic field in geoeffective parameters and in the prediction of geomagnetic activity

We suggest geoeffective independent parameters that can be calculated on the basis of conventional measurements of the solar wind, which allows them to be used to forecast space weather. We present the results of our analysis of the ground variations in planetary geomagnetic activity (K _p ) and geoeffective parameters calculated on the basis of solar wind and interplanetary magnetic field measurements in the Earth’s orbit for the period 1964–1996 by taking into account the change in the orientation of the geomagnetic moment during the Earth’s diurnal and annual motions.     

Ionospheric currents induced by solar wind interaction with planetary atmospheres

In a steady-state model for the interaction of the solar wind with the atmosphere of a non-magnetic planet, the magnetized solar wind acts as a dynamo over the dayside of the planet and induces Ohmic currents in the planet's ionosphere. A model for the dynamo mechanism and for the induced current configuration is developed. Based on this model and assumed model atmospheres of Mars and Venus, the distribution of currents entering the ionosphere through the ionopause is calculated. The requirement that the total current be of such a magnitude as to cancel the shock-compressed interplanetary magnetic field fixes the ionopause altitude. The calculations for Venus are in reasonable agreement with observations. The calculations for Mars indicate the possibility of an observable ionopause in the altitude range from 325 to 425 km.     

Reconstructing CME Structures from IPS Observations Using a Phenomenological Model

We present an extension of the Tappin?–?Howard (TH) phenomenological model (Tappin and Howard, Space Sci. Rev. 147, 55, 2009) for coronal mass ejection reconstruction to use interplanetary scintillation g-map data. The necessary changes to the model are discussed. We then use the modified model to reconstruct two major interplanetary disturbances observed using the Cambridge 3.6 ha Array in September 1980. We find that despite the lower cadence of IPS observations compared with white-light imagers, a consistent reconstruction can be generated which is in agreement with in-situ measurements and solar observations.     

Cycling Changes in the Amplitudes of the 27-Day Variation of the Galactic Cosmic Ray Intensity

We study quasi-periodical changes in the amplitudes of the 27-day variation of the galactic cosmic ray (GCR) intensity, and the parameters of solar wind and solar activity. We have recently found quasi-periodicity of three to four Carrington rotation periods (3?–?4 CRP) in the amplitudes of the 27-day variation of the GCR intensity (Gil and Alania in J. Atmos. Solar-Terr. Phys. 73, 294, 2011). A similar recurrence is recognized in parameters of solar activity (sunspot number, solar radio flux) and solar wind (components of the interplanetary magnetic field, solar wind velocity). We believe that the 3?–?4 CRP periodicity, among other periodicities, observed in the amplitudes of the 27-day variation of the GCR intensity is caused by a specific cycling structure of the Sun’s magnetic field, which may originate from the turbulent nature of the solar dynamo.     

Heliolatitude and Time Variations of Solar Wind Structure from in situ Measurements and Interplanetary Scintillation Observations

The 3D structure of the solar wind and its evolution in time are needed for heliospheric modeling and interpretation of energetic neutral atoms observations. We present a model to retrieve the solar wind structure in heliolatitude and time using all available and complementary data sources. We determine the heliolatitude structure of solar wind speed on a yearly time grid over the past 1.5 solar cycles based on remote-sensing observations of interplanetary scintillations, in situ out-of-ecliptic measurements from Ulysses, and in situ in-ecliptic measurements from the OMNI 2 database. Since in situ out-of-ecliptic information on the solar wind density structure is not available apart from the Ulysses data, we derive correlation formulae between the solar wind speed and density and use the information on the solar wind speed from interplanetary scintillation observations to retrieve the 3D structure of the solar wind density. With the variations of solar wind density and speed in time and heliolatitude available, we calculate variations in solar wind flux, dynamic pressure, and charge-exchange rate in the approximation of stationary H atoms.     

On the Relationship of the 27-day Variations of the Solar Wind Velocity and Galactic Cosmic Ray Intensity in Minimum Epoch of Solar Activity

We study the relationship of the 27-day variations of the galactic cosmic ray intensity with similar variations of the solar wind velocity and the interplanetary magnetic field based on observational data for the Bartels rotation period # 2379 of 23 November 2007 – 19 December 2007. We develop a three-dimensional (3-D) model of the 27-day variation of galactic cosmic ray intensity based on the heliolongitudinally dependent solar wind velocity. A consistent, divergence-free interplanetary magnetic field is derived by solving Maxwell’s equations with a heliolongitudinally dependent 27-day variation of the solar wind velocity reproducing in situ observations. We consider two types of 3-D models of the 27-day variation of galactic cosmic ray intensity, i) with a plane heliospheric neutral sheet, and ii) with the sector structure of the interplanetary magnetic field. The theoretical calculations show that the sector structure does not significantly influence the 27-day variation of galactic cosmic ray intensity, as had been shown before, based on observational data. Furthermore, good agreement is found between the time profiles of the theoretically expected and experimentally obtained first harmonic waves of the 27-day variation of the galactic cosmic ray intensity (with a correlation coefficient of 0.98±0.02). The expected 27-day variation of the galactic cosmic ray intensity is inversely correlated with the modulation parameter ζ (with a correlation coefficient of −0.91±0.05), which is proportional to the product of the solar wind velocity V and the strength of the interplanetary magnetic field B (ζ∼VB). The high anticorrelation between these quantities indicates that the predicted 27-day variation of the galactic cosmic ray intensity mainly is caused by this basic modulation effect.     

Wavelet Analysis on Solar Wind Parameters and Geomagnetic Indices

The geomagnetic activity is the result of the solar wind–magnetosphere interaction. It varies following the basic 11-year solar cycle; yet shorter time-scale variations appear intermittently. We study the quasi-periodic behavior of the characteristics of solar wind (speed, temperature, pressure, density) and the interplanetary magnetic field (B _x , B _y , B _z , β, Alfvén Mach number) and the variations of the geomagnetic activity indices (D _ST, AE, A _p and K _p). In the analysis of the corresponding 14 time series, which span four solar cycles (1966?–?2010), we use both a wavelet expansion and the Lomb/Scargle periodograms. Our results verify intermittent periodicities in our time-series data, which correspond to already known solar activity variations on timescales shorter than the sunspot cycle; some of these are shared between the solar wind parameters and geomagnetic indices.     

Effects of interplanetary magnetic field and magnetospheric substorm variations on the dayside aurora

Photometric observations of dayside auroras are compared with simultaneous measurements of geomagnetic disturbances from meridian chains of stations on the dayside and on the nightside to document the dynamics of dayside auroras in relation to local and global disturbances. These observations are related to measurements of the interplanetary magnetic field (IMF) from the satellites ISEE-1 and 3. It is shown that the dayside auroral zone shifts equatorward and poleward with the growth and decay of the circum-oval/polar cap geomagnetic disturbance and with negative and positive changes in the north-south component of the interplanetary magnetic field (B_z). The geomagnetic disturbance associated with the auroral shift is identified as the DP2 mode. In the post-noon sector the horizontal disturbance vector of the geomagnetic field changes from southward to northward with decreasing latitude, thereby changing sign near the center of the oval precipitation region. Discrete auroral forms are observed close to or equatorward of the ΔH = 0 line which separates positive and negative H-component deflections. This reversal moves in latitude with the aurora and it probably reflects a transition of the electric field direction at the polar cap boundary. Thus, the discrete auroral forms observed on the dayside are in the region of sunward-convecting field lines. A model is proposed to explain the equatorward and poleward movement of the dayside oval in terms of a dayside current system which is intensified by a southward movement of the IMF vector. According to this model, the Pedersen component of the ionospheric current is connected with the magnetopause boundary layer via field-aligned current (FAC) sheets. Enhanced current intensity, corresponding to southward auroral shift, is consistent with increased energy extraction from the solar wind. In this way the observed association of DP2 current system variations and auroral oval expansion/contraction is explained as an effect of a global, ‘direct’ response of the electromagnetic state of the magnetosphere due to the influence of the solar wind magnetic field. Estimates of electric field, current, and the rate of Joule heat dissipation in the polar cap ionosphere are obtained from the model.     

Flapping motions of the tail plasma sheet induced by the interplanetary magnetic field variations

Flapping motions of the magnetotail with an amplitude of several earth radii are studied by analysing the observations made in the near (x = ?25 ~ ?30 R_E and the distant (x? ?60 R_E) tail regions. It is found that the flapping motions result from fluctuations in the interplanetary magnetic field, especially Alfvénic fluctuations, when the magnitude of the interplanetary magnetic field is larger than ~10 γ and they propagate behind the Earth with the solar wind flow. Flappings tend to be observed in early phases of the magnetospheric substorm, and they have two fundamental modes with periods of ~200 and ~500 sec. In some limited cases a good correspondence with the long period micropulsations (Pc5) in the polar cap region is observed. These observational results are explained by the model in which the Alfvénic fluctuations in the solar wind penetrate into the magnetosphere along the connected interplanetary-magnetospheric field lines. The characteristics of the flapping reveal that the geomagnetic tail is a good resonator for the hydromagnetic disturbances in the solar wind.     

Dependence of the amount of open magnetic flux on the direction of the interplanetary magnetic field

The power generated by the solar wind—magnetosphere dynamo is proportional to the amount of the open magnetic flux Φ. It is difficult to use this fact in determining observationally the dependence of Φ on the orientation of the interplanetary magnetic field vector. It is shown that, for a simple vacuum superposition of the earth's dipole field and a uniform magnetic field, Φ is very closely proportional to sin θ/2) for a wide range of the intensity of the uniform field, where θ denotes the polar angle of the interplanetary magnetic field vector in the Y-Z plane of solar-magnetospheric coordinates.