Search

The 2022 eruption of the Hunga Tonga-Hunga Ha'apai volcano excited an atmospheric Lamb wave, which induced a fast-traveling tsunami. This tsunami was driven by the pressure-forced wave traveling at the speed of the Lamb wave and, thus, was much faster than conventional tsunamis. This was the first case in which ocean bottom monitoring systems widely observed an air pressure-induced tsunami. We found that the pressure-forced waves split and generated ocean gravity waves after passing the Japan Trench based on the S-net data. Our simulations show that changes in water depth can amplify or decrease the pressure-forced wave. Simultaneously, an ocean gravity wave is generated due to the conservation of water volume. Because the ocean gravity wave was slower than the pressure-forced wave near Japan, it was separated from, and traveled behind, the pressure-forced wave. We explained the wave separation phenomenon and reproduced the waveforms of different splitting stages observed by the stations near the Japan Trench. Graphical Abstract

Using meteor radar, radiosonde observations and MERRA-2 reanalysis data from 12 August to 31 October 2006, we report a dynamical coupling from the tropical lower atmosphere to the mesosphere and lower thermosphere through a quasi-27-day intraseasonal oscillation (ISO). It is interesting that the quasi-27-day ISO is observed in the troposphere, stratopause and mesopause regions, exhibiting a three-layer structure. In the MLT, the amplitude in the zonal wind increases from about 4 ms−1 at 90 km to 15 ms−1 at 100 km, which is different from previous observations that ISOs occurs generally in winter with an amplitude peak at about 80–90 km, and then are rapidly weakened with increasing height. Outgoing longwave radiation (OLR) and specific humidity demonstrate that there is a quasi-27-day periodicity in convective activity in the tropics, which causes the ISO of the zonal wind and gravity wave (GW) activity in the troposphere. The upward propagating GWs are further modulated by the oscillation in the troposphere and upper stratosphere. As the GWs propagate to the MLT, the quasi-27-day oscillation in the wind field is induced with a clear phase opposite to that in the lower atmosphere through instability and dissipation of these modulated GWs. Wavelet analysis shows that the quasi-27-day variability in the MLT appears as a case event rather than a persistent phenomenon, and has not a clear corresponding relation with the solar rotation effect within 1 year of observations.

In order to reveal solar activity dependence of the medium-scale traveling ionospheric disturbances (MSTIDs) at midlatitudes, total electron content (TEC) data obtained from a Global Positioning System (GPS) receiver network in Japan during 22 years from 1998 to 2019 were analyzed. We have calculated the detrended TEC by subtracting the 1-h running average from the original TEC data for each satellite and receiver pair, and made two-dimensional TEC maps of the detrended TEC with a spatial resolution of 0.15° × 0.15° in longitude and latitude. We have investigated MSTID activity, defined as $$\delta I/\overline{I}$$ δ I / I ¯ , where $$\delta I$$ δ I and $$\overline{I}$$ I ¯ are standard deviation of the detrended TEC and the average vertical TEC within the area of 133.0°–137.0° E and 33.0°–37.0° N for 1 h, respectively. From each 2-h time series of the detrended TEC data within the same area as the MSTID activity, auto-correlation functions (ACFs) of the detrended TEC were calculated to estimate the horizontal propagation velocity and direction of the MSTIDs. Statistical results of the MSTID activity and propagation direction of MSTIDs were consistent with previous studies and support the idea that daytime MSTIDs could be caused by atmospheric gravity waves, and that nighttime MSTIDs were caused by electro-dynamical forces, such as the Perkins instability. From the current long-term observations, we have found that the nighttime MSTID activity and occurrence rate increased with decreasing solar activity. For the daytime MSTID, the occurrence rate increased with decreasing solar activity, whereas the MSTID activity did not show distinct solar activity dependence. These results suggest that the secondary gravity waves generated by dissipation of the primary gravity waves propagating from below increase under low solar activity conditions. The mean horizontal phase velocity of the MSTIDs during nighttime did not show a distinct solar activity dependence, whereas that during daytime showed an anticorrelation with solar activity. The horizontal phase velocity of the daytime MSTIDs was widely distributed from 40 to 180 m/s under high solar activity conditions, whereas it ranged between 80 and 200 m/s, with a maximum occurrence at 130 m/s under low solar activity conditions, suggesting that gravity waves with low phase velocity could be dissipated by high viscosity in the thermosphere under low solar activity conditions.

The reason for stationary gravity waves at Venus’ cloud top to appear mostly at low latitudes in the afternoon is not understood. Since a neutral layer exists in the lower part of the cloud layer, the waves should be affected by the neutral layer before reaching the cloud top. To what extent gravity waves can propagate vertically through the neutral layer has been unclear. To examine the possibility that the variation of the neutral layer thickness is responsible for the dependence of the gravity wave activity on the latitude and the local time, we investigated the sensitivity of the vertical propagation of gravity waves on the neutral layer thickness using a numerical model. The results showed that stationary gravity waves with zonal wavelengths longer than 1000 km can propagate to the cloud-top level without notable attenuation in the neutral layer with realistic thicknesses of 5–15 km. This suggests that the observed latitudinal and local time variation of the gravity wave activity should be attributed to processes below the cloud. An analytical approach also showed that gravity waves with horizontal wavelengths shorter than tens of kilometers would be strongly attenuated in the neutral layer; such waves should originate in the altitude region above the neutral layer.

Propagation of medium-scale gravity waves (GWs) in the thermosphere/ionosphere is observed remotely, using multi-frequency and multi-point continuous Doppler sounding system located in the western part of Czechia. Reflection heights of the sounding radio waves are determined from a nearby ionosonde. Phase velocity vectors of GWs are calculated from time/phase delays between signals corresponding to different transmitter–receiver pairs that reflect in the ionosphere at different locations. As various frequencies reflect at different heights, reflection points of radio signals are separated both horizontally and vertically, and the investigation of GW propagation in the ionosphere is performed in three dimensions. Results obtained for two 1-year periods representing the solar maximum (July 2014–June 2015) and current solar minimum (September 2018–August 2019) are presented. It is shown that GWs in the ionosphere usually propagated with wave vectors directed obliquely downward. A statistical distribution of wave vector elevation angles is presented. A model of neutral winds is used to estimate the wave characteristics in the wind-rest frame. It is found that the distribution of elevation angles is narrower in the wind-rest frame than in the Earth frame. Seasonal and diurnal changes of propagation directions and attenuations of GWs are discussed. The wind-rest frame wavelengths of the analyzed GWs were usually from ~ 80 to 300 km, and the propagation velocities were mostly between ~ 100 and ~ 220 m/s.

Polar mesosphere winter echoes (PMWE) were observed at 70 km over Tromsø, Norway, on 8 January 2014 using the tristatic configuration of the European incoherent scatter VHF radar. For the interval 11:00–13:00 UT where the strongest patch of PMWE of about 6-min duration was detected, the spectra of the received signal were analysed for the Tromsø site and altitude profiles of spectral parameters were derived. For the remote sites Kiruna and Sodankylä, the Doppler velocities and their vertical shear were determined by using the measured autocorrelation functions. Ducted gravity waves with periods of 5–10 min were found in the vertical wind velocity between 66 and 81 km altitudes. The duct might be formed around 70 and 77 km altitude where horizontal wind maxima were observed with the Kiruna receiver. However, we did not find any close relation between wind shear at 70 km altitude and PMWE at the same height: the wind shear was present for 2 h, but PMWE for only 6 min. Enhanced spectral width in the vertical Tromsø beam was observed for the PMWE patch. We discussed these experimental findings in relation to the winter echo generation mechanism. Our conclusion is that the presence of patchy negatively charged small-sized dust might explain the observations although a gravity wave breaking mechanism cannot be completely rejected.

Numerical simulations are performed to simulate atmospheric perturbations observed at ionospheric heights just after the 2011 off the Pacific coast of Tohoku Earthquake. A time-dependent, two-dimensional, nonlinear, non-hydrostatic, compressible and neutral, numerical model is developed to reproduce the atmospheric perturbations. An impulsive upward surface motion is assumed as the source of the perturbations. Simulated atmospheric perturbations at 300-km altitude show remarkable agreement with oscillations observed in the ionospheric total electron content (TEC) when the source width is about 250 km. In the vicinity of the source, the acoustic resonance modes between the ground surface and the lower thermosphere are dominant. They have three dominant frequencies for the interval between 20 and 60 min after the impulsive input. The perturbation with the maximum amplitude has a frequency of 4.4 mHz. The other dominant modes have frequencies of 3.6 and 5.1 mHz. The beats between the dominant modes are also seen. In the distance, the gravity modes are dominant. The horizontal phase velocities are about 220 to 300 m/s, and the horizontal wavelengths are about 200 to 400 km. The good agreement between the simulation and the observations indicates that ionospheric oscillations generated by the earthquake are mainly due to the motion of the neutral atmosphere.

Accepted: 2011-06-11, 資料番号: SA1003052000

We performed an extensive experimental campaign (the spread F Experiment, or SpreadFEx) from September to November 2005 to attempt to define the role of neutral atmosphere dynamics, specifically wave motions propagating upward from the lower atmosphere, in seeding equatorial spread F and plasma bubbles extending to higher altitudes. Campaign measurements focused on the Brazilian sector and included ground-based optical, radar, digisonde, and GPS measurements at a number of fixed and temporary sites. Related data on convection and plasma bubble structures were also collected by GOES 12 and the GUVI instrument aboard the TIMED satellite. Initial results of our analyses of SpreadFEx and related data indicate 1) extensive gravity wave (GW) activity apparently linked to deep convection predominantly to the west of our measurement sites, 2) the presence of small-scale GWactivity confined to lower altitudes, 3) larger-scaleGWactivity apparently penetrating to much higher altitudes suggested by electron density and TEC fluctuations in the E and F regions, 4) substantial GW amplitudes implied by digisonde electron densities, and 5) apparent direct links of these perturbations in the lower F region to spread F and plasma bubbles extending to much higher altitudes. Related efforts with correlative data are defining 6) the occurrence and locations of deep convection, 7) the spatial and temporal evolutions of plasma bubbles, the 8) 2D (height-resolved) structures of plasma bubbles, and 9) the expected propagation of GWs and tides from the lower atmosphere into the thermosphere and ionosphere.

Doppler radar observations have made a large contribution towards improving our understanding of middle- and upper-atmosphere dynamics. This radar echo is mainly due to radiowave scattering by atmospheric turbulence, which occurs almost universally throughout the atmosphere as a result of gravity wave (GW) breaking. In the strictness sense, the radar pulse is back-scattered by radio refractive index (RRI) perturbations caused by turbulence whose size along the radar beam is half of the radar wavelength. Understanding of the validity of the MST (mesosphere, stratosphere and troposphere)-radar observation requires a precise knowledge of how RRI perturbations are formed and behave. A basic property of the RRI of the atmosphere is that it has a static vertical gradient in terms of distribution and the turbulence velocity is divergence-free. The RRI depends on three atmospheric properties—humidity, air density and electron density. These three elements each have a static vertical gradient along which turbulence transports each element, hereby perturbing the density distribution of each element; that is to say, it produces the RRI irregularity for scattering the radar pulse. Since the turbulence motion and winds are divergence-free, PRI perturbation co-moves with the total air flow, i.e. turbulence motion plus winds. This can be true even in the mesosphere where the perturbation is controlled electromagnetically. Co-movement of turbulence with local winds has been shown from a comparison of observations with radars and radiosondes. In addition to tracking turbulence as wind-tracers, MST-radar observations provide important data in the study of atmosphere turbulence dynamics.

Significant attention is being paid in recent times by several observational and modeling studies to quantify the spatial and temporal variabilities of diurnal tide in the mesosphere and lower thermosphere (MLT) region. These variabilities are ascribed to spatial and temporal variations in the tidal forcing or interactions between the propagating tides and background wind, planetary waves or gravity waves. The present work makes use of simultaneous ground-based radar wind observations of different durations from five equatorial/low latitude sites in the Indian, Indonesian and Pacific sectors: Tirunelveli (8.7°N, 77.8°E), Jakarta (6.4°S, 106.7°E), Pontianak (0.03°N, 109°E), Kauai (22°N, 160°W) and Christmas Island (2°N, 157°W). This study delineates the longitudinal differences in the tidal characteristics in (i) interannual time scales over Tirunelveli and Kauai during 1993–2002, (ii) seasonal time scales over Christmas Island, Jakarta and Tirunelveli for the years 1993–1997 and (iii) shorter than seasonal time scales over Christmas Island, Pontianak and Tirunelveli during 1996–1997. An important observational feature noticed in this work is the differing behavior of the long-term tidal fields over Tirunelveli and Kauai. The monthly tidal amplitudes over Tirunelveli reveal a strong QBO signature whereas a similar, strong QBO signal could not be traced in the long-term observations from Kauai.

We review measurements of nighttime atmospheric/ionospheric waves in the upper atmosphere in Japan, Indonesia, and Australia, using all-sky airglow imagers of optical mesosphere thermosphere imagers (OMTIs). The imagers observe two-dimensional patterns of airglow emissions from oxygen (wavelength: 557.7 nm) and hydroxyl (OH) (near-infrared band) in the mesopause region (80–100 km) and from oxygen (630.0 nm) in the thermosphere/ionosphere (200–300 km). Several statistical studies were done to investigate propagation characteristics of small-scale (less than 100 km) gravity waves in the mesopause region and medium-scale traveling ionospheric disturbances (MSTIDs, ∼100–1,000 km) in the thermosphere/ionosphere. Clear seasonal variations of occurrence and propagation directions were reported for these waves. The propagation directions in the mesopause region are controlled by wind filtering, ducting processes and relative location to the wave sources in the troposphere. Poleward-propagating waves tend to be observed in the summer in the mesopause region at several stations, suggesting that mesospheric gravity waves are generated by intense convective activity in the equatorial troposphere. On the other hand, systematic equatorward and westward motions were observed for all seasons for nighttime MSTIDs in the midlatitude ionosphere with geomagnetic conjugacy between the northern and southern hemispheres. Ionospheric instabilities may play important role for the generation and propagation of these MSTIDs. We also give an example of simultaneous observation of quasi-periodic southward-moving waves in the mesopause region and in the thermosphere at the geographic equator. From these results, we discuss mean wind acceleration by mesospheric gravity waves and penetration of gravity waves from the mesosphere to the thermosphere.

Using Global Positioning System (GPS) satellites, we have been conducting equatorial ionospheric scintillation observations at Kototabang, Indonesia since January 2003. Scintillations caused by equatorial plasma bubbles appear between 2000 and 0100 LT in equinoctial months with a seasonal asymmetry, and their activity decreases with decreasing solar activity. A comparison between scintillation index (S 4) and Earth’s brightness temperature (T bb) variations suggests that the scintillation activity can be related to tropospheric disturbances over the Indian Ocean to the west of Kototabang. To understand better the reasons of day-to-day variability of S 4, we analyze S 4, T bband lower thermospheric neutral wind ("Equation missing") data. The results show that S 4 fluctuates with periods of about 2.5, 5, 8, 14 and 25 days, possibly due to atmospheric waves from below and that similar periods are also found in the T bband "Equation missing" variations. Using a general circulation model, we made numerical simulations to determine the behavior of neutral wind in the equatorial thermosphere. The results indicate the following: (1) 2- to 20-day waves dissipate rapidly above about an altitude of 125 km, and 0.5- to 3-hour waves become predominant above 100 km, (2) zonal winds above 200 km altitude are, on the whole, eastward during sunset-sunrise, (3) zonal wind patterns due to short-period (1–4 h) atmospheric gravity waves (AGWs) above 120 km altitude change day by day, exhibit wavy structures with scale lengths of about 30–1000 km and, as a whole, move eastward at about 100−1 while changing patterns over time. These simulations suggest that the Rayleigh-Taylor instability responsible for plasma bubble generation can be seeded by AGWs with short periods of about 0.5–3 h, and that background conditions necessary for this instability are modulated by planetary-scale atmospheric waves propagating up to an altitude of about 120 km from below.

Using a general circulation model that contains the region from the ground surface to the upper thermosphere, we have examined characteristics of gravity waves in the equatorial thermosphere. At an altitude of 150 km, the dominant periods of gravity waves for zonal wave number 20 (zonal wavelength λ x ≈ 2000 km), 40 (λ x ≈ 1000 km) and 80 (λ x ≈ 500 km) are 3, 1.5 and 1 h, respectively. For individual zonal wave numbers, the corresponding dominant period becomes shorter at higher altitudes due to dissipation processes in the thermosphere, such as molecular viscosity and ion drag force, indicating that gravity waves with a larger horizontal phase velocity (larger vertical wavelength) can penetrate into the thermosphere. The longitudinal variation of gravity wave activity in the equatorial thermosphere and upward propagation of gravity waves from the lower atmosphere were also studied. The longitudinal distribution of gravity wave activity in the thermosphere is quite similar to that of gravity wave activity in the lower atmosphere and the cumulus convective activity in the tropical troposphere. Our results indicate that the strong energy flux due to gravity waves from the enhanced cumulus convective activity propagates upward into the upper thermosphere. The relation between the wind fluctuation associated with gravity waves and the ionospheric variation is discussed. Fluctuations of the neutral zonal wind with periods of 1–2 h are significant in the 200- to 300-km height region, and its amplitude sometimes exceeds 50 m s−1. These results suggest that upward propagating gravity waves can affect the ionospheric variation in the F-region.

In the Sumatra earthquake that occurred on December 26, 2004, significant ionospheric variations were detected immediately after the earthquake in both the TEC (total electron content) data of GPS (Global Positioning System) and the ionosonde data. A magnetic pulsation with a period of about 4 min was also observed in Phimai, Thailand. Recent studies have suggested that these events are associated with acoustic waves excited by a sudden large-scale displacement of the sea surface around the epicenter. In order to study these phenomena quantitatively, a time-dependent two-dimensional nonhydrostatic compressible atmosphere-ionosphere model has been used and compared with relevant data. By modeling in sea surface perturbation, we were able to reproduce an atmospheric oscillation with a period of about 4 min in the upper atmosphere above the epicenter. The electron density variations observed by GPS/TEC and by ionosondes were also reproduced fairly well. We found that the observed TIDs (traveling ionospheric disturbances) with long periods are caused by the ducted thermospheric gravity waves produced in the thermosphere through acoustic pulse from the epicenter. The good overall agreement between the simulation results and observations indicates that numerical simulation with the nonhydrostatic compressible atmosphere-ionosphere model could be a useful tool to investigate the relationship between variations in the upper atmosphere and various sources of disturbances in the lower atmosphere.

In the Antarctic upper mesosphere and lower thermosphere around 90 km, meteor radar observations at the South Pole have detected a significant semidiurnal wind component in summer which is found to be non-migrating with zonal wavenumber s = 1. It has been surmised that this component might possibly be excited through the non-linear interaction of the migrating semidiurnal tide with stationary planetary waves with zonal wavenumber s = 1 prevailing at stratospheric heights. The Kyushu University GCM has been successful in elucidating very unambiguously this conjecture. In the present paper, linearized steady tidal modeling is carried out in this connection to reproduce, by a fairly simplified but explicit model, trans-equatorial propagation of non-migrating semidiurnal tide forced in the opposite winter hemisphere and to compare latitudinal structures of migrating and non-linearly excited intermittent tides in view of polar latitudes where non-migrating tide tends to dominate over migrating tide. It is also shown that the re-analysis meteorological data for almost 10 years clearly supports the well-known N-S asymmetry in stationary planetary wave activity in winter polar stratospheric regions, possibly due to the difference in surface topography between two hemispheres. We suggest that there might be significant asymmetry in the summertime enhancement of semidiurnal non-migrating tide between both polar regions. This phenomenon may be in a context similar with the N-S asymmetry of gravity wave activity in the polar regions, possibly giving rise to N-S disparity of cold summertime mesopause temperature. Clarification awaits intensive bi-polar studies by coordinated radar and optical observations which are running both in the Arctic and Antarctic regions.

The global traveling ionospheric disturbances (TIDs) during the drastic magnetic storms of October 29–31, 2003 were analyzed using the Global Position System (GPS) total electron content (TEC) data observed in the Asian-Australian, European and North American sectors. We collected the most comprehensive set of the TEC data from more than 900 GPS stations on the International GNSS Services (IGS) website and introduce here a strategy that combines polynomial fitting and multi-channel maximum entropy spectral analysis to obtain TID parameters. The results of our study are summarized as follows: (1) large-scale TIDs (LSTIDs) and medium-scale TIDs (MSTIDs) were detected in all three sectors after the sudden commencement (SC) of the magnetic storm, and their features showed longitudinal and latitudinal dependences. The duration of TIDs was longer at higher latitudes than at middle latitudes, with a maximum of about 16 h. The TEC variation amplitude of LSTIDs was larger in the North American sector than in the two other sectors. At the lower latitudes, the ionospheric perturbations were more complicated, and their duration and amplitude were relatively longer and larger. (2) The periods and phase speeds of TIDs were different in these three sectors. In Europe, the TIDs propagated southward; in North America and Asia, the TIDs propagated southwestward; in the near-equator region, the disturbances propagated with the azimuth (the angle of the propagation direction of the LSTIDs measured clockwise from due north with 0°) of 210° showing the influence of Coriolis force; in the Southern Hemisphere, the LSTIDs propagated conjugatedly northwestward. Both the southwestward and northeastward propagating LSTIDs are found in the equator region. These results mean that the Coriolis effect cannot be ignored for the wave propagation of LSTIDs and that the propagation direction is correlated with polar magnetic activity.

For the Western-Pacific region spread-F has been found to occur with delays after geomagnetic activity (GA) ranging from 5 to 10 days as station groups are considered from low midlatitudes to equatorial regions. The statistical (superposed-epoch) analyses also indicate that at the equator the spread-F, and therefore associated medium-scale traveling ionospheric disturbances (MS-TIDs) occur with additional delays around 16, 22 and 28 days representing a 6-day modulation of the delay period. These results are compared with similar delays, including the modulation, for D-region enhanced hydroxyl emission (Shefov, 1969). It is proposed that this similarity may be explained by MS-TIDs influencing both the F and D regions as they travel. Long delays of over 20 days are also found near the equator for airglow-measured MS-TIDs (Sobral et al., 1997). These are recorded infrequently and have equatorward motions, while normally eastward motions are measured at the equator. Also in midlatitudes D-region absorption events have been shown (statistically) to have similar long delays after GA. It is suggested that atmospheric gravity waves and associated MS-TIDs may be generated by some of the precipitations responsible for the absorption. The recording of the delayed spread-F events depends on the GA being well below the average levels around sunset on the nights of recording. This implies that lower upper-atmosphere neutral particle densities are necessary.

Southwestward-propagating medium-scale traveling ionospheric disturbances (MSTIDs) observed over Shigaraki (34.85°N, 136.10°E) in Japan on the night of May 22, 1998 are analyzed in detail. The MSTIDs were detected with a 630.0 nm (OI) all-sky imager at Shigaraki and a large number of GPS (Global Positioning System) receivers distributed around Shigaraki. Each GPS receiver provided total electron content (TEC) between the GPS altitude (20,200 km) and the ground. MSTID amplitudes varied in space and time, and showed decay and enhancement during the southwestward propagation, suggesting that amplitudes of atmospheric gravity waves and the interaction process between gravity waves and F region plasma were highly variable. It is found that spatial and temporal fluctuations of the 630 nm intensity are well correlated with those of GPS-TEC except for a certain period of time. The Scheffield University Plasmasphere Ionosphere Model (SUPIM) is used to obtain theoretical relationships between the 630 nm airglow intensity and GPS-TEC and between their fluctuation amplitudes. The results indicate that the fluctuation amplitudes observed in weak airglow regions are caused by an electron density fluctuation of about ±20% occurring around an altitude of 250 km, where the 630 nm emission rate reaches a maximum, below the F layer peak altitude. Highly enhanced 630 nm intensity and GPS-TEC within a bright airglow region are due to an electron density enhancement of about 150% occurring at altitudes below 300 km.

Three-dimensional nonlinear breaking acoustic-gravity waves (AGWs) propagating from the Earth's surface to the upper atmosphere are simulated numerically. Horizontally moving periodical structures of vertical velocity on the Earth's surface are used as AGW sources in the model. The 3D algorithm for hydrodynamic equation solution uses finite-difference analogues of basic conservation laws. This approach allows us to select physically correct generalized wave solutions of hydrodynamic equations. The numerical simulation covers altitudes from the ground up to 500 km. Vertical profiles of the mean temperature, density, molecular viscosity, and thermal conductivity are specified from standard models of the atmosphere. Atmospheric waves in a few minutes can propagate to high altitudes above 100 km after activation of the surface wave forcing. Surfaces of constant phases are quasi-vertical first, and then become inclined to the horizon below about 100 km after some transition time interval. Vertical wavelengths decrease with time and tend to theoretically predicted values after times longer than several periods of the wave forcing. Decrease in vertical wavelengths and increase in AGW amplitudes can lead to wave instabilities, accelerations of the mean flow and wave-induced jet streams at altitudes above 100 km. AGWs may transport amplitude modulation of atmospheric wave sources in horizontal directions up to very high levels. Low wave amplitudes in the beginning of transition processes after activation of atmospheric wave sources could be additional reasons for slower amplitude grows with height compared to the nondissipative exponential growth predicted for stationary linear AGWs. Production of wave-induced mean jets and their superposition with nonlinear unstable dissipative AGWs can produce strong narrow peaks of horizontal speed in the upper atmosphere. This may increase the role of transient nonstationary waves in effective energy transport and variations of atmospheric parameters and gas admixtures in a broad altitude range.

An all sky optical imaging system was operated from Sriharikota rocket range (SHAR) (14°N, 80°E, 5.5°N dip latitude) during January–March, 1993 to observe ionospheric plasma depletions through 630 nm and 777.4 nm night glow emissions. Strong plasma depletions were observed only on four nights viz., 14, 17, 19 and 21 February, 1993. Except the 17 February, which was a magnetically disturbed day, all the other nights pertained to magnetically quiet period. A number of plasma depletion parameters such as, degree of depletion, east-west extent, tilt with respect to the geomagnetic field, inter-depletion distance, drift velocity and plasma enhancements or brightness patterns were estimated. Some of the important results are: (a) It was found that the east-west extent of plasma depletions varied with the degree of depletion; for the 630 nm images the degree of depletion ranged between 6–9% per 100 km east-west extent and for 777.4 nm images it was 3% per 100 km east-west extent, (b) The average inter-depletion distance (IDD) was in the range of 1500 ± 100 km during the magnetically disturbed period and 740 ± 60 km during quiet periods. This is suggestive of gravity wave modulation of the bottom side of the F-region. While the large scale gravity waves (1500 ± 100 km) of auroral origin could be responsible during magnetically disturbed period, smaller scale gravity waves (740 ± 60 km) having their origin in the lower atmosphere could produce initial perturbation in the bottom side of the F-region, (c) Plasma depletions are observed to have an eastward tilt in the range of 10–15° with respect to the geomagnetic field. It has been suggested here that these tilts are associated with the variation of plasma drift with altitude, (d) plasma depletions are observed to be moving eastwards with drift velocities in the range of 40–190 ms−1, and (e) Strong plasma enhancements or brightness patterns were observed on three nights. The degree of enhancement was by a factor of 1.4–3.8. These enhancements lasted for more than 15 minutes. Although, prima facie, these observations look similar to the transient brightness wave reported by Mendillo et al. (1997a), the high degree of enhancement and an extended duration of more than 15 minutes, observed in the present case, need to be understood.

This paper focuses on the effect of diabatic processes due to photochemical heating on long-period gravity waves in the stratosphere, mesosphere and lower thermosphere. A linear diabatic gravity wave model is established and compared to a model of pure dynamical adiabatic gravity waves. The results indicate that the photochemistry has a damping effect on gravity waves in most regions of the stratosphere and mesosphere. However, the photochemistry has a destabilizing effect on gravity waves in the mesopause region. The photochemical heating process can induce a comparatively strong enhancement of gravity waves at the mesopause for lower temperatures. In the summer polar mesopause region, this growth rate may be greater by about one order of magnitude than the growth rate of gravity waves at other seasons and locations.

In order to clarify the horizontal structure of the wavelike oscillation frequently observed in the night time sodium density profile with a small Gaussian half-width in the middle of the night and a broad distribution at dusk, simultaneous observations with two lidars at Shigaraki (34.9°N, 136.1°E) and Hachioji (35.6°N, 139.4°E) and the MU radar at Shigaraki has been carried out. In the campaign of 35 nights, simultaneous observation was successful in four nights. On December 27–28, 1995, a large scale wave motion was observed by two lidars and the MU radar with meteor observation mode. Phase velocities were almost the same at the two sites and there were only slight differences in phase between the two sites. The wave motion was inconsistent with the component of the tidal wave and similar with the results of hodograph analyses. It is possible that the wave observed by lidars on December 27–28, 1995 was a gravity wave. The results of analyses suggest the possibility of gravity waves which were observed at fixed local time.

This paper reports first steps toward a computationally inexpensive spectral gravity wave parameterization scheme whose predictions approximate those of a full three-dimensional (in spectral space) spectral model of atmospheric gravity waves. A reduction to two dimensions, as proposed by Hines, requiring the neglect of Coriolis and non-hydrostatic effects, is explored on the basis of comparisons with a full three-dimensional power-spectral model that includes Coriolis and non-hydrostatic effects. The reduction tries to be more realistic in terms of spectral shapes, though simpler in terms of wave-breaking criteria. It works remarkably well in the absence of, but less well in the presence of, background shear. The reasons for the discrepancies are under investigation, as are the implications for two-dimensional schemes, including Hines’ as well as ours.

We present a new rocket instrument which measures total atmospheric density with great precision and resolution by Rayleigh scattering of infrared light. Comparison with a ground-based lidar shows: (a) both instruments measure the same physical parameters, even though they have different integration times and volumes. (b) The observed density structures change little over the course of an hour and the horizontal distance of 14 km. The rocket instrument has a basic vertical resolution of approximately 8 m and a number density precision of 4.2 · 1019 m−3. Above 56 km we integrate over an increasing vertical range, reaching 80 m at 70 km. The measured number density profile shows remarkable alternations between very stable layering and metastable layering (adiabatic lapse rate) in the atmosphere between 52 km and 71 km. Comparison with the hodograph of the horizontal wind profile measured by a falling sphere 29 minutes later shows that the metastable height regions coincide with height regions where the hodograph deviates from an ideal spiral. The observation is tentatively interpreted as a gravity wave that is saturating (or encountering a critical level) in these height regions. The comparison of the fine-scale neutral number density observations with measurements of ion density by electrostatic skin probes on board the same vehicle shows a number of ion density enhancements in the stably-layered height regions. With one exception out of six cases, these enhancements occur where the vertical gradient of the meridional wind would collect positive ions as in sporadic E layers. This may be the first observation of such ion density enhancements in the height region 55 to 70 km.

Generation of internal gravity waves (IGWs) by tropospheric convections and vertical propagation of the generated IGWs throughout the middle atmosphere are simulated by a non-hydrostatic compressible nonlinear two-dimensional numerical model. The present simulation demonstrates that (i) IGWs are generated by the tropospheric dry convections, (ii) zonal mean wind is decelerated by critical layer absorption of the generated IGWs in the upper tropospheric shear zone, (iii) secondary IGWs are radiated by the critical layer instability, and (iv) the secondary IGWs break down and accelerate zonal mean winds in the upper middle atmosphere. The detailed analyses show that (1) eastward propagating IGWs generated by the convection in tropospheric westerly with small wavelength of the order of 10 km and short period of the order of 10 min are dominant. It is found that the dominant waves are selected by a filtering effects of the prescribed westerly. (2) Several secondary IGWs with smaller horizontal wavelength than the primary IGW are radiated. The secondary IGWs propagate vertically in the form of wavepackets and break in the upper middle atmosphere due to local convective instability because of the exponential growth of the wave amplitudes with height. In the breaking region, the observed and theoretically predicted universal power law of the wind fluctuation, which states the m−3 dependence for the power spectra versus the vertical wavenumber m due to the wave saturation and breakdown, is also realized in the present model simulation.

Simultaneous observations of the night airglow OH (6, 2) band emission intensity and rotational temperature, by a sky scanning airglow spectrophotometer, and meteor winds, by a middle and upper atmosphere radar (MU radar), were carried out at Shigaraki (34.9°N, 136.1°E), Japan, from October 29 to November 11, 1994, as the first phase of a campaign, and from July 25 to July 31, 1995 as the second phase. Horizontal structures in the OH emission intensity and rotational temperature were monitored optically, together with the background wind and its wave induced fluctuations, measured by MU radar. Since the MU radar makes a direct measurement of the vertical wavelength, and the OH spectrophotometer makes a direct measurement of the horizontal wavelength, the two techniques are mutually complementary to determine intrinsic wave parameters. Gravity waves with intrinsic periods of 2 to 9 hours, horizontal wavelengths of 500 to 3000 km and vertical wavelengths of 12 to 75 km were identified. Between the two different observation techniques, there is a reasonable agreement in the inferred wave characteristics.

Observations of winds and gravity waves (GW) by MF radars from the Arctic to the Equator are used to provide frequency spectra and spectral variances of horizontal motions, and information on the predominant azimuthal directions of propagation for the waves. The years used are mainly 1993/4; the height layer 76–88 km; and the GW bands 10 100 min. and 1–6 hrs. The high/mid-latitude locations of Tromso, Saskatoon, London/Urbana, Yamagawa, generally demonstrate similar behaviour: the monthly spectra have slopes near −5/3 in winter months, but smaller (absolute) slopes at higher frequencies (