First storm‐time plasma velocity estimates from high‐resolution ionospheric data assimilation

This paper uses data assimilation to estimate ionospheric state during storm time at subdegree resolution. We use Ionospheric Data Assimilation Four‐Dimensional (IDA4D) to resolve the three‐dimensional time‐varying electron density gradients of the storm‐enhanced density poleward plume. By Estimating Model Parameters from Ionospheric Reverse Engineering (EMPIRE), we infer the three‐dimensional plasma velocity from the densities. EMPIRE estimates of ExB drift are made by correcting the Weimer 2000 electric potential model. This is the first time electron densities derived from GPS total electron content (TEC) data are being used to estimate field‐aligned and field‐perpendicular drifts at such high resolution, without reference to direct drift measurements. The time‐varying estimated electron densities are used to construct the ionospheric spatial decorrelation in vertical total electron content (TEC) on horizontal scales of less than 100 km. We compare slant TEC (STEC) estimates to actual STEC GPS observations, including independent unassimilated data. The IDA4D density model of the extreme ionospheric storm on 20 November 2003 shows STEC delays of up to 210 TEC units, comparable to the STEC of the GPS ground stations. Horizontal drifts from EMPIRE are predicted to be northwestward within the storm‐enhanced density plume and its boundary, turning northeast at high latitudes. These estimates compare favorably to independent Assimilative Mapping of Ionospheric Electrodynamics‐assimilated high‐latitude ExB drift estimates. Estimated and measured Defense Meteorological Satellite Program in situ drifts differ by a factor of 2–3 and in some cases have incorrect direction. This indicates that significant density rates of change and more accurate accounting for production and loss may be needed when other processes are not dominant.

[1]  Per Enge,et al.  Assessment of Ionosphere Spatial Decorrelation for Global Positioning System-Based Aircraft Landing Systems , 2007 .

[2]  Anthony J. Mannucci,et al.  Automated daily processing of more than 1000 ground‐based GPS receivers for studying intense ionospheric storms , 2005 .

[3]  Per Enge,et al.  A Comprehensive Ionosphere Storm Data Analysis Method to Support LAAS Threat Model Development , 2005 .

[4]  Y. Otsuka,et al.  Physical mechanisms of the ionospheric storms at equatorial and higher latitudes during the recovery phase of geomagnetic storms , 2013 .

[5]  G. J. Bailey,et al.  Response of the ionosphere to super geomagnetic storms: Observations and modeling , 2008 .

[6]  D. Weimer,et al.  An improved model of ionospheric electric potentials including substorm perturbations and application to the Geospace Environment Modeling November 24, 1996, event , 2001 .

[7]  M. Yamamoto,et al.  A statistical study of the response of the dayside equatorial F2 layer to the main phase of intense geomagnetic storms as an indicator of penetration electric field , 2011 .

[8]  Thomas L. Gaussiran,et al.  Ionospheric Data Assimilation Three‐Dimensional (IDA3D): A global, multisensor, electron density specification algorithm , 2004 .

[9]  Gary S. Bust,et al.  History, current state, and future directions of ionospheric imaging , 2008 .

[10]  R. W. Meggs,et al.  Four‐dimensional GPS imaging of space weather storms , 2007 .

[11]  G. Bust,et al.  Error Propagation in Ionospheric Image-based Parameter Estimation , 2010 .

[12]  H. Yeh,et al.  Satellite observations of electric fields in the South Atlantic anomaly region during the July 2000 magnetic storm , 2005 .

[13]  C. Valladares,et al.  Modeling ionospheric super‐fountain effect based on the coupled TIMEGCM‐SAMI3 , 2013 .

[14]  G. Bust,et al.  Tracking of polar cap ionospheric patches using data assimilation , 2007 .

[15]  G. Bust,et al.  Neutral wind estimation from 4‐D ionospheric electron density images , 2009 .

[16]  Jonathan J. Makela,et al.  Combined Ionospheric Campaign 1: Ionospheric tomography and GPS total electron count (TEC) depletions , 2000 .

[17]  G. Bust,et al.  Deducing storm time F region ionospheric dynamics from 3‐D time‐varying imaging , 2011 .

[18]  Per Enge,et al.  Ionospheric Threat Parameterization for Local Area Global-Positioning-System-Based Aircraft Landing Systems , 2010 .

[19]  G. Bust,et al.  GNSS Imaging-derived Dynamics of Ionospheric Storm Transition Regions , 2011 .

[20]  A. Coster,et al.  Ionospheric and thermospheric variations associated with prompt penetration electric fields , 2012 .

[21]  Larry J. Paxton,et al.  Global thermosphere‐ionosphere response to onset of 20 November 2003 magnetic storm , 2006 .

[22]  R. W. Schunk,et al.  Solar-terrestrial energy program : handbook of ionospheric models , 1996 .

[23]  S. Wing,et al.  A new magnetic coordinate system for conjugate studies at high latitudes , 1989 .

[24]  Per Enge,et al.  Impact and mitigation of ionospheric anomalies on ground‐based augmentation of GNSS , 2009 .

[25]  S. Kawamura,et al.  A physical mechanism of positive ionospheric storms at low latitudes and midlatitudes , 2010 .

[26]  Cathryn N. Mitchell,et al.  Observations of the F region height redistribution in the storm‐time ionosphere over Europe and the USA using GPS imaging , 2006 .

[27]  G. Bust,et al.  Mapping the Time‐Varying Distribution of High‐Altitude Plasma During Storms , 2013 .

[28]  K. Shiokawa,et al.  Super plasma fountain and equatorial ionization anomaly during penetration electric field , 2009 .

[29]  W. Rideout,et al.  Multiradar observations of the polar tongue of ionization , 2005 .

[30]  D. Weimer,et al.  Improved Ionospheric Electrodynamic Models and Application to Calculating Joule Heating Rates , 2005 .