scispace - formally typeset
Search or ask a question

Showing papers by "Goddard Space Flight Center published in 1987"


Journal ArticleDOI
TL;DR: In this article, Orthogonally rotated principle component analysis (RPCA) was used to identify and describe the seasonality and persistence of the major modes of interannual variability.
Abstract: Orthogonally rotated principle component analysis (RPCA) of Northern Hemisphere 1-month mean 700 mb heights is used to identify and describe the seasonality and persistence of the major modes of interannual variability. The analysis is detailed and comprehensive, in that 1) a high resolution, approximately equal-area 358-point grid is used for the virtually maximum possible 35-year period of record, 2) a positive bias in the NMC data base in the early 1950s in the subtropics is largely eliminated for the first time, and 3) homogeneous, separate analyses of each month of the year are carried out, detailing the mouth-to-month changes in the dominant circulation patterns. Winter results are similar to those of other recent RPCA and teleconnection studies except that some less obvious patterns are identified and further detail of the better-known patterns is provided. Two north-south dipole patterns are found over the Pacific Ocean (West Pacific Oscillation and East Pacific pattern) and over the Atla...

3,330 citations


Journal ArticleDOI
TL;DR: The MSIS-86 empirical model of thermospheric temperature, density and composition as discussed by the authors uses new temperature and composition data from the Dynamics Explorer satellite to improve the representation of polar region morphology over that in theMSIS-83 model.
Abstract: The MSIS-86 empirical model of thermospheric temperature, density and composition uses new temperature and composition data from the Dynamics Explorer satellite to improve the representation of polar region morphology over that in the MSIS-83 model. Terms were added or changed to better represent seasonal variations in the polar regions under both quiet and magnetically disturbed conditions. Local time variations in the magnetic activity effect were added. In addition a new species, atomic nitrogen, was added to the previous list of N2, O2, He, O, H, and Ar covered by the model.

1,699 citations


Journal ArticleDOI
TL;DR: In this paper, large-scale revisions of the OGO 6 dawn-dusk measurement models are made, showing that the deformations of the two-cell patterns lead to sunward convection in dayside polar regions, while maintaining the integrity of the night-side convection pattern.
Abstract: The present analysis of electric field measurements from the Dynamics Explorer 2 satellite, which extends previous empirical models, emcompasses much data from polar crossings entering and exiting the high latitudes in all magnetic local time zones. The goal is to represent the typical distributions of convective electric fields with a minimum number of characteristic patterns. Significant large-scale revisions of the OGO 6 dawn-dusk measurement models are made. The deformations of the two-cell patterns lead to sunward convection in dayside polar regions, while maintaining the integrity of the nightside convection pattern.

1,060 citations


Journal ArticleDOI
TL;DR: In this article, the authors compared the accuracy of the three products, namely, SMMR, NOAA/NESDIS and USAFGWC, and concluded that the results of the comparisons are good for uniform snow covered areas, such as the Canadian high plains and the Russian steppes.
Abstract: Snow covers about 40 million km2 of the land area of the Northern Hemisphere during the winter season. The accumulation and depletion of snow is dynamically coupled with global hydrological and climatological processes. Snow covered area and snow water equivalent are two essential measurements. Snow cover maps are produced routinely by the National Environmental Satellite Data and Information Service of the National Oceanic and Atmospheric Administration (NOAA/NESDIS) and by the US Air Force Global Weather Center (USAFGWC). The snow covered area reported by these two groups sometimes differs by several million km2, Preliminary analysis is performed to evaluate the accuracy of these products. Microwave radiation penetrating through clouds and snowpacks could provide depth and water equivalent information about snow fields. Based on theoretical calculations, snow covered area and snow water equivalent retrieval algorithms have been developed. Snow cover maps for the Northern Hemisphere have been derived from Nimbus-7 SMMR data for a period of six years (1978–1984). Intercomparisons of SMMR, NOAA/NESDIS and USAFGWC snow maps have been conducted to evaluate and assess the accuracy of SMMR derived snow maps. The total snow covered area derived from SMMR is usually about 10% less than the other two products. This is because passive microwave sensors cannot detect shallow, dry snow which is less than 5 cm in depth. The major geographic regions in which the differences among these three products are the greatest are in central Asia and western China. Future study is required to determine the absolute accuracy of each product. Preliminary snow water equivalent maps have also been produced. Comparisons are made between retrieved snow water equivalent over large area and available snow depth measurements. The results of the comparisons are good for uniform snow covered areas, such as the Canadian high plains and the Russian steppes. Heavily forested and mountainous areas tend to mask out the microwave snow signatures and thus comparisons with measured water equivalent are poorer in those areas.

650 citations


Journal ArticleDOI
TL;DR: A theory of the origin of intraseasonal oscillations of the tropical atmosphere is presented and tested by simple model experiments in this paper, where the authors validate the basic theory against key features of the observed 40-50 day oscillation.
Abstract: A theory of the origin of intraseasonal oscillations of the tropical atmosphere is presented and tested by simple model experiments. This study forces on the validation of the basic theory against key features of the observed 40–50 day oscillation. It is shown that the observed eastward propagation of intraseasonal oscillation in the tropical atmosphere arises as an intrinsic mode of oscillation resulting from an interaction of convection and dynamics via the so-called “mobile” wave-CISK mechanism. Through this mechanism, the heat source feeds on the east-west asymmetry of forced equatorial waves. As a result, Kelvin waves are selectively amplified, which in turn causes the heat source to propagate eastward. This mechanism also prevents small-scale waves from immediate destabilization, contrary to the results of traditional wave-CISK theory. The “mobile” wave-CISK establishes a new dynamics equilibrium state between convection and the wind field to form a wave packet or collective motion with rel...

463 citations


Journal ArticleDOI
TL;DR: In this article, several currently available broadband parameterizations for longwave and shortwave radiation have been combined to produce a computationally fast radiation parameterization that is well suited for atmospheric circulation models.
Abstract: Several currently available broadband parameterizations for longwave and shortwave radiation have been combined to produce a computationally fast radiation parameterization that is well suited for atmospheric circulation models. The main features of the parameterization are the ability to include overlapping partly cloudy layers in the longwave, the use of a delta-Eddington technique to treat clouds in the shortwave, and a computational structure that is amenable to vectorization on supercomputers. Selected results of off-line one-dimensional computations using the code have been compared with more rigorous methods as part of an international intercomparison program and found to be quite accurate.

461 citations


Journal ArticleDOI
TL;DR: In this article, Hour-averaged data from the Helios and Voyager spacecraft are used to investigate the origin and evolution of low-frequency interplanetary fluctuations from 0.3 to 20 AU.
Abstract: Hour-averaged data from the Helios and Voyager spacecraft are used here to investigate the origin and evolution of low-frequency interplanetary fluctuations from 0.3 to 20 AU. The previously observed evolution toward a less purely Alfvenic state with increasing heliocentric distance is shown to occur more rapidly in the inner heliosphere and in low-speed as well as high-speed streams. It is concluded that outward-traveling flucutations are predominantly generated by the sun, but that in situ turbulence, most likely due to stream shear, generates fluctuations with both inward and outward senses of correlation.

393 citations


Journal ArticleDOI
TL;DR: In this article, the satellite data were combined with field data on soil respiration and a global map of net primary productivity to obtain the seasonal exchange of CO2 between the atmosphere and the terrestrial biosphere.
Abstract: Normalized difference vegetation indices derived from radiances measured by the Advanced Very High Resolution Radiometer were used to prescribe the phasing of terrestrial photosynthesis. The satellite data were combined with field data on soil respiration and a global map of net primary productivity to obtain the seasonal exchange of CO2 between the atmosphere and the terrestrial biosphere. The monthly fluxes of CO2 thus obtained were employed as source/sink functions in a global three-dimensional atmospheric tracer transport model to simulate the annual oscillations of CO2 in the atmosphere. The results demonstrate that satellite data of high spatial and temporal resolution can be used to provide quantitative information about seasonal and longer-term variations of photosynthetic activity on a global scale.

373 citations


Journal ArticleDOI
TL;DR: In this article, a two-stream approximation to the radiative transfer equation is used to calculate the vegetation indices (simple ratio and normalized difference), the fraction of incident photosynthetically active radiation (PAR) absorbed by the canopy, and the daily mean canopy net photosynthesis under clear-sky conditions.

350 citations


Journal ArticleDOI
TL;DR: In this paper, the properties of interplanetary magnetic field and plasma data acquired by Voyager between 1 and 11 AU are used to investigate the properties and evolution of inter-planetary MHD fluctuations and to answer several related questions regarding the Alfvenicity of solar wind fluctuations.
Abstract: The magnetic field and plasma data acquired by Voyager between 1 and 11 AU are used to investigate the properties of interplanetary MHD fluctuations and to attempt to answer several related questions regarding the Alfvenicity of solar wind fluctuations. These questions are: (1) the extent to which the interplanetary fluctuations are Alfvenic; (2) whether there is any evolution in propagation direction of the Alfvenic fluctuations; (3) whether the degree and evolution of Alfvenicity are correlated with structure; and (4) the importance and evolution of compressive contributions to the fluctuations. The conclusions on these points are summarized, and the results are related to theoretical work.

308 citations


Journal ArticleDOI
TL;DR: In this paper, a net-radiation-based empirical model for soil heat flux (G) is analyzed for inclusion in a canopy-temperature-based energy balance equation to estimate evaporation from a growing wheat crop.

Journal ArticleDOI
TL;DR: In this article, data from ISEE 3, Helios A, and Helios B were used to identify the components of two compound streams and to determine their configurations, each of which produced an unusually large geomagnetic storm, on April 3, 1979 and on April 25, 1979, respectively.
Abstract: Data from ISEE 3, Helios A, and Helios B were used to identify the components of two compound streams and to determine their configurations. In one case, ejecta containing a magnetic cloud associated with a disappearing quiescent filament were interacting with a corotating stream. In the second case, ejecta containing a magnetic cloud associated with a 2B flare were overtaking ejecta from a different source. Each of these compound streams produced an unusually large geomagnetic storm, on April 3, 1979, and on April 25, 1979, respectively. The largest geomagnetic storm in the period 1968-1986, which occurred on July 13, 1982, was associated with a compound stream. Thirty geomagnetic storms with A(p) greater than 90 occurred between 1972 and 1983, and there are interplanetary magnetic field and plasma data for 17 of these events. The data suggest that most large geomagnetic storms are associated with compound streams and/or magnetic clouds.

Journal ArticleDOI
TL;DR: In this article, the effect of spectral emissivity on the measurement of land surface (LS) temperature from satellite radiances using the split-window method is investigated, and formulas are derived to relate the LS temperatures to two brightness temperatures measured from space in the AVHRR 4 and AVHRr 5 channels, and their accuracies are discussed.
Abstract: The effect of emissivity on the measurement of land surface (LS) temperature from satellite radiances using the split-window method is investigated. Formulas are derived to relate the LS temperatures to two brightness temperatures measured from space in the AVHRR 4 and AVHRR 5 channels, and their accuracies are discussed. Results indicate that to determine LS temperatures from a satellite, the spectral emissivity must be known to an accuracy on the order of 0.005 for the average and 0.0007 for the difference in order to obtain an error of the order of 0.5 K.

Journal ArticleDOI
TL;DR: In this article, a connection between statistical and dynamical methods of description and prediction of persistent anomalies is established by computing and analyzing the empirical orthogonal functions (EOFs) in a simple deterministic model, and in Southern Hemisphere geopotential heights, on the other.
Abstract: Persistent anomalies with recurrent spatial patterns play an important role in the atmosphere's low-frequency variability. A connection between statistical and dynamical methods of description and prediction of persistent anomalies is established by computing and analyzing the empirical orthogonal functions (EOFs) in a simple deterministic model, on the one hand, and in Southern Hemisphere geopotential heights, on the other. The dynamical model is governed by the fully nonlinear, equivalent-barotropic vorticity equation on the sphere, with simplified forcing, dissipation and topography. The Southern Hemisphere data consist in gridded daily maps of 500 mb heights from June 1972 to July 1983. Two types of persistent anomalies appear in this time series, both having a strong wavenumber-three component; they differ by the value of the constant phase of this wave and by the strength of the wavenumber-one component. The first two EOFs bear a striking resemblance to these two patterns. It is concluded that the dynamical interpretation of EOFs is their pointing from the time mean to the most populated regions of the system's phase space. Pursuing this interpretation, a Markov-chain formulation of transitions from one persistent anomaly regime to another is introduced, and the implications for long-range forecasting are discussed.

Journal ArticleDOI
TL;DR: In this article, the solar exoatmospheric irradiances for the TM bands 1, 2, 3, and 4 have been calculated based on the earlier published data of Neckel and Labs.
Abstract: Based on solar irradiance data published by Neckel and Labs (1984) and Iqbal (1983), the solar exoatmospheric irradiances for Thematic Mapper (TM) bands 1, 2, 3, and 4 have been calculated. Results vary by up to 1 percent from previous published values, which were based on the earlier data of Neckel and Labs. For TM bands 5 and 7, integrated solar exoatmospheric irradiances have also been recalculated using solar irradiance data published by Labs and Neckel (1968), Arvesen et al. (1969), and Iqbal (1983). These irradiances vary by up to 6 percent from previously published results, which were based on data published by Thekaekara (1972).

Journal ArticleDOI
TL;DR: In this paper, three main approaches to land cover classification are considered, namely, principal components transformed images, the use of a characteristic curves procedure based on NDVI values plotted against time, and finally application of the maximum likelihood rule to multitemporal data sets.
Abstract: Various methods are compared for carrying out land cover classifications of South America using multitemporal Advanced Very High Resolution Radiometer data. Fifty-two images of the normalized difference vegetation index (NDVI) from a 1-year period are used to generate multitemporal data sets. Three main approaches to land cover classification are considered, namely the use of the principal components transformed images, the use of a characteristic curves procedure based on NDVI values plotted against time, and finally application of the maximum likelihood rule to multitemporal data sets. Comparison of results from training sites indicates that the last approach yields the most accurate results. Despite the reliance on training site figures for performance assessment, the results are nevertheless extremely encouraging, with accuracies for several cover types exceeding 90 per cent.

Journal ArticleDOI
TL;DR: In this paper, the authors use multitemporal satellite data to estimate the photosynthetic capacity of the terrestrial surface and record microwave surface brightness at the 37 GHz frequency for detecting and quantifying areas experiencing drought.

Journal ArticleDOI
TL;DR: In this article, a model of the spatial and temporal distribution of rainfall is described that produces random spatial rainfall patterns with these characteristics: (1) the model is defined on a grid with each grid point representing the average rain rate over the surrounding grid box, (2) rain occurs at any one grid point, on average, a specified percentage of the time and has a lognormal probability distribution, and (3) spatial correlation of the rainfall can be arbitrarily prescribed, and time stepping is carried out so that large scale features persist longer than small-scale features.
Abstract: A model of the spatial and temporal distribution of rainfall is described that produces random spatial rainfall patterns with these characteristics: (1) the model is defined on a grid with each grid point representing the average rain rate over the surrounding grid box, (2) rain occurs at any one grid point, on average, a specified percentage of the time and has a lognormal probability distribution, (3) spatial correlation of the rainfall can be arbitrarily prescribed, and (4) time stepping is carried out so that large-scale features persist longer than small-scale features. Rain is generated in the model from the portion of a correlated Gaussian random field that exceeds a threshold. The portion of the field above the threshold is rescaled to have a lognormal probability distribution. Sample output of the model designed to mimic radar observations of rainfall during the Global Atmospheric Research Program Atlantic Tropical Experiment (GATE), is shown. The model is intended for use in evaluating sampling strategies for satellite remote-sensing of rainfall and for development of algorithms for converting radiant intensity received by an instrument from its field of view into rainfall amount.

Journal ArticleDOI
27 Mar 1987-Science
TL;DR: The correlation between the satellite-derived green vegetation index and the ecological parameters associated with RVF virus suggested that satellite data may become a forecasting tool for RVF in Kenya and, perhaps, in other areas of sub-Saharan Africa.
Abstract: Data from the advanced very high resolution radiometer on board the National Oceanic and Atmospheric Administration's polar-orbiting meteorological satellites have been used to infer ecological parameters associated with Rift Valley fever (RVF) viral activity in Kenya. An indicator of potential viral activity was produced from satellite data for two different ecological regions in Kenya, where RVF is enzootic. The correlation between the satellite-derived green vegetation index and the ecological parameters associated with RVF virus suggested that satellite data may become a forecasting tool for RVF in Kenya and, perhaps, in other areas of sub-Saharan Africa.

Journal ArticleDOI
TL;DR: For all interplanetary (IP) shocks which generated so-called IP type II events, the associated solar events involved fast coronal mass ejections (CMEs).
Abstract: The interplanetary shocks which generate detectable low-frequency radio emission, represent as a group, the most energetic shocks produced by the sun. For all interplanetary (IP) shocks which generated so-called IP type II events, the associated solar events involved fast coronal mass ejections (CMEs). In comparison with the set of all CMEs detected by the Solwind coronagraph, the CMEs associated with IP type II events are the most massive and energetic. The majority belong to the structural classes described by the Solwind researchers as 'curved front' or 'halo'.

Journal ArticleDOI
TL;DR: In this article, the relationship between the evapotranspiration and soil moisture during the drying, supply-limited phase is studied, and the effects of natural soil heterogeneities on evapOTranspiration computed from the model are investigated.
Abstract: The relationship between the evapotranspiration and soil moisture during the drying, supply-limited phase is studied. A second scaling parameter, based on the evapotranspirational supply and demand concept of Federer (1982), is defined; the parameter, referred to as the threshold evapotranspiration, occurs in vegetation-covered surfaces just before leaf stomata close and when surface tension restricts moisture release from bare soil pores. A simple model for evapotranspiration is proposed. The effects of natural soil heterogeneities on evapotranspiration computed from the model are investigated. It is observed that the natural variability in soil moisture, caused by the heterogeneities, alters the relationship between regional evapotranspiration and the area average soil moisture.

Journal ArticleDOI
TL;DR: In this article, a spherically symmetric earth model and centroid-moment tensor solutions for earthquake sources were used to compute the earthquake-induced changes in the earth's rotation and low-degree harmonics of the gravitational field for the period 1977-1985.
Abstract: Analytical formulas based on the normal-mode theory are used together with a spherically symmetric earth model and the centroid-moment tensor solutions for earthquake sources to compute the earthquake-induced changes in the earth's rotation and low-degree harmonics of the gravitational field for the period 1977-1985. Spectral and statistical analyses are conducted on these changes. It is found that the earthquake-induced changes are two orders of magnitude smaller than those observed; most of these changes show strong evidence of nonrandomness either in their polarity or in their directions.

Journal ArticleDOI
TL;DR: In this article, the wave equation for a dissipative, compressible plasma is derived from the linearized magnetohydrodynamic equations for a plasma with transverse Alfven speed gradients.
Abstract: An improved method for calculating the resonance absorption heating rate is discussed and the results are compared with observations in the solar corona. To accomplish this, the wave equation for a dissipative, compressible plasma is derived from the linearized magnetohydrodynamic equations for a plasma with transverse Alfven speed gradients. For parameters representative of the solar corona, it is found that a two-scale description of the wave motion is appropriate. The large-scale motion, which can be approximated as nearly ideal, has a scale which is on the order of the width of the loop. The small-scale wave, however, has a transverse scale much smaller than the width of the loop, with a width of about 0.3-250 km, and is highly dissipative. These two wave motions are coupled in a narrow resonance region in the loop where the global wave frequency equals the local Alfven wave frequency. Formally, this coupling comes about from using the method of matched asymptotic expansions to match the inner and outer (small and large scale) solutions. The resultant heating rate can be calculated from either of these solutions. A formula derived using the outer (ideal) solution is presented, and shown to be consistent with observations of heatingmore » and line broadening in the solar corona. 34 references.« less

Journal ArticleDOI
TL;DR: In this article, a target calibration procedure for obtaining surface albedo from satellite data is presented, which addresses two key issues, the calibration of remotely-sensed, discrete wavelength, digital data and the derivation of an albedos measurement (defined over the solar short wave spectrum) from spectrally limited observations.
Abstract: A target calibration procedure for obtaining surface albedo from satellite data is presented. The methodology addresses two key issues, the calibration of remotely-sensed, discrete wavelength, digital data and the derivation of an albedo measurement (defined over the solar short wave spectrum) from spectrally limited observations. Twenty-seven Landsat observations, calibrated with urban targets (building roof-tops and parking lots), are used to derive spatial and seasonal patterns of surface reflectance and albedo for four land cover types: city, suburb, farm and forest.

Journal ArticleDOI
TL;DR: In this paper, the Compton-getting interpretation of the dipole anisotropy in the microwave background is presented as a diagnostic of the weak asymmetry exhibited by the cosmic X-ray background.

Journal ArticleDOI
TL;DR: In this article, an ecological model is developed to estimate annual net primary productivity (NPP) in 12 North American biomes, combining existing models which address canopy photosynthesis in response to light, temperature, and moisture availability, and account for respiration.

Journal ArticleDOI
TL;DR: In this article, the theory of gamma-ray and neutron production is discussed and results of calculations are compared to gamma-rays, neutron, and charged-particle observations from solar flares.
Abstract: Nuclear processes and particle acceleration in solar flares are reviewed. The theory of gamma-ray and neutron production is discussed and results of calculations are compared to gamma-ray, neutron, and charged-particle observations from solar flares. The implications of these comparisons on particle energy spectra, total numbers, anisotropies, electron-to-proton ratios, as well as on acceleration mechanisms and the interaction site, are presented. The information on elemental and isotopic abundances derived from gamma-ray observations is compared to abundances obtained from escaping accelerated particles and other sources.

Journal ArticleDOI
TL;DR: In this article, a computer simulation and theoretical study of the physical characteristics of the explosive coalescence of current-carrying loops is presented, which includes a large impulsive increase of the kinetic energies of electrons and ions, the simultaneous heating and acceleration of electrons in high and low energy ranges, and a break in the energy spectra.
Abstract: A computer simulation and theoretical study of the physical characteristics of the explosive coalescence of current-carrying loops is presented. Characteristics of the explosive coalescence include a large impulsive increase of the kinetic energies of electrons and ions, the simultaneous heating and acceleration of electrons and ions in high and low energy ranges, and a break in the energy spectra of electrons and ions. A characteristic double subpeak structure is found in the quasi-periodic oscillations found in the time profiles of the solar flares of June 7, 1980 and November 26, 1982 which can be explained in terms of the coalescence instability of two current loops.

Journal ArticleDOI
TL;DR: In this paper, the interaction of transverse vertical circulations associated with two separate jet steak/trough systems is found to be a common feature of cyclogenetic events which produce heavy snow along the East Coast of the United States.
Abstract: The interaction of transverse vertical circulations associated with two separate jet steak/trough systems is found to be a common feature of cyclogenetic events which produce heavy snow along the East Coast of the United States. The transverse circulations are identified for eight cases that span the period from 1960 to 1987 utilizing an isentropic analysis of the operational radiosonde data. The analyses depict the interaction of 1) a direct circulation located within the confluent entrance region of an upper-level jet streak over the northeastern United States or southeastern Canada with 2) an indirect circulation in the diffluent exit region of a jet streak associated with a trough nearing the East Coast. This interaction contributes to differential moisture and temperature advections and vertical motions necessary to produce heavy snowfall along the coast. It is suggested that the circulation patterns associated with the jet streak establish an environment within which boundary layer processe...

Journal ArticleDOI
01 Oct 1987-Tellus A
TL;DR: In this article, the authors study the forecast error growth in the 100-day ECMWF data set of 10-day forecasts previously utilized by Lorenz, separating the square of the error into systematic and random components.
Abstract: We study the forecast error growth in the 100-day ECMWF data set of 10-day forecasts previously utilized by Lorenz, separating the square of the error into systematic and random components. The random error variance is approximately equally distributed among all zonal wavenumbers m corresponding to the same total wavenumber n. Therefore, we combine the random error variance for all zonal wavenumbers and study its dependence on total (2-dimensional) wavenumber rather than on zonal wavenumber. We extend the Lorenz model for global r.m.s. error growth by including the effect of errors associated with deficiencies in the forecast model , and apply this new parameterization to the error variance of 10-day ECMWF operational forecasts obtaining an excellent fit. We point out that the commonly used parameter “doubling time of small errors” is not a good measure of error growth because it has to be determined by extrapolation to small errors and it is very sensitive to the method of extrapolation. On the other hand, the error growth at finite times is a better measure because it is well defined by the data. The results of the error fit to each total wavenumber n are the basis for the main conclusions of the paper. In the northern hemisphere winter data set, the error growth rateα increases almost monotonically with wavenumber, from about 0.3 day -1 at long waves to about 0.45 day -1 at medium and short waves. The saturation error variance V∞ is about 30% larger than the error variance at day 10 for long waves, which have not yet reached saturation. The scaled source of external error S / V ∞ (due to model deficiencies) varies from about 3% day -1 at long waves to about 20% day -1 at short waves. This increase of relative model error may be due to the second-order finite differences used in the ECMWF model during 1980/1981. Finally, we estimate the theoretical limit of dynamical predictability as a function of total wavenumber n . We define it as the time at which the error variance reaches 95% of the saturation value. This time decreases monotonically with n. In the winter, waves n ≪ 10 have not reached saturation at 10 days. In a perfect model with error growth similar to the present model, the predictability limit would be extended by about a week. In the summer, error growth rates and model deficiencies are larger than in winter, and the limit of predictability is close to 10 days even for long waves. DOI: 10.1111/j.1600-0870.1987.tb00322.x