Identifying efficient ensemble perturbations for initializing subseasonal-to-seasonal prediction

Abstract

The prediction of the weather at subseasonal-to-seasonal (S2S) timescales is dependent on both initial and boundary conditions. An open question is how to best initialize a relatively small-sized ensemble of numerical model integrations to produce reliable forecasts at these timescales. Reliability in this case means that the statistical properties of the ensemble forecast are consistent with the actual uncertainties about the future state of the geophysical system under investigation. In the present work, a method is introduced to construct initial conditions that produce reliable ensemble forecasts by projecting onto the eigenfunctions of the Koopman or the Perron-Frobenius operators, which describe the time-evolution of observables and probability distributions of the system dynamics, respectively. These eigenfunctions can be approximated from data by using the Dynamic Mode Decomposition (DMD) algorithm. The effectiveness of this approach is illustrated in the framework of a low-order ocean-atmosphere model exhibiting multiple characteristic timescales, and is compared to other ensemble initialization methods based on the Empirical Orthogonal Functions (EOFs) of the model trajectory and on the backward and covariant Lyapunov vectors of the model dynamics. Projecting initial conditions onto a subset of the Koopman or Perron-Frobenius eigenfunctions that are characterized by time scales with fast-decaying oscillations is found to produce highly reliable forecasts at all lead times investigated, ranging from one week to two months. Reliable forecasts are also obtained with the adjoint covariant Lyapunov vectors, which are the eigenfunctions of the Koopman operator in the tangent space. The advantages of these different methods are discussed.

arXiv commentSubmitted to JAMES, 35 pages, 12 figures. 9 pages supplementary material at the end

References in corpus (19)

• A Data-Driven Approximation of the Koopman Operator: Extending Dynamic Mode Decomposition
• On Dynamic Mode Decomposition: Theory and Applications
• Data-driven model reduction and transfer operator approximation
• On the numerical approximation of the Perron-Frobenius and Koopman operator
• Theory and computation of covariant Lyapunov vectors
• Study of dynamics in post-transient flows using Koopman mode decomposition
• Hyperbolicity and the effective dimension of spatially-extended dissipative systems
• Applied Koopman Operator Theory for Power Systems Technology
• Predictability of large-scale atmospheric motions: Lyapunov exponents and error dynamics
• Statistical and Dynamical Properties of Covariant Lyapunov Vectors in a Coupled Atmosphere-Ocean Model - Multiscale Effects, Geometric Degeneracy, and Error Dynamics
• Low-frequency variability and heat transport in a low-order nonlinear coupled ocean-atmosphere model
• Reduced-Order Models for Coupled Dynamical Systems: Data-driven Methods and the Koopman Operator
• Degenerate Kalman filter error covariances and their convergence onto the unstable subspace
• Crisis of the Chaotic Attractor of a Climate Model: A Transfer Operator Approach
• The Modular Arbitrary-Order Ocean-Atmosphere Model: MAOOAM v1.0
• Resonances in a Chaotic Attractor Crisis of the Lorenz Flow
• On temporal scale separation in coupled data assimilation with the ensemble Kalman filter
• Stochastic parameterization of subgrid-scale processes in coupled ocean-atmosphere systems: Benefits and limitations of response theory
• On the use of near-neutral Backward Lyapunov Vectors to get reliable ensemble forecasts in coupled ocean-atmosphere systems