629_647 Jodogne 35_1.pdf ANNALS OF GEOPHYSICS, VOL. 45, N. 5, October 2002 629 Ionosphere-plasmasphere response to geomagnetic storms studied with the RMI-Dourbes comprehensive database Jean-Claude Jodogne and Stanimir M. Stankov Royal Meteorological Institute of Belgium, Brussels, Belgium Abstract Presented is a review of the ionospheric storm research carried out at the Royal Meteorological Institute (RMI) of Belgium. The studies are based on the opportunities offered by the long-term measurements and the accumulated comprehensive database managed by the RMI Geophysics Centre at Dourbes (4.6°E, 50.1°N). Reported are case studies of major storms, and also results in developing forecasting and density reconstruction methods. 1. Introduction he ionospheric storm, defined as the set of large perturbations in the geomagnetic field and the ionosphere-plasmasphere system charac- teristics induced by solar flares and coronal mass ejections, is still a challenging subject of in- vestigation. Observations of the ionosphere- plasmasphere system response to storms have been well documented for several decades. Nevertheless, many questions remain unan- swered and many early studies are now being questioned and improved by new discoveries. Recently, particular importance has been attri- buited to studies of «space weather» (all geospace Mailing address: Dr. Stanimir M. Stankov, Royal Meteorological Institute, Av. Circulaire 3, B-1180 Brussels, Belgium; e-mail: sstankov@oma.be Key words ionospheric storms – synthetic geomag- netic index – electron profile reconstruction phenomenae caused by the radiation and atomic particles emitted by the Sun and other stars). Traditionally, the geomagnetic storm effects have been investigated by means of vertical radio-sounding measurement techniques (ionosondes) providing long series of data for the ionised atmospheric layers. Later, with advances in satellite technology, in situ ob- servations of the plasma composition and behaviour have helped the research. Most recently, measurements of the Total Electron Content (TEC), both ground- and space-based, contributed significantly to the knowledge of the storms’ global manifestation. In addition to the above methods, numerical modelling approaches allow the combination of various observations, better coupling with theoretical investigations, and spreading the information on essential processes and parameters over areas which are not directly observed. The purpose of this paper is to present the ionospheric storm research carried out at the Royal Meteorological Institute (RMI) of Bel- gium. The studies are based on the opportunities offered by the long-term measurements and the accumulated comprehensive database managed by the RMI Geophysics Centre at Dourbes 630 Jean-Claude Jodogne and Stanimir M. Stankov (4.6°E, 50.1°N). Reported are case studies of two major storms in view of some important topics of interest, and preliminary results in developing new forecasting and electron density recon- struction methods. 2. RMI Dourbes Geophysics Centre: history and observations The RMI Geophysics Centre at Dourbes is a complex observational site (fig. 1) consisting of several observatories – ionosphere sounding, atmospheric, geomagnetic, cosmic rays, TEC, etc. – all connected with optical-fibre commu- nication lines. A fast link allows commu- nication with the RMI in Brussels, the Royal Observatory of Belgium (ROB) and the Bel- gian Institute for Space Aeronomy (ISA). Some information is shared and available via the world wide web (WWW). First measurements at Dourbes started more than half a century ago and a vast amount of data have been accu- mulated during the years. 2.1. Vertical ionospheric sounding The ionospheric soundings are carried out by a digital sounder (DGS-256) developed by the University of Massachusetts Lowell and fea- turing the operational characteristics given in table I. The DGS-256 ionosonde automatically calculates the electron density profile assuming a model valey and outputs the coefficients of a Chebyshev polynomials expansion. It is very important because it provides an opportunity for comparison of the ionosonde TEC (computed using ionosonde measurements) with the GPS TEC (GPS-based observations of TEC). 2.2. GPS-based total electron content Some of the major advantages of using GPS to derive TEC are: the use of high frequencies (so that ionospheric absorption and geomagnetic field effects on the radio signal are small), elec- tron content evaluated up to a very high altitude (20 000 km, providing the plasmasphere con- tribution to the content), operational global coverage (e.g., using the International GPS geodetic network), cost efficient observations (continuous signals and commercially available receivers). At the RMI-Dourbes, a GPS receiver is col- located with the digital ionospheric sounder. Using the GPS signals on two coherent carrier frequencies (L1/L2 = 1575.42/1227.6 MHz), the TEC computation procedure (Warnant, 1997, 1998) is based on the «geometry-free» combi- nations of GPS code (P i p,GF ) and phase (F i p,GF ) measurements (2.1) where P i p is the code measurement made by receiver p on i – th satellite, F i p is the carrier phase measurement made by receiver p on the i – th P P P F F f f Fp GFi p Li p Li p GFi p Li L L p Li, , , , , ,/= − = − ( )1 2 1 1 2 2 , Fig. 1. The measurement and communication network of the Royal Meteorological Institute of Belgium. SAFIR – Thunderstorm alert system using radio interferometry. A 2p connection consists of two cables with two optic fibres in each cable. 631 Ionosphere-plasmasphere response to geomagnetic storms studied with the RMI-Dourbes comprehensive database satellite, and f L1 , f L2 – the frequencies on the L1, L2 carriers respectively. Rewritten as functions of TEC, the above equations read (2.2) (2.3) where N i p is the phase ambiguity, TEC i p is the slant electron content (along the i – th satellite raypath) in TECU, Di, D p are the i – th satellite and p− th re- ceiver differential group delays. The ambiguity is eliminated by the following combination of «geometry-free» code and phase measurements: (2.4) where λ L1 is the L 1 carrier wavelength. The formula requires the estimation of the receiver and satellite group delays, which are obtained from (2.2) after modelling TEC by means of a simple polynomial depending on latitude and local time. The conversion to vertical TEC is performed by assuming that the ionosphere is a layer of infinitesimal thickness located at a «mean ionospheric» height of 350 km and using a simple cosine function of the zenith angle at the «ionospheric point» (the raypath’s point at the mean ionospheric height). Finally, the TEC value is calculated from (2.3). To obtain a TEC value, representative of the ionosphere above a given observing station, the following is applied: first, selected are all TEC values within a lati- tude difference of 1.5° from the latitude of the observing station, and second, computed is the mean of these TEC values over 15 min periods. Measurements have been conducted in Dourbes on a regular basis since 1994 and a large TEC database has been created for the best part of the current solar-activity cycle. P D D p GF i p i p i , .= − × + −( )−1 05 10 17 TEC P F D D Np GFi L p GFi p i L p GFi, , ,− = −( ) −λ λ1 1 F Np GFi pi p GFi, ,.= − × +−5 52 10 17 TEC Operational parameter Range DB049 ionogram settings Frequency range 0.5 - 30 [MHz] 1-16 Frequency scale Linear or logarithmic Linear Frequency steps 5, 10, 25, 50, 100, 200 [kHz] 100 [kHz] Range resolution 2.5, 5.0, 10.0 [km] 5.0 [km] Number of range pixels 128, 256 128 Range start Selectable 60 [km] Pulse repetition range 50, 100, 200 [s– 1] 100 [s– 1] Pulse width 66, 133 [µs] 66 [µs] Transmitter RF power ≤ 10 [kW] 5 [kW] Duration of ionogram Selectable 4 [min] Phase code Interpulse pseudo-random biphase coding Digitization 12 bit linear Amplitude resolution 0.25 [dB] Phase resolution 1.4° Doppler resolution 4 [Hz] or 12 [Hz] 4 [Hz] Wave polarization O/X Data display Video and paper Data storage Disk, cartridge Table I. The Dourbes (DB049) digital ionosonde – operational parameters and settings (after Reinisch, 1996). 632 Jean-Claude Jodogne and Stanimir M. Stankov 2.3. Geomagnetic measurements From the beginning, the goal of the designers was the best measurements of the geomagnetic induction. The observatory building was erected without any ferromagnetic material even for the external paint. Double-core fluxgate magnet- ometers are used for measurements (Primdahl, 1979; Jankowski and Sucksdorff, 1996). These magnetometers are robust, have reliable elec- tronics, and are capable of recording not only the field variations but its full strength as well. Absolute declination, inclination and induction values are manually determined with Theodolite, Diflux and Overhauser protons magnetometer respectively. An optical pumping Cesium mag- netometer provides data each second. The precision level is 0.1 nT for the induction and 0.001° for the angles. Data are sent each hour to the METEOSAT satellite in the INTERMAGNET format. The large data bank on magnetic support started in 1957. A team uses also cryogenic mag- netometers to make paleo- and archeo-magnetic studies. The horizontal component is deduced from measurements of the total field (B) and inclination (I ), and their trigonometric relation H = B.cosI. Local magnetic observations can be very helpful in many aspects. One particular advan- tage, concerning the present study, is the op- portunity for near real-time detection of geo- magnetic storm onset. Such detection can be performed by identifying the SSC using obser- vations of the geomagnetic field’s horizontal component variations. The apparent shape of SSC(H) is not a simple increase but rather mani- fests some variety. Morphological SSC studies (Matsushita, 1967) reveal that the shape of SSC(H) depends on local time and latitude, and can be sorted in the following major types: i) a small negative impulse preceding the main positive impulse; ii) a main positive impulse alone; iii) an increase (lasting up to about six minutes) followed by a decrease to a level lower than the initial pre-SSC level. However, the SSC detection is not always easy, particularly when performed on a real-time basis. A similar abrupt disturbance without succeeding large changes, called «sudden impulse», can be mistaken for a SSC. 2.4. Cosmic rays The interplanetary space is continuously traversed by very energetic nuclear particles called cosmic rays. Their origin is not definitely determined − some certainly originate from the Sun and other stars, but others are associated with various energetic processes in the galaxy. Cosmic rays consist mostly of protons (~ 94%) and helium nuclei (~ 5.5%), with few nuclei of heavier elements (galactic component of the cosmic rays). The particles have a broad energy range, from less than 1 GeV (109eV) to more than 1011 GeV (1020eV). In interplanetary space, the cosmic ray flux is in the order of 0.6 particles cm– 2s– 1, the energy flux is 5 GeV cm– 2s– 1, and the space density – approximately 0.8 × 10– 10 particles cm– 3 (Akasofu and Chapman, 1972). The high energy of the galactic cosmic rays means that they also have a very high penetrating power. Near the Earth, their motion is influenced by the geomagnetic field. When penetrating deep into the atmosphere, the cosmic rays collide with the atmospheric constituents and the results are different according to the energy input. Protons with energies of about 0.1 GeV isotropically scatter the products of the collisions, such as protons, α-particles and neutrons. The paths of these products appear as radial lines from the collision point, forming a «star». If the energy of the protons is of higher order (several GeV), then π mesons are also produced, which are scattered forward and their paths appear as a «shower». When the energy of the incident protons is more than 103 GeV, the products are scattered into a narrow cone called a «jet». The protons produced collide with new atmospheric particles generating similar reactions. Thus, the newly-generated rays/particles permeate the atmosphere and can even reach ground level. Eleven-year cycle variation of the cosmic rays is detected, affecting the lower-energy cosmic rays much more than the higher-energy ones. The interplanetary shocks and their drivers (the fast coronal mass ejections) play an im- portant role in the modulation of the cosmic rays. Since the earliest observations of the cosmic rays (Forbush, 1938) it has been recognized that transient decreases in the cosmic ray intensity (Forbush decreases) occur at the times of sudden 633 Ionosphere-plasmasphere response to geomagnetic storms studied with the RMI-Dourbes comprehensive database commencement geomagnetic storms and that the cosmic ray intensity varies regularly and inver- sely to the sunspot cycle, with the ray flux most pronouncedly decreased at high solar activity. Therefore, the Forbush decreases (produced by large-scale solar-wind disturbances) can be successfully used for investigating ionospheric storms and related phenomena (Barnden, 1973; Barouch and Burlaga, 1975; Zhang and Burlaga, 1988; Cane et al., 1993, 1994; Shadrina et al., 2002; Gonzalez et al., 2002). At the RMI-Dourbes, a 9NM-64 neutrons monitor provides measurements of the secondary neutron component of the cosmic rays on the ground. After pressure correction this component faithfully follows the primary cosmic rays inten- sity. A large data bank of the hourly intensities has been existing since 1967. 3. Case studies of storms In a broad sense, the «ionospheric storm» is defined as the set of disturbances that appear in the ionosphere immediately following a geo- magnetic storm commencement (Hargreaves, 1992). The mechanisms accounting for the storm effects are composition changes, neutral winds, electric fields, etc. The individual case studies of geomagnetic storms and effects are considered as powerful instruments deepening our under- standing of the inter-related processes of energy and momentum transfer, and of the kinetic, chemical, and electrodynamic coupling in the chain of solar-terrestrial relations. The main reason for this is the fact that during ionospheric storms the ionospheric processes are intensified and a phenomenon normally occurring over a long period is observed during a relatively short time period. Two major topics were investigated: the neutral winds effect on the ionosphere and the electrodynamic mechanism. The global manifestation of the neutral wind effect starts at high latitudes where the per- turbations in the interplanetary medium and the magnetosphere yield increased energetic particle precipitation, convection electric fields and associated Joule heating. The Joule heating raises the temperature of the thermosphere at auroral latitudes resulting in enhanced equatorward winds and travelling ionospheric disturbances in all local time sectors. On the night-side the pole- to-equator wind field is enhanced. On the day- side, the equatorwards wind counteracts (and often overpowers) the solar-driven poleward winds. As a result of the air drag, the day-time F-region experiences upward motion and an ion density increase. This contributes to the «positive phase» in the critical frequency and total electron content values, whose increase is usually ob- served on the first day of the storm. Large-scale structures of the Earth’s plasma environment are determined by electrodynamic processes of magnetospheric and solar wind origin. For example, the electrodynamic drift can redistribute plasma to regions of either high or low chemical loss. Moreover, the E × B effect of the electric fields mapped down the geomagnetic field lines during storms is instantaneous. Unan- swered questions are: the exact interplay of processes during the positive phase, the degree to which magnetospheric convection effect can penetrate to midlatitudes, coupling of atmos- pheric regions and processes, etc. (Szuszczewicz et al., 1998; Namgaladze et al., 2000; Foerster and Jakowski, 2000). Geomagnetic storms are predominantly in- vestigated using magnetic, ionosonde, and GPS- based measurements, all available at the Dourbes Geophysics Centre. The geomagnetic activity characterisation is based on evaluation of the geomagnetic disturbances with two major types of indices in use – the K and Dst indices (Men- vielle and Berthelier, 1991). The ionosonde measurements are important not only with the critical frequency ( f 0 F 2 ) measurements but also with the peak density (N m F 2 ), peak height (h m F 2 ), and propagation factor M3000F2: a) The f0F2 values respond quickly to the storm, and together with TEC measurements, they can provide a crucial information on the slab thick-ness. b) The N m F2 values will be used for calculat- ing the slab thickness τ, where τ [km] = 806.405 TEC [1016m– 2]/f 0 F 2 2 [MHz]. This parameter can provide information on the plasma uplifting processes, plasmasphere fluxes, etc. c) The h m F2 will give more information on the profile and the F 2 layer response to the storm, particularly in connection to the uplifting process. 634 Jean-Claude Jodogne and Stanimir M. Stankov The TEC measurements are used together with the ionosonde measurements. No storm studies seem possible nowadays without the TEC value, especially after the introduction of the GPS system providing opportunities for a global cov- erage. 3.1. The storm on 10th January 1997 As a first example, we present the study of the storm on 10th January 1997, which has been well documented on a global scale because of the pre-observed and announced large coronal mass ejections on 6th January 1997 (Jakowski et al., 1999; Mikhailov and Foerster, 1999). This storm is characterised by a pronounced positive response from both TEC and f 0 F 2 parameters – for example, the noon TEC values are more than 100% higher than the monthly median (fig. 2c). The basic storm features are summarised in ta- ble II. Several aspects of the ionospheric response can be investigated from Dourbes and it will be demonstrated on this particular storm. The positive storm phase is generally ac- cepted to be induced by the strong southward thermospheric winds. The most significant consequence is the plasma uplifting effect – the strong meridional winds push the F2-layer plasma upward, leading to reduced F2-layer peak density and increased h m F 2 and TEC. This effect can be easily detected from measurements through the slab thickness shape parameter τ; increasing slab thickness (e.g., rising TEC accompanied by constant or decreasing f0F2) indicates plasma uplifting. A clear demonstration of this process is detected in the 09:00-10:00 UT period when TEC is rising, f0F2 is almost unchanged, and hmF2 also increases (fig. 2c,d). The wind-induced uplifting takes place along the geomagnetic field lines, and is therefore most efficient near the 45° geomagnetic inclination. Situated at 50.1°N, the Dourbes Geophysics Centre provides an excel- lent database for investigating the above-men- tioned effect. Other storm effects, the plasma fluxes from the plasmasphere, are due to the compression of the geomagnetic field during the onset phase of the storm (Foerster and Jakowski, 2000) and/or inter-hemispheric flows (Menietti et al., 1988). Such plasmaspheric fluxes can be easily detected: e.g., a high slab thickness value during night- time is certainly an indication of an enhanced plasma influx from above, and this is exactly the situation observed near 02:00UT (fig. 2d). The next effect of interest is associated with the generated convection electric field. The sharp Fig. 2a-d. Ionosphere-plasmasphere response to the 10th January 1997 geomagnetic storm, as observed at the site of the Dourbes ionosonde station DB049 (4.6°E, 50.1°N). a b c d 635 Ionosphere-plasmasphere response to geomagnetic storms studied with the RMI-Dourbes comprehensive database decrease of the slab thickness after 02:00 UT shows the reaction to this field. There is no unique explanation for the decrease; a possible expla- nation is the occurrence of zonal disturbance electric fields. A westward electric field moves the plasma from higher to lower L-shells with smaller volumes. Thus, the plasma pressure in the flux tubes is increased, causing an enhanced plasma flux from the plasmasphere down to F2 region altitudes. An eastward electric field results in the opposite effect of plasma decompres- sion so that the plasmaspheric flux decreases (Mikhailov and Foerster, 1999). The ionospheric storms can be much better investigated on a global scale. For the purpose, the RMI and ROB jointly upgrade the existing GPS TEC measurements to provide the TEC values at practically every point in the European region. The propagation velocity of the per- turbation is one of the important storm charac- teristics that can be investigated in this way. Some authors (Hernandez and Roble, 1978; Proelss et al., 1991; Proelss, 1997; Jakowski et al., 1999; Foerster and Jakowski, 2000) estimate this velocity in the order of 400-900 ms– 1. Theo- retically, due to the reduced ion drag on the night side, the propagation should be much faster than on the day side. It is possible to estimate the time delayed ionospheric reaction using simultaneous GPS TEC and ionosonde data. For example, if we observe the exact time of the day-time TEC peak (fig. 2c) at various geomagnetic latitudes near 0°E, we find that the peak is observed at 65°N in 11:30 UT, at 50°N (Dourbes) in 12:30 UT, at 35°N in 13:00 UT. Thus, it obviously takes about 90 min for the plasma uplifting to move from 65°N to 35°N, which corresponds to a wind velocity of approximately 600 ms– 1. It is of crucial importance to have higher resolution data (more frequent measurements) during storm times. For example, if we work with hourly data, the above peak would be observed in 11:00, 12:00, 13:00 UT, hence the time delay would be 120 min resulting in a different wind velocity estimate. At the Dourbes Geophysics Centre we are able to provide 15 min sampling rate for GPS TEC and even better rate for the digital iono- sonde. Thus, the TEC peak is observed at 12:30 UT which ensures better estimation of the propagation velocity. 3.2. The storm on 15th July 2000 Another example is the geomagnetic storm observed in the period 15th-17th July 2000. The SSC is recorded at 14:35 UT and follows a disturbed period starting approximately 24 h earlier. It is a very intense storm (fig. 3a,b) showing a pronounced negative response in both the TEC Table II. Summary of the 10th January 1997 storm characteristics. Feature Value Comments Storm onset (SSC) 01:00 UT 01:00 UT, 10/1/1997 Dst (min) – 78 nT 09:00 UT, 10/1/1997 Total storm period 10 h 08:00 UT (10/1/1997) – 17:00 UT (10/1/1997) (Dst ≤ – 50 nT) Initial plus main phase 8 h 01:00 UT (10/1/1997) – 09:00 UT (10/1/1997) (SSC → Dst (min)) Main phase 1 h 08:00 UT (10/1/1997) – 09:00 UT (10/1/1997) ( – 50 nT → Dst (min)) Recovery phase 9 h 09:00 UT (10/1/1997) – 18:00 UT (10/1/1997) (Dst (min) → – 50 nT) 636 Jean-Claude Jodogne and Stanimir M. Stankov and f 0 F 2 values. Around 03:00 UT on 16th July 2000, the relative differences increase to more than – 80% for TEC and – 50% for the critical frequency (fig. 3d). The most important storm characteristics are given in table III. At high latitudes, the negative phase is attri- buted to the increased molecular concentrations of the thermosphere, enhancing the recom- bination processes and reducing the electron density. The standard explanation for a negative phase at middle and low latitudes (Szuszczewicz et al., 1998) is the enhanced equatorward wind system advecting the modified thermospheric composition to lower latitudes. This equatorward penetration is deeper at night-time, and theo- retically, the corotation, coriolis forces, and zonal winds transport the modified thermospheric composition also to the day-time hemisphere and maintain the negative storm phase. Elements which complicate this simple picture involve: penetrating electric fields of magnetospheric origin, dynamo-driven fields, vertical and hori- zontal interactions between the wind system, thermospheric composition, ionospheric plasma, storm intensity, storm onset time (and associated local times in different longitude sectors), variations in pre-storm conditions, evolution of the controlling interplanetary-magnetospheric inputs, etc. This storm is interesting in that it reveals another phenomenon, the so-called «ionogram G condition». Under G conditions, a measure- ment is influenced or impossible because the ionization density of the layer is too small to deduce the measurement value accurately. Such conditions occur during ionospheric disturbances associated with geomagnetic storms. Then, the electron density in the F 2 layer becomes equal to or lower than the electron density in the F 1 layer. In the result, f0F2 is expressed by the value of f0F1 with the descriptive letter G. This is exactly the case on 16th July 2002 (fig. 4a): f 0 F 2 is equal to f 0 E during the day, except for a few hours (15:00-17:00 LT) in the afternoon. However, it is important to separate the G condition from two other confusing cases. The first case is when the F 2 trace disappears due to either interference or attenuation. The other case is when the trace cannot be recorded because of its extremely high virtual height exceeding the height range limit (descriptive letter W ). This is obviously what happens from 21:00 UT on 15th July 2000 (h’F2 = 520 km) to 03:00 UT on 16th July 2000 (h’F2 = 513 km) when h’F2 is above these values and above the height range limit (fig. 4b). In all of the above cases, the automatic ionogram scaling may fail to properly identify the f0F2 and hmF2 values, may not produce elec- tron density profiles, or deliver profiles which can be quite different from the real ones. The Fig. 3a-d. Ionosphere-plasmasphere response to the 15th-17th July 2000 geomagnetic storm as observed at the site of the Dourbes ionosonde station DB049 (4.6°E, 50.1°N). a b c d 637 Ionosphere-plasmasphere response to geomagnetic storms studied with the RMI-Dourbes comprehensive database associated ITEC (ionosonde-based TEC) calcu- lations, derived from such erratic profiles, can be very misleading. Concurrent GPS-based TEC measurements can be extremely useful. On the one hand, TEC measurements provide oppor- tunities for identification and correction of ionosonde measurements under such G/W conditions (fig. 3c,d). On the other hand, the TEC values can be used to reconstruct the electron profiles using the reconstruction method described in the following sections. 4. Synthetic index of local response to the geomagnetic activity The response of a given ionosphere charac- teristic F (like f 0 F 2 , M 3000 F 2 , or TEC) is better analysed through the relative deviation, F rel = = (Fmeas – Fmed)/Fmed, of the actual observed char- acteristic from its median value, Fmed. The re- sponse strongly depends, in a highly non-linear fashion, on geomagnetic indices such as the planetary indices A p and K p (Menvielle and Ber- thelier, 1991). For developing geomagnetically- correlated short-term forecasting methods, a new synthetic index is required (Stankov et al., 2001; Muhtarov et al., 2002; and the references there- in). This new synthetic index is actually a func- tion showing the Frel dependance on the geo- magnetic index (say K p ) used for the forecast. However, the new synthetic index (function) Fig. 4a,b. Ionosonde measurements during the 15th - 17th July 2000 geomagnetic storm recorded at the Dourbes ionosonde station DB049 (4.6°E, 50.1°N). Table III. Summary of the 15th July - 17th July 2000 storm characteristics. Feature Value Comments Storm onset (SSC) 14:35 UT 14:35 UT, 15/7/2000 Dst (min) – 300 nT 21:00 UT, 15/7/2000 Total storm period 49 h 17:00 UT (15/7/2000) – 18:00 UT (17/7/2000) (Dst ≤ – 50 nT) Initial plus main phase 6.5 h 14:35 UT (15/7/2000) – 18:00 UT (15/7/2000) ( SSC → Dst (min)) Main phase 4 h 17:00 UT (15/7/2000) – 21:00 UT (15/7/2000) (– 50 nT → Dst (min)) Recovery phase 45 h 21:00 UT (15/7/2000) – 18:00 UT (17/7/2000) (Dst (min) → – 50 nT) a b 638 Jean-Claude Jodogne and Stanimir M. Stankov demonstrates strong seasonal and spatial varia- bility and should be determined for each meas- urement site separately. At Dourbes, there are two types of measurements that can be used for the purpose – ionosonde ( f 0 F 2 ) and GPS TEC measurements. 4.1. The f 0 F 2 response to geomagnetic activity The f0F2 response to geomagnetic activity was first deduced from f0F2 measurements carried out at Dourbes over the period from July 1957 to June 1989. It is calculated in the following manner. All values of the f0F2 relative deviations in a given month of the year are sorted according to the three-hour values of the planetary index Kp, which index is recorded in step values of 0.00, 0.33, 0.67, 1.00, 1.33,..., 8.67 and 9.00. Then, for each level of K p and each month, the basic statistics are calculated – average, standard deviation, data number, etc. Taking into account the expected seasonal variations in the iono- spheric response to geomagnetic storms (Fuller- Rowell et al., 2000), calculations were made and the results presented in fig. 5. The top four panels show the relative deviations of f 0 F 2 for the months of January, April, July, and September; the contour plot at the bottom shows the detailed annual behaviour of the averages. The first important observation was the distribution of the zero value of the f 0 F 2 relative deviation ( f0F2rel = 0), i.e. when f0F2meas = f0F2med. This is the value representing the «average» magnetically quiet conditions. It is supported by the fact that the majority of the observations cluster at and around the median. It becomes obvious from the bottom contour plot that Frel = 0 is observed for K p = 2.67 except for the winter months when F rel varies between the higher values of 3.00 and 4.00. Considering the defi- nition of a quiet magnetic activity (0 < K p < 4), it follows that even within this interval a trend exists in the observations of Frel: a generally positive response for K p < 2.67 and negative response for Kp > 2.67. The second important fact is that the re- sponse of f0F2 under intense geomagnetic activity (K p > 4) is generally negative, i.e. during storms the measurements are below their monthly medians. This negative response is particularly strong near the equinoxes (April and September) when the decrease reaches 25-30% during severe geomagnetic storms (K p > 7). On the other hand, the maximum positive values of F rel are also observed during April and September for small values of K p . 4.2. The GPS TEC response to geomagnetic activity It is interesting to compare the critical frequency response with the response of the GPS TEC. The GPS TEC time-series data available at Dourbes (July 1994 - October 2001) are not so long as the f0F2 series, and data scattering is rather large, but some important differences have already been detected. The TEC relative de- viations (TEC rel) have been sorted in the same way as the f0F2 deviations and the results are provided in fig. 6. Again, the top four panels show the TEC relative deviations for the months of January, April, July, September, while the contour plot details the entire annual behaviour. For the GPS-based TEC measurements, our calculations show that the TEC response is generally stronger and much more complex than the f0F2-based observations. First, the «quiet» behaviour, TEC rel = 0, is detected at lower K p values of around 2.0, which is a sign for a stron- ger response to the geomagnetic activity. Second, it is obvious that, for quiet magnetic conditions (K p < 4), the dependence of the relative TEC on K p is similar to that of f 0 F 2 . Third, for higher values of the K p index, the relative TEC stops decreasing and in many cases increases signi- ficantly particularly during winter months. At this stage of investigation, the nature of this positive response is attributed mainly to the prolonged positive phase in the TEC storm-time behaviour. It is also clear that in November and December the increase of TECrel during storms can be as much as 20-25%. There are indications that such an increase can be expected also during January and even February. Alternatively, pro- nounced decreases are observed in both April and September for values of K p > 4, although not as strong as in the f 0 F 2 case. 639 Ionosphere-plasmasphere response to geomagnetic storms studied with the RMI-Dourbes comprehensive database Fig. 5. The f 0 F 2 relative deviations ( f 0 F 2rel ) response to geomagnetic activity for the site of ionosonde station DB049 (4.6°E, 50.1°N) based on data from the 1957-1989 period. The average values of f 0 F 2rel (dashes) are approximated by 3−rd degree polynomials (solid lines). The vertical bars represent the standard deviations (doubled), and circles represent the data number (scaled on the right-hand axis). 640 Jean-Claude Jodogne and Stanimir M. Stankov Fig. 6. The GPS TEC relative deviations (TEC rel ) response to geomagnetic activity for the site of ionosonde station DB049 (4.6°E, 50.1°N) based on data from the 1994-2001 period. The average values of TEC rel (dashes) are approximated by 3−rd degree polynomials (solid lines). The vertical bars represent the standard deviations (doubled), and circles represent the data number (scaled on the right-hand axis). 641 Ionosphere-plasmasphere response to geomagnetic storms studied with the RMI-Dourbes comprehensive database However, more GPS TEC measurement data and further analysis are required to better evaluate the TEC response to stronger geomagnetic activity, including the local time variation of the storm effects (Proelss, 1993). Statistical studies show (Davies, 1990) that, especially in middle latitudes, a geomagnetic storm follows different patterns depending on the local time of the Sudden Storm Commencement (SSC). On the one hand, a geomagnetic storm in which the SSC occurs during daytime generally produces a TEC increase. On the other hand, the TEC during a storm in which the SSC occurs at night-time either goes directly to the negative stage, or increases for a short period during the afternoon of the following day. Nevertheless, the case GPS TEC Degree Month 0 1 2 3 January 4.8198 × 10– 2 – 2.2835 × 10– 2 – 0.3204 × 10– 2 0.1971 × 10– 2 February – 1.9767 × 10– 2 4.0903 × 10– 2 – 1.4147 × 10– 2 0.1828 × 10– 2 March 6.5163 × 10– 2 – 7.4277 × 10– 2 2.5464 × 10– 2 – 0.2299 × 10– 2 April 3.8189 × 10– 2 0.9663 × 10– 2 – 1.4048 × 10– 2 0.1399 × 10– 2 May 1.7888 × 10– 2 – 1.0004 × 10– 2 – 0.3160 × 10– 2 0.0603 × 10– 2 June – 3.4079 × 10– 2 8.8947 × 10– 2 – 3.7321 × 10– 2 0.3786 × 10– 2 July 7.1668 × 10– 2 – 7.1640 × 10– 2 2.6913 × 10– 2 – 0.2537 × 10– 2 August 1.3113 × 10– 2 – 1.9644 × 10– 2 1.1797 × 10– 2 – 0.2037 × 10– 2 September 4.8579 × 10– 2 – 4.0954 × 10– 2 1.3645 × 10– 2 – 0.1542 × 10– 2 October 7.6576 × 10– 2 – 1.5230 × 10– 1 6.0279 × 10– 2 – 0.5746 × 10– 2 November – 4.0179 × 10– 2 – 1.4980 × 10– 2 2.1028 × 10– 2 – 0.1998 × 10– 2 December – 0.8146 × 10– 2 2.1755 × 10– 2 – 0.2223 × 10– 2 0.0500 × 10– 2 f 0 F 2 Degree Month 0 1 2 3 January 2.4165 × 10– 2 – 2.5847 × 10– 2 1.2986 × 10– 2 – 1.8028 × 10– 3 February – 1.0463 × 10– 3 2.4421 × 10– 2 – 6.4418 × 10– 3 – 1.4525 × 10– 4 March 2.4913 × 10– 2 1.0363 × 10– 2 – 8.2630 × 10– 3 3.6672 × 10– 4 April 4.7779 × 10– 2 – 9.7067 × 10– 3 – 2.2601 × 10– 3 – 2.1101 × 10– 4 May 4.5731 × 10– 2 – 1.6177 × 10– 2 – 1.4128 × 10– 3 – 1.8441 × 10– 4 June 4.3341 × 10– 2 – 3.1992 × 10– 2 7.2534 × 10– 3 – 1.0285 × 10– 3 July 3.6014 × 10– 2 – 1.8669 × 10– 2 2.7797 × 10– 3 – 5.9729 × 10– 4 August 4.2705 × 10– 2 – 2.9214 × 10– 2 7.4882 × 10– 3 – 1.1622 × 10– 3 September 6.0539 × 10– 2 – 4.4463 × 10– 2 1.1422 × 10– 2 – 1.5466 × 10– 3 October 3.6506 × 10– 2 – 2.6692 × 10– 2 7.4757 × 10– 3 – 1.2045 × 10– 3 November – 1.1878 × 10– 2 3.6587 × 10– 2 – 1.0752 × 10– 2 4.7759 × 10– 4 December – 5.0906 × 10– 3 2.2692 × 10– 2 – 6.7083 × 10– 3 3.3266 × 10– 4 Table IV. Coefficients for the 3– rd degree polynomial approximation of the average monthly response to geomagnetic activity; based on f0F2 (top table) and GPS TEC (bottom table) measurements at the site of the Dourbes ionosonde station DB049 (4.6°E,50.1°N). 642 Jean-Claude Jodogne and Stanimir M. Stankov studies presented here reveal that there are storms which do not follow the common pattern proving the necessity of further investigations of such a complex issue. 4.3. Synthetic index – approximation formulae for the local ionospheric response Considering the above analysis, it is possible to deduce simple analytical expressions of the f0F2 and TEC average response to the activity on a local basis. Such expressions can be used as a new synthetic index replacing the traditional A p and Kp indices. The purpose of creating such index is the necessity to «linearize» the relation- ship between the relative deviations (Frel) of the ionospheric characteristic and the standard geomagnetic index (e.g., Kp), which is absolutely necessary for a reliable forecast based on the same characteristic. It is known that the mutual correlation between two random variables is highest when they are linearly dependent (Muhta- rov et al., 2002). New analytical expressions, using third- degree polynomials, have been deduced for the synthetic index of local ionospheric response to geomagnetic activity. The polynomial coef- ficients are given in table IV for the site of the Dourbes ionosonde station. 5. Reconstruction of electron density profiles Many ionosphere-plasmasphere processes, like the ionospheric storms and propagation delays in GNSS, can be estimated and corrected much more easily if the (vertical) electron density profile is available at a given location on a real- time basis. As the theoretical ionospheric models are not sufficiently good for the purpose, actual on-line measurements immediately available should be considered. 5.1. Electron profiles deduced from ionosonde measurements The DGS-256 ionosonde automatically calculates the electron density profile assuming a model valey and outputs the coefficients of a Chebyshev polynomials expansion (Huang and Reinisch, 1996, 2001). The deduced electron density profile is very important as it provides opportunity to compare the TEC computed using ionosonde measurements with the GPS TEC available at the same site. The Ionosonde Total Electron Content (ITEC) is computed in the following manner (Warnant and Jodogne, 1998). The bottomside vertical electron density profile is obtained after standardized inversion of the h’( f ) traces on vertical incidence ionograms with an accuracy of about 10 km as long as the profiles are monotonic. As the topside profile cannot be measured by the ground-based ionosonde, it is modelled with the β-Chapman function from the peak height up to the GPS satellite height of 20 200 km (5.1) where z = (h – h max )/H T and H T is the top-side scale height. The topside scale height is cal- culated as the mean of the scale heights computed at three points of the bottom-side profile located at 88%, 90%, and 92% of the Nmax. The ITEC is obtained after numerically integrating the electron profile. For the low solar activity period (1995-1996), the estimated mean and standard deviations of the (TEC-ITEC) distribution are 0.46 and 1.72 [1016m– 2], respectively. 5.2. Electron profiles deduced from ionosonde and GPS TEC measurements A new technique has been recently developed (Stankov and Muhtarov, 2001, Stankov et al., 2002a) for retrieval of electron density profiles from concurrent GPS TEC and ionosonde data. The ionosonde measurements are used primarily for obtaining the bottom-side profile. Ionosonde data, together with simultaneously-measured TEC, are required for the topside electron profile determination. The unique determination of the top-side profile also depends on the availability of the O+-H+ transition level (upper transition level, UTL) value. The following «recon- struction» formula is proposed for calculation of the top-side (h > h m F 2 ) electron density profile: N h N z z T ( ) = − − −( )( )maxexp exp1 643 Ionosphere-plasmasphere response to geomagnetic storms studied with the RMI-Dourbes comprehensive database (5.2) where H O+ and H H+ are the top-side O+ and H+ scale heights, N O+ (h m F 2 ) and N H+ (h m F 2 ) are the O+ and H+ densities at the F2-layer peak height (h m F2 ), and sech(h) = 1/cosh(h), cosh(h) = = 0.5(exp(h) + exp(– h)). Several assumptions have been made: first, the height of the O+ density maximum is equal to the height of the H+ density maximum; second, the neutral particles and ions have the same scale height; third, the hydrogen and oxygen ion temperatures are equal. There are four unknown variables in the proposed formula: the ion densities N O+(h m F 2 ) and NH+(hmF2), and the ion scale heights HO+ and HH+. However, the number of unknown variables can be decreased by considering that H H+ ≈ 16.H O+ following the definition of scale height. Therefore, the next system is constructed to determine the remaining three unknowns (5.3) (5.4) (5.5) The first equation represents the principle of plasma quasi-neutrality and the fact that at the peak height, the only ion densities that really count are the oxygen and hydrogen ion densities. The second equation denotes the fact that the hydrogen and oxygen ion densities are equal at the O+-H+ transition level (h tr ); the level is obtained from a model based on satellite and rocket measurements (Stankov, 2002). The third equation is obtained after integrating the proposed N e (h) «reconstruction» formula (5.2) from h m F 2 to infinity. The above system (5.3)- (5.5) is solved numerically. 5.3. Exemplary reconstruction results The reconstruction is demonstrated on actual GPS TEC and digital ionosonde data obtained at the RMI Geophysics Centre and using its meas- urements and communication network (fig. 1). The deduced electron density distribution is plotted in fig. 7 (top panel). The required input values of GPS TEC, f0F2, f0E, M3000F2, UTL and h m F 2 are given in the middle and bottom panels of fig. 7. The TEC and f 0 F 2 show strongly- correlated diurnal behaviour; both quantities increase sharply in the early morning, reach their absolute maximum just before noon, and then start gradually decreasing. Relatively high values are maintained throughout the after- noon, followed by a 10% increase in the early evening. After that, both TEC and f0F2 fall rapidly to their corresponding absolute minima observed at 02:00-03:00 LT. In the bottom panel, the F 2 - peak-density and O+-H+ ion transition level are also provided. The level, starting from 680 km at midnight, increases up to 1000 km at noon and then decreases in a symmetrical fashion during the second half of the day. The h m F2 has highest values at midnight (around 400 km) and lowest values during day (varying between 300 and 330 km). Notice the detailed density dis- tribution above the h m F 2 . It is easy to detect the changes in the calculated topside scale height and resulting density distribution as they develop in time. 5.4. A note on the operational reconstruction Recently, an operational procedure has been developed for real-time reconstruction of the vertical electron profile at a single ionosonde location (Stankov et al., 2002b). In general, it is a stand-by procedure; its exe- cution is triggered by either time control system N h N h F h h F H e m m( ) = ( ) − ++ + O O sech2 2 2 2 ( ) . N h F h h F H m m+ ( ) − H H sech+ + 2 2 2 2 ( ) . TEC O O O Ht m mH N h F H N h F= ( ) + ( )+ + + +2 322 2. . . . . N h F N h F N Fm m mO H+ +( ) + ( ) =2 2 2 N h F h h F H m tr m O O sech+ + ( ) − =2 2 2 2 ( ) . N h F h h F H m tr m H O sech+ + = ( ) − 2 2 2 32 ( ) . 644 Jean-Claude Jodogne and Stanimir M. Stankov Fig. 7. Diurnal behaviour reconstruction: reconstructed vertical electron density (log scale, m– 3) (top), actual GPS TEC and ionosonde measurements (middle), and model-based UTL and H m F 2 (bottom). 645 Ionosphere-plasmasphere response to geomagnetic storms studied with the RMI-Dourbes comprehensive database or new measurements. Thus, it relies heavily on regular discrete influx of (ionosonde and TEC) measurement data. Both types of observations should be synchronized, so representative data are obtained for given location and time. Highest flexibility, in terms of time resolution, is offered by the digital ionosonde – a new block of meas- urement data is available within a delay of about 5 min. A longer delay is expected for receiving the GPS TEC value because the TEC derivation procedure (Warnant, 1997, 1998) requires time and a sufficient number of measurements to evaluate the hardware (receiver plus satellite) biases and detect possible irregularities (TID and scintillation effects). In practice, a TEC value can be obtained every 15 min, which is sufficient for most applications. Several distinct stages are observed in the operational reconstruction procedure: trans- mission of measurement data and retrieval of input parameters, construction of the bottom- and top- side electron profile, backup and display of results. The measurement data are transmitted using the File Transfer Protocol; the UTL values are provided by the empirical formulae (Kutiev et al., 1994) incorporated into the reconstruction software. If the TEC value is not available on time, it is possible to use the ionosonde-based TEC value. The final stage of the operational procedure is to store and display the results. As said before, the next round of calculations can be triggered by time control or the arrival of new measurements. 6. Conclusions Presented was a review of some important ionospheric response studies carried out at the Dourbes Geophysics Centre of the Royal Meteor- ological Institute. The studies were carried out with the help of the local comprehensive data- base. Three major topics of importance to the ion- ospheric response investigations were covered – case studies of storms, development of local geomagnetic activity indices for use in the short- term forecasting, and reconstruction of the vertical electron density distribution. According to each of the above subjects, the following conclusions can be made: Ionospheric storm studies: – The storm studies are much more reliable and efficient if combined ionosonde and GPS TEC measurements are used. – Higher resolution data (e.g., every 15 min) are required for a better understanding of the dynamic ionospheric storm behaviour on a global scale. – The local measurements of the horizontal geomagnetic field component can be used successfully as a reference for storm onset and subsequent main and recovery phases. Short-term forecasting of f 0 F 2 and TEC: – The development of local synthetic geo- magnetic activity indices is of crucial impor- tance to the f0F2 and TEC short-term forecasting. Such indices are local storm-response char- acteristics of high precision. – Differences are observed between the f0F2 and GPS TEC behaviour response to intense magnetic activity (K p > 4), mainly with the detected positive response of TEC during winter in contrast with the all-negative f0F2 relative deviations. Electron profile reconstruction: – A recently-developed electron density recon- struction technique can be successfully utilized in investigating the local storm-time ionosphere development. – A new operational procedure is being devel- oped, which will provide the reconstructed electron density on a near real-time basis. – An advantage of this method is its ap- plicability on a global scale through the ever- growing GPS TEC and ionosonde measurements network. In summary, the overall evaluation of the Dourbes measurement database proved a reliable research reference that will help to accomplish the above-stated goals. Acknowledgements This research was financially supported by the Belgian Federal Office for Scientific, Tech- nical and Cultural Affairs and by the Royal Me- 646 Jean-Claude Jodogne and Stanimir M. Stankov teorological Institute of Belgium. The authors thank Dr. R. Warnant for the useful discussions and GPS TEC data processing. REFERENCES AKASOFU, S.I. and S. CHAPMAN (1972): Solar-Terrestrial Physics (Oxford University Press, Oxford), pp. 901. BARNDEN L.R. (1973): The large-scale magnetic field configuration associated with Forbush decrease, in Proceedings XIII International Cosmic Rays Conference, Denver, 2, 1277-1282. BAROUCH, E. and L.F. BURLAGA (1975): Causes of Forbush decreases and other cosmic ray variations, J. Geophys. Res., 80, 449. CANE, H.V., I.G. RICHARDSON and T.T. VON ROSENVINGEN (1993): Cosmic ray decreases and particle acceleration in 1978-1982 and the associated solar wind structures, J. Geophys. Res., 98 (A8), 13,295-13,302. CANE, H.V., I.G. RICHARDSON, T.T. VON ROSENVINGEN and G. WIBBERENZ (1994): Cosmic ray decreases and shock structure: a multispacecraft study, J. Geophys. Res., 99 (A11), 21,429-21,441. DAVIES, K. (1990): Ionospheric Radio (Peter Peregrinus Ltd., London), pp. 580. FORBUSH, S.E. (1938): On the world-wide changes in cosmic- ray intensity, Phys. Rev., 54, 975. FULLER-ROWELL, T.J., M.V. CODRESCU, H. RISHBETH, R.J. MOFFETT and S. QUEGAN (2000): On the seasonal response of the thermosphere and ionosphere to geomagnetic storms, J. Geophys. Res., 101, 2343-2353. FOERSTER, M. and N. JAKOWSKI (2000): Geomagnetic storm effects on the topside ionosphere and plasmasphere, Surv. Geophys., 21, 47-87. GONZALEZ, W.D., B.T. TSURUTANI, R.P. LEPPING and R. SCHWEEN (2002): Interplanetary phenomena associated with very intense geomagnetic storms, J. Atmos. Sol. Terr. Phys., 64 (2), 173-181. HA R G R E AV E S, J.K. (1992): The Solar-Terrestrial Environment (Cambridge University Press, Cambridge), pp. 420. HERNANDEZ, G. and R.G. ROBLE (1978): Observations of large-scale thermospheric waves during geomagnetic storms, J. Geophys. Res., 83, 5531-5538. HUANG, X. and B. REINISCH (1996): Vertical electron density profiles from the Digisonde network, Adv. Space Res., 18 (6), 121-129. HUANG, X. and B. REINISCH (2001): Vertical total electron content from ionograms in real time, Radio Sci., 36, 335-342. JAKOWSKI, N., S. SCHLUTER and E. SARDON (1999): Total electron content of the ionosphere during the geomagnetic storm on 10 January 1997, J. Atmos. Sol. Terr. Phys., 61, 299-307. JANKOWSKI, J. and C. SUCKSDORFF (1996): IAGA Guide for Magnetic Measurements and Observatory Practice, IAGA, Warsaw, Poland. KUTIEV, I., S.M. STANKOV and P. MARINOV (1994): Analytical expression of O+-H+ transition surface for use in IRI, Adv. Space Res., 17 (12), 135-138. MATSUSHITA, S. (1967): Geomagnetic disturbances and storms, in Physics of the Geomagnetic Phenomena, edited by S. MATSUSHITA and W.H. CAMPBELL (Academic Press, NewYork), 1153-1202. MENIETTI, J.D., J.L. BURCH and D.L. GALLAGHER (1988): Statistical study of ion flows in the dayside and nightside plasmasphere, Planet. Space Sci., 36 (7), 693-702. MENVIELLE, M. and A. BERTHELIER (1991): The K-derived planetary indices: description and availability, Rev. Geophys., 29 (3), 415-432. MIKHAILOV, A.V. and M. FOERSTER (1999): Some F2-layer effects during the January 06-11, 1997 CEDAR storm period as observed with the Millstone Hill Incoherent Scatter Facility, J. Atmos. Sol. Terr. Phys., 61, 249-261. MU T H A R O V, P., I. KU T I E V and L. CA N D E R (2002): Geomagnetically correlated autoregression model for short-term prediction of ionospheric parameters, Inverse Problems, 18, 49-65. NAMGALADZE, A.A., M. FOERSTER and R.Y. YURIK (2000): Analysis of the positive ionospheric response to a moderate geomagnetic storm using a global numerical model, Ann. Geophysicae, 18, 461-477. PRIMDHAL, F. (1979): The Fluxgate Magnetometer, Geophysical Papers, R-55, Danish Meteorological Institute. PROELSS, G.W. (1993): On explaining the local time variation of ionospheric storm effect, Ann. Geophysicae, 11 (1), 1-9. PR O E L S S, G.W. (1997): Magnetic storm associated perturbations of the upper atmosphere, in Magnetic Storms, edited by B.T. TSURUTANI, W.D. GONZALEZ, Y. KAMIDE and J.K. ARBALLO, Geophys. Monogr. Ser., Am. Geophys. Un., Washington D.C., vol. 98, 227-241. PROELSS, G.W., L.H. BRACE, H.G. MAYR, G.R. CARIGNAN, T.L.KILLEEN and J.A. KLOBUCHAR (1991): Ionospheric storm effects at subauroral latitude: a case study, J. Geophys. Res., 96, 1275-1288. REINISCH, B.W. (1996): Ionosonde, in The Upper Atmos− phere - Data Analysis and Interpretation, edited by W. DIEMINGER, G.K. HARTMANN and R. LEITINGER, Springer, Berlin, 370-381. SHADRINA, L.P., S.A. STARODUBTSEV and I.Y. PLOTINIKOV (2002): Geomagnetic storms and Forbush decreases during the passage of the Earth through the flanks of large-scale solar wind disturbances, Int. J. Geomagn. Aeron., 3 (1), 45-50. STANKOV, S.M. (2002): Empirical modeling of ion transition levels based on satellite in situ measurements, C. R. Acad. Bulg. Sci., 55 (1), 35-40. STANKOV, S.M. and P.M. MUHTAROV (2001): Reconstruction of the electron density profile from the total electron content using upper transition level and vertical incidence sounding measurements, C. R. Acad. Bulg. Sci., 54 (9), 45-48. STANKOV, S.M., I.S. KUTIEV, N. JAKOWSKI and S. HEISE (2002a): Electron density profiles deduced from GPS TEC, O+-H+ transition height and ionosonde data, Acta Geod. Geophys. Hung., 37 (2-3), 171-181. STANKOV, S.M., R. WARNANT and J.C. JODOGNE (2002b): A real-time operational procedure for GPS TEC – based 647 Ionosphere-plasmasphere response to geomagnetic storms studied with the RMI-Dourbes comprehensive database reconstruction of the electron profile at a single ionosonde location, in Proceedings EGS XXVII General Assembly, 21-26 April 2002, Nice, France, Paper n. EGS02–A-05554. STANKOV, S.M., I.S. KUTIEV, N. JAKOWSKI and A. WEHRENPFENNIG (2001): A new method for total electron content forecasting using global positioning system measurements, in Proceedings ESTEC Workshop on Space Weather, 17-19 December 2001, Noordwijk, The Netherlands. SZUSZCZEWICZ, E.P., M. LESTER, P. WILKINSON, P. BLANCHARD, M. ABDU, R. HANBABA, K. IGARASHI, S. PULINETS and B.M. REDDY (1998): A comparative study of global ionospheric responses to intense magnetic storm conditions, J. Geophys. Res., 103 (A6), 11,665- 11,684. WARNANT, R. (1997): Reliability of the TEC computed using GPS measurements – the problem of hardware biases, Acta Geod. Geophys. Hung., 32 (3-4), 451-459. WARNANT, R. (1998): Detection of irregularities in the total electron content using GPS measurements – application to a mid-latitude station, Acta Geod. Geophys. Hung., 33 (1), 121-128. WARNANT, R. and J.C. JODOGNE (1998): A comparison between the TEC computed using GPS and ionosonde measurements, Acta Geod. Geophys. Hung., 33 (1), 147-153. ZHANG, G. and L.F. BURLAGA (1988): Magnetic clouds, geomagnetic disturbances, and cosmic ray decreases, J. Geophys. Res., 93, 2511. (received April 10, 2002; accepted October 15, 2002)