No abstract

[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.

[1] Recent observations have revealed large F-region electron density perturbations ($100%) and total electron content (TEC) perturbations ($30%) that appear to be correlated with tsunamis. The characteristic speed and horizontal wavelength of the disturbances are $200 m/s and $400 km. We describe numerical simulations using our spectral full-wave model (SFWM) of the upward propagation of a spectrum of gravity waves forced by a tsunami, and the interaction of these waves with the F-region ionosphere. The SFWM describes the propagation of linear, steady-state acoustic-gravity waves in a nonisothermal atmosphere with the inclusion of eddy and molecular diffusion of heat and momentum, ion drag, Coriolis force, and height-dependent mean winds. The tsunami is modeled as a deformation of our model lower boundary traveling at the shallow water wave speed of 200 m/s with a maximum vertical displacement of 50 cm and described by a modified Airy function in the horizontal direction. The derived vertical velocity spectrum at the surface describes the forcing at the lower boundary of the SFWM. A steady-state 1-D ionospheric perturbation model is used to calculate the electron density and TEC perturbations. The molecular diffusion strongly damps the waves in the topside (>300-km altitude) ionosphere. In spite of this, the F-region response is large, with vertical displacements of $2 to 5 km and electron density perturbations of $100%. Mean winds have a profound effect on the ability of the waves to propagate into the F-region ionosphere.Citation: Hickey, M. P., G. Schubert, and R. L. Walterscheid (2009), Propagation of tsunami-driven gravity waves into the thermosphere and ionosphere,

Abstract. From analyzing nine months of airglow imaging observations of atmospheric gravity waves (AGWs) over Adelaide, Australia (35øS) [Walterscheid et al., 1999] have proposed that many of the quasi-monochromatic waves seen in the images were primarily thermally ducted. Here are presented 15 months of observations, from February 1996 to May 1997, for AGW frequency and propagation direction from a northern latitude site, Urbana Illinois (40øN). As Adelaide, Urbana is geographically distant from large orographic features. Similar to what was found in Adelaide, the AGWs seem to originate from a preferred location during the time period around summer solstice. In conjunction with these airglow data there exists MF radar data to provide winds in the 90 km region and near-simultaneous lidar data which provide a temperature climatology. The temperature data have previously been analyzed by States and Gardner [2000]. The temperature and wind data are used here in a full wave model analysis to determine the characteristics of the wave ducting and wave reflection during the 15 month observation period. This model analysis is applied to this and another existing data set recently described by Nakamura et al. [1999]. It is shown that the existence of a thermal duct around summer solstice can plausibly account for our observations. However, the characteristics of the thermal duct and the ability of waves to be ducted is also greatly dependent on the characteristics of the background wind. A simple model is constructed to simulate the trapping of these waves by such a duct. It is suggested that the waves seen over Urbana originate no more than a few thousand kilometers from the observation site. IntroductionThe ducting or trapping of atmospheric gravity waves (AGWs) has long been a subject of interest [e.g., Pitteway and

[1] Ground-based Global Positioning System (GPS) measurements of ionospheric total electron content (TEC) show variations consistent with atmospheric internal gravity waves caused by ocean tsunamis following two recent seismic events: the Samoa earthquake of 29 September 2009 and the Chile earthquake of 27 February 2010. Both earthquakes produced ocean tsunamis that were destructive to coastal communities near the epicenters, and both were observed in tidal gauge and buoy measurements throughout the Pacific Ocean. We observe fluctuations in TEC correlated in time, space, and wave properties with these tsunamis using the Jet Propulsion Laboratory's Global Ionospheric Mapping software. These TEC measurements were band-pass filtered to remove ionospheric TEC variations with wavelengths and periods outside the typical range for tsunamis. Observable variations in TEC appear correlated with the tsunamis in some locations (Hawaii and Japan), but not in others (Southern California or near the epicenters). Where variations are observed, the typical amplitude tends to be ∼0.1-0.2 TEC units for these events, on the order of ∼1% of the background TEC value. These observations are compared to estimates of expected tsunami-driven TEC variations produced by Embry Riddle Aeronautical University's Spectral Full Wave Model, an atmosphere-ionosphere coupled model, and are found to be in good agreement. Significant TEC variations are not always seen when a tsunami is present, but in these two events the regions where a strong ocean tsunami was observed coincided with clear TEC observations, while a lack of clear TEC observations coincided with smaller sea surface height amplitudes. There exists the potential to apply these detection techniques to real-time GPS TEC data, providing estimates of tsunami speed and amplitude that may be useful for early warning systems.

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.

We use a nonlinear, fully compressible, two-dimensional numerical model to study the effects of dissipation on gravity wave packet spectra in the thermosphere. Numerical simulations are performed to excite gravity wave packets using either a time-dependent vertical body forcing at the bottom boundary or a specified initial wave perturbation. Three simulation case studies are performed to excite (1) a steady state monochromatic wave, (2) a spectrally broad wave packet, and (3) a quasi-monochromatic wave packet. In addition, we analyze (4) an initial condition simulation with an isothermal background. We find that, in cases where the wave is not continually forced, the dominant vertical wavelength decreases in time, predominantly due to a combination of refraction from the thermosphere and dissipation of the packets' high frequency components as they quickly reach high altitude. In the continually forced steady state case, the dominant vertical wavelength remains constant once initial transients have passed. The vertical wavelength in all simulations increases with altitude above the dissipation altitude (the point at which dissipation effects are greater than the wave amplitude growth caused by decreasing background density) at any fixed time. However, a shift to smaller vertical wavelength values in time is clearly exhibited as high-frequency, long vertical wavelength components reach high altitudes and dissipate first, to be replaced by slower waves of shorter vertical wavelength. Results suggest that the dispersion of a packet significantly determines its spectral evolution in the dissipative thermosphere.

No abstract