Interaction of energetic particles with HM-waves in the magnetosphere

Energetic particle response in electromagnetic fields of ULF HM-waves in the magnetosphere is reviewed. Pc4–5 geomagnetic pulsations observed at the synchronous altitude are classified into three types, in respect to their major magnetic field polarization in different directions, local time dependence, and different characteristics of accompanied flux modulations of energetic particles, i.e., two nearly transverse waves with the azimuthal and the radial polarization, and the compressional stormtime pulsations. Firstly, we formulate the drift kinetic theory of particle flux modulations under the constraint of the magnetic moment conservation. A generalized energy integral of the particle motion interacting with a ULF-wave with the three-dimensional structure propagating to the azimuthal direction is obtained in the L-shell coordinate of a mirror magnetic field. Its linearized form is reduced to the same form as the previously derived energy change, including the bounce-drift resonant interaction. It is shown that the perturbed guiding center distribution function of energetic particles consists of four contributions, the adiabatic mirror effect corresponding to pitch-angle change, the kinetic effects due to energy change and the accompanying L-shell displacement, and the bounceaveraged drift phase bunching. Secondly, the basic HM-wave modes constitutingcoupling ULF oscillations in non-uniform plasmas are discussed in different models of approach for different plasma states. The diamagnetic drift Alfvén wave and the compressional drift wave with a larger azimuthal mode number in a high-beta plasma are candidates for the stormtimes pulsations. The former is intrinsically a guided localized mode, while the latter is a non-localized mode. By making use of the above preparation, we apply the developed drift kinetic theory to interpret the phase relationships between the ion flux modulation and the geomagnetic pulsation in some selected examples of observations, demonstrating a fair agreement in theoretical results with the observations.     

Observations of Pg pulsations in the Northern Auroral zone and at lower latitude conjugate regions

An analysis is made of giant pulsation (Pg) data recorded at ground stations in the Northern Auroral Zone in Scandanavia (mainly at Tromsø, L = 6.4 and Kiruna, L = 5.5) during the period September 1976 to December 1977. They are shown to have a meridional variation of amplitude and polarization consistent with a field line resonance structure and their vertical component behaviour suggests that they also have a rapid azimuthal phase variation. Limited data from conjugate stations at L = 4.4 are used to show that Pg's are odd mode oscillations of the field line. Pg's are equated to the observation of a unique compressional wave in space at synchronous orbit and it is suggested that they result from the drift wave instability of the compressional Alfven wave at the outer edge of the quiet time ring current.     

Observations of a quiet-time Pc5 wave in the outer magnetosphere

We report an observation of the radial profile of a Pc5 magnetic pulsation and the associated energetic electron flux oscillations from 10 to 18 Re, recorded by the IMP-5 satellite at 19.00 M.L.T. on 21 March 1970. The Pc5 pulsation was mainly compressional and occurred during extremely quiet geomagnetic conditions. Fluxes of energetic electrons detected above three energy thresholds (18, 45, and 80 keV) were found to oscillate out of phase with magnetic field intensity. One new result is that both the wave amplitude and the wave period increased with radial distance. Second, the electron flux oscillation amplitude was roughly proportional to magnetic field fluctuation amplitude and wave period. The wave event is found to be interpreted better as an ion drift wave because of lack of polarization reversal. The characteristics of energetic electron flux oscillations are shown to agree qualitatively with theoretical calculations of the kinetic perturbation of distribution functions by compressional waves.     

Evolution of latitude zonal structure of the large-scale magnetic field in solar cycles

Properties of a latitude zonal component of the large-scale solar magnetic field are analyzed on the basis of H charts for 1905–1982. Poleward migration of prominences is used to determine the time of reversal of the polar magnetic field for 1870–1905. It is shown that in each hemisphere the polar, middle latitude and equatorial zones of the predominant polarity of large-scale magnetic field can be detected by calculating the average latitude of prominence samples referred to one boundary of the large-scale magnetic field. The cases of a single and three-fold polar magnetic field reversal are investigated. It is shown that prominence samples referred to one boundary of the large-scale magnetic field do not have any regular equatorward drift. They manifest a poleward migration with a variable velocity up to 30 m s^-1 depending on the phase of the cycle. The direction of migration is the same for both low-latitude and high-latitude zones. Two different time intervals of poleward migration are found. One lasts from the beginning of the cycle to the time of polar magnetic field reversal and the other lasts from the time of reversal to the time of minimum activity. The velocity of poleward migration of prominences during the first period is from 5 m s^-1 to 30 m s^-1 and the second period is devoid of regular latitude drift.     

The effect of damping on geomagnetic pulsation amplitude and phase at ground observatories

We present some results from a model of forced oscillations of the magnetosphere. The purpose of this work is to examine the effects and consequences of damping on geomagnetic pulsations as observed on the ground. The aim of the current work is to quantify the amount of damping applicable to geomagnetic pulsation waveforms. Ionospheric conductivities vary with latitude and time of day and this variation will effect the damping of geomagnetic pulsations. The variations in ionospheric conductivities are taken into account to predict the changes in amplitude and phase of geomagnetic pulsations over an extended latitudinal array of ground observatories. Three situations are modelled where the damping factor γ/ω_n, which is related to the amplitude loss per cycle, is different: (i) γ/ω_n approximately equal to 0.01, this corresponds to the ionospheric Joule damping of Newton et al. (1978); (ii) λ/ω_n equal to 0.1, this value is consistent with the empirically determined day-time damping factors from the observed latitude-dependent transient decays of the pulsation single effect events discussed by Siebert (1964). The value of 0.1 as the damping factor is taken as typical of day-time conditions and its effect on amplitude and phase for continuous pulsations is considered; and (iii) λ/ω_n is latitude-dependent; three different levels of damping are used appropriate for the night-time conditions associated with the auroral electrojet, plasmatrough and plasmasphere.The results from the model suggest that observationally determined damping factors are greater than those computed from ionospheric Joule damping alone. The model also illustrates the broadening of the latitudinal resonance width with increasing damping and the reducing of the phase change across resonance to less than 180°. The model also successfully reproduces features of pulsation single effect events and Pi2 pulsations.     

Measuring Stellar Differential rotation with asteroseismology

The variation of rotation with latitude is poorly known on stars other than the Sun. Several indirect techniques, photometric and spectroscopic, have been used to search for departure from rigid rotation for sufficiently fast rotators. Here we investigate the possibility of measuring stellar differential rotation for solar-type stars through asteroseismology. Rotationally split frequencies of global oscillation provide information about rotation at different latitudes depending on the azimuthal order, m, of the mode of pulsation. We present a method to estimate differential rotation based on the realization that the m = ±1 and m = ±2 components of quadrupole oscillations can be observed simultaneously in asteroseismology. Rotational frequency splittings can be inverted to provide an estimate of the difference in stellar angular velocity between the equator and 45° latitude. The precision of the method, assessed through Monte Carlo simulations, depends on the value of the mean rotation and on the inclination angle between the rotation axis and the line of sight.     

A comparison between simultaneous I.P.D.P. groundbased observations and observations of energetic protons obtained by satellites

A study of simultaneous groundbased observations of I.P.D.P. (intervals of pulsation of diminishing period) magnetic field fluctuation events and satellite observations of energetic protons have been performed. Some of our results are as follows. (1) The region of I.P.D.P. occurrence is always located equatorward of the isotropic proton precipitation. (2) The I.P.D.P. generation is not connected with the poleward leap of the aurora and the poleward expansion of the precipitating protons. (3) In the evening to afternoon sector enhanced pitch angle scattering is found near L = 4 during I.P.D.P. events, earlier shown to be associated with ion cyclotron resonance. (4) I.P.D.P. events seem to be associated with increased fluxes of (40–60) keV protons injected during substorms near the plasmapause in the equatorial plane.In order to explain the observations we invoke the following model: at substorm onset ring current protons are injected deep into the nightside magnetosphere covering a certain region in L and L.T., with the inner edge of the proton population following McIlwain's injection boundary. The protons drift azimuthally westward and generate ion cyclotron waves in a certain L interval at or inside the plasmapause. By taking into account the shape and position of the plasmapause and the injection boundary, the exterrt and position of the wave generating region can be determined. The frequency-time dispersion of the I.P.D.P. is largely attributed to the L-dependent drift velocity of protons in a narrow energy band. The model is able to explain the observations during several individual events. Also, the model predicts the general trends that have been found by statistical analysis of I.P.D.P. events and accounts for the constant frequency observed by satellites during I.P.D.P. events.     

High latitude,outer zone boundary observations of electrons and protons

The diurnal variation of the high latitude outer zone boundary at 1400km has been determined for electrons ?140keV electrons, and for protons in two energy intervals: 0.56?E?1.1 MeV, 1.1?E?3.2 MeV, from detectors aboard the NOAA-2 satellite. The dependence of the 140 keV electron boundary on D_st has been examined as well?. A wel?l defined correlation of boundary position with D_st is found to exist during the main phase of disturbances, together with an evident local time dependence. All the boundaries are found to be consistent witn the supposition of adiabatic drift and demonstrate the stability of the boundary position over approximately ten years of comparable observation. No statistically significant hemispheric differences in boundary location were observed to occur.     

On the diurnal period variation of midlatitude ULF pulsations

Magnetometer studies of the periods of mid-latitude ULF pulsations have produced conflicting results on the variation of the pulsation periods with both latitude and local time. Since the mid-latitude geomagnetic field is not expected to be significantly distorted by the solar wind, the observed diurnal period variations should be determined by changes in the ambient plasma density. We have applied a physically realistic plasmasphere model to the determination of pulsation eigenperiods over a 24-h interval at L=2.3 (appropriate to Wellington, New Zealand). The resulting model pulsation eigenperiods are largest during the day, with minimum and maximum values at 05.00 and 18.00 L.T. respectively. The model predicts a general increase in the eigenperiods during the replenishment of the protonosphere after a period of geomagnetic activity.     

Motion of charged particles in the magnetosphere

The adiabatic motion of charged particles in the magnetosphere has been investigated using Mead-Fairfield magnetospheric field model (Mead and Fairfield, 1975). Since the motion of charged particles in a dipolar field geometry is well understood, we bring out in this paper some important features in characteristic motion due to non-dipolar distortions in the field geometry. We look at the tilt averaged picture of the field configuration and estimate theoretically the parameters like bounce period, longitudinal invariant and the bounce averaged drift velocities of the charged particle in the Mead-Fairfield field geometry. These parameters are evaluated as a function of pitch angle and azimuthal position in the region of ring current (5 to 7 Earth radii from the centre of the Earth) for four ranges of magnetic activity. At different longitudes the non-dipolar contribution as a percentage of dipole value in bounce period and longitudinal invariant show maximum variation for particles close to 90° pitch angles. For any low pitch angle, these effects maximize at the midnight meridian. The radial component of the bounce averaged drift velocity is found to be greatest at the dawn-dusk meridians and the contribution vanishes at the day and midnight meridians for all pitch angles. In the absence of tilt-dependent terms in the model, the latitudinal component of the drift velocity vanishes. On the other hand, the relative non-dipolar contribution to bounce averaged azimuthal drift velocity is very high as compared to similar contribution in other characteristic parameters of particle motion. It is also shown that non-dipolar contribution in bounce period, longitudinal invariant and bounce averaged drift velocities increases in magnitude with increase in distance and magnetic activity.     

Modulation Cycles of Galactic Cosmic Ray Diurnal Anisotropy Variation

The diurnal variation of the galactic cosmic ray (GCR) count rates measured by a ground-based neutron monitor (NM) station represents an anisotropic flow of GCR at 1 AU. The variation of the local time of GCR maximum intensity (we call the phase) is thought in general to have a period of two sunspot cycles (22 years). However, other interpretations are also possible. In order to determine the cyclic behavior of GCR anisotropic variation more precisely, we have carried out a statistical study on the diurnal variation of the phase. We examined 54-year data of Huancayo (Haleakala), 40-year data from Rome, and 43-year data from Oulu NM stations using the ‘pile-up’ method and the F-test. We found that the phase variation has two components: of 22-year and 11-year cycles. All NM stations show mainly the 22-year phase variation controlled by the drift effect due to solar polar magnetic field reversal, regardless of their latitudinal location (cut-off rigidity). However, the lower the NM station latitude is (the higher the cut-off rigidity is), the higher is the contribution from the 11-year phase variation controlled by the diffusion effect due to the change in strength of the interplanetary magnetic fields associated with the sunspot cycle.     

Meridional drift in the large-scale solar magnetic field pattern

A method of two-dimensional correlation functions has been applied to a sequence of synoptic maps of the large-scale magnetic field to obtain the meridional drift pattern of field structures. The meridional drift profile obtained is antisymmetric about the equator. The meridional drift is directed from the equator to the poles at latitudes below 45°. A maximum drift velocity of 11–13 m s^–1 is attained in the latitude range 30°. A picture of the space-time distribution of meridional drift is also obtained, which may be interpreted as resulting from the effect of azimuthal convective rolls (3 rolls per hemisphere) on the large-scale magnetic field. Rolls originate at high latitudes following the cycle maximum, and migrate equatorwards until the minimum of the next cycle. The picture in the equatorial region can correspond to convective rolls with lifetimes of about two years, or to the process of interaction of rolls from two hemispheres.     

VLF goniometer observations at halley bay,Antarctica—II. Magnetospheric structure deduced from whistler observations

On 26 July 1967, a magnetically quiet day (ΣK_p = 12?) with high whistler activity at Halley Bay, it was found possible, by measurement of whistler nose-frequency and dispersion and the bearings of the whistler exit points, to make a detailed study of the magnetospheric structure associated with the whistler ducts.During the period 0509–2305 UT most of the exit points of whistlers inside the plasmasphere were situated along a strip about 100km wide passing through Halley Bay in an azimuthal direction 30°E of N between 57° and 62° invariant latitude. A mechanism which can give rise to such a well-defined locus which co-rotates with the Earth is not clear. Nevertheless, it does appear that the locus coincides with the contour of solar zenith angle 102° at 1800 UT 25 July. This was also the time of occurrence of a sub-storm and it is suggested that the magnetospheric structure was initiated by proton precipitation along the solar zenith angle 102° contour.At mid-day knee-whistlers observed outside the plasmapause had exit points which were closely aligned along an L-shell at an invariant latitude of 62.5°. They exhibited a marked variation (~ 3:1) in electron tube content over about 12° of invariant longitude and a drift of about 8 msec^?1 to lower L-shells.Throughout the period of observation the plasmapause lay about 2° polewards of the mean position found by Carpenter (1968) for moderately disturbed days.     

Measurement of Kodaikanal white-light images – IV. Axial Tilt Angles of Sunspot Groups

The Kodaikanal sunspot data set, covering the interval 1906–1987, is used in conjunction with the similar Mount Wilson sunspot data set, covering the interval 1917–1985, to examine characteristics of sunspot group axial tilt angles. Good agreement is demonstrated between various results derived from the two independent data sets. In particular, the tendency for sunspot groups near the average tilt angle to be larger than those far from the average tilt angle is confirmed. Similarly the faster residual rotation rate for groups near the average tilt angle is also confirmed. Other confirmations are made for the relationships between latitude drift of sunspot groups and tilt angle, polarity separations, and axial expansion. Evidence is presented that tilt angles averaged over these long time intervals differ between the north and south hemispheres by about 1.4 deg. It is suggested that residual tilt angles show a slight systematic variation with phase in the activity cycle.     

Ionospheric Joule dissipation as a damping mechanism for high latitude ULF pulsations: Observational evidence

On 9 January 1979 an SI-excited pulsation event was observed by the Scandinavian Magnetometer Array. The pulsation period shows a clear variation with latitude which suggests decoupled oscillations of individual magnetic field shells. The pulsation amplitudes exhibit an e-fold decay with the damping rate γ varying both in longitudinal and latitudinal directions. Assuming Joule heating in the ionosphere as the dominant damping mechanism (and thus γ ～ Σ^?1_p) approximate height-integrated Pedersen conductivity profiles were calculated which fit well with previously observed Σ_p distributions. This is interpreted as observational evidence for ionosopheric Joule dissipation as the major damping mechanism for high-latitude ULF-pulsations.     

The horizontal propagation of Pc1 pulsations in the ionosphere

The propagation properties of the various elements of the plane-wave angular spectrum of a Pc1 pulsation signal in the ionosphere are determined by a full-wave numerical analysis. A spectral component is characterized by the wave-vector azimuthal direction, and the Snell constant S. The isotropic R-mode transmission coefficient to ground is fairly flat for S ? 400, but thereafter (S > 500) drops rapidly with increasing S. Coupling of energy from the field-guided L-mode to the R-mode occurs along the entire length of the L-mode trajectory within the ionospheric duct in which the R-mode can propagate. Within the duct, the R-mode attenuation is determined largely by R to L-mode coupling, which is larger for E-W than for N-S azimuths, especially for steep angles of incidence (S < 100). This should lead to enhanced injection of energy into E-W high altitude, high velocity paths, but to higher E-W attenuation at oblique angles. For oblique propagation (S ? 200) horizontal group velocities are slightly higher than the Alfvén phase velocity at the F-layer peak, but about twice as high for steep angles (S ≈ 100).     

A Comparative Study of Jovian Anticyclone Properties from a Six-Year (1994-2000) Survey

We have used the Hubble Space Telescope archived images of Jupiter for the period 1994-2000, complemented by ground-based telescopic observations, to study in detail the long-term properties of synoptic-scale anticyclonic vortices (size > 1500 km, lifetime > months). We have also analyzed a set of Voyager 1 and 2 images obtained in 1979 to compare anticyclone properties from the two different periods. The latitudinal range covered by this survey spans 75°N to 75°S, encompassing 12 anticyclonic zones. We present data on vortex size, aspect ratio, number, latitude location, lifetime, motion, interaction, and morphology for more than 100 vortices. We study empirically the relation between these properties and the mean ambient zonal flow.We show that most of these properties are not related to latitude and location relative to the jet pattern. However, a significant linear anticorrelation is found when plotting vortex relative speed (vortex propagation speed minus zonal flow velocity) against the zonal flow velocity at the central latitude of each vortex. As the mean eastward flow increases its velocity within a given anticyclonic domain, vortex velocity becomes more westward. This relation holds for all anticyclones except those moving at high velocities (at latitudes 20°S and 23°N), whose origin appears to be of a different nature. Moreover, there is also some evidence that the drift rate could be connected to the planetary minus flow vorticity gradient (most conspicuous at 40 and 45°N). We present simple dynamical interpretations of these observations.     

A theoretical evaluation of the lunar tidal variations in the ionospheric F2-layer

Time-varying solutions of the full continuity equation for electrons in the F2-region are obtained. The effects of production, loss, diffusion and electrodynamic ‘E × B’ drift are taken into account. The ‘E × B’ drift term consists of a solar and a lunar component. The solar component of drift is assumed diurnal with 14.6m/sec maximum upward speed at mid-day. The lunar component is assumed sinusoidal with period of half lunar day and amplitude one tenth of the solar drift; the phase is assumed to remain constant in lunar time, in accordance with Chapman's phase law.The results show that the lunar variations in the F2-region are markedly dependent on solar time and latitude. It is also shown that the average semi-diurnal lunar variations in N_mF2 and h_mF2 at any particular lunar time are almost opposite in phase to each other (i.e. out of phase by 6 hr) in the magnetic equatorial zone, and out of phase by 2 hr at moderate latitudes. The phase of δh_mF2 is 10 hr at low latitudes and 9 hr at moderate latitudes. The phase of δN_mF2 is 4 hr at low latitudes and 11 lunar hr at moderate latitudes.The results also show that the phase of the lunar semi-monthly oscillations in N_mF2 undergoes a rapid shift of about 5 lunar hr in going from 8 to 12° and the so called phase reversal occurs at about 10° lat at which the amplitude of N_mF2. becomes extremely small.These and other results are in good agreement with observations. Thus it is shown that the main features of the observed lunar tidal variations of the F2-region within 20° of the magnetic equator can be explained satisfactorily by the superposition of a small lunar drift on a large solar drift.     

14 July 2000, a near-global coronal event and its association with energetic electron events detected in the interplanetary medium

On 14 July 2000, the LASCO coronagraphs showed a very fast halo coronal mass ejection in association with the radio bursts seen shortly after 10:00 UT. Radio imaging observations by the Nançay radioheliograph (NRH) of these bursts showed a very complex event that can be regarded as global: the sources encompassed all the visible range in longitude and a huge span in latitude. Another interesting feature of the radio event is its recurrent nature: after the most intense phase shortly after 10:00 UT, two other strong outbursts are detected, one at about 12:50 UT and another at about 13:48 UT. All of these sub-events showed similar development and likely evidence for CMEs. The launch of a CME in association with the 14:00 UT sub-event is inferred from WIND/WAVES, with interplanetary type II signatures in the hectometric wavelength range at that time. These later events were not detected by LASCO due to energetic particles hitting the CCD. During the Bastille Day event, energetic particle observations measured in situ by ACE/EPAM are dominated by energetic electrons. Changes in anisotropy and energy spectrum of the ～38–350 keV electrons suggest a good correlation with the coronal radio observations. In addition to the three main radio events and particle observations, the NRH data reveal moving features in the southern hemisphere. These moving features, located at about 45 deg south and with an angular extent of about 45 deg, are illuminated by non-thermal electrons and are seen at distances up to 2.5 solar radii from the Sun center. More generally, we interpret the global and recurrent coronal activity, revealed by the radio data, as responsible for populating the interplanetary medium with energetic electrons.     

Plasmatrough ion mass densities determined from ULF pulsation eigenperiods

Auroral radar studies of ULF pulsations have proved useful in determining the spatial characteristics of resonant oscillations. A particular class of ringing or transient pulsations has been identified in the radar data as toroidal mode eigenoscillations. We have considered a total of 64 events of this type recorded by either the STARE radar in Scandinavia, or the Slope Point radar in New Zealand, giving a combined latitudinal coverage of approx. 12°. These events are interpreted as toroidal mode eigenoscillations; the periods for individual events and the mean periods increase with geomagnetic latitude. Use of hydromagnetic resonance theory allows the equatorial ion mass density to be determined. The densities obtained are appropriate to the plasmatrough and range from ～ 10 to 100 a.m.u. cm^?3 near geosynchronous orbit. The radial variation in the equatorial plane is typically R^?5 in the midnight-noon sector and R^?3 in the noon-midnight sector. To reconcile these pulsation periods with in situ electron density measurements implies that H⁺ ion densities in the range ～ 1–10 cm^?3 and ～50% O⁺ ions are required.