skip to main content


Title: Storm time neutral density assimilation in the thermosphere ionosphere with TIDA
To improve Thermosphere–Ionosphere modeling during disturbed conditions, data assimilation schemes that can account for the large and fast-moving gradients moving through the modeled domain are necessary. We argue that this requires a physics based background model with a non-stationary covariance. An added benefit of using physics-based models would be improved forecasting capability over largely persistence-based forecasts of empirical models. As a reference implementation, we have developed an ensemble Kalman Filter (enKF) software called Thermosphere Ionosphere Data Assimilation (TIDA) using the physics-based Coupled Thermosphere Ionosphere Plasmasphere electrodynamics (CTIPe) model as the background. In this paper, we present detailed results from experiments during the 2003 Halloween Storm, 27–31 October 2003, under very disturbed ( K p  = 9) conditions while assimilating GRACE-A and B, and CHAMP neutral density measurements. TIDA simulates this disturbed period without using the L1 solar wind measurements, which were contaminated by solar energetic protons, by estimating the model drivers from the density measurements. We also briefly present statistical results for two additional storms: September 27 – October 2, 2002, and July 26 – 30, 2004, to show that the improvement in assimilated neutral density specification is not an artifact of the corrupted forcing observations during the 2003 Halloween Storm. By showing statistical results from assimilating one satellite at a time, we show that TIDA produces a coherent global specification for neutral density throughout the storm – a critical capability in calculating satellite drag and debris collision avoidance for space traffic management.  more » « less
Award ID(s):
1651459
NSF-PAR ID:
10335694
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Space Weather and Space Climate
Volume:
12
ISSN:
2115-7251
Page Range / eLocation ID:
16
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Inter-hemispheric asymmetry (IHA) in Earth’s ionosphere–thermosphere (IT) system can be associated with high-latitude forcing that intensifies during storm time, e.g., ion convection, auroral electron precipitation, and energy deposition, but a comprehensive understanding of the pathways that generate IHA in the IT is lacking. Numerical simulations can help address this issue, but accurate specification of high-latitude forcing is needed. In this study, we utilize the Active Magnetosphere and Planetary Electrodynamics Response Experiment-revised fieldaligned currents (FACs) to specify the high-latitude electric potential in the Global Ionosphere and Thermosphere Model (GITM) during the October 8–9, 2012, storm. Our result illustrates the advantages of the FAC-driven technique in capturing high-latitude ion drift, ion convection equatorial boundary, and the storm-time neutral density response observed by satellite. First, it is found that the cross-polar-cap potential, hemispheric power, and ion convection distribution can be highly asymmetric between two hemispheres with a clear Bydependence in the convection equatorial boundary. Comparison with simulation based on mirror precipitation suggests that the convection distribution is more sensitive to FAC, while its intensity also depends on the ionospheric conductance-related precipitation. Second, the IHA in the neutral density response closely follows the IHA in the total Joule heating dissipation with a time delay. Stronger Joule heating deposited associated with greater high-latitude electric potential in the southern hemisphere during the focus period generates more neutral density as well, which provides some evidences that the high-latitude forcing could become the dominant factor to IHAs in the thermosphere when near the equinox. Our study improves the understanding of storm-time IHA in high-latitude forcing and the IT system.

     
    more » « less
  2. Abstract

    Ionospheric storm enhanced density (SED) has been extensively investigated using total electron content deduced from GPS ground and satellite‐borne receivers. However, dayside in situ electron density measurements have not been analyzed in detail for SEDs yet. We report in situ electron density measurements of a SED event in the Northern Hemisphere (NH) at the noon meridian plane measured by the Challenging Minisatellite Payload (CHAMP) polar‐orbiting satellite at about 390 km altitude during the 20 November 2003 magnetic storm. The CHAMP satellite measurements render rare documentation about the dayside SED's life cycle at a fixed magnetic local time (MLT) through multiple passes. Solar wind drivers triggered the SED onset and controlled its lifecycle through its growth and retreat phases. The SED electron density enhancement extended from the equatorial ionization anomaly to the noon cusp. The midlatitude electron density increased to a maximum at the end of the growth phase. Afterward, the dayside SED region retreated gradually to lower magnetic latitudes. The observations showed a hemisphere asymmetry, with the NH electron density exhibiting a more significant enhancement. The simulations using the Thermosphere Ionosphere Electrodynamic General Circulation model show a good agreement with the CHAMP observations. The simulations indicate that the dayside midlatitude electron density enhancement has a complicated dependence on vertical ion drift, neutral wind, magnetic latitude, MLT, and the height of the F2 layer. Finally, we discuss the notion of using the mean cross‐polar cap electric field as a proxy for assessing the effects of solar wind drivers on producing midlatitude electron density enhancement.

     
    more » « less
  3. Abstract

    The dayside equatorial ionospheric electrodynamics exhibit strong variability driven simultaneously by highly changeable external forcings that originate from the solar extreme ultraviolet (EUV), magnetosphere, and lower atmosphere. We investigate this variability by carrying out comprehensive data‐driven ensemble modeling using a coupled model of the thermosphere and ionosphere, with the focus on the verticalE×Bdrift variability during a solar minimum and minor storm period. The variability of verticalE×Bdrift in response to the changes and uncertainty of primary forcings (i.e., solar EUV, high‐latitude plasma convection and auroral particle precipitation, and lower‐atmospheric tide and wave forcing) is investigated by ensemble forcing sensitivity experiments that incorporate data‐driven stochastic perturbations of these forcings into the model. Second, the impact of assimilating FORMOsa SATellite‐3/Constellation Observing System for Meteorology, Ionosphere, and Climate (FORMOSAT‐3/COSMIC) electron density profiles (EDPs) on the reduction of uncertainty of the modeled verticalE×Bdrift variability resulting from inadequately specified external forcing is revealed. The Communication and Navigation Outage Forecasting System (C/NOFS) ion drift velocity observations are used for validation. The validation results support the importance of the use of a data‐driven forcing perturbation methods in ensemble modeling and data assimilation. In conclusion, the solar EUV dominates the global‐scale day‐to‐day variability, while the lower atmosphere tide and wave forcing is critical to determining the regional variability. The modeled verticalE×Bdrift is also sensitive to the magnetospheric forcing. The ensemble data assimilation of FORMOSAT‐3/COSMIC EDPs helps to reduce the uncertainty and improves agreement of the modeled verticalE×Bdrifts with C/NOFS observations.

     
    more » « less
  4. Abstract

    Thermospheric mass density perturbations are commonly observed during geomagnetic storms and fundamental to upper atmosphere dynamics, but the sources of these perturbations are not well understood. Large neutral density perturbations during storms greatly affect the drag experienced by low Earth orbit. We investigated the thermospheric density perturbations at all latitudes observed along the CHAMP and GRACE satellite trajectories during the August 24–25, 2005 geomagnetic storm. Observations show that large neutral density enhancements occurred not only at high latitudes, but also globally. Large density perturbations were seen in the equatorial regions away from the high‐latitude, magnetospheric energy sources. We used the high‐resolution Multiscale Atmosphere Geospace Environment (MAGE) model to simulate consecutive neutral density changes observed by satellites during the storm. The MAGE simulation, which resolved mesoscale high‐latitude convection electric fields and field‐aligned currents, and included physics‐based specification of auroral precipitation, was contrasted with a standalone ionosphere‐thermosphere simulation driven by a high‐latitude electrodynamics empirical model. The comparison demonstrates that first‐principles representations of highly dynamic and localized Joule heating events in a fully coupled whole geospace model is critical to accurately capture both generation and propagation of traveling atmospheric disturbances (TADs) that produce neutral density perturbations globally. The MAGE simulation shows that larger density peaks in the equatorial region observed by CHAMP and GRACE are the result of TADs generated at high‐latitudes in both hemispheres, and intersect at low‐latitudes. This study reveals the importance of investigating thermospheric density variations at all latitudes in a fully coupled geospace model with sufficiently high resolving power.

     
    more » « less
  5. Abstract

    The upper boundary height of the traditional community general circulation model of the ionosphere‐thermosphere system is too low to be applied to the topside ionosphere/thermosphere study. In this study, the National Center for Atmospheric Research Thermosphere‐Ionosphere‐Electrodynamics General Circulation Model (NCAR‐TIEGCM) was successfully extended upward by four scale heights from 400–600 km to 700–1,200 km depending on solar activity, named TIEGCM‐X. The topside ionosphere and thermosphere simulated by TIEGCM‐X agree well with the observations derived from a topside sounder and satellite drag data. In addition, the neutral density, temperature, and electron density simulated by TIEGCM‐X are morphologically consistent with the NCAR‐TIEGCM simulations before extension. The latitude‐altitude distribution of the equatorial ionization anomaly derived from TIEGCM‐X is more reasonable. During geomagnetic storm events, the thermospheric responses of TIEGCM‐X are similar to NCAR‐TIEGCM. However, the ionospheric storm effects in TIEGCM‐X are stronger than those in NCAR‐TIEGCM and are even opposites at some middle and low latitudes due to the presence of more closed magnetic field lines. Defense Meteorological Satellite Program observations prove that the ionospheric storm effect of TIEGCM‐X is more reasonable. The well‐validated TIEGCM‐X has significant potential applications in ionospheric/thermospheric studies, such as the responses to storms, low‐latitude dynamics, and data assimilation.

     
    more » « less