The authors argue that the concept of snapshot attractors and of their natural probability distributions are the only available tools by means of which mathematically sound statements can be made about averages, variances, etc., for a given time instant in a changing climate. A basic advantage of the snapshot approach, which relies on the use of an ensemble, is that the natural distribution and thus any statistics based on it are independent of the particular ensemble used, provided it is initiated in the past earlier than a convergence time. To illustrate these concepts, a tutorial presentation is given within the framework of a low-order model in which the temperature contrast parameter over a hemisphere decreases linearly in time. Furthermore, the averages and variances obtained from the snapshot attractor approach are demonstrated to strongly differ from the traditional 30-yr temporal averages and variances taken along single realizations. The authors also claim that internal variability can be quantified by the natural distribution since it characterizes the chaotic motion represented by the snapshot attractor. This experience suggests that snapshot-attractor-based calculations might be appropriate to be evaluated in any large-scale climate model, and that the application of 30-yr temporal averages taken along single realizations should be complemented with this more appealing tool for the characterization of climate changes, which seems to be practically feasible with moderate ensemble sizes.

We study the teleconnection between El Niño–Southern Oscillation (ENSO) and the Indian summer monsoon (IM) in large ensemble simulations, the Max Planck Institute Earth System Model (MPI-ESM), and the Community Earth System Model (CESM1). We characterize ENSO by the June–August Niño-3 box-average SST and the IM by the June–September average precipitation over India, and define their teleconnection in a changing climate as an ensemble-wise correlation. To test robustness, we also consider somewhat different variables that can characterize ENSO and the IM. We utilize ensembles converged to the system’s snapshot attractor for analyzing possible changes in the teleconnection. Our main finding is that the teleconnection strength is typically increasing on the long term in view of appropriately revised ensemble-wise indices. Indices involving a more western part of the Pacific reveal, furthermore, a short-term but rather strong increase in strength followed by some decrease at the turn of the century. Using the station-based Southern Oscillation index (SOI) as opposed to area-based indices leads to the identification of somewhat more erratic trends, but the turn-of-the-century “bump” is well detectable with it. All this is in contrast, if not in contradiction, to the discussion in the literature of a weakening teleconnection in the late twentieth century. We show here that this discrepancy can be due to any of three reasons: 1) ensemble-wise and temporal correlation coefficients used in the literature are different quantities; 2) the temporal moving correlation has a high statistical variability but possibly also persistence; or 3) MPI-ESM does not represent the Earth system faithfully.

A time series resulting from a single initial condition is shown to be insufficient for quantifying the internal variability in a climate model, and thus one is unable to make meaningful climate projections based on it. The authors argue that the natural distribution, obtained from an ensemble of trajectories differing solely in their initial conditions, of the snapshot attractor corresponding to a particular forcing scenario should be determined in order to quantify internal variability and to characterize any instantaneous state of the system in the future. Furthermore, as a simple measure of internal variability of any particular variable of the model, the authors suggest using its instantaneous ensemble standard deviation. These points are illustrated with the intermediate-complexity climate model Planet Simulator forced by a CO2 scenario, with a 40-member ensemble. In particular, the leveling off of the time dependence of any ensemble average is shown to provide a much clearer indication of reaching a steady state than any property of single time series. Shifts in ensemble averages are indicative of climate changes. The dynamical character of such changes is illustrated by hysteresis-like curves obtained by plotting the ensemble average surface temperature versus the CO2 concentration. The internal variability is found to be the most pronounced on small geographical scales. The traditionally used 30-yr temporal averages are shown to be considerably different from the corresponding ensemble averages. Finally, the North Atlantic Oscillation (NAO) index, related to the teleconnection paradigm, is also investigated. It is found that the NAO time series strongly differs in any individual realization from each other and from the ensemble average, and climatic trends can be extracted only from the latter.

In an incompressible flow, fluid density remains invariant along fluid element trajectories. This implies that the spatial distribution of non-interacting noninertial particles in such flows cannot develop density inhomogeneities beyond those that are already introduced in the initial condition. However, in certain practical situations, density is measured or accumulated on (hyper-) surfaces of dimensionality lower than the full dimensionality of the flow in which the particles move. An example is the observation of particle distributions sedimented on the floor of the ocean. In such cases, even if the initial distribution of noninertial particles is uniform within a finite support in an incompressible flow, advection in the flow will give rise to inhomogeneities in the observed density. In this paper we analytically derive, in the framework of an initially homogeneous particle sheet sedimenting towards a bottom surface, the relationship between the geometry of the flow and the emerging distribution. From a physical point of view, we identify the two processes that generate inhomogeneities to be the stretching within the sheet, and the projection of the deformed sheet onto the target surface. We point out that an extreme form of inhomogeneity, caustics, can develop for sheets. We exemplify our geometrical results with simulations of particle advection in a simple kinematic flow, study the dependence on various parameters involved, and illustrate that the basic mechanisms work similarly if the initial (homogeneous) distribution occupies a more general region of finite extension rather than a sheet. a) gabor@ifisc.uib-csic.es 1 Sedimentation of small particles in complex flows is an outstanding problem in science and technology. In particular, the sinking of biogenic particles from the marine surface to the bottom is a fundamental process of the biological carbon pump, playing a key role in the global carbon cycle. A complete understanding of this problem is still lacking. It has been recently shown that despite fluid incompressibility, sedimenting particles moving as noninertial particles in the ocean on large scales show density inhomogeneities when accumulated on some bottom surface. Here, we analytically derive the relation between the geometry of the flow and the emerging distribution for an initially homogeneous sheet of particles. From a physical point of view, we identify the two processes that generate inhomogeneities to be the stretching within the sheet, and the projection of the deformed sheet onto the target surface. We point out conditions under which an extreme form of inhomogeneity, caustics, can develop.

We study the forced response of the teleconnection between the El Niño–Southern Oscillation (ENSO) and the Indian summer monsoon (IM) in the Max Planck Institute Grand Ensemble, a set of Earth system ensemble simulations under historical and Representative Concentration Pathway (RCP) forcing. The forced response of the teleconnection, or a characteristic of it, is defined as the time dependence of a correlation coefficient evaluated over the ensemble. We consider the temporal variability of spatial averages and that with respect to dominant spatial modes in the sense of Maximal Covariance Analysis, Canonical Correlation Analysis and Empirical Orthogonal Function analysis across the ensemble. A further representation of the teleconnection that we define here takes the point of view of the predictability of the spatiotemporal variability of the Indian summer monsoon. We find that the strengthening of the ENSO-IM teleconnection is robustly or consistently featured in view of various teleconnection representations, whether sea surface temperature (SST) or sea level pressure (SLP) is used to characterize ENSO, and both in the historical period and under the RCP8.5 forcing scenario. It is found to be associated dominantly with the principal mode of ENSO variability. Concerning representations that involve an autonomous characterisation of the Pacific, in terms of a linear regression model, the main contributor to the strengthening is the regression coefficient, which can outcompete even a declining ENSO variability when it is represented by SLP. We also find that the forced change of the teleconnection is typically nonlinear by 1) formally rejecting the hypothesis that ergodicity holds, i.e., that expected values of temporal correlation coefficients with respect to the ensemble equal the ensemble-wise correlation coefficient itself, and also showing that 2) the trivial contributions of the forced changes in means and standard deviations are insignificant here. We also provide, in terms of the test statistics, global maps of the degree of nonlinearity/nonergodicity of the forced change of the teleconnection between local precipitation and ENSO.

In nonautonomous dynamical systems, like in climate dynamics, an ensemble of trajectories initiated in the remote past defines a unique probability distribution, the natural measure of a snapshot attractor, for any instant of time, but this distribution typically changes in time. In cases with an aperiodic driving, temporal averages taken along a single trajectory would differ from the corresponding ensemble averages even in the infinite-time limit: ergodicity does not hold. It is worth considering this difference, which we call the nonergodic mismatch, by taking time windows of finite length for temporal averaging. We point out that the probability distribution of the nonergodic mismatch is qualitatively different in ergodic and nonergodic cases: its average is zero and typically nonzero, respectively. A main conclusion is that the difference of the average from zero, which we call the bias, is a useful measure of nonergodicity, for any window length. In contrast, the standard deviation of the nonergodic mismatch, which characterizes the spread between different realizations, exhibits a power-law decrease with increasing window length in both ergodic and nonergodic cases, and this implies that temporal and ensemble averages differ in dynamical systems with finite window lengths. It is the average modulus of the nonergodic mismatch, which we call the ergodicity deficit, that represents the expected deviation from fulfilling the equality of temporal and ensemble averages. As an important finding, we demonstrate that the ergodicity deficit cannot be reduced arbitrarily in nonergodic systems. We illustrate via a conceptual climate model that the nonergodic framework may be useful in Earth system dynamics, within which we propose the measure of nonergodicity, i.e., the bias, as an order-parameter-like quantifier of climate change.

Teleconnections are striking features of the Earth climate system which appear as statistically correlated climate-related patterns between remote geographical regions of the globe. In a changing climate, however, the strength of teleconnections might change, and an appropriate characterization of these correlations and their change (more appropriate than detrending the time series) is lacking in the literature. Here we present a novel approach, based on the theory of snapshot attractors, corresponding in our context to studying parallel climate realizations. Imagining an ensemble of parallel Earth systems, instead of the single one observed (i.e., the real Earth), the ensemble, after some time, characterizes the appropriate probabilities of all options permitted by the climate dynamics, reflecting the internal variability of the climate. We claim that the relevant quantities for characterizing teleconnections in a changing climate are correlation coefficients taken over the temporally evolving ensemble in any time instant. As a particular example, we consider the teleconnections of the North Atlantic Oscillation (NAO). In a numerical climate model, we demonstrate that this approach provides the only statistically correct characterization, in contrast to commonly used temporal correlations evaluated along single detrended time series. The teleconnections of the NAO are found to survive the climate change, but their strength might be time-dependent.

Based on the theory of "snapshot/pullback attractors", we show that important features of the climate change that we are observing can be understood by imagining many replicas of Earth that are not interacting with each other. Their climate systems evolve in parallel, but not in the same way, although they all obey the same physical laws, in harmony with the chaotic-like nature of the climate dynamics. These parallel climate realizations evolving in time can be considered as members of an ensemble. We argue that the contingency of our Earth's climate system is characterized by the multiplicity of parallel climate realizations rather than by the variability that we experience in a time series of our observed past. The natural measure of the snapshot attractor enables one to determine averages and other statistical quantifiers of the climate at any instant of time. In this paper, we review the basic idea for climate changes associated with monotonic drifts, and illustrate the large number of possible applications. Examples are given in a low-dimensional model and in numerical climate models of different complexity. We recall that systems undergoing climate change are not ergodic, hence temporal averages are generically not appropriate for the instantaneous characterization of the climate. In particular, teleconnections, i.e. correlated phenomena of remote geographical locations are properly characterized only by correlation coefficients evaluated with respect to the natural measure of a given time instant, and may also change in time. Physics experiments dealing with turbulent-like phenomena in a changing environment are also worth being interpreted in view of the attractor-based ensemble approach. The possibility of the splitting of the snapshot attractor to two branches, near points where the corresponding time-independent system undergoes bifurcation as a function of the changing parameter, is briefly mentioned. This can lead in certain climate-change scenarios to the coexistence of two distinct sub-ensembles representing dramatically different climatic options. The problem of pollutant spreading during climate change is also discussed in the framework of parallel climate realizations. Keywords Climate dynamics • Nonautonomous systems • Ensembles • Snapshot attractors • Natural measures Communicated by Valerio Lucarini.