A diagnostic algorithm, based on the empirical normal mode decomposition technique, is proposed as a diagnostic tool in studies of the atmospheric variability. It begins by analyzing the transient eddies in terms of empirical modes that are orthogonal with respect to wave activities. Time-dependent amplitudes together with wave activity spectra are used to classify the modes and compute their propagation properties. The algorithm is applied to a sequence of four Northern Hemisphere winters taken from the National Centers for Environmental Prediction reanalyses, with a focus on the upper troposphere and lower stratosphere, giving a set of empirical modes of wind, pressure, specific volume, and potential vorticity. Results indicate that most of the wave activity is carried by large-scale, eastward-propagating modes centered at middle and high latitudes. Some properties of the leading modes, such as their average phase speeds, are in good agreement with the predictions of linear dynamics. Characteristics of the leading wavenumber-5 mode, such as its dipolar pressure pattern near the summer hemisphere tropopause, its propagation speed of 12 m s−1 and decay rate of 3 days, can be explained by the theory of quasi modes, defined as superpositions of singular modes sharply peaked in the phase speed domain. Other large-scale, midlatitude modes also show properties compatible with the quasi-modal description, suggesting that quasi modes play an important role in the upper-troposphere dynamics.