Search

The energy surplus resulting from radiative forcing causes warming of the Earth system. This initial warming drives a myriad of changes including in sea surface temperatures (SSTs), leading to different radiative feedbacks. The relationship between the radiative feedbacks and the pattern of SST changes is referred to as the "pattern effect". The current approach to study the pattern effect relies on diagnosing the response of atmosphere-only models to perturbations in the SST boundary condition. Here, we argue that the fluctuation-dissipation relation (FDR), together with coarse-graining procedures, is a computationally cheap and theoretically grounded alternative to model experiments. We introduce a protocol to study the pattern effect and present its application in a state-of-the-art coupled climate model. By focusing on the coupled dynamics, we unveil the role of the slow ocean component in setting the pattern effect. We present a new "sensitivity map", representing a first, qualitative prediction of the response of the average top-of-the-atmosphere (TOA) radiative flux to perturbations in the SST field. We find negative sensitivity throughout the tropics, in contrast to the current understanding of a positive-negative dipole of sensitivity in the tropical Pacific. Considering only the shortest time scales, the response is dominated by the fast atmospheric variability and we recover results in qualitative agreement with the literature. Therefore, the difference between our results and previous studies, largely comes from including the atmosphere-ocean coupling. The framework offers a conceptually novel perspective on the pattern effect: feedbacks in the coupled system are encoded in a temporally and spatially dependent response operator, rather than time-independent maps as for previous studies.

We propose a data-driven framework to simplify the description of spatiotemporal climate variability into few entities and their causal linkages. Given a high-dimensional climate field, the methodology first reduces its dimensionality into a set of regionally constrained patterns. Time-dependent causal links are then inferred in the interventional sense through the fluctuation-response formalism, as shown in Baldovin et al. (2020). These two steps allow to explore how regional climate variability can influence remote locations. To distinguish between true and spurious responses, we propose a novel analytical null model for the fluctuation-dissipation relation, therefore allowing for uncertainty estimation at a given confidence level. Finally, we select a set of metrics to summarize the results, offering a useful and simplified approach to explore climate dynamics. We showcase the methodology on the monthly sea surface temperature field at global scale. We demonstrate the usefulness of the proposed framework by studying few individual links as well as "link maps", visualizing the cumulative degree of causation between a given region and the whole system. Finally, each pattern is ranked in terms of its "causal strength", quantifying its relative ability to influence the system's dynamics. We argue that the methodology allows to explore and characterize causal relationships in high-dimensional spatiotemporal fields in a rigorous and interpretable way.

Studies agree on a significant global mean sea level rise in the 20th century and its recent 21st century acceleration in the satellite record. At regional scale, the evolution of sea level probability distributions is often assumed to be dominated by changes in the mean. However, a quantification of changes in distributional shapes in a changing climate is currently missing. To this end, we propose a novel framework quantifying significant changes in probability distributions from time series data. The framework first quantifies linear trends in quantiles through quantile regression. Quantile slopes are then projected onto a set of four $orthogonal$ polynomials quantifying how such changes can be explained by $independent$ shifts in the first four statistical moments. The framework proposed is theoretically founded, general and can be applied to any climate observable with close-to-linear changes in distributions. We focus on observations and a coupled climate model (GFDL-CM4). In the historical period, trends in coastal daily sea level have been driven mainly by changes in the mean and can therefore be explained by a shift of the distribution with no change in shape. In the modeled world, robust changes in higher order moments emerge with increasing CO2 concentration. Such changes are driven in part by ocean circulation alone and get amplified by sea level pressure fluctuations, with possible consequences for sea level extremes attribution studies.

The threat of global warming and the demand for reliable climate predictions pose a formidable challenge being the climate system multiscale, high-dimensional and nonlinear. Spatiotemporal recurrences of the system hint to the presence of a low-dimensional manifold containing the high-dimensional climate trajectory that could make the problem more tractable. Here we argue that reproducing the geometrical and topological properties of the low-dimensional attractor should be a key target for models used in climate projections. In doing so, we propose a general data-driven framework to characterize the climate attractor and showcase it in the tropical Pacific ocean using a reanalysis as observational proxy and two state-of-the-art models. The analysis spans four variables simultaneously over the periods 1979-2019 and 2060-2100. At each time t, the system can be uniquely described by a state space vector parameterized by N variables and their spatial variability. The dynamics is confined on a manifold with dimension lower than the full state space that we characterize through manifold learning algorithms, both linear and nonlinear. The local geometry and local stability of the high-dimensional, multi-variable climate attractor are quantified through the local dimension and persistence metrics. Model biases that hamper climate predictability are identified and found to be similar in the multivariate attractor of the two models during the historical period while diverging under the warming scenario considered. Finally, the relationships between different sub-spaces (univariate fields), and therefore among climate variables, are evaluated. The proposed framework provides a comprehensive, physically based, test for assessing climate feedbacks and opens new avenues for improving their model representation.

One of the greatest sources of uncertainty in future climate projections comes from limitations in modelling clouds and in understanding how different cloud types interact with the climate system. A key first step in reducing this uncertainty is to accurately classify cloud types at high spatial and temporal resolution. In this paper, we introduce Cumulo, a benchmark dataset for training and evaluating global cloud classification models. It consists of one year of 1km resolution MODIS hyperspectral imagery merged with pixel-width 'tracks' of CloudSat cloud labels. Bringing these complementary datasets together is a crucial first step, enabling the Machine-Learning community to develop innovative new techniques which could greatly benefit the Climate community. To showcase Cumulo, we provide baseline performance analysis using an invertible flow generative model (IResNet), which further allows us to discover new sub-classes for a given cloud class by exploring the latent space. To compare methods, we introduce a set of evaluation criteria, to identify models that are not only accurate, but also physically-realistic. CUMULO can be download from https://www.dropbox.com/sh/i3s9q2v2jjyk2it/AACxXnXfMF5wuIqLXqH4NJOra?dl=0 .