Search

ESA’s Swarm constellation mission consisting of three identical satellites, successfully launched on 22 November 2013, provides the excellent opportunity of reliable electric current density estimates in the ionosphere based on multi-satellite magnetic field measurements. The final constellation with Swarm A/C orbiting the Earth at about 470 km (flying side-by-side) and Swarm B at about 520 km altitude was achieved on 17 April 2014. Swarm satellites need a period of about 130 days to cover all local times; therefore the coverage of all local times and local seasons can be achieved in 5 years from the final constellation when the dual-satellite current density estimations have been started. We present the climatology of field-aligned currents (FACs) at polar cap, cusp and auroral regions, radial, interhemispheric (IHFACs) and F region dynamo currents at mid-latitude and equatorial regions, using the single and dual satellite approaches for the current density estimations. Also, the higher-flying Swarm B satellite provides from time to time well-matched observations with the 50 km lower Swarm A/C pair that gives the possibility for zonal current estimates by using Ampère’s ring integral for determining the mean current density passing through the encircled area.

The neutral density of the thermosphere is one of the key parameters used, in particular, to compute atmospheric drag acting on the Earth’s orbiting satellites and objects at an altitude below 1000 km. Various empirical and physical models of the thermospheric density have been derived since 1961. The most widely used empirical models are CIRA86, NRLMSISE00, JB2008, and DTM2013. A new empirical model (CHTherm-2018) has been recently developed at GFZ using 9-year accelerometer measurements from the CHAMP satellite from August 2000 to July 2009 when it flew at an altitude range from 460 to 310 km. A new method for estimating scale factors of integrated thermospheric density has been developed at DGFI-TUM using Satellite Laser Ranging (SLR) observations to spherical Earth orbiting satellites, in combination with a full precise orbit determination of these satellites. In this method, the total drag coefficient is computed using the Sentman model by averaging over all thermospheric species. Using this approach, we have estimated scale factors of the thermospheric density provided by CIRA86, NRLMSISE00, JB2008, DTM2013 and CH-Therm-2018 using SLR observations to the ANDE-P satellite from 16 August 2009 to 3 October 2009, the ANDE-C satellite from 16 August 2009 to 26 March 2010 and SpinSat from 29 December 2014 to 21 June 2015, i.e. at the periods of low and high solar activity, and the altitude range from 248 to 425 km. In this presentation, we describe the methods and models used and discuss the scaled factors obtained for various models of thermospheric density.

In this paper we estimate zonal plasma drift in the equatorial ionospheric F region without counting on ion drift meters. From June 2001 to June 2004 zonal plasma drift velocity is estimated from electron, neutral, and magnetic field observations of Challenging Mini-satellite Payload (CHAMP) in the 09:00–20:00 LT sector. The estimated velocities are validated against ion drift measurements by the Republic of China Satellite-1/Ionospheric Plasma and Electrodynamics Instrument (ROCSAT-1/IPEI) during the same period. The correlation between the CHAMP (altitude ~ 400 km) estimates and ROCSAT-1 (altitude ~ 600 km) observations is reasonably high (R ≈ 0.8). The slope of the linear regression is close to unity. However, the maximum westward drift and the westward-to-eastward reversal occur earlier for CHAMP estimates than for ROCSAT-1 measurements. In the equatorial F region both zonal wind and plasma drift have the same direction. Both generate vertical currents but with opposite signs. The wind effect (F region wind dynamo) is generally larger in magnitude than the plasma drift effect (Pedersen current generated by vertical E field), thus determining the direction of the F region vertical current.

Total electron content (TEC) between Low-Earth-Orbit (LEO) satellites and the Global Navigation Satellite System (GNSS) satellites can be used to constrain the three-dimensional morphology of equatorial plasma bubbles (EPBs). In this study we investigate TEC measured onboard the Challenging Minisatellite Payload (CHAMP) from 2001 to 2005. We only use TEC data obtained when CHAMP passed through EPBs: that is, when in situ plasma density measurements at CHAMP altitude also show EPB signatures. The observed TEC gradient along the CHAMP track is strongest when the corresponding GNSS satellite is located equatorward and westward of CHAMP with elevation angles of about 40–60°. These elevation and azimuth angles are in agreement with the angles expected from the morphology of the plasma depletion shell proposed by Kil et al.(2009).

The atmospheric lunar tide is one known source of ionospheric variability. The subject received renewed attention as recent studies found a link between stratospheric sudden warmings and amplified lunar tidal perturbations in the equatorial ionosphere. There is increasing evidence from ground observations that the lunar tidal influence on the ionosphere depends on longitude. We use magnetic field measurements from the CHAMP satellite during July 2000 to September 2010 and from the two Swarm satellites during November 2013 to February 2017 to determine, for the first time, the complete seasonal- longitudinal climatology of the semidiurnal lunar tidal variation in the equatorial electrojet intensity. Significant longitudinal variability is found in the amplitude of the lunar tidal variation, while the longitudinal variability in the phase is small. The amplitude peaks in the Peruvian sector (similar to 285 degrees E) during the Northern Hemisphere winter and equinoxes, and in the Brazilian sector (similar to 325 degrees E) during the Northern Hemisphere summer. There are also local amplitude maxima at similar to 55 degrees E and similar to 120 degrees E. The longitudinal variation is partly due to the modulation of ionospheric conductivities by the inhomogeneous geomagnetic field. Another possible cause of the longitudinal variability is neutral wind forcing by nonmigrating lunar tides. A tidal spectrum analysis of the semidiurnal lunar tidal variation in the equatorial electrojet reveals the dominance of the westward propagating mode with zonal wave number 2 (SW2), with secondary contributions by westward propagating modes with zonal wave numbers 3 (SW3) and 4 (SW4). Eastward propagating waves are largely absent from the tidal spectrum. Further study will be required for the relative importance of ionospheric conductivities and nonmigrating lunar tides.

We present in a statistical study a comparison of thermospheric mass density enhancements (ρrel) with electron temperature (Te), small-scale field-aligned currents (SSFACs), and vertical ion velocity (Vz) at high latitudes around noon magnetic local time (MLT). Satellite data from CHAMP (CHAllenging Minisatellite Payload) and DMSP (Defense Meteorological Satellite Program) sampling the Northern Hemisphere during the years 2002–2005 are used. In a first step we investigate the distribution of the measured quantities in a magnetic latitude (MLat) versus MLT frame. All considered variables exhibit prominent peak amplitudes in the cusp region. A superposed epoch analysis was performed to examine causal relationship between the quantities. The occurrence of a thermospheric relative mass density anomaly, ρrel >1.2, in the cusp region is defining an event. The location of the density peak is taken as a reference latitude (Δ MLat = 0°). Interestingly, all the considered quantities, SSFACs, Te, and Vz are co-located with the density anomaly. The amplitudes of the peaks exhibit different characters of seasonal variation. The average relative density enhancement of the more prominent density peaks considered in this study amounts to 1.33 during all seasons. As expected, SSFACs are largest in summer with average amplitudes equal to 2.56 μA m−2, decaying to 2.00 μA m−2 in winter. The event related enhancements of Te and Vz are both largest in winter (Δ Te =730 K, Vz =136 m s−1) and smallest in summer (Δ Te = 377 K, Vz = 57 m s−1. Based on the similarity of the seasonal behaviour we suggest a close relationship between these two quantities. A correlation analysis supports a linear relation with a high coefficient greater than or equal to 0.93, irrespective of season. Our preferred explanation is that dayside reconnection fuels Joule heating of the thermosphere causing air upwelling and at the same time heating of the electron gas that pulls up ions along affected flux tubes.

By using the accelerometer measurements from CHAMP and GRACE satellites, the tidal signatures of the thermospheric mass density and zonal wind at midlatitudes have been analyzed in this study. The results show that the mass density and zonal wind at southern midlatitudes are dominated by a longitudinal wave-1 pattern. The most prominent tidal components in mass density and zonal wind are the diurnal tides D0 and DW2 and the semidiurnal tides SW1 and SW3. This is consistent with the tidal signatures in the F region electron density at midlatitudes as reported by Xiong and Lühr (2014). These same tidal components are observed both in the thermospheric and ionospheric quantities, supporting a mechanism that the non-migrating tides in the upper atmosphere are excited in situ by ion–neutral interactions at midlatitudes, consistent with the modeling results of Jones Jr. et al. (2013). We regard the thermospheric dynamics as the main driver for the electron density tidal structures. An example is the in-phase variation of D0 between electron density and mass density in both hemispheres. Further research including coupled atmospheric models is probably needed for explaining the similarities and differences between thermospheric and ionospheric tidal signals at midlatitudes.