Relationship of the diffuse aurora and SAR arc dynamics to substorms and storms

Investigation results of a diffuse aurora (DA) and stable auroral red (SAR) arc dynamics based on spectrophotometric observations at the Yakutsk meridian (199°E geomagnetic longitude) are presented. The relationship of an equatorward extension of DA in the 557.7 nm emission to a substorm growth phase during the magnetospheric convection intensification after the turn of IMF B_Z to the south is shown. The formation of SAR arc during the substorm expansion phase is investigated. The association of SAR arc dynamics with the development of asymmetric ring current (substorm injection) during the main phase of a storm is analyzed. It is shown how the pulsating precipitations of energetic ring current particles develop in the outer plasmasphere based on photometric observations.     

Magnetospheric Boundary Layer Structure and Dynamics as Seen From Cluster and Double Star Measurements

In this review, we discuss the structure and dynamics of the magnetospheric Low-Latitude Boundary Layer (LLBL) based on recent results from multi-satellite missions Cluster and Double Star. This boundary layer, adjacent to the magnetopause on the magnetospheric side, usually consists of a mixture of plasma of magnetospheric and magnetosheath origins, and plays an important role in the transfer of mass and energy from the solar wind into the magnetosphere and subsequent magnetospheric dynamics. During southward Interplanetary Magnetic Field (IMF) conditions, this boundary layer is generally considered to be formed as a result of the reconnection process between the IMF and magnetospheric magnetic field lines at the dayside magnetopause, and the structure and plasma properties inside the LLBL can be understood in terms of the time history since the reconnection process. During northward IMF conditions, the LLBL is usually thicker, and has more complex structure and topology. Recent observations confirm that the LLBL observed at the dayside can be formed by single lobe reconnection, dual lobe reconnection, or by sequential dual lobe reconnection, as well as partially by localized cross-field diffusion. The LLBL magnetic topology and plasma signatures inside the different sub-layers formed by these processes are discussed in this review. The role of the Kelvin-Helmholtz instability in the formation of the LLBL at the flank magnetopause is also discussed. Overall, we conclude that the LLBL observed at the flanks can be formed by the combination of processes, (dual) lobe reconnection and plasma mixing due to non-linear Kelvin-Helmholtz waves.      

Impulsive plasma waves observed by DE-1 in the magnetosphere

Impulsive plasma waves (1–9 kHz) with durations less than 100 msec have been found in DE-1 wide-band electric field data (650 Hz – 40 kHz) received at Kashima, Japan. The waves are associated with a strong narrow-band ELF hiss, and were observed at geocentric distances from 3.1 to 4.9 Re (earth's radius) in the low-latitude nightside magnetosphere. Local electron densities and plasmapause locations estimated suggest that the waves were observed outside the nightside plasmapause. The waves are discussed in terms of Landau resonant trapping of magnetospheric electrons by the associated whistler-mode ELF hiss.     

Case study of September 24–26, 1998 magnetic storm

We compute global magnetospheric parameters based upon solar wind data obtained from the WIND spacecraft upstream. Using the paraboloid magnetospheric model, calculations of the dynamic global magnetospheric current systems have been made. The solar wind dynamic pressure, the interplanetary magnetic field, the strength of the tail current, and the ring current control the polar cap and auroral oval size and location during the magnetic storm. The model calculations demonstrate that the polar cap and the auroral oval areas are mainly controlled by the tail current. The substorm onset at 0630 UT on September 25, 1998 happened near the minimum in the main phase field depression. The substorm expansion onset time is also marked by a sudden enhancement in the solar wind dynamic pressure and an enhancement in the tail current. The magnetic signatures of these two effects cancel each other, which explains why the D_st profile shows no strong time variation during the substorm. Evidence for the substorm expansion includes not only the signature in the AL index but also the strong asymmetry of the low latitude magnetic disturbances (substorm positive bay signature). Model calculations were checked by comparison with the GOES 8 and 10 magnetic field measurements.     

Interconnection of high-latitude and low-latitude boundary layers when IMF BY is dominant

The antiparallel merging model places the location of the reconnection region for a dominant interplanetary magnetic field (IMF) B_Y at high latitudes at the dayside magnetopause and predicts that the low-latitude boundary layer (LLBL) is located on open field lines of the magnetospheric flanks. Interball-1 data obtained in the wide local time range near the low-latitude magnetopause makes it possible to analyze the LLBL plasma population and to find a link between possible reconnection at high latitudes and LLBL occurrence. We found that no boundary layer was observed in the regions which have no topological connection with the merging site. All cases of LLBL observations are located downstream from a specific boundary. This boundary coincides with the first magnetospheric field line touching the reconnection region and can be located in a wide local time region depending on the instant IMF direction. Even the LLBL on closed field lines shows the tendency to be concentrated in the vicinity of this boundary. Thus we show that all types of observed LLBLs are linked to reconnection sites predicted by the antiparallel merging model.     

MHD surface waves – a source of global magnetospheric modes (resonances)?

Various mechanisms have been proposed for explanation of the global magnetospheric modes whose frequency peaks are in the frequency range 1–4 mHz. Recent papers claim: basic characteristics of the 1–4 mHz activity events observed on ground give evidences for an existence of MHD surface mode excited on the Earth magnetopause. The discrete frequencies of such MHD surface wave modes suggest an emergence of standing wave structures along the magnetic field lines lying on the magnetopause. Such discrete frequencies of MHD surface waves on magnetopause however, are not stable, at all. Contrariwise, MHD surface wave modes supported by the two plasma boundaries – the magnetopause and the plasmapause, are in accordance with existing experimental facts: discrete set of almost stable frequencies, field amplitude peaks and positions, energy dissipation, and field distribution from high to low latitudes. Mechanisms of the global magnetospheric mode resonance are pointed out as well as tools for their identification and discrimination.     

Magnetopause poleward of the cusp: Comparison of plasma and magnetic signature of the boundary for southward and northward directed interplanetary magnetic field

A coherent data set of high-latitude dayside magnetopause encounters by old (Heos 2, Hawkeye, Prognoz 7, 8) and new (Polar, Interball Tail, Cluster) spacecraft is needed to build a realistic model of the magnetopause (MP) including an indentation in the cusp. In building such a coherent data set a caution is necessary as the dayside magnetopause at high-latitudes may be less clearly defined than in the case of observations at low latitudes. It is due to expected presence of bundles of newly-reconnected magnetic field lines forming an extended boundary layer on the magnetosheath (MS) side of the magnetopause in the cusp region. Moreover, numerical magnetohydrodynamic (MHD) models of the solar wind-magnetosphere interaction predict that under northward interplanetary magnetic field (IMF) an additional thin current sheet should form inside the magnetopause at high latitudes on the dayside (e.g., Wu, 1983; Palmroth et al., 2001). Such a thin currect sheet is absent in empirical magnetosphere models. This internal current sheet, if a real one, may be mistaken for the magnetopause if magnetic field data are only taken into account and/or plasma data are unavailable. The Interball-Tail orbit allows for a full transition of magnetopause boundary layers at high-latitudes. We compare plasma and magnetic field signatures of the magnetopause poleward of the cusp for southward and northward IMF. The distance between the magnetic signature of the magnetopause (the current layer) and a cold and laminarly antisunward flowing MS plasma (so called free-flow MS) was found to be 0.5 to 1 R_E, at least. These observations were made under nominal solar wind of v350 km/s and p_dyn=1 to 4 nPa. We also observed several transient magnetic field reversals in the cusp related to pulses of solar wind dynamic pressure and/or the IMF discontinuity arrival. These transient reversals occurred at the same distance to the model MP as well defined full MP crossing, so most probably they represent just short encounters with the magnetopause current layer. Our analysis suggests that an indentation of the magnetopause with a subtle structure dependent on the local magnetic shear would explain and allow to predict the magnetic configuration in the high-altitude cusp.     

Occurrence rate of SAR arcs during the 23rd solar activity cycle

The occurrence rate of SAR arcs during 1997–2007 has been analyzed based on the photometric observations at the Yakutsk meridian (Maimaga station, corrected geomagnetic coordinates: 57°N, 200°E). SAR arcs appeared in 114 cases (∼500 h) during ∼370 nights of observations (∼3170 h). The occurrence frequency of SAR arcs increases to 27% during the growth phase of solar activity and has a clearly defined maximum at a decline of cycle 23. The SAR arc registration probability corresponds to the variations in geomagnetic activity in this solar cycle. The dates, intervals of UT, and geomagnetic latitudes of SAR arc observations at the Yakutsk meridian are presented.     

The automatic detection of geomagnetic-storm sudden commencements

Storm Sudden Commencements (SSCs) are sharp offsets in the global geomagnetic field H component which accompany some geomagnetic storms. SSC signatures, which are a strong function of geographic location and local time, are usually explained as the result of magnetospheric and ionospheric responses to a compression caused by propagating discontinuities in the solar wind. The Space Environment Services Center (SESC) in Boulder, Colorado, has recently instituted an SSC ‘alarm’ to alert both forecasters and customers of potential geomagnetic storm conditions. The present software design looks for rapid change, coincidental within one minute, in the H-components among any two of the three magnetometers received in real-time at the SESC. These are the Boulder magnetometer, and the magnetometers on board the GOES geosynchronous satellites located at 75° and 135° West Longitude. The results are not yet satisfactory and offer some clues toward a better understanding of the physics of SSCs.     

The thermosphere as a sink of magnetospheric energy: A review of recent observations of dynamics

The study of the dynamics and thermodynamics of the earth's upper atmosphere has made significant progress over the past few years owing to the availability of new global-scale data sets from the Dynamics Explorer satellites. The thermospheric wind and temperature fields at high latitude have been observed to depend strongly on forcing processes of magnetospheric origin. A key momentum source is due to the drag effect of ions convecting in response to electric fields mapped down on the ionosphere from magnetospheric boundary regions. Likewise, an important heat source derives from Joule or frictional dissipation due to ion/neutral difference velocities governed, in turn, by magnetospheric forcing. In this paper we discuss the progress made over the last 2–3 years initiated by the new satellite measurements and we review published data on ion and neutral motions in the context of the energy and momentum coupling between the magnetosphere and the ionosphere/neutral upper atmosphere. The observations indicate the existence of a “flywheel effect” which implies direct feedback from the neutral thermosphere to the magnetosphere via the release of energy and momentum previously “stored” in the neutral thermosphere.     

Bow shock: Power aspects

It is clear that the primary energy source for magnetospheric processes is the solar wind, but the process of energy transfer from the solar wind into the magnetosphere, or rather, to convecting magnetospheric plasma, appears to be rather complicated. Bow shock is a powerful transformer of the solar wind kinetic energy into the gas dynamic and electromagnetic energy. A jump of the magnetic field tangential component at front crossing means that the front carries an electric current. The solar wind kinetic energy partly transforms to gas kinetic and electromagnetic energy during its passage through the bow shock front. The transition layer (magnetosheath) can use part of this energy for accelerating of plasma, but can conversely spend part its kinetic energy on the electric power generation, which afterwards may be used by the magnetosphere. Thereby, transition layer can be both consumer (sink) and generator (source) of electric power depending upon special conditions. The direction of the current behind the bow shock front depends on the sign of the IMF Bz-component. It is this electric current which sets convection of plasma in motion.     

Ionospheric response to solar and magnetospheric protons during January 15–22, 2005: EAGLE whole atmosphere model results

We present an analysis of the ionosphere and thermosphere response to Solar Proton Events (SPE) and magnetospheric proton precipitation in January 2005, which was carried out using the model of the entire atmosphere EAGLE. The ionization rates for the considered period were acquired from the AIMOS (Atmospheric Ionization Module Osnabrück) dataset. For numerical experiments, we applied only the proton-induced ionization rates of that period, while all the other model input parameters, including the electron precipitations, corresponded to the quiet conditions. In January 2005, two major solar proton events with different energy spectra and proton fluxes occurred on January 17 and January 20. Since two geomagnetic storms and several sub-storms took place during the considered period, not only solar protons but also less energetic magnetospheric protons contributed to the calculated ionization rates. Despite the relative transparency of the thermosphere for high-energy protons, an ionospheric response to the SPE and proton precipitation from the magnetotail was obtained in numerical experiments. In the ionospheric E layer, the maximum increase in the electron concentration is localized at high latitudes, and at heights of the ionospheric F2 layer, the positive perturbations were formed in the near-equatorial region. An analysis of the model-derived results showed that changes in the ionospheric F2 layer were caused by a change in the neutral composition of the thermosphere. We found that in the recovery phase after both solar proton events and the enhancement of magnetospheric proton precipitations associated with geomagnetic disturbances, the TEC and electron density in the F region and in topside ionosphere/plasmasphere increase at low- and mid-latitudes due to an enhancement of atomic oxygen concentration. Our results demonstrate an important role of magnetospheric protons in the formation of negative F-region ionospheric storms. According to our results, the topside ionosphere/plasmasphere and bottom-side ionosphere can react to solar and magnetospheric protons both with the same sign of disturbances or in different way. The same statement is true for TEC and foF2 disturbances. Different disturbances of foF2 and TEC at high and low latitudes can be explained by topside electron temperature disturbances.     

Variation of energetic particle precipitation with local magnetic time

The detailed study of the precipitation of magnetospheric particles into the atmosphere is complicated by the rather complex spatial configuration of the precipitation region and its variability with geomagnetic activity. In this paper we will introduce polar oval coordinates and apply them to POES observations of 30 keV to 2.5 MeV electrons and comparable protons to illustrate the dependence of particle precipitation on local time and geomagnetic activity. These coordinates also allow an easy separation of the spatial precipitation patterns of solar and magnetospheric particles. The results indicate that (a) the spatial precipitation pattern of energetic magnetospheric electrons basically follows the pattern of the field parallel Birkeland currents up to MeV energies and (b) at least in the mesosphere the influence of magnetospheric electrons is comparable to the one of solar electrons. Implications for modeling of atmospheric chemistry will be sketched.     

On the long-term variability of the heliosphere–magnetosphere environment

Annual means of measured and reconstructed solar, heliospheric, and magnetospheric parameters are used to infer solar activity signatures at the Hale and Gleissberg cycles timescales. Available open solar flux, modulation strength, cosmic ray flux, total solar irradiance data, reconstructed back to 1700, solar wind parameters (speed and density) and the magnitude of the heliospheric magnetic field at 1 AU, reconstructed back to 1870, as well as the time series of geomagnetic activity indices (aa, IDV, IHV), going back to 1870, have been considered. Simple filtering procedures (successive 11-, 22-, and 88-year running averages and differences between them) and scaling by the standard deviation from the average value for the common interval covered by the data show that the long-discussed variation in the 20th century (a pronounced increase since ∼1900, followed by a depression in the ‘60s and a new, slower, increase) seen in the 11-year averages of parameters such as geomagnetic activity indices and reconstructed heliospheric magnetic field strength, solar wind speed, open solar flux, is a result of the superposition in data of solar activity signatures at Hale and Gleissberg cycles timescales. The Hale and Gleissberg signals were characterized and similarities and differences in the temporal behavior of the analyzed parameters at these timescales are discussed. The similarities in the studied parameters point to a common pacing source, the solar dynamo.     

Mapping mean lake surface from satellite altimetry and GPS kinematic surveys

Lake water height is a key variable in water cycle and climate change studies, which is achievable using satellite altimetry constellation. A method based on data processing of altimetry from several satellites has been developed to interpolate mean lake surface (MLS) over a set of 22 big lakes distributed on the Earth. It has been applied on nadir radar altimeters in Low Resolution Mode (LRM: Jason-3, Saral/AltiKa, CryoSat-2) in Synthetic Aperture Radar (SAR) mode (Sentinel-3A), and in SAR interferometric (SARin) mode (CryoSat-2), and on laser altimetry (ICESat). Validation of the method has been performed using a set of kinematic GPS height profiles from 18 field campaigns over the lake Issykkul, by comparison of altimetry’s height at crossover points for the other lakes and using the laser altimetry on ICESat-2 mission. The precision reached ranges from 3 to 7 cm RMS (Root Mean Square) depending on the lakes. Currently, lake water level inferred from satellite altimetry is provided with respect to an ellipsoid. Ellipsoidal heights are converted into orthométric heights using geoid models interpolated along the satellite tracks. These global geoid models were inferred from geodetic satellite missions coupled with absolute and regional anomaly gravity data sets spread over the Earth. However, the spatial resolution of the current geoid models does not allow capturing short wavelength undulations that may reach decimeters in mountaineering regions or for rift lakes (Baikal, Issykkul, Malawi, Tanganika). We interpolate in this work the geoid height anomalies with three recent geoid models, the EGM2008, XGM2016 and EIGEN-6C4d, and compare them with the Mean Surface of 22 lakes calculated using satellite altimetry. Assuming that MLS mimics the local undulations of the geoid, our study shows that over a large set of lakes (in East Africa, Andean mountain and Central Asia), short wavelength undulations of the geoid in poorly sampled areas can be derived using satellite altimetry. The models used in this study present very similar geographical patterns when compared to MLS. The precision of the models largely depends on the location of the lakes and is about 18 cm, in average over the Earth. MLS can serve as a validation dataset for any future geoid model. It will also be useful for validation of the future mission SWOT (Surface Water and Ocean Topography) which will measure and map water heights over the lakes with a high horizontal resolution of 250 by 250 m.     

Emission heights of coronal bright points on Fe XII radiance map

The study of coronal bright points (BPs) is important for understanding coronal heating and the origin of the solar wind. Previous studies indicated that coronal BPs have a highly significant tendency to coincide with magnetic neutral lines in the photosphere. Here we further studied the emission heights of the BPs above the photosphere in the bipolar magnetic loops that are apparently associated with them. As BPs are seen in projection against the disk their true emission heights are unknown. The correlation of the BP locations on the Fe XII radiance map from EIT with the magnetic field features (in particular neutral lines) was investigated in detail. The coronal magnetic field was determined by an extrapolation of the photospheric field (derived from 2-D magnetograms obtained from the Kitt Peak observatory) to different altitudes above the disk. It was found that most BPs sit on or near a photospheric neutral line, but that the emission occurs at a height of about 5 Mm. Some BPs, while being seen in projection, still seem to coincide with neutral lines, although their emission takes place at heights of more than 10 Mm. Such coincidences almost disappear for emissions above 20 Mm. We also projected the upper segments of the 3-D magnetic field lines above different heights, respectively, on to the tangent x–y plane, where x is in the east–west and y in the south–north direction. The shape of each BP was compared with the respective field-line segment nearby. This comparison suggests that most coronal BPs are actually located on the top of their associated magnetic loops. Finally, we calculated for each selected BP region the correlation coefficient between the Fe XII intensity enhancement and the horizontal component of the extrapolated magnetic field vector at the same x–y position in planes of different heights, respectively. We found that for almost all the BP regions we studied the correlation coefficient, with increasing height, increases to a maximal value and then decreases again. The height corresponding to this maximum was defined as the correlation height, which for most bright points was found to range below 20 Mm.     

SAMBO-GEOS: On three-dimensional substorm dynamics — A case study for 4 March 1979

For investigations of the isolated magnetospheric substorm the ground-based and satellite geophysical data are discussed. The active phase of the substorm has two parts: an active-convective phase and a classical active phase. The analysed data show that the breakup of the substorm takes place at closed and not very much stretched magnetic field lines.     

Magnetospheric electric field variations caused by storm-time shock fronts

On January 20, 2005 there was an X 7.1 solar flare at 0636 UT with an accompanied halo coronal mass ejection (CME). The resultant interplanetary shock impacted earth ∼36 h later. Near earth, the Advanced Composition Explorer (ACE) spacecraft observed two impulses with a staircase structure in density and pressure. The estimated earth-arrival times of these impulses were 1713 UT and 1845 UT on January 21, 2005. Three MINIature Spectrometer (MINIS) balloons were aloft on January 21st; one in the northern polar stratosphere and two in the southern polar stratosphere. MeV relativistic electron precipitation (REP) observed by all three balloons is coincident (<3 min) with the impulse arrivals and magnetospheric compression observed by both GOES 10 and 12. Balloon electric field data from the southern hemisphere show no signs of the impulse electric field directly reaching the ionosphere. Enhancement of the balloon-observed convection electric field by as much as 40 mV/m in less than 20 min during this time period is consistent with typical substorm growth. Precipitation-induced ionospheric conductivity enhancements are suggested to be (a) the result of both shock arrival and substorm activity and (b) the cause of rapid (<6 min) decreases in the observed electric field (by as much as 40 mV/m). There is poor agreement between peak cross polar cap potential in the northern hemisphere calculated from Super Dual Auroral Radar Network (SuperDARN) echoes and horizontal electric field at the MINIS balloon locations in the southern hemisphere. Possible reasons for this poor agreement include (a) a true lack of north–south conjugacy between measurement sites, (b) an invalid comparison between global (SuperDARN radar) and local (MINIS balloon) measurements and/or (c) radar absorption resulting from precipitation-induced D-region ionosphere density enhancements.     

Simultaneous observations of the plasma density on the same field line by the CPMN ground magnetometers and the Cluster satellites

By applying the cross-phase method and the amplitude-ratio method to magnetic field data obtained from two ground stations located close to each other, we can determine the frequency of the field line resonance (FLR), or the field line eigenfrequency, for the field line running through the midpoint of the two stations. From thus identified FLR frequency we can estimate the equatorial plasma mass density (ρ)$(\rho )$ by using the T05s magnetospheric field model [Tsyganenko, N.A., Sitnov, M.I. Modeling the dynamics of the inner magnetosphere during strong geomagnetic storms, J. Geophys. Res. 110, A03208, 2005] and the equation of Singer et al. [Singer, H.J., Southwood, D.J., Walker, R.J., Kivelson, M.G. Alfven wave resonances in a realistic magnetospheric magnetic field geometry, J. Geophys. Res. 86 (A6) 4589–4596, 1981].