A statistical procedure is developed to obtain the parameters of spectral harmonics of internal gravity waves (IGW) with horizontal wavelengths larger than ~ 60 km using measurements of wind velocity in the middle atmosphere from the middle and upper atmosphere (MU) radar. The method is applied to find sets of IGW parameters for the winter and summer seasons of 1987-1988 for 6-km-thick layers centered at 70, 75, and 80 km at Shigaraki, Japan (34.9øN, 139.4øE). The numbers of selected IGW harmonics having periods from 5 min up to 6 hours vary between 360 and 2000 depending on the season and altitude. Mean vertical wave numbers are 0.6-0.9 km -• in winter and 0.5-1 km -• in summer, and horizontal wave numbers are (3-4)x10 -2 km -• and (2-4)x10 -2 km -•, respectively. Predominant azimuths of wave vectors and momentum fluxes concentrate near 640 in winter and 900 in summer at all altitudes. In general, the difference between the azimuths of the wave vector and momentum flux is close to zero or 1800 . Turbulent and instrumental noise can contribute 20-30% of the observed intensity of the wind variations at altitudes of 70-85 km. The IGW momentum fluxes are 0.2-2.7 m•s -: in winter and 0.2-1.9 m:s -• in summer. The wave drag of the mean flow is 0.5-45 m2s -• in winter and 2-25 m2s -• in summer. In most cases the momentum flux and wave drag are directed oppositely to the mean wind. 1. 1976; Reid and Vincent, 1987; Gardner and Voe!z, 1987; Gavrilov and Kalov, 1987; Kashcheyev and Oleynikov, 1994], directions of IGW propagation [Frills and Vin-Unfortunately, most experiments are not able to measure all IGW parameters simultaneously. Also, IGW characteristics are usually obtained for a relatively broad range of frequencies or for cases when one IGW harmonic dominates over others, and visual inspection reveals clear wavelike structures or quasi-elliptical wind hodographs.In this paper we reevaluate the methods used to determine IGW characteristics and modify them to obtain parameters for each spectral harmonic of the complicated noisy wave spectrum present in the atmosphere. Formulae are obtained that allow the calculation of the parameters of IGW harmonics using measurements of wind velocity in the middle atmosphere with the MU radar in five directions of the antenna beam (one vertical and four inclined to the north, west, south, and east from zenith). The statistical method includes spectral analysis of the data for the five directions of the antenna beam, calculation of the axes and the tilt angle of the velocity hodograph ellipse, amplitudes and momentum flux, vertical and horizontal wave numbers, and azimuths of wave vector and momentum flux for each spectral harmonic. Spectra of wind variations in the middle and upper atmosphere are very complicated, and may include (in addition to IGWs) harmonics caused by turbulence, trapped waves, impacts of solar events, etc. Spectral analysis of short data series is notorious for giving "computational" peaks of spectral density that do not exist in nature (see, for example, Gavri...