MST radar studies of wind and turbulence in the middle atmosphere

A survey is presented of recent developments in the observation of wind and turbulence in the stratosphere and mesosphere using MST radars. One of the highlights of these developments is the growing recognition that the MST/ST radar is a valuable tool for routine monitoring of the atmospheric wind field. Furthermore, preliminary observations have shown the feasibility of monitoring atmospheric turbulence as well. Recent observations of mesospheric turbulence support theoretical models that emphasize the role of propagating waves in coupling the lower and middle atmospheres. Scientific groups in several countries are now planning or constructing MST radars so that within a few years observations should be available from diverse geographical locations spanning the globe.     

MF Doppler and spaced antenna radar measurements of upper middle atmosphere winds

Observations of upper mesospheric and lower thermospheric wind velocities obtained simultaneously over six days with MF Doppler and Spaced Antenna (SA) radars at Adelaide, Australia in November 1980 are presented. To obtain these measurements, the large (~ 1 km diameter) Buckland Park MF array was run in a dual beam Doppler radar configuration, and a portable SA radar was operated adjacent to the main array. Hourly mean values of wind velocity show considerable consistency, with cross correlation coefficients of about 0.6–0.8 for the entire observational period. However, agreement between the magnitudes of the wind velocities as measured by each technique is found to be significantly improved when the effect of the aspect sensitivity of the backscattering irregularities on the effective beam pointing angle of the Doppler radar beams is taken into account. This is also found to be true for SA and Doppler radar observations obtained in adjacent periods of 2–5 days over two years with the Buckland Park facility operating alternatively as a Doppler and SA radar. Some representative examples of these results are also presented and discussed. A preliminary comparison between MF Doppler and SA radar derived vertical wind velocities is also briefly considered.     

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.     

VHF radar observations of nonlinear interactions in the summer polar mesosphere

VHF radar measurements of velocities and echo power in the summer polar mesosphere have been analysed using maximum entropy, bispectral and cross-spectral methods in order to study wave-wave interactions. The results show nonlinear interactions of second and even third order between the diurnal and semidiurnal tides and planetary waves with periods of 2 and 3 days in the velocity field. Similar analyses of time series of echo power suggest corresponding variations of the temperature field.     

A study of the VHF radar data via a random variance model

Data sets from some VHF radars have been analysed. Gaussian distributions with random variance are proposed for the signal's quadrature components. The suggested distributions explain the data sets satisfactorily, especially as the length of the data series increases. Non-stationarity of the signals will also be interpreted using the proposed model. Moreover, a χ²-goodness of fit test for the proposed model has been conducted and its results are persuasive. We suggest that it is better to use the proposed distribution for the quadrature components than to use the Nakagami distribution for the amplitude distribution or the regular Gaussian distribution for the quadrature components. In addition, the sampling time should be less than 4 min to guarantee the stationarity of the data.     

Some aspects of Doppler radar measurements of the mean and fluctuating components of the wind field in the upper middle atmosphere

The basic assumptions made when a Doppler radar is used to measure the mean and fluctuating components of the wind field in the middle atmosphere with various beam configurations are examined. Particular reference is made to the measurement of the various components of the Reynolds stress tensor associated with short period internal gravity waves. It is shown that it is not generally possible to measure the upward flux of horizontal momentum with the conventional Doppler radar beam configuration in the upper middle atmosphere and that an optimum beam configuration is that in which beams are directed at +θ,0 and − θ to the zenith in both the zonal and meridional planes. This allows five of the six components of the Reynolds stress tensor (all those except the horizontal transport of momentum) to be obtained directly from the mean square radial velocities. In addition, the mean wind components and, in principle, the horizontal divergence and stretching deformations may be obtained. The power spectrum of the horizontal velocity may also be calculated using only the assumption that the statistics of the motions are horizontally homogeneous.     

Foil chaff clouds as a tool for in-situ measurements of atmospheric motions in the middle atmosphere: their flight behaviour and implications for radar tracking

The aerodynamic behaviour of foil chaff (rectangular plates) used for in-situ measurements of atmospheric motions in the middle atmosphere (up to 100 km altitude) can be described by Stokes' Law in which the corrections are applied to the coefficient of viscosity. The results obtained with this approximation are in good agreement with observations and allow us to explain in detail certain peculiarities occasionally seen in the tracking of chaff clouds by radar.     

Climatology of gravity waves in the middle atmosphere

An attempt is made at the statistical analysis of small-scale disturbances in the stratosphere and mesosphere with the aid of meteorological rocket observations at many stations from 77°N to 8°S for several years.By applying a high-pass filter to daily rocket data in the height range 20–65 km, wind and temperature fluctuations with characteristic vertical scales close to or less than 10 km are obtained, which are considered to be due to internal gravity waves. Results are expressed in terms of parameters which tend to emphasize smallscale vertical fluctuations and which should provide qualitative measures of gravity wave activity.It is found that the gravity wave activity shows a notable annual cycle in higher latitudes with the maximum in wintertime, while it shows a semiannual cycle in lower latitudes with the maxima around equinoxes. It is also found from the standard deviation around the monthly mean that the temporal variability of gravity waves is very large.     

Climatology of the middle atmosphere of the Southern hemisphere

Main features of spatial distribution and time variations of meteorological parameters in the Southern hemisphere at altitudes 25–80 km are reviewed on the basis of zonal empirical models revised in 1982. Meridional distribution and seasonal variations are analysed. For comparison purposes with the Northern hemisphere, a model developed by Cole and Kantor in 1978 is used. It is revealed that the compilation of new models of the Southern hemisphere atmosphere has not resulted in substantial revision of hemispheric-structure outlined earlier in studies conducted in the Central Aerological Observatory. Meridional distribution of the parameters in summer is characterized by higher values of temperature, pressure and density gradients in the stratosphere of the Southern hemisphere than in that of the Northern hemisphere. This resulted in greater advancement of the core of the summer-time easterly (low towards the equator in the Southern hemisphere than in its northern counterpart. In winter, meridional temperature gradients in the middle stratosphere are greater in the Southern hemisphere than those in the Northern hemisphere, however, rapid attenuation of the gradients with height is observed in the Southern hemisphere, and above 35–40 km they become negative near 50–60°S, in contrast to positive values typical for the Northern hemisphere stratosphere. In the wind field, specific features of the Southern hemisphere westerly flow are high intensity and relatively low altitude of the maximum speed (as compared to the Northern hemisphere).The phases of the annual temperature wave at low latitudes are similar south and north of the equator; south of 30°S a reversal of the phase is observed. The semi-annual oscillation of temperature and wind is less pronounced in middle and high latitudes of the Southern hemisphere than in the Northern hemisphere.The origin of the primary differences between the hemispheres is related mainly to lower intensity of large-scale stratospheric processes in the Southern hemisphere as compared to those in the Northern hemisphere, which is confirmed by values of the standard deviation of the parameters in the two hemispheres.In summer, temperature and pressure fields based on satellite data are symmetric relative to the poles in both hemispheres. In winter, the distortion of the mean pressure field in the mesosphere may be as great in the Southern as in the Northern hemisphere.     

High resolution turbulence observations in the middle and lower atmosphere by the MU radar with fast beam steerability: preliminary results

Refractive index fluctuations or turbulence in the mesosphere, stratosphere and troposphere are observed with the aid of the fast beam steerability of the MU (middle and upper atmosphere) radar which operates at 46.5 MHz with 1 MW peak radiation power and 8330 m² antenna aperture. Morphology of the mesospheric and stratospheric turbulence is studied by making use of the high altitude and time resolutions. Sixteen beam observations based on the fast beam steerability reveal advection properties and spatial variability of echoing regions in the troposphere. These results demonstrate new possibilities for this system in the investigation of three dimensional structures of turbulence.     

First results with the Adelaide VHF radar: spaced antenna studies of tropospheric winds

The first results from a VHF radar of the ST type located at Buckland Park near Adelaide, Australia (35°S, 138°E), are presented. The radar is designed to be versatile and can be used to measure velocities in the lower atmosphere using both the spaced antenna (SA) and Doppler beam-swinging (DBS) techniques. Here studies of irregularities and motions made with the spaced antenna technique are discussed. It is shown that the scale of the diffraction pattern formed by the backscattered radiation varies with altitude, with the mean pattern scale being smaller in the troposphere than in the stratosphere. The observations are consistent with the backscattered energy decreasing as a function of off-vertical angle by 1.5 dB per degree in the troposphere and by about 2.8 dB per degree in the lower stratosphere. An intercomparison of zonal velocities measured with the SA and DBS methods shows good agreement. In May and August 1984 an extensive comparison was made between the velocities measured by the SA method and winds determined from over 80 balloon-borne radiosondes released from Adelaide Airport, situated some 36 km to the south of the radar. The velocities were compared on a statistical basis and showed excellent agreement, although the SA speeds tended to be 1–2 m s⁻¹ smaller in magnitude than the radiosonde velocities. Overall, the rms differences between the two sets of measurements was only 3–4ms⁻¹ throughout the troposphere, a result which is consistent with the random errors inherent in each technique, as well as the spatial separation between the radar and balloon observations. The utility of the SA method for meteorological observations is illustrated by a study of both the horizontal and vertical wind fields during the passage of a cold front made in November 1984. The high time resolution available with the radar allows detailed studies of the development of the pre-frontal jet, the wind convergence into the front and associated vertical motions.     

The characteristics of VHF radar echoes from the tropopause region near a jet stream

Measurements carried out with the SOUSY radar system to the south-east of a jet stream during a night in November 1981 are described. Particular attention is paid to the association between power and spectral width of echoes returned from the vicinity of the tropopause located near 12 km. A negative association is found between power and spectral width of echoes from 150 m height channels between 10.65 km and 11.70 km, both with the radar beam directed vertically and at 7° to the vertical in the northerly and easterly directions. Evidence for a positive association is only found for height channels centred at 11.85 and 12.15 km and then chiefly with the beam directed vertically. A partial correlation analysis suggests that the possible influences of the finite beamwidth of the radar system and of wind shear on the spectral width do not affect the associations found. The results are discussed in terms of the scattering and partial reflection mechanisms giving rise to the echoes.     

Electrodynamic response of the middle atmosphere to auroral pulsations

On the nights of 21 and 28 October 1987, two Nike Orion payloads (NASA 31.066 and 31.067) were launched from Andøya, Norway, as part of the MAC/EPSILON campaign, to study the effect of auroral energetics on the middle atmosphere. Each payload carried detectors to measure relativistic electrons from 0.1 to 1.0MeV in 12 differential energy channels, and bremsstrahlung X-rays from >5 to >80keV in 5 integral channels. In addition, instrumentation to measure bulk ion properties and electric fields was also carried by these and/or near simultaneous flights. Flight 31.066 was launched during the recovery phase of a moderate magnetic substorm, during relatively stable auroral conditions. Flight 31.067 was launched during highly active post-break-up conditions during which Pc 5 pulsations (> 150s period) were in progress. The energetic radiation of the first event was composed almost entirely of relativistic electrons below 200 keV with negligible contributions from bremsstrahlung X-rays, while the radiation of the second event was dominated by much softer electrons ( < 100 kcV), which produced high X-ray fluxes that exceeded the cosmic ray background as an ionizing source down to altitudes below 30 km. Simultaneous conductivity measurements during both events show consistency with the ionizing radiations, with the pulsation event producing free electrons down to 55 km. far below their expected altitude range during night-time. These comparisons are discussed to evaluate the impact of such events on the middle atmosphere.     

Possible effects of solar variability on the middle atmosphere

This paper discusses the possible effects on the middle atmosphere of

• 1.(i) solar variability in emission in the u.v. and far u.v.,
• 2.(ii) solar proton events,
• 3.(iii) relativistic electron precipitation events,
• 4.(iv) corpuscular heating in auroras,
• 5.(v) joule heating by the auroral electrojet,
• 6.(vi) auroral NO production and
• 7.(vii) gravity wave emission by the auroral electrojet. It is important to establish the maximum depth in the atmosphere to which these effects go, in order to be able to determine whether the reflection of planetary scale waves from the troposphere is substantially altered by them.

     

Absolute reflectivities and aspect sensitivities of VHF radio wave scatterers measured with the SOUSY radar

By accurately calibrating the SOUSY radar in West Germany it has been made possible to measure absolute values of effective reflection coefficients and turbulence structure constants. Some typical values of these parameters as a function of altitude are presented. Such profiles are presented for both a vertically directed beam, and also for two beams directed 7° off-vertical. Comparisons of powers on the vertical and off-vertical beams show that scatter became more aspect sensitive at the tropopause and in the lower stratosphere, but, unexpectedly, scatter was observed to become considerably more isotropic in the higher regions of the stratosphere (above 15–18 km) on this occasion. An enhancement of signal from the tropopause occurred not only on the vertical beam, but also on the off-vertical beams.Comparisons of signal strengths scattered from the mesosphere and measured with the vertical and off-vertical beams showed that for the present observations mesospheric scatter was close to isotropic. The backscatter cross-sections at VHF have been compared with other measurements at medium and high frequencies at other locations, and these comparisons help set some limits on the scales of turbulent and specular scatterers in the mesosphere.     

Fractal analysis of gravity wave spectra in the middle atmosphere

Wind fluctuations in the middle atmosphere behave like colored noise processes. They have a continuum of scales without dominant features and a power spectrum density (PSD) that often decays with frequency ƒ as ƒ^−β. Spectral index β is generally obtained through least-square fit to PSD estimated by Fourier methods. Graphs of colored noise have fractal plane-filling properties depending on β. An efficient method for finding β using the fractal dimension (D), based on analysis of 1/ƒ noise in galactic X-ray luminosities by McHardy I. and Czerny B., (1987, Nature325, 696), is described. An empirical relation is found between D and β and its validity is confirmed in limiting cases. Then D is obtained from power-law dependence of a length metric L(μ) on scale μ. The method is applied to middle-atmospheric velocity data from the Poker Flat radar in Alaska. Variations of D follow those in β, from an earlier analysis by Bemraet al., (1986, Handbook for MAP20, 216), but show an offset of 0.1–0.2 even after corrections for outliers, gaps, and additive noise. Usefulness of this method for screening data as an aid to spectral analysis is examined.     

A search for a solar cycle-temperature relationship in the middle atmosphere over Volgograd

The variation of temperature in the middle atmosphere (15–80 km) at Volgograd (49°N, 44°E) during an 11-year solar cycle (1971–1982) has been studied. The temperature of the stratosphere did not show any significant influence of the sunspot cycle, but the temperatures of the mesosphere showed a strong in-phase relationship with the solar cycle. Computed correlation and regression coefficients were positive and highly significant in this region. At 60 and 70km the temperature variations were almost linearly related to the sunspot number. Seasonal studies indicated that solar activity has a much stronger influence on temperature during the winter than during the summer.     

Simulation of gravity wave effects under solstice conditions using a 3-D circulation model of the middle atmosphere

Vertically propagating gravity waves can transport momentum and energy from the troposphere up to the mesosphere and thus modify the circulation of the middle atmosphere. The effects of regional gravity wave sources, together with temporal changes of gravity wave activity, are studied under solstice conditions in a 3-D circulation model using a simplified parameterization scheme for the gravity momentum deposition. In this way we can reproduce the reversal of the mean zonal wind with height and very low temperatures at the summer mesopause region. Using a stochastic forcing by taking the gravity wave parameters at random, characteristic oscillations are found with periods in the planetary scale range (2, 4 and 5 days) and in the tidal range (1 day, 16 h and 12 h).     

A simple model of gravity-wave momentum and energy fluxes transferred through the middle atmosphere to the upper atmosphere

Vertical fluxes of momentum and energy through the middle atmosphere are calculated by using a simple semi-empirical model of quasi-monochromatic internal gravity waves with dominant vertical wavenumbers. In this model those dominant gravity waves are assumed to saturate and break at each observational altitude by an effective critical-layer mechanism. The dominant value of the vertical wave-number is expressed by an exponential function of altitude, decreasing upward with a scale height of 34 km. This expression gives the momentum and energy flux densities decreasing upward with scale heights of 12 and 18 km, respectively, and typical values at 100 km altitude are estimated as 4 × 10⁻⁵ Pa and 4 × 10⁻³ W/m². A heat flux induced by wavebreaking turbulence also has an order of magnitude similar to that of the wave energy flux. Variabilities around these values and comparisons with other momentum and heat inputs to the upper atmosphere are only briefly discussed.     

A diurnal study of the electrical structure of the equatorial middle atmosphere

Electrical parameters measured from 115 km down to below 20 km during the Project Condor campaign at the Punta Lobos Rocket Range near Lima, Peru, are presented Ten rocket-launched payloads measured electrical conductivity. A strong diurnal influence due to solar ultraviolet radiation is shown. Nine of the payloads also measured electric fields. No large mesospheric vertical electric fields are found in the data. A calculation of the d.c. global conduction current density at 18 km is smaller than previously measured at low latitudes and does not show the conventional diumal variation.