skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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
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. 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
  2. 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
  3. Low Earth orbit (LEO) radio occultation|radio occultations (RO) constellations can provide global electron density profiles (EDPs) to better specify and forecast the ionosphere‐thermosphere (I‐T) system. To inform future RO constellation design, this study uses comprehensive Observing System Simulation Experiments (OSSEs) to assess the ionospheric specification impact of assimilating synthetic EDPs into a coupled I‐T model. These OSSEs use 10 different sets of RO constellation configurations containing 6 or 12 LEO satellites with base orbit parameter combinations of 520 or 800 km altitude, and 24° or 72° inclination. The OSSEs are performed using the Ensemble Adjustment Kalman Filter implemented in the data assimilation (DA) Research Testbed and the Thermosphere‐Ionosphere‐Electrodynamics General Circulation Model (TIEGCM). A different I‐T model is used for the nature run, the Whole Atmosphere Model‐Ionosphere Plasmasphere Electrodynamics (WAM‐IPE), to simulate the period of interest is the St. Patrick's Day storm on March 13–18, 2015. Errors from models and EDP retrieval are realistically accounted for in this study through distinct I‐T models and by retrieving synthetic EDPs through an extension Abel inversion algorithm. OSSE assessment, using multiple metrics, finds that greater EDP spatial coverage leading to improved specification at altitudes 300 km and above, with the 520 km altitude constellations performing best due to yielding the highest observation counts. A potential performance limit is suggested with two 6‐satellite constellations. Lastly, close examination of Abel inversion error impacts highlights major EDP limitations at altitudes below 200 km and dayside equatorial regions with large horizontal gradients and low electron density magnitudes. 
    more » « less
  4. Abstract The geospace environment is volatile and highly driven. Space weather has effects on Earth's magnetosphere that cause a dynamic and enigmatic response in the thermosphere, particularly on the evolution of neutral mass density. Many models exist that use space weather drivers to produce a density response, but these models are typically computationally expensive or inaccurate for certain space weather conditions. In response, this work aims to employ a probabilistic machine learning (ML) method to create an efficient surrogate for the Thermosphere Ionosphere Electrodynamics General Circulation Model (TIE‐GCM), a physics‐based thermosphere model. Our method leverages principal component analysis to reduce the dimensionality of TIE‐GCM and recurrent neural networks to model the dynamic behavior of the thermosphere much quicker than the numerical model. The newly developed reduced order probabilistic emulator (ROPE) uses Long‐Short Term Memory neural networks to perform time‐series forecasting in the reduced state and provide distributions for future density. We show that across the available data, TIE‐GCM ROPE has similar error to previous linear approaches while improving storm‐time modeling. We also conduct a satellite propagation study for the significant November 2003 storm which shows that TIE‐GCM ROPE can capture the position resulting from TIE‐GCM density with <5 km bias. Simultaneously, linear approaches provide point estimates that can result in biases of 7–18 km. 
    more » « less
  5. Abstract A new TEC‐based ionospheric data assimilation system (TIDAS) over the continental US and adjacent area (20°–60°N, 60°–130°W, and 100–600 km) has been developed through assimilating heterogeneous ionospheric data, including dense ground‐based Global Navigation Satellite System (GNSS) Total Electron Content (TEC) from 2,000+ receivers, Constellation Observing System for Meteorology, Ionosphere, and Climate radio occultation data, JASON satellite altimeter TEC, and Millstone Hill incoherent scatter radar measurements. A hybrid Ensemble‐Variational scheme is utilized to reconstruct the regional 3‐D electron density distribution: a more realistic and location‐dependent background error covariance matrix is calculated from an ensemble of corrected NeQuick outputs, and a three‐dimensional variational (3DVAR) method is adopted for measurement updates to obtain an optimal state estimation. The spatial‐temporal resolution of the reanalyzed 3‐D electron density product is as high as 1° × 1° in latitude and longitude, 20 km in altitude, and 5 min in universal time, which is sufficient to reproduce ionospheric fine structure and storm‐time disturbances. The accuracy and reliability of data assimilation results are validated using ionosonde and other measurements. TIDAS reanalyzed electron density is able to successfully reconstruct the 3‐D morphology and dynamic evolution of the storm‐enhanced density (SED) plume observed during the St. Patrick's day geomagnetic storm on 17 March 2013 with high fidelity. Using TIDAS, we found that the 3‐D SED plume manifests as a ridge‐like high‐density channel that predominantly occurred between 300 and 500 km during 19:00–21:00 UT for this event, with the F2 region peak height being raised by 40–60 km and peak density enhancement of 30%–50%. 
    more » « less