Study of the magnetic turbulence in a corotating interaction region in the interplanetary medium

We study the geometry of magnetic fluctuations in a CIR observed by Pioneer 10 at 5 AU between days 292 and 295 in 1973. We apply the methodology proposed by Bieber et al. to make a comparison of the relative importance of two geometric arrays of vector propagation of the magnetic field fluctuations: slab and two-dimensional (2D). We found that inside the studied CIR this model is not applicable due to the restrictions imposed on it. Our results are consistent with Alfvenic fluctuations propagating close to the radial direction, confirming Mavromichalaki et al.’s findings. A mixture of isotropic and magnetoacoustic waves in the region before the front shock would be consistent with our results, and a mixture of slab/2D and magnetoacoustic waves in a region after the reverse shock. We base the latter conclusions on the theoretical analysis made by Kunstmann. We discuss the reasons why the composite model can not be applied in the CIR studied although the fluctuations inside it are two dimensional.     

Interannual signals in length of day and atmospheric angular momentum

Atmospheric angular momentum (AAM) and length of day (LOD) series are investigated for their characteristics on interannual time scales during the half-century period 1949 to 1998. During this epoch, the interannual variability in LOD can be separated naturally into three bands: a quasi-biennial, a triennial-quadrennial and one at six-seven years. The atmosphere appears to excite the first two bands, while it does not contribute to the last. Considering the quasi-biennial (QB) band alone, the atmosphere appears to excite most of its signal in LOD, but it arises from separate fluctuations with stratospheric and tropospheric origin. Thus, although close in frequency, stratospheric and tropospheric processes differ in their amplitude and phase variability. The time shift can be noted especially during the strong El Niño events of 1982/83 and 1997/98 when both processes have positive phase and thus combine to help produce particularly strong peak in AAM and LOD. In addition, we have reconfirmed the downward propagation in the stratosphere and upward propagation in the troposphere of AAM observed in earlier studies for other variables. In the triennial-quadrennial (TQ) band, time-variable spectral analyses reveal that LOD and AAM contain strong variability, with periods shorter than four years before 1975 and longer thereafter. This signal originates mainly within the troposphere and propagates upwards from the lower to the higher layers of the troposphere. According to a zonal analysis, an equatorial poleward mode, strongly linked to the SOI, explains more than 60% of the total variability at these ranges. In addition, this study also indicates that an equatorward mode, originating within polar latitudes, explains, on average, more than 15% of the triennial-quadrennial oscillation (TQO) variability in AAM, and up to 30% at certain epochs. Finally, a six year period in LOD noted in earlier studies, as well as in lengthier series covering much of the century, is found to be absent in atmospheric excitations, and it is thus likely to arise from mantle/core interactions.     

Detection of zero anisotropy at 5.2 AU during the November 1998 solar particle event: Ulysses Anisotropy Telescopes observations

For the first time during the mission, the Anisotropy Telescopes instrument on board the Ulysses spacecraft measured constant zero anisotropy of protons in the 1.3/2.2 MeV energy range, for a period lasting more than three days. This measurement was made during the energetic particle event taking place at Ulysses between 25 November and 15 December 1998, an event characterised by constant high proton fluxes within a region delimited by two interplanetary forward shocks, at a distance of 5.2 AU from the Sun and heliographic latitude of 17°S. We present the ATs results for this event and discuss their possible interpretation and their relevance to the issue of intercalibration of the two telescopes.     

Climate hypersensitivity to solar forcing?

We compare the equilibrium climate responses of a quasi-dynamical energy balance model to radiative forcing by equivalent changes in CO₂, solar total irradiance (S_tot) and solar UV (S_UV). The response is largest in the S_UV case, in which the imposed UV radiative forcing is preferentially absorbed in the layer above 250 mb, in contrast to the weak response from global-columnar radiative loading by increases in CO₂ or S_tot. The hypersensitive response of the climate system to solar UV forcing is caused by strongly coupled feedback involving vertical static stability, tropical thick cirrus ice clouds and stratospheric ozone. This mechanism offers a plausible explanation of the apparent hypersensitivity of climate to solar forcing, as suggested by analyses of recent climatic records. The model hypersensitivity strongly depends on climate parameters, especially cloud radiative properties, but is effective for arguably realistic values of these parameters. The proposed solar forcing mechanism should be further confirmed using other models (e.g., general circulation models) that may better capture radiative and dynamical couplings of the troposphere and stratosphere.     

Cross-field gradients: general concept, importance of multi-spacecraft measurements and study at 1 AU of the source intensity gradient for E>30 keV solar event electrons

Three main physical processes (and associated properties) are currently used to describe the flux and anisotropy time profiles of solar energetic particle events, called SEP profiles. They are (1) the particle scattering (due to magnetic waves), (2) the particle focusing (due to the decrease of the amplitude of the interplanetary magnetic field (IMF) with the radial distance to the Sun) and (3) the finite injection profile at the source. If their features change from one field line to another, i.e. if there is a cross IMF gradient (CFG), then the shape of the SEP profiles will depend, at onset time, on the relative position of the spacecraft to the IMF and might vary significantly on small distance scale (e.g. 10⁶ km). One type of CFG is studied here. It is called intensity CFG and considers variations, at the solar surface, only of the intensity of the event. It is shown here that drops of about two orders of magnitude over distances of 10⁴ km at the Sun (1° of angular distance) can influence dramatically the SEP profiles at 1 AU. This CFG can lead to either an under or overestimation of both the parallel mean free path and of the injection parameters by factor up to, at least, 2/3 and 18, respectively. Multi-spacecraft analysis can be used to identify CFG. Three basic requirements are proposed to identify, from the observation, the type of the CFG being measured.     

ERNE observations of energetic particles associated with Earth-directed coronal mass ejections in April and May, 1997

Two Earth-directed coronal mass ejections (CMEs), which were most effective in energetic (1–50 MeV) particle acceleration during the first 18 months since the Solar and Heliospheric Observatory (SOHO) launch, occurred on April 7 and May 12, 1997. In the analysis of these events we have deconvoluted the injection spectrum of energetic protons by using the method described by Anttila et al. In order to apply the method developed earlier for data of a rotating satellite (Geostationary Operational Environmental Satellites, GOES), we first had to develop a method to calculate the omnidirectional energetic particle intensities from the observations of Energetic and Relativistic Nuclei and Electrons (ERNE), which is an energetic particle detector onboard the three-axis stabilized SOHO spacecraft. The omnidirectional intensities are calculated by fitting an exponential pitch angle distribution from directional information of energetic protons observed by ERNE. The results of the analysis show that, compared to a much faster and more intensive CMEs observed during the previous solar maximum, the acceleration efficiency decreases fast when the shock propagates outward from the Sun. The particles injected at distances <0.5 AU from the Sun dominate the particle flux during the whole period, when the shock propagates to the site of the spacecraft. The main portion of particles injected by the shock during its propagation further outward from the Sun are trapped around the shock, and are seen as an intensity increase at the time of the shock passage.     

Observations of interplanetary scintillation during the 1998 Whole Sun Month: a comparison between EISCAT, ORT and Nagoya data

Observations of interplanetary scintillation (IPS) allow accurate solar wind velocity measurements to be made at all heliographic latitudes and at a range of distances from the Sun. The data may be obtained with either single, double or multiple antennas, each requiring a different method of analysis. IPS data taken during the 1998 whole sun month (30th July–31st August 1998) by EISCAT, the ORT (Ooty Radio Telescope), India, and the Nagoya IPS system, Japan, allow the results of individual methods of analysis to be compared. Good agreement is found between the velocity measurements using each method, and when combined an improved understanding of the structure of the solar wind can be obtained.On leave from the Physical Research Laboratory, Ahmedabad 380 009, India     

Ionospheric measurements of relative coronal brightness during the total solar eclipses of 11 August, 1999 and 9 July, 1945

Swept-frequency (1/10 MHz) ionosonde measurements were made at Helston, Cornwall (50°06N, 5°18W) during the total solar eclipse on August 11, 1999. Soundings were made every three minutes. We present a method for estimating the percentage of the ionising solar radiation which remains unobscured at any time during the eclipse by comparing the variation of the ionospheric E-layer with the behaviour of the layer during a control day. Application to the ionosonde date for 11 August, 1999, shows that the flux of solar ionising radiation fell to a minimum of 25±2% of the value before and after the eclipse. For comparison, the same technique was also applied to measurements made during the total solar eclipse of 9 July, 1945, at Sörmjöle (63°68N, 20°20E) and yielded a corresponding minimum of 16 ± 2%. Therefore the method can detect variations in the fraction of solar emissions that originate from the unobscured corona and chromosphere. We discuss the differences between these two eclipses in terms of the nature of the eclipse, short-term fluctuations, the sunspot cycle and the recently-discovered long-term change in the coronal magnetic field.     

Simulating solar MHD

Two aspects of solar MHD are discussed in relation to the work of the MHD simulation group at KIS. Photospheric magneto-convection, the nonlinear interaction of magnetic field and convection in a strongly stratified, radiating fluid, is a key process of general astrophysical relevance. Comprehensive numerical simulations including radiative transfer have significantly improved our understanding of the processes and have become an important tool for the interpretation of observational data. Examples of field intensification in the solar photosphere (‘convective collapse’) are shown. The second line of research is concerned with the dynamics of flux tubes in the convection zone, which has far-reaching implications for our understanding of the solar dynamo. Simulations indicate that the field strength in the region where the flux is stored before erupting to form sunspot groups is of the order of 10⁵ G, an order of magnitude larger than previous estimates based on equipartition with the kinetic energy of convective flows.     

Long-duration high-energy proton events observed by GOES in October 1989

We consider the prolonged injection of the high-energy (> 10 MeV) protons during the three successive events observed by GOES in October 1989. We apply a solar-rotation-stereoscopy approach to study the injection of the accelerated particles from the CME-driven interplanetary shock waves in order to find out how the effectiveness of the particle acceleration and/or escape depends on the angular distance from the shock axis. We use an empirical model for the proton injection at the shock and a standard model of the interplanetary transport. The model can reproduce rather well the observed intensity-time profiles of the October 1989 events. The deduced proton injection rate is highest at the nose of the shock; the injection spectrum is always harder near the Sun. The results seem to be consistent with the scheme that the CME-driven interplanetary shock waves accelerate a seed particle population of coronal origin.     

Probing the magnetic topology of coronal mass ejections by means of Ulysses/HI-SCALE energetic particle observations

In this work, solar flare energetic particle fluxes (E_e 42 keV) observed by the HI-SCALE instrument onboard Ulysses, a spacecraft that is probing the heliosphere in 3-D, are utilized as diagnostics of the large-scale structure and topology of the interplanetary magnetic field (IMF) embedded within two well-identified interplanetary coronal mass ejection (ICME) structures. On the basis of the energetic solar flare particle observations firm conclusions are drawn on whether the detected ICMEs have been detached from the solar corona or are still magnetically anchored to it when they arrive at 2.5 AU. From the development of the angular distributions of the particle intensities, we have inferred that portions of the ICMEs studied consisted of both open and closed magnetic field lines. Both ICMEs present a filamentary structure comprising magnetic filaments with distinct electron anisotropy characteristics. Subsequently, we studied the evolution of the anisotropies of the energetic electrons along the magnetic field loop-like structure of one ICME and computed the characteristic decay time of the anisotropy which is a measure of the amount of scattering that the trapped electron population underwent after injection at the Sun.     

Observations of the response time of high-latitude ionospheric convection to variations in the interplanetary magnetic field using EISCAT and IMP-8 data

We have combined ∼300 h of tristatic measurements of the field-perpendicular F region ionospheric flow measured overhead at Tromsø by the EISCAT UHF radar, with simultaneous IMP-8 measurements of the solar wind and interplanetary magnetic field (IMF) upstream of the Earth’s magnetosphere, in order to examine the response time of the ionospheric flow to changes in the north-south component of the IMF (B_z). In calculating the flow response delay, the time taken by field changes observed by the spacecraft to first effect the ionosphere has been carefully estimated and subtracted from the response time. Two analysis methods have been employed. In the first, the flow data were divided into 2 h-intervals of magnetic local time (MLT) and cross-correlated with the “half-wave rectifier” function V²B_s, where V is the solar wind speed, and B_s is equal to IMF B_z if the latter is negative, and is zero otherwise. Response delays, determined from the time lag of the peak value of the cross-correlation coefficient, were computed versus MLT for both the east-west and north-south components of flow. The combined data set suggests minimum delays at ∼1400 MLT, with increased response times on the nightside. For the 12-h sector centred on 1400 MLT, the weighted average response delay was found to be 1.3 ± 0.8 min, while for the 12-h sector centred on 0200 MLT the weighted average delay was found to increase to 8.8 ± 1.7 min. In the second method we first inspected the IMF data for sharp and enduring (at least ∼5 min) changes in polarity of the north-south component, and then examined concurrent EISCAT flow data to determine the onset time of the corresponding enhancement or decay of the flow. For the case in which the flow response was timed from whichever of the flow components responded first, minimum response delays were again found at ∼1400 MLT, with average delays of 4.8 ± 0.5 min for the 12-h sector centred on 1400 MLT, increasing to 9.2 ± 0.8 min on the nightside. The response delay is thus found to be reasonably small at all local times, but typically ∼6 min longer on the nightside compared with the dayside. In order to make an estimate of the ionospheric information propagation speed implied by these results, we have fitted a simple theoretical curve to the delay data which assumes that information concerning the excitation and decay of flow propagates with constant speed away from some point on the equatorward edge of the dayside open-closed field line boundary, taken to lie at 77° magnetic latitude. For the combined cross-correlation results the best-fit epicentre of information propagation was found to be at 1400 MLT, with an information propagation phase speed of 9.0 km s⁻¹. For the combined event analysis, the best-fit epicentre was also found to be located at 1400 MLT, with a phase speed of 6.8 km s⁻¹.     

Transport of thermal plasma above the auroral ionosphere in the presence of electrostatic ion-cyclotron turbulence

The electron component of intensive electric currents flowing along the geomagnetic field lines excites turbulence in the thermal magnetospheric plasma. The protons are then scattered by the excited electromagnetic waves, and as a result the plasma is stable. As the electron and ion temperatures of the background plasma are approximately equal each other, here electrostatic ion-cyclotron (EIC) turbulence is considered. In the nonisothermal plasma the ion-acoustic turbulence may occur additionally. The anomalous resistivity of the plasma causes large-scale differences of the electrostatic potential along the magnetic field lines. The presence of these differences provides heating and acceleration of the thermal and energetic auroral plasma. The investigation of the energy and momentum balance of the plasma and waves in the turbulent region is performed numerically, taking the magnetospheric convection and thermal conductivity of the plasma into account. As shown for the quasi-steady state, EIC turbulence may provide differences of the electric potential of δ V ≈ 1–10 kV at altitudes of 500 < h < 10 000 km above the Earth’s surface. In the turbulent region, the temperatures of the electrons and protons increase only a few times in comparison with the background values.     

Heliospheric current sheet inclinations at Venus and Earth

We investigate the inclinations of heliospheric current sheet at two sites in interplanetary space, which are generated from the same solar source. From the data of solar wind magnetic fields observed at Venus (0.72 AU) and Earth (1 AU) during December 1978–May 1982 including the solar maximum of 1981, 54 pairs of candidate sector boundary crossings are picked out, of which 16 pairs are identified as sector boundaries. Of the remainder, 12 pairs are transient structures both at Venus and Earth, and 14 pairs are sector boundaries at one site and have transient structures at the other site. It implies that transient structures were often ejected from the coronal streamer belt around the solar maximum. For the 16 pairs of selected sector boundaries, we determine their normals by using minimum variance analysis. It is found that most of the normal azimuthal angles are distributed between the radial direction and the direction perpendicular to the spiral direction both at Venus and Earth. The normal elevations tend to be smaller than ≈45° with respect to the solar equatorial plane, indicating high inclinations of the heliospheric current sheet, in particular at Earth. The larger scatter in the azimuth and elevation of normals at Venus than at Earth suggests stronger effects of the small-scale structures on the current sheet at 0.72 AU than at 1 AU. When the longitude difference between Venus and Earth is small (<40° longitudinally), similar or the same inclinations are generally observed, especially for the sector boundaries without small-scale structures. This implies that the heliospheric current sheet inclination tends to be maintained during propagation of the solar wind from 0.72 AU to 1 AU. Detailed case studies reveal that the dynamic nature of helmet streamers causes variations of the sector boundary structure.     

Geomagnetic storm effects at low latitudes

The geomagnetic horizontal (H) field from the chain of nine observatories in India are used to study the storm-time and disturbance daily variations. The peak decrease in storm-time variation in H showed significant enhancements at the equatorial electrojet stations over and above the normally expected decrease due to the ring current effects corrected for geomagnetic latitudes. The disturbance daily variation of H at equatorial stations showed a large decrease around midday hours over and above the usual dawn-maximum and dusk-minimum seen at any mid-latitude stations around the world. These slow and persistent additional decreases of H of disturbance daily variation at equatorial latitudes could be the effect of a westward electric field due to the Disturbance Ionospheric dynamo coupled with abnormally large electrical conductivities in the E region over the equator.     

Application of nonlinear forecasting techniques for meteorological modeling

A nonlinear forecasting method was used to predict the behavior of a cloud coverage time series several hours in advance. The method is based on the reconstruction of a chaotic strange attractor using four years of cloud absorption data obtained from half-hourly Meteosat infrared images from Northwestern Spain. An exhaustive nonlinear analysis of the time series was carried out to reconstruct the phase space of the underlying chaotic attractor. The forecast values are used by a non-hydrostatic meteorological model ARPS for daily weather prediction and their results compared with surface temperature measurements from a meteorological station and a vertical sounding. The effect of noise in the time series is analyzed in terms of the prediction results.     

Measurements of the direction of the solar wind using interplanetary scintillation

EISCAT observations of the interplanetary scintillation of a single source were made over an extended period of time, during which the orientation of the baselines between the two observing sites changed significantly. Assuming that maximum correlation between the scintillations observed at the two sites occurs when the projected baseline is parallel to the direction of plasma flow, this technique can be used to make a unique determination of the direction of the solar wind. In the past it has usually been assumed that the plasma flow is radial, but measurements of eleven sources using this technique have indicated conclusively that in at least six cases observed at mid or high heliocentric latitude there is a significant non-radial component directed in four cases towards the heliocentric equator and in two cases towards the pole.     

Effects of substorms on the stormtime ring current index Dst

There has been some discussion in recent times regarding whether or not substorm expansive phase activity plays any role of importance in the formation of the stormtime ring current. I explore this question using the K_p index as a proxy for substorm expansive phase activity and the D_st index as a proxy for symmetric ring current strength. I find that increases in D_st are mildly related to the strength of substorm expansive phase activity during the development of the storm main phase. More surprisingly, I find that the strength of D_st during the storm recovery phase is positively correlated with the strength of substorm expansive phase activity. This result has an important bearing on the question of how much the Dst index reflects activity other than that of the stormtime symmetric ring current strength for which it is supposed to be a proxy.     

A simple oblique dip model for geomagnetic micropulsations

It is pointed out that simple models adopted so far have tended to neglect the obliquity of the magnetic field lines entering the Earths surface. A simple alternative model is presented, in which the ambient field lines are straight, but enter wedge shaped boundaries at half a right-angle. The model is illustrated by assuming an axially symmetric, compressional, impulse type disturbance at the outer boundary, all other boundaries being assumed to be perfectly conducting. The numerical method used is checked from the instant the excitation ceases, by an analytical method. The first harmonic along field lines is found to be of noticeable size, but appears to be mainly due to coupling with the fundamental, and with the first harmonic across field lines.     

The response of the azimuthal component of the ionospheric electric field to auroral arc brightening

We have analyzed the response of azimuthal component of the ionospheric electric field to auroral arc activity. We have chosen for analysis three intervals of coordinated EISCAT and TV observations on 18 February, 1993. These intervals include three kinds of arc activity: the appearance of a new auroral arc, the gradual brightening of the existing arc and variations of the arc luminosity. The arcs were mostly east-west aligned. In all cases, the enhancement of arc luminosity is accompanied by a decrease in the westward component of the ionospheric electric field. In contrast, an increase of that component seems to be connected with arc fading. The observed response is assumed to have the same nature as the short circuit of an external electric field by the conductor. The possible consequence of this phenomenon is discussed.