Annales Geophysicae

Abstract. The Cluster mission, ESA’s first cornerstone project, together with the SOHO mission, dating back to the first proposals in 1982, was finally launched in the summer of 2000. On 16 July and 9 August, respectively, two Russian Soyuz rockets blasted off from the Russian cosmodrome in Baikonour to deliver two Cluster spacecraft, each into their proper orbit. By the end of August 2000, the four Cluster satellites had reached their final tetrahedral constellation. The commissioning of 44 instruments, both individually and as an ensemble of complementary tools, was completed five months later to ensure the optimal use of their combined observational potential. On 1 February 2001, the mission was declared operational. The main goal of the Cluster mission is to study the small-scale plasma structures in three dimensions in key plasma regions, such as the solar wind, bow shock, magnetopause, polar cusps, magnetotail and the auroral zones. With its unique capabilities of three-dimensional spatial resolution, Cluster plays a major role in the International Solar Terrestrial Program (ISTP), where Cluster and the Solar and Heliospheric Observatory (SOHO) are the European contributions. Cluster’s payload consists of state-of-the-art plasma instrumentation to measure electric and magnetic fields from the quasi-static up to high frequencies, and electron and ion distribution functions from energies of nearly 0 eV to a few MeV. The science operations are coordinated by the Joint Science Operations Centre (JSOC), at the Rutherford Appleton Laboratory (UK), and implemented by the European Space Operations Centre (ESOC), in Darmstadt, Germany. A network of eight national data centres has been set up for raw data processing, for the production of physical parameters, and their distribution to end users all over the world. The latest information on the Cluster mission can be found at http://sci.esa.int/cluster/.

Abstract. With the recent increase in the satellite-based navigation applications, the ionospheric total electron content (TEC) and the L-band scintillation measurements have gained significant importance. In this paper we present the temporal and spatial variations in TEC derived from the simultaneous and continuous measurements made, for the first time, using the Indian GPS network of 18 receivers located from the equator to the northern crest of the equatorial ionization anomaly (EIA) region and beyond, covering a geomagnetic latitude range of 1° S to 24° N, using a 16-month period of data for the low sunspot activity (LSSA) years of March 2004 to June 2005. The diurnal variation in TEC at the EIA region shows its steep increase and reaches its maximum value between 13:00 and 16:00 LT, while at the equator the peak is broad and occurs around 16:00 LT. A short-lived day minimum occurs between 05:00 to 06:00 LT at all the stations from the equator to the EIA crest region. Beyond the crest region the day maximum values decrease with the increase in latitude, while the day minimum in TEC is flat during most of the nighttime hours, i.e. from 22:00 to 06:00 LT, a feature similar to that observed in the mid-latitudes. Further, the diurnal variation in TEC show a minimum to maximum variation of about 5 to 50 TEC units, respectively, at the equator and about 5 to 90 TEC units at the EIA crest region, which correspond to range delay variations of about 1 to 8 m at the equator to about 1 to 15 m at the crest region, at the GPS L1 frequency of 1.575 GHz. The day-to-day variability is also significant at all the stations, particularly during the daytime hours, with maximum variations at the EIA crest regions. Further, similar variations are also noticed in the corresponding equatorial electrojet (EEJ) strength, which is known to be one of the major contributors for the observed day-to-day variability in TEC. The seasonal variation in TEC maximizes during the equinox months followed by winter and is minimum during the summer months, a feature similar to that observed in the integrated equatorial electrojet (IEEJ) strength for the corresponding seasons. In the Indian sector, the EIA crest is found to occur in the latitude zone of 15° to 25° N geographic latitudes (5° to 15° N geomagnetic latitudes). The EIA also maximizes during equinoxes followed by winter and is not significant in the summer months in the LSSA period, 2004–2005. These studies also reveal that both the location of the EIA crest and its peak value in TEC are linearly related to the IEEJ strength and increase with the increase in IEEJ.

Abstract. The magnetic structure and geomagnetic response of 73 magnetic clouds (MC) observed by the WIND and ACE satellites in solar cycle 23 are examined. The results have been compared with the surveys from the previous solar cycles. The preselected candidate MC events were investigated using the minimum variance analysis to determine if they have a flux-rope structure and to obtain the estimation for the axial orientation (θC, φC). Depending on the calculated inclination relative to the ecliptic we divided MCs into "bipolar" (θC<45°) and "unipolar" (θC>45°). The number of observed MCs was largest in the early rising phase, although the halo CME rate was still low. It is likely that near solar maximum we did not identify all MCs at 1AU, as they were crossed far from the axis or they had interacted strongly with the ambient solar wind or with other CMEs. The occurrence rate of MCs at 1AU is also modified by the migration of the filament sites on the Sun towards the poles near solar maximum and by the deflection of CMEs towards the equator due to the fast solar wind flow from large polar coronal holes near solar minimum. In the rising phase nearly all bipolar MCs were associated with the rotation of the magnetic field from the south at the leading edge to the north at the trailing edge. The results for solar cycles 21-22 showed that the direction of the magnetic field in the leading portion of the MC starts to reverse at solar maximum. At solar maximum and in the declining phase (2000-2003) we observed several MCs with the rotation from the north to the south. We observed unipolar (i.e. highly inclined) MCs frequently during the whole investigated period. For solar cycles 21-22 the majority of MCs identified in the rising phase were bipolar while in the declining phase most MCs were unipolar. The geomagnetic response of a given MC depends greatly on its magnetic structure and the orientation of the sheath fields. For each event we distinguished the effect of the sheath fields and the MC fields. All unipolar MCs with magnetic field southward at the axis were geoeffective (Dst<-50nT) while those with the field pointing northward did not cause magnetic storms at all. About half of the all identified MCs were not geoffective or the sheath fields preceding the MC caused the storm. MCs caused more intense magnetic storms (Dst<-100nT) than moderate magnetic storms (-50nT ≥Dst≥-100nT).

Abstract. On 3 July 2001, the four Cluster satellites traversed along the dawnside magnetospheric flank and observed large variations in all plasma parameters. The estimated magnetopause boundary normals were oscillating in the z-direction and the normal component of the magnetic field showed systematic  2–3 min bipolar variations for 1 h when the IMF had a small positive bz-component and a Parker-spiral orientation in the x,y-plane. Brief  33 s intervals with excellent deHoffman Teller frames were observed satisfying the Walén relation. Detailed comparisons with 2-D MHD simulations indicate that Cluster encountered rotational discontinuities generated by Kelvin-Helmholtz instability. We estimate a wave length of  6 RE and a wave vector with a significant z-component.

Abstract. In this work, cluster and principal component analysis are used to divide Spain in a limited number of climatically homogeneous zones, based on seasonal rainfall for 32 Spanish localities for the period 1912–2000. Using the hierarchical technique of clustering Ward's method, three clusters have been obtained in winter and spring, and four clusters have been obtained in summer and autumn. Results are similar to those obtained by applying principal component analysis. Centroid series of each cluster and principal component series of each EOF have been compared to analyze the temporal patterns. The comparison of both methods indicates that cluster analysis is suitable to establish spatio-temporal patterns of seasonal rainfall distribution in Spain. Key words. Meteorology and atmospheric dynamics (climatology; precipitation; general or miscellaneous)

Abstract. CHAMP (CHAllenging Minisatellite Payload) and GRACE (Gravity Recovery And Climate Experiment) formed a satellite configuration for precise atmospheric sounding during the first activation of the GPS (Global Positioning System) radio occultation experiment aboard GRACE on 28 and 29 July 2004. 338 occultations were recorded aboard both satellites, providing globally distributed vertical profiles of refractivity, temperature and specific humidity. The combined set of CHAMP and GRACE profiles shows excellent agreement with meteorological analysis. Almost no refractivity bias is observed between 5 and 30km, the standard deviation is between 1 and 2% within this altitude interval. The GRACE satellite clock stability is significantly improved in comparison with CHAMP. This allows for the application of a zero difference technique for precise analysis of the GRACE occultation data.

Abstract. The strong increase in the flux of relativistic electrons during the recovery phase of magnetic storms and during other active periods is investigated with the help of Hamiltonian formalism and simulations of test electrons which interact with whistler waves. The intensity of the whistler waves is enhanced significantly due to injection of 10-100 keV electrons during the substorm. Electrons which drift in the gradient and curvature of the magnetic field generate the rising tones of VLF whistler chorus. The seed population of relativistic electrons which bounce along the inhomogeneous magnetic field, interacts resonantly with the whistler waves. Whistler wave propagating obliquely to the magnetic field can interact with energetic electrons through Landau, cyclotron, and higher harmonic resonant interactions when the Doppler-shifted wave frequency equals any (positive or negative) integer multiple of the local relativistic gyrofrequency. Because the gyroradius of a relativistic electron may be the order of or greater than the perpendicular wavelength, numerous cyclotron, harmonics can contribute to the resonant interaction which breaks down the adiabatic invariant. A similar process diffuses the pitch angle leading to electron precipitation. The irreversible changes in the adiabatic invariant depend on the relative phase between the wave and the electron, and successive resonant interactions result in electrons undergoing a random walk in energy and pitch angle. This resonant process may contribute to the 10-100 fold increase of the relativistic electron flux in the outer radiation belt, and constitute an interesting relation between substorm-generated waves and enhancements in fluxes of relativistic electrons during geomagnetic storms and other active periods.Key words. Magnetospheric physics (energetic particles · trapped; plasma waves and instabilities; storms and substorms)

Abstract. A physical model of the coupled thermosphere and ionosphere has been used to determine the accuracy of model predictions of the ionospheric response to geomagnetic activity, and assess our understanding of the physical processes. The physical model is driven by empirical descriptions of the high-latitude electric field and auroral precipitation, as measures of the strength of the magnetospheric sources of energy and momentum to the upper atmosphere. Both sources are keyed to the time-dependent TIROS/NOAA auroral power index. The output of the model is the departure of the ionospheric F region from the normal climatological mean. A 50-day interval towards the end of 1997 has been simulated with the model for two cases. The first simulation uses only the electric fields and auroral forcing from the empirical models, and the second has an additional source of random electric field variability. In both cases, output from the physical model is compared with F-region data from ionosonde stations. Quantitative model/data comparisons have been performed to move beyond the conventional "visual" scientific assessment, in order to determine the value of the predictions for operational use. For this study, the ionosphere at two ionosonde stations has been studied in depth, one each from the northern and southern mid-latitudes. The model clearly captures the seasonal dependence in the ionospheric response to geomagnetic activity at mid-latitude, reproducing the tendency for decreased ion density in the summer hemisphere and increased densities in winter. In contrast to the "visual" success of the model, the detailed quantitative comparisons, which are necessary for space weather applications, are less impressive. The accuracy, or value, of the model has been quantified by evaluating the daily standard deviation, the root-mean-square error, and the correlation coefficient between the data and model predictions. The modeled quiet-time variability, or standard deviation, and the increases during geomagnetic activity, agree well with the data in winter, but is low in summer. The RMS error of the physical model is about the same as the IRI empirical model during quiet times. During the storm events the RMS error of the model improves on IRI, but there are occasionally false-alarms. Using unsmoothed data over the full interval, the correlation coefficients between the model and data are low, between 0.3 and 0.4. Isolating the storm intervals increases the correlation to between 0.43 and 0.56, and by smoothing the data the values increases up to 0.65. The study illustrates the substantial difference between scientific success and a demonstration of value for space weather applications.Key words: Ionosphere (ionospheric disturbances; mid-latitude ionosphere; modeling and forecasting)

Abstract. We have used time-delay feed-forward neural networks to compute the geomagnetic-activity index Dst one hour ahead from a temporal sequence of solar-wind data. The input data include solar-wind density n, velocity V and the southward component Bz of the interplanetary magnetic field. Dst is not included in the input data. The networks implement an explicit functional relationship between the solar wind and the geomagnetic disturbance, including both direct and time-delayed non-linear relations. In this study we especially consider the influence of varying the temporal size of the input-data sequence. The networks are trained on data covering 6600 h, and tested on data covering 2100 h. It is found that the initial and main phases of geomagnetic storms are well predicted, almost independent of the length of the input-data sequence. However, to predict the recovery phase, we have to use up to 20 h of solar-wind input data. The recovery phase is mainly governed by the ring-current loss processes, and is very much dependent on the ring-current history, and thus also the solar-wind history. With due consideration of the time history when optimizing the networks, we can reproduce 84% of the Dst variance.