We present a new MHD model for simulating the large-scale structure of the solar corona and solar wind under "steady state" conditions stemming from the Wang-Sheeley-Arge empirical model. The processes of turbulent heating in the solar wind are parameterized using a phenomenological, thermodynamical model with a varied polytropic index. We employ the Bernoulli integral to bridge the asymptotic solar wind speed with the assumed distribution of the polytropic index on the solar surface. We successfully reproduce the mass flux from Sun to Earth, the temperature structure, and the large-scale structure of the magnetic field. We reproduce the solar wind speed bimodal structure in the inner heliosphere. However, the solar wind speed is in a quantitative agreement with observations at 1 AU for solar maximum conditions only. The magnetic field comparison demonstrates that the input magnetogram needs to be multiplied by a scaling factor in order to obtain the correct magnitude at 1 AU.

[1] We observe ionospheric perturbations caused by the Tohoku earthquake and tsunami of March 11, 2011. Perturbations near the epicenter were found in measurements of ionospheric total electron content (TEC) from 1198 GPS receivers in the Japanese GEONET network. For the first time for this event, we compare these observations with the estimated magnitude and speed of a tsunami-driven atmospheric gravity wave, using an atmosphere-ionosphere-coupling model and a tsunami model of sea-surface height, respectively. Traveling ionospheric disturbances (TIDs) were observed moving away from the epicenter at approximate speeds of 3400 m/s, 1000 m/s and 200-300 m/s, consistent with Rayleigh waves, acoustic waves, and gravity waves, respectively. We focus our analysis on gravity waves moving south and east of the epicenter, since tsunamis propagating in the deep ocean have been shown to produce gravity waves detectable in ionospheric TEC in the past. Observed southeastward gravity wave perturbations, seen $60 min after the earthquake, are mostly between 0.5 to 1.5 TECU, representing up to $5% of the background vertical TEC (VTEC). Comparisons of observed TID gravity waves with the modeled tsunami speed in the ocean and the predicted VTEC perturbation amplitudes from an atmosphere-ionosphere-coupling model show the measurements and models to be in close agreement. Due to the dense GPS network and high earthquake magnitude, these are the clearest observations to date of the effect of a major earthquake and tsunami on the ionosphere near the epicenter. Such observations from a future realtime GPS receiver network could be used to validate tsunami models, confirm the existence of a tsunami, or track its motion where in situ buoy data is not available.

Recent advances in GPS data processing have demonstrated that ground-based GPS receivers are capable of detecting ionospheric TEC perturbations caused by surface-generated Rayleigh, acoustic and gravity waves. There have been a number of publications discussing TEC perturbations immediately following the M 9.0 Tohoku earthquake in Japan on March 11, 2011. Most investigators have focused on the ionospheric responses up to a few hours following the earthquake and tsunami. In our research, in addition to March 11, 2011 we investigate global ionospheric TEC perturbations a day before and after the event. We also compare indices of geomagnetic activity on all three days with perturbations in TEC, revealing strong geomagnetic storm conditions that are also apparent in processed GEONET TEC observations. In addition to the traveling ionospheric disturbances (TIDs) produced by the earthquake and tsunami, we also detect "regular" TIDs across Japan about 5 hours following the Tohoku event, concluding these are likely due to geomagnetic activity. The variety of observed TEC perturbations are consistent with tsunami-generated gravity waves, auroral activity, regular TIDs and equatorial fluctuations induced by increased geomagnetic activity. We demonstrate our capabilities to monitor TEC fluctuations using JPL's real-time Global Assimilative Ionospheric Model (GAIM) system. We show that a real-time global TEC monitoring network is able to detect the acoustic and gravity waves generated by the earthquake and tsunami. With additional real-time stations deployed, this new capability has the potential to provide real-time monitoring of TEC perturbations that could potentially serve as a plug-in to enhance existing early warning systems.

1] Objective quantification of model performance based on metrics helps us evaluate the current state of space physics modeling capability, address differences among various modeling approaches, and track model improvements over time. The Coupling, Energetics, and Dynamics of Atmospheric Regions (CEDAR) Electrodynamics Thermosphere Ionosphere (ETI) Challenge was initiated in 2009 to assess accuracy of various ionosphere/thermosphere models in reproducing ionosphere and thermosphere parameters. A total of nine events and five physical parameters were selected to compare between model outputs and observations. The nine events included two strong and one moderate geomagnetic storm events from GEM Challenge events and three moderate storms and three quiet periods from the first half of the International Polar Year (IPY) campaign, which lasted for 2 years, from March 2007 to March 2009. The five physical parameters selected were NmF2 and hmF2 from ISRs and LEO satellites such as CHAMP and COSMIC, vertical drifts at Jicamarca, and electron and neutral densities along the track of the CHAMP satellite. For this study, four different metrics and up to 10 models were used. In this paper, we focus on preliminary results of the study using ground-based measurements, which include NmF2 and hmF2 from Incoherent Scatter Radars (ISRs), and vertical drifts at Jicamarca. The results show that the model performance strongly depends on the type of metrics used, and thus no model is ranked top for all used metrics. The analysis further indicates that performance of the model also varies with latitude and geomagnetic activity level.Citation: Shim, J. S., et al. (2011), CEDAR Electrodynamics Thermosphere Ionosphere (ETI) Challenge for systematic assessment of ionosphere/thermosphere models: NmF2, hmF2, and vertical drift using ground-based observations, Space Weather, 9, S12003,

In an effort to quantitatively assess the current capabilities of Ionosphere/Thermosphere (IT) models, an IT model validation study using metrics was performed. This study is a main part of the CEDAR Electrodynamics Thermosphere Ionosphere (ETI) Challenge, which was initiated at the CEDAR workshop in 2009 to better comprehend strengths and weaknesses of models in predicting the IT system, and to trace improvements in ionospheric/thermospheric specification and forecast. For the challenge, two strong geomagnetic storms, four moderate storms, and three quiet time intervals were selected. For the selected events, we obtained four scores (i.e., RMS error, prediction efficiency, ratio of the maximum change in amplitudes, and ratio of the maximum amplitudes) to compare the performance of models in reproducing the selected physical parameters such as vertical drifts, electron and neutral densities, NmF2, and hmF2. In this paper, we present the results from comparing modeled values against space‐based measurements including NmF2 and hmF2 from the CHAMP and COSMIC satellites, and electron and neutral densities at the CHAMP satellite locations. It is found that the accuracy of models varies with the metrics used, latitude and geomagnetic activity level.