Search Results

You are looking at 11 - 20 of 61 items for

  • Author or Editor: Michael Ghil x
  • Refine by Access: All Content x
Clear All Modify Search
Cécile Penland and Michael Ghil

Abstract

Multivariate linear prediction based on single-lag inverse modeling is developed further and critically examined. The method is applied to the National Meteorological Center analyses of Northern Hemisphere 700-mb geopotential height anomalies, which have been filtered to eliminate periods shorter than 10 days. Empirically derived normal modes of the randomly forced linear system are usually correlated, even at zero lag, suggesting that combinations of modes should be used in predictions. Due to nonlinearities in the dynamics and the neglect of interactions with other pressure levels, the lag at which the analysis is performed is crucial; best predictions obtain when the autocovariances involved in the analysis are calculated at a lag comparable to the exponential decay times of the modes. Errors in prediction have a significant seasonal dependence, indicating that the annual cycle affects the higher-order statistics of the field. Optimized linear predictions using this method are useful for about half a day longer than predictions made by persistence.

Conditional probabilities are much more efficiently calculated using normal-mode parameters than from histograms, and yield similar results. Maps of the model's Fourier spectra—integrated over specified frequency intervals and consistent with the assumptions made in a linear analysis—agree with maps obtained from fast Fourier transforms of the data.

Full access
Masahide Kimoto and Michael Ghil

Abstract

This paper presents an observational analysis of recurrent flow patterns in the Northern Hemisphere (NH) winter, based on a 37-year series of daily 700-mb height anomalies. Large-scale anomaly patterns that appear repeatedly and persist beyond synoptic time scales are identified by searching for local maxima of probability density in a phase subspace, which is spanned by the leading empirical orthogenal functions (EOFs).

By using an angular probability density function (PDF), we focus on the shape, not magnitude, of the anomaly patterns. The PDF estimate is nonparametric; that is, our algorithm makes no a priori assumption on symmetry with respect to the climatological mean as in one-point correlation and rotated EOF analyses. The local density maxima are searched by iterative bump hunting.

Based on observed partial decoupling between the Pacific (PAC) and the Atlantic-Eurasian (ATL) sectors, the classification algorithm is applied separately to each of the two. Seven PAC and six ATL patterns are obtained. Anomaly maps that belong to the neighborhood of each PDF peak are associated with distinct flow regimes. These include regional blocked and zonal flows, and wave train-like anomaly patterns, some of them well known from previous studies, others revealed by our analysis for the first time.

Successive appearances of flow regimes are generally separated by unclassifiable, transient periods. A Markov chain describes transitions between different flow regimes; highly likely, as well as unlikely routes of transition exist. Chains of preferred transitions may be related to the existence of oscillatory modes in the NH extratropics.

A synoptic characterization of onsets and breaks for the flow regimes obtained is given by compositing. In situ evolutions of anomaly patterns, slow westward shifts of high-latitude anomaly centers, and successive down-stream increase of anomaly magnitudes are the typical signatures of such events.

Full access
Nathan Paldor and Michael Ghil

Abstract

The linear instability of a zonal geostrophic jet with a cosh−2 meridional profile on an f plane is investigated in a reduced-gravity, shallow-water model. The stability theory developed here extends classic quasigeostrophic theory to cases where the change of active-layer depth across the jet is not necessarily small. A shooting method is used to integrate the equations describing the cross-stream structure of the alongstream wave perturbations. The phase speeds of these waves are determined by the boundary conditions of regularity at infinity. Regions exist in parameter space where the waves that propagate along the jet will grow exponentially with time. The wavelength of the most unstable waves is 2π R, where R is the internal deformation radius on the deep side, and their e-folding time is about 25 days.

The upper-layer thickness of the basic state in the system has a spatial structure resembling that of the isopycnals across the Gulf Stream. The unstable waves obtained in the present analysis have a wavelength that is in agreement with some recent observations—based on infrared imaging of the sea surface temperature field—of the fastest- growing meanders’ wavelength. Calculated growth rates fall toward the low end of the range of values obtained from these infrared observations on the temporal evolution of Gulf Stream meanders.

Full access
Fei Chen and Michael Ghil

Abstract

An idealized North Atlantic Ocean model is forced by climatological wind stress, restoring temperature, and a diagnosed salinity flux. Both centennial and interdecadal oscillations are sustained in the model if the diagnosed salinity flux is characterized by net evaporation in high latitudes. To investigate further the role of salinity fluxes two different linear profiles are imposed: one has net evaporation in high latitudes and the other net precipitation. The first salinity flux induces a purely interdecadal oscillation in the model, while the second one causes a millennial and a decadal-to-interdecadal oscillation. Next, the authors consider a boundary condition for temperature expressed as the sum of a fixed heat flux and a restoring term. Constant heat flux characterized by net cooling in high latitudes leads to an interdecadal oscillation similar to the one caused by net evaporation.

Both the decadal-to-interdecadal and the purely interdecadal oscillation are upper-ocean phenomena. Inter-decadal anomalies are mainly confined to high latitudes, with their center moving anticlockwise near the north-west corner of the model domain; they are amplified and sink in that region. Decadal-to-interdecadal anomalies are mainly confined to midlatitudes, advected eastward by the mean flow, and disappear near the cast coast.

The physical mechanisms for the two oscillations are different. The interdecadal oscillation is caused by surface-density variations in northern high latitudes; variations are due to either net evaporation from the applied salinity flux or constant cooling from the applied heat flux. The decadal-to-interdecadal oscillation is a by-product of deep-water warming, due to the strong braking effect of salinity forcing on thermal forcing: surface saline water from the subtropics overlies continuously warming intermediate water to provide a favorable environment for the decadal-to-interdecadal oscillation. Further analysis implies that in a fully coupled ocean-atmosphere situation the decadal-to-interdecadal oscillation is less likely to exist.

Full access
Michael Ghil and Kingtse Mo

Abstract

We have examined systematically oscillatory modes in the Northern Hemisphere and in the tropics. The 700 mb heights were used to analyze extratropical oscillations, and the outgoing longwave radiation to study tropical oscillations in convection. All datasets were band-pass filtered to focus on the intraseasonal (IS) band of 10–120 days. Leading spatial patterns of variability were obtained by applying EOF analysis to these IS data. The leading principal components (PCs) were subjected to singular spectrum analysis (SSA). SSA is a statistical technique related to EOF analysis, but in the time domain, rather than the spatial domain. It helps identify nonlinear oscillations in short and noisy time series.

In the Northern Hemisphere, there are two important modes of oscillation with periods near 48 and 23 days, respectively. The 48-day mode is the most important of the two. It has both traveling and standing components, and is dominated by a zonal wavenumber two. The 23-day mode has the spatial structure and propagation properties described by Branstator and by Kushnir.

In the tropics, the 40–50 day oscillation documented by Madden and Julian, Weickmann, Lau, their colleagues, and many other authors dominates the Indian and Pacific oceans from 60°E to the date line. From 170°W to 90°W, however, a 24–28 day oscillation is equally strong. The extratropical modes are often independent of, and sometimes lead, the tropical modes.

Full access
Andreas Groth and Michael Ghil

Abstract

Singular spectrum analysis (SSA) along with its multivariate extension (M-SSA) provides an efficient way to identify weak oscillatory behavior in high-dimensional data. To prevent the misinterpretation of stochastic fluctuations in short time series as oscillations, Monte Carlo (MC)–type hypothesis tests provide objective criteria for the statistical significance of the oscillatory behavior. Procrustes target rotation is introduced here as a key method for refining previously available MC tests. The proposed modification helps reduce the risk of type-I errors, and it is shown to improve the test’s discriminating power. The reliability of the proposed methodology is examined in an idealized setting for a cluster of harmonic oscillators immersed in red noise. Furthermore, the common method of data compression into a few leading principal components, prior to M-SSA, is reexamined, and its possibly negative effects are discussed. Finally, the generalized Procrustes test is applied to the analysis of interannual variability in the North Atlantic’s sea surface temperature and sea level pressure fields. The results of this analysis provide further evidence for shared mechanisms of variability between the Gulf Stream and the North Atlantic Oscillation in the interannual frequency band.

Full access
Masahide Kimoto and Michael Ghil

Abstract

Recurrent and persistent flow patterns are identified by examining multivariate probability density functions (PDFs) in the phase space of large-scale atmospheric motions. This idea is pursued systematically here in the hope of clarifying the extent to which intraseasonal variability can be described and understood in terms of multiple flow regimes.

Bivariate PDFs of the Northern Hemisphere (NH) wintertime anomaly heights at 700 mb are examined in the present paper, using a 37-year dataset. The two-dimensional phase plane is defined by the two leading empirical orthogonal functions (EOFs) of the anomaly fields. PDFs on this plane exhibit synoptically intriguing and statistically significant inhomogeneities on the periphery of the distribution. It is shown that these inhomogeneities are due to the existence of persistent and recurrent anomaly patterns, well-known as dominant teleconnection patterns; that is, the Pacific/North American (PNA) pattern, its reverse, and zonal and blocked phases of the North Atlantic Oscillation (NAO). It is argued that the inhomogeneities are obscured when PDFs are examined in a smaller-dimensional subspace than dynamically desired.

Full access
Ricardo Todling and Michael Ghil

Abstract

Sequential data assimilation schemes approaching true optimality for sizable atmospheric models are becoming a reality. The behavior of the Kalman filter (KF) under difficult conditions needs therefore to be understood. In this two-part paper we implement a KF for a two-dimensional shallow-water model, with one or two layers. The model is linearized about a basic flow that depends on latitude; this permits the one-layer (1-L) case to be barotropically unstable. Constant vertical shear in the two-layer (2-L) case induces baroclinic instability.

A model-error covariance matrix for the KF simulations is constructed based on the hypothesis that an ensemble of slow modes dominates the errors. In the 1-L case, the system is stable for a meridionally constant basic flow. Assuming equipartition of energy in the construction of the model-error covariance matrix has a deleterious effect on the process of data assimilation in both the stable and unstable cases. Estimation errors are found to be smaller for a model-error spectrum that decays exponentially with wavenumber than an equipartition spectrum. Then the model-error covariance matrix for the 2-L model is also obtained using a decaying-energy spectrum.

The barotropically unstable 1-L case is studied for a basic velocity profile that has a cosine-square shape. Given this linear instability, forecast errors grow exponentially when no observations are present. The KF keeps the errors bounded, even when very few observations are available. The best placement of a single observation is determined in this simple situation and shown to be where the instability is strongest. The 2-L case and a comparison with the performance of a currently operational data assimilation scheme will appear in Part II.

Full access
Hisanori Itoh and Michael Ghil

Abstract

Numerical experiments are performed to clarify the excitation mechanism of mixed Rossby-gravity waves (Yanai waves) in the tropical troposphere, as well as the selection of zonal wavenumbers 4–5 and of the five-day period. The model used is governed by the primitive equations on an equatorial β-plane. Moisture budgets are calculated explicitly.

A nonlinear wave-CISK mechanism produces Yanai waves with the same spectral peaks in wavenumber and frequency as observed. In the absence of antisymmetric lateral forcing, these peaks do not appear distinctly, because the symmetric equatorially trapped modes, i.e., Kelvin-like waves having different spectral peaks, are dominant. It is the lateral antisymmetric forcing which puts the peaks characterizing the antisymmetric Yanai waves in evidence.

It appears that Yanai waves of very small wavenumbers (1–3) cannot have large amplitudes because their frequencies are too large for moisture to be effectively supplied for the convection associated with these waves. Symmetric Kelvin modes are dominant in the absence of forcing asymmetries due at least in part to the difference in the nature of heating between symmetric and antisymmetric modes: precipitation, and hence heating, is not normally distributed. Given a strongly skewed distribution of heating, it can be shown that symmetric modes are excited more effectively. Finally, our results indicate that the vertical wavenumber, and hence the period of Yanai waves are selected by the height of cumulus convection, while the lateral forcing selects the horizontal wavenumber within a certain band provided by the nonlinear wave-CISK mechanism.

Full access
Kingtse C. Mo and Michael Ghil

Abstract

Persistent anomalies with recurrent spatial patterns play an important role in the atmosphere's low-frequency variability. We establish a connection between statistical and dynamical methods of description and prediction of persistent anomalies. This is done 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. Model solutions exhibit persistent anomalies identifiable with blocked, zonal and wave-train anomalies in Northern Hemisphere atmospheric data. Flow structures similar to the patterns above occur as high-variance EOFs of this nonlinear model.

The Southern Hemisphere data we analyze 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.

We conclude 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, we introduce a Markov-chain formulation of transitions from one persistent anomaly regime to another, and discuss the implications for long-range forecasting.

Full access