skip to main content


Title: Explaining Solar Flare‐Induced Ionospheric Ion Upflow at Millstone Hill (42.6°N)
Abstract

Previous studies have shown that solar flares can significantly affect Earth's ionosphere and induce ion upflow with a magnitude of ∼110 m/s in the topside ionosphere (∼570 km) at Millstone Hill (42.61°N, 71.48°W). We use simulations from the Thermosphere‐Ionosphere‐Electrodynamics General Circulation Model (TIEGCM) and observations from Incoherent Scatter Radar (ISR) at Millstone Hill to reveal the mechanism of ionospheric ion upflow near the X9.3 flare peak (07:16 LT) on 6 September 2017. The ISR observed ionospheric upflow was captured by the TIEGCM in both magnitude and morphology. The term analysis of the F‐region ion continuity equation during the solar flare shows that the ambipolar diffusion enhancement is the main driver for the upflow in the topside ionosphere, while ion drifts caused by electric fields and neutral winds play a secondary role. Further decomposition of the ambipolar diffusive velocity illustrates that flare‐induced changes in the vertical plasma density gradient is responsible for ion upflow. The changes in the vertical plasma density gradient are mainly due to solar extreme ultraviolet (EUV, 15.5–79.8 nm) induced electron density and temperature enhancements at the F2‐region ionosphere with a minor and indirectly contribution from X‐ray (0–15.5 nm) and ultraviolet (UV, 79.8–102.7 nm).

 
more » « less
Award ID(s):
2033787 1952737
NSF-PAR ID:
10444595
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Space Physics
Volume:
127
Issue:
7
ISSN:
2169-9380
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The high‐resolution thermosphere‐ionosphere‐electrodynamics general circulation model has been used to investigate the response ofF2region electron density (Ne) at Millstone Hill (42.61°N, 71.48°W, maximum obscuration: 63%) to the Great American Solar Eclipse on 21 August 2017. Diagnostic analysis of model results shows that eclipse‐induced disturbance winds causeF2region Ne changes directly by transporting plasma along field lines, indirectly by producing enhanced O/N2ratio that contribute to the recovery of the ionosphere at and below theF2peak after the maximum obscuration. Ambipolar diffusion reacts to plasma pressure gradient changes and modifies Ne profiles. Wind transport and ambipolar diffusion take effect from the early phase of the eclipse and show strong temporal and altitude variations. The recovery ofF2region electron density above theF2peak is dominated by the wind transport and ambipolar diffusion; both move the plasma to higher altitudes from below theF2peak when more ions are produced in the lowerF2region after the eclipse. As the moon shadow enters, maximizes, and leaves a particular observation site, the disturbance winds at the site change direction and their effects on theF2region electron densities also vary, from pushing plasma downward during the eclipse to transporting it upward into the topside ionosphere after the eclipse. Chemical processes involving dimming solar radiation and changing composition, wind transport, and ambipolar diffusion together cause the time delay and asymmetric characteristic (fast decrease of Ne and slow recovery of the eclipse effects) of the topside ionospheric response seen in Millstone Hill incoherent scatter radar observations.

     
    more » « less
  2. This paper studies the three-dimensional (3-D) ionospheric electron density variation over the continental US and adjacent regions during the August 2017 Great American Solar Eclipse event, using Millstone Hill incoherent scatter radar observations, ionosonde data, the Swarm satellite measurements, and a new TEC-based ionospheric data assimilation system (TIDAS). The TIDAS data assimilation system can reconstruct a 3-D electron density distribution over continental US and adjacent regions, with a spatial–temporal resolution of 1∘× 1∘ in latitude and longitude, 20 km in altitude, and 5 min in universal time. The combination of multi-instrumental observations and the high-resolution TIDAS data assimilation products can well represent the dynamic 3-D ionospheric electron density response to the solar eclipse, providing important altitude information and fine-scale details. Results show that the eclipse-induced ionospheric electron density depletion can exceed 50% around the F2-layer peak height between 200 and 300 km. The recovery of electron density following the maximum depletion exhibits an altitude-dependent feature, with lower altitudes exhibiting a faster recovery than the F2 peak region and above. The recovery feature was also characterized by a post-eclipse electron density enhancement of 15–30%, which is particularly prominent in the topside ionosphere at altitudes above 300 km.

     
    more » « less
  3. Abstract

    Sudden enhancement in high‐frequency absorption is a well‐known impact of solar flare‐driven Short‐Wave Fadeout (SWF). Less understood, is a perturbation of the radio wave frequency as it traverses the ionosphere in the early stages of SWF, also known as the Doppler flash. Investigations have suggested two possible sources that might contribute to it’s manifestation: first, enhancements of plasma density in the D‐and lower E‐regions; second, the lowering of the F‐region reflection point. Our recent work investigated a solar flare event using first principles modeling and Super Dual Auroral Radar Network (SuperDARN) HF radar observations and found that change in the F‐region refractive index is the primary driver of the Doppler flash. This study analyzes multiple solar flare events observed across different SuperDARN HF radars to determine how flare characteristics, properties of the traveling radio wave, and geophysical conditions impact the Doppler flash. In addition, we use incoherent scatter radar data and first‐principles modeling to investigate physical mechanisms that drive the lowering of the F‐region reflection points. We found, (a) on average, the change in E‐ and F‐region refractive index is the primary driver of the Doppler flash, (b) solar zenith angle, ray’s elevation angle, operating frequency, and location of the solar flare on the solar disk can alter the ionospheric regions of maximum contribution to the Doppler flash, (c) increased ionospheric Hall and Pedersen conductance causes a reduction of the daytime eastward electric field, and consequently reduces the vertical ion‐drift in the lower and middle latitude ionosphere, which results in lowering of the F‐region ray reflection point.

     
    more » « less
  4. Abstract

    E‐region models have traditionally underestimated the ionospheric electron density. We believe that this deficiency can be remedied by using high‐resolution photoabsorption and photoionization cross sections in the models. Deep dips in the cross sections allow solar radiation to penetrate deeper into the E‐region producing additional ionization. To validate our concept, we perform a study of model electron density profiles (EDPs) calculated using the Atmospheric Ultraviolet Radiance Integrated Code (AURIC; D. Strickland et al., 1999,https://doi.org/10.1016/s0022-4073(98)00098-3) in the E‐region of the terrestrial ionosphere. We compare AURIC model outputs using new high‐resolution photoionization and photoabsorption cross sections, and solar spectral irradiances during low solar activity with incoherent scatter radar (ISR) measurements from the Arecibo and Millstone Hills observatories, Constellation Observing System for Meteorology Ionosphere and Climate (COSMIC‐1) observations, and outputs from empirical models (IRI‐2016 and FIRI‐2018). AURIC results utilizing the new high‐resolution cross sections reveal a significant difference to model outputs calculated with the low‐resolution cross sections currently used. Analysis of AURIC EDPs using the new high‐resolution data indicate fair agreement with ISR measurements obtained at various times at Arecibo but very good agreement with Millstone Hills ISR observations from ∼96–140 km. However, discrepancies in the altitude of the E‐region peak persist. High‐resolution AURIC calculations are in agreement with COSMIC‐1 observations and IRI‐2016 model outputs between ∼105 and 140 km while FIRI‐2018 outputs underestimate the EDP in this region. Overall, AURIC modeling shows increased E‐region electron densities when utilizing high‐resolution cross sections and high‐resolution solar irradiances, and are likely to be the key to resolving the long standing data‐model discrepancies.

     
    more » « less
  5. Abstract

    Ion upflow in theFregion and topside ionosphere can greatly influence the ion density and fluxes at higher altitudes and thus has significant impact on ion outflow. We investigated the statistical characteristics of ion upflow and downflow using a 3‐year (2011–2013) data set from the Poker Flat Incoherent Scatter Radar (PFISR). Ion upflow is twice more likely to occur on the nightside than on the dayside in PFISR observations, while downflow events occur more often in the afternoon sector. Upflow and downflow on the dayside tend to occur at altitudes ~500 km, higher than those on the nightside. Both upflow and downflow occur more frequently as ion convection speed increases. Upflow observed from 16 to 6 magnetic local time through midnight is associated with temperature and density enhancements. Occurrence rates of upflow on the nightside and downflow on the dayside increase with geomagnetic activity level. On the nightside, occurrence rate of ion upflow increases with enhanced solar wind and interplanetary magnetic field (IMF) drivers as well as southwestward local magnetic perturbations. The lack of correlation of upflow on the dayside with the solar wind and IMF parameters is because PFISR is usually equatorward of the dayside auroral zone. Occurrence rate of downflow does not show strong dependence on the solar wind and IMF conditions. However, it occurs much more frequently on the dayside when the IMFBy > 10 nT and the IMFBz < −10 nT, which we suggest is associated with the decaying of the dayside storm‐enhanced density (SED) and the SED plume.

     
    more » « less