ALBERT

Description: Magnetic field and dynamic pressure ULF fluctuations in coronal-mass-ejection-driven sheath regions Annales Geophysicae, 31, 1559-1567, 2013 Author(s): E. K. J. Kilpua, H. Hietala, H. E. J. Koskinen, D. Fontaine, and L. Turc Compressed sheath regions form ahead of interplanetary coronal mass ejections (ICMEs) that are sufficiently faster than the preceding solar wind. The turbulent sheath regions are important drivers of magnetospheric activity, but due to their complex internal structure, relatively little is known on the distribution of the magnetic field and plasma variations in them. In this paper we investigate ultra low frequency (ULF) fluctuations in the interplanetary magnetic field (IMF) and in dynamic pressure ( P dyn ) using a superposed epoch analysis of 41 sheath regions observed during solar cycle 23. We find strongest fluctuation power near the shock and in the vicinity of the ICME leading edge. The IMF and P dyn ULF power have different profiles within the sheath; the former is enhanced in the leading part of the sheath, while the latter is increased in the trailing part of the sheath. We also find that the ICME properties affect the level and distribution of the ULF power in sheath regions. For example, sheath regions associated with strong or fast ICMEs, or those that are crossed at intermediate distances from the center, have strongest ULF power and large variation in the power throughout the sheath region. The weaker or slower ICMEs, or those that are crossed centrally, have in general considerably weaker ULF power with relatively smooth profiles. The strong and abrupt decrease of the IMF ULF power at the ICME leading edge could be used to distinguish the ICME from the preceding sheath plasma.

Description: Cluster observation of few-hour-scale evolution of structured plasma in the inner magnetosphere Annales Geophysicae, 31, 1569-1578, 2013 Author(s): M. Yamauchi, I. Dandouras, H. Rème, R. Lundin, and L. M. Kistler Using Cluster Ion Spectrometry (CIS) data from the spacecraft-4 perigee traversals during the 2001–2006 period (nearly 500 traversals after removing those that are highly contaminated by radiation belt particles), we statistically examined the local time distribution of structured trapped ions at sub- to few-keV range as well as inbound–outbound differences of these ion signatures in intensities and energy–latitude dispersion directions. Since the Cluster orbit during this period was almost constant and approximately north–south symmetric at nearly constant local time near the perigee, inbound–outbound differences are attributed to temporal developments in a 1–2 h timescale. Three types of structured ions at sub- to few keV range that are commonly found in the inner magnetosphere are examined: – Energy–latitude dispersed structured ions at less than a few keV, – Short-lived dispersionless ion stripes at wide energy range extending 0.1–10 keV, – Short-lived low-energy ion bursts at less than a few hundred eV. The statistics revealed that the wedge-like dispersed ions are most often observed in the dawn sector (60% of traversals), and a large portion of them show significant enhancement during the traversals at all local times. The short-lived ion stripes are predominantly found near midnight, where most stripes are significantly enhanced during the traversals and are associated with substorm activities with geomagnetic AL 〈 −300 nT. The low-energy bursts are observed at all local times and under all geomagnetic conditions, with moderate peak of the occurrence rate in the afternoon sector. A large portion of them again show significant enhancement or decay during the traversals.

Description: Influence of the terrestrial magnetic field geometry on the cutoff rigidity of cosmic ray particles Annales Geophysicae, 31, 1637-1643, 2013 Author(s): K. Herbst, A. Kopp, and B. Heber Studies of the propagation of charged energetic particles in the Earth's magnetic field go back to Carl Størmer. In the end, his investigations finally lead to the definition of the so-called cutoff rigidity R C ; that is, the minimum momentum per charge a particle must have in order to reach a certain geographical location. Employing Monte Carlo simulations with the PLANETOCOSMICS code we investigate the correlation between the geomagnetic field structure and the cutoff rigidity. We show that the geometry of the magnetic field has a considerable influence on the resulting cutoff rigidity distribution. Furthermore, we will present a simple geometry-based parameter, δ B , which is able to reflect the location-dependent cutoff rigidity. We show that this correlation is also visible in the temporal evolution of the Earth's magnetic field, at least over the last 100 yr. Using latitude scans with neutron monitors, changes of the relative counting rates at different positions are calculated, showing small variations for, e.g., Kiel and Moscow, while large ones occur at Mexico City as well as on the British Virgin Islands.

Description: Electron acceleration at Jupiter: input from cyclotron-resonant interaction with whistler-mode chorus waves Annales Geophysicae, 31, 1619-1630, 2013 Author(s): E. E. Woodfield, R. B. Horne, S. A. Glauert, J. D. Menietti, and Y. Y. Shprits Jupiter has the most intense radiation belts of all the outer planets. It is not yet known how electrons can be accelerated to energies of 10 MeV or more. It has been suggested that cyclotron-resonant wave-particle interactions by chorus waves could accelerate electrons to a few MeV near the orbit of Io. Here we use the chorus wave intensities observed by the Galileo spacecraft to calculate the changes in electron flux as a result of pitch angle and energy diffusion. We show that, when the bandwidth of the waves and its variation with L are taken into account, pitch angle and energy diffusion due to chorus waves is a factor of 8 larger at L -shells greater than 10 than previously shown. We have used the latitudinal wave intensity profile from Galileo data to model the time evolution of the electron flux using the British Antarctic Survey Radiation Belt (BAS) model. This profile confines intense chorus waves near the magnetic equator with a peak intensity at ∼5° latitude. Electron fluxes in the BAS model increase by an order of magnitude for energies around 3 MeV. Extending our results to L = 14 shows that cyclotron-resonant interactions with chorus waves are equally important for electron acceleration beyond L = 10. These results suggest that there is significant electron acceleration by cyclotron-resonant interactions at Jupiter contributing to the creation of Jupiter's radiation belts and also increasing the range of L -shells over which this mechanism should be considered.

Description: The coupling between the solar wind and proton fluxes at GEO Annales Geophysicae, 31, 1631-1636, 2013 Author(s): R. J. Boynton, S. A. Billings, O. A. Amariutei, and I. Moiseenko The relationship between the solar wind and the proton flux at geosynchronous Earth orbit (GEO) is investigated using the error reduction ratio (ERR) analysis. The ERR analysis is able to search for the most appropriate inputs that control the evolution of the system. This approach is a black box method and is able to derive a mathematical model of a system from input-output data. This method is used to analyse eight energy ranges of the proton flux at GEO from 80 keV to 14.5 MeV. The inputs to the algorithm were solar wind velocity, density and pressure; the Dst index; the solar energetic proton (SEP) flux; and a function of the interplanetary magnetic field (IMF) tangential magnitude and clock angle. The results show that for lowest five energy channels (80 to 800 keV) the GEO proton fluxes are controlled by the solar wind velocity with a lag of two to three days. However, above 350 keV, the SEP fluxes, accounts for a significant portion of the GEO proton flux variance. For the highest three energy channels (0.74 to 14.5 MeV), the SEPs account for the majority of the ERR. The results also show an anisotropy of protons with gyrocenters inside GEO and outside GEO, where the protons inside GEO are controlled partly by the Dst index and also an IMF-clock angle function.

Description: A statistical study on O + flux in the dayside magnetosheath Annales Geophysicae, 31, 1005-1010, 2013 Author(s): R. Slapak, H. Nilsson, and L. G. Westerberg Studies on terrestrial oxygen ion (O + ) escape into the interplanetary space have considered a number of different escape paths. Recent observations however suggest a yet insufficiently investigated additional escape route for hot O + : along open magnetic field lines in the high altitude cusp and mantle. Here we present a statistical study on O + flux in the high-latitude dayside magnetosheath. The O + is generally seen relatively close to the magnetopause, consistent with observations of O + flowing primarily tangentially to the magnetopause. We estimate the total escape flux in this region to be ~ 7 × 10 24 s −1 , implying this escape route to significantly contribute to the overall total O + loss into interplanetary space.

Description: Simulation of the influence of historical land cover changes on the global climate Annales Geophysicae, 31, 995-1004, 2013 Author(s): Y. Wang, X. Yan, and Z. Wang In order to estimate biogeophysical effects of historical land cover change on climate during last three centuries, a set of experiments with a climate system model of intermediate complexity (MPM-2) is performed. In response to historical deforestation, the model simulates a decrease in annual mean global temperature in the range of 0.07–0.14 °C based on different grassland albedos. The effect of land cover changes is most pronounced in the middle northern latitudes with maximum cooling reaching approximately 0.6 °C during northern summer. The cooling reaches 0.57 °C during northern spring owing to the large effects of land surface albedo. These results suggest that land cover forcing is important for study on historical climate change and that more research is necessary in the assessment of land management options for climate change mitigation.

Description: IMF effect on the polar cap contraction and expansion during a period of substorms Annales Geophysicae, 31, 1021-1034, 2013 Author(s): A. T. Aikio, T. Pitkänen, I. Honkonen, M. Palmroth, and O. Amm The polar cap boundary (PCB) location and motion in the nightside ionosphere has been studied by using measurements from the EISCAT radars and the MIRACLE magnetometers during a period of four substorms on 18 February 2004. The OMNI database has been used for observations of the solar wind and the Geotail satellite for magnetospheric measurements. In addition, the event was modelled by the GUMICS-4 MHD simulation. The simulation of the PCB location was in a rather good agreement with the experimental estimates at the EISCAT longitude. During the first three substorm expansion phases, neither the local observations nor the global simulation showed any poleward motions of the PCB, even though the electrojets intensified. Rapid poleward motions of the PCB took place only in the early recovery phases of the substorms. Hence, in these cases the nightside reconnection rate was locally higher in the recovery phase than in the expansion phase. In addition, we suggest that the IMF B z component correlated with the nightside tail inclination angle and the PCB location with about a 17-min delay from the bow shock. By taking the delay into account, the IMF northward turnings were associated with dipolarizations of the magnetotail and poleward motions of the PCB in the recovery phase. The mechanism behind this effect should be studied further.

Description: A comparison of bow shock models with Cluster observations during low Alfvén Mach number magnetic clouds Annales Geophysicae, 31, 1011-1019, 2013 Author(s): L. Turc, D. Fontaine, P. Savoini, H. Hietala, and E. K. J. Kilpua Magnetic clouds (MCs) are very geoeffective solar wind structures. Their properties in the interplanetary medium have been extensively studied, yet little is known about their characteristics in the Earth's magnetosheath. The Cluster spacecraft offer the opportunity to observe MCs in the magnetosheath, but before MCs reach the magnetosphere, their structure is altered when they interact with the terrestrial bow shock (BS). The physics taking place at the BS strongly depends on Θ Bn , the angle between the shock normal and the interplanetary magnetic field. However, in situ observations of the BS during an MC's crossing are seldom available. In order to relate magnetosheath observations to solar wind conditions, we need to rely on a model to determine the shock's position and normal direction. Yet during MCs, the models tend to be less accurate, because the Alfvén Mach number ( M A ) is often significantly lower than in regular solar wind. On the contrary, the models are generally optimised for high M A conditions. In this study, we compare the predictions of four widely used models available in the literature (Wu et al., 2000; Chapman and Cairns, 2003; Jeřáb et al., 2005; Měrka et al., 2005b) to Cluster's dayside BS crossings observed during five MC events. Our analysis shows that the Θ Bn angle is well predicted by all four models. On the other hand, the Jeřáb et al. (2005) model yields the best estimates of the BS position during low M A MCs. The other models locate the BS either too far from or too close to Earth. The results of this paper can be directly used to estimate the BS parameters in all studies of MC interaction with Earth's magnetosphere.

Description: Investigation of radiative effects of the optically thick dust layer over the Indian tropical region Annales Geophysicae, 31, 647-663, 2013 Author(s): S. K. Das, J.-P. Chen, M. Venkat Ratnam, and A. Jayaraman Optical and physical properties of aerosols derived from multi-satellite observations (MODIS-Aqua, OMI-Aura, MISR-Terra, CALIOP-CALIPSO) have been used to estimate radiative effects of the dust layer over southern India. The vertical distribution of aerosol radiative forcing and heating rates are calculated with 100 m resolution in the lower atmosphere, using temperature and relative humidity data from balloon-borne radiosonde observations. The present study investigates the optically thick dust layer of optical thickness 0.18 ± 0.06 at an altitude of 2.5 ± 0.7 km over Gadanki, transported from the Thar Desert, producing radiative forcing and heating rate of 11.5 ± 3.3 W m −2 and 0.6 ± 0.26 K day −1 , respectively, with a forcing efficiency of 43 W m −2 and an effective heating rate of 4 K day −1 per unit dust optical depth. Presence of the dust layer increases radiative forcing by 60% and heating rate by 60 times at that altitude compared to non-dusty cloud-free days. Calculation shows that the radiative effects of the dust layer strongly depend on the boundary layer aerosol type and mass loading. An increase of 25% of heating by the dust layer is found over relatively cleaner regions than urban regions in southern India and further 15% of heating increases over the marine region. Such heating differences in free troposphere may have significant consequences in the atmospheric circulation and hydrological cycle over the tropical Indian region.

Description: Low level jet intensification by mineral dust aerosols Annales Geophysicae, 31, 625-632, 2013 Author(s): O. Alizadeh Choobari, P. Zawar-Reza, and A. Sturman Modification of the intensity of a low level jet (LLJ) and near-surface wind speed by mineral dust is important as it has implications for dust emission and its long-range transport. Using the Weather Research and Forecasting with Chemistry (WRF/Chem) regional model, it is shown that direct radiative forcing by mineral dust reduces temperature in the lower atmosphere, but increases it in the layers aloft. The surface cooling is shown to be associated with a reduction of turbulent kinetic energy (TKE) and hence vertical mixing of horizontal momentum. Changes in the vertical profile of temperature over the regions that are under the influence of a LLJ are shown to result in an intensification of the LLJ and near-surface wind speed, but a decrease of winds aloft. These changes in the wind speed profile differ from results of previous research which suggested a decrease of wind speed in the lower atmosphere and its increase in the upper boundary layer.

Description: Parametric validations of analytical lifetime estimates for radiation belt electron diffusion by whistler waves Annales Geophysicae, 31, 599-624, 2013 Author(s): A. V. Artemyev, D. Mourenas, O. V. Agapitov, and V. V. Krasnoselskikh The lifetimes of electrons trapped in Earth's radiation belts can be calculated from quasi-linear pitch-angle diffusion by whistler-mode waves, provided that their frequency spectrum is broad enough and/or their average amplitude is not too large. Extensive comparisons between improved analytical lifetime estimates and full numerical calculations have been performed in a broad parameter range representative of a large part of the magnetosphere from L ~ 2 to 6. The effects of observed very oblique whistler waves are taken into account in both numerical and analytical calculations. Analytical lifetimes (and pitch-angle diffusion coefficients) are found to be in good agreement with full numerical calculations based on CRRES and Cluster hiss and lightning-generated wave measurements inside the plasmasphere and Cluster lower-band chorus waves measurements in the outer belt for electron energies ranging from 100 keV to 5 MeV. Comparisons with lifetimes recently obtained from electron flux measurements on SAMPEX, SCATHA, SAC-C and DEMETER also show reasonable agreement.

Description: Coherent amplitude modulation of electron-beam-driven Langmuir waves Annales Geophysicae, 31, 633-638, 2013 Author(s): K. Baumgärtel A linear approach to the phenomenon of irregular amplitude modulation of beam-driven Langmuir waves, developed in a previous paper, is extended to explain periodic modulation as well. It comes about by beating of the fastest growing mode of the instability with beam-aligned plasma oscillations. They are naturally generated in a uniform domain of beam–plasma interaction prior to the onset of the instability. Particle-in-cell (PIC) simulations support the results of the linear analysis.

Description: Seasonal variation of the ion upflow in the topside ionosphere during SAPS (subauroral polarization stream) periods Annales Geophysicae, 31, 1521-1534, 2013 Author(s): H. Wang and H. Lühr A statistical study has been performed by using two years of DMSP (Defense Meteorological Satellite Program) plasma observations to investigate the seasonal effect of SAPS (subauroral polarization stream) on the ion upflow in the duskside ionosphere of the Northern Hemisphere. There are obvious upflows occurring in the topside ionosphere around the SAPS region, exceeding 200 m s −1 at winter solstice, indicating an important relationship between SAPS and the local plasma upward motion. Both SAPS and ion upward velocities show similar seasonal variations, largest in winter and smallest in summer, irrespective of geomagnetic activity. A good correlation is found and a linear relationship is derived between SAPS and the ion upflow velocities. During December solstice the average upflow flux can reach about 2 × 10 8 cm −2 s −1 for more disturbed periods, which is comparable to the typical upflow flux in the dayside cusp region. The depression of the ion temperatures around the peak SAPS region can be understood in terms of the adiabatic cooling. The hot ion cools down when expanding into the low ion concentration region. The electron temperature elevates around the SAPS region because of the reduced Coulomb cooling in the low ion density region. Both the changes of ion and electron temperatures are larger in winter than in summer, however, for Kp 〈 4 the electron temperatures are almost seasonably independent. The present work highlights the important role of the SAPS-related frictional heating at mid-latitudes on the local formation of the strong upward flow, which might provide a direct ionospheric ion source for the ring current and plasmasphere in the duskside sector.

Description: Magnetic field generation in a jet-sheath plasma via the kinetic Kelvin-Helmholtz instability Annales Geophysicae, 31, 1535-1541, 2013 Author(s): K.-I. Nishikawa, P. Hardee, B. Zhang, I. Duţan, M. Medvedev, E. J. Choi, K. W. Min, J. Niemiec, Y. Mizuno, A. Nordlund, J. T. Frederiksen, H. Sol, M. Pohl, and D. H. Hartmann We have investigated the generation of magnetic fields associated with velocity shear between an unmagnetized relativistic jet and an unmagnetized sheath plasma. We have examined the strong magnetic fields generated by kinetic shear (Kelvin–Helmholtz) instabilities. Compared to the previous studies using counter-streaming performed by Alves et al. (2012), the structure of the kinetic Kelvin–Helmholtz instability (KKHI) of our jet-sheath configuration is slightly different, even for the global evolution of the strong transverse magnetic field. In our simulations the major components of growing modes are the electric field E z , perpendicular to the flow boundary, and the magnetic field B y , transverse to the flow direction. After the B y component is excited, an induced electric field E x , parallel to the flow direction, becomes significant. However, other field components remain small. We find that the structure and growth rate of KKHI with mass ratios m i / m e = 1836 and m i / m e = 20 are similar. In our simulations saturation in the nonlinear stage is not as clear as in counter-streaming cases. The growth rate for a mildly-relativistic jet case (γ j = 1.5) is larger than for a relativistic jet case (γ j = 15).

Description: Ionospheric response to total solar eclipse of 22 July 2009 in different Indian regions Annales Geophysicae, 31, 1549-1558, 2013 Author(s): S. Kumar, A. K. Singh, and R. P. Singh The variability of ionospheric response to the total solar eclipse of 22 July 2009 has been studied analyzing the GPS data recorded at the four Indian low-latitude stations Varanasi (100% obscuration), Kanpur (95% obscuration), Hyderabad (84% obscuration) and Bangalore (72% obscuration). The retrieved ionospheric vertical total electron content (VTEC) shows a significant reduction (reflected by all PRNs (satellites) at all stations) with a maximum of 48% at Varanasi (PRN 14), which decreases to 30% at Bangalore (PRN 14). Data from PRN 31 show a maximum of 54% at Kanpur and 26% at Hyderabad. The maximum decrement in VTEC occurs some time (2–15 min) after the maximum obscuration. The reduction in VTEC compared to the quiet mean VTEC depends on latitude as well as longitude, which also depends on the location of the satellite with respect to the solar eclipse path. The amount of reduction in VTEC decreases as the present obscuration decreases, which is directly related to the electron production by the photoionization process. The analysis of electron density height profile derived from the COSMIC (Constellation Observing System for Meteorology, Ionosphere & Climate) satellite over the Indian region shows significant reduction from 100 km altitude up to 800 km altitude with a maximum of 48% at 360 km altitude. The oscillatory nature in total electron content data at all stations is observed with different wave periods lying between 40 and 120 min, which are attributed to gravity wave effects generated in the lower atmosphere during the total solar eclipse.

Description: Application of nonlinear autoregressive moving average exogenous input models to geospace: advances in understanding and space weather forecasts Annales Geophysicae, 31, 1579-1589, 2013 Author(s): R. J. Boynton, M. A. Balikhin, S. A. Billings, and O. A. Amariutei The nonlinear autoregressive moving average with exogenous inputs (NARMAX) system identification technique is applied to various aspects of the magnetospheres dynamics. It is shown, from an example system, how the inputs to a system can be found from the error reduction ratio (ERR) analysis, a key concept of the NARMAX approach. The application of the NARMAX approach to the Dst (disturbance storm time) index and the electron fluxes at geostationary Earth orbit (GEO) are reviewed, revealing new insight into the physics of the system. The review of studies into the Dst index illustrate how the NARMAX approach is able to find a coupling function for the Dst index from data, which was then analytically justified from first principles. While the review of the electron flux demonstrates how NARMAX is able to reveal new insight into the physics of the acceleration and loss processes within the radiation belt.

Description: Energy exchange and wave action conservation for magnetohydrodynamic (MHD) waves in a general, slowly varying medium Annales Geophysicae, 32, 1495-1510, 2014 Author(s): A. D. M. Walker Magnetohydrodynamic (MHD) waves in the solar wind and magnetosphere are propagated in a medium whose velocity is comparable to or greater than the wave velocity and which varies in both space and time. In the approximation where the scales of the time and space variation are long compared with the period and wavelength, the ray-tracing equations can be generalized and then include an additional first-order differential equation that determines the variation of frequency. In such circumstances the wave can exchange energy with the background: wave energy is not conserved. In such processes the wave action theorem shows that the wave action, defined as the ratio of the wave energy to the frequency in the local rest frame, is conserved. In this paper we discuss ray-tracing techniques and the energy exchange relation for MHD waves. We then provide a unified account of how to deal with energy transport by MHD waves in non-uniform media. The wave action theorem is derived directly from the basic MHD equations for sound waves, transverse Alfvén waves, and the fast and slow magnetosonic waves. The techniques described are applied to a number of illustrative cases. These include a sound wave in a medium undergoing a uniform compression, an isotropic Alfvén wave in a steady-state shear layer, and a transverse Alfvén wave in a simple model of the magnetotail undergoing compression. In each case the nature and magnitude of the energy exchange between wave and background is found.

Description: On the influence of CMEs on the global 3-D coronal electron density Annales Geophysicae, 29, 1019-1028, 2011 Author(s): M. Kramar, J. Davila, H. Xie, and S. Antiochos In order to analize the influence of a Coronal Mass Ejection (CME) on the coronal streamer belt, we made 3-D reconstructions of the electron density in the corona at heliospheric distances from 1.5 to 4 R ⊙ for periods before and after a CME occured. The reconstructions were performed using a tomography technique. We studied two CME cases: (i) a slow CME on 1 June 2008; (ii) two fast CMEs on 31 December 2007 and 2 January 2008. For the first case of slow CME, it was found: (i) the potential magnetic field configuration in the CME initiation region before the CME does not agree with the coronal density structure while after the CME the agreement between the field and density is much better. This could be manifistation of that that the field was non-potential before the CME and after the CME the field relaxes towards a more potential state. (ii) It was shown that the dimming caused by the slow CME is not due to rotation of the corona and a line-of-sight (LOS) effect but a streamer blow out effect took place.

Description: Longitudinal study of the ionospheric response to the geomagnetic storm of 15 May 2005 and manifestation of TADs Annales Geophysicae, 29, 1063-1070, 2011 Author(s): S. Sharma, P. Galav, N. Dashora, and R. Pandey Response of low latitude ionosphere to the geomagnetic storm of 15 May 2005 has been studied using total electron content (TEC) data, obtained from three GPS stations namely, Yibal, Udaipur and Kunming situated near the northern crest of equatorial ionization anomaly at different longitudes. Solar wind parameters, north-south component of the interplanetary magnetic field (IMF B z ) and AE index data have been used to infer the strength of the geomagnetic storm. A large value of eastward interplanetary electric field at 06:15 UT, during the time of maximum southward IMF B z has been used to infer the transmission of an eastward prompt penetration electric field (PPEF) which resulted in a peak in TEC at 07:45 UT due to the local uplift of plasma in the low latitudes near the anomaly crest over a wide range of longitudes. Wave-like modulations superposed over the second enhancement in TEC between 09:15 UT to 10:30 UT have been observed at all the three stations. The second enhancement in TEC along with the modulations of up to 5 TECU have been attributed to the combined effect of super plasma fountain and traveling atmospheric disturbances (TAD). Observed large enhancements in TEC are a cause of concern for satellite based navigation and ground positioning. Increased [O/N 2 ] ratio between 09:15 UT to 10:15 UT when modulations in TEC have been also observed, confirms the presence of TADs over a wide range of longitudes.

Description: Volume cross section of auroral radar backscatter and RMS plasma fluctuations inferred from coherent and incoherent scatter data: a response on backscatter volume parameters Annales Geophysicae, 29, 1081-1092, 2011 Author(s): M. V. Uspensky, P. Janhunen, A. V. Koustov, and K. Kauristie Norway and Finland STARE radar measurements in the eastward auroral electrojet are combined with EISCAT CP-1 measurements of the electron density and electric field vector in the common scattering volume to investigate the variation of the auroral radar volume cross section (VCS) with the flow angle of observations (radar look direction with respect to the E × B electron drift). The data set available consists of ~6000 points for flow angles of 40–85° and electron drifts between 500 and 2000 m s −1 . The EISCAT electron density N(h) -profile data are used to estimate the effective electron density, aspect angle and thickness of the backscattering layer. It is shown that the flow angle variation of the VCS is rather weak, only ~5 dB within the range of the considered flow angles. The VCS values themselves respond almost linearly to the square of both the electron drift velocity magnitude and the effective electron density. By adopting the inferred shape of the VCS variation with the flow angle and the VCS dependence upon wavelength, the relative amplitude of electrostatic electron density fluctuations over all scales is estimated. Inferred values of 2–4 percent react nearly linearly to the electron drift velocity in the range of 500–1000 m s −1 but the rate of increase slows down at electron drifts 〉1000 m s −1 and density fluctuations of ~5.5 percent due to, perhaps, progressively growing nonlinear wave losses.

Description: Flux quanta, magnetic field lines, merging – some sub-microscale relations of interest in space plasma physics Annales Geophysicae, 29, 1121-1127, 2011 Author(s): R. A. Treumann, R. Nakamura, and W. Baumjohann We clarify the notion of magnetic field lines in plasma by referring to sub-microscale (quantum mechanical) particle dynamics. It is demonstrated that magnetic field lines in a field of strength B carry single magnetic flux quanta Φ 0 = h/e . The radius of a field line in the given magnetic field B is calculated. It is shown that such field lines can merge and annihilate only over the length ℓ ∥ of their strictly anti-parallel sections, for which case we estimate the power generated. The length ℓ ∥ becomes a function of the inclination angle θ of the two merging magnetic flux tubes (field lines). Merging is possible only in the interval 1 2 πθ≤π. This provides a sub-microscopic basis for "component reconnection" in classical macro-scale reconnection. We also find that the magnetic diffusion coefficient in plasma appears in quanta D 0 m = e Φ 0 / m e = h/m e . This lets us conclude that the bulk perpendicular plasma resistivity is limited and cannot be less than η 0⊥ =μ 0 e Φ 0 / m e =μ 0 h/m e ~10 −9 Ohm m. This resistance is an invariant.

Description: Could periodic patterns in human mortality be sensitive to solar activity? Annales Geophysicae, 29, 1113-1120, 2011 Author(s): R. Díaz-Sandoval, R. Erdélyi, and R. Maheswaran Seasonal behaviour of human diseases have been observed and reported in the literature for years. Although the Sun plays an essential role in the origin and evolution of life on Earth, it is barely taken into account in biological processes for the development of a specific disease. Higher mortality rates occur during the winter season in the Northern Hemisphere for several diseases, particularly diseases of the cardiovascular and respiratory systems. This increment has been associated with seasonal and social causes. However, is there more behind these correlations, in particular in terms of solar variability? In this paper we attempt to make a first step towards answering this question. A detailed wavelet analysis of periodicities for diseases from England and Wales seem to reveal that mortality periodicities (3 days to half a year) could be due to the Earth's position around the Sun. Moreover, crosswavelet and wavelet coherence analysis show common features between medical diseases and solar proxies around solar maximum activity suggesting that this relation, if any, has to be searched in times of high solar activity.

Description: Excitation of zero-frequency magnetic field-aligned currents by ionospheric heating Annales Geophysicae, 29, 1147-1152, 2011 Author(s): A. V. Streltsov and T. R. Pedersen Time-dependent, three-dimensional numerical simulations of the reduced MHD model describing shear Alfvén waves in the magnetosphere provide an interesting prediction superficially similar to results of several ionospheric heating experiments conducted at high altitudes. In these experiments, heating of the ionospheric F-region with a constant/zero-frequency beam of HF waves causes luminous structures in the ionosphere in the form of a ring or a solid spot with a characteristic size comparable to the size of the heated spot. Simulations suggest that spots/rings or similar optical appearance might be associated with a magnetic field-aligned current system produced by the ionospheric heating. Two of the most interesting features of this current system are (1) strong localization across the ambient magnetic field and (2) distinctive non-symmetrical luminous signatures (ring/spot) in magnetically conjugate locations in the ionosphere.

Description: Comment on "The PC index: review of methods", by McCreadie and Menvielle (2010) Annales Geophysicae, 29, 1137-1146, 2011 Author(s): P. Stauning The Polar Cap (PC) index is a controversial topic within the IAGA scientific community. Since 1997 discussions of the validity of the index to be endorsed as an official IAGA index have ensued. The article: "The PC index: review of methods", written by two members of the IAGA PC index committee, H. McCreadie and M. Menvielle, holds a critical review of some aspects of the methods used to derive PC index values. However, a number of incorrect statements and factual errors have been found and shall be called attention to and discussed in this commentary. Further critical comments concern the Corrigendum issued by the same authors and published in Ann. Geophys., 29, 813–814, 2011.

Description: The role of the bow shock in solar wind-magnetosphere coupling Annales Geophysicae, 29, 1129-1135, 2011 Author(s): R. E. Lopez, V. G. Merkin, and J. G. Lyon In this paper we examine the role of the bow shock in coupling solar wind energy to the magnetosphere using global magnetohydrodynamic simulations of the solar wind-magnetosphere interaction with southward IMF. During typical solar wind conditions, there are two significant dynamo currents in the magnetospheric system, one in the high-latitude mantle region tailward of the cusp and the other in the bow shock. As the magnitude of the (southward) IMF increases and the solar wind becomes a low Mach number flow, there is a significant change in solar wind-magnetosphere coupling. The high-latitude magnetopause dynamo becomes insignificant compared to the bow shock and a large load appears right outside the magnetopause. This leaves the bow shock current as the only substantial dynamo current in the system, and the only place where a significant amount of mechanical energy is extracted from the solar wind. That energy appears primarily as electromagnetic energy, and the Poynting flux generated at the bow shock feeds energy back into the plasma, reaccelerating it to solar wind speeds. Some small fraction of that Poynting flux is directed into the magnetosphere, supplying the energy needed for magnetospheric dynamics. Thus during periods when the solar wind flow has a low Mach number, the main dynamo in the solar wind-magnetosphere system is the bow shock.

Description: Measurements of aerosol intensive properties over Visakhapatnam, India for 2007 Annales Geophysicae, 29, 973-985, 2011 Author(s): K. Niranjan, B. Spandana, T. Anjana Devi, V. Sreekanth, and B. L. Madhavan Information on the aerosol intensive properties like Single Scattering Albedo (SSA) and asymmetry parameter are very limited, particularly over the peninsular India, though extensive reports are available on the aerosol bulk properties. In view of the importance of these parameters in evaluating the aerosol radiative forcing, we present for the first time the temporal variation in SSA with measurements on aerosol absorption and scattering coefficients over Visakhapatnam (17.72° N, 83.32° E; located on the east coast of India) for the year 2007. The inferred SSA was in the range of 0.65 and 0.9 with an annual mean of 0.76 ± 0.013 and with a probable value of 0.80, indicating a marginal atmospheric warming over the region. The mixed layer contribution to column Aerosol Optical depth is found to be 35 % in summer while it is well above 35 % in winter, indicating the confinement of aerosol within the boundary layer during winter. The asymmetry parameter which represents the angular scattering in radiative forcing estimation is found to be around 0.65 ± 0.1 for the location. The day to day variability in SSA is found to be well correlated with the variations in surface BC mass concentrations and/or the relative dominance of the fine/coarse mode aerosol. The results are discussed in light of the aerosol physical and optical properties and the asymmetry parameter.

Description: Corrigendum to "Downward auroral currents from the reconnection Hall-region", published in Ann. Geophys., 29, 679–685, 2011 Annales Geophysicae, 29, 1061-1061, 2011 Author(s): R. A. Treumann, R. Nakamura, and W. Baumjohann No abstract available.

Description: Hybrid models of solar wind plasma heating Annales Geophysicae, 29, 1071-1079, 2011 Author(s): L. Ofman, A.-F. Viñas, and P. S. Moya Remote sensing and in-situ observations show that solar wind ions are often hotter than electrons, and the heavy ions flow faster than the protons by up to an Alfvén speed. Turbulent spectrum of Alfvénic fluctuations and shocks were detected in solar wind plasma. Cross-field inhomogeneities in the corona were observed to extend to several tens of solar radii from the Sun. The acceleration and heating of solar wind plasma is studied via 1-D and 2-D hybrid simulations. The models describe the kinetics of protons and heavy ions, and electrons are treated as neutralizing fluid.The expansion of the solar wind is considered in 1-D hybrid model. A spectrum of Alfvénic fluctuations is injected at the computational boundary, produced by differential streaming instability, or initial ion temperature anisotropy, and the parametric dependence of the perpendicular heating of H + -He ++ solar wind plasma is studied. It is found that He ++ ions are heated efficiently by the Alfvénic wave spectrum below the proton gyroperiod.

Description: On the effect of IMF turning on ion dynamics at Mercury Annales Geophysicae, 29, 987-996, 2011 Author(s): D. C. Delcourt, T. E. Moore, and M.-C. H. Fok We investigate the effect of a rotation of the Interplanetary Magnetic Field (IMF) on the transport of magnetospheric ion populations at Mercury. We focus on ions of planetary origin and investigate their large-scale circulation using three-dimensional single-particle simulations. We show that a nonzero B X component of the IMF leads to a pronounced asymmetry in the overall circulation pattern. In particular, we demonstrate that the centrifugal acceleration due to curvature of the E × B drift paths is more pronounced in one hemisphere than the other, leading to filling of the magnetospheric lobes and plasma sheet with more or less energetic material depending upon the hemisphere of origin. Using a time-varying electric and magnetic field model, we investigate the response of ions to rapid (a few tens of seconds) re-orientation of the IMF. We show that, for ions with gyroperiods comparable to the field variation time scale, the inductive electric field should lead to significant nonadiabatic energization, up to several hundreds of eVs or a few keVs. It thus appears that IMF turning at Mercury should lead to localized loading of the magnetosphere with energetic material of planetary origin (e.g., Na + ).

Description: Instability of combined gravity-inertial-Rossby waves in atmospheres and oceans Annales Geophysicae, 29, 997-1003, 2011 Author(s): J. F. McKenzie The properties of the instability of combined gravity-inertial-Rossby waves on a β-plane are investigated. The wave-energy exchange equation shows that there is an exchange of energy with the background stratified medium. The energy source driving the instability lies in the background enthalpy released by the gravitational buoyancy force. It is shown that if the phase speed of the westward propagating low frequency-long wavelength Rossby wave exceeds the Poincaré-Kelvin (or "equivalent" shallow water) wave speed, instability arises from the merging of Rossby and Poincaré modes. There are two key parameters in this instability condition; namely, the equatorial/rotational Mach (or Froude) number M and the latitude θ 0 of the β-plane. In general waves equatorward of a critical latitude for given M can be driven unstable, with corresponding growth rates of the order of a day or so. Although these conclusions may only be safely drawn for short wavelengths corresponding to a JWKB wave packet propagating internally and located far from boundaries, nevertheless such a local instability may play a significant role in atmosphere-ocean dynamics.

Description: The correlation between solar and geomagnetic activity – Part 3: An integral response model Annales Geophysicae, 29, 1005-1018, 2011 Author(s): Z. L. Du An integral response model is proposed to describe the relationship between geomagnetic activity ( aa index) and solar activity (represented by sunspot number R z ): The aa at a given time t is the integral of R z at past times ( t '≤ t ) multiplied by an exponential decay factor of the time differences ( e −(t−t')/τ ), where τ is the decay time scale (~40 months). The correlation coefficient of aa with the reconstructed series based on this model ( r f =0.85) is much higher than that of aa with R z ( r 0 =0.61). If this model is applied to each solar cycle, the correlation coefficient will be higher ( r f =0.95). This model can naturally explain some phenomena related to aa and R z , such as (i) the significant increase in the aa index (and its baseline) over the twentieth century; (ii) the longer lag times of aa to R z at solar cycle maxima than at minima; and (iii) the variations in the correlations related to solar and Hale cycles. These results demonstrate that aa depends not only on the present R z but also on past values. The profile of aa can be better predicted from R z by this model than by point-point correspondence.

Description: Cluster observations of trapped ions interacting with magnetosheath mirror modes Annales Geophysicae, 29, 1049-1060, 2011 Author(s): J. Soucek and C. P. Escoubet Mirror modes are among the most intense low frequency plasma wave phenomena observed in the magnetosheaths of magnetized planets. They appear as large amplitude non-propagating fluctuations in the magnetic field magnitude and plasma density. These structures are widely accepted to represent a non-linear stage of the mirror instability, dominant in plasmas with large ion beta and a significant ion temperature anisotropy T ⊥ / T ∥ 〉1. It has long been recognized that the mirror instability both in the linear and non-linear stage is a kinetic process and that the behavior of resonant particles at small parallel velocities is crucial for its development and saturation. While the dynamics of the instability and the effect of trapped particles have been studied extensively in theoretical models and numerical simulations, only spurious observations of the trapped ions were published to date. In this work we used data from the Cluster spacecraft to perform the first detailed experimental study of ion velocity distribution associated with mirror mode oscillations. We show a conclusive evidence for the predicted cooling of resonant ions at small parallel velocities and heating of trapped ions at intermediate pitch angles.

Description: MHD waves generated by high-frequency photospheric vortex motions Annales Geophysicae, 29, 1029-1035, 2011 Author(s): V. Fedun, S. Shelyag, G. Verth, M. Mathioudakis, and R. Erdélyi In this paper, we discuss simulations of MHD wave generation and propagation through a three-dimensional open magnetic flux tube in the lower solar atmosphere. By using self-similar analytical solutions for modelling the magnetic field in Cartesian coordinate system, we have constructed a 3-D magnetohydrostatic configuration which is used as the initial condition for non-linear MHD wave simulations. For a driver we have implemented a high-frequency vortex-type motion at the footpoint region of the open magnetic flux tube. It is found that the implemented swirly source is able to excite different types of wave modes, i.e. sausage, kink and torsional Alfvén modes. Analysing these waves by magneto-seismology tools could provide insight into the magnetic structure of the lower solar atmosphere.

Description: Average properties and small-scale variations of the mesospheric Na and Fe layers as observed simultaneously by two closely colocated lidars at 30° N Annales Geophysicae, 29, 1037-1048, 2011 Author(s): L. Chen and F. Yi We report the average properties and small-scale variation features of the mesospheric Na and Fe layers at 30° N from extensive simultaneous and common-volume Na and Fe lidar measurements at Wuhan, China. The annual mean Na and Fe density profiles are derived in terms of an averaging method taken from an early literature. The mean Na and Fe profiles preserve the sharp gradients present in most individual density profiles near the layer bottom. Near the bottommost of the layers the mean Na and Fe scale heights are respectively −0.42 and −0.30 km. The mean layer parameters coincide well with the previous report. The Na and Fe densities in the lowest several kilometers of the layers consistently exhibit nearly the same time variations. A clear-cut distinction between the Na and Fe time variations always appears in an altitude range near 90 km. A relatively weak positive correlation between them persistently occurs also in an altitude range near 100 km. The mean increase and decrease rates for both Na and Fe are altitude dependent and have a single-peak structure. The time constant of the layer variation is ~0.07–2.0 h for Na and ~0.02–1.7 h for Fe, suggesting that the variability is dominated by small-scale processes. However, there is also a slow net increase in each of the annual mean column abundances (Na and Fe) during night.

Description: Polarisation in the auroral red line during coordinated EISCAT Svalbard Radar/optical experiments Annales Geophysicae, 29, 1101-1112, 2011 Author(s): M. Barthélémy, J. Lilensten, F. Pitout, C. Simon Wedlund, R. Thissen, D. Lorentzen, F. Sigernes, J. Moen, G. Gronoff, I. McCrea, H. Rothkael, H. Ménager, and A. Aruliah The polarisation of the atomic oxygen red line in the Earth's thermosphere is observed in different configurations with respect to the magnetic field line at high latitude during several coordinated Incoherent Scatter radar/optical experiment campaigns. When pointing northward with a line-of-sight nearly perpendicular to the magnetic field, we show that, as expected, the polarisation is due to precipitated electrons with characteristic energies of a few hundreds of electron Volts. When pointing toward the zenith or southward with a line-of-sight more parallel to the magnetic field, we show that the polarisation practically disappears. This confirms experimentally the predictions deduced from the recent discovery of the red line polarisation. We show that the polarisation direction is parallel to the magnetic field line during geomagnetic activity intensification and that these results are in agreement with theoretical work.

Description: New insights into the stratospheric and mesosphere-lower thermospheric ozone response to the abrupt changes in solar forcing Annales Geophysicae, 29, 1093-1099, 2011 Author(s): Karanam Kishore Kumar, K. V. Subrahmanyam, and Sherine Rachel John Using a unique set of satellite based observations of the vertical distribution of ozone during the recent annular solar eclipse of 15 January 2010, we demonstrate for the first time, a complete picture of the response of stratospheric ozone to abrupt changes in solar forcing. The stratospheric ozone decreased after the maximum obscuration of the Sun and then gradually increased with time. A dramatic increase in stratospheric ozone of up to 4 ppmv is observed 3 h after the maximum obscuration of the Sun. The present study also reports for the first time the mesosphere-lower thermospheric ozone response to solar eclipse. Thus it is envisaged that the present results will have important implications in understanding the ozone response to abrupt changes in solar forcing and time-scales involved in such response.

Description: Corrigendum to "The PC index: review of methods" published in Ann. Geophys., 28, 1887–1903, 2010 Annales Geophysicae, 29, 813-814, 2011 Author(s): H. McCreadie and M. Menvielle No abstract available.

Description: Pre-monsoon aerosol characteristics over the Indo-Gangetic Basin: implications to climatic impact Annales Geophysicae, 29, 789-804, 2011 Author(s): A. K. Srivastava, S. Tiwari, P. C. S. Devara, D. S. Bisht, Manoj K. Srivastava, S. N. Tripathi, P. Goloub, and B. N. Holben Sun/sky radiometer observations over the Indo-Gangetic Basin (IGB) region during pre-monsoon (from April–June 2009) have been processed to analyze various aerosol characteristics in the central and eastern IGB region, represented by Kanpur and Gandhi College, respectively, and their impacts on climate in terms of radiative forcing. Monthly mean aerosol optical depth (AOD at 500 nm) and corresponding Angstrom Exponent (AE at 440–870 nm, given within the brackets) was observed to be about 0.50 (0.49) and 0.51 (0.65) in April, 0.65 (0.74) and 0.67 (0.91) in May and 0.69 (0.45) and 0.77 (0.71) in June at Kanpur and Gandhi College, respectively. Results show a positive gradient in AOD and AE from central to eastern IGB region with the advancement of the pre-monsoon, which may be caused due to diverse geographical location of the stations having different meteorological conditions and emission sources. Relatively lower SSA was observed at the eastern IGB (0.89) than the central IGB (0.92) region during the period, which suggests relative dominance of absorbing aerosols at the eastern IGB as compared to central IGB region. The absorbing aerosol optical properties over the station suggest that the atmospheric absorption over central IGB region is mainly due to dominance of coarse-mode dust particles; however, absorption over eastern IGB region is mainly due to dominance of fine-particle pollution. The derived properties from sun/sky radiometer during pre-monsoon period are used in a radiative-transfer model to estimate aerosol radiative forcing at the top-of-the atmosphere (TOA) and at the surface over the IGB region. Relatively large TOA and surface cooling was observed at the eastern IGB as compared to the central IGB region. This translates into large heating of the atmosphere ranging from 0.45 to 0.55 K day −1 at Kanpur and from 0.45 to 0.59 K day −1 at Gandhi College.

Description: Annual and semiannual variations of vertical total electron content during high solar activity based on GPS observations Annales Geophysicae, 29, 865-873, 2011 Author(s): M. P. Natali and A. Meza Annual, semiannual and seasonal variations of the Vertical Total Electron Content (VTEC) have been investigated during high solar activity in 2000. In this work we use Global IGS VTEC maps and Principal Component Analysis to study spatial and temporal ionospheric variability. The behavior of VTEC variations at two-hour periods, at noon and at night is analyzed. Particular characteristics associated with each period and the geomagnetic regions are highlighted. The variations at night are smaller than those obtained at noon. At noon it is possible to see patterns of the seasonal variation at high latitude, and patterns of the semiannual anomaly at low latitudes with a slow decrease towards mid latitudes. At night there is no evidence of seasonal or annual anomaly for any region, but it was possible to see the semiannual anomaly at low latitudes with a sudden decrease towards mid latitudes. In general, the semiannual behavior shows March–April equinox at least 40 % higher than September one. Similarities and differences are analyzed also with regard to the same analysis done for a period of low solar activity.

Description: Modelling of spacecraft spin period during eclipse Annales Geophysicae, 29, 875-882, 2011 Author(s): E. Georgescu, F. Plaschke, U. Auster, K.-H. Fornaçon, and H. U. Frey The majority of scientific satellites investigating the Earth magnetosphere are spin stabilized. The attitude information comes usually from a sun sensor and is missing in the umbra; hence, the accurate experimental determination of vector quantities is not possible during eclipses. The spin period of the spacecraft is generally not constant during these times because the moment of inertia changes due to heat dissipation. The temperature dependence of the moment of inertia for each spacecraft has a specific signature determined by its design and distribution of mass. We developed an "eclipse-spin" model for the spacecraft spin period behaviour using magnetic field vector measurements close to the Earth, where the magnetic field is dominated by the dipole field, and in the magnetospheric lobes, where the magnetic field direction is mostly constant. The modelled spin periods give us extraordinarily good results with accumulated phase deviations over one hour of less than 10 degrees. Using the eclipse spin model satellite experiments depending on correct spin phase information can deliver science data even during eclipses. Two applications for THEMIS B, one in the lobe and the other in the lunar wake, are presented.

Description: Apparent temperature anisotropies due to wave activity in the solar wind Annales Geophysicae, 29, 909-917, 2011 Author(s): D. Verscharen and E. Marsch The fast solar wind is a collisionless plasma permeated by plasma waves on many different scales. A plasma wave represents the natural interplay between the periodic changes of the electromagnetic field and the associated coherent motions of the plasma particles. In this paper, a model velocity distribution function is derived for a plasma in a single, coherent, large-amplitude wave. This model allows one to study the kinetic effects of wave motions on particle distributions. They are by in-situ spacecraft measured by counting, over a certain sampling time, the particles coming from various directions and having different energies. We compare our results with the measurements by the Helios spacecraft, and thus find that by assuming high wave activity we are able to explain key observed features of the measured distributions within the framework of our model. We also address the recent discussions on nonresonant wave–particle interactions and apparent heating. The applied time-averaging procedure leads to an apparent ion temperature anisotropy which is connected but not identical to the intrinsic temperature of the underlying distribution function.

Description: The numerical simulation on ionospheric perturbations in electric field before large earthquakes Annales Geophysicae, 32, 1487-1493, 2014 Author(s): S. F. Zhao, X. M. Zhang, Z. Y. Zhao, and X. H. Shen Many observational results have shown electromagnetic abnormality in the ionosphere before large earthquakes. The theoretical simulation can help us to understand the internal mechanism of these anomalous electromagnetic signals resulted from seismic regions. In this paper, the horizontal and vertical components of electric and magnetic field at the topside ionosphere are simulated by using the full wave method that is based on an improved transfer matrix method in the lossy anisotropic horizontally stratified ionosphere. Taken account into two earthquakes with electric field perturbations recorded by the DEMETER satellite, the numerical results reveal that the propagation and penetration of ULF (ultra-low-frequency) electromagnetic waves into the ionosphere is related to the spatial distribution of electron and ion densities at different time and locations, in which the ion density has less effect than electron density on the field intensity. Compared with different frequency signals, the minimum values of electric and magnetic field excited by earthquakes can be detected by satellite in current detection capability have also been calculated, and the lower frequency wave can be detected easier.

Description: A comparison between VEGA 1, 2 and Giotto flybys of comet 1P/Halley: implications for Rosetta Annales Geophysicae, 32, 1441-1453, 2014 Author(s): M. Volwerk, K.-H. Glassmeier, M. Delva, D. Schmid, C. Koenders, I. Richter, and K. Szegö Three flybys of comet 1P/Halley, by VEGA 1, 2 and Giotto, are investigated with respect to the occurrence of mirror mode waves in the cometosheath and field line draping in the magnetic pile-up region around the nucleus. The time interval covered by these flybys is approximately 8 days, which is also the approximate length of an orbit or flyby of Rosetta around comet 67P/Churyumov–Gerasimenko. Thus any significant changes observed around Halley are changes that might occur for Rosetta during one pass of 67P/CG. It is found that the occurrence of mirror mode waves in the cometosheath is strongly influenced by the dynamical pressure of the solar wind and the outgassing rate of the comet. Field line draping happens in the magnetic pile-up region. Changes in nested draping regions (i.e. regions with different B x directions) can occur within a few days, possibly influenced by changes in the outgassing rate of the comet and thereby the conductivity of the cometary ionosphere.

Description: Magnetic curvature driven Rayleigh-Taylor instability revisited Annales Geophysicae, 29, 411-413, 2011 Author(s): O. A. Pokhotelov and O. G. Onishchenko The problem of incomplete finite ion Larmor radius (FLR) stabilization of the magnetic curvature driven Rayleigh-Taylor instability (RTI) in low beta plasma with homogeneous ion temperature is investigated. For this purpose a model hydrodynamic description of nonlinear flute waves with arbitrary spatial scales compared to the ion Larmor radius is developed. It is shown that the RTI is not stabilized by FLR effects in a plasma with cold electrons when the ratio of characteristic spatial scale of the plasma inhomogeneity to local effective radius of curvature of the magnetic field lines is larger than 1/4. The crucial role in the absence of the complete FLR stabilization plays the contribution of the compressibility of the polarization part of the ion velocity.

Description: Forecasting relativistic electron flux using dynamic multiple regression models Annales Geophysicae, 29, 415-420, 2011 Author(s): H.-L. Wei, S. A. Billings, A. Surjalal Sharma, S. Wing, R. J. Boynton, and S. N. Walker The forecast of high energy electron fluxes in the radiation belts is important because the exposure of modern spacecraft to high energy particles can result in significant damage to onboard systems. A comprehensive physical model of processes related to electron energisation that can be used for such a forecast has not yet been developed. In the present paper a systems identification approach is exploited to deduce a dynamic multiple regression model that can be used to predict the daily maximum of high energy electron fluxes at geosynchronous orbit from data. It is shown that the model developed provides reliable predictions.

Description: Airglow observations over the equatorial ionization anomaly zone in Taiwan Annales Geophysicae, 29, 749-757, 2011 Author(s): J. Y. Liu, P. K. Rajesh, I. T. Lee, and T. C. Chow Airglow imaging at mid-latitude stations often show intensity modulations associated with medium scale travelling ionospheric disturbances (MSTID), while those carried out near the equatorial regions reveal depletions caused by equatorial plasma bubbles (EPB). Two all sky cameras are used to observe plasma depletions in the 630.0 nm emission over the equatorial ionization anomaly (EIA) region, Taiwan (23° N, 121° E; 13.5° N Magnetic) during 1998–2002 and 2006–2007. The results show EPB and MSTID depletions in different solar activity conditions. Several new features of the EPB depletions such as bifurcation, secondary structure on the walls, westward tilt, etc., are discussed in this paper. Evidence of tilted depletions with secondary structures developing on the eastern wall that later evolve to appear as bifurcations, are presented for the first time. Moreover, detail investigations are carried out using International Reference Ionosphere (IRI) model as well as the electron density from Ionosonde and Global Positioning System (GPS) Occultation Experiment (GOX) onboard FORMOSAT-3/COSMIC satellite, to understand the conditions that favor the propagation of MSTID to the latitude of Taiwan.

Description: 3-D mesoscale MHD simulations of a cusp-like magnetic configuration: method and first results Annales Geophysicae, 29, 759-770, 2011 Author(s): E. Adamson and K. Nykyri We present a local mesoscale model of the magnetospheric cusp region with high resolution (up to 300 km). We discuss the construction and implementation of the initial configuration and give a detailed description of the numerical simulation. An overview of simulation results for the case of strongly northward interplanetary magnetic field (IMF) is then presented and compared with data from Cluster 2 spacecraft from 14 February 2003. Results show a cusp diamagnetic cavity (CDC) with depth normal to the magnetospheric boundary on the order of 1–2 R E and a much larger extent of ~5–9 R E tangential to the boundary, bounded by a gradual inner boundary with the magnetospheric lobe and a more distinct exterior boundary with the magnetosheath. These results are qualitatively consistent with observational data.

Description: Polar cap ion beams during periods of northward IMF: Cluster statistical results Annales Geophysicae, 29, 771-787, 2011 Author(s): R. Maggiolo, M. Echim, J. De Keyser, D. Fontaine, C. Jacquey, and I. Dandouras Above the polar caps and during prolonged periods of northward IMF, the Cluster satellites detect upward accelerated ion beams with energies up to a few keV. They are associated with converging electric field structures indicating that the acceleration is caused by a quasi-static field-aligned electric field that can extend to altitudes higher than 7 R E (Maggiolo et al., 2006; Teste et al., 2007). Using the AMDA science analysis service provided by the Centre de Données de la Physique des Plasmas, we have been able to extract about 200 events of accelerated upgoing ion beams above the polar caps from the Cluster database. Most of these observations are taken at altitudes lower than 7 R E and in the Northern Hemisphere. We investigate the statistical properties of these ion beams. We analyze their geometry, the properties of the plasma populations and of the electric field inside and around the beams, as well as their dependence on solar wind and IMF conditions. We show that ~40 % of the ion beams are collocated with a relatively hot and isotropic plasma population. The density and temperature of the isotropic population are highly variable but suggest that this plasma originates from the plasma sheet. The ion beam properties do not change significantly when the isotropic, hot background population is present. Furthermore, during one single polar cap crossing by Cluster it is possible to detect upgoing ion beams both with and without an accompanying isotropic component. The analysis of the variation of the IMF B Z component prior to the detection of the beams indicates that the delay between a northward/southward turning of IMF and the appearance/disappearance of the beams is respectively ~2 h and 20 min. The observed electrodynamic characteristics of high altitude polar cap ion beams suggest that they are closely connected to polar cap auroral arcs. We discuss the implications of these Cluster observations above the polar cap on the magnetospheric dynamics and configuration during prolonged periods of northward IMF.

Description: Diurnal variation of short-period (20–120 min) gravity waves in the equatorial Mesosphere and Lower Thermosphere and its relation to deep tropical convection Annales Geophysicae, 29, 623-629, 2011 Author(s): N. Venkateswara Rao, Y. Shibagaki, and T. Tsuda We study short period gravity waves (20–120 min) in the equatorial Mesosphere and Lower Thermosphere (MLT) using a Medium Frequency (MF) radar at Pameungpeuk (7.4° S, 107.4° E), Indonesia. In particular, we study local time and seasonal variation of the gravity wave variance and its relation to tropical convection. The gravity wave variance at 88 km enhances between 20:00 LT and 07:00 LT, with a peak at 02:00–03:00 LT. The enhancement is mainly observed during February–April and September–October and shows inter-annual variability. Convective activity over the same location persists from 16:00–21:00 LT with a peak activity ~18:00 LT and enhances between November–April. Time delay between the peak of convection and that of gravity wave activity ranges 1–15 h, which is consistent with theoretical calculations and previous reports based on reverse ray tracing analysis.

Description: Observations of AGW/TID propagation across the polar cap: a case study Annales Geophysicae, 29, 1355-1363, 2011 Author(s): H. T. Cai, F. Yin, S. Y. Ma, and I. W. McCrea In this paper, we present observational evidence for the trans-polar propagation of large-scale Traveling Ionospheric Disturbances (TIDs) from their nightside source region to the dayside. On 13 February 2001, the 32 m dish of EISCAT Svalbard Radar (ESR) was directing toward the geomagnetic pole at low elevation (30°) during the interval 06:00–12:00 UT (MLT ≈ UT + 3 h), providing an excellent opportunity to monitor the ionosphere F-region over the polar cap. The TIDs were first detected by the ESR over the dayside north polar cap, propagating equatorward, and were subsequently seen by the mainland UHF radar at auroral latitudes around geomagnetic local noon. The propagation properties of the observed ionization waves suggest the presence of a moderately large-scale TIDs, propagating across the northern polar cap from the night-time auroral source during substorm conditions. Our results agree with the theoretical simulations by Balthazor and Moffett (1999) in which poleward-propagating large-scale traveling atmospheric disturbances were found to be self-consistently driven by enhancements in auroral heating.

Description: Statistical evidence for O + energization and outflow caused by wave-particle interaction in the high altitude cusp and mantle Annales Geophysicae, 29, 945-954, 2011 Author(s): M. Waara, R. Slapak, H. Nilsson, G. Stenberg, M. André, and I. A. Barghouthi We present a statistical study of the low (

Description: O + heating associated with strong wave activity in the high altitude cusp and mantle Annales Geophysicae, 29, 931-944, 2011 Author(s): R. Slapak, H. Nilsson, M. Waara, M. André, G. Stenberg, and I. A. Barghouthi We use the Cluster spacecraft to study three events with intense waves and energetic oxygen ions (O + ) in the high altitude cusp and mantle. The ion energies considered are of the order 1000 eV and higher, observed above an altitude of 8 earth radii together with high wave power at the O + gyrofrequency. We show that heating by waves can explain the observed high perpendicular energy of O + ions, using a simple gyroresonance model and 25–45% of the observed wave spectral density at the gyrofrequency. This is in contrast to a recently published study where the wave intensity was too low to explain the observed high altitude ion energies. Long lasting cases (〉10 min) of high perpendicular-to-parallel temperature ratios are sometimes associated with low wave activity, suggesting that high perpendicular-to-parallel temperature ratio is not a good indicator of local heating. Using multiple spacecraft, we show that the regions of enhanced wave activity are at least one order of magnitude larger than the gyroradius of the heated ions.

Description: Simulation of aerosol optical properties over a tropical urban site in India using a global model and its comparison with ground measurements Annales Geophysicae, 29, 955-963, 2011 Author(s): D. Goto, K. V. S. Badarinath, T. Takemura, and T. Nakajima Aerosols have great impacts on atmospheric environment, human health, and earth's climate. Therefore, information on their spatial and temporal distribution is of paramount importance. Despite numerous studies have examined the variation and trends of BC and AOD over India, only very few have focused on their spatial distribution or even correlating the observations with model simulations. In the present study, a three-dimensional aerosol transport-radiation model coupled with a general circulation model. SPRINTARS, simulated atmospheric aerosol distributions including BC and aerosol optical properties, i.e., aerosol optical thickness (AOT), Ångström Exponent (AE), and single scattering albedo (SSA). The simulated results are compared with both BC measurements by aethalometer and aerosol optical properties measured by ground-based skyradiometer and by satellite sensor, MODIS/Terra over Hyderabad, which is a tropical urban area of India, for the year 2008. The simulated AOT and AE in Hyderabad are found to be comparable to ground-based measured ones. The simulated SSA tends to be higher than the ground-based measurements. Both these comparisons of aerosol optical properties between the simulations with different emission inventories and the measurements indicate that, firstly the model uncertainties derived from aerosol emission inventory cannot explain the gaps between the simulations and the measurements and secondly the vertical transport of BC and the treatment of BC-containing particles can be the main issue in the global model to solve the gap.

Description: Energy and flux variations across thin auroral arcs Annales Geophysicae, 29, 1699-1712, 2011 Author(s): H. Dahlgren, B. Gustavsson, B. S. Lanchester, N. Ivchenko, U. Brändström, D. K. Whiter, T. Sergienko, I. Sandahl, and G. Marklund Two discrete auroral arc filaments, with widths of less than 1 km, have been analysed using multi-station, multi-monochromatic optical observations from small and medium field-of-view imagers and the EISCAT radar. The energy and flux of the precipitating electrons, volume emission rates and local electric fields in the ionosphere have been determined at high temporal (up to 30 Hz) and spatial (down to tens of metres) resolution. A new time-dependent inversion model is used to derive energy spectra from EISCAT electron density profiles. The energy and flux are also derived independently from optical emissions combined with ion-chemistry modelling, and a good agreement is found. A robust method to obtain detailed 2-D maps of the average energy and number flux of small scale aurora is presented. The arcs are stretched in the north-south direction, and the lowest energies are found on the western, leading edges of the arcs. The large ionospheric electric fields (250 mV m −1 ) found from tristatic radar measurements are evidence of strong currents associated with the region close to the optical arcs. The different data sets indicate that the arcs appear on the boundaries between regions with different average energy of diffuse precipitation, caused by pitch-angle scattering. The two thin arcs on these boundaries are found to be related to an increase in number flux (and thus increased energy flux) without an increase in energy.

Description: Storm-time total electron content and its response to penetration electric fields over South America Annales Geophysicae, 29, 1765-1778, 2011 Author(s): P. M. de Siqueira, E. R. de Paula, M. T. A. H. Muella, L. F. C. Rezende, M. A. Abdu, and W. D. Gonzalez In this work the response of the ionosphere due to the severe magnetic storm of 7–10 November 2004 is investigated by analyzing GPS Total Electron Content (TEC) maps constructed for the South America sector. In order to verify the disturbed zonal electric fields in South America during the superstorm, ionospheric vertical drift data obtained from modeling results are used in the analysis. The vertical drifts were inferred from Δ H magnetometer data (Jicamarca-Piura) following the methodology presented by Anderson et al. (2004). Also used were vertical drifts measured by the Jicamarca ISR. Data from a digisonde located at São Luís, Brazil (2.33° S, 44.2° W, dip latitude 0.25°) are presented to complement the Jicamarca equatorial data. Penetration electric fields were observed by the comparison between the equatorial vertical drifts and the Interplanetary Electric Field (IEF). The TEC maps obtained from GPS data reflect the ionospheric response over the South America low-latitude and equatorial region. They reveal unexpected plasma distributions and TEC levels during the main phase of the superstorm on 7 November, which is coincident with the local post-sunset hours. At this time an increase in the pre-reversal enhancement was expected to develop the Equatorial Ionization Anomaly (EIA) but we observed the absence of EIA. The results also reveal well known characteristics of the plasma distributions on 8, 9, and 10 November. The emphasized features are the expansion and intensification of EIA due to prompt penetration electric fields on 9 November and the inhibition of EIA during post-sunset hours on 7, 8, and 10 November. One important result is that the TEC maps provided a bi-dimensional view of the ionospheric changes offering a spatial description of the electrodynamics involved, which is an advantage over TEC measured by isolated GPS receivers.

Description: Climatology of rapid geomagnetic variations at high latitudes over two solar cycles Annales Geophysicae, 29, 1783-1792, 2011 Author(s): A. Viljanen and E. Tanskanen We investigate the characteristics of rapid geomagnetic variations at high latitudes based on the occurrence of large time derivatives of the horizontal magnetic field ( d H / dt exceeding 1 nT s −1 ). Analysis of IMAGE magnetometer data from North Europe in 1983–2010, covering more than two solar cycles, confirms and specifies several previous findings. We show that d H / dt activity is high around the midnight and early morning hours, and nearly vanishes at noon and early afternoon. This happens during all seasons, although the midnight maximum is nearly invisible during summer. As indicated by modelled ionospheric equivalent currents, large d H / dt values occur predominantly during westward ionospheric electrojets. Before and around midnight, d H / dt tends to be north-south oriented, whereas in the morning hours, its direction is more west-east directed. d H / dt tends to be more strictly north-south oriented during winter than other seasons. The seasonal occurrence of large d H / dt values is similar to the variation of the maximum amplitude of westward equivalent currents. The yearly fraction of east-west directed large d H / dt vectors at the Kilpisjärvi station (MLAT 65.88) varies from 31 to 47 % without any clear correlation with the general geomagnetic activity nor with the yearly averages of solar wind parameters.

Description: Is thermospheric long-term cooling due to CO 2 or O 3 ? Annales Geophysicae, 29, 1779-1782, 2011 Author(s): P. L. Walsh and W. L. Oliver While greenhouse gases trap heat emanating from the Earth and thereby heat the surface atmosphere, they act as emitters in the high atmosphere and cool the air there. In 1989 Roble and Dickinson (1989) estimated the cooling that would occur in the thermosphere (250–500 km altitude) due to a doubling of greenhouse gas densities. Ever since, long-term data bases have been scoured for evidence of this thermospheric "global cooling." Here we show evidence that the thermosphere did indeed cool over the period 1966–1987, but the data suggest that the cooling accelerated at a "breakpoint year" around 1979 to a rate far larger than may be attributed to greenhouse cooling. This 1979 breakpoint year appears to coincide with a breakpoint year in ozone (O 3 ) column density. Further, the cooling was confined largely to the daytime thermosphere while the nighttime showed only a small trend. These results suggest, first, that the greenhouse cooling of the thermosphere may well not be detectable with current data sets and, second, that the long-term cooling that is clearly seen may be due largely to O 3 depletion.

Description: Characteristics of the terrestrial field-aligned current system Annales Geophysicae, 29, 1713-1729, 2011 Author(s): J. W. Gjerloev, S. Ohtani, T. Iijima, B. Anderson, J. Slavin, and G. Le We present the first ever comprehensive statistical study of the spatiotemporal characteristics of field-aligned currents in the terrestrial magnetosphere-ionosphere system using multi point measurements. We determine how the FAC density, variability and scale size are coupled. The three ST 5 satellites were in a pearls-on-a-string formation making measurements of the magnetic field with variable inter-spacecraft separations ranging from a few seconds to about 10 min. More than 4700 sets of satellite passes are analyzed using a robust correlation analysis aimed at determining the variability of the FAC system as a function of scale size and satellite spacing. We find significant differences between the FAC characteristics on the dayside and on the nightside in terms of dynamics of the current systems. On the dayside the FAC characteristics are found to be independent of IMF B z and geomagnetic activity while the nightside indicates increased variability during disturbed conditions. The boundary separating highly and poorly correlated FACs can be fitted by a linear line for satellite separations shorter than 60 s (dayside) and 160 s (nightside). We interpret this as the dayside and nightside magnetospheric reconfiguration times respectively. For times exceeding this the FAC characteristics are suggested to be controlled by the solar wind (dayside) and plasma sheet (nightside) dynamics. Finally, the characteristics of FAC system with scale sizes larger than ~200 km (at ionospheric altitude) appear to be stable and repeatable on time scales of the order of a minute (i.e. comparable to the low-altitude orbiting satellite's traverse time across the auroral belt). In this sense, our results effectively validate the Iijima and Potemra (1978) assumption that on average the large-scale currents with scale sizes of the Region1 and Region2 are quasi-persistently significant in the transport of energy and momentum between the magnetosphere and the ionosphere.

Description: Whistler-Langmuir oscillitons and their relation to auroral hiss Annales Geophysicae, 29, 1739-1753, 2011 Author(s): K. Sauer and R. D. Sydora A new type of oscilliton (soliton with superimposed spatial oscillations) is described which arises in plasmas if the electron cyclotron frequency Ω e is larger than the electron plasma frequency ω e , which is a typical situation for auroral regions in planetary magnetospheres. Both high-frequency modes of concern, the Langmuir and the whistler wave, are completely decoupled if they propagate parallel to the magnetic field. However, for oblique propagation two mixed modes are created with longitudinal and transverse electric field components. The lower mode (in the literature commonly called the whistler mode, e.g. Gurnett et al., 1983) has whistler wave characteristics at small wave numbers and asymptotically transforms into the Langmuir mode. As a consequence of the coupling between these two modes, with different phase velocity dependence, a maximum in phase velocity appears at finite wave number. The occurrence of such a particular point where phase and group velocity coincide creates the condition for the existence of a new type of oscillating nonlinear stationary structure, which we call the whistler-Langmuir (WL) oscilliton. After determining, by means of stationary dispersion theory, the parameter regime in which WL oscillitons exist, their spatial profiles are calculated within the framework of cold (non-relativistic) fluid theory. Particle-in-cell (PIC) simulations are used to demonstrate the formation of WL oscillitons which seem to play an important role in understanding electron beam-excited plasma radiation that is observed as auroral hiss in planetary magnetospheres far away from the source region.

Description: Density and temperature of energetic electrons in the Earth's magnetotail derived from high-latitude GPS observations during the declining phase of the solar cycle Annales Geophysicae, 29, 1755-1763, 2011 Author(s): M. H. Denton and T. E. Cayton Single relativistic-Maxwellian fits are made to high-latitude GPS-satellite observations of energetic electrons for the period January 2006–November 2010; a constellation of 12 GPS space vehicles provides the observations. The derived fit parameters (for energies ~0.1–1.0 MeV), in combination with field-line mapping on the nightside of the magnetosphere, provide a survey of the energetic electron density and temperature distribution in the magnetotail between McIlwain L-values of L =6 and L =22. Analysis reveals the characteristics of the density-temperature distribution of energetic electrons and its variation as a function of solar wind speed and the Kp index. The density-temperature characteristics of the magnetotail energetic electrons are very similar to those found in the outer electron radiation belt as measured at geosynchronous orbit. The energetic electron density in the magnetotail is much greater during increased geomagnetic activity and during fast solar wind. The total electron density in the magnetotail is found to be strongly correlated with solar wind speed and is at least a factor of two greater for high-speed solar wind ( V SW =500–1000 km s −1 ) compared to low-speed solar wind ( V SW =100–400 km s −1 ). These results have important implications for understanding (a) how the solar wind may modulate entry into the magnetosphere during fast and slow solar wind, and (b) if the magnetotail is a source or a sink for the outer electron radiation belt.

Description: Lower hybrid resonances stimulated by the four CLUSTER relaxation sounders deep inside the plasmasphere: observations and inferred plasma characteristics Annales Geophysicae, 29, 2003-2018, 2011 Author(s): S. Kougblénou, G. Lointier, P. M. E. Décréau, J.-G. Trotignon, J.-L. Rauch, X. Vallières, P. Canu, A. Masson, and J. Pickett The frequency range of the WHISPER relaxation sounder instrument on board CLUSTER, 4–80 kHz, has been chosen so as to encompass the electron gyro-frequency, F ce , and the electron plasma frequency, F p , in most regions to be explored. Measurement of those frequencies, which are triggered as resonances by the sounder, provides a direct estimation of in situ fundamental plasma characteristics: electron density and magnetic field intensity. In the late mission phase, CLUSTER penetrated regions deep inside the plasmasphere where F ce and F p are much higher than the upper frequency of the sounder's range. However, they are of the right order of magnitude as to place the lower hybrid frequency, F lh , in the 4–15 kHz band. This characteristic frequency, placed at a resonance of the medium, is triggered by the sounder's transmitter and shows up as an isolated peak in the received spectrum, not present in spectra of naturally occuring VLF waves. This paper illustrates, from analysis of case events, how measured F lh values give access to a plasma diagnostic novel of its kind. CLUSTER, travelling along its orbit, encounters favourable conditions where F ce is increasing and F p decreasing, such that F ce / F p increases from values below unity to values above unity. Measured F lh values thus give access, in turn, to the effective mass, M eff , indicative of plasma ion composition, and to the core plasmasphere electron density value, a parameter difficult to measure. The analysed case events indicate that the estimated quantities ( M eff in the 1.0–1.4 range, N e in the 5 × 10 2 –10 4 cm −3 range) are varying with external factors (altitude, L value, geomagnetic activity) in a plausible way. Although covering only a restricted region (mid-latitude, low altitude inner plasmasphere), these measurements are available, since late 2009, for all CLUSTER perigee passes not affected by eclipses (on average, roughly a third of a total of ~200 passes per year) and offer multipoint observations previously unavailable in this region.

Description: Observations of electron vorticity in the inner plasma sheet Annales Geophysicae, 29, 1517-1527, 2011 Author(s): C. Gurgiolo, M. L. Goldstein, A. F. Viñas, W. H. Matthaeus, and A. N. Fazakerley From a limited number of observations it appears that vorticity is a common feature in the inner plasma sheet. With the four Cluster spacecraft and the four PEACE instruments positioned in a tetrahedral configuration, for the first time it is possible to directly estimate the electron fluid vorticity in a space plasma. We show examples of electron fluid vorticity from multiple plasma sheet crossings. These include three time periods when Cluster passed through a reconnection ion diffusion region. Enhancements in vorticity are seen in association with each crossing of the ion diffusion region.

Description: First detection of the O III 495.8911 and 500.6843 nm lines in the Earth's upper atmosphere Annales Geophysicae, 29, 2235-2238, 2011 Author(s): O. Witasse, T. Slanger, and R. Thissen We report the first detection of two emission lines of the atomic oxygen doubly-charged ion at 495.8911 and 500.6843 nm in the terrestrial upper atmosphere. They correspond to the transitions 1 D 2 - 3 P 1 and 1 D 2 - 3 P 2 of the O ++ ion, respectively. The measurements were performed on 30 October 2003 during the "Halloween" storms, with the Ultraviolet and Visual Echelle Spectrograph (UVES) mounted on the Very Large Telescope (VLT) in Chile. The intensities of these emissions are ~70 mRayleigh, and ~260 mRayleigh, respectively. These emissions constitute a new diagnostic of the state of the ionosphere.

Description: Solar wind and geomagnetism: toward a standard classification of geomagnetic activity from 1868 to 2009 Annales Geophysicae, 30, 421-426, 2012 Author(s): J. L. Zerbo, C. Amory Mazaudier, F. Ouattara, and J. D. Richardson We examined solar activity with a large series of geomagnetic data from 1868 to 2009. We have revisited the geomagnetic activity classification scheme of Legrand and Simon (1989) and improve their scheme by lowering the minimum Aa index value for shock and recurrent activity from 40 to 20 nT. This improved scheme allows us to clearly classify about 80% of the geomagnetic activity in this time period instead of only 60% for the previous Legrand and Simon classification.

Description: A statistical study of the performance of the Hakamada-Akasofu-Fry version 2 numerical model in predicting solar shock arrival times at Earth during different phases of solar cycle 23 Annales Geophysicae, 30, 405-419, 2012 Author(s): S. M. P. McKenna-Lawlor, C. D. Fry, M. Dryer, D. Heynderickx, K. Kecskemety, K. Kudela, and J. Balaz The performance of the Hakamada Akasofu-Fry, version 2 (HAFv.2) numerical model, which provides predictions of solar shock arrival times at Earth, was subjected to a statistical study to investigate those solar/interplanetary circumstances under which the model performed well/poorly during key phases (rise/maximum/decay) of solar cycle 23. In addition to analyzing elements of the overall data set (584 selected events) associated with particular cycle phases, subsets were formed such that those events making up a particular sub-set showed common characteristics. The statistical significance of the results obtained using the various sets/subsets was generally very low and these results were not significant as compared with the hit by chance rate (50%). This implies a low level of confidence in the predictions of the model with no compelling result encouraging its use. However, the data suggested that the success rates of HAFv.2 were higher when the background solar wind speed at the time of shock initiation was relatively fast. Thus, in scenarios where the background solar wind speed is elevated and the calculated success rate significantly exceeds the rate by chance, the forecasts could provide potential value to the customer. With the composite statistics available for solar cycle 23, the calculated success rate at high solar wind speed, although clearly above 50%, was indicative rather than conclusive. The RMS error estimated for shock arrival times for every cycle phase and for the composite sample was in each case significantly better than would be expected for a random data set. Also, the parameter "Probability of Detection, yes" (PODy) which presents the Proportion of Yes observations that were correctly forecast (i.e. the ratio between the shocks correctly predicted and all the shocks observed), yielded values for the rise/maximum/decay phases of the cycle and using the composite sample of 0.85, 0.64, 0.79 and 0.77, respectively. The statistical results obtained through detailed analysis of the available data provided insights into how changing circumstances on the Sun and in interplanetary space can affect the performance of the model. Since shock arrival predictions are widely utilized in making commercially significant decisions re. protecting space assets, the present detailed archival studies can be useful in future operational decision making during solar cycle 24. It would be of added value in this context to use Briggs-Rupert methodology to estimate the cost to an operator of acting on an incorrect forecast.

Description: Circulation changes in the winter lower atmosphere and long-lasting solar/geomagnetic activity Annales Geophysicae, 30, 1719-1726, 2012 Author(s): J. Bochníček, H. Davídkovová, P. Hejda, and R. Huth The paper describes the association between high long-lasting solar/geomagnetic activity and geopotential height (GPH) changes in the winter lower atmosphere, based on their development in the Northern Hemisphere in the winter periods (December–March) of 1950–1969 and 1970–2002. Solar/geomagnetic activity is characterised by the 60-day mean of the sunspot number R/by the 60-day mean of the daily sum of the Kp index. The GPH distributions in the lower atmosphere are described by 60-day anomalies from their long-term daily average at 20 hPa/850 hPa. The data have been adopted from the NCEP/NCAR reanalysis. The 60-day mean values of solar/geomagnetic activity and GPH anomalies were calculated in five-day steps over the whole winter period. The analysis was carried out using composite maps which represent their distribution of the GPH anomalies during high solar activity ( R ≥ 100) and high geomagnetic activity (ΣKp ≥ 20). Analysis has shown that the distribution of GPH anomalies depends on solar activity, geomagnetic activity and the phase of winter period (early or late winter). The nature of this relationship then depends on the time interval involved, i.e. 1950–1969 or 1970–2002. Positive anomalies in the polar stratosphere (20 hPa) were detected during the whole winter periods of the years 1950–1969. Significant anomalies were detected in the lower troposphere (850 hPa) during the second half of the winter period. The distribution of GPH anomalies on the maps compiled with regard to solar activity was similar to the distribution on maps compiled with regard to geomagnetic activity. In the interval 1970–2002, significant negative GPH anomalies were detected in the stratosphere at high latitudes, and positive anomalies were detected in the region of low latitudes. The distribution of GPH anomalies in the lower troposphere was substantially affected by situations in which, together with high solar activity, also high geomagnetic activity occurred.

Description: Support vector machines for TEC seismo-ionospheric anomalies detection Annales Geophysicae, 31, 173-186, 2013 Author(s): M. Akhoondzadeh Using time series prediction methods, it is possible to pursue the behaviors of earthquake precursors in the future and to announce early warnings when the differences between the predicted value and the observed value exceed the predefined threshold value. Support Vector Machines (SVMs) are widely used due to their many advantages for classification and regression tasks. This study is concerned with investigating the Total Electron Content (TEC) time series by using a SVM to detect seismo-ionospheric anomalous variations induced by the three powerful earthquakes of Tohoku (11 March 2011), Haiti (12 January 2010) and Samoa (29 September 2009). The duration of TEC time series dataset is 49, 46 and 71 days, for Tohoku, Haiti and Samoa earthquakes, respectively, with each at time resolution of 2 h. In the case of Tohoku earthquake, the results show that the difference between the predicted value obtained from the SVM method and the observed value reaches the maximum value (i.e., 129.31 TECU) at earthquake time in a period of high geomagnetic activities. The SVM method detected a considerable number of anomalous occurrences 1 and 2 days prior to the Haiti earthquake and also 1 and 5 days before the Samoa earthquake in a period of low geomagnetic activities. In order to show that the method is acting sensibly with regard to the results extracted during nonevent and event TEC data, i.e., to perform some null-hypothesis tests in which the methods would also be calibrated, the same period of data from the previous year of the Samoa earthquake date has been taken into the account. Further to this, in this study, the detected TEC anomalies using the SVM method were compared to the previous results (Akhoondzadeh and Saradjian, 2011; Akhoondzadeh, 2012) obtained from the mean, median, wavelet and Kalman filter methods. The SVM detected anomalies are similar to those detected using the previous methods. It can be concluded that SVM can be a suitable learning method to detect the novelty changes of a nonlinear time series such as variations of earthquake precursors.

Description: Magnetosheath dynamic pressure enhancements: occurrence and typical properties Annales Geophysicae, 31, 319-331, 2013 Author(s): M. O. Archer and T. S. Horbury The first comprehensive statistical study of large-amplitude (〉 100%) transient enhancements of the magnetosheath dynamic pressure reveals events of up to ~ 15 times the ambient dynamic pressure with durations up to 3 min and an average duration of around 30 s, predominantly downstream of the quasi-parallel shock. The dynamic pressure transients are most often dominated by velocity increases along with a small fractional increase in the density, though the velocity is generally only deflected by a few degrees. Superposed wavelet transforms of the magnetic field show that, whilst most enhancements exhibit changes in the magnetosheath magnetic field, the majority are not associated with changes in the Interplanetary Magnetic Field (IMF). However, there is a minority of enhancements that do appear to be associated with solar wind discontinuities which cannot be explained simply by random events. In general, it is found that during periods of magnetosheath dynamic pressure enhancements the IMF is steadier than usual. This suggests that a stable foreshock and hence foreshock structures or processes may be important in the generation of the majority of magnetosheath dynamic pressure enhancements.

Description: New type of ensemble of quasi-periodic, long-lasting VLF emissions at the auroral zone Annales Geophysicae, 30, 1655-1660, 2012 Author(s): J. Manninen, N. G. Kleimenova, and O. V. Kozyreva A new type of the series of quasi-periodic (QP) very low frequency (VLF) emissions in frequency range of 1–5 kHz, and not associated with geomagnetic pulsations, has been discovered at auroral latitudes ( L = 5.3) during the Finnish VLF campaign (held in December 2011). At least five unusually spectacular events, each with a duration of several hours, have been observed during the night under conditions of quiet geomagnetic activity (Kp = 0–1), although QPs usually occur during the daytime. Contrary to the QP emissions typically occurring during the day, the spectral structure of these QP events represented an extended, complicated sequence of repeated discrete rising VLF signals. Their duration was about 2–3 min each, with the repetition periods ranging from ~1 min to ~10 min. Two such nighttime non-typical events are reported in this paper. The fine structure of the separated QP elements may represent a mixture of the different frequency band signals, which seem to have independent origins. It was found that the periodic signals with lower frequency appear to trigger the strong dispersive upper frequency signals. The temporal dynamics of the spectral structure of the QPs studied were significantly controlled by some disturbances in the solar wind and interplanetary magnetic field (IMF). This finding is very important for future theoretical investigations because the generation mechanism of this new type of QP emissions is not yet understood.

Description: Spatial characteristics of wave-like structures in diffuse aurora obtained using optical observations Annales Geophysicae, 30, 1693-1701, 2012 Author(s): K. Axelsson, T. Sergienko, H. Nilsson, U. Brändström, Y. Ebihara, K. Asamura, and M. Hirahara We present the results of a statistical study using optical images from ALIS (Auroral Large Imaging System) to investigate the spatial and temporal variations of structures in diffuse aurora. Analysis of conjugate Reimei data shows that such fine structures are a result of modulation of high-energy precipitating electrons. Pitch angle diffusion into the loss cone due to interaction of whistler mode waves with plasma sheet electrons is the most feasible mechanism leading to high-energy electron precipitation. This suggests that the fine structure is an indication of modulations of the efficiency of the wave–particle interaction. The scale sizes and variations of these structures, mapped to the magnetosphere, can give us information about the characteristics of the modulating wave activity. We found the scale size of the auroral stripes and the spacing between them to be on average 13–14 km, which corresponds to 3–4 ion gyro radii for protons with an energy of 7 keV. The structures move southward with a speed close to zero in the plasma convection frame.

Description: Response of polar mesosphere summer echoes to geomagnetic disturbances in the Southern and Northern Hemispheres: the importance of nitric oxide Annales Geophysicae, 31, 333-347, 2013 Author(s): S. Kirkwood, E. Belova, P. Dalin, M. Mihalikova, D. Mikhaylova, D. Murtagh, H. Nilsson, K. Satheesan, J. Urban, and I. Wolf The relationship between polar mesosphere summer echoes (PMSE) and geomagnetic disturbances (represented by magnetic K indices) is examined. Calibrated PMSE reflectivities for the period May 2006–February 2012 are used from two 52.0/54.5 MHz radars located in Arctic Sweden (68° N, geomagnetic latitude 65°) and at two different sites in Queen Maud Land, Antarctica (73°/72° S, geomagnetic latitudes 62°/63°). In both the Northern Hemisphere (NH) and the Southern Hemisphere (SH) there is a strong increase in mean PMSE reflectivity between quiet and disturbed geomagnetic conditions. Mean volume reflectivities are slightly lower at the SH locations compared to the NH, but the position of the peak in the lognormal distribution of PMSE reflectivities is close to the same at both NH and SH locations, and varies only slightly with magnetic disturbance level. Differences between the sites, and between geomagnetic disturbance levels, are primarily due to differences in the high-reflectivity tail of the distribution. PMSE occurrence rates are essentially the same at both NH and SH locations during most of the PMSE season when a sufficiently low detection threshold is used so that the peak in the lognormal distribution is included. When the local-time dependence of the PMSE response to geomagnetic disturbance level is considered, the response in the NH is found to be immediate at most local times, but delayed by several hours in the afternoon sector and absent in the early evening. At the SH sites, at lower magnetic latitude, there is a delayed response (by several hours) at almost all local times. At the NH (auroral zone) site, the dependence on magnetic disturbance is highest during evening-to-morning hours. At the SH (sub-auroral) sites the response to magnetic disturbance is weaker but persists throughout the day. While the immediate response to magnetic activity can be qualitatively explained by changes in electron density resulting from energetic particle precipitation, the delayed response can largely be explained by changes in nitric oxide concentrations. Observations of nitric oxide concentration at PMSE heights by the Odin satellite support this hypothesis. Sensitivity to geomagnetic disturbances, including nitric oxide produced during these disturbances, can explain previously reported differences between sites in the auroral zone and those at higher or lower magnetic latitudes. The several-day lifetime of nitric oxide can also explain earlier reported discrepancies between high correlations for average conditions (year-by-year PMSE reflectivities and K indices) and low correlations for minute-to-day timescales.

Description: Observations of poleward-propagating large-scale traveling ionospheric disturbances in southern China Annales Geophysicae, 31, 377-385, 2013 Author(s): F. Ding, W. Wan, B. Ning, B. Zhao, Q. Li, Y. Wang, L. Hu, R. Zhang, and B. Xiong We report here on two cases of poleward-propagating large-scale traveling ionospheric disturbances (LSTIDs) in China during a medium-scale storm between 27 May and 1 June 2011. The observations were conducted by making use of the Global Positioning System network and ionosondes in China and Southeast Asia. One northeastward-propagating LSTID occurred on the morning of 30 May, while the other was observed during the nighttime of 1 June. Both poleward-traveling LSTIDs occurred during the storm's recovery phase in southern China's low-latitude region (geomagnetic latitude ~ 7.3–24° N) and experienced severe dissipation during their propagation from south to north. Although the initial relative amplitude of the nighttime LSTID was ~ 60% larger than that of the morning event, the nighttime event dissipated more quickly than the morning event because of a strong nighttime enhancement in background total electronic content (TEC) during storm time, which led to strong ion-drag dissipation during the evening. The poleward-propagating LSTIDs exhibit a narrower latitudinal range, a smaller amplitude, and a slightly higher elevation compared with the equatorward-moving LSTIDs observed in the same region. Given these features, the poleward-propagating LSTIDs were likely excited by some local source near southern China. Excitation of secondary LSTIDs during the dissipation of some primary medium-scale disturbances from the lower atmosphere is a possible mechanism.

Description: A new method for solving the MHD equations in the magnetosheath Annales Geophysicae, 31, 419-437, 2013 Author(s): C. Nabert, K.-H. Glassmeier, and F. Plaschke We present a new analytical method to derive steady-state magnetohydrodynamic (MHD) solutions of the magnetosheath in different levels of approximation. With this method, we calculate the magnetosheath's density, velocity, and magnetic field distribution as well as its geometry. Thereby, the solution depends on the geomagnetic dipole moment and solar wind conditions only. To simplify the representation, we restrict our model to northward IMF with the solar wind flow along the stagnation streamline. The sheath's geometry, with its boundaries, bow shock and magnetopause, is determined self-consistently. Our model is stationary and time relaxation has not to be considered as in global MHD simulations. Our method uses series expansion to transfer the MHD equations into a new set of ordinary differential equations. The number of equations is related to the level of approximation considered including different physical processes. These equations can be solved numerically; however, an analytical approach for the lowest-order approximation is also presented. This yields explicit expressions, not only for the flow and field variations but also for the magnetosheath thickness, depending on the solar wind parameters. Results are compared to THEMIS data and offer a detailed explanation of, e.g., the pile-up process and the corresponding plasma depletion layer, the bow shock and magnetopause geometry, the magnetosheath thickness, and the flow deceleration.

Description: Investigation of a mesospheric bore event over northern China Annales Geophysicae, 31, 409-418, 2013 Author(s): Q. Li, J. Xu, J. Yue, X. Liu, W. Yuan, B. Ning, S. Guan, and J. P. Younger A mesospheric bore event was observed using an OH all-sky airglow imager (ASAI) at Xinglong (40.2° N, 117.4° E), in northern China, on the night of 8–9 January 2011. Simultaneous observations by a Doppler meteor radar, a broadband sodium lidar, and TIMED/SABER OH intensity and temperature measurements are used to investigate the characteristics and environment of the bore propagation and the possible relations with the Na density perturbations. The bore propagated from northeast to southwest and divided the sky into bright and dark halves. The calculations show that the bore has an average phase velocity of 68 m s −1 . The crests following the bore have a horizontal wavelength of ~ 22 km. These parameters are consistent with the hydraulic jump theory proposed by Dewan and Picard, as well as the previous bore reports. Simultaneous wind measurements from the Doppler meteor radar at Shisanling (40.3° N, 116.2° E) and temperature data from SABER on board the TIMED satellite are used to characterize the propagating environment of the bore. The result shows that a thermal-Doppler duct exists near the OH layer that supports the horizontal propagation of the bore. Simultaneous Na lidar observations at Yanqing (40.4° N, 116.0° E) suggest that there is a downward displacement of Na density during the passage of the mesospheric bore event.

Description: Near real-time estimation of water vapour in the troposphere using ground GNSS and the meteorological data Annales Geophysicae, 30, 1379-1391, 2012 Author(s): J. Bosy, J. Kaplon, W. Rohm, J. Sierny, and T. Hadas The near real-time (NRT) high resolution water vapour distribution models can be constructed based on GNSS observations delivered from Ground Base Augmentation Systems (GBAS) and ground meteorological data. Since 2008 in the territory of Poland, a GBAS system called ASG-EUPOS (Active Geodetic Network) has been operating. This paper addresses the problems concerning construction of the NRT model of water vapour distribution in the troposphere near Poland. The first section presents all available GNSS and ground meteorological stations in the area of Poland and neighbouring countries. In this section, data feeding scheme is discussed, together with timeline and time resolution. The high consistency between measured and interpolated temperature value is shown, whereas some discrepancy in the pressure is observed. In the second section, the NRT GNSS data processing strategy of ASG-EUPOS network is discussed. Preliminary results show fine alignment of the obtained Zenith Troposphere Delays (ZTDs) with reference data from European Permanent Network (EPN) processing center. The validation of NRT troposphere products against daily solution shows 15 mm standard deviation of obtained ZTD differences. The last section presents the first results of 2-D water vapour distribution above the GNSS network and application of the tomographic model to 3-D distribution of water vapour in the atmosphere. The GNSS tomography model, working on the simulated data from numerical forecast model, shows high consistency with the reference data (by means of standard deviation 4 mm km −1 or 4 ppm), however, noise analysis shows high solution sensitivity to errors in observations. The discrepancy for real data preliminary solution (measured as a mean standard deviation) between reference NWP data and tomography data was on the level of 9 mm km −1 (or 9 ppm) in terms of wet refractivity.

Description: Variations of the ionospheric TEC using simultaneous measurements from the China Crustal Movement Observation Network Annales Geophysicae, 30, 1423-1433, 2012 Author(s): Y. W. Wu, R. Y. Liu, B. C. Zhang, Z. S. Wu, J. S. Ping, J. M. Liu, and Z. J. Hu Variations of the ionospheric Total Electron Content (TEC) over China are investigated using the TEC data obtained from China Crustal Movement Observation Network in the year 2004. The results show a single-peak occurred in post-noon for the diurnal variation and two peaks exit around two equinox points, respectively, for the seasonal variation. Overall, the values of TEC increased from the north to the south of China. There were small but clear longitudinal differences in both sides of the longitudes with zero magnetic declination. The intensity of the day-to-day variation of TEC was not a monotonic change along the latitudes. It was usually weaker in the middle of China than that in the north or south. Comparing with the maximum F-layer electron density ( Nm F2) derived from the ionosonde stations in China, it is found that the day-to-day variation of TEC was less significant than that of Nm F2, and that the northern crest of the equatorial anomaly identified from the Nm F2 data can reach Guangzhou-region. While, the TEC crest was hardly observed in the same location. This is probably caused by the tilt of topside ionosphere near the northern anomaly crest region at lower latitudes.

Description: Collisionless reconnection: magnetic field line interaction Annales Geophysicae, 30, 1515-1528, 2012 Author(s): R. A. Treumann, W. Baumjohann, and W. D. Gonzalez Magnetic field lines are quantum objects carrying one quantum Φ 0 = 2π h/e of magnetic flux and have finite radius λ m . Here we argue that they possess a very specific dynamical interaction. Parallel field lines reject each other. When confined to a certain area they form two-dimensional lattices of hexagonal structure. We estimate the filling factor of such an area. Anti-parallel field lines, on the other hand, attract each other. We identify the physical mechanism as being due to the action of the gauge potential field, which we determine quantum mechanically for two parallel and two anti-parallel field lines. The distortion of the quantum electrodynamic vacuum causes a cloud of virtual pairs. We calculate the virtual pair production rate from quantum electrodynamics and estimate the virtual pair cloud density, pair current and Lorentz force density acting on the field lines via the pair cloud. These properties of field line dynamics become important in collisionless reconnection, consistently explaining why and how reconnection can spontaneously set on in the field-free centre of a current sheet below the electron-inertial scale.

Description: Observation of a tropopause fold by MARA VHF wind-profiler radar and ozonesonde at Wasa, Antarctica: comparison with ECMWF analysis and a WRF model simulation Annales Geophysicae, 30, 1411-1421, 2012 Author(s): M. Mihalikova, S. Kirkwood, J. Arnault, and D. Mikhaylova Tropopause folds are one of the mechanisms of stratosphere–troposphere exchange, which can bring ozone rich stratospheric air to low altitudes in the extra-tropical regions. They have been widely studied at northern mid- or high latitudes, but so far almost no studies have been made at mid- or high southern latitudes. The Moveable Atmospheric Radar for Antarctica (MARA), a 54.5 MHz wind-profiler radar, has operated at the Swedish summer station Wasa, Antarctica (73° S, 13.5° W) during austral summer seasons from 2007 to 2011 and has observed on several occasions signatures similar to those caused by tropopause folds at comparable Arctic latitudes. Here a case study is presented of one of these events when an ozonesonde successfully sampled the fold. Analysis from European Center for Medium Range Weather Forecasting (ECMWF) is used to study the circumstances surrounding the event, and as boundary conditions for a mesoscale simulation using the Weather Research and Forecasting (WRF) model. The fold is well resolved by the WRF simulation, and occurs on the poleward side of the polar jet stream. However, MARA resolves fine-scale layering associated with the fold better than the WRF simulation.

Description: Chorus wave-normal statistics in the Earth's radiation belts from ray tracing technique Annales Geophysicae, 30, 1223-1233, 2012 Author(s): H. Breuillard, Y. Zaliznyak, V. Krasnoselskikh, O. Agapitov, A. Artemyev, and G. Rolland Discrete ELF/VLF (Extremely Low Frequency/Very Low Frequency) chorus emissions are one of the most intense electromagnetic plasma waves observed in radiation belts and in the outer terrestrial magnetosphere. These waves play a crucial role in the dynamics of radiation belts, and are responsible for the loss and the acceleration of energetic electrons. The objective of our study is to reconstruct the realistic distribution of chorus wave-normals in radiation belts for all magnetic latitudes. To achieve this aim, the data from the electric and magnetic field measurements onboard Cluster satellite are used to determine the wave-vector distribution of the chorus signal around the equator region. Then the propagation of such a wave packet is modeled using three-dimensional ray tracing technique, which employs K. Rönnmark's WHAMP to solve hot plasma dispersion relation along the wave packet trajectory. The observed chorus wave distributions close to waves source are first fitted to form the initial conditions which then propagate numerically through the inner magnetosphere in the frame of the WKB approximation. Ray tracing technique allows one to reconstruct wave packet properties (electric and magnetic fields, width of the wave packet in k -space, etc.) along the propagation path. The calculations show the spatial spreading of the signal energy due to propagation in the inhomogeneous and anisotropic magnetized plasma. Comparison of wave-normal distribution obtained from ray tracing technique with Cluster observations up to 40° latitude demonstrates the reliability of our approach and applied numerical schemes.

Description: High-precision measurement of satellite velocity using the EISCAT radar Annales Geophysicae, 30, 1555-1565, 2012 Author(s): T. Nygrén, J. Markkanen, A. Aikio, and M. Voiculescu This paper presents a method of measuring the velocity of a hard target using radar pulses reflected from the target flying through the radar beam. The method has two stages. First, the Doppler shifts of the echo pulses are calculated at a high accuracy with an algorithm which largely improves the accuracy given by the Fourier transform. The algorithm also calculates the standard deviations of the Doppler frequencies with Monte Carlo simulation. The second step is to fit the results from a sequence of radar pulses to a velocity model allowing linear variation of the second time derivative of target range. The achieved accuracies are demonstrated using radio pulses reflected by a satellite passing through the beam of the EISCAT UHF radar working at 930-MHz frequency. At high SNR levels, the standard deviations of the frequency from a single pulse reach typically down to 0.2 Hz. The best standard deviations of velocity fit are below 5 mm s −1 while those of the second time derivative of range are below 1 cm s −2 .

Description: Equatorial plasma bubbles and L-band scintillations in Africa during solar minimum Annales Geophysicae, 30, 675-682, 2012 Author(s): V. V. Paznukhov, C. S. Carrano, P. H. Doherty, K. M. Groves, R. G. Caton, C. E. Valladares, G. K. Seemala, C. T. Bridgwood, J. Adeniyi, L. L. N. Amaeshi, B. Damtie, F. D'Ujanga Mutonyi, J. O. H. Ndeda, P. Baki, O. K. Obrou, B. Okere, and G. M. Tsidu We report on the longitudinal, local time and seasonal occurrence of equatorial plasma bubbles (EPBs) and L band (GPS) scintillations over equatorial Africa. The measurements were made in 2010, as a first step toward establishing the climatology of ionospheric irregularities over Africa. The scintillation intensity is obtained by measuring the standard deviation of normalized GPS signal power. The EPBs are detected using an automated technique, where spectral analysis is used to extract and identify EPB events from the GPS TEC measurements. Overall, the observed seasonal climatology of the EPBs as well as GPS scintillations in equatorial Africa is adequately explained by geometric arguments, i.e., by the alignment of the solar terminator and local geomagnetic field, or STBA hypothesis (Tsunoda, 1985, 2010a). While plasma bubbles and scintillations are primarily observed during equinoctial periods, there are longitudinal differences in their seasonal occurrence statistics. The Atlantic sector has the most intense, longest lasting, and highest scintillation occurrence rate in-season. There is also a pronounced increase in the EPB occurrence rate during the June solstice moving west to east. In Africa, the seasonal occurrence shifts towards boreal summer solstice, with fewer occurrences and shorter durations in equinox seasons. Our results also suggest that the occurrence of plasma bubbles and GPS scintillations over Africa are well correlated, with scintillation intensity depending on depletion depth. A question remains about the possible physical mechanisms responsible for the difference in the occurrence phenomenology of EPBs and GPS scintillations between different regions in equatorial Africa.

Description: The ultra-fast Kelvin waves in the equatorial ionosphere: observations and modeling Annales Geophysicae, 31, 209-215, 2013 Author(s): A. N. Onohara, I. S. Batista, and H. Takahashi The main purpose of this study is to investigate the vertical coupling between the mesosphere and lower thermosphere (MLT) region and the ionosphere through ultra-fast Kelvin (UFK) waves in the equatorial atmosphere. The effect of UFK waves on the ionospheric parameters was estimated using an ionospheric model which calculates electrostatic potential in the E-region and solves coupled electrodynamics of the equatorial ionosphere in the E- and F-regions. The UFK wave was observed in the South American equatorial region during February–March 2005. The MLT wind data obtained by meteor radar at São João do Cariri (7.5° S, 37.5° W) and ionospheric F-layer bottom height ( h 'F) observed by ionosonde at Fortaleza (3.9° S; 38.4° W) were used in order to calculate the wave characteristics and amplitude of oscillation. The simulation results showed that the combined electrodynamical effect of tides and UFK waves in the MLT region could explain the oscillations observed in the ionospheric parameters.

Description: Payload charging events in the mesosphere and their impact on Langmuir type electric probes Annales Geophysicae, 31, 187-196, 2013 Author(s): T. A. Bekkeng, A. Barjatya, U.-P. Hoppe, A. Pedersen, J. I. Moen, M. Friedrich, and M. Rapp Three sounding rockets were launched from Andøya Rocket Range in the ECOMA campaign in December 2010. The aim was to study the evolution of meteoric smoke particles during a major meteor shower. Of the various instruments onboard the rocket payload, this paper presents the data from a multi-Needle Langmuir Probe (m-NLP) and a charged dust detector. The payload floating potential, as observed using the m-NLP instrument, shows charging events on two of the three flights. These charging events cannot be explained using a simple charging model, and have implications towards the use of fixed bias Langmuir probes on sounding rockets investigating mesospheric altitudes. We show that for a reliable use of a single fixed bias Langmuir probe as a high spatial resolution relative density measurement, each payload should also carry an additional instrument to measure payload floating potential, and an instrument that is immune to spacecraft charging and measures absolute plasma density.

Description: Electric potential differences across auroral generator interfaces Annales Geophysicae, 31, 251-261, 2013 Author(s): J. De Keyser and M. Echim Strong localized high-altitude auroral electric fields, such as those observed by Cluster, are often associated with magnetospheric interfaces. The type of high-altitude electric field profile (monopolar, bipolar, or more complicated) depends on the properties of the plasmas on either side of the interface, as well as on the total electric potential difference across the structure. The present paper explores the role of this cross-field electric potential difference in the situation where the interface is a tangential discontinuity. A self-consistent Vlasov description is used to determine the equilibrium configuration for different values of the transverse potential difference. A major observation is that there exist limits to the potential difference, beyond which no equilibrium configuration of the interface can be sustained. It is further demonstrated how the plasma densities and temperatures affect the type of electric field profile in the transition, with monopolar electric fields appearing primarily when the temperature contrast is large. These findings strongly support the observed association of monopolar fields with the plasma sheet boundary. The role of shear flow tangent to the interface is also examined.

Description: Variability of ionospheric TEC during solar and geomagnetic minima (2008 and 2009): external high speed stream drivers Annales Geophysicae, 31, 263-276, 2013 Author(s): O. P. Verkhoglyadova, B. T. Tsurutani, A. J. Mannucci, M. G. Mlynczak, L. A. Hunt, and T. Runge We study solar wind–ionosphere coupling through the late declining phase/solar minimum and geomagnetic minimum phases during the last solar cycle (SC23) – 2008 and 2009. This interval was characterized by sequences of high-speed solar wind streams (HSSs). The concomitant geomagnetic response was moderate geomagnetic storms and high-intensity, long-duration continuous auroral activity (HILDCAA) events. The JPL Global Ionospheric Map (GIM) software and the GPS total electron content (TEC) database were used to calculate the vertical TEC (VTEC) and estimate daily averaged values in separate latitude and local time ranges. Our results show distinct low- and mid-latitude VTEC responses to HSSs during this interval, with the low-latitude daytime daily averaged values increasing by up to 33 TECU (annual average of ~20 TECU) near local noon (12:00 to 14:00 LT) in 2008. In 2009 during the minimum geomagnetic activity (MGA) interval, the response to HSSs was a maximum of ~30 TECU increases with a slightly lower average value than in 2008. There was a weak nighttime ionospheric response to the HSSs. A well-studied solar cycle declining phase interval, 10–22 October 2003, was analyzed for comparative purposes, with daytime low-latitude VTEC peak values of up to ~58 TECU (event average of ~55 TECU). The ionospheric VTEC changes during 2008–2009 were similar but ~60% less intense on average. There is an evidence of correlations of filtered daily averaged VTEC data with Ap index and solar wind speed. We use the infrared NO and CO 2 emission data obtained with SABER on TIMED as a proxy for the radiation balance of the thermosphere. It is shown that infrared emissions increase during HSS events possibly due to increased energy input into the auroral region associated with HILDCAAs. The 2008–2009 HSS intervals were ~85% less intense than the 2003 early declining phase event, with annual averages of daily infrared NO emission power of ~ 3.3 × 10 10 W and 2.7 × 10 10 W in 2008 and 2009, respectively. The roles of disturbance dynamos caused by high-latitude winds (due to particle precipitation and Joule heating in the auroral zones) and of prompt penetrating electric fields (PPEFs) in the solar wind–ionosphere coupling during these intervals are discussed. A correlation between geoeffective interplanetary electric field components and HSS intervals is shown. Both PPEF and disturbance dynamo mechanisms could play important roles in solar wind–ionosphere coupling during prolonged (up to days) external driving within HILDCAA intervals.

Description: Case study of stratospheric gravity waves of convective origin over Arctic Scandinavia – VHF radar observations and numerical modelling Annales Geophysicae, 31, 239-250, 2013 Author(s): A. Réchou, J. Arnault, P. Dalin, and S. Kirkwood Orography is a well-known source of gravity and inertia-gravity waves in the atmosphere. Other sources, such as convection, are also known to be potentially important but the large amplitude of orographic waves over Scandinavia has generally precluded the possibility to study such other sources experimentally in this region. In order to better understand the origin of stratospheric gravity waves observed by the VHF radar ESRAD (Esrange MST radar) over Kiruna, in Arctic Sweden (67.88° N, 21.10° E), observations have been compared to simulations made using the Weather Research and Forecasting model (WRF) with and without the effects of orography and clouds. This case study concerns gravity waves observed from 00:00 UTC on 18 February to 12:00 UTC on 20 February 2007. We focus on the wave signatures in the static stability field and vertical wind deduced from the simulations and from the observations as these are the parameters which are provided by the observations with the best height coverage. As is common at this site, orographic gravity waves were produced over the Scandinavian mountains and observed by the radar. However, at the same time, southward propagation of fronts in the Barents Sea created short-period waves which propagated into the stratosphere and were transported, embedded in the cyclonic winds, over the radar site.

Description: Long-term trends in the ionospheric F2 region with different solar activity indices Annales Geophysicae, 31, 291-303, 2013 Author(s): J. Mielich and J. Bremer A new comprehensive data collection by Damboldt and Suessmann (2012a) with monthly fo F2 and M(3000)F2 median values is an excellent basis for the derivation of long-term trends in the ionospheric F2 region. Ionospheric trends have been derived only for stations with data series of at least 22 years (124 stations with fo F2 data and 113 stations with M(3000)F2 data) using a twofold regression analysis depending on solar and geomagnetic activity. Three main results have been derived: Firstly, it could be shown that the solar 10.7 cm radio flux F10.7 is a better index for the description of the solar activity than the relative solar sunspot number R as well as the solar EUV proxy E10.7. Secondly, the global mean fo F2 and hm F2 trends derived for the interval between 1948 and 2006 are in surprisingly good agreement with model calculations of an increasing atmospheric greenhouse effect (Rishbeth and Roble, 1992). Thirdly, during the years 2007 until 2009, the hm F2 values and to a smaller amount the fo F2 values strongly decrease. The reason for this effect is a reduction of the thermospheric density and ionization due to a markedly reduced solar EUV irradiation and extremely small geomagnetic activity during the solar cycle 23/24 minimum.

Description: Simultaneous measurement of OI 557.7 nm, O 2 (0, 1) Atmospheric Band and OH (6, 2) Meinel Band nightglow at Kolhapur (17° N), India Annales Geophysicae, 31, 197-208, 2013 Author(s): N. Parihar, A. Taori, S. Gurubaran, and G. K. Mukherjee Near-simultaneous measurements of OI 557.7 nm, O 2 (0, 1) Atmospheric Band and OH (6, 2) Meinel Band nightglow were carried out at Kolhapur (17° N), India during February–March 2007. Atmospheric temperatures around 87 and 94 km were derived from the knowledge of intensity measurements of spectral features OH (6, 2) Meinel Band and O 2 Atmospheric Band, respectively. An account of the behaviour of derived temperatures has been presented. The nocturnal behaviour of OH and O 2 temperatures is governed by the waves of tidal origin, whereas the signatures of planetary wave-like oscillations is noted in the night-to-night variation of two temperatures. This is probably the first report of planetary waves observed in nightglow temperatures from the Indian subcontinent.

Description: Refilling process in the plasmasphere: a 3-D statistical characterization based on Cluster density observations Annales Geophysicae, 31, 217-237, 2013 Author(s): G. Lointier, F. Darrouzet, P. M. E. Décréau, X. Vallières, S. Kougblénou, J. G. Trotignon, and J.-L. Rauch The Cluster mission offers an excellent opportunity to investigate the evolution of the plasma population in a large part of the inner magnetosphere, explored near its orbit's perigee, over a complete solar cycle. The WHISPER sounder, on board each satellite of the mission, is particularly suitable to study the electron density in this region, between 0.2 and 80 cm −3 . Compiling WHISPER observations during 1339 perigee passes distributed over more than three years of the Cluster mission, we present first results of a statistical analysis dedicated to the study of the electron density morphology and dynamics along and across magnetic field lines between L = 2 and L = 10. In this study, we examine a specific topic: the refilling of the plasmasphere and trough regions during extended periods of quiet magnetic conditions. To do so, we survey the evolution of the ap index during the days preceding each perigee crossing and sort out electron density profiles along the orbit according to three classes, namely after respectively less than 2 days, between 2 and 4 days, and more than 4 days of quiet magnetic conditions (ap ≤ 15 nT) following an active episode (ap 〉 15 nT). This leads to three independent data subsets. Comparisons between density distributions in the 3-D plasmasphere and trough regions at the three stages of quiet magnetosphere provide novel views about the distribution of matter inside the inner magnetosphere during several days of low activity. Clear signatures of a refilling process inside an expended plasmasphere in formation are noted. A plasmapause-like boundary, at L ~ 6 for all MLT sectors, is formed after 3 to 4 days and expends somewhat further after that. In the outer part of the plasmasphere ( L ~ 8), latitudinal profiles of median density values vary essentially according to the MLT sector considered rather than according to the refilling duration. The shape of these density profiles indicates that magnetic flux tubes are not fully replenished after 6 days of quiet conditions. In addition, the outer plasmasphere in the night and dawn sectors (22:00 to 10:00 MLT range) maintains an overall clear deficit of ionospheric population, when compared to the situation in the noon and dusk sectors (10:00 to 22:00 MLT range).

Description: Spatiotemporal variability analysis of diffuse radiation in China during 1981–2010 Annales Geophysicae, 31, 277-289, 2013 Author(s): X. L. Ren, H. L. He, L. Zhang, L. Zhou, G. R. Yu, and J. W. Fan Solar radiation is the primary driver of terrestrial plant photosynthesis and the diffuse component can enhance canopy light use efficiency (LUE), which in turn influences the carbon balance of terrestrial ecosystems. In this study we calculated the spatial data of diffuse radiation in China from 1981 to 2010, using a radiation decomposition model and spatial interpolation method based on observational data. Furthermore, we explored the spatiotemporal characteristics of diffuse radiation using GIS and trend analysis techniques. The results show the following: (1) The spatial patterns of perennial average of annual diffuse radiation during 1981–2010 are complex and inhomogeneous in China, generally lower in the north and higher in the south and west. The perennial average ranges from 1730.20 to 3064.41 MJ m −2 yr −1 across the whole country. (2) There is an increasing trend of annual diffuse radiation in China from 1981 to 2010 on the whole, with mean increasing amplitude of 7.03 MJ m −2 yr −1 per decade. Whereas a significant downtrend was observed in the first 10 years, distinct anomalies in 1982, 1983, 1991 and 1992 occurred due to the eruptions of El Chinchon and Pinatubo. (3) The spatial distribution of the temporal variability of diffuse radiation showed significant regional heterogeneity in addition to the seasonal differences. Northwestern China has the most evident downtrend, with highest decreasing rate of 6% per decade, while the Tibetan Plateau has the most evident uptrend, with highest increasing rate of up to 9% per decade. Such quantitative spatiotemporal characteristics of diffuse radiation are essential in regional scale modeling of terrestrial carbon dynamics.

Description: Comment on "Storming the Bastille: the effect of electric fields on the ionospheric F-layer" by Rishbeth et al. (2010) Annales Geophysicae, 31, 145-150, 2013 Author(s): B. T. Tsurutani, A. J. Mannuccci, O. P. Verkhoglyadova, and G. S. Lakhina No abstract available.

Description: Geosynchronous magnetic field responses to fast solar wind dynamic pressure enhancements: MHD field model Annales Geophysicae, 30, 1285-1295, 2012 Author(s): T. R. Sun, C. Wang, N. L. Borodkova, and G. N. Zastenker We performed global MHD simulations of the geosynchronous magnetic field in response to fast solar wind dynamic pressure ( P d ) enhancements. Taking three P d enhancement events in 2000 as examples, we found that the main features of the total field B and the dominant component B z can be efficiently predicted by the MHD model. The predicted B and B z varies with local time, with the highest level near noon and a slightly lower level around mid-night. However, it is more challenging to accurately predict the responses of the smaller component at the geosynchronous orbit (i.e., B x and B y ). In contrast, the limitations of T01 model in predicting responses to fast P d enhancements are presented.

Description: Low-frequency magnetic field fluctuations in Earth's plasma environment observed by THEMIS Annales Geophysicae, 30, 1271-1283, 2012 Author(s): L. Guicking, K.-H. Glassmeier, H.-U. Auster, Y. Narita, and G. Kleindienst Low-frequency magnetic wave activity in Earth's plasma environment was determined based on a statistical analysis of THEMIS magnetic field data. We observe that the spatial distribution of low-frequency magnetic field fluctuations reveals highest values in the magnetosheath, but the observations differ qualitatively from observations at Venus presented in a previous study since significant wave activity at Earth is also observed in the nightside magnetosheath. Outside the magnetosheath the low-frequency wave activity level is generally very low. By means of an analytical streamline model for the magnetosheath plasma flow, we are able to investigate the spatial and temporal evolution of wave intensity along particular streamlines in order to characterise possible wave generation mechanisms. We observe a decay of wave intensity along the streamlines, but contrary to the situation at Venus, we obtain good qualitative agreement with the theoretical concept of freely evolving/decaying turbulence. Differences between the dawn region and the dusk region can be observed only further away from the magnetopause. We conclude that wave generation mechanisms may be primarily attributed to processes at or in the vicinity of the bow shock. The difference with the observations of the Venusian magnetosheath we interpret to be the result of the different types of solar wind interaction processes since the Earth possesses a global magnetic field while Venus does not, and therefore the observed magnetic wave activities may be caused by diverse magnetic field controlled characteristics of wave generation processes.

Description: Variations of ionospheric profile parameters during solar maximum and comparison with IRI-2007 over Chung-Li, Taiwan Annales Geophysicae, 30, 1249-1257, 2012 Author(s): Y. J. Chuo This paper studies the seasonal changes in the diurnal variation of ionospheric bottomside slab thickness (B0), based on observations during high solar activities at the equatorial ionization anomaly (EIA) area station of Chung-Li (121.1° E, 24.9° N), Taiwan. The data examined in this investigation are derived from ionograms recorded at Chung-Li in 1999, and are compared with International Reference Ionosphere (IRI-2007) model values. In our data set B0 shows largest values and biggest changes during the daytime (06:00–12:00 LT) particularly in the summer. Moreover, the diurnal variation of B0 shows an abnormal peak during the pre-sunrise period, especially in the winter. The variation in the F-peak height ( hm F2) is related to a thermospheric wind traveling toward the equator, which also enhances B0 during the pre-sunrise period. The results of the comparison with the IRI model show that B0 is overestimated, in both the B0-table and the Gulyaeva option, after noon LT in the equinox (B0-table and Gulyaeva average values for the overestimation are 11 and 47 km, respectively) and summer (B0-table and Gulyaeva average values for the overestimation are 23 and 71 km, respectively) periods. Furthermore, the modeled values are underestimated at approximately 31 and 14 km for the table and Gulyaeva option during the daytime in the winter, respectively. The F2-layer maximum electron density ( Nm F2) data show reasonably favorable agreement with the model for a high correlation coefficient of approximately 0.97, with the major difference observed at approximately noon in the equinox and winter seasons. Regarding the hm F2 data, the model shows agreement with the observed values, and the largest discrepancy (average value is 39 km) was observed in the summer and the smallest (average value is 11 km) in the equinox season. This paper provides a comprehensive discussion on the relationship among B0, the Nm F2 and the hm F2 for geomagnetic storm events.

Description: Features of highly structured equatorial plasma irregularities deduced from CHAMP observations Annales Geophysicae, 30, 1259-1269, 2012 Author(s): C. Xiong, H. Lühr, S. Y. Ma, C. Stolle, and B. G. Fejer In this study five years of CHAMP (Challenging Mini-satellite Payload) fluxgate magnetometer (FGM) data is used to investigate the characteristics of Equatorial Plasma Bubbles (EPBs). We filtered the FGM data by using band-passes with four different cut-off periods to get the EPBs with different maximum spatial scale sizes in the meridional plane ranging from 76–608 km. Associated with the EPB observations at about 400 km, the typical altitude of CHAMP during the year 2000–2005, we also investigate the post-sunset equatorial vertical plasma drift data from ROCSAT-1 (Republic of China Satellite 1). Since the height of the F-layer is highly correlated with the vertical plasma drift and solar flux, we sorted the ROCSAT-1 data into different groups by F10.7. From the integrated vertical drift we have estimated the post-sunset uplift of the ionosphere. By comparing the properties of EPB occurrence for different scale sizes with the global distribution of plasma vertical uplift, we have found that EPBs reaching higher altitudes are more structured than those which are sampled by CHAMP near the top side of the depleted fluxtube. Such a result is in accord with 3-D model simulations (Aveiro and Hysell, 2010). Small-scale EPB structures are observed by CHAMP when the irregularities reach apex heights of 800 km and more. Such events are encountered primarily in the Brazilian sector during the months around November, when the post-sunset vertical plasma drift is high.

Description: Intensification of dayside diffuse auroral precipitation: contribution of dayside Whistler-mode chorus waves in realistic magnetic fields Annales Geophysicae, 30, 1297-1307, 2012 Author(s): R. Shi, D. Han, B. Ni, Z.-J. Hu, C. Zhou, and X. Gu Compared to the recently improved understanding of nightside diffuse aurora, the mechanism(s) responsible for dayside diffuse aurora remains poorly understood. While dayside chorus has been thought as a potential major contributor to dayside diffuse auroral precipitation, quantitative analyses of the role of chorus wave scattering have not been carefully performed. In this study we investigate a dayside diffuse auroral intensification event observed by the Chinese Arctic Yellow River Station (YRS) all-sky imagers (ASI) on 7 January 2005 and capture a substantial increase in diffuse auroral intensity at the 557.7 nm wavelength that occurred over almost the entire ASI field-of-view near 09:24 UT, i.e., ~12:24 MLT. Computation of bounce-averaged resonant scattering rates by dayside chorus emissions using realistic magnetic field models demonstrates that dayside chorus scattering can produce intense precipitation losses of plasma sheet electrons on timescales of hours (even approaching the strong diffusion limit) over a broad range of both energy and pitch angle, specifically, from ~1 keV to 50 keV with equatorial pitch angles from the loss cone to up to ~85° depending on electron energy. Subsequent estimate of loss cone filling index indicates that the loss cone can be substantially filled, due to dayside chorus driven pitch angle scattering, at a rate of ≥0.8 for electrons from ~500 eV to 50 keV that exactly covers the precipitating electrons for the excitation of green-line diffuse aurora. Estimate of electron precipitation flux at different energy levels, based on loss cone filling index profile and typical dayside electron distribution observed by THEMIS spacecraft under similar conditions, gives a total precipitation electron energy flux of the order of 0.1 erg cm −2 s −1 with ~1 keV characteristic energy (especially when using T01s), which can be very likely to cause intense green-line diffuse aurora activity on the dayside. Therefore, dayside chorus scattering in the realistic magnetic field can greatly contribute to the YRS ASI observed intensification of dayside green-line aurora. Besides wave induced scattering and changes in the ambient magnetic field, variations in associated electron flux can also contribute to enhanced diffuse aurora emissions, the possibility of which we cannot exactly rule out due to lack of simultaneous observations of magnetospheric particles. Since the geomagnetic activity level was rather low during the period of interest, it is reasonable to infer that changes in the associated electron flux in the magnetosphere should be small, and consequently its contribution to the observed enhanced diffuse auroral activity should be small as well. Our results support the scenario that dayside chorus could play a major role in the production of dayside diffuse aurora, and also demonstrate that changes in magnetospheric magnetic field should be considered to reasonably interpret observations of dayside diffuse aurora.

Description: Oxygen ion energization by waves in the high altitude cusp and mantle Annales Geophysicae, 30, 1309-1314, 2012 Author(s): M. Waara, H. Nilsson, R. Slapak, M. André, and G. Stenberg We present a comparative study of low frequency electric field spectral densities and temperatures observed by the Cluster spacecraft in the high altitude cusp/mantle region. We compare the relation between the O + temperature and wave intensity at the oxygen gyrofrequency at each measurement point and find a clear correlation. The trend of the correlation agrees with the predictions by both an asymptotic mean-particle theory and a test-particle approach. The perpendicular to parallel temperature ratio is also consistent with the predictions of the asymptotic mean-particle theory. At times the perpendicular temperature is significantly higher than predicted by the models. A simple study of the evolution of the particle distributions (conics) at these altitudes indicates that enhanced perpendicular temperatures would be observed over many R E after heating ceases. Therefore, sporadic intense heating is the likely explanation for cases with high temperature and comparably low wave activity. We observe waves of sufficient amplitude to explain the highest observed temperatures, while the theory in general overestimates the temperature associated with the highest observed wave activity, indicating that such high wave activity is very sporadic.

Description: The "step feature" of suprathermal ion distributions: a discriminator between acceleration processes? Annales Geophysicae, 30, 1315-1319, 2012 Author(s): H. J. Fahr and H. Fichtner The discussion of exactly which process is causing the preferred build-up of v −5 -power law tails of the velocity distribution of suprathermal particles in the solar wind is still ongoing. Criteria allowing one to discriminate between the various suggestions that have been made would be useful in order to clarify the physics behind these tails. With this study, we draw the attention to the so-called "step feature" of the velocity distributions and offer a criterion that allows one to distinguish between those scenarios that employ velocity diffusion, i.e. second-order Fermi processes, which are prime candidates in the present debate. With an analytical approximation to the self-consistently obtained velocity diffusion coefficient, we solve the transport equation for suprathermal particles. The numerical simulation reveals that this form of the diffusion coefficient naturally leads to the step feature of the velocity distributions. This finding favours – at least in regions of the appearance of the step feature (i.e. for heliocentric distances up to about 11 AU and at lower energies) – the standard velocity diffusion as a consequence of the particle's interactions with the plasma wave turbulence as opposed to that caused by velocity fluctuation-induced compressions and rarefactions.

Description: Plasma parameter analysis of the Langmuir decay process via Particle-in-Cell simulations Annales Geophysicae, 30, 1169-1183, 2012 Author(s): M. A. Diaz, M. Zettergren, J. L. Semeter, and M. Oppenheim The beam-plasma mechanism, based on the Langmuir decay process, has been proposed to explain naturally enhanced ion-acoustic lines (NEIALs), which are spectral distortions in incoherent scatter radar (ISR) data frequently observed in the vicinity of auroral arcs. In this work the effect of the Langmuir decay process on the ISR spectrum is studied and compared with an analytical model for different plasma parameters by using an electrostatic parallel particle-in-cell (EPPIC) code. Simulations show that the code is working in accordance with theory for a wide range of beam and plasma values and that the features of the spectrum are sensitive to changes of those values. These results suggest that the EPPIC code might be used to build a spectrum-plasma parameter model which will allow estimation of beam and plasma parameters from observed spectra. Simulations also confirm that background electron density ( n e ) plays an important role in determining the maximum detectable wavenumber of the enhancement. Specifically, results demonstrate that an increase in n e makes the enhancements of the ion acoustic more likely line at large wavenumbers, a finding consistent with statistical studies showing more frequent NEIAL occurrence near solar maximum. Finally, the simulations expose some inaccuracies of the current theoretical model in quantifying the energy passed from the beam to the Langmuir waves as well as with the range of enhanced wavenumbers. These differences may be attributable to the weak Langmuir turbulent regime assumption used in the theory.