Prepared Feb. 15, 2019 by M. Burleigh (mburleig@umich.edu)


Data Archive to support related publication:
Burleigh, M., Zettergren, M., Lynch, K., Lessard, M., Moen, J., Clausen, L., Kenward, D., Hysell, D., and Liemohn, M. (2019). Transient ionospheric upflow driven by poleward moving auroral forms observed during the Rocket Experiment for Neutral Upwelling 2 (RENU2) campaign. Geophysical Research Letters. (Submitted).


Overview: 
The Geospace Environment Model of Ion-Neutral Interactions with Transverse Ion Acceleration (GEMINI-TIA) has been used here to simulate poleward moving auroral forms (PMAFs) in the high-latitude ionosphere as observed during the Rocket Experiment for Neutral Upwelling 2 (RENU2) campaign. Brightness from an all-sky imager has been used as a spatiotemporal constraint for auroral inputs selected from in situ rocket measurements which drive the model. This method allows for realistic ionospheric forcing that is not captured in traditional "on-off" methods of describing PMAFs. Transient forcing and steady forcing simulations have been generated for comparison and discussion in the related publication.


Software:
The .mat files were generated using MATLAB 2018a and tested for compatibility with MATLAB versions 2012, 2014, 2016, 2017a,b, and 2018a. This file format should also be compatible with Octave (a free alternative to MATLAB).


File Description:
Model outputs are archived here and include:

dataRENU2_SS_all.mat  --  This is the control simulation where the model was simply turned on and allowed to run. GEMINI-TIA has been developed with dynamic background conditions that include ambient effects such as photoionization, collisional cooling, etc. that realistically evolve over time generating natural diurnal variations. 

dataRENU2_steady_cusp_all.mat  --  This simulation used a representative form of steady cusp forcing. Model inputs include total energy flux (0.75 mW/m^2) and characteristic energy (100 eV) which were applied constantly for 20 minutes from 6:45 to 7:05 UT using a latitudinal Gaussian envelope, centered on ~77 {degrees}, with a half-width of ~0.6 {degrees} to create latitudinal structure. 

dataRENU2_fullstrength_all.mat -- This simulation used all-sky imager normalized brightness data combined with DC electric field (8 mV/m), total energy flux (0.75 mW/m^2), and characteristic energy (100 eV) model inputs to provide realistic spatiotemporal variability.


Variables in use in these files:
nst       Density (m^-3) [vertical (x1) x horizontal (x2) x specie x time]
vsx1t     Field-aligned velocity (m/s) [vertical x horizontal x specie x time]
Tspart    Parallel temperature (K) [vertical x horizontal x specie x time]
Tsperpt   Perpendicular temperature (K) [vertical x horizontal x specie x time] 
simdate   Each row contains a date/time organized as [Day,Month,Year,Hour (in UT),Min,Sec]
timeUT    Time in UT
xg        Structure containing model grid information. Structure members include: xg.x1 magnetic field-aligned position variable, xg.x2 horizontal position variable, xg.x3 horizontal position constant, xg.h1, xg.h2, and xg.h3 are the corresponding the metric factors, lx stores the number of grid points, xg.dx** are the finite differences (1,2,3 corresponds to the dimension, f,b,h indicates forward, backward, or half finite difference), xg.glat,glon are the latitudes and longitudes (degrees geographic) of each grid point, xg.alt is the altitude of each grid point, xg.r,theta,phi for each grid point are the radial distance (from the center of the Earth), magnetic colatitude (in radians), and magnetic longitude (in radians), nullpts is a logical that indicates what portion of the grid (if necessary) is treated as a bottom boundary to avoid locations with negative altitude, and Bmag is the magnitude of the local magnetic field.

The specie order is O^+, NO^+, N_2^+, O_2^+, N^+, H^+, e^-


For more information on GEMINI-TIA please see: 
Burleigh, M., & Zettergren, M. (2017). Anisotropic fluid modeling of ionospheric upflow: Effects of low-altitude anisotropy and thermospheric winds. Journal of Geophysical Research: Space Physics. doi: 10.1002/2016JA023329

If you would like to use GEMINI-TIA please contact M. Burleigh (mburleig@umich.edu). Effort is currently underway to open source this model. 


Use and Access:
This data set is made available under a Attribution-NonCommercial 4.0 International license (CC BY-NC 4.0).


Citation for this Data Archive:
Burleigh, M. (2019). Modeled transient ionospheric upflow as driven by poleward moving auroral forms observed during the Rocket Experiment for Neutral Upwelling 2 (RENU2) campaign. [Data set]. University of Michigan Deep Blue Data Repository. https://doi.org/10.7302/pcrt-dj64