Search NASA⌕ Search

Engineering topics

Ghil, Michael

Publications and source records attributed to Ghil, Michael.

17 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Advanced data assimilation in strongly nonlinear dynamical systems

Advanced data assimilation methods are applied to simple but highly nonlinear problems. The dynamical systems studied here are the stochastically forced double well and the Lorenz model. In both systems, linear approximation of the dynamics about the critical points near which regime transitions occur is not always sufficient to track their occurrence or nonoccurrence. Straightforward application of the extended Kalman filter yields mixed results. The ability of the extended Kalman filter to track transitions of the double-well system from one stable critical point to the other depends on the frequency and accuracy of the observations relative to the mean-square amplitude of the stochastic forcing. The ability of the filter to track the chaotic trajectories of the Lorenz model is limited to short times, as is the ability of strong-constraint variational methods. Examples are given to illustrate the difficulties involved, and qualitative explanations for these difficulties are provided. Three generalizations of the extended Kalman filter are described. The first is based on inspection of the innovation sequence, that is, the successive differences between observations and forecasts; it works very well for the double-well problem. The second, an extension to fourth-order moments, yields excellent results for the Lorenz model but will be unwieldy when applied to models with high-dimensional state spaces. A third, more practical method--based on an empirical statistical model derived from a Monte Carlo simulation--is formulated, and shown to work very well. Weak-constraint methods can be made to perform satisfactorily in the context of these simple models, but such methods do not seem to generalize easily to practical models of the atmosphere and ocean. In particular, it is shown that the equations derived in the weak variational formulation are difficult to solve conveniently for large systems.

Miller, Robert N.↗
NASA NTRS1993

Intraseasonal variability in a barotropic model with seasonal forcing

It has recently been suggested that oscillatory topographic instability could contribute to low-frequency variability over the Northern Hemisphere midlatitudes. A barotropic potential vorticity model, with a hierarchy of forcing and topography configurations on the sphere, is used to investigate the nature of low-frequency oscillations induced by such instabilities. Steady-state solutions of the model include multiple unstable equilibria that sustain oscillatory instabilities with periods of 10 - 15 days, and 150 - 180 days, for a realistic forcing pattern, Time-dependent solutions exhibit chaotic behavior with episodic oscillations, featuring both the intraseasonal (35 - 50 day) and biweekly (10 - 15 day) modes. The former is dominated by standing spatial patterns, the latter by traveling wave patterns. The phases of the intraseasonal oscillation are robust for all cases, exhibiting a clear oscillatory exchange of atmospheric angular momentum with the solid earth via mountain troque. It is demonstrated, through linear stability analysis on the sphere, that the intraseasonal oscillations are induced by topographic instabilities. The role of the seasonal cycle is studied by prescribing an annual cycle in the forcing. In this case, the winter forcing is more favorable than the summer for the occurrence of episodic intraseasonal oscillations. Recent observations are consistent with this model result.

Strong, Christopher M.↗
NASA NTRS1993

Multiple flow regimes in the Northern Hemisphere winter. I - Methodology and hemispheric regimes. II - Sectorial regimes and preferred transitions

Multivariate pdfs in the phase space of large-scale atmospheric motions are examined in order to identify recurrent and persistent flow patterns. The extent to which intraseasonal variablity can be described and understood in terms of multiple flow regimes is clarified. It is shown that the synoptically intriguing and statistically significant inhomogeneities exhibited by the 2D phase plane on the periphery of the distribution are due to the existence of persistent and recurrent anomaly patterns. It is argued that these inhomgeneities are obscured when PDFs are examined in a smaller-dimensional subspace than dynamically desired. A synoptic characterization of onsets and breaks for the flow regimes is given by compositing. In situ evolutions of anomaly patterns, slow westward shifts of high-latitude anomaly centers, and successive downstream increase of anomaly magnitudes are the typical signatures of such events.

Kimoto, Masahide↗
NASA NTRS1992

Adaptive filtering and prediction of the Southern Oscillation index

Singular spectrum analysis (SSA), a variant of principal component analysis, is applied to a time series of the Southern Oscillation index (SOI). The analysis filters out variability unrelated to the Southern Oscillation and separates the high-frequency, 2- to 3-year variability, including the quasi-biennial oscillation, from the lower-frequency 4- to 6-year El Nino cycle. The maximum entropy method (MEM) is applied to forecasting the prefiltered SOI. Prediction based on MEM-associated autoregresive models has useful skill for 30-36 months. A 1993-1994 La Nina event is predicted based on data through February 1992.

Keppenne, Christian L.↗
NASA NTRS1992

Aspects Of 40- to 50-Day Oscillations In LOD And AAM

Report presents study of fluctuations in rotation of Earth, focusing on irregular intraseasonal oscillations in length of day (LOD) and atmospheric angular momentum (AAM) with periods varying from 40 to 50 days. Study draws upon and extends results of prior research.

Dickey, Jean O.↗
NASA NTRS1992

Extreme weather events

Models that try to forecast the detailed geographical distribution of oceanic and atmospheric variables affected by the ENSO cycle are briefly discussed. Combinations of singular-spectrum analysis and the maximum entropy method that hold promise for predicting the ENSO cycle 2-3 yrs in advance are addressed.

Keppenne, Christian L.↗
NASA NTRS1991

Adaptive filtering and maximum entropy spectra with application to changes in atmospheric angular momentum

The spectral resolution and statistical significance of a harmonic analysis obtained by low-order MEM can be improved by subjecting the data to an adaptive filter. This adaptive filter consists of projecting the data onto the leading temporal empirical orthogonal functions obtained from singular spectrum analysis (SSA). The combined SSA-MEM method is applied both to a synthetic time series and a time series of AAM data. The procedure is very effective when the background noise is white and less so when the background noise is red. The latter case obtains in the AAM data. Nevertheless, reliable evidence for intraseasonal and interannual oscillations in AAM is detected. The interannual periods include a quasi-biennial one and an LF one, of 5 years, both related to the El Nino/Southern Oscillation. In the intraseasonal band, separate oscillations of about 48.5 and 51 days are ascertained.

Penland, Cecile↗
NASA NTRS1991

Intraseasonal oscillations in the global atmosphere. I - Northern Hemisphere and tropics

Oscillatory modes in the Northern Hemisphere and in the tropics were examined systematically. 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 empirical orthogonal function analysis to these IS data. The leading principal components were subjected to singular spectrum analysis. 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 (1987). In the tropics, the 40-50 day oscillation documented by Madden and Julian (1972), Weickmann (1983), Lau, and their colleagues, dominates the Indian and Pacific oceans from 60 deg E to the date line. From 170 deg W to 90 deg W, however, a 24-28 day oscillation is equally strong. The extratropical modes are often independent of, and sometimes lead, the tropical modes.

Ghil, Michael↗
NASA NTRS1991

Intraseasonal oscillations in the global atmosphere. II - Southern Hemisphere

A systematic examination of oscillatory modes in the global atmosphere is completed by studying 12 years of 500 mb geopotential heights in the Southern Hemisphere. The data were band-pass filtered to focus on intraseasonal (IS) phenomena, and spatial empirical orthogonal functions (EOFs) were obtained. The leading principal components were subjected to singular spectrum analysis in order to identify nonlinear IS oscillations with high statistical confidence. In the Southern Hemisphere, the dominant mode has a period of 23 days, with spatial patterns carried by the second and third winter EOF of the IS band. It has a zonal wavenumber-four structure. The 40-day mode is second, and is dominated by wavenumbers three and four, while a 16-day mode is too weak to separate its spatial behavior from the previous two. The IS dynamics in the Southern Hemisphere is more complex and dominated by shorter wavenumbers than the Northern Hemisphere. No statistically significant correlations between the Southern Hemisphere and the tropics or the Northern Hemisphere are apparent in the IS band.

Ghil, Michael↗
NASA NTRS1990

Statistical significance test for transition matrices of atmospheric Markov chains

Low-frequency variability of large-scale atmospheric dynamics can be represented schematically by a Markov chain of multiple flow regimes. This Markov chain contains useful information for the long-range forecaster, provided that the statistical significance of the associated transition matrix can be reliably tested. Monte Carlo simulation yields a very reliable significance test for the elements of this matrix. The results of this test agree with previously used empirical formulae when each cluster of maps identified as a distinct flow regime is sufficiently large and when they all contain a comparable number of maps. Monte Carlo simulation provides a more reliable way to test the statistical significance of transitions to and from small clusters. It can determine the most likely transitions, as well as the most unlikely ones, with a prescribed level of statistical significance.

Vautard, Robert↗
NASA NTRS1989

Empirical orthogonal functions and multiple flow regimes in the Southern Hemisphere winter

A new empirical orthogonal function (EOF) analysis of winter 500 mb geopotential height anomalies in the Southern Hemisphere is presented. An earlier EOF analysis prefiltered the anomalies to exclude wavenumbers 5 and higher; the present analysis does not. The different preprocessing of data affects the results. All three distinct planetary flow regimes identified in the winter circulation of the Southern Hemisphere by a pattern correlation method are captured by the new set of EOFs; only two of those regimes were captured by the earlier set. The new results, therefore, lend further support to the idea that EOFs point to distinct planetary flow regimes.

Farrara, John D.↗
NASA NTRS1988

Cluster analysis of multiple planetary flow regimes

A modified cluster analysis method developed for the classification of quasi-stationary events into a few planetary flow regimes and for the examination of transitions between these regimes is described. The method was applied first to a simple deterministic model and then to a 500-mbar data set for Northern Hemisphere (NH), for which cluster analysis was carried out in the subspace of the first seven empirical orthogonal functions (EOFs). Stationary clusters were found in the low-frequency band of more than 10 days, while transient clusters were found in the band-pass frequency window between 2.5 and 6 days. In the low-frequency band, three pairs of clusters determined EOFs 1, 2, and 3, respectively; they exhibited well-known regional features, such as blocking, the Pacific/North American pattern, and wave trains. Both model and low-pass data exhibited strong bimodality.

Mo, Kingtse↗
NASA NTRS1988

The generation mechanism of mixed Rossby-gravity waves in the equatorial troposphere

Numerical experiments are performed to clarify the excitation mechanism of mixed Rossby-gravity (Yanai) waves in the tropical troposphere, as well as the selection of zonal wavenumbers four through five and of the five-day period. The results show that the nonlinear wave-CISK mechanism can select the same wavenumber and period of Yanai waves as observed, that in the absence of asymmetries in lateral forcing, symmetric Kelvin modes are dominant over antisymmetric Yanai modes, and that if lateral antisymmetric forcing with arbitrary periods is given, the correct spectral peak of Yanai waves becomes apparent. It is shown how the nonlinear wave-CISK mechanism selects the observed spectral peaks in zonal wavenumber and frequency. The question of why symmetric modes dominate antisymmetric Yanai waves is addressed, as is the reason why the tropical atmosphere selects the periods of Yanai waves independently of those of lateral forcing.

Itoh, Hisanori↗
NASA NTRS1987

Statistics and dynamics of persistent anomalies

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.

Mo, Kingtse C.↗
NASA NTRS1987

Cluster analysis of multiple planetary flow regimes

A modified cluster analysis method was developed to identify spatial patterns of planetary flow regimes, and to study transitions between them. This method was applied first to a simple deterministic model and second to Northern Hemisphere (NH) 500 mb data. The dynamical model is governed by the fully-nonlinear, equivalent-barotropic vorticity equation on the sphere. Clusters of point in the model's phase space are associated with either a few persistent or with many transient events. Two stationary clusters have patterns similar to unstable stationary model solutions, zonal, or blocked. Transient clusters of wave trains serve as way stations between the stationary ones. For the NH data, cluster analysis was performed in the subspace of the first seven empirical orthogonal functions (EOFs). Stationary clusters are found in the low-frequency band of more than 10 days, and transient clusters in the bandpass frequency window between 2.5 and 6 days. In the low-frequency band three pairs of clusters determine, respectively, EOFs 1, 2, and 3. They exhibit well-known regional features, such as blocking, the Pacific/North American (PNA) pattern and wave trains. Both model and low-pass data show strong bimodality. Clusters in the bandpass window show wave-train patterns in the two jet exit regions. They are related, as in the model, to transitions between stationary clusters.

Mo, Kingtse↗
NASA NTRS1986

Dynamics, statistics and predictability of planetary flow regimes

Consideration is given to regimes of low-frequency variability in large-scale atmospheric dynamics. The model utilized is the fully-nonlinear, equivalent-baratropic vorticity equation on the sphere, with simplified forcing, dissipation and topography. It is found that certain limited regions in the system's phase space are visited repeatedly and for extended periods by model solutions. Flow patterns coupled with these regions correspond to synoptically-defined zonal and blocked Northern Hemisphere midlatitude flows. It is shown that the system's macrodynamics can be depicted by two or more planetary flow regimes, the expected residence time in each regime, and the transition properties from one regime to another. These model-derived ideas are also applied to a time series of atmospheric data from the Southern Hemisphere.

Ghil, Michael↗