041_053 adg v–5 n01.pdf ANNALS OF GEOPHYSICS, VOL. 45, N. 1, February 2002 41 GPS detection of the instantaneous response of the global ionosphere to strong magnetic storms with sudden commencement Edward L. Afraimovich, Eugene A. Kosogorov, Ludmila A. Leonovich, Oleg S. Lesyuta and Igor I. Ushakov Institute of Solar-Terrestrial Physics SD, Russian Academy of Science, Irkutsk, Russia Abstract Using a new technology for global GPS detection of ionospheric disturbances, GLOBDET, it has been established that a drastic increase in the time derivative of the magnetic field strength during magnetic storms is accompanied by an almost simultaneous decrease in mid-latitude total electron content on the entire dayside. The corresponding correlation coefficient is not below −0.8; the delay with respect to the time of a magnetic storm sudden commencement is about 3-10 min. This is most pronounced for magnetic storms with a well- marked sudden storm commencement. The sudden storm commencements presented in the paper were observed during the initial storm phase. The analysis reported here was made for a set of from 90 to 300 GPS stations for 10 days in 1998-2001 with a different level of geomagnetic activity (Dst from −6 nT to −295 nT, and K p from 0 to 9). The «simultaneous» total electron content response for the events under consideration was 0.1-0.4 TECU, and the travel velocity of the disturbance from the dayside to the nightside was in the order of 10-20 km/s. Results obtained are consistent with earlier ionospheric parameter measurements obtained using high temporal resolution methods. 1. Introduction The ionospheric response to a Sudden Storm Commencement (SSC) should be set off from the large number of phenomena observed during different phases of magnetospheric storm development. Many morphological research efforts have been made to puzzle out the dynamic nature of the processes accompanying SSC. SSC response characteristics in mid-latitude VLF- emissions were considered in Mullayarov and Muzlov (2000). Observations of post-SSC ionospheric fluxes and particle precipitation events were presented in Nielsen and Honary (2000). The solar wind velocity and density effects on the cyclic behavior of the SSC occurrence frequency were investigated in Makarov (1994). Distinctive decreases in maximum ionization (N m ) and in ionization at fixed levels (N h ) after SSC were observed for a strong magnetic storm of July 11-19, 1959 (Potapova and Shapiro, Mailing address: P r o f . E d w a r d L . A f r a i m o v i c h , Institute of Solar-Terrestrial Physics SD, Russian Academy of Sciences, P.O. Box 4026, Irkutsk, 664033, Russia; e-mail: afra@iszf.irk.ru Key words global ionospheric response – sudden storm commencement – magnetic storms 42 Edward L. Afraimovich, Eugene A. Kosogorov, Ludmila A. Leonovich, Oleg S. Lesyuta and Igor I. Ushakov 1961). A valley in the diurnal variation of the F2-layer critical frequency (f0 F2), and a drastic increase in the height of the F2-layer maximum (h max F 2 ) occurred almost simultaneously with SSC on the entire day-side (for mid-latitude stations). The valleys during evening SSC are hidden by an abrupt diurnal change in the parameter f0 F2; therefore, a scatter of the delay time τ between the f 0 F 2 minimum and the SSC time was observed. The duration (∆T ) of the f0 F2 response to the SSC varied from 1 to 22 min. Table I presents the following charac- teristics: magnetic storm date; reference to a relevant publication; τ - delay between the SSC time and the ionospheric response; ∆T - duration of the ionospheric response to the SSC; ∆ | A| - amplitude increment relative to quiet values; V - displacement velocity of disturbances, and n - largest number of stations involved in the investigation. The post-SSC ionospheric response for the magnetic storm of June 15, 1965 was in- vestigated in great depth by Shashunkina (1968, 1972), based on data from a number of ion- ospheric stations located at different sites of the globe. Using data on f 0 F 2 , the F 2 -layer maximum electron density (N max F 2 ), the maximum height hmaxF2, the F2-layer thickness y m , and electron density distribution profiles with height (N h ), it was found that the delay of the maximum of the ionospheric response to SSC varied from 0 to 90 min. The response duration according to f0 F2 data varied from 90 (low solar activity) to 180 min (high solar Table I. Characteristics of the ionospheric response to SSC. N Date Reference Technique Parameters τ, min ∆T, min ∆| A| V, m/s n 1 11-19.07-1959 Potapova and Sphiro I f0 F2, N(h) 0-138 1-22 0.2-2 MHz − 14 (1961) NmF2, Nh 2 15.06.1965 Shashunkina I f0 F2, N(h) 0-90 90-180 1-1.3 MHz 0.72-2.5 20 (1968, 1972) hmF2, NmF2, h′F2, ym 3 25-26.05.1990 Buonsanto FRLPS, TEC, NmF2, 108-300 102-360 10-32 TECU, − 4, (1995) ISR h m F 2 180-300 150 km 1 4 25-26.09.1958 Rüster (1965) I hmF2, f0 F2 150 120 100 km − 1 5 13-14.09.1959 Rüster (1965) I hmF2 180 120 150 km − 1 6 3-4.10.1959 Rüster (1965) I f0 F2 90 90 100 km − 1 7 8.07.1958 Bowman I ∆h′F 3-3.5 180-540 20-100 km 0.37-0.72 5 (1965) hmF2 8 12-14.03.1989 Huang and Cheng I f0 F2, TEC 0-0.6 30-60 5 MHz, − 1, (1991) FRLPS 65 TECU 3 9 15.02.1967; Huang et al. FDS ∆f 0-0.07 1-4 0.5-2 Hz − 2 13.03.1970 (1973) 10 6.01.1998; − GLOBDET TEC 4-10 15-30 0.08-0.38 10-20 90-300 23.04.1998; 6.04.2000; 8.06.2000; 13-15.07.2000; 31.03.2001; 4-11.04.2001 I = ionosonde; FRLPS = Faraday rotation of the linearly-polarized signal; ISR = incoherent scatter radar; FDS = frequency Doppler shift measurements along the HF path; 1TECU = 1016m−2. 43 GPS detection of the instantaneous response of the global ionosphere to strong magnetic storms with sudden commencement activity) depending on the station’s latitude. The delay of the ionospheric response to SSC increased from high to low latitudes. Some publications (Huang and Cheng, 1991; Buonsanto, 1995) contain estimates of variations of ionospheric parameters after SSC made on the basis of data from Incoherent Scatter Stations (ISS) and by determining the value of the Total Electron Content (TEC) from Faraday effect measurements of the signal from geostationary satellites (GOES2, GOES7). The delay of the ionospheric response to SSC for the magnetic storm of May 25-26, 1990 was 108-300 min (according to TEC measurements) and 180-300 min (according to data from ISS) (Buonsanto, 1995). The response duration ∆T was 102-360 min. For the magnetic storm of March 12-14, 1989, τ ranged from 0 to 36 min according to TEC and f 0 F 2 data. The response duration ∆T was 30-60 min (Huang and Cheng, 1991). Information on the ionospheric response to SSC for the magnetic storms of September 25- 26, 1958, September 13-14, 1959, October 3-4, 1959, and July 8, 1958 according to ionosonde data is reported by Bowman (1965) and Rüster (1965). The delay τ for these storms averaged 90-210 min. The ionospheric response to an instantaneous change in the magnetic field lasted 90-540 min. The ionospheric response to a SSC of the magnetic storms of February 15 and September 19, 1967; October 29, 1968; May 14, 1969; April 28, 1969, and of March 31, 1970, was inves- tigated by measuring the frequency Doppler shift along an HF path (Huang et al., 1973). The delay τ between a maximum frequency Doppler shift and the SSC time was 0-4 min; the response duration was 1-4 min. Of all the above- mentioned data these results are distinguished for the highest time resolution and hence for the highest reliability of determination of the response delay τ. Investigations of the ionospheric response to SSC made by a number of authors to date, contributed significantly to the study of the phenomena described above. A demerit char- acteristic for most investigations is the low spatial and temporal resolution of measuring facilities (time steps of 5 or 15 min; the largest number of ionospheric stations involved in the study of the same event does not exceed 20). For that reason, in spite of the large body of experimental evidence available, it is impossible to obtain reasonably reliable estimates of the main parameters of ionospheric response to SSC. Because of this, the results presented in table I are most likely to refer to different manifestations of geomagnetic disturbances in the mid-latitude ionosphere. Thus, in some cases (lines 2 and 7 of table I) the ionospheric parameter variations were caused by large-scale traveling ionospheric disturbances (LS TID) with typical time periods of 1-2 h and wavelengths in the order of 1000 km (Bowman, 1965; Shashunkina, 1968, 1972). It is generally believed that LS TIDs are manifestations of Acoustic-Gravity Waves (AGW) whose generation regions lie in equatorial zones of the northern and southern hemispheres. These concepts received further support through recent experiments using GPS data (Ho et al., 1998; Afraimovich et al., 2000b). The LS TID response delay of about 1-2 h is consistent with the idea that TIDs are generated in auroral regions and travel equatorward at a velocity of about 300-400 m/s. In other cases (lines 9 and 10 of table I) the response delay is too small to be consistent with the above mechanism associated with the generation and propagation of AGW. A global spatial averaging of TEC variation, obtained using data from the international GPS network, seems to have been used for the first time by Afraimovich et al. (2000a) in investigating the global ionospheric response to SSC. An «instantaneous» ionospheric response to a sudden commencement, SSC, of a strong magnetic storm of April 6, 2000 was detected; it appeared as a single negative TEC disturbance with a duration of about 20 min, with about 200- s delay with respect to the SSC at mid-latitudes, almost simultaneously on the entire dayside. The study reported by Afraimovich et al. (2000a) used the new technology, GLOBDET, developed at the ISTP SD RAS for global detection of ionospheric disturbances of natural and artificial origins using data from the international network of two-frequency multi-channel receivers of the navigation GPS system which substantially improves the sensitivity and space-time resolution of experiment (Afraimovich, 2000). 44 Edward L. Afraimovich, Eugene A. Kosogorov, Ludmila A. Leonovich, Oleg S. Lesyuta and Igor I. Ushakov The objective of this paper is to obtain statistically more significant estimates of the main parameters of the response to SSC for magnetic storms with a different level of geomagnetic activity. The SSC presented in the paper were observed during the initial storm phase. 2. The method for determining the global ionospheric response to a sudden storm commencement We now outline the procedure of processing GPS data which is directly used in this paper. The input data in our experiment are represented by series of «oblique» TEC values along the beam to the GPS satellite at selected receiving points, as well as the corresponding series of values of the elevation θ (t) and azimuth α (t) of the beam calculated from our developed CONVTEC program to convert standard (for the GPS system) RINEX-files available through the Internet (Gurtner, 1993). The variations of the «oblique» TEC I (t) are determined on the basis of phase measurements at each of the spatially separated two-frequency GPS receivers using the known formula from Hofmann-Wellenhof et al. (1992) (2.1) where L 1 λ 1 and L 2 λ 2 are the increments of the radio signal phase path caused by the phase delay in the ionosphere (m); L1 and and L2 stand for the number of complete phase rotations, and λ 1 and λ 2 are the wavelengths (m) for the frequencies f1 = 1575.42 MHz and f2 = 1227.60 MHz, respectively; const is some unknown initial phase path (m); and nL is the error in determining the phase path (m). The sensitivity of phase measurements in the GPS system permits irregularities and wave processes to be detected over a wide range of amplitudes (as large as 10–3 of the diurnal TEC variation) and periods (from several days to a few minutes). Standard data available via the Internet have a time resolution of 30 s. Ionospheric data provide their data with a time resolution of 5-15 min. Hence the GPS data have higher time resolution compared with ionosondes and incoherent scatter stations. To determine the characteristics associated with disturbances involves selecting several series of TEC I(t) measurements with a duration of 1 h at least. To exclude the influence of the signal reception conditions, we used only observations with satellite elevations θ (t) larger than 30°. For the purpose of eliminating the variations of a regular ionosphere as well as trends introduced by the satellite motion, use is made of the procedure of removing the linear trend by preliminarily smoothing the initial series with the selected time window of a duration in the order of 30 min. The GLOBDET software package used in this study performs a global coherent accu- mulation of TEC variations over the time interval chosen for the analysis (Afraimovich, 2000) (2.2) where dI i is the filtered TEC series; and n is the number of beam. The correction coefficient K i is required for converting the «oblique» TEC to an equivalent «vertical» value (Klobuchar, 1997) (2.3) where R z is Earth’s radius; and hmax is the height of the ionospheric F2-layer maximum. 3. Geometry of the experiment The investigation reported in this paper is based on data from the global network of GPS receiving stations available via the Internet. At present (as of February 2001), data from no fewer than 800 GPS stations can be accessed through the Internet. Figure 1a illustrates the geometry of part of the global network used in nL+ + ]const I f f f f L L= − −( ) +[ 1 40 308 1 2 2 2 1 2 2 2 1 1 2 2. λ λ S t dI Ki i i n = × = ∑ 1 ( ) K R R h i z z i= +         θcos arcsin cos max ( ) 45 GPS detection of the instantaneous response of the global ionosphere to strong magnetic storms with sudden commencement 09:10 17:10 01:10 09:10 LT 0 60 120 180 240 300 360 90 60 30 0 -30 -60 -90 La tit ud e, d eg . Longitude, deg. 0 60 120 180 240 300 360 90 60 30 0 -30 -60 -90 La tit ud e, d eg . a b Fig. 1a,b. The geometry of the part of the global GPS network used in this study when analyzing the ionospheric response to the SSC of the strong magnetic storms of June 8, 2000 (263 stations) (a). Dots correspond to the location of GPS stations. The upper scale indicates the local time LT, corresponding to the time 09:10 UT. Panel (b) shows the coordinates of subionospheric points for the height of the F 2 -layer maximum h max = 300 km for all satellites visible at the SSC time June 8, 2000, for each of the GPS stations marked on panel a (a total of 1325 LOS). 46 Edward L. Afraimovich, Eugene A. Kosogorov, Ludmila A. Leonovich, Oleg S. Lesyuta and Igor I. Ushakov this case study of the magnetic storm of June 8, 2000 (263 stations). As is evident from fig. 1a, the set of stations we selected from the part of the global GPS network available to use, covers rather densely North America and Europe, and to a much lesser extent Asia. GPS stations are more sparsely distributed on the Pacific and Atlantic Oceans. However, such coverage of the surface of the globe makes it possible, even today, to tackle the problem of global detection of disturbances with hither to unprecedented spatial accumulation. This ensures the number of statistically independent series which is, as a minimum, by two orders of magnitude larger than could be realized when taking records with other measuring facilities. Thus, in the Western hemisphere the corresponding number of stations can reach no less than 500 today, and the number of beams to the satellite can be no less than 2000-3000. The upper scale in fig. 1a shows the Local Time (LT), corresponding to 09:10 UT (SSC time for the magnetic storm of June 8, 2000). For a variety of reasons, slightly differing sets of GPS stations were selected for different events to be analyzed; however, the geometry of experiment for all events was virtually identical. The station coordinates are not given here for reasons of space. This information may be found at the electronic address http://lox.ucsd.edu/cgi- bin/allCoords.cgi? Figure 1b presents the locations subionospheric points for the height hmax = 300 km for the whole set of GPS satellites during the SSC of June 8, 2000, for each of the GPS stations marked off in fig. 1a (the number of beams is 1325). 4. Ionospheric response to a sudden commencement of the strong magnetic storm of April 6, 2000 Figure 2a-d plots the TEC time dependence I(t) on the dayside for April 6, 2000 (a), and the dI(t) variations obtained upon removing the linear trend from the I(t) series and smoothing with the 30-min time window (b) for stations ALGO (PRN02, thick line) and OHIG (PRN01, thin line). The same applies for the nightside for stations IRKT (PRN21, thick line), KERG (PRN22, thin line) (c) and (d). As is evident from fig. 2b, the decrease of dI(t) relative to the SSC time, which is characteristic for the ALGO and OHIG stations, occurs with a delay τ of about 5 min. For the IRKT and KERG stations (fig. 2d), there is a delay of about 32-37 min between the valley in S(t) and the SSC time. The geographic coordinates of GPS stations OHIG and ALGO are shown in fig. 2b, and those for stations IRKT and KERG are indicated in fig. 2d. For the time interval 16:00-18:00 UT, April 6, 2000, fig. 3a plots the variations in magnetic flux at geostationary orbit of the GOES10 satellite (135°W), fig. 3b the variations of the H - component of the magnetic field H(t), and fig. 3c the time derivative of the H-component (thin line) for station Irkutsk (52.20°N; 104.30°E). Results of a global summation of S(t) are presented for the dayside (472 trajectories) in fig. 3d. The dependence S(t) for the nightside (245 trajectories) is plotted in fig. 3e. For comparison, fig. 3c plots the inverted curve S(t) for the dayside (thick line). The SSC time (16:39 UT) in fig. 3a-e is shown by the vertical dashed line. An analysis of the TEC data for our entire chosen set of GPS points for the time interval 16:30-17:30 UT, including the SSC time, shows a negative disturbance of a duration of about 20 min. The delay σ between a maximum of the dH(t)/dt and a maximum of the inverted S(t) series was 6 min. The delay τ between the SSC time and a maximum of the S(t) curve was 7 min. Figure 4a plots the time of a minimum t min of the filtered TEC series versus latitude, obtained for each beam; the amplitude of the TEC response to the SSC of the strong magnetic storm of April 6, 2000 on the dayside; and the distribution of t min and of the response amplitude A as a function of Local Time (LT) (for the time interval 16:24-17:36 UT). The SSC time in fig. 4a,c is shown by the thin horizontal line. The scatter of the time t min is due to the fact that there exist background TEC fluctuations with periods similar to periods of oscillations caused by the SSC; therefore, to improve the signal/noise ratio we apply a global coherent accumulation of dI(t) series. 47 GPS detection of the instantaneous response of the global ionosphere to strong magnetic storms with sudden commencement Fig. 2a-h. Time dependencies of TEC I(t) on the dayside on April 6, 2000 (a), and variations of the dI(t) with the linear trend removed and a smoothing with the 30-min time window (b) for stations ALGO (PRN02 - thick line) and OHIG (PRN01 - thin line). Same on the nightside for the station IRKT (PRN21 - thick line) and KERG (PRN22 - thin line) (c) and (d). Time dependencies of TEC I(t) on the dayside on June 8, 2000 (e), and variations of the dI(t) (f ) for stations PFAN (PRN03 - thick line). Same on the nightside for the station MEM2 (PRN02 - thick line) and RIOG (PRN01 - thin line) (g) and (h). All panels show GPS station names, PRN numbers of the GPS satellites and geographical coordinates of stations. 0 1 0 2 0 3 0 4 0 I ( t ) , T E C U O H I G A L G O - 1 .8 - 1 .2 - 0 .6 0 0 .6 1 .2 d I ( t ) , T E C U 0 1 2 3 4 I ( t) , T E C U 1 6 :3 0 1 7 :0 0 1 7 :3 0 T i m e , U T - 0 .4 - 0 .2 0 0 .2 0 .4 d I ( t ) , T E C U I R K T K E R G 0 4 8 1 2 I ( t ) , T E C U - 0 .1 6 - 0 .0 8 0 d I ( t ) , T E C U P F A N S S C ( - 6 3° N ; 3 0 2 ° E ) ( 4 6° N ; 2 8 2 ° E ) ( - 4 9 .3° N ; 7 0 ° E ) ( 5 2° N ; 1 0 4 .3° E ) 0 2 4 6 I ( t ) , T E C U M E M 2 R I O G 0 9 :0 0 0 9 :3 0 1 0 :0 0 T i m e , U T - 0 .2 - 0 .1 0 0 .1 0 .2 d I ( t ) , T E C U ( 3 5 .4 6 ° N ; 2 7 0 ° E ) ( - 5 3 .8° N ; 2 9 2 .2° E ) S S C (47.5° N; 9.8° E) April 6, 2000 June 8, 2000 a b c e d g h f 48 Edward L. Afraimovich, Eugene A. Kosogorov, Ludmila A. Leonovich, Oleg S. Lesyuta and Igor I. Ushakov An examination of figs. 3a-d and 4a-d, suggests the conclusion that on the dayside of the Earth the largest value of the TEC response amplitude was 0.4 TECU (1-2% of the background value of the TEC I0; 1TECU = 10 16 m–2). The value of I0 was determined from standard IONEX-files (Mannucci et al., 1998), available via the Internet, by averaging the TEC values for the latitude, longitude and time ranges of our interest. The duration of the ionospheric response to the SSC was 20 min. The delay τ between the largest value of the TEC response amplitude and the SSC time, in mid-latitudes on the dayside, was 6 min. In high latitudes the delay τ lasts for more than 15 min. On the nightside, A = 0.1 TECU, and τ = 45 min. 5. Ionospheric response to other magnetic storms Similar results were also obtained for some other magnetic storms: of January 6 and April 23, 1998; June 8 and July 13, 14 and 15, 2000; and of March 31 and April 4, 11, 2001. Data on the magnetic field were received from the GOES spacecraft (the satellite numbers are 9 and 10 for April 23, 1998 and June 8, 2001; March 31 and April 4, 11, 2001, respectively). Figure 2 (right) plots the time dependencies of the TEC I(t) on the dayside for June 8, 2000 (e) and the dI(t) variations (f ) for GPS station PFAN (PRN03). The same applies for the nightside for stations MEM2 (PRN02, thick curve) and RIOG (PRN01, thin curve). It is apparent from fig. 2f that for station PFAN lying on the dayside, the delay τ between the SSC onset time and a minimum dI(t) is about 3 min. On the nightside, τ = 7-9 min for stations MEM2 and RIOG. The geographic coordinates of GPS station PFAN and for stations MEM2 and RIOG are shown in fig. 2f,h respectively. Figure 4e-h shows the same parameters as in fig. 4a-d, but for the magnetic storm of June 8, 2000. The response amplitude A in fig. 4h is plotted as a function of Local Time (LT) for the time interval 08:42-09:42 UT. Figure 5a-e plots the same dependencies as in fig. 2a-h, but for the time interval 08:00-10:00 UT, June 8, 2000. The result of a global coherent summation for the dayside (fig. 5d) is shown for 406 beams. Such a result for S(t) on the nightside (632 beam) is presented in fig. 5e. Figure 6a-e plots the same dependencies as in fig. 2a-h, but for the time interval 15:30-17:30 UT, April 11, 2001. The result of a global coherent summation for the dayside (fig. 6d) is shown for 403 beams. Such a result for S(t) on the nightside (97 beam) is presented in fig. 6e. Fig. 3a-e. The variations in magnetic flux (a) at geostationary orbit of the GOES10 station (135°W), the variations of the H-component of the magnetic field H(t) (b), and of the time derivative dH(t)/dt (c) (thin line) for the time interval 16:00-18:00 UT, April 6, 2000, at station Irkutsk (52.2°N, 104.3°E). The result of a global spatial averaging of S(t) (d) for 472 LOS on the dayside. A similar result for 245 LOS on the nightside (e). For comparison, panel (c) presents the inverted curve S(t) (heavy line) on the dayside. The time of SSC is shown on panels (a-e) by a vertical dashed line. a b c d e 49 GPS detection of the instantaneous response of the global ionosphere to strong magnetic storms with sudden commencement Fig. 4a-h. The latitudinal dependencies (obtained for each LOS) of the tmin (a) and of the amplitude A (b) of the ionospheric response to the magnetic storm SSC on April 6, 2000, on the dayside; the total distributions of tmin (c) and A (d) as a function of local time LT (for the time interval 16:24-17:36 UT). The latitudinal dependencies of the t min (e) and of the amplitude A (f ) of the ionospheric response to the magnetic storm SSC on June 8, 2000, on the dayside; the total distributions of t min (g) and (h) as a function of local time LT (for the time interval 08:42-09:42 UT). The SSC times is shown on panels (a, c, e, g) by a horizontal straight line. - 8 0 - 4 0 0 4 0 8 0 1 6 :3 0 1 7 :0 0 1 7 :3 0 t m i n , U T 6 8 1 0 1 2 1 6 :3 0 1 7 :0 0 1 7 :3 0 t m i n , U T - 8 0 - 4 0 0 4 0 8 0 L a ti t u d e , d e g . - 1 .5 - 1 .0 - 0 .5 0 .0 A , T E C U 6 8 1 0 1 2 T i m e , L T - 2 .0 - 1 .5 - 1 .0 - 0 .5 0 .0 A , T E C U S S C S S C - 8 0 - 4 0 0 4 0 8 0 0 9 :0 0 0 9 :1 1 0 9 :2 3 0 9 :3 5 0 9 :4 7 t m i n , U T 9 1 2 1 5 0 9 :0 0 0 9 :1 1 0 9 :2 3 0 9 :3 5 0 9 :4 7 t m in , U T - 8 0 - 4 0 0 4 0 8 0 L a t i t u d e , d e g . - 0 .4 - 0 .2 0 .0 A ,T E C U 9 1 2 1 5 2 T im e , L T - 0 .6 - 0 .4 - 0 .2 0 .0 A ,T E C U S S C S S C June 8, 2000April 6, 2000 1 8 2 1 1 4 1 6 1 4 1 6 1 8 2 1 a b c e d g h f 50 Edward L. Afraimovich, Eugene A. Kosogorov, Ludmila A. Leonovich, Oleg S. Lesyuta and Igor I. Ushakov An examination of figs. 4a-h and 5a-h suggests that the ionospheric response to the SSC lasted 20 min. The delay σ between a maximum of the inverted curve S(t) and a maximum of the curve dH/dt is 3.5 min. The delay τ between the SSC time and a maximum of the S(t) curve was 4.5 min. On the dayside of the Earth the largest value of the response in amplitude and the delay τ with respect to the SSC time were 0.1 TECU and 2 min, respectively. On the nightside, these values were 0.06 TECU and 8 min, respectively. The travel velocity of the disturbance from mid- to high-latitudes and from the dayside to the nightside averaged 10-20 km/s. Fig. 5a-e. The same as fig. 2a-h, but for the time interval 08:00-10:00 UT, June 8, 2000. The result of a global spatial averaging of S(t) for 406 LOS on the dayside (d). A similar result for 632 LOS on the nightside (e). Fig. 6a-e. The same as fig. 2a-h, but for the time interval 15:30-17:30 UT, April 11, 2001. The result of a global spatial averaging of S(t) for 403 LOS on the dayside (d). A similar result for 97 LOS on the nightside (e). The main values characterizing the magnetic storms being analyzed are presented in table II: Dst at the SSC time (according to the data from the GOES satellite); the SSC time; max dH/dt - maximum value of the derivative H(t); ∆T - time interval under investigation; m - number of averagings (of beams) in the case of a global coherent accumulation;A∑- modulus of a maximum response amplitude for the S(t) series (it is shown on a percentage basis how much our measured value of the TEC amplitude differs from the mean background value of I - 0 calculated from standard IONEX-files); σ - delay between maxima of the series S(t) (inverted values) and 1 0 0 1 2 0 1 4 0 1 6 0 1 8 0 n T - 1 0 1 2 d H (t )/ d t, n T /s ec - 0 .1 0 - 0 .0 5 0 .0 0 S (t ), T E C U S S C n = 4 0 6 0 8 : 0 0 0 8 : 3 0 0 9 : 0 0 0 9 : 3 0 1 0 : 0 0 T i m e , U T - 0 .0 4 0 0 .0 4 0 .0 8 S (t ), T E C U n = 6 3 2 - 4 0 0 4 0 8 0 H (t ), n T S ( t ) I r k u t s k June 8, 2000 G O E S 1 0 ( 1 3 5 0 W ) a b c d e a b c d e n = 403 n = 97 Irkutsk 51 GPS detection of the instantaneous response of the global ionosphere to strong magnetic storms with sudden commencement dH/dt; τ - delay between a minimum of the S(t) series and the SSC time; and τ is the correlation coefficient between S(t) and dH/dt. The distinguishing features of the storms of April 6; June 8, July 15, 2000 and April 4, 2001 are the relatively high maximum values of dH/ dt at the SSC time - 0.2, 0.4 nT/s, respectively. For these storms, high values of the correlation coefficient of S(t) and dH/dt, – 1 and – 0.9, were obtained. It is evident from the table that the global negative TEC disturbance is quite well correlated with the derivative of the magnetic field strength (the correlation coefficient not below – 0.8) but is delayed with respect to these variations by about 3-9.5 min. 6. Discussion and conclusions Using the new technology for global GPS detection of ionospheric disturbances, GLOBDET, it has been ascertained that a drastic increase in the time derivative of the magnetic field strength at the time of magnetic storms is accompanied by an almost simultaneous Table II. Statistics of experiments. N Date DST, nT Time of max dH/dt, ∆T, UT m | A∑|, σ, τ, r SSC, UT nT/s TECU min min 1 6.01.1998 102.25 14:16 0.1 14 - 15 395 0.1 −6 −6.5 −0.8 2 23.04.1998 148.84 18:25 0.1 18 - 18:50 401 0.1 −3 −4 −0.8 0.4 I - 0 = 14.5 3 6.04.2000 144.82 16:39 0.2 16:30 - 18 472 2.6% −6 −7 −1 0.1 I - 0 = 35.3 4 8.06.2000 114 9:10 0.2 8:55 - 9:52 406 0.3% −3.5 −4.5 −0.9 0.2 I - 0 = 35.3 5 13.07.2000 88.32 9:42 0.1 9:19 - 10:32 227 0.3% −3 −13.5 −0.9 0.1 I - 0 = 37.3 6 14.07.2000 101.45 15:32 0.1 15 - 16 350 0.4% −4.5 −6 −0.8 0.1 I - 0 = 35.3 7 15.07.2000 166.5 14:37 0.4 14 - 15:21 203 0.3% −9.5 −10 −0.9 0.2 I - 0 = 25.9 8 31.03.2001 143 00:52 0.15 00:22 - 01:19 103 0.8% −4.5 −7.5 −0.9 0.1 I - 0 = 39.6 9 04.04.2001 132 14:55 0.15 14:32 - 15:17 91 0.2% −6 −3.5 −0.9 0.3 I - 0 = 35.3 10 11.04.2001 30.6 15:47 0.2 15:30 - 17:30 403 0.8% −9.0 −51 −0.9 52 Edward L. Afraimovich, Eugene A. Kosogorov, Ludmila A. Leonovich, Oleg S. Lesyuta and Igor I. Ushakov decrease of the mid-latitude TEC throughout the dayside. The corresponding correlation coefficient is not below – 0.8; the delay with respect to the time of the sudden magnetic storm commencement is about 3-10 min. This is most pronounced for magnetic storms with a clearly-defined SSC. The amplitude of the «instantaneous» TEC response for the events under consideration was 0.1-0.4 TECU, and the travel velocity of the disturbance from the dayside to the nightside was in the order of 10-20 km/s. These results are in reasonably good agreement with data obtained by recording the ionospheric response to a sudden com- mencement of strong magnetic storms, based on measurements of the frequency Doppler shift along a HF path (Huang et al., 1973). Thus the disturbances which we have analyzed do not fit in the model of Acoustic- Gravity Waves (AGW), and the appropriate physical mechanisms should be pursued when modeling the electromagnetic set of phenomena accompanying a strong geomagnetic disturbance. Effects caused by electric fields in the iono- sphere influence the behavior of its parameters and manifest themselves differently with a different latitude, which is due primarily to the geomagnetic field geometry. Under quiet conditions, the electric field has no marked effect on the behavior of ionospheric parameters as a function of time. However, during magnetic disturbances when the electric field is enhanced significantly (Kozelova et al., 2001), large variations in the height of the F2-layer maximum are observed, which can be caused by the vertical component of the E→ × Β→ - drift of ionospheric plasma (Rüster, 1965). According to Bryunelli and Namgaladze (1988), ionospheric effects of magnetospheric electric fields in mid-latitudes are caused by the divergence of the vertical ion flux produced by the horizontal electric field. The electric field-induced transport during a disturbance in the F2-region with the E → × Β→ - drift velocity is effective both on day and night. The electric field can reach 10-15 mV/m in mid- latitudes (Ogawa et al., 1975); the unperturbed values of the field are smaller by a factor of 2-3. The influence of the vertical component of the E→ × Β→ - drift on the height distribution of electron density was investigated by different methods (Tanaka and Irao, 1973). As the E→ × Β→ - drift velocity is virtually independent of the F 2 - layer height, the profile of the electron density distribution is displaced, under the action of the electric field, as a single whole without changing its form substantially (Bryunelli and Namgaladze, 1988). The eastward-directed electric field causes the height of the F2-layer maximum to rise, while the westward-directed field is responsible for its lowering. The h m F 2 variations are quite well correlated with zonal electric field variations, with a small delay, in the order of 20 min, at night and with a still smaller delay during the daytime. N m F 2 variations at the time of an instantaneous change of the field are poorly pronounced and show up as a decrease in critical frequencies by about 0.1-0.3 MHz. As suggested in Park (1974), such effects corresponded to electric fields with an amplitude of 5-10 mV/m. Unambiguous interpretation of the iono- spheric effects in terms of the action of electric fields in earlier studies was made difficult by the fact that similar effects can also be caused by AGW generated by auroral electrojets during substorms. The deciding criterion for separating the effects of electric fields and AGW is the delay between geomagnetic and ionospheric effects. In the former case, where the influence of electric fields is predominant, the delay is a few minutes (Bryunelli and Namgaladze, 1988; Volkov and Namgaladze, 2001). In the latter case, where ionospheric effects are triggered by AGW, the delay averages 45-60 min (Shashunkina, 1968, 1972; Afraimovich et al., 2000b). The analysis made in this paper which describes the occurrence of ionospheric effects on abrupt changes of the electric field, and also the data on TEC variations together with dH/dt obtained in this study, point to the penetration of electric fields of magnetospheric origin with typical nonstationarity times of 30-60 min, into the mid-latitude ionosphere during periods of drastic magnetic disturbances. Acknowledgements The author is grateful to N.N. Klimov and E.A. Ponomarev for their encouraging interest in this study and active participation in discus- 53 GPS detection of the instantaneous response of the global ionosphere to strong magnetic storms with sudden commencement sions. We are grateful to S.A. Nechayev and V.V. Kharchenko for providing data from magnetic observatory Irkutsk. Thanks are also due to V.G. Mikhalkovsky for his assistance in preparing the English version of the T E X manuscript. We greatly appreciate the Referee’s efforts to improve the submitted manuscript. This work was done with support from both the Russian foundation for Basic Research 00-05-72026 and RFBR grant of leading scientific schools of the Russian Federation 00-15-98509. REFERENCES AFRAIMOVICH, E.L. (2000): The GPS global detection of the ionospheric response to solar flares, Radio Sci., 35, 1417-1424. AFRAIMOVICH, E.L., E.A. KOSOGOROV, L.A. LEONOVICH, O.S. L E S Y U TA and I.I. U S H A KOV (2000a): Instantaneous Response of the Ionosphere to a Sudden Commencement of the Strong Magnetic Storm of April 6, 2000, http://xxx.lanl.gov/abs/physics/0007006. AFRAIMOVICH, E.L., E.A. KOSOGOROV, L.A. LEONOVICH, K.S. PALAMARCHOUK, N.P. PEREVALOVA and O.M. PIROG (2000b): Determining parameters of large-scale traveling ionospheric disturbances of auroral origin using GPS-arrays, J. Atmos. Sol.-Terr. Phys., 62, 553-565. BOWMAN, G.G. (1965): Travelling disturbances associated with ionospheric storms, J. Atmos. Terr. Phys., 27, 1247- 1261. BUONSANTO, M.J. (1995): A case study of the ionospheric storm dusk effect, J. Geophys. Res., 100, 23,857- 23,869. BRYUNELLI, B.E. and A.A. NAMGALADZE (1988): The Physics of the Ionosphere (Moscow: Nauka), pp. 528 (in Russian). GURTNER, W. (1993): RINEX: the Receiver Independent Exchange Format Version 2, http://igscb.jpl.nasa.gov/ igscb/data/format/rinex2.txt. HO, C.M., B.A. IIJIMA, X.P. LINDQWISTER, A.J. MANNUCCI, L. SPARKS, M.J. REYES and B.D. WILSON (1998): Ionospheric total electron content perturbations monitored by the GPS global network during two Northern Hemisphere winter storms, J. Geophys. Res., 103, 26,409-26,420. HOFMANN-WELLENHOF, B., H. LICHTENEGGER and J. COLLINS (1992): Global Positioning System: Theory and Practice (Springer-Verlag, Wien, New York), p. 327. H UA N G , Y.-N. and K. C H E N G (1991): Ionospheric disturbances at the equatorial anomaly crest region during the March 1989 magnetic storm, J. Geophys. Res., 96, 13,953-13,965. HUANG, Y.-N., K. NAJITA and P. YUEN (1973): The ionospheric effects of geomagnetic sudden commencements as measured with an HF Doppler sounder at Hawaii, J. Atmos. Terr. Phys., 35, 173-183. KLOBUCHAR, J.A. (1997): Real-time ionospheric science: the new reality, Radio Sci., 32, 1943-1952. KOZELOVA, T.V., B.V. KOZELOV and L.L. LAZUTIN (2001): Electric fields during the substorm according to data from the CRRES satellite, Geomagn. Aeron., 41, 39-49 (in Russian). MAKAROV, G.A. (1994): Sudden geomagnetic storm commencements and conditions in the solar wind, Geomagn. Aeron., 34, 1-6 (in Russian). MANNUCCI, A.J., C.M. HO, U.J. LINDQWISTER, T.F. RUNGE, B.D. WILSON and D.N. YUAN (1998): A global mapping technique for GPS-driven ionospheric TEC measurements, Radio Sci., 33, 565-582. MULLAYAROV, V.A. and E.O. MUZLOV (2000): SC effects in VLF emission, Geomagn. Aeron., 40, 102-106 (in Russian). NIELSEN, N. and F. HONARY (2000): Observations of ionospheric flows and particle precipitation following a sudden commencement, Ann. Geophys., 18, 908- 917. OGAWA, T., Y. TANAKA, A. HUZITA and M. YASUHARA (1975): Horizontal electric fields in middle latitude, Planet. Space Sci., 23, 825-830. PARK, C.G. (1974): A morphological study of substorm- associated disturbances in the ionosphere, J. Geophys. Res., 79, 2821-2827. POTAPOVA, N.I. and D.B. SHAPIRO (1961): Height distribution of ion density during the June 1959 ionospheric distur-bances, Ionos. Issled., 6, 23-28 (in Russian). RÜSTER, R. (1965): Height variations of the F2-layer above Tsumeb during geomagnetic bay-disturbances, J. Atmos. Terr. Phys., 27, 1229-1245. SHASHUNKINA, V.M. (1968). Ionospheric effect of a sudden magnetic storm commencement in the years of maximum and minimum solar activity, Geomagn. Aeron., 8, 184-187 (in Russian). SHASHUNKINA, V.M. (1972): Research results on the sudden magnetic storm commencement effect, Ionos. Issled., 20, 154-165 (in Russian). TANAKA, T. and K. HIRAO (1973): Effects of an electric field on the dynamical behavior of the ionospheres and application to the storm time disturbance of the F-layer, J. Atmos. Terr. Phys., 35, 1443-1452. VOLKOV, M.A. and A.A. NAMGALADZE (2001): Influence of the magnetospheric convection electric field on the formation of field-aligned currents of the substorm break-up phase, Geomagn. Aeron., 41, 33-38 (in Russian).