This contribution introduces the dynamics of shallow and rotating flows that characterizes large-scale motions of the atmosphere and oceans. It then focuses on an important aspect of climate dynamics on interannual and interdecadal scales, namely the wind-driven ocean circulation. Studying the variability of this circulation and slow changes therein is treated as an application of the theory of nonautonomous dynamical systems. The contribution concludes by discussing the relevance of these mathematical concepts and methods for the highly topical issues of climate change and climate sensitivity.

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.

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. 

Decline in the Arctic sea ice extent (SIE) is an area of active scientific research with profound socio-economic implications. Of particular interest are reliable methods for SIE forecasting on subseasonal time scales, in particular from early summer into fall, when sea ice coverage in the Arctic reaches its minimum. Here, we apply the recent data-adaptive harmonic (DAH) technique of Chekroun and Kondrashov, (2017), Chaos, 27 for the description, modeling and prediction of the Multisensor Analyzed Sea Ice Extent (MASIE, 2006–2016) data set. The DAH decomposition of MASIE identifies narrowband, spatio-temporal data-adaptive modes over four key Arctic regions. The time evolution of the DAH coefficients of these modes can be modelled and predicted by using a set of coupled Stuart–Landau stochastic differential equations that capture the modes’ frequencies and amplitude modulation in time. Retrospective forecasts show that our resulting multilayer Stuart–Landau model (MSLM) is quite skilful in predicting September SIE compared to year-to-year persistence; moreover, the DAH–MSLM approach provided accurate real-time prediction that was highly competitive for the 2016–2017 Sea Ice Outlook.

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. 

We study the influence of interannual climate variability on the economy of several countries in the Sahel region. In the agricultural sector, we are able to identify coupled climate-economic modes that are statistically significant on interannual time scales. In particular, precipitation is a key climatic factor for agriculture in this semi-arid region. Locality and diversity characterize the Sahel's climatic and economic system, with the coupled climate-economic patterns exhibiting substantial differences from country to country. Large-scale atmospheric patterns — like the El Niño–Southern Oscillation and its quasi-biennial and quasi-quadrennial oscillatory modes — have quite limited influence on the economies, while more location-specific rainfall patterns play an important role.

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.