The propagation of compressional Pi 2 waves in the inner magnetosphere is investigated by analyzing the onset delay times between the ground and the geosynchronous altitude. We use the compressional component (northward) of magnetic data from low-latitude stations and the geosynchronous satellite ETS-VIII (GMLat. = -10.8°, GMLon. = 217.5°). The onset delays are determined by a cross-correlation analysis, and we analyzed the events with high waveform correlations (correlation coefficient greater than 0.75). Some of these high-correlation events have the properties of propagating waves; Pi 2 waveforms at the ground stations and the satellite were synchronized with each other when the data were shifted by onset delays. The results of the statistical analysis show that 87% of the Pi 2 onsets at a ground station (Kuju, GMLat. = 26.13°, GMLon. = 202.96°) were delayed from the Pi 2 onsets at ETS-VIII, and the average of the delay times was 29 sec. This clearly shows Pi 2 onsets (initial perturbations of Pi 2) propagated from the geosynchronous altitude to the low-latitude ground. The delay times tended to be larger around the midnight sector than around the dawn and dusk sectors. These results are consistent with two-dimensional propagation of fast waves estimated by the model of Uozumi et al. (J Geophys Res 114:A11207, 2009). The delay times are nearly identical to the travel time of fast waves from geosynchronous altitude to the low-latitude ground, and the local time variation of the delay shows the azimuthal propagation along the geosynchronous orbit. We conclude that the initial compressional perturbations of Pi 2 waves propagate radially and longitudinally as a fast wave in the inner magnetosphere.

Pi 2 geomagnetic pulsations are defined as impulsive geomagnetic oscillations with a period of 40 to 150 s (Jacobs et al.1964). Pi 2 pulsations occur during substorms and other auroral activities such as pseudo breakups and poleward boundary intensifications (Saito1969; Voronkov et al.2003; Lyons et al.1999). Although it is accepted that Pi 2 pulsations (or their source disturbances) are generated in the nightside magnetosphere, their generation and propagation mechanisms have not been determined conclusively.

Simultaneous observations in the magnetosphere and on the ground have greatly improved our understanding of low-latitude Pi 2 pulsations. Using data from the geosynchronous ATS-6 satellite and mid-latitude ground stations, (Sakurai and McPherron1983) found that Pi 2 onsets were almost simultaneous (within 1 min). (Yumoto et al.1989) reported that a magnetic disturbance in the current disruption region (approximately 9 RE) was observed about 1 min earlier than Pi 2 pulsations observed at low-latitude stations. They interpreted the delay time as the travel time of the compressional disturbance. (Takahashi et al.1995) found that Pi 2 pulsations in the inner magnetosphere had high coherence with H-component Pi 2s at a low-latitude ground station. The high-coherence Pi 2s were detected only in the nightside inner magnetosphere (L ≤ 5). They suggested that low-latitude compressional Pi 2 pulsations (their oscillations appear in the compressional component) had a radial standing structure in the plasmasphere due to a cavity mode resonance (e.g., Sutcliffe and Yumoto1989; Yeoman and Orr1989; Nosé2010). (Takahashi et al.2003) analyzed electric and magnetic field perturbations of Pi 2s observed by the Combined Release and Radiation Effects Satellite (CRRES) in the inner magnetosphere. They found the radial standing structure of the fast wave not only inside the plasmasphere but also just outside the plasmasphere. On the other hand, (Kepko and Kivelson1999) and (Kepko et al.2001) compared bursty bulk flows (BBFs) (Baumjohann et al.1990; Angelopoulos et al.1992) with Pi 2 pulsations at geosynchronous orbit and at ground stations. The waveforms of each compressional Pi 2 wave and earthward BBF were well correlated when the data were shifted by the delay times. These studies argued that low-latitude Pi 2 pulsations are generated by the forced oscillation of Pi 2 periodicities associated with the braking of BBFs. The relationship between BBFs and Pi 2 pulsations was subsequently investigated statistically (Kim et al.2007; Hsu et al.2012). Recently, (Murphy et al.2011) showed an event, where a low-latitude Pi 2 was observed approximately 60 s earlier than a periodic BBF, and their waveforms were well correlated. They speculated that the fast wave and the BBF were simultaneously launched from the reconnection region, and the fast wave propagated toward Earth faster than the BBF.

Statistical studies of timing relationships have revealed the propagation morphology of ground Pi 2 pulsations. (Uozumi et al.2000) studied latitudinal propagation characteristics of Pi 2s by examining delays of wave packets. They found that the arrival time of Pi 2s at low-latitude stations are earlier than at a station near the plasmapause and later than at a station in the polar region. (Uozumi et al.2004) investigated the magnetic local time (MLT) dependences of the delay time of Pi 2s detected at two pairs of latitudinally and longitudinally separated high-latitude stations. They concluded that the most probable center of azimuthal propagation is 22.5 MLT in the high-latitude region. (Uozumi et al.2009) showed that the timing of Pi 2s in the high-latitude region varied with vector components of magnetic data. However, there have only been a few statistical studies investigating the timing of Pi 2 onsets between satellites and ground stations.

Although it is believed that Pi 2 pulsations or their source disturbances propagate as fast-mode waves in the low-latitude magnetosphere, direct evidence of Pi 2 propagation and its morphology have not been statistically established. In this paper, we statistically investigate onset delay times of compressional Pi 2s between the geosynchronous satellite ETS-VIII and the low-latitude ground stations using cross-correlation analyses. We also used a numerical model developed by (Uozumi et al.2009) to estimate the travel time of the fast-mode wave in the magnetosphere. From both the observations and model estimation, we discuss the propagation of compressional Pi 2 waves in the low-latitude magnetosphere.

Analysis procedure

Data description

We used the magnetic data from the geosynchronous ETS-VIII satellite and the low-latitude ground observatories. The ETS-VIII was launched on December 18, 2006 and placed in geosynchronous orbit at 146°E geographic longitude (Koga and Obara2008). The projection point of ETS-VIII has been located in the geomagnetic southern hemisphere (GMLat. = -10.8°, GMLon. = 217.5°). The coordinate system of the ETS-VIII magnetic field data is orthogonal, consisting of P, E, and N coordinates; P is parallel to the rotational axis, E points radially earthward, and N is azimuthally eastward. We also used ground magnetic data from the low-latitude magnetic data acquisition system (MAGDAS) stations (Yumoto and The MAGDAS Group2006; Yumoto2007) and from station Honolulu (HON) operated by the US Geological Survey (USGS). All the data have temporal resolutions of 1 s. The station codes and the coordinates are listed in Table1. In this study, since we were interested in the compressional component of Pi 2s, we used the P (parallel to the rotational Earth axis) component of the satellite data and the H (geomagnetic northward) component of the ground data to compare the onset timings.

Table 1

Station codes and coordinates of the ground observatories

Station name

Station code

Geographic

Geographic

Geomagnetic

Geomagnetic

L

latitude

longitude

latitude

longitude

Popov Island

PPI

42.98

131.73

36.27

203.74

1.56

Kuju

KUJ

33.06

131.23

26.13

202.96

1.24

Muntinlupa

MUT

14.37

121.02

6.79

192.25

1.01

Yap

YAP

9.50

138.08

1.49

209.06

1.00

Darwin

DAW

-12.41

130.92

-21.91

202.81

1.16

Learmonth

LMT

-22.22

114.10

-33.59

185.11

1.46

Townsville

TWV

-19.63

146.86

-28.73

220.30

1.30

Melbourne

MLB

-38.36

145.18

-49.46

222.51

2.37

Honolulu

HON

21.32

202.00

21.64

269.74

1.16

Fayum

FYM

29.18

35.50

25.76

112.65

1.23

Event selection

We used H-component magnetic data at the ground station Kuju (KUJ) to identify events used in this study. To compare onset timings, we require an undisturbed background before Pi 2 onsets. Therefore, data from KUJ were processed into time-series data of Pi 2 wave energy (Uozumi et al.2000), defined as an envelope of ΔH2 (ΔH: filtered H with a range of 40 to 150 s) to evaluate the background disturbance and to estimate the maximum energy value and time of each Pi 2. We defined the background disturbance as the average Pi 2 wave energy 5 to 10 min before the time of maximum energy. Potential Pi 2 events at KUJ were identified, where the maximum energy was greater than 0.01 nT2 and the background disturbance was smaller than 10% of the maximum energy. We further removed events with large noise and other types of disturbances by visual examination. A total of 344 Pi 2 events were selected during the period from December 1, 2009 to December 1, 2010, at local time of ETS-VIII (LT = UT + 9.7 h) between 19 and 04 LT. Of the 344 events, 15%, 48%, 27%, 6%, and 4% had Kp values of 0, 1, 2, 3, and 4, respectively.

Calculation of delay times

To calculate onset delay times, we used a cross-correlation analysis. The calculation procedures were as follows. First, we applied a bandpass filter (filter range set to 40 to 150 s). Next, we set a fixed segment of 5 min from Tmax - 120 s to Tmax + 180 s in the filtered data, where Tmax is the time when the Pi 2 wave energy is maximum. Finally, we calculated a cross-correlation function (R(τ)) for a range of time shifts, -90 ≤ τ ≤ 90 s. We defined τ when R(τ) is a maximum peak as a delay time (ΔT). To avoid analytical errors, the delay time was only used in the analysis when R(ΔT) was greater than 0.75, referred to as a ‘high-correlation event’ in this paper.

Event examples

This section provides examples of Pi 2 events. Figure1a shows the filtered P-component data from ETS-VIII and the H-component data from the low-latitude ground stations (at 13:20 UT on August 14, 2010), and the satellite and the ground stations around the 135°E meridian (except HON and Fayum (FYM)) were in the premidnight sector. The first positive peak of the Pi 2 oscillation at ETS-VIII appeared at 13:19:30 UT (red vertical line), whereas the first simultaneous positive peaks at all the ground stations appeared at 13:20:10 UT (black vertical line). The correlation coefficients and the delay times between the KUJ data and other stations are shown in Table2. A positive delay time indicates that a Pi 2 wave was delayed from that of KUJ. As the delay times between KUJ and other ground stations are between -6 and 5 s, the difference between the ground stations is within 11 s. The oscillations on the ground after the onsets are almost in-phase. The small onset time delay and phase delay among the nightside low latitudes are consistent with previous studies (Kitamura et al.1988; Shinohara et al.1997; Nosé et al.2006). Figure1b shows the filtered H-component at KUJ without a time shift and the P-component data at ETS-VIII with a 42-s positive time shift. The onset timings for both Pi 2s are simultaneous, and the waveforms were well correlated (R(ΔT) = 0.97).

Figure 1

Plots of geomagnetic data around 13:20 UT on August 14, 2010.(a) Filtered P-component data from ETS-VIII (ETS P) and H-component data from low-latitude ground stations around 13:20 UT on August 14, 2010 (filter range of 40 to 150 s). (b) Filtered H-component data from KUJ without a time shift and P-component data from ETS-VIII with a positive time shift of 42 s.

Table 2

Correlation coefficients and delay times between KUJ and other stations and ETS-VIII

Station/satellite

R(ΔT)

ΔT

ETS-VIII

0.97

-42

PPI

0.97

4

MUT

0.95

2

YAP

0.99

5

DAW

0.99

5

LMT

0.99

5

TWV

0.97

2

MLB

0.96

-1

HON

0.93

-4

FYM

0.88

-6

A positive delay time indicates that Pi 2s were delayed from the Pi 2 at KUJ.

Figure2 shows the waveform examples of eight high-correlation events. Each event was observed at different local times. The upper panel of each example shows the non-filtered magnetic variation (without any time shift), and the lower panels show the filtered magnetic variations (the filter range of 40 to 150 s) with a time shift of ΔT. Except for event A, Pi 2s at ETS-VIII start 20 to 56 s earlier than Pi 2s at KUJ. It is evident from the lower panel (with time shifts) of all the events that Pi 2 waveforms at ETS-VIII and KUJ are well synchronized. The oscillations in the upper panel (without time shifts) of events D, E, and F appear to be out of phase. It should be noted that one cannot identify from a two-point magnetic observation whether these are due to resonance structures of standing waves or propagation delays (this will be discussed later).

Figure 2

Waveform examples of eight high-correlation events observed at different local times. Upper panels of each example show the non-filtered magnetic variation without a time shift, and lower panels show the filtered magnetic variation (the filter range of 40 to 150 s) with a time shift of ΔT. Horizontal arrows in the lower panels indicate the shift directions, with time shifts provided on the arrows.

Statistical analysis

We statistically investigated delay times between the P-component Pi 2s at ETS-VIII and the H-component Pi 2s at KUJ for high-correlation events (0.75 ≤ R(ΔT)). As a result of the cross-correlation analyses, 131 events (38%) have R(ΔT) values greater than 0.75. Of the 131 high-correlation events, 24%, 54%, 18%, 4%, and 0% had Kp values of 0, 1, 2, 3, and 4, respectively. Of the high-correlation events, 78% occurred during periods of very low geomagnetic activity (Kp = 0,1), whereas only 53% of the low-correlation events (0.75 ≥ R(ΔT)) occurred during these periods. Therefore, the high-correlation events tend to occur during low geomagnetic activity.

Figure3a shows the distribution of local times at ETS-VIII (LT = UT + 9.7 h) when high-correlation events occurred. Most events occurred between 23 and 24 LT, with more events seen in the premidnight sector than in the postmidnight sector. Figure3b shows the distribution of delay times of high-correlation events. The average delay time of all the events was -29 s (the median is -35 s), with 87% of events showing negative delay times that indicate that the Pi 2 pulsations at ETS-VIII occurred before the Pi 2 at KUJ. Figure3c shows the local time dependence of the delay times between ETS-VIII and KUJ. Because there are some outliers, we show not only the average values and standard deviations but also the median values in each hourly bin. The delay times tend to be longer around the midnight sector than around the dawn and dusk sectors. The dashed red line represents the delay time as a function of the local time calculated by the numerical model of (Uozumi et al.2009).

Figure 3

Analysis of 131 high-correlation events.(a) Distribution of local time at ETS-VIII (LT = UT + 9.7 h). (b) Distribution of delay times. (c) Local time dependence of the delay times.

Model estimation of fast-wave travel time

We estimated the fast-wave travel time using the model presented in (Uozumi et al.2009). This model is based on a three-dimensional magnetohydrodynamic (MHD) propagation model for Pi 2 pulsations developed by (Uozumi et al.2007) and applied to a realistic plasmapause model which was developed by (Gallagher et al.2000). The propagation path of MHD waves is based on the mode conversion from fast waves to Alfvén waves (Tamao1964; Itonaga and Yoshikawa1996). The Alfvén wave reaching the ground station is excited at the equatorial intersection of the field line by inward fast waves reaching the intersection. We focused on the two-dimensional propagation of the fast waves from a local source in the magnetic equatorial plane.

In this study, the wave source is set at a radial distance of r = 10 RE and the longitudinal width of LT = 22 to 24 h. (McPherron et al.2011) have suggested that plasma flows are clustered around 23 LT and likely decelerated at approximately 10 RE. It is believed that the position of aurora breakups is related to the Pi 2 source (or energy source) (e.g., Kuwashima1978; Rostoker and Samson1981). Breakup regions are located around 23 MLT, and their longitudinal width is about 2 MLT on average (Liou et al.2002). For these reasons, we assumed that the Pi 2 source is not a point source but a longitudinal wide source. We set Kp equal to 1 because this is where most of the high-correlation events occurred. He+ and O+ densities were set to 1%. Figure4 shows the isochron map of wave fronts of the fast wave in the XGSM-YGSM plane (ZGSM = 0) estimated from the model. The wave fronts expand anisotropically, depending on the Alfvén velocity distribution. The propagation speed dramatically decreases in the plasmasphere due to the large mass density. In the lower altitude region (r ≤ 2 RE), the speed dramatically increases due to strong magnetic fields. On the ground on the nightside (r = 1 RE), the travel time of the wave front to longitudinally separated points is not significant (< 6 s).

Figure 4

Isochron map showing wave fronts of fast-mode waves illustrated by the model of (Uozumi et al.2009). The wave source is set at r = 10 RE and LT = 22 to 24 h. Blue and red contours are time intervals of 2 and 10 s, respectively. The dashed circle indicates geosynchronous orbit. The green ellipse illustrates the plasmapause.

The delay time as a function of local time at ETS-VIII calculated by the model of (Uozumi et al.2009) (ΔTU 09 (LT)) is as follows:

ΔTU09(LT)=TETS(LT)-TKUJ(LT),

(1)

where TETS (LT) and TKUJ (LT) are the travel time from the source to ETS-VIII and KUJ, respectively, represented by the dashed line in Figure3c. The minimum peak of ΔTU 09 (LT) appears at 23.5 LT and its value is -42 s. The value of ΔTU 09 (LT) decreases more notably on the dusk side than on the dawn side.

The observational results show Pi 2 onsets (initial perturbations of Pi 2) propagated from ETS-VIII to the low-latitude ground stations because in most cases (87%), Pi 2s at ETS-VIII started before Pi 2s at KUJ. Assuming that the location of the Pi 2 source is greater than the geosynchronous altitude, the waves reach the ground through ETS-VIII when the ETS-VIII is in the midnight sector. Therefore, the delay time in the midnight sector is probably equal to the wave travel time from the geosynchronous altitude to the low-latitude ground. The delay time in the midnight sector is calculated to be -37 s, the average within the 23 to 01 LT. We obtained an average propagation velocity of 920 km/s. This is in the same order of the typical Alfvén velocity in the inner magnetosphere. In addition, only fast-mode waves can directly propagate across the ambient magnetic field. Thus, we consider that initial perturbations of Pi 2s propagate as fast-mode waves.

The LT dependence of delay times would be the result of the longitudinal propagation of the fast waves in the magnetosphere. As the isochron map (Figure4) shows, the arrival times of the fast waves are not very different at each longitudinal ground point, whereas at geosynchronous altitude, the arrival times vary with the local time due to a large longitudinal spatial displacement. When the satellite is located at the magnetospheric flanks, the distance from the source to the satellite becomes longer and comparable to the distance from the source to the ground. In that case, it takes more time for the fast waves to travel from the source to the satellite, and consequently, the difference in the arrival times between the satellite and the ground is smaller. We suggest that the longitudinal propagation at the geosynchronous altitude contributes to the LT dependence of delay times.

Many factors appear to be responsible for the propagation delay of fast waves. The shape of the plasmasphere would be the most important factor. Since the Alfvén velocity dramatically decelerates inside the plasmasphere, the delay time becomes near zero or even positive (earlier at ground stations) when the satellite is inside the plasmasphere. The ratio of heavy ions to proton mass densities is also an important factor. The ratio of O+ ions varies depending on geomagnetic activity levels (Lennartsson and Sharp1982). Because O+ ions are about 16 times heavier than H+ ions, their contribution to plasma mass density is large, and consequently, the Alfvén velocity can be dramatically altered. Equally important, all components that affect the distribution of magnetic fields and plasma mass densities can be responsible for the delay of the fast waves. To better estimate the source region and the travel time for individual events, more detailed information of the magnetospheric conditions is required.

Propagation of initial Pi 2 perturbations or Pi 2 pulsations themselves involves two scenarios of low-latitude Pi 2 generation, a plasmaspheric resonance, and the BBF-driven mechanism. In the case of the plasmaspheric resonance, an initial Pi 2 perturbation would be a source disturbance triggering the plasmaspheric resonance. This source disturbance has been recognized as the impulsive compressional MHD disturbance (Yumoto et al.1989; Fujita et al.2000; Fujita et al.2002). (Fujita et al.2000) conducted a simulation of the plasmaspheric virtual resonance (PVR) and showed that the initial perturbation of magnetic fields was delayed between the inside and outside of the plasmasphere. Both observations and simulations of PVR have shown that a phase reversal occurs near the plasmapause. The 74% of high-correlation events with delay times within half a period of the Pi 2 range (20 to 75 s) found in this study implies that onset delays can contribute to phase delays of approximately 180° between geosynchronous altitude and low latitudes. Events D, E, and F in Figure2 are good examples showing onset delays that correspond to those phase delays of approximately 180°. On the other hand, Pi 2s driven by BBFs are expected to show the characteristics of propagating forced waves. Propagating waves observed at two different points should have high correlation when data are shifted by a propagation delay. Our observations show that there are high-correlation events which have time delays consistent with the fast-wave travel times. Also, there are several events whose waveforms are synchronized, including their irregular features when the data are shifted by the onset delays (for example, Figure1b and events F and G in Figure2). It is plausible that such Pi 2s propagated as a forced fast-mode waves. As stated before, since the delay time corresponds to approximately half the period of Pi 2 pulsations, it is possible that propagation delays contributed to out-of-phase oscillations. However, one cannot conclude from magnetic observation of only two points whether resonance structures of standing waves or propagation delays contribute to these out-of-phase oscillations. Observations by multiple satellites in elliptical orbits (e.g., THEMIS) are potentially useful to reveal the radial structure of propagating delays. The distance from the magnetic equatorial plane to elliptical orbit satellites is variable even if the local time and distance from Earth are the same. Differences in local time and latitude between ground stations and elliptical orbit satellites changes by the hour. Spatial coverage of elliptical orbit satellites is much larger than that of geosynchronous satellites, but the amount of data in each spatial grid is small by necessity. As propagation delay times vary between events, the error relating to a difference in conditions in each case cannot be reduced without enormous data sets.

We have presented the statistical analysis of onset delay times of compressional Pi 2s which have a high correlation between the low-latitude ground stations and the geosynchronous satellite. The results showed that there were events with quite similar waveforms when the data are shifted by the onset delays. While onset timings at all low-latitude stations were simultaneous, onset timings at the geosynchronous satellite were clearly earlier than at the ground stations. Of the Pi 2s at KUJ, 87% were delayed from Pi 2s at ETS-VIII. The delay time tends to be shorter when the satellite is located at the magnetospheric flanks, which indicates that Pi 2s (or initial Pi 2 perturbations) propagate longitudinally from the midnight source region. Propagation of initial Pi 2 perturbations or Pi 2 pulsations themselves involves both plasmaspheric resonance and BBF-driven mechanisms. Although we cannot determine whether Pi 2 pulsations themselves are standing waves or propagating waves, we conclude that, in either case, initial perturbations of compressional Pi 2s propagate as fast waves radially and azimuthally.

Acknowledgements

This work was supported by the JSPS Core-to-Core Program, B. Asia-Africa Science Platforms. MAGDAS magnetic data from MUT and YAP are supported by Mr. Commodore Romeo I. Ho and Dr. David Aranug, respectively. The Honolulu magnetic field data were provided by the USGS Geomagnetism Program. Kp-index data was provided by the World Data Center for Geomagnetism, Kyoto, Website. Financial support was provided by the Ministry of Education, Science and Culture of Japan in the form of Grants-in-Aid for Overseas Scientific Survey (18253005 and 22253007).

Authors' original submitted files for images

Below are the links to the authors’ original submitted files for images.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SI analyzed the data and drafted the manuscript. KY conceived of the study and helped interpret the results. TU provided the simulation results of the numerical model. HK helped revise the draft manuscript. SA and AI helped analyze the data. KK, HM, and TO provided the ETS-VIII magnetic data. RM, VAA, AM, and AL made substantial contributions to acquisition of MAGDAS magnetic data. AY made substantial contributions to conception and interpretation of this study. All authors read and approved the final manuscript.

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.