Source author record

Michael Ghil

Michael Ghil appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
16topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2023arXiv

A multi-model ensemble Kalman filter for data assimilation and forecasting

Data assimilation (DA) aims to optimally combine model forecasts and observations that are both partial and noisy. Multi-model DA generalizes the variational or Bayesian formulation of the Kalman filter, and we prove that it is also the minimum variance linear unbiased estimator. Here, we formulate and implement a multi-model ensemble Kalman filter (MM-EnKF) based on this framework. The MM-EnKF can combine multiple model ensembles for both DA and forecasting in a flow-dependent manner; it uses adaptive model error estimation to provide matrix-valued weights for the separate models and the observations. We apply this methodology to various situations using the Lorenz96 model for illustration purposes. Our numerical experiments include multiple models with parametric error, different resolved scales, and different fidelities. The MM-EnKF results in significant error reductions compared to the best model, as well as to an unweighted multi-model ensemble, with respect to both probabilistic and deterministic error metrics.

preprint2022arXiv

PaleoJump: A database for abrupt transitions in past climates

Tipping points (TPs) in the Earth system have been studied with growing interest and concern in recent years due to the potential risk of anthropogenic forcing causing abrupt, and possibly irreversible, climate transitions. Paleoclimate records are essential for identifying TPs in Earth's past and for properly understanding the climate system's underlying nonlinearities and bifurcation mechanisms. Due to the variations in quality, resolution, and dating methods, it is crucial to select the records that give the best representation of past climates. Furthermore, as paleoclimate time series vary in their origin, time spans, and periodicities, an objective, automated methodology is crucial for identifying and comparing TPs. To reach this goal, we present here the PaleoJump database of carefully selected, high-resolution records originating in ice cores, marine sediments, speleothems, terrestrial records, and lake sediments. These records describe climate variability on centennial, millennial, or longer time scales and cover all the continents and ocean basins. We provide an overview of their spatial distribution and discuss the gaps in coverage. Our statistical methodology includes an augmented Kolmogorov-Smirnov test and Recurrence Quantification Analysis; it is applied here to selected records to automatically detect abrupt transitions therein and to investigate the presence of potential tipping elements. These transitions are shown in the PaleoJump database together with other essential information, including location, temporal scale and resolution, along with temporal plots. This database represents, therefore, a valuable resource for researchers investigating TPs in past climates.

preprint2022arXiv

Random templex encodes topological tipping points in noise-driven chaotic dynamics

Random attractors are the time-evolving pullback attractors of stochastically perturbed, deterministically chaotic dynamical systems. These attractors have a structure that changes in time, and that has been characterized recently using {\sc BraMAH} cell complexes and their homology groups. This description has been further improved for their deterministic counterparts by endowing the cell complex with a directed graph, which encodes the order in which the cells in the complex are visited by the flow in phase space. A templex is a mathematical object formed by a complex and a digraph; it provides a finer description of deterministically chaotic attractors and permits their accurate classification. In a deterministic framework, the digraph of the templex connects cells within a single complex for all time. Here, we introduce the stochastic version of a templex. In a random templex, there is one complex per snapshot of the random attractor and the digraph connects the generators or ``holes'' of successive cell complexes. Tipping points appear in a random templex as drastic changes of its holes in motion, namely their birth, splitting, merging, or death. This paper introduces and computes the random templex for the noise-driven Lorenz system's random attractor (LORA).

preprint2020arXiv

The Physics of Climate Variability and Climate Change

The climate system is a forced, dissipative, nonlinear, complex and heterogeneous system that is out of thermodynamic equilibrium. The system exhibits natural variability on many scales of motion, in time as well as space, and it is subject to various external forcings, natural as well as anthropogenic. This paper reviews the observational evidence on climate phenomena and the governing equations of planetary-scale flow, as well as presenting the key concept of a hierarchy of models as used in the climate sciences. Recent advances in the application of dynamical systems theory, on the one hand, and of nonequilibrium statistical physics, on the other, are brought together for the first time and shown to complement each other in helping understand and predict the system's behavior. These complementary points of view permit a self-consistent handling of subgrid-scale phenomena as stochastic processes, as well as a unified handling of natural climate variability and forced climate change, along with a treatment of the crucial issues of climate sensitivity, response, and predictability.

preprint2015arXiv

DADA: Data Assimilation for the Detection and Attribution of Weather- and Climate-related Events

We describe a new approach allowing for systematic causal attribution of weather and climate-related events, in near-real time. The method is purposely designed to facilitate its implementation at meteorological centers by relying on data treatments that are routinely performed when numerically forecasting the weather. Namely, we show that causal attribution can be obtained as a by-product of so-called data assimilation procedures that are run on a daily basis to update the meteorological model with new atmospheric observations; hence, the proposed methodology can take advantage of the powerful computational and observational capacity of weather forecasting centers. We explain the theoretical rationale of this approach and sketch the most prominent features of a "data assimilation-based detection and attribution" (DADA) procedure. The proposal is illustrated in the context of the classical three-variable Lorenz model with additional forcing. Several theoretical and practical research questions that need to be addressed to make the proposal readily operational within weather forecasting centers are finally laid out.

preprint2015arXiv

Low-Dimensional Galerkin Approximations of Nonlinear Delay Differential Equations

This article revisits the approximation problem of systems of nonlinear delay differential equations (DDEs) by a set of ordinary differential equations (ODEs). We work in Hilbert spaces endowed with a natural inner product including a point mass, and introduce polynomials orthogonal with respect to such an inner product that live in the domain of the linear operator associated with the underlying DDE. These polynomials are then used to design a general Galerkin scheme for which we derive rigorous convergence results and show that it can be numerically implemented via simple analytic formulas. The scheme so obtained is applied to three nonlinear DDEs, two autonomous and one forced: (i) a simple DDE with distributed delays whose solutions recall Brownian motion; (ii) a DDE with a discrete delay that exhibits bimodal and chaotic dynamics; and (iii) a periodically forced DDE with two discrete delays arising in climate dynamics. In all three cases, the Galerkin scheme introduced in this article provides a good approximation by low-dimensional ODE systems of the DDE's strange attractor, as well as of the statistical features that characterize its nonlinear dynamics.

preprint2015arXiv

Low-frequency variability and heat transport in a low-order nonlinear coupled ocean-atmosphere model

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.

preprint2014arXiv

Data-driven non-Markovian closure models

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 very broad class of MSM applications, a class that includes non-polynomial predictors and nonlinearities that do not necessarily preserve quadratic energy invariants. The approach is applied to conceptual nonlinear models borrowed from climate dynamics and population dynamics. In both cases, it is shown that the resulting closure models are able to capture the main statistical features of the dynamics, even in presence of weak time-scale separation.

preprint2010arXiv

A delay differential model of ENSO variability, Part 2: Phase locking, multiple solutions, and dynamics of extrema

We consider a highly idealized model for El Nino/Southern Oscillation (ENSO) variability, as introduced in an earlier paper. The model is governed by a delay differential equation for sea surface temperature in the Tropical Pacific, and it combines two key mechanisms that participate in ENSO dynamics: delayed negative feedback and seasonal forcing. We perform a theoretical and numerical study of the model in the three-dimensional space of its physically relevant parameters: propagation period of oceanic waves across the Tropical Pacific, atmosphere-ocean coupling, and strength of seasonal forcing. Phase locking of model solutions to the periodic forcing is prevalent: the local maxima and minima of the solutions tend to occur at the same position within the seasonal cycle. Such phase locking is a key feature of the observed El Nino (warm) and La Nina (cold) events. The phasing of the extrema within the seasonal cycle depends sensitively on model parameters when forcing is weak. We also study co-existence of multiple solutions for fixed model parameters and describe the basins of attraction of the stable solutions in a one-dimensional space of constant initial model histories.

preprint2010arXiv

Another Look at Climate Sensitivity

We revisit a recent claim that the Earth's climate system is characterized by sensitive dependence to parameters; in particular, that the system exhibits an asymmetric, large-amplitude response to normally distributed feedback forcing. Such a response would imply irreducible uncertainty in climate change predictions and thus have notable implications for climate science and climate-related policy making. We show that equilibrium climate sensitivity in all generality does not support such an intrinsic indeterminacy; the latter appears only in essentially linear systems. The main flaw in the analysis that led to this claim is inappropriate linearization of an intrinsically nonlinear model; there is no room for physical interpretations or policy conclusions based on this mathematical error. Sensitive dependence nonetheless does exist in the climate system, as well as in climate models -- albeit in a very different sense from the one claimed in the linear work under scrutiny -- and we illustrate it using a classical energy balance model (EBM) with nonlinear feedbacks. EBMs exhibit two saddle-node bifurcations, more recently called "tipping points", which give rise to three distinct steady-state climates, two of which are stable. Such bistable behavior is, furthermore, supported by results from more realistic, nonequilibrium climate models. In a truly nonlinear setting, indeterminacy in the size of the response is observed only in the vicinity of tipping points. We show, in fact, that small disturbances cannot result in a large-amplitude response, unless the system is at or near such a point. We discuss briefly how the distance to the bifurcation may be related to the strength of Earth's ice-albedo feedback.

preprint2010arXiv

Averaging of time-periodic systems without a small parameter

In this article, we present a new approach to averaging in non-Hamiltonian systems with periodic forcing. The results here do not depend on the existence of a small parameter. In fact, we show that our averaging method fits into an appropriate nonlinear equivalence problem, and that this problem can be solved formally by using the Lie transform framework to linearize it. According to this approach, we derive formal coordinate transformations associated with both first-order and higher-order averaging, which result in more manageable formulae than the classical ones. Using these transformations, it is possible to correct the solution of an averaged system by recovering the oscillatory components of the original non-averaged system. In this framework, the inverse transformations are also defined explicitly by formal series; they allow the estimation of appropriate initial data for each higher-order averaged system, respecting the equivalence relation. Finally, we show how these methods can be used for identifying and computing periodic solutions for a very large class of nonlinear systems with time-periodic forcing. We test the validity of our approach by analyzing both the first-order and the second-order averaged system for a problem in atmospheric chemistry.

preprint2010arXiv

Climate dynamics and fluid mechanics: Natural variability and related uncertainties

The purpose of this review-and-research paper is twofold: (i) to review the role played in climate dynamics by fluid-dynamical models; and (ii) to contribute to the understanding and reduction of the uncertainties in future climate-change projections. To illustrate the first point, we focus on the large-scale, wind-driven flow of the mid-latitude oceans which contribute in a crucial way to Earth's climate, and to changes therein. We study the low-frequency variability (LFV) of the wind-driven, double-gyre circulation in mid-latitude ocean basins, via the bifurcation sequence that leads from steady states through periodic solutions and on to the chaotic, irregular flows documented in the observations. This sequence involves local, pitchfork and Hopf bifurcations, as well as global, homoclinic ones. The natural climate variability induced by the LFV of the ocean circulation is but one of the causes of uncertainties in climate projections. Another major cause of such uncertainties could reside in the structural instability in the topological sense, of the equations governing climate dynamics, including but not restricted to those of atmospheric and ocean dynamics. We propose a novel approach to understand, and possibly reduce, these uncertainties, based on the concepts and methods of random dynamical systems theory. As a very first step, we study the effect of noise on the topological classes of the Arnol'd family of circle maps, a paradigmatic model of frequency locking as occurring in the nonlinear interactions between the El Nino-Southern Oscillations (ENSO) and the seasonal cycle. It is shown that the maps' fine-grained resonant landscape is smoothed by the noise, thus permitting their coarse-grained classification. This result is consistent with stabilizing effects of stochastic parametrization obtained in modeling of ENSO phenomenon via some general circulation models.

preprint2009arXiv

Transport on river networks: A dynamical approach

This study is motivated by problems related to environmental transport on river networks. We establish statistical properties of a flow along a directed branching network and suggest its compact parameterization. The downstream network transport is treated as a particular case of nearest-neighbor hierarchical aggregation with respect to the metric induced by the branching structure of the river network. We describe the static geometric structure of a drainage network by a tree, referred to as the static tree, and introduce an associated dynamic tree that describes the transport along the static tree. It is well known that the static branching structure of river networks can be described by self-similar trees (SSTs); we demonstrate that the corresponding dynamic trees are also self-similar. We report an unexpected phase transition in the dynamics of three river networks, one from California and two from Italy, demonstrate the universal features of this transition, and seek to interpret it in hydrological terms.

preprint2007arXiv

A delay differential model of ENSO variability: Parametric instability and the distribution of extremes

We consider a delay differential equation (DDE) model for El-Nino Southern Oscillation (ENSO) variability. The model combines two key mechanisms that participate in ENSO dynamics: delayed negative feedback and seasonal forcing. We perform stability analyses of the model in the three-dimensional space of its physically relevant parameters. Our results illustrate the role of these three parameters: strength of seasonal forcing $b$, atmosphere-ocean coupling $κ$, and propagation period $τ$ of oceanic waves across the Tropical Pacific. Two regimes of variability, stable and unstable, are separated by a sharp neutral curve in the $(b,τ)$ plane at constant $κ$. The detailed structure of the neutral curve becomes very irregular and possibly fractal, while individual trajectories within the unstable region become highly complex and possibly chaotic, as the atmosphere-ocean coupling $κ$ increases. In the unstable regime, spontaneous transitions occur in the mean ``temperature'' ({\it i.e.}, thermocline depth), period, and extreme annual values, for purely periodic, seasonal forcing. The model reproduces the Devil's bleachers characterizing other ENSO models, such as nonlinear, coupled systems of partial differential equations; some of the features of this behavior have been documented in general circulation models, as well as in observations. We expect, therefore, similar behavior in much more detailed and realistic models, where it is harder to describe its causes as completely.