We present and apply a novel method of describing and modeling complex multivariate datasets in the geosciences and elsewhere. Data-adaptive harmonic (DAH) decomposition identifies narrow-banded, spatio-temporal modes (DAHMs) whose frequencies are not necessarily integer multiples of each other. The evolution in time of the DAH coefficients (DAHCs) of these modes can be modeled using a set of coupled Stuart-Landau stochastic differential equations that capture the modes’ frequencies and amplitude modulation in time and space. This methodology is applied first to a challenging synthetic dataset and then to Arctic sea ice concentration (SIC) data from the US National Snow and Ice Data Center (NSIDC). The 36-year (1979–2014) dataset is parsimoniously and accurately described by our DAHMs. Preliminary results indicate that simulations using our multilayer Stuart-Landau model (MSLM) of SICs are stable for much longer time intervals, beyond the end of the twenty-first century, and exhibit interdecadal variability consistent with past historical records. Preliminary results indicate that this MSLM is quite skillful in predicting September sea ice extent. 

The multiscale variability of the ocean circulation due to its nonlinear dynamics remains a big challenge for theoretical understanding and practical ocean modeling. This paper demonstrates how the data-adaptive harmonic (DAH) decomposition and inverse stochastic modeling techniques introduced in (Chekroun and Kondrashov, (2017), Chaos, 27), allow for reproducing with high fidelity the main statistical properties of multiscale variability in a coarse-grained eddy-resolving ocean flow. This fully-data-driven approach relies on extraction of frequency-ranked time-dependent coefficients describing the evolution of spatio-temporal DAH modes (DAHMs) in the oceanic flow data. In turn, the time series of these coefficients are efficiently modeled by a family of low-order stochastic differential equations (SDEs) stacked per frequency, involving a fixed set of predictor functions and a small number of model coefficients. These SDEs take the form of stochastic oscillators, identified as multilayer Stuart–Landau models (MSLMs), and their use is justified by relying on the theory of Ruelle–Pollicott resonances. The good modeling skills shown by the resulting DAH-MSLM emulators demonstrates the feasibility of using a network of stochastic oscillators for the modeling of geophysical turbulence. In a certain sense, the original quasiperiodic Landau view of turbulence, with the amendment of the inclusion of stochasticity, may be well suited to describe turbulence. 

This chapter considers the sub-seasonal–to–seasonal (S2S) prediction problem as intrinsically more difficult than either short-range weather prediction or interannual–to–multidecadal climate prediction. The difficulty arises from the comparable importance of atmospheric initial states and of parameter values in determining the atmospheric evolution on the S2S time scale. The chapter relies on the theoretical framework of dynamical systems and the practical tools this framework helps provide to low-order modeling and prediction of S2S variability. The emphasis is on mid-latitude variability and the complementarity of the nonlinear-waves vs. multiple-regime points of view in understanding this variability. Empirical model reduction and the forecast skill of the models thus produced in real-time prediction are reviewed.

Watch the video online

Spectral analyses of the North Atlantic temperature field in the Simple Ocean Data Analysis (SODA) reanalysis identify prominent and statistically significant interannual oscillations along the Gulf Stream front and in large regions of the North Atlantic. A 7–8-yr oscillatory mode is characterized by a basin-wide southwest-to-northeast–oriented propagation pattern in the sea surface temperature (SST) field. This pattern is found to be linked to a seesaw in the meridional-dipole structure of the zonal wind stress forcing (TAUX). In the subpolar gyre, the SST and TAUX fields of this mode are shown to be in phase opposition, which suggests a cooling effect of the wind stress on the upper ocean layer. Over all, this mode’s temperature field is characterized by a strong equivalent-barotropic component, as shown by covariations in SSTs and sea surface heights, and by phase-coherent behavior of temperature layers at depth with the SST field. Recent improvements of multivariate singular spectrum analysis (M-SSA) help separate spatio-temporal patterns. This methodology is developed further and applied to studying the ocean’s response to variability in the atmospheric forcing. Statistical evidence is shown to exist for other mechanisms generating oceanic variability of similar 7–8-yr periodicity in the Gulf Stream region; the latter variability is likewise characterized by a strongly equivalent-barotropic component. Two other modes of biennial variability in the Gulf Stream region are also identified, and it is shown that interannual variability in this region cannot be explained by the ocean’s response to similar variability in the atmospheric forcing alone.