scispace - formally typeset
Search or ask a question

Showing papers in "Annales Geophysicae in 2013"


Journal ArticleDOI
TL;DR: In this article, a two-dimensional structure of medium-scale traveling ionospheric disturbances (MSTIDs) over Europe has been revealed, for the first time, by using maps of the total electron content (TEC) obtained from more than 800 GPS receivers of the European GPS receiver networks.
Abstract: . Two-dimensional structures of medium-scale traveling ionospheric disturbances (MSTIDs) over Europe have been revealed, for the first time, by using maps of the total electron content (TEC) obtained from more than 800 GPS receivers of the European GPS receiver networks. From statistical analysis of the TEC maps obtained 2008, we have found that the observed MSTIDs can be categorized into two groups: daytime MSTID and nighttime MSTID. The daytime MSTID frequently occurs in winter. Its maximum occurrence rate in monthly and hourly bin exceeds 70% at lower latitudes over Europe, whereas it is approximately 45% at higher latitudes. Since most of the daytime MSTIDs propagate southward, we speculate that they could be caused by atmospheric gravity waves in the thermosphere. The nighttime MSTIDs also frequently occur in winter but most of them propagate southwestward, in a direction consistent with the theory that polarization electric fields play an important role in generating the nighttime MSTIDs. The nighttime MSTID occurrence rate shows distinct latitudinal difference: The maximum of the occurrence rate in monthly and hourly bin is approximately 50% at lower latitudes in Europe, whereas the nighttime MSTID was rarely observed at higher latitudes. We have performed model calculations of the plasma density perturbations caused by a gravity wave and an oscillating electric field to reproduce the daytime and nighttime MSTIDs, respectively. We find that TEC perturbations caused by gravity waves do not show dip angle dependencies, while those caused by the oscillating electric field have a larger amplitude at lower latitudes. These dip angle dependencies of the TEC perturbation amplitude could contribute to the latitudinal variation of the MSTID occurrence rate. Comparing with previous studies, we discuss the longitudinal difference of the nighttime MSTID occurrence rate, along with the E- and F-region coupling processes. The seasonal variation, of the nighttime MSTID occurrence rate in Europe, is not consistent with the theory that the longitudinal and seasonal variations of the nighttime MSTID occurrence could be attributed to those of the Es layer occurrence.

168 citations


Journal ArticleDOI
TL;DR: In this paper, high-speed jets are identified as intervals when the anti-sunward component of the dynamic pressure in the subsolar magnetosheath exceeds half of its upstream solar wind value.
Abstract: . Using 2008–2011 data from the five Time History of Events and Macroscale Interactions during Substorms (THEMIS) spacecraft in Earth's subsolar magnetosheath, we study high-speed jets identified as intervals when the anti-sunward component of the dynamic pressure in the subsolar magnetosheath exceeds half of its upstream solar wind value. Based on our comprehensive data set of 2859 high-speed jets, we obtain the following statistical results on jet properties and favorable conditions: high-speed jets occur predominantly downstream of the quasi-parallel bow shock, i.e., when interplanetary magnetic field cone angles are low. Apart from that, jet occurrence is only very weakly dependent (if at all) on other upstream conditions or solar wind variability. Typical durations and recurrence times of high-speed jets are on the order of tens of seconds and a few minutes, respectively. Relative to the ambient magnetosheath, high-speed jets exhibit higher speed, density and magnetic field intensity, but lower and more isotropic temperatures. They are almost always super-Alfvenic, often even super-magnetosonic, and typically feature 6.5 times as much dynamic pressure and twice as much total pressure in anti-sunward direction as the surrounding plasma does. Consequently, they are likely to have significant effects on the magnetosphere and ionosphere if they impinge on the magnetopause.

101 citations


Journal ArticleDOI
TL;DR: 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 as mentioned in this paper.
Abstract: . 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.

86 citations


Journal ArticleDOI
TL;DR: In this article, the main advances in the area of data-based modeling of the Earth's distant magnetic field achieved during the last two decades are reviewed. But the main goal of the approach is to extract maximum information from available data, using physically realistic and flexible mathematical structures, parameterized by the most relevant and routinely accessible observables.
Abstract: . This paper reviews the main advances in the area of data-based modelling of the Earth's distant magnetic field achieved during the last two decades. The essence and the principal goal of the approach is to extract maximum information from available data, using physically realistic and flexible mathematical structures, parameterized by the most relevant and routinely accessible observables. Accordingly, the paper concentrates on three aspects of the modelling: (i) mathematical methods to develop a computational "skeleton" of a model, (ii) spacecraft databases, and (iii) parameterization of the magnetospheric models by the solar wind drivers and/or ground-based indices. The review is followed by a discussion of the main issues concerning further progress in the area, in particular, methods to assess the models' performance and the accuracy of the field line mapping. The material presented in the paper is organized along the lines of the author Julius-Bartels' Medal Lecture during the General Assembly 2013 of the European Geosciences Union.

76 citations


Journal ArticleDOI
TL;DR: In this article, the authors analyzed 79 events between 1996 and 2009 reported in existing ICME/MC catalogs (Wind magnetic cloud list and the Richardson and Cane ICME list) using near-Earth observations by ACE (Advanced Composition Explorer) and Wind.
Abstract: . The relationship of magnetic clouds (MCs) to interplanetary coronal mass ejections (ICMEs) is still an open issue in space research. The view that all ICMEs would originate as magnetic flux ropes has received increasing attention, although near the orbit of the Earth only about one-third of ICMEs show clear MC signatures and often the MC occupies only a portion of the more extended region showing ICME signatures. In this work we analyze 79 events between 1996 and 2009 reported in existing ICME/MC catalogs (Wind magnetic cloud list and the Richardson and Cane ICME list) using near-Earth observations by ACE (Advanced Composition Explorer) and Wind. We perform a systematic comparison of cases where ICME and MC signatures coincided and where ICME signatures extended significantly beyond the MC boundaries. We find clear differences in the characteristics of these two event types. In particular, the events where ICME signatures continued more than 6 h past the MC rear boundary had 2.7 times larger speed difference between the ICME's leading edge and the preceding solar wind, 1.4 times higher magnetic fields, 2.1 times larger widths and they experienced three times more often strong expansion than the events for which the rear boundaries coincided. The events with significant mismatch in MC and ICME boundary times were also embedded in a faster solar wind and the majority of them were observed close to the solar maximum. Our analysis shows that the sheath, the MC and the regions of ICME-related plasma in front and behind the MC have different magnetic field, plasma and charge state characteristics, thus suggesting that these regions separate already close to the Sun. Our study shows that the geometrical effect (the encounter through the CME leg and/or far from the flux rope center) does not contribute much to the observed mismatch in the MC and ICME boundary times.

65 citations


Journal ArticleDOI
TL;DR: In this paper, 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.
Abstract: . 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.

64 citations


Journal ArticleDOI
TL;DR: In this paper, GPS-TEC data were observed at the same local time at two equatorial stations on both longitudes: Lagos (6.52° N, 3.4° E, 4.25° N magnetic latitude), Nigeria; and Pucallpa (8.38° S, 74.57° W,4.04° S magnetic latitude) during the minimum (2009, 2010) and ascending (2011) phases of solar cycle 24.
Abstract: . GPS-TEC data were observed at the same local time at two equatorial stations on both longitudes: Lagos (6.52° N, 3.4° E, 3.04° S magnetic latitude), Nigeria; and Pucallpa (8.38° S, 74.57° W, 4.25° N magnetic latitude), Peru during the minimum (2009, 2010) and ascending (2011) phases of solar cycle 24. These data were grouped into daily, seasonal and solar activity sets. The day-to-day variations in vertical TEC (VTEC) recorded the maximum during 14:00–16:00 LT and minimum during 04:00–06:00 LT at both longitudes. Seasonally, during solar minimum, maximum VTEC values were observed during March equinox and minimum during solstices. However, during the ascending phase of the solar activity, the maximum values were recorded during the December solstice and minimum during the June solstice. VTEC also increased with solar activity at both longitudes. On longitude by longitude comparison, the African GPS station generally recorded higher VTEC values than the American GPS station. Furthermore, harmonic analysis technique was used to extract the annual and semi-annual components of the amplitudes of the TEC series at both stations. The semi-annual variations dominated the TEC series over the African equatorial station, while the annual variations dominated those over the American equatorial station. The GPS-TEC-derived averages for non-storm days were compared with the corresponding values derived by the IRI-2007 with the NeQuick topside option. The NeQuick option of IRI-2007 showed better performance at the American sector than the African sector, but generally underestimating TEC during the early morning hours at both longitudes.

60 citations


Journal ArticleDOI
TL;DR: In this article, the role of disturbance dynamos caused by high-latitude winds (due to particle precipitation and Joule heating in the auroral zones) and prompt penetrating electric fields (PPEFs) in the solar wind-ionosphere coupling during these intervals are discussed.
Abstract: . 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 CO2 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 × 1010 W and 2.7 × 1010 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.

56 citations


Journal ArticleDOI
TL;DR: In this article, the authors used the Sodankyla ion and neutral chemistry model to calculate ionic production and loss rates of neutral H ox and NOy species, imposed by particle precipitation, for a range of atmospheric conditions and levels of ionization.
Abstract: . In the polar regions, precipitation of solar high-energy protons and electrons affects the neutral composition of the middle atmosphere. Here we use the Sodankyla Ion and Neutral Chemistry model to calculate ionic production and loss rates of neutral HOx and NOy species, imposed by particle precipitation, for a range of atmospheric conditions and levels of ionization. We also analyse in detail the ionic reaction sequences leading to the HOx and NOy changes. Our results show that particle impact ionization and positive ion chemistry cause net production of N, NO, HNO2, H, andOH from N2 and H2O. On the other hand, negative ion chemistry redistributes the NOy species, without net production or loss, so that NO, NO2, and N2O5 are converted to HNO3 and NO3. Based on the model results, we provide tables of so-called P/Q numbers (i.e. production and loss rates of neutral species divided by ionization rates) at altitudes between 20 and 90 km. These numbers can be easily used to parameterise the ion chemistry effects when modelling atmospheric response to particle precipitation. Compared to earlier studies, our work is the first to consider in detail the NOy effect of negative ion chemistry, and the diurnal and seasonal variability of the P/Q numbers.

55 citations


Journal ArticleDOI
TL;DR: In this paper, it is shown that direct radiative forcing by mineral dust reduces temperature in the lower atmosphere, but increases it in the layers aloft, and that surface cooling is associated with a reduction of turbulent kinetic energy (TKE) and hence vertical mixing of horizontal momentum.
Abstract: . 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.

54 citations


Journal ArticleDOI
TL;DR: Damboldt and Suessmann as mentioned in this paper derived long-term trends in the ionospheric F2 region using a twofold regression analysis depending on solar and geomagnetic activity.
Abstract: . A new comprehensive data collection by Damboldt and Suessmann (2012a) with monthly foF2 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 foF2 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 foF2 and hmF2 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 hmF2 values and to a smaller amount the foF2 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.

Journal ArticleDOI
TL;DR: In this article, the authors present occurrence rate statistics for nighttime medium-scale traveling ionospheric disturbances (MSTIDs) in the Central Pacific and South American sectors using data collected by 630.0 nm filtered CCD imaging systems.
Abstract: . We present occurrence rate statistics for nighttime medium-scale traveling ionospheric disturbances (MSTIDs) in the Central Pacific and South American sectors using data collected by 630.0 nm filtered CCD imaging systems. The data were collected from September 2006 through December 2012. In general, the statistics are in good agreement with the basic linear theory of MSTIDs, with observations coinciding with low F10.7A values, representative of solar minimum. Overall, MSTIDs are observed in approximately 68% of the usable nights near the solstices at mid-latitudes and approximately 20% of the usable nights for equinox periods. Observations closer to the geomagnetic equator yielded a maximum occurrence rate of about 10–20% during the solstices and about 0–3% during the equinoxes. The lower number of MSTID observations near the low latitudes is attributed to limitations of MSTID growth rate, propagation, and/or geometrical observational effects. The relatively large number of MSTID occurrences during the solstices can be accounted for by the neutral wind contribution to the MSTID growth rate either at the local or magnetic conjugate point.

Journal ArticleDOI
TL;DR: In this paper, the authors investigated the regional variability of raindrop size distribution (DSD) along the Equator through a network of Parsivel disdrometers in Indonesia and found that the DSD at PT has more large drops than at the other three sites.
Abstract: . Regional variability of raindrop size distribution (DSD) along the Equator was investigated through a network of Parsivel disdrometers in Indonesia. The disdrometers were installed at Kototabang (KT; 100.32° E, 0.20° S), Pontianak (PT; 109.37° E, 0.00° S), Manado (MN; 124.92° E, 1.55° N) and Biak (BK; 136.10° E, 1.18° S). It was found that the DSD at PT has more large drops than at the other three sites. The DSDs at the four sites are influenced by both oceanic and continental systems, and majority of the data matched the maritime-like DSD that was reported in a previous study. Continental-like DSDs were somewhat dominant at PT and KT. Regional variability of DSD is closely related to the variability of topography, mesoscale convective system propagation and horizontal scale of landmass. Different DSDs at different sites led to different Z–R relationships in which the radar reflectivity at PT was much larger than at other sites, at the same rainfall rate.

Journal ArticleDOI
TL;DR: In this paper, the authors present the observed equatorial electrojet spectrum and its relation to tidal driving, showing that a large fraction of the resulting spectrum is related to the switch on/off of the EEJ between day and night.
Abstract: . Based on 10 yr of magnetic field measurements by the CHAMP satellite we draw a detailed picture of the equatorial electrojet (EEJ) tidal variations. For the first time the complete EEJ spectrum related to average solar tides has been compiled. A large fraction of the resulting spectrum is related to the switch on/off of the EEJ between day and night. This effect has carefully been considered when interpreting the results. As expected, largest amplitudes are caused by the migrating tides representing the mean diurnal variation. Higher harmonics of the daily variations show a 1/f fall-off in amplitude. Such a spectrum is required to represent the vanishing of the EEJ current at night. The migrating tidal signal exhibits a distinct annual variation with large amplitudes during December solstice and equinox seasons but a depression by a factor of 1.7 around June–July. A rich spectrum of non-migrating tidal effects is deduced. Most prominent is the four-peaked longitudinal pattern around August. Almost 90% of the structure can be attributed to the diurnal eastward-propagating tide DE3. In addition the westward-propagating DW5 is contributing to wave-4. The second-largest non-migrating tide is the semi-diurnal SW4 around December solstice. It causes a wave-2 feature in satellite observations. The three-peaked longitudinal pattern, often quoted as typical for the December season, is significantly weaker. During the months around May–June a prominent wave-1 feature appears. To first order it represents a stationary planetary wave SPW1 which causes an intensification of the EEJ at western longitudes beyond 60° W and a weakening over Africa/India. In addition, a prominent ter-diurnal non-migrating tide TW4 causes the EEJ to peak later, at hours past 14:00 local time in the western sector. A particularly interesting non-migrating tide is the semi-diurnal SW3. It causes largest EEJ amplitudes from October through December. This tidal component shows a strong dependence on solar flux level with increasing amplitudes towards solar maximum. We are not aware of any previous studies mentioning this behaviour of SW3. The main focus of this study is to present the observed EEJ spectrum and its relation to tidal driving. For several of the identified spectral components we cannot offer convincing explanations for the generation mechanisms.

Journal ArticleDOI
TL;DR: In this article, a combined iterative and non-iterative reconstruction approach for the 3D water vapour inversion using GPS observations and COSMIC profiles is presented, which is evaluated using 10-day GPS observations in Hong Kong and radio occultation data.
Abstract: . Traditionally, balloon-based radiosonde soundings are used to study the spatial distribution of atmospheric water vapour. However, this approach cannot be frequently employed due to its high cost. In contrast, GPS tomography technique can obtain water vapour in a high temporal resolution. In the tomography technique, an iterative or non-iterative reconstruction algorithm is usually utilised to overcome rank deficiency of observation equations for water vapour inversion. However, the single iterative or non-iterative reconstruction algorithm has their limitations. For instance, the iterative reconstruction algorithm requires accurate initial values of water vapour while the non-iterative reconstruction algorithm needs proper constraint conditions. To overcome these drawbacks, we present a combined iterative and non-iterative reconstruction approach for the three-dimensional (3-D) water vapour inversion using GPS observations and COSMIC profiles. In this approach, the non-iterative reconstruction algorithm is first used to estimate water vapour density based on a priori water vapour information derived from COSMIC radio occultation data. The estimates are then employed as initial values in the iterative reconstruction algorithm. The largest advantage of this approach is that precise initial values of water vapour density that are essential in the iterative reconstruction algorithm can be obtained. This combined reconstruction algorithm (CRA) is evaluated using 10-day GPS observations in Hong Kong and COSMIC profiles. The test results indicate that the water vapor accuracy from CRA is 16 and 14% higher than that of iterative and non-iterative reconstruction approaches, respectively. In addition, the tomography results obtained from the CRA are further validated using radiosonde data. Results indicate that water vapour densities derived from the CRA agree with radiosonde results very well at altitudes above 2.5 km. The average RMS value of their differences above 2.5 km is 0.44 g m−3.

Journal ArticleDOI
TL;DR: In this paper, the amplitude and phase scintillation indices are customarily obtained by specialised GPS Ionospheric Scintillation and TEC Monitors (GISTMs) from L1 signal recorded at the rate of 50 Hz.
Abstract: . The amplitude and phase scintillation indices are customarily obtained by specialised GPS Ionospheric Scintillation and TEC Monitors (GISTMs) from L1 signal recorded at the rate of 50 Hz. The scintillation indices S4 and σΦ are stored in real time from an array of high-rate scintillation receivers of the Canadian High Arctic Ionospheric Network (CHAIN). Ionospheric phase scintillation was observed at high latitudes during a moderate geomagnetic storm (Dst = −61 nT) that was caused by a moderate solar wind plasma stream compounded with the impact of two coronal mass ejections. The most intense phase scintillation (σΦ ~ 1 rad) occurred in the cusp and the polar cap where it was co-located with a strong ionospheric convection, an extended tongue of ionisation and dense polar cap patches that were observed with ionosondes and HF radars. At sub-auroral latitudes, a sub-auroral polarisation stream that was observed by mid-latitude radars was associated with weak scintillation (defined arbitrarily as σΦ 0.1 rad and DPR > 2 mm s−1, both mapped as a function of magnetic latitude and magnetic local time, are very similar.

Journal ArticleDOI
TL;DR: Wang et al. as mentioned in this paper investigated the effect of total cloud cover (TCC) and sunshine duration (SSD) in the variation of diurnal temperature range (DTR) in China during 1954-2009.
Abstract: . This study aims to investigate the effect of total cloud cover (TCC) and sunshine duration (SSD) in the variation of diurnal temperature range (DTR) in China during 1954–2009. As expected, the inter-annual variation of DTR was mainly determined by TCC. Analysis of trends of 30-year moving windows of DTR and TCC time series showed that TCC changes could account for that of DTR in some cases. However, TCC decreased during 1954–2009, which did not support DTR reduction across China. DTRs under sky conditions such as clear, cloudy and overcast showed nearly the same decreasing rate that completely accounted for the overall DTR reduction. Nevertheless, correlation between SSD and DTR was weak and not significant under clear sky conditions in which aerosol direct radiative effect should be dominant. Furthermore, 30–60% of DTR reduction was associated with DTR decrease under overcast conditions in south China. This implies that aerosol direct radiative effect appears not to be one of the main factors determining long-term changes in DTR in China.

Journal ArticleDOI
TL;DR: In this article, the authors investigate the correlation between the geomagnetic field structure and the cutoff rigidity, and present a simple geometry-based parameter, δB, which is able to reflect the location-dependent cutoff Rigidity.
Abstract: . Studies of the propagation of charged energetic particles in the Earth's magnetic field go back to Carl Stormer. In the end, his investigations finally lead to the definition of the so-called cutoff rigidity RC; 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.

Journal ArticleDOI
TL;DR: In this paper, the authors investigated the distribution of the magnetic field and plasma variations in the compressed sheath regions of the interplanetary magnetic field (IMF) and in dynamic pressure (Pdyn) during solar cycle 23.
Abstract: . 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 (Pdyn) 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 Pdyn 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.

Journal ArticleDOI
TL;DR: In this paper, the authors used an automated recognition of substorm and storm phases in the auroral electrojet (AL) and ring current (Dst) index data from 1995-2009.
Abstract: . Substorm properties during different storm phases have been studied using an automated recognition of substorm and storm phases in the auroral electrojet (AL) and ring current (Dst) index data from 1995–2009. The large number of events (about 500 storms and 15 000 substorms) provides statistically reliable distributions, average behaviour and long time series of simple parameters, such as durations and intensities. The phases of storms and substorms have been examined independently. Substorm phases have been further combined to single and multi-cycle events. The former consist of one growth, one expansion and one recovery phase, while the latter include multiple expansion and recovery phases after one growth phase. Our findings show that most substorms take place during non-storm times, and substorms during storm initial phases resemble isolated non-storm time substorms. Both during storm initial phases and non-storm times, the substorm growth phases may last longer than the other substorm phases. Substorm recovery phase is typically the longest phase but its duration also varies most. The longest substorm recovery phase duration was observed during multi-cycle substorms. The longest substorm expansion and storm main phases were found during the years close to the solar maximum. The shortest substorm events (the shortest phase durations) are the single-cycle substorms. The period of expansion onsets during multi-cycle substorms varied hugely for events with a small number of expansion phases. For events with a larger number of expansions, a clearer periodicity of about one hour (median value) was suggested.

Journal ArticleDOI
TL;DR: In this paper, a global climatology of Pc1 pulsations as observed by the CHAMP satellite from 2000 to 2010 is presented, where the authors present a global analysis of the occurrence rate of pc1 events.
Abstract: . We present a global climatology of Pc1 pulsations as observed by the CHAMP satellite from 2000 to 2010. The Pc1 center frequency and bandwidth are about 1 and 0.5 Hz, respectively. The ellipticity is mostly linear with the major axis almost aligned with the magnetic zonal direction. The diurnal variation of Pc1 occurrences shows a primary maximum early in the morning and a secondary maximum during pre-midnight hours. The annual variations of the occurrence rates exhibit a clear preference for local summer. The solar cycle dependence of the occurrence rate reveals a maximum at the declining phase (2004–2005). Neither magnetic activity nor solar wind velocity controls the Pc1 occurrence rate significantly. Pc1 occurrence rate peaks at subauroral latitudes, but the steep cutoff towards higher latitudes is due to auroral field-aligned currents masking the Pc1 pulsations. The center frequency of Pc1 pulsations does not show a clear dependence on latitude. The global distribution of Pc1 exhibits highest occurrence rates near the longitude sector of the South Atlantic Anomaly. Pc1 events at auroral latitudes, although they are rarely detected, show a clear occurrence peak around local noon. A majority of the auroral Pc1 events are observed during solar minimum years.

Journal ArticleDOI
TL;DR: In this paper, the authors examined the Alfvenicity of a set of 188 solar wind directional discontinuities (DDs) identified in the Cluster data from 2003 by Knetter (2005), with the objective of separating rotational discontinuity (RDs) from tangential ones (TDs).
Abstract: . We examine the Alfvenicity of a set of 188 solar wind directional discontinuities (DDs) identified in the Cluster data from 2003 by Knetter (2005), with the objective of separating rotational discontinuities (RDs) from tangential ones (TDs). The DDs occurred over the full range of solar wind velocities and magnetic shear angles. By performing the Walen test in the de Hoffmann–Teller (HT) frame, we show that 77 of the 127 crossings for which a good HT frame was found had plasma flow speeds exceeding 80% of the Alfven speed at an average angular deviation of 7.7°; 33 cases had speeds exceeding 90% of the Alfven speed at an average angle of 6.4°. We show that the angular deviation between flow velocity (in the HT frame) and the Alfven velocity can be obtained from a reduced form of the Walen correlation coefficient. The corresponding results from the Walen test expressed in terms of jumps in flow speed and corresponding jumps in Alfven speed are similar: 66 of the same 127 cases had velocity jumps exceeding 80% with average angular deviation of 5.8°, and 22 exceeding 90% of the jump in Alfven speed, with average angular deviation 6.2°. We conclude that a substantial fraction of the 127 events can be identified as RDs. We present further evidence for coupling across the DDs by showing that, for most of the 127 crossings, the HT frame velocities, evaluated separately on the two sides of the DD, are nearly the same – a result required for RDs but not for TDs. We also show that the degree of Alfvenicity is nearly the same for the DDs and fluctuations in which the DDs are embedded. Whatever process causes deviations from ideal Alfvenicity appears to operate equally for the DDs as for the surrounding fluctuations. Finally, our study has established a unique relation between the strahl electron pitch angle and the sign of the Walen slope, implying antisunward propagation in the plasma frame for all 127 cases.

Journal ArticleDOI
TL;DR: In this paper, an experimental confirmation of the plasmasphere wind has been obtained using ion measurements acquired in the outer plasmosphere by the CIS experiment onboard the four Cluster spacecraft.
Abstract: . Plumes, forming at the plasmapause and released outwards, constitute a well-established mode for plasmaspheric material release to the Earth's magnetosphere. They are associated to active periods and the related electric field change. In 1992, Lemaire and Shunk proposed the existence of an additional mode for plasmaspheric material release to the Earth's magnetosphere: a plasmaspheric wind, steadily transporting cold plasmaspheric plasma outwards across the geomagnetic field lines, even during prolonged periods of quiet geomagnetic conditions. This has been proposed on a theoretical basis. Direct detection of this wind has, however, eluded observation in the past. Analysis of ion measurements, acquired in the outer plasmasphere by the CIS experiment onboard the four Cluster spacecraft, provide now an experimental confirmation of the plasmaspheric wind. This wind has been systematically detected in the outer plasmasphere during quiet and moderately active conditions, and calculations show that it could provide a substantial contribution to the magnetospheric plasma populations outside the Earth's plasmasphere. Similar winds should also exist on other planets, or astrophysical objects, quickly rotating and having an atmosphere and a magnetic field.

Journal ArticleDOI
TL;DR: In this paper, Lockwood et al. presented a new composite of geomagnetic activity which is designed to be as homogeneous in its construction as possible, by only combining data that, by virtue of the locations of the source observatories used, have similar responses to solar wind and IMF (interplanetary magnetic field) variations.
Abstract: . We present a new composite of geomagnetic activity which is designed to be as homogeneous in its construction as possible. This is done by only combining data that, by virtue of the locations of the source observatories used, have similar responses to solar wind and IMF (interplanetary magnetic field) variations. This will enable us (in Part 2, Lockwood et al., 2013a) to use the new index to reconstruct the interplanetary magnetic field, B, back to 1846 with a full analysis of errors. Allowance is made for the effects of secular change in the geomagnetic field. The composite uses interdiurnal variation data from Helsinki for 1845–1890 (inclusive) and 1893–1896 and from Eskdalemuir from 1911 to the present. The gaps are filled using data from the Potsdam (1891–1892 and 1897–1907) and the nearby Seddin observatories (1908–1910) and intercalibration achieved using the Potsdam–Seddin sequence. The new index is termed IDV(1d) because it employs many of the principles of the IDV index derived by Svalgaard and Cliver (2010), inspired by the u index of Bartels (1932); however, we revert to using one-day (1d) means, as employed by Bartels, because the use of near-midnight values in IDV introduces contamination by the substorm current wedge auroral electrojet, giving noise and a dependence on solar wind speed that varies with latitude. The composite is compared with independent, early data from European-sector stations, Greenwich, St Petersburg, Parc St Maur, and Ekaterinburg, as well as the composite u index, compiled from 2–6 stations by Bartels, and the IDV index of Svalgaard and Cliver. Agreement is found to be extremely good in all cases, except two. Firstly, the Greenwich data are shown to have gradually degraded in quality until new instrumentation was installed in 1915. Secondly, we infer that the Bartels u index is increasingly unreliable before about 1886 and overestimates the solar cycle amplitude between 1872 and 1883 and this is amplified in the proxy data used before 1872. This is therefore also true of the IDV index which makes direct use of the u index values.

Journal ArticleDOI
TL;DR: In this paper, simultaneous LAD1/2-observations from TWINS1 and TWINS2 between June 2008 and June 2010 were used to study seasonal variations in the H-geocorona.
Abstract: . The Lyman-α Detectors (LAD) on board the two TWINS 1/2-satellites allow for the simultaneous stereo imaging of the resonant emission glow of the H-geocorona from very different orbital positions. Terrestrial exospheric atomic hydrogen (H) resonantly scatters solar Lyman-α (121.567 nm) radiation. During the past solar minimum, relevant solar parameters that influence these emissions were quite stable. Here, we use simultaneous LAD1/2-observations from TWINS1 and TWINS2 between June 2008 and June 2010 to study seasonal variations in the H-geocorona. Data are combined to produce two datasets containing (summer) solstice and (combined spring and fall) equinox emissions. In the range from 3 to 10 Earth radii (RE), a three-dimensional (3-D) mathematical model is used that allows for density asymmetries in longitude and latitude. At lower geocentric distances (

Journal ArticleDOI
TL;DR: In this article, the authors presented a new reconstruction of the interplanetary magnetic field (IMF, B) for 1846-2012 with a full analysis of errors, based on the homogeneously constructed IDV(1d) composite of geomagnetic activity presented in Part 1 (Lockwood et al., 2013a).
Abstract: . We present a new reconstruction of the interplanetary magnetic field (IMF, B) for 1846–2012 with a full analysis of errors, based on the homogeneously constructed IDV(1d) composite of geomagnetic activity presented in Part 1 (Lockwood et al., 2013a). Analysis of the dependence of the commonly used geomagnetic indices on solar wind parameters is presented which helps explain why annual means of interdiurnal range data, such as the new composite, depend only on the IMF with only a very weak influence of the solar wind flow speed. The best results are obtained using a polynomial (rather than a linear) fit of the form B = χ · (IDV(1d) − β)α with best-fit coefficients χ = 3.469, β = 1.393 nT, and α = 0.420. The results are contrasted with the reconstruction of the IMF since 1835 by Svalgaard and Cliver (2010).

Journal ArticleDOI
Xiaoli Ren1, Honglin He1, Li Zhang1, Lei Zhou1, Guirui Yu1, Jiangtao Fan1 
TL;DR: In this paper, 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, was calculated using GIS and trend analysis techniques.
Abstract: . 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.

Journal ArticleDOI
TL;DR: In this paper, a modified composite day analysis for all-sky interferometric meteor radars when measuring mean winds and tides is described. But the results of this analysis are often contaminated by the tidal fields and suggest a method of empirical correction derived from a simple simulation model.
Abstract: . The advantages of using a composite day analysis for all-sky interferometric meteor radars when measuring mean winds and tides are widely known. On the other hand, problems arise if this technique is applied to Hocking's (2005) gravity wave analysis for all-sky meteor radars. In this paper we describe how a simple change in the procedure makes it possible to use a composite day in Hocking's analysis. Also, we explain how a modified composite day can be constructed to test its ability to measure gravity wave momentum fluxes. Test results for specified mean, tidal, and gravity wave fields, including tidal amplitudes and gravity wave momentum fluxes varying strongly with altitude and/or time, suggest that the modified composite day allows characterization of monthly mean profiles of the gravity wave momentum fluxes, with good accuracy at least at the altitudes where the meteor counts are large (from 89 to 92.5 km). In the present work we also show that the variances measured with Hocking's method are often contaminated by the tidal fields and suggest a method of empirical correction derived from a simple simulation model. The results presented here greatly increase our confidence because they show that our technique is able to remove the tide-induced false variances from Hocking's analysis.

Journal ArticleDOI
TL;DR: In this article, the authors study the momentum deposition in the thermosphere from the dissipation of small amplitude gravity waves (GWs) within a wave packet using a fully nonlinear two-dimensional compressible numerical model.
Abstract: We study the momentum deposition in the ther- mosphere from the dissipation of small amplitude gravity waves (GWs) within a wave packet using a fully nonlinear two-dimensional compressible numerical model. The model solves the nonlinear propagation and dissipation of a GW packet from the stratosphere into the thermosphere with re- alistic molecular viscosity and thermal diffusivity for various Prandtl numbers. The numerical simulations are performed for GW packets with initial vertical wavelengths ( z) rang- ing from 5 to 50 km. We show that z decreases in time as a GW packet dissipates in the thermosphere, in agreement with the ray trace results of Vadas and Fritts (2005) (VF05). We also find good agreement for the peak height of the mo- mentum flux (zdiss) between our simulations and VF05 for GWs with initial z 2H in an isothermal, windless back- ground, where H is the density scale height. We also confirm that zdiss increases with increasing Prandtl number. We in- clude eddy diffusion in the model, and find that the momen- tum deposition occurs at lower altitudes and has two sep- arate peaks for GW packets with small initial z. We also simulate GW packets in a non-isothermal atmosphere. The net z profile is a competition between its decrease from vis- cosity and its increase from the increasing background tem- perature. We find that the wave packet disperses more in the non-isothermal atmosphere, and causes changes to the mo- mentum flux and z spectra at both early and late times for GW packets with initial z 10 km. These effects are caused by the increase in T in the thermosphere, and the decrease in T near the mesopause.

Journal ArticleDOI
TL;DR: In this article, the dispersion properties of short-wavelength magnetic turbulence during a rare high-speed solar wind event with a flow velocity of about 700 km s−1 using magnetic field and ion data from the Cluster spacecraft were investigated.
Abstract: . Physical processes of the energy transport in solar wind turbulence are a subject of intense studies, and different ideas exist to explain them. This manuscript describes the investigation of dispersion properties in short-wavelength magnetic turbulence during a rare high-speed solar wind event with a flow velocity of about 700 km s−1 using magnetic field and ion data from the Cluster spacecraft. Using the multi-point resonator technique, the dispersion relations (i.e., frequency versus wave-number values in the solar wind frame) of turbulent magnetic fluctuations with wave numbers near the inverse ion inertial length are determined. Three major results are shown: (1) the wave vectors are uniformly quasi-perpendicular to the mean magnetic field; (2) the fluctuations show a broad range of frequencies at wavelengths around the ion inertial length; and (3) the direction of propagation at the observed wavelengths is predominantly in the sunward direction. These results suggest the existence of high-frequency dispersion relations partly associated with normal modes on small scales. Therefore nonlinear energy cascade processes seem to be acting that are not described by wave–wave interactions.