An Mw 9.3 earthquake originated in the Indian Ocean off the western coast of northern Sumatra at 00:58:53 Universal Time (UT) on 26 December 2004. Two giant ionospheric disturbances at 01:19 and 04:10 UT are observed by a network of digital Doppler sounders in Taiwan. The first disturbance excited mainly by Rayleigh waves, which consists of a packet of short‐period Doppler shift variations, results in vertical ionospheric fluctuations with a maximum velocity of about 70 m/s and displacement of about 200 m. The second disturbance, in a W‐shaped pulse propagating at a horizontal speed of 360 ± 70 m/s, is attributable to coupling of the atmospheric gravity waves (AGW) excited by broad crustal uplift together with the following big tsunami waves around the earthquake source zone. The accompanying ionosonde data suggest that the AGW in the atmosphere may have caused the ionosphere to move up and down by about 40 km.

This paper presents our effort to assimilate FORMOSAT‐3/COSMIC (F3/C) GPS Occultation Experiment (GOX) observations into the National Center for Atmospheric Research (NCAR) Thermosphere Ionosphere Electrodynamics General Circulation Model (TIE‐GCM) by means of ensemble Kalman filtering (EnKF). The F3/C electron density profiles (EDPs) uniformly distributed around the globe which provide an excellent opportunity to monitor the ionospheric electron density structure. The NCAR TIE‐GCM simulates the Earth's thermosphere and ionosphere by using self‐consistent solutions for the coupled nonlinear equations of hydrodynamics, neutral and ion chemistry, and electrodynamics. The F3/C EDP are combined with the TIE‐GCM simulations by EnKF algorithms implemented in the NCAR Data Assimilation Research Testbed (DART) open‐source community facility to compute the expected value of electron density, which is ‘the best’ estimate of the current ionospheric state. Assimilation analyses obtained with real F3/C electron density profiles are compared with independent ground‐based observations as well as the F3/C profiles themselves. The comparison shows the improvement of the primary ionospheric parameters, such as NmF2 and hmF2. Nevertheless, some unrealistic signatures appearing in the results and high rejection rates of observations due to the applied outlier threshold and quality control are found in the assimilation experiments. This paper further discusses the limitations of the model and the impact of ensemble member creation approaches on the assimilation results, and proposes possible methods to avoid these problems for future work.

[1] This paper presents a statistical study of the pre-earthquake ionospheric anomaly by using the total electron content (TEC) data from the global ionosphere map. A total of 736 M ≥ 6.0 earthquakes in the global area during 2002-2010 are selected. The anomaly day is first defined. Then the occurrence rates of abnormal days for both the days within 1-21 days prior to the earthquakes (P E ) and the background days (P N ) are calculated. The results show that the values of P E depend on the earthquake magnitude, the earthquake source depth, and the number of days prior to the earthquake. The P E is larger for earthquakes with greater magnitude and lower depth and for days closer to the earthquakes. The results also show that the occurrence rate of anomaly within several days before the earthquakes is overall larger than that during the background days, especially for the large-magnitude and low-depth earthquakes. These results indicate that the anomalous behavior of TEC within just a few days before the earthquakes is related with the forthcoming earthquakes with high probability.

We present a statistical analysis of the occurrence probability of equatorial spread F irregularities measured by the Communication/Navigation Outage Forecasting System satellite during 2008-2012. We use different criteria (plasma density perturbations, ΔN, and relative density perturbations, ΔN/N 0 ) to identify the occurrence of ionospheric irregularities. The purpose of this study is to determine whether the occurrence probability of irregularities is the same for different criteria, whether the patterns of irregularity occurrence vary with solar activity and with local time, and how the patterns of irregularity occurrence are correlated with ionospheric scintillation. It is found that the occurrence probability of irregularities and its variation with local time are significantly different when different identification criteria are used. The occurrence probability based on plasma density perturbations is high in the evening sector and becomes much lower after midnight. In contrast, the occurrence probability based on relative density perturbations is low in the evening sector but becomes very high after midnight in the June solstice. We have also compared the occurrence of ionospheric irregularities with scintillation. The occurrence pattern of the S4 index and its variation with local time are in good agreement with the irregularity occurrence based on plasma density perturbations but are significantly different from those based on relative density perturbations. This study reveals that the occurrence pattern of equatorial ionospheric irregularities varies with local time and that only the occurrence probability of irregularities based on plasma density perturbations is consistent with the occurrence of scintillation at all local times.