Inverse Modeling Using a Novel Scintillation Model Sigma to Characterize High Latitude Irregularities During a Storm and a Substorm

Monday, 15 December 2014
Kshitija Deshpande1, Gary S Bust2, C. Robert Clauer3, Hyomin Kim1 and Daniel R Weimer4, (1)Virginia Polytechnic Institute and State University, Blacksburg, VA, United States, (2)JHU Applied Physics Lab, Laurel, MD, United States, (3)National Institute of Aerospace, Hampton, VA, United States, (4)Virginia Tech, Department of Electrical and Computer Engineering, Blacksburg, VA, United States
We have developed a high fidelity ``Satellite-beacon Ionospheric-scintillation Global Model of the upper Atmosphere" (SIGMA) which is a full 3D EM wave propagation model to simulate GPS scintillations globally. We demonstrate in this work that the results from this model can form a basic framework on the use of inverse method to understand the physics of high latitude irregularities using GPS scintillations. We are using SIGMA and an inverse method to understand the physics of the irregularities observed with GPS receivers from six different inter-hemispheric high latitude stations during a geomagnetic storm on 9 March 2012, and from Autonomous Adaptive Low-Power Instrument Platform (AAL-PIP) Antarctic stations during a substorm on 9 January 2014. We utilize ancillary observations from SuperDARN, ISRs, riometers etc. to obtain some of the input parameters of SIGMA. Further, we implement a uniform-grid SIGMA simulation or a non-linear optimization of the model to obtain the rest of the unknowns that give us the best fit with data. The input parameters of SIGMA thus derived represent the physical and propagation parameters related to the physics of the irregularity that produced those GNSS scintillations. Some of our findings from this investigation include that the spectral indices and outer scales for ionospheric heights of 350~km are higher than those at 120~km. The best fits we obtained from our inverse method mostly agree with the observations except for some cases, which might be because the spectral model we use is insufficient to describe irregularity physics. We need more auxiliary data in order to facilitate the possibility of accomplishing a unique solution to the inverse problem.