Inertio-gravity waves and subtropical multiple tropopauses: Vertical wavenumber spectra of wind and temperature observed by the MU radar,radiosondes and operational rawinsonde network

We have carried out continuous observations of the tropopause region over Japan for three weeks during the Baiu (early summer rain “in Japan”) season in 1991, by using a VHF Doppler radar (the MU radar), radiosondes launched at the radar site and operational rawinsondes at five meteorological stations. Based on these observations, we try to examine the hypothesis that the multiple tropopauses and the dominant inertio-gravity waves are one and the same feature, and obtain some interesting results that are not inconsistent with this hypothesis. First, vertical wavenumber spectra and hodographs analyzed from the radar wind data in the tropopause region suggest that inertio-gravity waves with vertical wavelengths of ∼ 2 km are quasi-monochromatically dominant (with 2–3 day scale variabilities of 10–20%), and are in accordance with activities of the subtropical jet stream and mesoscale cyclone-front system activities observed by the operational network. Second, striking (potential) temperature fluctuations are detected simultaneously by the radiosondes and rawinsondes, which appear as multiple tropopauses in meridional cross-section analysis. Third, vertical wavenumber spectra analyzed from the radiosonde temperature data are consistent with the radar wind spectrum, if we assume that both wind and temperature fluctuations are mainly induced by the dominant inertio-gravity waves. Finally, we confirm that the dominant interio-gravity waves can be barely detected also from routine rawinsonde (1.5-km running-mean) wind data if the amplitude is larger than 1.5 m/s. However, the monochromatic wave structures are generally quite localized in space and time.     

A study of the vertical motion field near the high-latitude summer mesopause during MAC/SINE

We present the results of high-resolution observations of the vertical velocity field obtained with the EISCAT and SOUSY VHF radars near the high-latitude summer mesopause during the MAC/SINE campaign in northern Norway in 1987. The data reveal an energetic motion field with maximum amplitudes of ~ 10 m/s and characteristic periods of ~5–30 min. Motions exhibit a high degree of vertical coherence and a quasi-periodic structure, with typical durations of 5–10 cycles. Estimates of the mean vertical velocity are downward at lower levels and are near zero or positive at greater heights. The mean vertical velocity variance is found to be ~5 m²/s², consistent with other high-latitude measurements. Frequency spectra computed for each radar are found to exhibit considerable variability, while vertical wavenumber spectra are seen to be somewhat variable in amplitude and to have slopes approaching −3 at lower wavenumbers. These results are suggestive of an energetic spectrum of gravity wave motions near the mesopause that has a large vertical flux of wave energy, that may have observed wave frequencies differing significantly from intrinsic frequencies due to Doppler shifting by large horizontal winds, and that is consistent with the separability of the frequency and wavenumber dependence of the motion spectrum and with gravity wave saturation at sufficiently small vertical scales.     

Dominant vertical scales of gravity waves in the middle atmosphere observed with the MU radar and rocketsondes

We have simultaneously observed wind motions in the altitude range of 5–90 km by means of the MU radar, rocketsondes and radiosondes. Dominant vertical scales of wind fluctuations due to gravity waves were 2–5 km in the lower stratosphere, about 5–15 km in the upper stratosphere and longer than 15 km in the mesosphere. The increase in the vertical scale with altitude is interpreted in terms of the saturation of upward propagating gravity waves. In the stratosphere, the observed vertical wavenumber spectra showed smaller amplitudes and more gradual slopes than the model values. Furthermore, the wind velocity variance in the stratosphere increases exponentially with an e-folding height of about 9 km, implying that the gravity waves were not fully saturated. On the other hand, the spectra in the upper stratosphere and mesosphere agreed fairly well with the model spectra. The variance in the mesosphere seems to cease increase of the wave amplitudes and agrees reasonably well with the model value.     

Observational evidence of local turbulence production due to gravity wave saturation in the summer mesosphere

Five foil chaff and two falling sphere rockets flown during the MAC/SINE Campaign on 15 July 1987 at Andenes, Northern Norway (69°17′N). From these rocket measurements, turbulent energy dissipation rates, vertical wind shears and Richardson numbers as functions of height were derived in the range from 82 to 92km. Turbulent energy dissipation rates generally range from 1.4 × 10⁻⁵ to 2.0 × 10⁻²W/kg and are consistent with other experiments performed at the same latitude. Strong wind shears of the order of 50–90 m/s/km are observed at various heights. Good correspondence between turbulence intensity peaks, regions of strong wind shear and low Richardson number is found. Vertical wavenumber spectra of the five scalar winds measured by the foil chaff rockets indicate that there is an excellent agreement with the saturation hypothesis, suggesting that the turbulence intensity peaks measured in this salvo are linked directly to the saturation of gravity wave motions via dynamical instabilities.     

Equatorial spread-F and neutral atmospheric turbulence : a review and a comparative anatomy

A number of satellite and rocket plasma density spectra obtained during equatorial spread-F conditions are presented and discussed in the light of similar measurements in the neutral atmosphere. We discuss this comparison in some detail and find both distinct similarities and subtle differences. The horizontal spectral measurements show a peak at an outer scale quite similar to the scale of the undulations caused by gravity wave interactions with the ionosphere. This feature is similar to a buoyancy subrange but it is easy to show that the amplitudes of the plasma fluctuations are too large to be directly driven by the neutral atmosphere. At intermediate scales the plasma fluctuations have a one-dimensional horizontal spectrum with a power law well described by a ($\text{−5}\text{3}$) slope. Once again it can be shown that the neutral fluid cannot be similarly structured at 400 km altitude due to the high viscosity coefficient. A plasma cascade process seems to be operating but it is not at all clear how the spectrum is formed. Furthermore, vertical power spectra seem to run the gamut in spectral form with slopes (n) varying in the range from −1 to −3. So the horizontal spectra are near universal in form, while the vertical spectra are quite variable.     

Variability of winds and temperatures in the lower thermosphere

High-resolution daytime incoherent scatter radar measurements of plasma temperatures and drifts in the ionospheric E-region above Millstone Hill (42.6°N, 71.5°W) have been used to derive horizontal neutral winds and temperatures in the lower thermosphere (105–130 km) during five multi-day campaigns in 1987–1991. The underlying semi-diurnal tidal component has been determined from the observations, with characteristic average amplitudes of 50 ± 15 m/s and 30 ± 10 K. Phase propagation with altitude follows the expected structure of semi-diurnal tidal modes, but reveals complex coupling of tidal modes, particularly above 115 km. Day-to-day variability in the winds and temperatures is large, and the deviations from the semi-diurnal harmonic can exceed 40 m/s and 50 K. No strong correlations have so far been found with geophysical parameters to explain the observed variability.     

Gravity wave activity in the upper mesosphere over Urbana,Illinois: lidar observations and analysis of gravity wave propagation models

We analyze 375 h of Na Wind/Temperature lidar measurements of the mesopause region (≈ 80–105 km) Na density and temperature profiles on 57 nights distributed over 2 yr at Urbana, Illinois. These observations yield a high-resolution seasonal data set of gravity wave activity in the upper mesosphere. From this data, we present measurements of the Brunt-Väisälä period, the relative atmospheric density perturbations and their spectra, and the parameters of 143 quasi-monochromatic gravity waves. The direct measurement of the Brunt-Väisälä period allows accurate calculation of the horizontal velocity perturbations and vertical displacement perturbations from the density measurements. The horizontal velocity and vertical displacement vertical wave number spectrum magnitudes and indices show considerable seasonal and nightly variability. The gravity wave amplitudes, wavelengths, and observed periods exhibit systematic relationships similar to those found in previous studies, and are consistent with the MU radar measurements of intrinsic gravity wave parameters. Here, we present a detailed analysis of the observations in terms of Diffusive-Filtering Theory models of gravity wave propagation. The magnitudes of the vertical wave number spectrum, the form of the joint vertical wave number and frequency spectrum, and the systematic relationships between the monochromatic gravity wave parameters are consistent with the Diffusive-Filtering model. We compare these results with a variety of radar, lidar, and airglow observations from other sites. This observational study suggests that the complex nonlinear interactions of the gravity wave field may be modeled successfully as a diffusive damping process, where the effective diffusivity is a function of the total wave variance.     

Characteristics of ionospheric irregularities capable of producing quasi-horizontal traces on ionograms

This paper presents simulated ionograms calculated for a parabolic ionospheric layer containing irregularities in the form of small amplitude waves. With small amplitudes, perturbation techniques can be used enabling results for the irregular ionospheres to be calculated from the results for smooth ionospheres. This approach is relatively straightforward and avoids having to ray trace new paths each time the irregularity parameters are changed. It is, however, restricted to irregularities which do not cause multiple echoes. Irregularities with vertical wavelengths of up to a few kilometres can produce significant changes in the ionosphere over height intervals smaller than those involved in reflecting a single pulse. Consequently, in the simulation procedure, it is essential to consider not just the carrier frequency but the complete frequency spectrum of the pulse. Irregularities with vertical wavelengths of the order of 10 km or more can produce ripples in an ionogram trace. These will, of course, be more evident on ionograms with high frequency resolution. Irregularities with vertical wavelengths of up to several kilometres and amplitudes up to a few per cent can produce significant pulse spreading and splitting. The actual effects depend not just on the irregularity properties but also on the ionosonde pulse width, gain and frequency and height resolutions. Some simulations show trace splitting and quasi-horizontal traces similar in many respects to effects observed by Bowman (1987, J. atmos. terr. Phys. 49, 1007) and Bowmanet al. (1988, J. atmos. terr. Phys. 50, 797). Consequently it is suggested that, at least in some cases, small amplitude (≤3%) and small scale (≤4 km) irregularities produce the spread-ifF reported by these authors.     

Estimates of gravity wave momentum fluxes in the winter and summer high mesosphere over northern Scandinavia

As part of the MAP/WINE campaign (winter 1983–1984) and the MAC/SINE campaign (summer 1987) high resolution wind profiles were obtained in the upper mesosphere using the foil cloud technique. Vertical winds were derived from the fall rate of the foil clouds and are used for estimating the momentum fluxes associated with vertical wavelengths shorter than about 10 km. From the ensemble average of 15 observations over an altitude range of 74–89 km we calculate a zonal net momentum flux of +12.6 ± 4.5 m²s⁻² in summer. The average of 14 measurements in winter between 73 and 85 km indicates a zonal net momentum flux of −3.7 ± 2.4 m²2 s⁻².     

Tidal winds from the MLT global radar network during the first LTCS campaign—September 1987

Winds and tides were measured by a number of MLT (Mesosphere, Lower Thermosphere) radars with locations varying from 43–70°N, 35–68°S, during the first LTCS (Lower Thermosphere Coupling Study) Campaign, 21–25 September 1987. The mean winds were globally westerly, consistent with early winter-like (NH) and late winter (SH) circulations.The semi-diurnal tide had vertical wavelengths near or less than 100 km at most locations, with some latitudinal variation (longer/shorter at lower latitudes in the NH/SH)—amplitudes decreased at high latitudes. The global tide was closer to anti-symmetric, with northward components being in phase at 90 km. Numerical model calculations [Forbes and Vial (1989), J. atmos. lerr. Phys. 51, 649] for September have rather similar wavelengths and amplitudes; however, the global tide was closer to symmetric, and detailed latitudinal trends differed from observed.The diurnal tide had similar wavelengths in each hemisphere, with short values (~30 km) at 35°, long (evanescence) at 68–70°, and irregular phase structures at mid-latitudes. The tide was neither symmetric nor anti-symmetric. Model calculations for the equinox [Forbes. S and Hagan (1988), Planet. Space Sci. 36, 579] were by nature symmetric, and showed the short wavelengths extending to mid-latitudes (43–52°). Southern hemisphere phases were significantly (6–8 h) different from observations. Amplitudes decreased at high latitudes in model and observation profiles.     

Tidal wind observations with incoherent scatter radar and meteorological rockets during MAP/WINE

Measurements of winds in the mesosphere and lower thermosphere were carried out during the main phase of the MAP/WINE project in January and February 1984 with the EISCAT UHF incoherent scatter radar near Tromsö, Norway, and with meteorological rockets launched from the Andøya Rocket Range, Norway. The radar measurements yield wind profiles between the altitudes of about 80 km and 105 km and the rockets between about 60 km and 90 km. Results from both techniques are combined to yield mean profiles which are particularly evaluated in terms of tidal variations. It is found that the semidiurnal tide constitutes an essential wind contribution between 85 km and 105 km. Whereas the tidal amplitudes are below 5 m s⁻¹ at about 80 km, they increase to 20–30 m s⁻¹ at 100 km. The average vertical wavelength of 35 km points to the S₄² mode, but coupling and superposition of different modes cannot be excluded.     

Reconstruction of horizontally-inhomogeneous ionospheric structure from oblique-incidence backscatter experiments

Fridman and Fridman [(1994) J. atmos. terr. Phys. 56, 115] suggested a method of reconstructing the horizontally-inhomogeneous ionospheric structure using vertical- and oblique-incidence backscatter sounding (OBS) ionograms measured at a single location. In the present paper this technique has been used to analyze experimental data and tested against independent vertical sounding (VS) measurements. By using the OBS and VS ionograms measured at Irkutsk as source data for the method we reconstructed ionization profiles over Tomsk (1050 km to the west of Irkutsk). We found that the reconstructed profiles are in reasonable agreement with the profiles obtained from VS measurements at Tomsk.     

Wave characteristics obtained from OH rotational temperatures and 557.7 nm airglow intensities

Variations of OH rotational temperature and 557.7 nm atomic oxygen intensity have been measured from Calgary, Alberta, Canada (51°10′N, 114°13′W) from 1985 to 1987. For three nights studied in detail the OH rotational temperature wave structure at 85 km was negatively correlated with the green line emission at 95 km, indicating that wave activity linked the two regions. The lower altitude region displayed high and low frequency wave structure, but by 95 km the high frequency component had disappeared. Temperature data from 16 nights during which there was obvious wave activity yielded horizontal wavelengths from about 5–100 km and inferred vertical wavelengths from 0.7 to 8 km. The horizontal and vertical angles of propagation imply a statistical source to the observed waves as being located south of Calgary along the Rocky Mountain range. There appeared to be very few, if any, wave structures propagating towards the southeast indicating a probable filtering mechanism by the background winds between the Earth's surface and 85 km.     

Rocket-borne measurements of atmospheric infrared fluxes

In the Energy Budget Campaign two rockets, one from Andøya Rocket Range, Norway, and one from Esrange, Sweden, each carrying a liquid helium cooled infrared spectrometer, were simultaneously launched as part of salvo B. The launches occurred during the recovery phase of the last of four auroral magnetic events after a Joule heating criteria was exceeded. At Andøya, zenith radiance altitude profiles were obtained of nitric oxide (NO) near 5.4μm from 70 to 185 km (rocket apogee), of ozone (O₃) near 9.6 μm from 70 to 105 km (instrument sensitivity) and of carbon dioxide (CO₂) near 15 μm from 70 to 150 km (instrument sensitivity). Measured CO₂ spectra at 72 km are shown to compare favorably to those calculated for local thermodynamic equilibrium conditions and instrument resolution. By comparing Andøya and Esrange CO₂ radiance profiles it is shown that there is evidence for spatial variation in the emission. Further, it is shown that the very disturbed conditions of salvo B prior to and during these launches appears to have significantly changed the O₃ 9.6μm radiance profiles compared to previous rocket measurements in polar disturbed and quiet atmospheres. Using the nitric oxide radiance profiles and spectrum, previous rocket results and computed models it is shown that no radiance increase could be detected from prompt auroral energy deposition. The results support the thesis that the NO density in auroral regions is significantly enhanced over mid-latitude values and that for weak auroras, the reaction NO(v = 0) + O → NO(r = 1) + O is the dominant radiation mechanism.     

Some further results on long term mesospheric and lower thermospheric wind observations by the Arecibo UHF radar

A second series of long term mesospheric and lower thermospheric wind observations was conducted at Arecibo (18.4°N, 66.8°W) between 6 and 20 March 1981 using the UHF Doppler radar, following the first observations in August 1980 (Hirota et al., 1983). Zonal and meridional wind velocities were measured during the morning (8–10 LT) and afternoon (13–15 LT) periods. The mean wind profile averaged over the entire observational period shows the predominance of the diurnal tide. The fluctuating wind vector rotates clockwise relative to height with a characteristic vertical scale of about 10 km. The phase difference inferred by a cross correlation analysis between morning and afternoon profiles indicates that the dominant period is about 20–30 h. This oscillation is discussed in relation to internal inertia-gravity waves observed by the same radar in the lower stratosphere. On the other hand, wind fluctuation with a vertical scale larger than 20 km shows a substantial day-to-day variation with a period of 5–8 days. This long period oscillation shows a good correlation with the global scale geopotential height anomalies at 1 mb (46–48 km) observed by the Tiros-N satellite at 20°N. Our evidence suggests that westward travelling planetary-scale waves with zonal wavenumber one may propagate up to the lower thermosphere.     

Electric field measurements during the MAC/EPSILON campaign

As part of the MAC/EPSILON campaign in northern Norway during October and November 1987, five rocket-borne payloads made three-axis electric field measurements of the middle atmosphere. Flights 31.066 and 31.067 consisted of large multi-experiment packages, while the other three flights (30.036, 30.037 and 30.038) were devoted primarily to electrical measurements. Simultaneous measurements of the horizontal electric field made by flights 31.066 and 30.036 were in general agreement in their limited altitude region of overlap. A simultaneous small temporal feature was observed in both datasets. The relatively more extensive horizontal E-field datasets from 31.066 and 31.067 both exhibited a decreasing mapping function with decreasing altitude, which is an indication of the observation of fields from a local auroral patch. Small-scale variations in the horizontal fields of the flights were similar to observed wave-like variations in the neutral wind field. No unusual features were observed at high altitudes in the measured E₁ field. Two of the payloads observed small vertical layer structures between 40 and 50 km. No electric field structure was observed in association with the presence of a sudden sodium layer.     

Radar and optical observations of mesospheric wave activity during the lunar eclipse of 6 July 1982

Simultaneous measurements were made using a 2.66 MHz interferometer radar, infrared photometers, and imaging systems during the total lunar eclipse of 6 July 1982. The radar data showed that a series of six discrete scatterers passed overhead at 103 km with an average spacing of 54 min, and two passed overhead at 88 km, also 54 min apart. The 88 km events were approximately 27 min out of phase with those at 103 km. One of the 88 km events was examined in detail; the radar returns appeared to come from a single scatterer or a few clustered scatterers, with a velocity of 135 m s⁻¹ almost due south, at 6° below the horizontal. The speed and period give a horizontal wavelength of 440 km, and the phase shift between 88 and 103 km activity suggests a 30 km vertical wavelength, in agreement with values for typical medium-scale traveling ionospheric disturbances (TIDs). Infrared images were made in the near infrared, and photometric measurements were made on and off the 8−3 band of OH. These observations, made from one site near the radar and a second site 575 km south, showed wavelike structures appearing first over the radar, then further south until they filled most of the sky. The speed of development of the infrared structure pattern in the sky is consistent with the 135 m s⁻¹ southward wave speed observed by the radar, but the structures themselves appeared in place, then drifted slowly northward at 10 m s⁻¹. The photographically determined wavelengths were 30–60 km, considerably shorter than the 440 km determined with the radar.     

Equatorial scintillations: advances since ISEA-6

Since the last equatorial aeronomy meeting in 1980, our understanding of the morphology of equatorial scintillations has advanced greatly due to more intensive observations at the equatorial anomaly locations in the different longitude zones. The unmistakable effect of the sunspot cycle in controlling irregularity belt width and electron concentration responsible for strong scintillation in the GHz range has been demonstrated. The fact that night-time F-region dynamics is an important factor in controlling the magnitude of scintillations has been recognized by interpreting scintillation observations in the light of realistic models of total electron content at various longitudes. A hypothesis based on the alignment of the solar terminator with the geomagnetic flux tubes as an indicator of enhanced scintillation occurrence and another based on the influence of a transequatorial thermospheric neutral wind have been postulated to describe the observed longitudinal variation.A distinct class of equatorial irregularities known as the bottomside sinusoidal (BSS) type has been identified. Unlike equatorial bubbles, these irregularities occur in very large patches, sometimes in excess of several thousand kilometers in the E-W direction and are associated with frequency spread on ionograms. Scintillations caused by such irregularities exist only in the VHF band, exhibit Fresnel oscillations in intensity spectra and are found to give rise to extremely long durations (~ several hours) of uninterrupted scintillations. These irregularities maximize during solstices, so that in the VHF range, scintillation morphology at an equatorial station is determined by considering occurrence characteristics of both bubble type and BSS type irregularities.The temporal structure of scintillations in relation to the in situ measurements of irregularity spatial structure within equatorial bubbles has been critically examined. A two-component irregularity spectrum with a shallow slope (p₁ ~ 1.5) at long scalelengths (> 1km) and steep slope (p₂ ~−3) at shorter scalelengths has been found in both vertical and horizontal spectra. Phase and intensity scintillation modelling was found to be consistent with this two-component irregularity spectrum.Finally, the information provided by the major experimental undertaking represented by Project Condor in the fields of night-time scintillations and zonal irregularity drifts with be briefly outlined.     

Role of vertical winds on the Rayleigh-Taylor mode instabilities of the night-time equatorial ionosphere

In view of the recent observations on the presence of vertical winds in the equatorial ionosphere in the evening and night-time, the role of vertical winds in the Rayleigh-Taylor (R-T) mode instability has been re-examined. The mathematical treatment of Chiu and Straus, earlier developd for a case of horizontal winds, is extended to evaluate the role of vertical winds in causing the R-T mode instability. It is shown that the vertical (downward) winds of small magnitude have a very significant effect on the instability growth rate in the. F-region. A downward wind of l m s⁻¹ can cause the same growth rate as a 200 m s⁻¹ eastward wind at 260 km altitude. Furthermore, a downward wind of 16m s⁻¹ at 300 km can be as effective as that due to the gravitational drift itself. Similarly, an upward wind can inhibit the instability on the bottomside of the F-region. It appears that the polarity of the vertical winds (upward or downward) at the base of the F-layer plays an important role in the growth of the R-T mode plasma instability in the equatorial ionosphere.     

Calculation of auroral Bahner volume emission height profiles in the upper atmosphere

Energetic protons entering the atmosphere will either travel as auroral protons or as neutral hydrogen atoms due to charge-exchange and excitation interactions with atmospheric constituents. Our objective is to develop a simple procedure to evaluate the Balmer excitation rates of H_α and H_β, and produce the corresponding volume emission rates vs height, using semi-empirical range relations in air, starting from proton spectra observed from rockets above the main collision region as measured by Reasoneret al. [(1968) J. geophys. Res.73, 4185] and Søbraaset al. [(1974) J. geophys. Res.79, 1851]. The main assumptions are that the geomagnetic field is parallel and vertical, and that the pitch angle of the proton/hydrogen atom is preserved in collisions with atmospheric constituents before being thermalized. Calculations show that the largest energy losses occur in the height interval between 100 and 125 km, and the corresponding volume emission rate vs height profiles have maximum values in this height interval. The calculted volume emission rate height profile of H_β compares favorably with that measured with a rocket-borne photometer.