skip to main content


Title: Importance of lower atmospheric forcing and magnetosphere-ionosphere coupling in simulating neutral density during the February 2016 geomagnetic storm
During geomagnetic storms a large amount of energy is transferred into the ionosphere-thermosphere (IT) system, leading to local and global changes in e.g., the dynamics, composition, and neutral density. The more steady energy from the lower atmosphere into the IT system is in general much smaller than the energy input from the magnetosphere, especially during geomagnetic storms, and therefore details of the lower atmosphere forcing are often neglected in storm time simulations. In this study we compare the neutral density observed by Swarm-C during the moderate geomagnetic storm of 31 January to 3 February 2016 with the Thermosphere-Ionosphere-Electrodynamics-GCM (TIEGCM) finding that the model can capture the observed large scale neutral density variations better in the southern than northern hemisphere. The importance of more realistic lower atmospheric (LB) variations as specified by the Whole Atmosphere Community Climate Model eXtended (WACCM-X) with specified dynamics (SD) is demonstrated by improving especially the northern hemisphere neutral density by up to 15% compared to using climatological LB forcing. Further analysis highlights the importance of the background atmospheric condition in facilitating hemispheric different neutral density changes in response to the LB perturbations. In comparison, employing observationally based field-aligned current (FAC) versus using an empirical model to describe magnetosphere-ionosphere (MI) coupling leads to an 7–20% improved northern hemisphere neutral density. The results highlight the importance of the lower atmospheric variations and high latitude forcing in simulating the absolute large scale neutral density especially the hemispheric differences. However, focusing on the storm time variation with respect to the quiescent time, the lower atmospheric influence is reduced to 1–1.5% improvement with respect to the total observed neutral density. The results provide some guidance on the importance of more realistic upper boundary forcing and lower atmospheric variations when modeling large scale, absolute and relative neutral density variations.  more » « less
Award ID(s):
2002574
NSF-PAR ID:
10423064
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Frontiers in Astronomy and Space Sciences
Volume:
9
ISSN:
2296-987X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  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. 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

    The ionospheric density displays hemispheric asymmetries in the polar region due to various hemispheric differences, for example, in the offset between geographic and geomagnetic poles and in the geomagnetic field strength. Using ground‐based ionospheric measurements from Vertical Incidence Pulsed Ionospheric Radar with Dynasonde analysis at Jang Bogo Station (JBS), Antarctica and from EISCAT Svalbard Radar (ESR) where both sites are located mostly in the polar cap, we investigate the hemispheric differences in the ionospheric density between the northern and southern hemispheres for geomagnetically quiet and solar minimum condition. The results are also compared with Thermosphere Ionosphere Electrodynamic Global Circulation Model (TIEGCM) simulations. The observations show larger density and stronger diurnal and seasonal variations at JBS in the southern hemisphere than at Svalbard in the northern hemisphere. The diurnal variations of the density peak height are also observed to be much larger at JBS. In both hemispheres, the ionospheric density is significantly reduced in winter due to the limited solar production at high geographic latitudes, but TIEGCM considerably overestimates winter density, which is even larger than summer density, especially in the northern hemisphere. Also existed are the differences in the equinoctial asymmetry between the observations and the simulations: the daytime F‐region density is observed to be larger in fall than in spring in both hemispheres, but TIEGCM shows the opposite. In general, most of the observed asymmetrical density are much weaker in the model simulation, which may result from lack of proper magnetospheric forcings and neutral dynamics in the model.

     
    more » « less
  5. Abstract

    The National Aeronautics and Space Administration Global‐scale Observations of the Limb and Disk (GOLD) mission observed a unique structure of thermospheric column density ratio of O to N2(∑O/N2) during a geomagnetic storm on day of year (DOY) 130 (May 10) to DOY 132 in 2019. The percentage difference of ∑O/N2between the storm time (DOY 131) and the quiet time (DOY 128) had a relatively enhanced ∑O/N2region sandwiched by two depleted regions over North America and the Atlantic Ocean in the Northern Hemisphere. This enhanced ∑O/N2region is called the neutral tongue here. The National Center for Atmospheric Research Thermosphere Ionosphere Electrodynamics General Circulation Model also predicted a similar ∑O/N2structure with the same spatial and temporal evolution as that seen in the GOLD observations on DOY 131. Model diagnostic analysis revealed that the neutral tongue was formed when the neutral wind in the region changed from mostly equatorward to poleward, which transported the ∑O/N2rich air from lower latitudes into this longitude section, forming the neutral tongue. The neutral tongue separated a large depletion region formed earlier into two smaller depleted regions, one on each side of the neutral tongue.

     
    more » « less