Publications by Author: Ghil, Michael

2015
Mukhin, Dmitry, Dmitri Kondrashov, Evgeny Loskutov, Andrey Gavrilov, Alexander Feigin, and Michael Ghil. “Predicting critical transitions in ENSO models. Part II: Spatially dependent models.” Journal of Climate 28, no. 5 (2015): 1962–1976. Abstract
The present paper is the second part of a two-part study on empirical modeling and prediction of climate variability. This paper deals with spatially distributed data, as opposed to the univariate data of Part I. The choice of a basis for effective data compression becomes of the essence. In many applications, it is the set of spatial empirical orthogonal functions that provides the uncorrelated time series of principal components (PCs) used in the learning set. In this paper, the basis of the learning set is obtained instead by applying multichannel singular-spectrum analysis to climatic time series and using the leading spatiotemporal PCs to construct a reduced stochastic model. The effectiveness of this approach is illustrated by predicting the behavior of the Jin–Neelin–Ghil (JNG) hybrid seasonally forced coupled ocean–atmosphere model of El Niño–Southern Oscillation. The JNG model produces spatially distributed and weakly nonstationary time series to which the model reduction and prediction methodology is applied. Critical transitions in the hybrid periodically forced coupled model are successfully predicted on time scales that are substantially longer than the duration of the learning sample.
Vannitsem, Stéphane, Jonathan Demaeyer, Lesley De Cruz, and Michael Ghil. “Low-frequency variability and heat transport in a low-order nonlinear coupled ocean–atmosphere model.” Physica D: Nonlinear Phenomena 309 (2015): 71–85. Abstract
We formulate and study a low-order nonlinear coupled ocean–atmosphere model with an emphasis on the impact of radiative and heat fluxes and of the frictional coupling between the two components. This model version extends a previous 24-variable version by adding a dynamical equation for the passive advection of temperature in the ocean, together with an energy balance model. The bifurcation analysis and the numerical integration of the model reveal the presence of low-frequency variability (LFV) concentrated on and near a long-periodic, attracting orbit. This orbit combines atmospheric and oceanic modes, and it arises for large values of the meridional gradient of radiative input and of frictional coupling. Chaotic behavior develops around this orbit as it loses its stability; this behavior is still dominated by the LFV on decadal and multi-decadal time scales that is typical of oceanic processes. Atmospheric diagnostics also reveals the presence of predominant low- and high-pressure zones, as well as of a subtropical jet; these features recall realistic climatological properties of the oceanic atmosphere. Finally, a predictability analysis is performed. Once the decadal-scale periodic orbits develop, the coupled system’s short-term instabilities–as measured by its Lyapunov exponents–are drastically reduced, indicating the ocean’s stabilizing role on the atmospheric dynamics. On decadal time scales, the recurrence of the solution in a certain region of the invariant subspace associated with slow modes displays some extended predictability, as reflected by the oscillatory behavior of the error for the atmospheric variables at long lead times.
Chang, C. P., Michael Ghil, M. Latif, and J. M. Wallace, ed. Climate Change: Multidecadal and Beyond. World Scientific Publ. Co./Imperial College Press, 2015.
Kondrashov, Dmitri, Mickaël D. Chekroun, and Michael Ghil. “Data-driven non-Markovian closure models.” Physica D: Nonlinear Phenomena 297 (2015): 33–55. Abstract

This paper has two interrelated foci: (i) obtaining stable and efficient data-driven closure models by using a multivariate time series of partial observations from a large-dimensional system; and (ii) comparing these closure models with the optimal closures predicted by the Mori–Zwanzig (MZ) formalism of statistical physics. Multilayer stochastic models (MSMs) are introduced as both a generalization and a time-continuous limit of existing multilevel, regression-based approaches to closure in a data-driven setting; these approaches include empirical model reduction (EMR), as well as more recent multi-layer modeling. It is shown that the multilayer structure of MSMs can provide a natural Markov approximation to the generalized Langevin equation (GLE) of the MZ formalism. A simple correlation-based stopping criterion for an EMR–MSM model is derived to assess how well it approximates the GLE solution. Sufficient conditions are derived on the structure of the nonlinear cross-interactions between the constitutive layers of a given MSM to guarantee the existence of a global random attractor. This existence ensures that no blow-up can occur for a broad class of MSM applications, a class that includes non-polynomial predictors and nonlinearities that do not necessarily preserve quadratic energy invariants. The EMR–MSM methodology is first applied to a conceptual, nonlinear, stochastic climate model of coupled slow and fast variables, in which only slow variables are observed. It is shown that the resulting closure model with energy-conserving nonlinearities efficiently captures the main statistical features of the slow variables, even when there is no formal scale separation and the fast variables are quite energetic. Second, an MSM is shown to successfully reproduce the statistics of a partially observed, generalized Lotka–Volterra model of population dynamics in its chaotic regime. The challenges here include the rarity of strange attractors in the model’s parameter space and the existence of multiple attractor basins with fractal boundaries. The positivity constraint on the solutions’ components replaces here the quadratic-energy–preserving constraint of fluid-flow problems and it successfully prevents blow-up.

Chavez, M., Michael Ghil, and J. Urrutia-Fucugauchi, ed. Extreme Events: Observations, Modeling and Economics. Geophysical Monographs. Vol. 214. Washington, DC: American Geophysical Union & Wiley, 2015. Publisher's Version Abstract

The monograph covers the fundamentals and the consequences of extreme geophysical phenomena like asteroid impacts, climatic change, earthquakes, tsunamis, hurricanes, landslides, volcanic eruptions, flooding, and space weather. This monograph also addresses their associated, local and worldwide socio-economic impacts. The understanding and modeling of these phenomena is critical to the development of timely worldwide strategies for the prediction of natural and anthropogenic extreme events, in order to mitigate their adverse consequences. This monograph is unique in as much as it is dedicated to recent theoretical, numerical and empirical developments that aim to improve: (i) the understanding, modeling and prediction of extreme events in the geosciences, and, (ii) the quantitative evaluation of their economic consequences. The emphasis is on coupled, integrative assessment of the physical phenomena and their socio-economic impacts. With its overarching theme, Extreme Events: Observations, Modeling and Economics will be relevant to and become an important tool for researchers and practitioners in the fields of hazard and risk analysis in general, as well as to those with a special interest in climate change, atmospheric and oceanic sciences, seismo-tectonics, hydrology, and space weather.

Groth, Andreas, Patrice Dumas, Michael Ghil, and Stéphane Hallegatte. “Impacts of natural disasters on a dynamic economy.” In Extreme Events : Observations, Modeling, and Economics, edited by Eric Chavez, Michael Ghil, and Jaime Urrutia-Fucugauchi, 343–360. American Geophysical Union and Wiley-Blackwell, 2015. Abstract

This chapter presents a modeling framework for macroeconomic growth dynamics; it is motivated by recent attempts to formulate and study “integrated models” of the coupling between natural and socioeconomic phe­ nomena. The challenge is to describe the interfaces between human activities and the functioning of the earth system. We examine the way in which this interface works in the presence of endogenous business cycle dynam­ ics, based on a nonequilibrium dynamic model. Recent findings about the macroeconomic response to natural disasters in such a nonequilibrium setting have shown a more severe response to natural disasters during expan­ sions than during recessions. These findings raise questions about the assessment of climate change damages or natural disaster losses that are based purely on long-term growth models. In order to compare the theoretical findings with observational data, we analyze cyclic behavior in the U.S. economy, based on multivariate singular spectrum analysis. We analyze a total of nine aggregate indicators in a 52 year interval (1954–2005) and demon­ strate that the behavior of the U.S. economy changes significantly between intervals of growth and recession, with higher volatility during expansions.

PDF
Groth, Andreas, Patrice Dumas, Michael Ghil, and Stéphane Hallegatte. “Impacts of natural disasters on a dynamic economy.” In Extreme Events : Observations, Modeling, and Economics, edited by Eric Chavez, Michael Ghil, and Jaime Urrutia-Fucugauchi, 343–360. American Geophysical Union and Wiley-Blackwell, 2015. Abstract

This chapter presents a modeling framework for macroeconomic growth dynamics; it is motivated by recent attempts to formulate and study “integrated models” of the coupling between natural and socioeconomic phe­ nomena. The challenge is to describe the interfaces between human activities and the functioning of the earth system. We examine the way in which this interface works in the presence of endogenous business cycle dynam­ ics, based on a nonequilibrium dynamic model. Recent findings about the macroeconomic response to natural disasters in such a nonequilibrium setting have shown a more severe response to natural disasters during expan­ sions than during recessions. These findings raise questions about the assessment of climate change damages or natural disaster losses that are based purely on long-term growth models. In order to compare the theoretical findings with observational data, we analyze cyclic behavior in the U.S. economy, based on multivariate singular spectrum analysis. We analyze a total of nine aggregate indicators in a 52 year interval (1954–2005) and demon­ strate that the behavior of the U.S. economy changes significantly between intervals of growth and recession, with higher volatility during expansions.

PDF
Ghil, Michael. “A mathematical theory of climate sensitivity or, How to deal with both anthropogenic forcing and natural variability?” In Climate Change: Multidecadal and Beyond, edited by C. P. Chang, Michael Ghil, Mojib Latif, and J. M. Wallace, 31–51. World Scientific Publ. Co./Imperial College Press, 2015. Abstract

Recent estimates of climate evolution over the coming century still differ by several degrees. This uncertainty motivates the work presented here. There are two basic approaches to apprehend the complexity of climate change: deterministically nonlinear and stochastically linear, i.e., the Lorenz and the Hasselmann approach. The grand unification of these two approaches relies on the theory of random dynamical systems. We apply this theory to study the random attractors of nonlinear, stochastically perturbed climate models. Doing so allows one to examine the interaction of internal climate variability with the forcing, whether natural or anthropogenic, and to take into account the climate system's non-equilibrium behavior in determining climate sensitivity. This non-equilibrium behavior is due to a combination of nonlinear and random effects. We give here a unified treatment of such effects from the point of view of the theory of dynamical systems and of their bifurcations. Energy balance models are used to illustrate multiple equilibria, while multi-decadal oscillations in the thermohaline circulation illustrate the transition from steady states to periodic behavior. Random effects are introduced in the setting of random dynamical systems, which permit a unified treatment of both nonlinearity and stochasticity. The combined treatment of nonlinear and random effects is applied to a stochastically perturbed version of the classical Lorenz convection model. Climate sensitivity is then defined mathematically as the derivative of an appropriate functional or other function of the system’s state with respect to the bifurcation parameter. This definition is illustrated by using numerical results for a model of the El Niño–Southern Oscillation. The concept of a hierarchy of models is the thread that runs across this chapter, and the robustness of elementary bifurcations across such a hierarchy is emphasized.

PDF
Ghil, Michael. “A mathematical theory of climate sensitivity or, How to deal with both anthropogenic forcing and natural variability?” In Climate Change: Multidecadal and Beyond, edited by C. P. Chang, Michael Ghil, Mojib Latif, and J. M. Wallace, 31–51. World Scientific Publ. Co./Imperial College Press, 2015. Abstract

Recent estimates of climate evolution over the coming century still differ by several degrees. This uncertainty motivates the work presented here. There are two basic approaches to apprehend the complexity of climate change: deterministically nonlinear and stochastically linear, i.e., the Lorenz and the Hasselmann approach. The grand unification of these two approaches relies on the theory of random dynamical systems. We apply this theory to study the random attractors of nonlinear, stochastically perturbed climate models. Doing so allows one to examine the interaction of internal climate variability with the forcing, whether natural or anthropogenic, and to take into account the climate system's non-equilibrium behavior in determining climate sensitivity. This non-equilibrium behavior is due to a combination of nonlinear and random effects. We give here a unified treatment of such effects from the point of view of the theory of dynamical systems and of their bifurcations. Energy balance models are used to illustrate multiple equilibria, while multi-decadal oscillations in the thermohaline circulation illustrate the transition from steady states to periodic behavior. Random effects are introduced in the setting of random dynamical systems, which permit a unified treatment of both nonlinearity and stochasticity. The combined treatment of nonlinear and random effects is applied to a stochastically perturbed version of the classical Lorenz convection model. Climate sensitivity is then defined mathematically as the derivative of an appropriate functional or other function of the system’s state with respect to the bifurcation parameter. This definition is illustrated by using numerical results for a model of the El Niño–Southern Oscillation. The concept of a hierarchy of models is the thread that runs across this chapter, and the robustness of elementary bifurcations across such a hierarchy is emphasized.

PDF
Groth, Andreas, and Michael Ghil. “Monte Carlo Singular Spectrum Analysis (SSA) revisited: Detecting oscillator clusters in multivariate datasets.” Journal of Climate 28, no. 19 (2015): 7873–7893. 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.

PDF
Groth, Andreas, Michael Ghil, Stéphane Hallegatte, and Patrice Dumas. “The Role of Oscillatory Modes in U.S. Business Cycles.” OECD Journal: Journal of Business Cycle Measurement and Analysis, no. 2015/1 (2015): 63–81. Abstract

We apply multivariate singular spectrum analysis to the study of U.S. business cycle dynamics. This method provides a robust way to identify and reconstruct oscillations, whether intermittent or modulated. We show such oscillations to be associated with comovements across the entire economy. The problem of spurious cycles generated by the use of detrending filters is addressed and we present a Monte Carlo test to extract significant oscillations. The behavior of the U.S. economy is shown to change significantly from one phase of the business cycle to another: the recession phase is dominated by a five-year mode, while the expansion phase exhibits more complex dynamics, with higher-frequency modes coming into play. We show that the variations so identified cannot be generated by random shocks alone, as assumed in ‘real’ business-cycle models, and that endogenous, deterministically generated variability has to be involved.

PDF
Ghil, Michael, and I. Zaliapin. “Understanding ENSO variability and its extrema: A delay differential equation approach.” In Extreme Events: Observations, Modeling and Economics, edited by M. Chavez, Michael Ghil, and J. Urrutia-Fucugauchi, 63–78. American Geophysical Union & Wiley, 2015. Abstract

The El-Nino/Southern-Oscillation (ENSO) phenomenon is the most prominent signal of seasonal-to-interannual climate variability. The past 30 years of research have shown that ENSO dynamics is governed, by and large, by the interplay of the nonlinear mechanisms, and that their simplest version can be studied in autonomous or forced delay differential equation (DDE) models. This chapter briefly reviews the results of Ghil et al., Zaliapin and Ghil, and Ghil and Zaliapin and pursues their DDE model analysis by focusing on multiple model solutions for the same parameter values and the dynamics of local extrema. It first introduces the DDE model of ENSO variability, reviews the main theoretical results concerning its solutions, and comments on the appropriate numerical integration methods. Novel results on multiple solutions and their extrema are reported and illustrated. After discussing the model's pullback attractor, the chapter explores parameter dependence in the model over its entire 3D parameter space.

PDF
Ghil, Michael, and I. Zaliapin. “Understanding ENSO variability and its extrema: A delay differential equation approach.” In Extreme Events: Observations, Modeling and Economics, edited by M. Chavez, Michael Ghil, and J. Urrutia-Fucugauchi, 63–78. American Geophysical Union & Wiley, 2015. Abstract

The El-Nino/Southern-Oscillation (ENSO) phenomenon is the most prominent signal of seasonal-to-interannual climate variability. The past 30 years of research have shown that ENSO dynamics is governed, by and large, by the interplay of the nonlinear mechanisms, and that their simplest version can be studied in autonomous or forced delay differential equation (DDE) models. This chapter briefly reviews the results of Ghil et al., Zaliapin and Ghil, and Ghil and Zaliapin and pursues their DDE model analysis by focusing on multiple model solutions for the same parameter values and the dynamics of local extrema. It first introduces the DDE model of ENSO variability, reviews the main theoretical results concerning its solutions, and comments on the appropriate numerical integration methods. Novel results on multiple solutions and their extrema are reported and illustrated. After discussing the model's pullback attractor, the chapter explores parameter dependence in the model over its entire 3D parameter space.

PDF
2014
Roques, Lionel, Mickaël D. Chekroun, Michel Cristofol, Samuel Soubeyrand, and Michael Ghil. “Parameter estimation for energy balance models with memory.” Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences 470, no. 2169 (2014). Publisher's Version Abstract
We study parameter estimation for one-dimensional energy balance models with memory (EBMMs) given localized and noisy temperature measurements. Our results apply to a wide range of nonlinear, parabolic partial differential equations with integral memory terms. First, we show that a space-dependent parameter can be determined uniquely everywhere in the PDE’s domain of definition D, using only temperature information in a small subdomain E⊂D. This result is valid only when the data correspond to exact measurements of the temperature. We propose a method for estimating a model parameter of the EBMM using more realistic, error-contaminated temperature data derived, for example, from ice cores or marine-sediment cores. Our approach is based on a so-called mechanistic-statistical model that combines a deterministic EBMM with a statistical model of the observation process. Estimating a parameter in this setting is especially challenging, because the observation process induces a strong loss of information. Aside from the noise contained in past temperature measurements, an additional error is induced by the age-dating method, whose accuracy tends to decrease with a sample’s remoteness in time. Using a Bayesian approach, we show that obtaining an accurate parameter estimate is still possible in certain cases.
Chekroun, Mickaël D., J. David Neelin, Dmitri Kondrashov, James C. McWilliams, and Michael Ghil. “Rough parameter dependence in climate models and the role of Ruelle-Pollicott resonances.” Proceedings of the National Academy of Sciences 111, no. 5 (2014): 1684-1690. Abstract

Despite the importance of uncertainties encountered in climate model simulations, the fundamental mechanisms at the origin of sensitive behavior of long-term model statistics remain unclear. Variability of turbulent flows in the atmosphere and oceans exhibits recurrent large-scale patterns. These patterns, while evolving irregularly in time, manifest characteristic frequencies across a large range of time scales, from intraseasonal through interdecadal. Based on modern spectral theory of chaotic and dissipative dynamical systems, the associated low-frequency variability may be formulated in terms of Ruelle-Pollicott (RP) resonances. RP resonances encode information on the nonlinear dynamics of the system, and an approach for estimating them—as filtered through an observable of the system—is proposed. This approach relies on an appropriate Markov representation of the dynamics associated with a given observable. It is shown that, within this representation, the spectral gap—defined as the distance between the subdominant RP resonance and the unit circle—plays a major role in the roughness of parameter dependences. The model statistics are the most sensitive for the smallest spectral gaps; such small gaps turn out to correspond to regimes where the low-frequency variability is more pronounced, whereas autocorrelations decay more slowly. The present approach is applied to analyze the rough parameter dependence encountered in key statistics of an El-Niño–Southern Oscillation model of intermediate complexity. Theoretical arguments, however, strongly suggest that such links between model sensitivity and the decay of correlation properties are not limited to this particular model and could hold much more generally.

2013
Ghil, Michael. “Lecture 1: Data Assimilation: How We Got Here and Where To Next?Workshop on Mathematics of Climate Change, Related Hazards and Risks, CIMAT, Guanajuato, Mexico, 2013. Abstract

Lecture 1: Data Assimilation: How We Got Here and Where To Next?
Ghil, Michael. “Lecture 2: Toward a Mathematical Theory of Climate Sensitivity.” Workshop on Mathematics of Climate Change, Related Hazards and Risks, CIMAT, Guanajuato, Mexico, 2013. Abstract

Lecture 2: Toward a Mathematical Theory of Climate Sensitivity
Ghil, Michael. “Lecture 3 : The Coupled Dynamics of Climate and Economics.” Workshop on Mathematics of Climate Change, Related Hazards and Risks, CIMAT, Guanajuato, Mexico, 2013. Abstract

Lecture 3 : The Coupled Dynamics of Climate and Economics
Sella, Lisa, Gianna Vivaldo, Andreas Groth, and Michael Ghil. “Economic Cycles and their Synchronization: A spectral survey.” Fondazione Eni Enrico Mattei (FEEM) 105, no. 105 (2013): 1. Publisher's Version Abstract

The present work applies several advanced spectral methods to the analysis of macroeconomic fluctuations in three countries of the European Union: Italy, The Netherlands, and the United Kingdom. We focus here in particular on singular-spectrum analysis (SSA), which provides valuable spatial and frequency information of multivariate data and that goes far beyond a pure analysis in the time domain. The spectral methods discussed here are well established in the geosciences and life sciences, but not yet widespread in quantitative economics. In particular, they enable one to identify and describe nonlinear trends and dominant cycles –- including seasonal and interannual components –- that characterize the deterministic behavior of each time series. These tools have already proven their robustness in the application on short and noisy data, and we demonstrate their usefulness in the analysis of the macroeconomic indicators of these three countries. We explore several fundamental indicators of the countries' real aggregate economy in a univariate, as well as a multivariate setting. Starting with individual single-channel analysis, we are able to identify similar spectral components among the analyzed indicators. Next, we consider combinations of indicators and countries, in order to take different effects of comovements into account. Since business cycles are cross-national phenomena, which show common characteristics across countries, our aim is to uncover hidden global behavior across the European economies. Results are compared with previous findings on the U.S. indicators \citepGroth.ea.FEEM.2012. Finally, the analysis is extended to include several indicators from the U.S. economy, in order to examine its influence on the European market.

PDF
de Viron, O., J. O. Dickey, and Michael Ghil. “Global modes of climate variability.” Geophysical Research Letters 40, no. 9 (2013): 1832-1837. Abstract

The atmosphere, hydrosphere and cryosphere form a fully coupled climate system. This system exhibits a number of large-scale phenomena, such as the El Nino Southern Oscillation (ENSO), the Asian Monsoon, the North Atlantic Oscillation (NAO), and the Madden-Julian Oscillation (MJO). While these modes of variability are not exactly periodic, they are oscillatory in character, and their state is monitored using so-called climate indices. Each of these scalar indices is a combination of several climate variables. Here, we use a comprehensive set of 25 climate indices for time intervals that range between 1948 and 2011, and estimate an optimal set of lags between these indices to maximize their correlation. We show that most of the index pairs drawn from this set present a significant correlation on interannual time scales. It is also shown that, on average, about two-thirds of the total variability in each index can be described by using only the four leading principal components of the entire set of lagged indices. Our index set's leading orthogonal modes exhibit several interannual frequencies and capture separately variability associated with the North Atlantic and the North Pacific. These modes are associated, in turn, with large-scale variations of sea surface temperatures.

PDF

Pages