Pc5 Waves Research Articles

Formation of Butterfly Pitch Angle Distributions of Relativistic Electrons in the Outer Radiation Belt With a Monochromatic Pc5 Wave

AbstractRadial transport of relativistic electrons in the inner magnetosphere has been considered one of the acceleration mechanisms of the outer radiation belt electrons and can be driven by the drift resonance with Pc5 ultralow‐frequency waves. In this study, we focus on the pitch angle dependences of the radial transport and investigate the pitch angle distributions (PADs) of the relativistic electrons after interaction with a monochromatic Pc5 wave, using two simulation models of the inner magnetosphere: Geomagnetic Environment Modeling System for Integrated Studies (GEMSIS)‐Ring Current and GEMSIS‐Radiation Belt models. The results show that butterfly‐like PADs (two peaks in small pitch angles) are produced at a fixed radial distance (L‐shell) and energy, where the electrons with small pitch angles satisfy the resonance condition. The radial range and timescale of the butterfly PAD appearance correspond to the maximum radial transport of the electron (i.e., resonant width) and its oscillation period in the radial direction, respectively. We analytically derive the resonant width and oscillation period of relativistic electrons interacting with a monochromatic Pc5 wave, assuming conservation of first and second adiabatic invariants, and compare them with the simulation result. The results show that electrons with small equatorial pitch angles at a fixed magnetic moment can interact with a monochromatic Pc5 wave in a wider L‐shell range than perpendicular electrons due to the pitch angle changes, in the course of the radial transport.

A Distributed Lag Autoregressive Model of Geostationary Relativistic Electron Fluxes: Comparing the Influences of Waves, Seed and Source Electrons, and Solar Wind Inputs

AbstractRelativistic electron flux at geosynchronous orbit depends on enhancement and loss processes driven by ultralow frequency (ULF) Pc5, chorus, and electromagnetic ion cyclotron (EMIC) waves, seed electron flux, magnetosphere compression, the “Dst effect,” and substorms, while solar wind inputs such as velocity, number density, and interplanetary magnetic field Bz drive these factors and thus correlate with flux. Distributed lag regression models show the time delay of highest influence of these factors on log10 high‐energy electron flux (0.7–7.8 MeV, Los Alamos National Laboratory satellites). Multiple regression with an autoregressive term (flux persistence) allows direct comparison of the magnitude of each effect while controlling other correlated parameters. Flux enhancements due to ULF Pc5 and chorus waves are of equal importance. The direct effect of substorms on high‐energy electron flux is strong, possibly due to injection of high‐energy electrons by the substorms themselves. Loss due to electromagnetic ion cyclotron waves is less influential. Southward Bz shows only moderate influence when correlated processes are accounted for. Adding covariate compression effects (pressure and interplanetary magnetic field magnitude) allows wave‐driven enhancements to be more clearly seen. Seed electrons (270 keV) are most influential at lower relativistic energies, showing that such a population must be available for acceleration. However, they are not accelerated directly to the highest energies. Source electrons (31.7 keV) show no direct influence when other factors are controlled. Their action appears to be indirect via the chorus waves they generate. Determination of specific effects of each parameter when studied in combination will be more helpful in furthering modeling work than studying them individually.

Transmission of a Magnetospheric Pc1 Wave Beam Through the Ionosphere to the Ground

AbstractA characteristic feature of the upper ionosphere is the occurrence of the ionospheric Alfvén resonator (IAR) and the ionospheric fast‐mode waveguide (IFW), which can trap electromagnetic wave energy in the Pc1 frequency range from fractions of a hertz to a few hertz. This wave trapping ensures the dependence of ionospheric transmission/reflective properties on frequency. We numerically model magnetospheric Alfvén wave transmission through the ionosphere to the ground based on solution of the magnetohydrodynamic full‐wave equations in a realistic ionosphere, whose parameters are reconstructed from the International Reference Ionosphere model. The spatial structure of an incident wave is modeled as a localized beam with a finite latitudinal scale λ⊥ and an azimuthally propagating wave. The IAR and IFW modes are coupled owing to the Hall conductivity and geomagnetic field inclination. Ground spatial and spectral structures of the Pc1 wave (0.1‐ to 6‐Hz band) have been calculated for summer day/night conditions at the Antarctic Halley observatory, though results may be qualitatively applied to any midlatitude site. An incident wave with azimuthal wave vector ky=1.7·10−3 km−1, and λ⊥=102 km has been considered. The model predicts that beneath the incident beam the ground magnetic response “duplicates” its structure after accounting for a π/2 rotation and some latitudinal shift/widening. The transmission has an oscillatory dependence on frequency, thus forming “transmission windows” at resonant frequencies of IAR and IFW modes. The spectra vary depending on distance from an incidence point. Interference between IFW modes is revealed in the nonmonotonic and frequency‐dependent character of latitudinal variations of wave amplitude.

Survey of ionospheric Pc3-5 ULF wave signatures in SuperDARN high time resolution data.

Ionospheric signatures of ultra-low frequency (ULF) wave in the Pc3-5 band (1.7-40.0 mHz) were surveyed using ~6 s resolution data from Super Dual Auroral Radar Network (SuperDARN) radars in the northern hemisphere from 2010 to 2016. Numerical experiments were conducted to derive wave period dependent thresholds for automated detection of ULF waves using the Lomb-Scargle periodogram technique. The spatial occurrence distribution, frequency characteristics, seasonal effects, solar wind condition and geomagnetic activity level dependence have been studied. Pc5 wave events were found to dominate at high and polar latitudes with a most probable frequency of 2.08 ± 0.07 mHz while Pc3-4 waves were relatively more common at midlatitudes on the nightside with a most probable frequency of 11.39 ± 0.14 mHz. At high latitudes, the occurrence rate of Pc4-5 waves maximizes in the dusk sector and during winter. These events tend to occur during low geomagnetic activity and northward interplanetary magnetic field (IMF). For the category of radially bounded but longitudinally extended Pc4 events in the duskside ionosphere, an internal driving source is suggested. At midlatitudes, the Pc3-4 occurrence rate maximizes premidnight and during equinox. This tendency becomes more prominent with increasing auroral electrojet (AE) index and during southward IMF, which suggests many of these events are Pi2 and Pc3-4 pulsations associated with magnetotail dynamics during active geomagnetic intervals. The overall occurrence rate of Pc3-5 wave events is lowest in summer, which suggests that the ionospheric conductivity plays a role in controlling ULF wave occurrence.

An initial ULF wave index derived from 2 years of Swarm observations

Abstract. The ongoing Swarm satellite mission provides an opportunity for better knowledge of the near-Earth electromagnetic environment. Herein, we use a new methodological approach for the detection and classification of ultra low-frequency (ULF) wave events observed by Swarm based on an existing time-frequency analysis (TFA) tool and utilizing a state-of-the-art high-resolution magnetic field model and Swarm Level 2 products (i.e., field-aligned currents – FACs – and the Ionospheric Bubble Index – IBI). We present maps of the dependence of ULF wave power with magnetic latitude and magnetic local time (MLT) as well as geographic latitude and longitude from the three satellites at their different locations in low-Earth orbit (LEO) for a period spanning 2 years after the constellation's final configuration. We show that the inclusion of the Swarm single-spacecraft FAC product in our analysis eliminates all the wave activity at high altitudes, which is physically unrealistic. Moreover, we derive a Swarm orbit-by-orbit Pc3 wave (20–100 MHz) index for the topside ionosphere and compare its values with the corresponding variations of solar wind variables and geomagnetic activity indices. This is the first attempt, to our knowledge, to derive a ULF wave index from LEO satellite data. The technique can be potentially used to define a new Level 2 product from the mission, the Swarm ULF wave index, which would be suitable for space weather applications. Keywords. Space plasma physics (waves and instabilities)

Discovery of 1 Hz Range Modulation of Isolated Proton Aurora at Subauroral Latitudes

AbstractIsolated proton aurora (IPA) is a manifestation of the wave‐particle interaction visible at subauroral latitudes, with activity on many timescales. We herein present the first observational evidence of rapid luminous modulation of IPA correlated with simultaneously observed Pc1 waves observed on the ground, which are equivalent to the electromagnetic ion cyclotron (EMIC) waves in the magnetosphere. The fastest luminous modulation of IPA was observed in the 1 Hz frequency range, which was twice the frequency of the related Pc1 waves. The time lag between variations of Pc1 wave power and the IPA luminosity suggests that the source regions of IPA are distributed near the magnetic equator, suggesting an EMIC wave‐energetic (a few tens of keV) proton or relativistic (MeV or sub‐MeV) electron interaction. The generation mechanism of this 1 Hz luminous modulation remains an open issue, but this study supports the importance of nonlinear pitch angle scattering via wave‐particle interactions.

Applying the new method of time-frequency transforms to the analysis of the characteristics of geomagnetic Рс5 pulsations

This study considers the possibility of using the new methods of time-frequency transforms, such as chirplet and warblet transforms, to analyze the digital observational data of geomagnetic pulsations of Pc5 type. For this purpose, necessary algorithms of calculation and appropriate software were developed. The chirplet transform method (CT) is used to analyze signals with a linear frequency modulation. A chirplet variation, the so-called warblet transform, is used to analyze signals with a nonlinear frequency modulation. Since, in studying geomagnetic pulsations, it is difficult to make assumptions on the character of the behavior of the instantaneous frequency of the signal, the special generalized warblet transform (GWT) was used for the analysis. The GWT has a high spatiotemporal resolution and was developed to analyze oscillations both with a periodic and nonperiodic change of the instantaneous frequency. The software developed for GWT calculation was used to study daytime geomagnetic Pc5 pulsations with durations of several hours that were detected via the network of ground-based magnetometers of the Scandinavian IMAGE profile during the magnetic storm of May 29–30, 2003. For the first time, temporal variations of the instantaneous frequency of geomagnetic pulsations are determined and their possible use in studying the fine spatial structure of Pc5 waves is shown.

Relativistic electron dynamics produced by azimuthally localized poloidal mode ULF waves: Boomerang‐shaped pitch angle evolutions

AbstractWe present an analysis of “boomerang‐shaped” pitch angle evolutions of outer radiation belt relativistic electrons observed by the Van Allen Probes after the passage of an interplanetary shock on 7 June 2014. The flux at different pitch angles is modulated by Pc5 waves, with equatorially mirroring electrons reaching the satellite first. For 90° pitch angle electrons, the phase change of the flux modulations across energy exceeds 180° and increasingly tilts with time. Using estimates of the arrival time of particles of different pitch angles at the spacecraft location, a scenario is investigated in which shock‐induced ULF waves interact with electrons through the drift resonance mechanism in a localized region westward of the spacecraft. Numerical calculations on particle energy gain with the modified ULF wavefield reproduce the observed boomerang stripes and modulations in the electron energy spectrogram. The study of boomerang stripes and their relationship to drift resonance taking place at a location different from the observation point adds new understanding of the processes controlling the dynamics of the outer radiation belt.

Nighttime Pc3 pulsations: MM100 and MAGDAS observations

In this paper, we present a statistical and case analysis of nighttime Pc3 pulsations observed from middle to equatorial latitudes during the year 2003. We found two groups of nighttime Pc3 pulsations. Pc3s of the first group are in fact the nightside counterpart of morning Pc3 pulsations with large azimuthal scales slowly attenuating toward midnight. Such night signatures of morning Pc3 waves are observed during the periods of fast solar wind (V>500,hbox {km/s}). The second type is the locally generated night Pc3 pulsations. They can be observed under moderate solar wind velocities. Maximal occurrence rates and amplitudes for these pulsations are recorded at middle geomagnetic latitudes near the local magnetic midnight. Probably, they are associated with auroral activations or local non-substorm bursty processes.Graphical abstractAn example of a nighttime Pc3 pulsationGraphical abstractEmpirical probability density function of the solar wind speed for the intervals with the two types of night Pc3 pulsations and for all the intervals analyzed.

Fast‐moving diffuse auroral patches: A new aspect of daytime Pc3 auroral pulsations

AbstractAuroral pulsations are a convenient diagnostic of wave‐particle interactions in the magnetosphere. A case study of a daytime Pc3 (22–100 mHz) auroral pulsation event, measured with a ~2 Hz sampling all‐sky camera at South Pole Station (74.4°S magnetic latitude) on 17 May 2012, is presented. The daytime Pc3 auroral pulsations were most active in a closed field line region where the aurora was dominated by diffuse green‐line emissions and within ±2 h of magnetic local noon. Usually, but not always, the corresponding periodic variations were recorded with a colocated search coil magnetometer. Of particular interest is the two‐dimensional auroral signature, indicating that the temporal luminosity variations at a given point were due to repeated formation and horizontal motion of faint, nonpulsating auroral patches with scale sizes of ~100 km. The individual patches propagated equatorward with speeds of 15 km s−1 up to 20–25 km s−1 one after another along the magnetic meridian through local magnetic zenith. These properties differ considerably from typical pulsating aurorae, being periodic on‐off luminosity variations in a particular auroral patch and drifting in accordance with the convection electric field in the magnetosphere. We speculate that such repetitive patterns of the fast‐moving auroral patches, being another aspect of the daytime Pc3 auroral pulsations, may be a visible manifestation of compressional Pc3 waves which propagate earthward and cause modulation of precipitating keV electron fluxes in the dayside outer magnetosphere.

In situ statistical observations of Pc1 pearl pulsations and unstructured EMIC waves by the Van Allen Probes

AbstractWe present here the first in situ statistical survey of structured Pc1 pearl pulsations compared with unstructured electromagnetic ion cyclotron (EMIC) waves observed by the Van Allen Probes spacecraft. This data set was compiled from observations spanning 8 September 2012 through 31 August 2015 and comprises over 1630 h of total EMIC wave activity, of which 291 h exhibited pearl structure. Additionally, we have identified 29 wave events demonstrating periodically oscillating wave packets, mostly about the magnetic equator, indicated by the reversal of Poynting flux along the background magnetic field. We have found several stark differences between Pc1 pearl pulsations and unstructured EMIC waves. While unstructured EMIC waves demonstrate the predicted behavior of a higher occurrence across the dayside with enhanced wave power at dusk, pearl pulsations occur uniformly across magnetic local time, with a small enhancement in the late morning sector. Pearl pulsations were more often observed during magnetospherically quiet periods, particularly in the late recovery period of geomagnetic storms. The mean excitation frequency of pearl pulsations was observed to be independent of the local ion cyclotron frequency, and individual wave investigations indicate that the modulation period also remained constant for the duration of the event over a finite range in L. We examine three possible generation mechanisms—the bouncing wave packet model, modulation by ultralow‐frequency Pc4 and Pc5 waves, and the formation of an ion cyclotron resonator—but are unable to definitively confirm the validity of any one model.

On the predictive potential of Pc5 ULF waves to forecast relativistic electrons based on their relationships over two solar cycles

AbstractA statistical study of relativistic electron (>2 MeV) fluence derived from geosynchronous satellites and Pc5 ultralow frequency (ULF) wave power computed from a ground magnetic observatory data located in Canada's auroral zone has been carried out. The ground observations were made near the foot points of field lines passing through the GOESs from 1987 to 2009 (cycles 22 and 23). We determine statistical relationships between the two quantities for different phases of a solar cycle and validate these relationships in two different cycles. There is a positive linear relationship between log fluence and log Pc5 power for all solar phases; however, the power law indices vary for different phases of the cycle. High index values existed during the descending phase. The Pearson's cross correlation between electron fluence and Pc5 power indicates fluence enhancement 2–3 days after strong Pc5 wave activity for all solar phases. The lag between the two quantities is shorter for extremely high fluence (due to high Pc5 power), which tends to occur during the declining phases of both cycles. Most occurrences of extremely low fluence were observed during the extended solar minimum of cycle 23. The precursory attribute of Pc5 power with respect to fluence and the enhancement of fluence due to rising Pc5 power both support the notion of an electron acceleration mechanism by Pc5 ULF waves. This precursor behavior establishes the potential of using Pc5 power to predict relativistic electron fluence.

Effects of ULF waves on local and global energetic particles: Particle energy and species dependences

AbstractAfter 06:13 UT on 24 August 2005, an interplanetary shock triggers large‐amplitude ultralow‐frequency (ULF) waves (|δB| ≥ 15 nT) in the Pc4–Pc5 wave band (1.6–9 mHz) near the noon geosynchronous orbit (6.6 RE). The local and global effects of ULF waves on energetic particles are observed by five Los Alamos National Laboratory satellites at different magnetic local times. The large‐amplitude ULF waves cause the synchronous oscillations of energetic electrons and protons (≥75 keV) at the noon geosynchronous orbit. When the energetic particles have a negative phase space density radial gradient, they undergo rapid outward radial diffusion and loss in the wave activity region. In the particle drift paths without strong ULF waves, only the rapidly drifting energetic electrons (≥225 keV) display energy‐dispersive oscillations and flux decays, whereas the slowly drifting electrons (<225 keV) and protons (75–400 keV) have no ULF oscillation and loss feature. When the dayside magnetopause is compressed to the geosynchronous orbit, most of energetic electrons and protons are rapidly lost because of open drift trajectories. The global and multicomposition particle measurements demonstrate that the effect of ULF waves on nonlocal particle flux depends on the particle energy and species, whereas magnetopause shadowing effect is independent of the energetic particle species. For the rapidly drifting outer radiation belt particles (≥225 keV), nonlocal particle loss/acceleration processes could also change their fluxes in the entire drift trajectory in the absence of “Dst effect” and substorm injection.

Solar wind‐driven Pc5 waves observed at a polar cap station and in the near cusp ionosphere

AbstractWe present the results of a comparative study conducted in Antarctica by using the ULF geomagnetic field measurements at Terra Nova Bay (Altitude Adjusted Corrected Geomagnetic Coordinates latitude 80°S) and simultaneous data from the Super Dual Auroral Radar Network radar at South Pole Station. Pc5 waves observed at Terra Nova Bay around local magnetic noon, when the station is close to the dayside cusp, can be interpreted as spatial integrated signals, produced by ionospheric currents associated to field line resonances at somewhat lower latitudes. The radar, providing the Doppler velocities of ionospheric plasma over a range of geomagnetic latitudes, allows to detect the occurrence of such field line resonances. In the reported case, our analysis shows evidence of resonant signals in the ionosphere at 75°S and 76°S that find correspondence in frequency and time with the geomagnetic signals observed at Terra Nova Bay around local noon. During the period of interest, oscillations of the solar wind dynamic pressure at the same frequency are detected by Geotail, just upstream of the morning flank of the bow shock. All the observations are consistent with the interpretation of the signals at Terra Nova Bay in terms of signatures of field line resonances occurring at lower latitudes, driven by solar wind oscillations transmitted into the magnetosphere. We discuss also the possibility of an additional contribution to the signals at Terra Nova Bay, due to the direct propagation of the solar wind waves along the local open field line.

Latitudinal amplitude-phase structure of MHD waves: STARE radar and image magnetometer observations and modeling

We have developed a numerical model that yields a steady-state distribution of field components of MHD wave in an inhomogeneous plasma box simulating the realistic magnetosphere. The problem of adequate boundary condition at the ionosphere–magnetosphere interface for coupled MHD mode is considered. To justify the model’s assumptions, we have derived the explicit inequality showing when the ionospheric inductive Hall effect can be neglected upon the consideration of Alfven wave reflection from the ionospheric boundaries. The model predicts a feature of the ULF spatial amplitude/phase distribution that has not been noticed by the field line resonance theory: the existence of a region with opposite phase delays on the source side of the resonance. This theoretical prediction is supported by the amplitude-phase latitudinal structures of Pc5 waves observed by STARE radar and IMAGE magnetometers. A gradual decrease in azimuthal wave number m at smaller L-shells was observed at longitudinally separated radar beams.

On the origin of the dawn‐dusk asymmetry of toroidal Pc5 waves

AbstractObservations in space and on the ground have shown that azimuthally polarized (toroidal) magnetospheric Pc5 waves (frequency = 1.7–6.7 mHz) are more prominent on the dawnside than on the duskside, but the origin of this dawn‐dusk asymmetry is unknown. To determine the origin, we conducted a statistical analysis of toroidal Pc5 waves at the Time History of Events and Macroscale Interactions during Substorms (THEMIS) D spacecraft using ion bulk velocity measured during 2008–2014. We found that both the detection rate of narrowband Pc5 waves, which are attributed to fundamental toroidal‐mode standing Alfvén waves, and the amplitude of velocity oscillations in the fixed Pc5 band exhibit a dawn‐dusk asymmetry consistent with previous studies regardless of whether the interplanetary magnetic field is directed along the Parker spiral or orthogonal to it. These results imply that the asymmetry originates within the magnetosphere, not from a dawn‐dusk asymmetry of the Kelvin‐Helmholtz instability on the magnetopause. Based on theoretical studies published previously and a global magnetohydrodynamic simulation, we propose that the asymmetry originates from the local time dependence of the radial mass density variation and the associated gradient of the frequency of standing Alfvén waves. The THEMIS data indicate that in the L range 8–12 the average gradient is smaller in the dusk sector than in the dawn sector. This will result in weaker field line resonance (FLR), and thus weaker toroidal Pc5 waves, in the dusk sector. Density structures associated with drainage plumes may also suppress FLR in the dusk sector.

Широтная амплитудно-фазовая структура МГД-волн: моделирование наблюдения на радаре STARE и магнитометрах сети IMAGE

We have developed a numerical model that yields a steady-state distribution of field components of MHD wave in an inhomogeneous plasma box simulating the realistic magnetosphere. The problem of adequate boundary condition at the ionosphere–magnetosphere interface for coupled MHD mode is considered. To justify the model’s assumptions, we have derived the explicit inequality showing when the ionospheric inductive Hall effect can be neglected upon the consideration of Alfven wave reflection from the ionospheric boundaries. The model predicts a feature of the ULF spatial amplitude/phase distribution that has not been noticed by the field line resonance theory: the existence of a region with opposite phase delays on the source side of the resonance. This theoretical prediction is supported by the amplitude-phase latitudinal structures of Pc5 waves observed by STARE radar and IMAGE magnetometers. A gradual decrease in azimuthal wave number m at smaller L-shells was observed at longitudinally separated radar beams.

Correspondence between the ULF wave power spatial distribution and auroral oval boundaries

The world-wide spatial distribution of the wave power in the Pc5 band during magnetic storms has been compared with auroral oval boundaries. The poleward and equatorward auroral oval boundaries are estimated using either the British Antarctic Survey database containing IMAGE satellite UV observations of the aurora or the OVATION model based on the DMSP particle data. The “epicenter” of the spectral power of broadband Pc5 fluctuations during the storm growth phase is mapped inside the auroral oval. During the storm recovery phase, the spectral power of narrowband Pc5 waves, both in the dawn and dusk sectors, is mapped inside the auroral oval or around its equatorward boundary. This observational result confirms previously reported effects: the spatial/temporal variations of the Pc5 wave power in the morning/pre-noon sector are closely related to the dynamics of the auroral electrojet and magnetospheric field-aligned currents. At the same time, narrowband Pc5 waves demonstrate typical resonant features in the amplitude-phase latitudinal structure. Thus, the location of the auroral oval or its equatorward boundary is the preferred latitude for magnetospheric field-line Alfven resonator excitation. This effect is not taken into account by modern theories of ULF Pc5 waves, but it could be significant for the development of more adequate models.

Соответствие между пространственным распределением мощности УНЧ-волн и границами аврорального овала

The location of the auroral oval boundaries has been mapped onto the Pc5 wave power spatial distribution. The poleward and equatorward boundaries of auroral oval are estimated using either the BAS database based on UV observations of the aurora by the IMAGE satellite or the OVATION model based on the DMSP particle data. The epicenter of the spectral power of broadband fluctuations in the Pc5 band during storm growth phase is mapped inside the auroral oval. During the storm recovery phase, the spectral power of narrowband Pc5 waves, both in the morning and dusk sector, is mapped inside the auroral oval or around its equatorward boundary. This observational result confirms the effects earlier reported: the spatial/temporal variations of the Pc5 wave power in the morning/pre-noon sector are closely related to the location of the auroral electrojet and magnetospheric field-aligned currents. At the same time, narrowband Pc5 waves demonstrate typical resonant features in the amplitude-phase latitudinal structure. Thus, the location of the auroral oval (or its equatorward border) is the preferred latitude of magnetospheric field-line Alfven resonator excitation. This effect is not taken into account by modern theories of ULF Pc5 waves, but it could be significant for development of more adequate models.

Modulation of the ionosphere by Pc5 waves observed simultaneously by GPS/TEC and EISCAT

Earlier studies demonstrated that the monitoring of the ionospheric total electron content (TEC) by global satellite navigation systems is a powerful method to study the propagation of transient disturbances in the ionosphere, induced by internal gravity waves. This technique has turned out to be sensitive enough to detect ionospheric signatures of magnetohydrodynamic waves as well. However, the effect of TEC modulation by ULF waves is not well examined as a responsible mechanism has not been firmly identified. During periods with intense Pc5 waves distinct pulsations with the same periodicity were found in the TEC data from high-latitude GPS receivers in Scandinavia. We analyze jointly responses in TEC variations and EISCAT ionospheric parameters to global Pc5 pulsations during the recovery phase of the strong magnetic storms on October 31, 2003. Comparison of periodic fluctuations of the electron density at different altitudes from EISCAT data shows that main contribution into TEC pulsations is provided by the lower ionosphere, up to ~150 km, that is the E-layer and lower F-layer. This observational fact favors the TEC modulation mechanism by field-aligned plasma transport induced by Alfven wave. Analytical estimates and numerical modeling support the effectiveness of this mechanism. Though the proposed hypothesis is basically consistent with the analyzed event, the correspondence between magnetic and ionospheric oscillations is not always perfect, so further studies need to be conducted to understand fully the TEC modulations associated with Pc5 pulsations.