Mickaël D. Chekroun

Predicting stochastic systems by noise sampling, and application to the El Niño-Southern Oscillation.

Interannual and interdecadal prediction are major challenges of climate dynamics. In this article we develop a prediction method for climate processes that exhibit low-frequency variability (LFV). The method constructs a nonlinear stochastic model from past observations and estimates a path of the “weather” noise that drives this model over previous finite-time windows. The method has two steps: (i) select noise samples—or “snippets”—from the past noise, which have forced the system during short-time intervals that resemble the LFV phase just preceding the currently observed state; and (ii) use these snippets to drive the system from the current state into the future. The method is placed in the framework of pathwise linear-response theory and is then applied to an El Niño–Southern Oscillation (ENSO) model derived by the empirical model reduction (EMR) methodology; this nonlinear model has 40 coupled, slow, and fast variables. The domain of validity of this forecasting procedure depends on the nature of the system’s pathwise response; it is shown numerically that the ENSO model’s response is linear on interannual time scales. As a result, the method’s skill at a 6- to 16-month lead is highly competitive when compared with currently used dynamic and statistic prediction methods for the Niño-3 index and the global sea surface temperature field.

Data-driven non-Markovian closure models

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.

Rough parameter dependence in climate models and the role of Ruelle-Pollicott resonances

Significance It is shown on a geophysical fluid model that the sensitive dependence on parametric variations of a given set of observations is related to the dominant modes of variability of the flow, as encoded through the Ruelle-Pollicott resonances. The latter are estimated by Markov modeling techniques. Such modeling and analysis methods can be used for any type of high-dimensional dissipative and chaotic dynamical systems. 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.

Parameter estimation for energy balance models with memory

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 PDEs 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 samples remoteness in time. Using a Bayesian approach, we show that obtaining an accurate parameter estimate is still possible in certain cases.

Stochastic Parameterizing Manifolds and Non-Markovian Reduced Equations

General Introduction.- Preliminaries.- Invariant Manifolds.- Pullback Characterization of Approximating, and Parameterizing Manifolds.- Non-Markovian Stochastic Reduced Equations.- On-Markovian Stochastic Reduced Equations on the Fly.- Proof of Lemma 5.1.-References.- Index.

Exploring the Pullback Attractors of a Low-Order Quasigeostrophic Ocean Model: The Deterministic Case

AbstractA low-order quasigeostrophic double-gyre ocean model is subjected to an aperiodic forcing that mimics time dependence dominated by interdecadal variability. This model is used as a prototype of an unstable and nonlinear dynamical system with time-dependent forcing to explore basic features of climate change in the presence of natural variability. The study relies on the theoretical framework of nonautonomous dynamical systems and of their pullback attractors (PBAs), that is, of the time-dependent invariant sets attracting all trajectories initialized in the remote past. The existence of a global PBA is rigorously demonstrated for this weakly dissipative nonlinear model. Ensemble simulations are carried out and the convergence to PBAs is assessed by computing the probability density function (PDF) of localization of the trajectories. A sensitivity analysis with respect to forcing amplitude shows that the PBAs experience large modifications if the underlying autonomous system is dominated by small-a...

Data-adaptive harmonic spectra and multilayer Stuart-Landau models.

Harmonic decompositions of multivariate time series are considered for which we adopt an integral operator approach with periodic semigroup kernels. Spectral decomposition theorems are derived that cover the important cases of two-time statistics drawn from a mixing invariant measure. The corresponding eigenvalues can be grouped per Fourier frequency and are actually given, at each frequency, as the singular values of a cross-spectral matrix depending on the data. These eigenvalues obey, furthermore, a variational principle that allows us to define naturally a multidimensional power spectrum. The eigenmodes, as far as they are concerned, exhibit a data-adaptive character manifested in their phase which allows us in turn to define a multidimensional phase spectrum. The resulting data-adaptive harmonic (DAH) modes allow for reducing the data-driven modeling effort to elemental models stacked per frequency, only coupled at different frequencies by the same noise realization. In particular, the DAH decomposition extracts time-dependent coefficients stacked by Fourier frequency which can be efficiently modeled-provided the decay of temporal correlations is sufficiently well-resolved-within a class of multilayer stochastic models (MSMs) tailored here on stochastic Stuart-Landau oscillators. Applications to the Lorenz 96 model and to a stochastic heat equation driven by a space-time white noise are considered. In both cases, the DAH decomposition allows for an extraction of spatio-temporal modes revealing key features of the dynamics in the embedded phase space. The multilayer Stuart-Landau models (MSLMs) are shown to successfully model the typical patterns of the corresponding time-evolving fields, as well as their statistics of occurrence.

Diversity, Nonlinearity, Seasonality, and Memory Effect in ENSO Simulation and Prediction Using Empirical Model Reduction

AbstractA suite of empirical model experiments under the empirical model reduction framework are conducted to advance the understanding of ENSO diversity, nonlinearity, seasonality, and the memory effect in the simulation and prediction of tropical Pacific sea surface temperature (SST) anomalies. The model training and evaluation are carried out using 4000-yr preindustrial control simulation data from the coupled model GFDL CM2.1. The results show that multivariate models with tropical Pacific subsurface information and multilevel models with SST history information both improve the prediction skill dramatically. These two types of models represent the ENSO memory effect based on either the recharge oscillator or the time-delayed oscillator viewpoint. Multilevel SST models are a bit more efficient, requiring fewer model coefficients. Nonlinearity is found necessary to reproduce the ENSO diversity feature for extreme events. The nonlinear models reconstruct the skewed probability density function of SST an...

Approximation of Stochastic Hyperbolic Invariant Manifolds

In this chapter, we consider situations where the subspace \(\fancyscript{H}^{\mathfrak {c}}\) spanned by the low modes contains a combination of critical modes and modes that remain stable as \(\lambda \) varies in some interval \(\varLambda \). We derive similar results to Theorem 6.1 and Corollary 6.1 concerning the leading-order approximation of the corresponding local stochastic invariant manifolds under suitable conditions. A possible relaxation of these conditions is also discussed. In particular, when the linear operator, \(L_\lambda \), is self-adjoint, some non-resonance (NR) conditions arise naturally. As explained in Volume II [37], these NR conditions play an important role in the construction of certain type of manifolds which may not necessarily approximate an invariant manifold, but still convey useful dynamical information; see [37, Chaps. 6–7]. This is formulated in Volume II via the concept of stochastic parameterizing manifolds; see [37, Sect. 4.2].

Pullback Attractor Crisis in a Delay Differential ENSO Model

We study the pullback attractor (PBA) of a seasonally forced delay differential model for the El Nino–Southern Oscillation (ENSO); the model has two delays, associated with a positive and a negative feedback. The control parameter is the intensity of the positive feedback and the PBA undergoes a crisis that consists of a chaos-to-chaos transition. Since the PBA is dominated by chaotic behavior, we refer to it as a strange PBA. Both chaotic regimes correspond to an overlapping of resonances but the two differ by the properties of this overlapping. The crisis manifests itself by a brutal change not only in the size but also in the shape of the PBA. The change is associated with the sudden disappearance of the most extreme warm (El Nino) and cold (La Nina) events, as one crosses the critical parameter value from below. The analysis reveals that regions of the strange PBA that survive the crisis are those populated by the most probable states of the system. These regions are those that exhibit robust foldings with respect to perturbations. The effect of noise on this phase-and-parameter space behavior is then discussed. It is shown that the chaos-to-chaos crisis may or may not survive the addition of small noise to the evolution equation, depending on how the noise enters the latter.