# Random templex encodes topological tipping points in noise-driven chaotic dynamics

Gisela Charó, Michael Ghil, Denisse Sciamarella

► To cite this version:

Gisela Charó, Michael Ghil, Denisse Sciamarella. Random templex encodes topological tipping points in noise-driven chaotic dynamics. Chaos: An Interdisciplinary Journal of Nonlinear Science, 2023, 33 (10), 10.1063/5.0140660. hal-04764307

HAL Id: hal-04764307

https://hal.science/hal-04764307v1

Submitted on 9 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Random templex encodes topological tipping points in noise-driven chaotic dynamics

Special Collection: **Theory-informed and Data-driven Approaches to Advance Climate Sciences**

Gisela D. Charó ⚙️; Michael Ghil ⚙️; Denisse Sciamarella 🐦

<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

Check for updates

Chaos 33, 103141 (2023)

https://doi.org/10.1063/5.0140660

<figure><img src="image_5.png" /></figure>

View
Online

<figure><img src="image_6.png" /></figure>

Export
Citation

## Articles You May Be Interested In

Templex-based dynamical units for a taxonomy of chaos

*Chaos* (November 2024)

A templex-based study of the Atlantic Meridional Overturning Circulation dynamics in idealized chaotic models

*Chaos* (January 2025)

Templex: A bridge between homologies and templates for chaotic attractors

*Chaos* (August 2022)

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>
# Random templex encodes topological tipping points in noise-driven chaotic dynamics

Cite as: Chaos 33, 103141 (2023); doi: 10.1063/5.0140660

Submitted: 29 December 2022 · Accepted: 2 October 2023 ·

Published Online: 30 October 2023

<figure><img src="image_9.png" /></figure>

View Online

<figure><img src="image_10.png" /></figure>

Export Citation

<figure><img src="image_11.png" /></figure>

CrossMark

Gisela D. Charó,<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> ID Michael Chil,<fnref n="4" />,<fnref n="5" />,<fnref n="6" /> ID and Denisse Sciamarella<fnref n="1" />,<fnref n="3" />,<fnref n="7" />,⁸ ID

## AFFILIATIONS

**Note:** This article is part of the Focus Issue, Theory-informed and Data-driven Approaches to Advance Climate Sciences.

ªAuthor to whom correspondence should be addressed: denisse.sciamarella@cnrs.fr

## ABSTRACT

Random attractors are the time-evolving pullback attractors of deterministically chaotic and stochastically perturbed dynamical systems. These attractors have a structure that changes in time and that has been characterized recently using Branched Manifold Analysis through Homologies cell complexes and their homology groups. This description has been further improved for their deterministic counterparts by endowing the cell complex with a directed graph (digraph), which encodes the order in which the cells in the complex are visited by the flow in phase space. A templex is a mathematical object formed by a cell complex and a digraph; it provides a finer description of deterministically chaotic attractors and permits their accurate classification. In a deterministic framework, the digraph of the templex connects cells within a single complex for all time. Here, we introduce the stochastic version of a templex. In such a random templex, there is one complex per snapshot of the random attractor and the digraph connects the generators or “holes” of successive cell complexes. Tipping points appear in a random templex as drastic changes of its holes in time, through their birth, splitting, merging, or death. This paper introduces random complexes and computes them for the noise-driven Lorenz system’s random attractor.

© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). https://doi.org/10.1063/5.0140660

Branched manifolds underlying chaotic attractors have topological properties that remain invariant in a deterministic framework and that can be characterized using homologies.<fnref n="1" />,<fnref n="2" /> A more complete description is obtained if the cell complex whose homologies are computed is endowed with a directed graph (digraph) that prescribes cell connections in terms of the flow direction. Such a topological description is given by a templex, which carries the information of the structure of the branched manifold, as well as information on the flow.<fnref n="3" /> This work revisits the templex concept in a stochastic framework. Stochastic attractors in the pullback approach—like

the LOrenz Random Attractor (LORA)<fnref n="4" />—include sharp transitions in a Branched Manifold Analysis through Homologies (BraMAH).<fnref n="5" />⁻<fnref n="7" /> These sharp transitions can be suitably described using what we call here a random templex, computed from a sequence of BraMAH cell complexes and a digraph. The BraMAH cell complexes are such that changes can be followed in terms of how the generators of the homology groups, i.e., the “holes” of these complexes, evolve. The nodes of the digraph are the generators of the homology groups, and its directed edges indicate the correspondence between holes from one snapshot to the next. Topological tipping points

<footnote n="1">¹Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, C1428EGA, Buenos Aires, Argentina</footnote>

<footnote n="2">²CONICET — Universidad de Buenos Aires, Centro de Investigaciones del Mar y la Atmósfera (CIMA), C1428EGA, Buenos Aires, Argentina</footnote>

<footnote n="3">³CNRS — IRD — CONICET — UBA, Instituto Franco-Argentino para el Estudio del Clima y sus Impactos (IRL 3351 IFAECI), C1428EGA, Buenos Aires, France</footnote>

<footnote n="4">⁴Geosciences Department and Laboratoire de Météorologie Dynamique (CNRS and IPSL), École Normale Supérieure and PSL University, 75231 Paris Cedex 05, France</footnote>

<footnote n="5">⁵Department of Atmospheric & Oceanic Sciences, University of California, Los Angeles, California 90095-1565, USA</footnote>

<footnote n="6">⁶Departments of Mathematics and of Finance, Imperial College London, London SW7 2AZ, United Kingdom</footnote>

<footnote n="7">⁷Centre National de la Recherche Scientifique, 75794 Paris CEDEX 16, France</footnote>
can be identified with the creation, destruction, splitting, or
merging of holes, through a definition in terms of the nodes in
the digraph.

I. INTRODUCTION

The topological characterization of noise-driven chaos is a
challenging issue that is crucial in the understanding of complex
systems, where part of the dynamics remains unresolved and is mod-
eled as noise. While additive noise in a system of equations will blur
the topological structure, multiplicative noise may radically change
it, as shown by Charó et al.⁷ These authors extended the concept of
a branched manifold to account for the integer-dimensional set in
phase space that robustly supports the system's invariant measure at
each instant.

Such a branched manifold, however, does not contain—as does
its deterministic counterpart—any information about the future or
the past of the invariant measure. In other words, the evolution of
the system is not completely described by the branched manifold,
which is now itself time-dependent. The latter requires, therefore,
additional information for its complete description.

The templex was introduced in the realm of deterministic
attractors in order to provide more topological information than
that contained in a cell complex.³ This missing information con-
cerns the flow around the branched manifold. This information can
be spelled out using a digraph⁹ that connects the cells of the complex
according to the flow.

But what about the flow in a cell complex representing the
invariant measure of a random attractor? One could try to pose it
in terms of connections between cells of different cell complexes,
but algebraic topology definitions are such that the number and
distribution of cells in a cell complex are arbitrary. It is not the indi-
vidual cells, but the topological properties of each cell complex that
characterize the changes from one instant to the next. Such topo-
logical properties are encoded by the generators of the homology
groups of the cell complex and also by the torsion groups. Herein,
we concentrate on the homology groups exclusively and leave tor-
sions for future work. The properties of interest are independent
of the particular cell decomposition that is adopted to build the
cell complex from data. Homologies hence enable us to connect a
cell complex of a random attractor at a certain instant, with a cell
complex corresponding to a different instant.

Homological properties can be computed at distinct times, and
they can also be tracked across sufficiently small time intervals, thus
helping us detect sudden changes in the topology. The information
on the mapping of the holes at a certain instant to the holes at the
next time step can be encoded using a digraph. This naturally leads
to the definition of a “random templex” as the mathematical object
that condenses the topological information regarding the evolution
of the system’s invariant measure in a finite time window. Topolog-
ical tipping points (TTPs) are contained in a random templex in the
form of creation, splitting, merging, and destruction of generators
or, equivalently, in certain simple characteristics of the nodes of the
digraph associated with a sequence of cell complexes.

This paper provides the theoretical background that leads to the
concept of a random templex and shows how to compute it for the

Lorenz Random Attractor (LORA). A brief explanation concerning
templexes in a deterministic and autonomous framework appears in
Sec. II. The construction of a random templex in a nonautonomous
and stochastically perturbed framework is discussed in Sec. III and
illustrated by application to LORA. Finally, in Sec. IV, using the
LORA templex in Sec. III, we define rigorously the TTPs introduced
in a merely intuitive fashion by Charó et al.⁷ and compute them for
LORA in a time window that contains the various types of TTPs.
Section V contains conclusions and perspectives, including possible
applications to the effects of global change on climatic subsystems,
often referred to recently as tipping elements.⁹

Four appendixes clarify the concepts and associated method-
ologies of pullback attractors, in which we define pullback, as
opposed to forward, attraction (Appendix A); random attractors,
in which the time-dependent forcing is random, as opposed to
deterministic (Appendix B); sample measures and their evolution
in time (Appendix C); and computing topological features from a
templex, on the algebraic calculations involved in determining the
homologies and the flow on a random templex (Appendix D).

II. THE DETERMINISTIC FRAMEWORK

Homology theory allows us to classify manifolds in terms of a
small number of topological properties.¹⁰ In order to study the topo-
logical structure of solutions of deterministic dynamical systems,
the characterization must be done from a set of discretely sampled
trajectories, i.e., a set of points with as many coordinates as time-
dependent variables. In the case of the Lorenz¹¹ chaotic attractor, for
instance, the points have three coordinates; we refer hereafter to this
model as L63. The phase space is three dimensional, but the points
lie on a butterfly shaped surface. How can the topological properties
of this surface be computed from the set of points? This is achieved
by building a cell complex.

A. Cell complexes

A cell complex K is a layered structure formed by a set of *k*-cells with *k* = 0, 1, ..., *d*. Each cell stands for an Euclidean closed set with a certain dimension: points are 0-cells, segments are 1-cells, filled polygons are 2-cells, and so forth. The highest cell dimension defines the dimension of the complex.

There exist methods to build a complex from a point cloud that
are quite different. Our choice will be to build a BraMAH complex,
i.e., a cell complex whose cells are formed gathering subsets of points
which can be locally approximated by a *d*-disk, with *d* being the local
dimension of the underlying manifold. Further details of this proce-
dure can be found in *Sciamarella and Mindlin*.⁶ A BraMAH complex
for the Lorenz attractor is shown in Fig. 1. Its highest dimensional
cell is a 2-cell and therefore *d* = 2. The polygon figures forming the
2-cells “pave” the surface of the attractor. The number of cells used
to pave this surface is not fixed. It depends on the criteria used to
define the sets of points. But this number has no special signifi-
cance, because the homologies finally obtained will not depend on
the particular size or distribution of the cells.¹²

In order to make homologies computable, *k*-cells with *k* > 1
must have a “direction” or “orientation.” In a 1-cell, this direction
<figure><img src="image_12.png" /><figcaption>FIG. 1. Construction of a BraMAH complex for the Lorenz attractor.¹¹ The point cloud is decomposed into clusters of points that are used to form the 2-complex K, with the 2-cells {$γ_i : i = 1, ..., 14$}.</figcaption></figure>

can be denoted by an arrow. For instance, if (1, 2) stands for instructions to travel along the path from the 0-cell (1) to the 0-cell (2), the expression −(1, 2) has a natural meaning: it means traveling in the opposite direction, i.e., from (2) to (1). A direction or orientation can also be assigned to 2-cells: clockwise or counterclockwise. A 2-complex K is said to be directed if each 1-cell and 2-cell is assigned a direction.

A directed complex leads to an algebra of chains. A k-chain is a formal linear combination of the k-cells in a cell complex. The boundary of a directed 2-cell, for instance, is the chain formed by the 1-cells on its boundary, with a positive sign if the direction of the edge is consistent with the direction of the 2-cell, and with a negative sign otherwise. It is, thus, the linear combination of 1-cells forming its boundary, while using the signs in accordance with the arrows. All k-chains in a complex form an abelian group.¹²

The purpose of introducing these algebraic concepts is to come up with something that will distinguish the important characteristics of a given cell complex K. The chains of 1-cells forming loops that are not the boundaries of the 2-cells will be important, as well as a chain of 2-cells enclosing a cavity, as in a torus or a sphere. These features are summarized by the homology groups of K, whose generators identify the “holes” at level k of the cell complex, namely, $H_k = [h_1^k, ..., h_{\beta_k}^k]$, where $h_i^k$ is the ith k-hole and $\beta_k$ is the Betti number, which counts the number of k-holes. It was Poincaré who proved the invariance of these numbers for a given set, independently of the details of the construction of a cell complex for the set. Moreover, $\beta_k = 0$ for $k \ge n$, where n is the dimension of the space in which the set is embedded.¹⁰,¹³

The significance of such homology groups can be understood level by level. At level 0, we have $H_0(K)$, which contains the 0-holes that identify the disconnected pieces of the cell complex K. L63 has a single connected component and therefore $\beta_0 = 1$.

Stepping up to $k = 1$, $H_1(K)$ lists the 1-generators, which are associated with the non-trivial loops or 1-holes in the complex. When the manifold underlying a complex has “handles,” the 1-generators encircle these handles forming closed sequences of 1-cells that can be traveled through sequentially. In L63, there are two such handles, $\beta_1 = 2$, one in each wing of the butterfly, as shown in Fig. 2.

These generators need not strictly contour the boundaries of the geometrical hole. A generator may wander around a hole, without tightly encircling the empty space. But the tight holes encircling handles can be retrieved algebraically, from the 2-complex itself, as long as it is uniformly oriented, i.e., the orientation of the 2-cells is propagated so that shared borders (1-cells) of two adjacent 2-cells are canceled out when the borders of all the 2-cells of the complex K are summed.

This information is contained in what we call the orientability chain ($O_1$). It is defined as the 1-chain obtained by applying the border operator $\partial_2$ to the sum of all the 2-cells⁶ (see Appendix D). This yields a sum of 1-cells, which includes boundaries and torsions,

$$O_1 = \partial_2 \left( \sum_i \gamma_i \right) = \sum_j a_j \sigma_j,$$

where $\gamma_i$ denotes the 2-cells, $\sigma_j$ denotes the 1-cells, and the $a_j$ are integers. The 1-cells in $O_1$ whose $a_j = \pm 1$ are the boundaries of the complex. In the case we are dealing with, $O_1$ can be used to replace the generators by their homologically equivalent “tight” holes. These 1-holes will be called minimal holes.

A graphic example comparing generators and minimal holes for the BraMAH complex of the L63 attractor is given in Figs. 2(a) and 2(b). Comparison of the two panels shows that the 1-generator in the right wing is minimal, while the 1-generator in the left wing is not.

Passing on to the next layer, generators of $H_2$ identify empty volumes enclosed by a surface, here called 2-holes. This is achieved by looking for chains of 2-cells forming the set of polygons that enclose the cavity. As there are no cavities enclosed by polygons in the L63 attractor, $\beta_2 = 0$, i.e., there are no 2-holes.

## B. Digraphs and complexes

What does a cell complex lack in order to fully characterize an attractor in phase space? The cell complex describes the shape of the attractor, but not the dynamics of the flow on the attractor. This missing information can be brought into the description by indicating the order in which the cells are visited by the flow. This can be done using a directed graph defined so that its nodes represent the highest dimensional cells of the complex—the 2-cells in the case of a 2-complex.

Two nodes are linked by a directed edge if the flow connects the cells in a given order. In the directed graph for L63 shown in Fig. 3, there is an edge connecting $\gamma_7$ to $\gamma_8$ and $\gamma_{13}$ to $\gamma_8$ because trajectories can flow from the 2-cells $\gamma_7$ and $\gamma_{13}$ to the 2-cell $\gamma_8$.

This leads to the definition of a dual object formed by a cell complex K and its digraph G. Charó et al.³ introduced this novel type of object and called it a **templex**, a contraction of template + complex.
<figure><img src="image_13.png" /><figcaption>FIG. 2. Complex K of Fig. 1 showing in color (a) the generators of $H_1(K)$ or 1-holes obtained in the homology computation and (b) the minimal 1-holes.</figcaption></figure>

**Definition 1:** A 2-templex $T = (K, G)$, is a templex of dimension 2, where K is a 2-complex, whose underlying structure is a branched 2-manifold associated with a deterministic dynamical system, and $G = (N, E)$ is a digraph, such that (i) the nodes N are the 2-cells of K and (ii) the edges E are the connections between the 2-cells associated with the flow.

The properties of a templex that characterize both the shape of the branched manifold and the flow upon it can be derived using the combined properties of the cell complex and the digraph. The result is expressed in terms of what we call **stripexes**. They have the same role and inherit their name from strips in templates.¹⁴ Stripexes identify the nonequivalent ways of circulating around the attractor according to the flow on it. The L63 attractor has four stripexes, which can be computed looking for the cycles in the digraph and retaining only the nonredundant ones. For further details on how to do these computations, the reader is referred to Charó *et al*.³

The nonequivalent paths along the complex representing the L63 attractor are

$$ \gamma_1 \to \gamma_2 \to \gamma_4 \to \gamma_6 \to \gamma_1, \quad (1a) $$

$$ \gamma_8 \to \gamma_9 \to \gamma_{11} \to \gamma_{13} \to \gamma_8, \quad (1b) $$

$$ \gamma_1 \to \gamma_3 \to \gamma_5 \to \gamma_7 \to \gamma_8, \quad (1c) $$

<figure><img src="image_14.png" /><figcaption>FIG. 3. The digraph G for the complex K of L63 shown in Fig. 2. The digraph indicates the connections associated with the flow between the cells. Together, the complex and the digraph form the templex $T = (K, G)$.</figcaption></figure>

$$ \gamma_8 \to \gamma_{10} \to \gamma_{12} \to \gamma_{14} \to \gamma_1. \quad (1d) $$

The stripexes of $T = (K, G)$ are the four sub-templexes associated with the distinct paths apparent in Fig. 3. They are denoted by $S_i$ ($i = 1, 2, 3, 4$). The first two stripexes are paths along either wing of the butterfly, according to Eqs. (1a) and (1d). The last two are complementary paths going from one wing to the other, according to Eqs. (1b) and (1c).

In a deterministic framework, the topological structure of an attractor can be accurately described by a templex because the branched manifold is invariant. But this is not the case for a random attractor. As shown by Charó *et al.*,¹⁷ the concept of branched manifold and cell complex can still be of help in the stochastic context. Section III will discuss how the templex theory can be extended to deal with a time-varying topological structure.

### III. THE STOCHASTIC FRAMEWORK

In a stochastic framework, ensembles of trajectories driven by the same noise path $\omega_t$ can be tracked and provide considerable insight into the overall behavior of the system, as shown by Chekroun *et al.*,⁴,¹⁵ and by Tél and associates¹⁶,¹⁷ in the climate sciences and by a vast literature on random dynamical systems.¹⁸,¹⁹ This pullback approach cancels out the well-known smoothing effect of noise and makes the fractal structure of the noise-driven chaotic dynamics emerge.

When the L63 model¹¹ is perturbed by a multiplicative noise in the Itô sense,¹⁹ with $W_t$ a Wiener process and $\sigma > 0$ the noise intensity, we get the stochastic Lorenz model of Chekroun *et al.*,⁴

$$ dx = s(y - x)dt + \sigma x dW_t, \quad (2a) $$

$$ dy = (rx - y - xz)dt + \sigma y dW_t, \quad (2b) $$

$$ dz = (-bz + xy)dt + \sigma z dW_t. \quad (2c) $$

The model's three parameters take the standard values for deterministically chaotic behavior, namely, $(s, r, b) = (28, 10, 8/3)$.
In the pullback approach, one is interested in the time-
dependent sample measures $\mu_t$ (see Appendix C) that are driven by
the same noise realization $\omega$, when starting from any bounded com-
pact set of initial data at a time $s$ in the remote past and watching
them converge at times $t$ as $s \to -\infty$.

The point clouds that are obtained in this way show the time-dependent stretching and folding mechanisms caused by the noise-driven nonlinearities. These point clouds evolve, combining the smoothness of the L63 deterministic convection with sudden deformations of the pullback attractor’s support. In the stochastic setting, these pullback attractors are called random attractors<sup>18,19</sup> or snapshot attractors.<sup>17,20</sup> See Appendixes A and B for details.

The stochastic Lorenz system (2) can be modeled using an RDS $(\varphi, \theta_t)$; see Appendix B for additional details on RDSs. The driving system $\theta_t$ is defined using the Wiener shift, as in Caravallo, Łukaszewicz, and Real,<sup>21</sup> and modeled using normally distributed random numbers. At time $t = N\Delta t$, where $\Delta t = 0.005$ is the time step and $N$ is the number of time steps, $\theta_t$ is modeled using MATLAB's function `randn(N, 3)`. Solutions of this stochastic model generate data sets in the form of three-dimensional point clouds.

The structure of a random attractor is captured by its sample measure $\hat{\mu}_t$; see Appendix C. A numerical estimation $\hat{\mu}_t$ of the time-dependent measure $\mu_t$ at time $t$ can be computed by letting a large set of $N_0$ initial points “flow” in phase space from the remote past until time $t$, for a fixed noise realization $\omega$. In practice, the convergence of the $\mu_t$'s approximation is observed for a set of $N_0 = 10^8$ initial points, as in Chekroun et al.<sup>4</sup>

Each initial point is forward-integrated following the stochastic
Lorenz system with all trajectories starting at the same initial time.
Each point in the random attractor at instant *t* is associated with a
value of *μ*ᵢₜ that is obtained by averaging over a volume encompass-
ing that point: higher or lower *μ*ᵢₜ-values correspond to more or less
populated regions of the random attractor.

Topological characterization of a snapshot using Branched Manifold Analysis through Homologies (BraMAH) is not straightforward, as already discussed in Charó et al.⁷ In order to characterize the topology of the evolving structure that is observed in phase space, we assume that the point cloud’s most populated regions at time *t* can be fitted by a branched manifold associated with a threshold *μ*ᵢₜ for the sample measure. Threshold selection is crucial for a judicious fit of the branched manifold to the sequence of point clouds as the attractor evolves. The topological analysis is, therewith, threshold dependent. This threshold is kept fixed for all time, so as to compare the topologies from one time to the next and to trace their changes.

For illustration purposes, Fig. 4 shows a sequence of the pro-
jections of point clouds $C^{t_j}$ (in black) onto the (y, z)-plane and their
corresponding sample measures (in color) obtained from LORA
simulations for three different consecutive time instants $0 < t_1$
$< t_2$, using the noise amplitude $\sigma = 0.3$. Each snapshot is the
result of an integration that is illustrated using the three “pullback
tunnels”—blue, red, and green. The yellow tunnel represents the
evolution of the theoretical random attractor from real time $t = 0$
through $t = t_1$ and on to $t = t_2$, as illustrated by the discrete snap-
shots $\mathcal{A}(0, \omega)$, $\mathcal{A}(t_1, \omega)$, and $\mathcal{A}(t_2, \omega)$. As in Fig. A2 of Chekroun et

<figure><img src="image_15.png" /><figcaption>FIG. 4. Diagram showing the projections onto the (y, z)-plane of nine point clouds {$C_{s=0}^{t_j} : s_i = -\tau_2, -\tau_1, 0; j = 0, 1, 2$} (plots of the numerical estimation of the sample measure $\hat{\mu}_t$) obtained from LORA simulations for three different consecutive time instants $0 < t_1 < t_2$, and plotted as the blue, red, and green tunnels. The yellow tunnel represents the theoretical random attractor and displays the snapshots $\mathcal{A}(0, \omega)$, $\mathcal{A}(t_1, \omega)$, and $\mathcal{A}(t_2, \omega)$ at the time instants $t = 0, t_1, t_2$. The sieved point clouds $\bar{C}^0, \bar{C}^{t_1}$, and $\bar{C}^{t_2}$ (plots of colored holes) are obtained by applying a threshold to the point clouds $C_{s=0}^0, C_{s=0}^{t_1}, C_{s=0}^{t_2}$, respectively. We show the projections onto the (y, z)-planes of the sieved point clouds $\bar{C}^{t_j}$ with the 1-holes in colors chosen for mutual contrast. The analysis of the branched manifolds approximated by cell complexes in terms of homologies yields $\mathcal{H}_1(\bar{C}^0) \sim \mathbb{Z}^{21}$, $\mathcal{H}_1(\bar{C}^{t_1}) \sim \mathbb{Z}^{18}$, and $\mathcal{H}_1(\bar{C}^{t_2}) \sim \mathbb{Z}^{19}$, respectively.</figcaption></figure>
<figure><img src="image_16.png" /><figcaption>FIG. 5. Phase portrait showing the point cloud (in gray) juxtaposed on the cell complex, and the minimal 1-holes (in colors) for a single snapshot of the Lorenz Random Attractor (LORA) at $t = 40.09$. There are ten 1-holes and their barycenters are labeled with numbers.</figcaption></figure>

al.,⁴ reproduced here as Fig. 11 in Appendix C, the tunnels occur in the pullback time $s$, while the random attractor $A_r(\omega)$ evolves in the real (forward) time $t$.

Figure 2 of Charó et al.⁷ shows how the threshold $\bar{\mu}_t$, denoted there by $\bar{n}$, affects the analysis of the topology of the branched manifold: $\bar{n}$-values that are too high will tend to disassemble the structure while values that are too low will fill up all the holes. Note that the choice of $\bar{n}$ was done for each snapshot separately, while here it is the same for the entire set of point clouds that were treated.

Finally, we only retain here the points whose projection onto the $(y, z)$-plane is surrounded by more than $n = 100$ points within a grid box of size $\delta y \times \delta z = 0.0684 \times 0.0587$.

For the three sieved clouds in Fig. 4—$\bar{C}^0$, $\bar{C}^{t_1}$, and $\bar{C}^{t_2}$—the analysis of the branched manifolds approximated by cell complexes in terms of homologies yields $\mathcal{H}_1(\bar{C}^0) \sim \mathbb{Z}^{21}$, $\mathcal{H}_1(\bar{C}^{t_1}) \sim \mathbb{Z}^{18}$, and $\mathcal{H}_1(\bar{C}^{t_2}) \sim \mathbb{Z}^{19}$, respectively.

The pullback invariance of the random attractor $\mathcal{A}(t, \omega)$ and the pullback attraction of random sets $B$ is illustrated in Ghil et al.,²², Fig. A2 reproduced here as Fig. 11 in Appendix B. The difference between this figure and our Fig. 4 is the topological description shown through the colored point clouds $\bar{C}^{t_i}$ supporting the sample measure that approximates the most populated regions of the random attractor at the times $t = 0, t_1$, and $t_2$, respectively. For further details on constructing a BRAMAH complex for a snapshot from a point cloud, see Charó et al.⁷

Defining a complex for a random attractor involves establishing a link between different snapshots, in order to incorporate time into the description, as done in the deterministic case. But how can the topology of the point cloud at a given instant be related to the topology at another instant? We can compute the approximations $\hat{\mu}_t$ for several time instants and build cell complexes for each snapshot. Establishing a correspondence between the cells in subsequent

<figure><img src="image_17.png" /><figcaption>FIG. 6. Minimal 1-hole tracking from a snapshot (a) at $t = 40.07$ to the next one (b) at $t = 40.075$. Colors and numbers are arbitrary and specific to this figure. The correspondences are shown in (c) and will be used to construct the edges of the digraph $\mathcal{D}$ in LORA's random template $\mathcal{R}$.</figcaption></figure>

cell complexes is not trivial, since, as we have explained in Sec. II A, the number of cells and their distribution is arbitrary.

The cell complex as a whole flows in phase space for a fixed noise realization $\omega$, following the trajectories of the system (2). Due to the fact that the evolution of the deterministic convection is smooth, the minimal 1-holes in the complexes—which are intrinsic to the snapshot—can be followed from one snapshot to the next. This can be used to check if the holes are simply displaced or if they are modified drastically, through splitting, merging, creation, or destruction events. An example of minimal holes of a LORA point cloud is provided in Fig. 5 for $t = 40.09$ and $\sigma = 0.3$.
<figure><img src="image_18.png" /><figcaption>**FIG. 8.** A merging event. In the top panel, five holes that are distinct at *t* = 40.07 are shown to merge into a single dark red hole at *t* = 40.075. (Bottom panels) Phase portraits of the snapshots with the distinctly colored holes that merge into a single dark-red hole.</figcaption></figure>

occur at specific times, which are associated with TTPs, as mentioned in Sec. I. A random templex will be shown to encode such TTPs.

We limit our description here to the homology groups through
the minimal 1-holes, without discussion of orientability proper-
ties, which are associated with the torsion groups of a cell com-
plex. Neglecting torsion properties works here because the LORA
complexes do not present distinct torsion groups, which is the
case for many other known deterministic attractors, when changing
parameters and hence, supposedly, for their corresponding ran-
dom attractors. But we do not exclude that torsion groups might
be relevant for other random attractors or that twists as defined in
Charó et al.³ might also play a role.

Figure 6 shows how holes are tracked from one snapshot to the next using two successive snapshots of LORA, at $t = 40.07$ and $t = 40.075$. Tracking is performed by searching for the minimal distance between the barycenters of the minimal holes at consecutive snapshots. The symmetry of LORA is used in the tracking process. In Fig. 6(a), $b_1$ and $b_2$ are symmetric at time $t = 40.07$, and they must also be mapped to symmetric 1-holes $c_1$ and $c_2$ in the subsequent snapshot at $t = 40.075$, Fig. 6(b).

Applying this procedure to every hole in all the snapshots, we
obtain the digraph of the random templex. The result for ten LORA
snapshots in a time window $T_w = [40.065, 40.11]$ is shown in Fig. 7.
The digraph $\mathcal{D}$ here is not singly connected: $\mathcal{D}$ has 15 connected
components, and each of these directed subgraphs tells the story of
one or several holes.

Even if there are only ten snapshots within $T_w$, a single con-
nected component in $\mathcal{D}$ may have more than ten nodes because of

<figure><img src="image_19.png" /><figcaption>**FIG. 9.** A splitting event. (Top panel) A hole at *t* = 40.100 splits into two different holes with different locations in phase space at *t* = 40.105. (Bottom panels) Phase portraits of the snapshots with the dark-red-colored hole that splits into a blue and a fire red one.</figcaption></figure>

the existence of merging and splitting events that enable connections
between the storylines of different holes. Important LORA prop-
erties that can be extracted from its random templex are given in
Sec. IV.

IV. TOPOLOGICAL TIPPING POINTS (TTPs)

TTPs, introduced by Charó et al.,⁷ can be now identified and
classified using the digraph *D* of LORA’s random templex.

**Definition 3:** A *Topological Tipping Point (TTP)* occurring at time $t^* \in [t_1, t_s]$ and at position $\bar{x}^*$ is encoded in the digraph $\mathcal{D} = (N, E)$ of a random templex $\mathcal{R}$ either (i) as a node that receives or emanates two or more edges or else (ii) as an initial or terminal node of a connected component in $\mathcal{D}$ that does not correspond to time $t_1$ or $t_s$, respectively.

This definition allows one to classify TTPs. The nodes that
receive two or more edges are merging TTPs. The nodes that
emanate more than one edge are splitting TTPs. Merging can be fol-
lowed by splitting. The initial nodes of a connected component of
the digraph that do not correspond to time $t_1$ are creation events.
Destruction events are terminal nodes of a connected component
that do not correspond to time $t_s$.

Each tipping point is associated with a particular snapshot, that
is, a particular time instant in the time window $T_w$, and with a par-
ticular location in phase space, as defined by the barycenter of the
corresponding hole. Figures 8 and 9 show, in detail, an example of a
merging node and of a splitting node, respectively. In both figures,
the top plot shows the location of the nodes in phase space, as identi-
fied by the holes’ barycenters, with an indication of the time instant,
and the two bottom plots show the phase portraits of the snapshots
at which the holes merge or split being colored.
<figure><img src="image_20.png" /><figcaption>FIG. 10. The 15 constellations that represent the essence of the evolution of LORA's 1-holes within the time window $T_w$. The plot shows the embedding of the digraph $\mathcal{D}$ of the random templex $\mathcal{R}$ into the L63 model's phase space, by using the coordinates of the barycenters of the nodes $\mathcal{N}$.</figcaption></figure>

The digraph of LORA's random templex $\mathcal{R}$ is shown in Fig. 7. This tree corresponds to $s = 10$ snapshots within $t_1 \le t \le t_s$ with $t_1 = 40.065$ and $t_s = 40.11$. Applying the definition of TTP to the nodes and edges in $\mathcal{D}$, one can detect them and classify them according to the type of event. We find 18 splitting TTPs (in red), 12 merging TTPs (in blue), 2 mergings followed by splittings (in magenta), 8 creation TTPs (in green), and 12 destruction TTPs (in black).

A better picture of how the holes are evolving in the system's phase space can be gained by using the coordinates of the barycenters of the holes for an embedding of $\mathcal{D}$ into this space. We call the embeddings of $\mathcal{D}$'s distinct connected components into phase space **constellations**.

**Definition 4:** A constellation $\mathcal{C}$ is the set of immersed nodes and edges forming a connected component in the digraph $\mathcal{D}$ of a random templex. Each node is immersed in the phase space using the coordinates of the corresponding hole's barycenter.

Figure 10 presents the 15 constellations characterizing the evolution of LORA in $T_w$. While Fig. 7 is a purely graph-theoretical representation of the evolution of LORA's random templex, Fig. 10 is a step in connecting it with a more geometric representation of this evolution. Such a representation could provide a bridge between the random templex—whose simplicity benefits from the invariance of the topology and of its breakdowns—and a more detailed description of the flow's dynamics in phase space.

V. CONCLUDING REMARKS

In this paper, we have introduced novel tools into algebraic topology and used them to provide insights into the behavior of random attractors, which combine deterministic chaos with stochastic perturbations.<sup>18,19</sup> We summarize our results in Subsection V A and discuss them in Subsection V B.

A. Summary

Section I reviewed the basic tools of this trade: branched manifolds, cell complexes, and homology groups<sup>12</sup> and it laid out the road plan for the paper.

In Sec. II, we presented the by now fairly well-known theory for deterministic strange attractors with time-independent forcing. In Sec. II A, we illustrated the connection between the branched manifold that supports the Lorenz (L63) strange attractor's<sup>11</sup> invariant measure and a cell complex $K$, constructed by the Branched Manifold Analysis through Homologies (BraMAH) of Sciamarella and Mindlin;<sup>5</sup> see Fig. 1. We emphasized the independence of the homology groups $\mathcal{H}_k$ and associated Betti numbers $\beta_k$ from the details of the construction of the cell complex.<sup>10,12,13</sup> In particular, we showed that 1-generators, i.e., loops around a 1-hole with $\beta = 1$, need not be minimal to still provide the homological information; see Fig. 2.

In Sec. II B, we recalled the directed graph (digraph) $G$ associated with the cell complex $K$, as introduced by Charó et al.<sup>3</sup> for the L63 attractor, as well as for the spiral and funnel Rössler attractors, among others. This digraph $G$, with its nodes $N$ that are the cells and its edges $E$ that are the connections between them, is shown here for L63 in Fig. 3; it provides the direction of the flow on the branched manifold from one cell to another. Together, $K$ and $G$ form the templex $T = (K, G)$ associated with a particular nonlinear and, possibly, chaotic dynamics.

In Sec. III, we turned to our paper's main focus, namely, extending the recent digraph and templex concepts and tools reviewed in Sec. II to the study of random attractors that evolve in time.<sup>18-20</sup> To do so, we had to somehow incorporate the time into the definition of the digraph, by relating the cell complexes from one snapshot to another. This was implemented by constructing and labeling the minimal holes of one snapshot and tracking them to the next one; see Figs. 5 and 6.

More precisely, the correspondence between snapshots was established by identifying the 1-generators of the homology groups, i.e., the 1-holes of two consecutive complexes. In order to implement hole tracking between snapshots, the minimal holes were defined using the algebraic procedure described in the appendix. These holes are the nodes $\mathcal{N}$ used to encode, through the edges $\mathcal{E}$ of the digraph $\mathcal{D}$, how the random attractor's invariant measure evolves in time within a certain time window $T_w$.

The random templex $\mathcal{R}$ is hence defined as the couple $\mathcal{R} = (\mathcal{K}, \mathcal{D})$. The object $\mathcal{K}$ is the indexed family of cell complexes over the time interval $T_w$ under consideration, and the object $\mathcal{D}$ is the tree illustrated in Fig. 7, with nodes $\mathcal{N}$ at each snapshot and edges $\mathcal{E}$ connecting one snapshot with another.

Each node of the digraph is a minimal hole at a certain time. Each hole has the possibility of connecting with itself in the subsequent snapshot, but also with other holes, through mergings or splittings. Holes can disappear as well from a snapshot to the next or be created at an intermediate snapshot. Each connected component of the digraph $\mathcal{D}$ tells the story of how a hole evolves and how it connects to other holes, as time progresses; see Fig. 7. Examples of merging and splitting of minimal 1-holes were given in Figs. 8 and 9, respectively.

Much of the motivation of the work herein had to do with the visually striking, rapid changes in the evolution of LORA's
invariant measure as described by Chekroun et al.⁷ These visually rapid changes in time were documented by changes in LORA's algebraic topology; see Fig. 4 in Charó et al.,⁷ in which $H_1$ changed from $Z^2$ to $Z^{10}$ and back to $Z^4$ in very small time steps of $\Delta t = 0.09$.

We thus addressed in Sec. IV, using the more sophisticated topological apparatus of Sec. III, the existence and nature of topological tipping points (TTPs) conjectured by Charó et al.⁷ A TTP occurring at time $t^*$ and at position $\bar{x}^*$ is encoded in the digraph $\mathcal{D} = (N, E)$ of a random complex $\mathcal{R}$ either (i) as a node that receives or emanates two or more edges or else (ii) as an initial or terminal node of a connected component in $\mathcal{D}$ that does not correspond to the beginning or end of the time interval $T_w$ under consideration.

First, by considering time steps of $\Delta t = 0.005$ that are even smaller compared to the L63 model's¹¹ characteristic time of order unity, we confirmed, at least numerically, that the holes actually appear, disappear, and merge or split instantaneously. While a truly rigorous mathematical proof that this is indeed so might be difficult, the numerical evidence is quite substantial. Another matter that is left open at this point is whether changes in torsion—such as between a Möbius band, with torsion, and a regular one, without it—might be as sudden or not.

## B. Discussion

Given the fairly novel and even surprising nature of the concepts and tools introduced herein, a number of other questions are worth mentioning. One concerns the possibility of establishing closer connections between the metric flow of the solutions of a dynamical system in phase space and the topological description provided herein. The constellations shown in Fig. 10 and discussed at the end of Sec. IV suggest such a possibility, since one uses the embedding of the digraph $\mathcal{D}$ of the random complex $\mathcal{R}$ into the L63 model's phase space, by using the coordinates of the barycenters of the nodes $\mathcal{N}$.

This representation emphasizes how a random complex contains the information of when—i.e., in which snapshot—and where, i.e., in which phase space location—topological tipping is taking place. Charó et al.⁷ emphasized already that BraMAH-based topological data analysis³,⁶ is not restricted to low-dimensional dynamical systems. Clearly, some forms of reduced-order modeling²³,²⁴ will be necessary to treat models with substantial geographical resolution. But it is at least conceivable to use the mixed localization approach of our constellations for the study of tipping in subsystems of a larger system, in the spirit of the rather popular “tipping elements” of Lenton et al.⁹

## ACKNOWLEDGMENTS

It is a pleasure to acknowledge stimulating discussions with M.D. Chekroun on extending the results of Charó et al.⁷ and improving the detection of topological tipping points. Two anonymous reviewers and one anonymous editor have made insightful comments that helped clarify certain items in the original version of this paper. This work has received funding from the ANR project TeMPlex ANR-23-CE56-0002 (D.S.). G.D.C. gratefully acknowledges her postdoctoral scholarship from CONICET. The present paper is TiPES Contribution No. 210; this project has

received funding from the EU Horizon 2020 Research and Innovation Program under Grant Agreement No. 820970, and it helps support the work of M.G.

## AUTHOR DECLARATIONS

### Conflict of Interest

The authors have no conflicts to disclose.

### Author Contributions

**Gisela D. Charó:** Conceptualization (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Software (equal); Validation (equal); Visualization (equal); Writing – original draft (equal); Writing – review & editing (equal). **Michael Ghil:** Conceptualization (equal); Funding acquisition (equal); Investigation (equal); Resources (equal); Supervision (equal); Writing – original draft (equal); Writing – review & editing (equal). **Denisse Sciamarella:** Conceptualization (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Software (equal); Supervision (equal); Validation (equal); Visualization (equal); Writing – original draft (equal); Writing – review & editing (equal).

## DATA AVAILABILITY

The data that support the findings of this study are available within the article. The code that computes complex properties is available at git.cima.fcen.uba.ar/sciamarella/complex-properties.git.

## APPENDIX A: PULLBACK ATTRACTORS (PBAS)

We only consider here the finite-dimensional case of nonautonomous dynamical systems (NDS) theory, as illustrated, for instance, by a differential equation of the form

$$ \dot{x} = F(t, x), \quad t \in \mathbb{R}, \qquad (A1) $$

where $F$ denotes a smooth time-dependent vector field that governs the time evolution of the state $x$ in a phase space $X$.²⁵,²⁶ Let us denote by $\Phi(t, s)x_0$ the solution map of (A1) at time $t$, where $x_0$ is the initial state at time $s \le t$. The operator $\Phi(t, s)$ provides a two-time description of the system's evolution: the time $s$ when the system was initialized and the time $t \ge s$ of the system's current state. One speaks of a two-parameter semigroup with parameters $s$ and $t$.

In the autonomous case in which the right-hand side of Eq. (A1) does not depend explicitly on $t$, $F = F(x)$, and a one-parameter semigroup suffices to determine the evolution of the system. In this case, only the time interval $t' = t - s$ matters and $\Phi(t, s)$ reduces to $\Phi(t')$, i.e., to a standard flow. Hence, forward-in-time attraction only is of interest, resulting in attractors that are fixed, time-independent objects such as fixed points, limit cycles, tori, and strange attractors. In the case of nonautonomous systems, such as the climate system,²⁷,²⁸ there is explicit time dependence of forcing or coefficients, $F = F(t, x)$, and we need to define pullback attraction and pullback attractors (PBAs) or snapshot attractors.⁴,¹⁹,²¹,²²,²⁹,³⁰

**Definition 5:** An indexed family of objects {$\mathcal{A}(t)}_{t \in \mathbb{R}}$ is a PBA iff satisfies the following three fundamental properties:
(i) Each snapshot $A(t)$ is a compact subset of $X$.

(ii) **Invariance:** $\{A(t)\}_{t \in \mathbb{R}}$ is invariant with respect to the dynamics: $\Phi(t, s)A(s) = A(t)$ for all $s \le t$.

(iii) **Pullback attraction:** For all bounded subsets $B$ of $X$,

$$ \lim_{s \to -\infty} \mathrm{dist}(\Phi(t, s)B, A(t)) = 0, $$

for any fixed $t \in \mathbb{R}$

In simple words, the family $\{A(t)\}_{t \in \mathbb{R}}$ is invariant under the system's dynamics and it attracts at each time $t$ in the pullback sense all compact initial subsets $B$ from the remote past.

## APPENDIX B: RANDOM ATTRACTORS

In natural systems, such as the climate system, random time-dependent forcing is present²⁸ and it is necessary to use stochastic differential equations (SDEs) to model the system.²⁸,³¹ In the theory of random dynamical systems (RDSs), when the time-dependent forcing is random, PBAs are known as random attractors.¹⁹,²⁵,²⁶

In this framework, to have a reasonably self-consistent description of the system's random aspects, the noise has to be modeled as a time-dependent stochastic process in a probability space $(\Omega, \mathcal{F}, \mathbb{P})$.⁴ In this probability space, $\Omega$ is the sample space; $\mathcal{F}$ is a $\sigma$-algebra of measurable subsets of $\Omega$; and $\mathbb{P}$ is a probability measure on $\mathcal{F}$.¹⁹ To track the noise effect in time, one defines a time-dependent shift $\theta_t: \Omega \to \Omega$:

**Definition 6:** Let $(\Omega, \mathcal{F}, \mathbb{P})$ be a probability space and a time-dependent mapping $\theta_t: \Omega \to \Omega$ with $t \in \mathbb{R}$ a flow on $\Omega$. The mapping $\theta_t$ satisfies

(i) $\theta_0\omega = \omega, \forall \omega \in \Omega$;

(ii) $\theta_{t+\tau}\omega = \theta_t \circ \theta_\tau\omega, \forall t, \tau \in \mathbb{R}$;

(iii) $\theta_t \mathbb{P} = \mathbb{P} \forall t \in \mathbb{R}$.

The mapping $\theta_t$ connects each realization $\omega$ at time 0 to its state after a time $t$ has elapsed and the condition (iii) means that the probability measure $\mathbb{P}$ is preserved by $\theta_t$. Then $(\Omega, \mathcal{F}, \mathbb{P}, \theta)$ is called a driving dynamical system.

In the case of a stochastic-dynamic system driven by a Wiener process $W$, $\theta_t$ is defined on $\Omega$ according to $W_s(\theta_t\omega) = W_{t+s}(\omega) - W_s(\omega)$.¹⁹ To define an RDS, one can use $\theta_t$, to let the noise realization $\omega$ evolve in time, and define a cocycle mapping $\varphi$ to describe the evolution of the state $\mathbf{x}$.

**Definition 7:** Let $(\Omega, \mathcal{F}, \mathbb{P})$ be a probability space. An RDS $(\varphi, \theta_t)$ on $X$ consists on a mapping $\varphi: \mathbb{R} \times \Omega \times X \to X$ with the cocycle property and $\theta_t$ a driving dynamical system acting on the probability space. The RDS $(\varphi, \theta_t)$ satisfies²¹

(i) $\varphi(0, \omega)x = x, \forall x \in X$ and $\omega \in \Omega$, and

(ii) $\varphi(t+s, \omega)x = \varphi(t, \theta_s(\omega)) \circ \varphi(s, \omega)x, \forall s, t \in \mathbb{R}$ and $\omega \in \Omega$.

The above equalities set (i) the initial condition for the cocycle; and (ii) the cocycle property. This cocycle property is satisfied for a broad class of SDEs like those of interest here; see Arnold.¹⁹ In fact, Fig. 4 and its discussion in Sec. III of the main text herein illustrate very well the difference between (a) the one-parameter semigroups of autonomous dynamics and their purely forward attractors; and (b) the two-parameter semigroups of nonautonomous dynamics and their attractors, which can be either forward or pullback.

<figure><img src="image_21.png" /><figcaption>**FIG. 11.** Schematic diagram showing the pullback invariance of the random attractor $A(\omega)$ and pullback attraction of random sets $B$. Taken from Ref. 15.</figcaption></figure>

The evolution of a stochastic dynamical system can be modeled using an RDS $(\varphi, \theta)$, and its random attractor $\{A(t; \omega)\}_{t \in \mathbb{R}}$ provides the natural extension of a PBA to the random framework.

**Definition 8:** Let $(\varphi, \theta)$ be an RDS on $X$ and $\mathcal{B}$ a family of random sets. A random attractor $\{A(t; \omega)\}_{t \in \mathbb{R}}$ satisfies two properties:

(1) **Invariance:** $\varphi(t, \omega)A(\omega) = A(\theta_t\omega) \forall t \ge 0$.

(2) **Pullback attraction:** For all subsets $B \in \mathcal{B}$,

$$ \lim_{s \to -\infty} \mathrm{dist}(\varphi(t, \theta_s\omega)B, A(\omega)) = 0. $$

The time-indexed family of random compact sets $\{A(t; \omega)\}_{t \in \mathbb{R}}$ provides a complete description of all the system's possible states that are likely to be observed at any particular time $t$. Figure 11 shows a schematic diagram of a random attractor and the properties of pullback invariance and pullback attraction of random sets. The latter properties are further emphasized in Fig. 4 of Sec. III herein.

Note that, in practice, one does not need to let $s \to -\infty$ to compute a good numerical approximation of a random attractor. In the case of LORA⁴,⁷ that we treated here, the convergence time of the sample measure estimate $\hat{\mu}_t$ seems to be of the order of $s^* \simeq 40$, in the nondimensional units of the deterministic Lorenz convection model.¹¹ As far as our numerical results show, this convergence is uniform in $t$.

## APPENDIX C: SAMPLE MEASURES

A fundamental object in any RDS $(\theta, \varphi)$ is its random invariant measure $\mu$.

**Definition 9:** Let $(\Omega, \mathcal{F}, \mathbb{P})$ be the probability space of Appendix B, $\mathcal{B}$ its $\sigma$-algebra of Borel sets in $X$, and consider the RDS $(\varphi, \theta)$. A probability measure $\mu$ on $(\Omega \times X, \mathcal{F} \otimes \mathcal{B})$ is an invariant measure for the RDS $\varphi$ if it satisfies these two conditions:¹⁹

(i) $\theta_t\mu = \mu, \forall t \in \mathbb{R}$; and

(ii) for $\pi_{\Omega}(\omega, \mathbf{x}) = \omega$ being the projection of $(\omega, \mathbf{x})$ onto $\Omega$, the marginal $\pi_{\Omega}\mu$ of $\mu$ is $\theta$-invariant, i.e.,

$$ \pi_{\Omega}\mu = \theta_t(\pi_{\Omega}\mu) = \mathbb{P}, \forall t \in \mathbb{R}. $$
The property (ii) means that $\mu$ “lifts” $\mathbb{P}$ into the extended phase space $\Omega \times X$.⁴ Any probability measure $\mu$ possesses a disintegration or factorization,¹⁹ also known as a sample measure $\mu_{\omega}$.⁴,³² Sample measures enable one to work in the phase space X, rather than in the extended phase space $\Omega \times X$.

**Definition 10:** A sample measure $\mu_{\omega} : \Omega \times \mathcal{B} \to [0, 1]$ is such that

(i) For any Borel set $B \in \mathcal{B}$, $\mu_{\omega}$ is $\mathcal{F}$-measurable;

(ii) $\exists \Omega' \in \Omega$ such that $\mathbb{P}(\Omega') = 1$ and $\mu_{\omega}$ is a probability measure on $(X, \mathcal{B})$ for all $\omega$ in $\Omega'$; and

(iii) $\forall f \text{ in } L_{\mu}^{1}(\Omega \times X) \text{ we have}$

$$ \int_{\Omega \times X} f(\omega, x) \mu(d\omega, dx) = \int_{\Omega} \left( \int_{X} f(\omega, x) \mu_{\omega}(dx) \right) \mathbb{P}(d\omega). \quad (C1) $$

The sample measure $\mu_{\theta_{t}\omega}$ supported by the random attractor $\mathcal{A}(t; \omega)$ at time $t$ and for the noise realization $\omega$ provides the spatiotemporal probability distributions of the parts of the phase space $X$ occupied by the RDS. The sample measure $\mu_{\omega}$ is unique $\mathbb{P}$-almost surely,³⁵ and the invariance property (i) of Definition 9 translates into $\varphi(t, \omega)\mu_{\omega} = \mu_{\theta_{t}\omega}$.

When an RDS $(\theta, \varphi)$ has a random attractor $\mathcal{A}(t; \omega)$, then $\mu_{\omega}(\mathcal{A}(\omega)) = 1$ for almost all $\omega$ in $\Omega$, i.e., it supports every invariant measure. For a given noise realization $\omega$, the random attractor determines the parts of the phase space $X$ onto which any bounded random set $B$ is mapped at time $t$, when $B$ is propagated by $\varphi$ from a sufficiently remote past. Physical measures are sample measures that are “naturally chosen” by the dynamics.⁴,²⁸

**Definition 11:** A physical measure $\mu$ is a probability measure which satisfies: For any continuous observable $\psi: X \to \mathbb{R}$, the time average equals the ensemble average for almost all initial data $x_0$ that lie in a Lebesgue-positive set $B_{\mu}$, called the basin of attraction of $\mu$.⁴

An example of a class of probability measures that also coincides with that of physical measures for many dynamical systems are the Sinai–Ruelle–Bowen (SRB) measures.³³⁻³⁵

**Definition 12:** A probability measure $\mu$ is an SRB measure if its conditional measures on the system's unstable manifolds are absolutely continuous with respect to Lebesgue measure.

Chekroun et al.⁴ have shown in their Appendix C that, when the Kolmogorov operator associated with their Lorenz stochastic model (2), is hypoelliptic³⁶ and its leading Lyapunov exponent is positive, LORA supports a random SRB measure. They also provided numerical evidence that the measure to which the LORA's point clouds converge must be physical, too. But they gave no indication of the uniqueness of LORA's random SRB measure or of its physical measure.

## APPENDIX D: COMPUTING TOPOLOGICAL FEATURES FROM A TEMPLEX

This appendix details the algebraic computations that enable one to extract the topological features—such as homologies, joining loci, and stripexes—from a templex. For the sake of simplicity, we consider the chaotic attractor produced by the Rössler system³⁷ with

<figure><img src="image_22.png" /><figcaption>FIG. 12. Chaotic attractor of the Rössler system (D1). The parameter values are $a = 0.43295$, $b = 2$, and $c = 4$.</figcaption></figure>

$$ a = 0.43295, b = 2 \text{ and } c = 4. $$

$$ \dot{x} = -y - z, \qquad (\text{D1a}) $$

$$ \dot{y} = x + ay, \qquad (\text{D1b}) $$

$$ \dot{z} = b + z(x - c). \qquad (\text{D1c}) $$

A point cloud of this spiral attractor appears in Fig. 12.

A cell complex $K'$ for this system is shown in Fig. 13(a) using a planar diagram of the complex, i.e., a diagram in which the cells that appear to be duplicated are in fact two copies of the same cell and must therefore be glued together and considered as one. The advantage of such a planar diagram over the attractor's three-dimensional representation is that one can see the whole structure, whereas if we plot the complex in three dimensions—juxtaposed on the point cloud as in Fig. 1—some parts of the plot are obscured by the perspective. In the case of $K'$, the 1-cells that must be glued together are drawn with heavy lines. If this planar diagram is drawn on paper, one can obtain a model of the attractor by gluing the heavy lines together.

The complex $K'$ is made up of six 2-cells,

$$ \gamma_1 = \langle 1, 3, 2, 0 \rangle, \gamma_2 = \langle 2, 3, 6, 5, 4 \rangle, \gamma_3 = \langle 4, 5, 8, 7 \rangle, $$

$$ \gamma_4 = \langle 5, 6, 9, 8 \rangle, \gamma_5 = \langle 7, 8, 1, 0 \rangle, \gamma_6 = \langle 1, 8, 9, 0 \rangle. $$

The 2-cell denoted by $\gamma_6$ is attached through the 1-cell $\langle 0, 1 \rangle$ to the cells $\gamma_1$ and $\gamma_5$ with a gluing direction that highlights the folding that is taking place.

The boundary operator $\partial_k$ computes the boundary of an oriented k-cell $\partial_k: C_k \to C_{k-1}$, with $C_k = \sum_i a_i \sigma_k^i$, where $a_i \in \mathbb{Z}$ and $\sigma_k^i$ is a k-cell indexed by $i \in \mathbb{N}$. The action of the boundary maps is
<figure><img src="image_23.png" /><figcaption>FIG. 13. A complex $T' = (K', G')$ for the Rössler attractor. (a) Planar diagram of the oriented complex $K'$. (b) Digraph $G'$ indicating the connections between the 2-cells in $K'$.</figcaption></figure>

the following:

$$ \begin{aligned} \partial_0(1) &= 0, & \partial_1(0, 1) &= (1) - (0), \\ \partial_2(1, 8, 9, 0) &= (1, 8) + (8, 9) - (0, 9) + (0, 1). \end{aligned} $$

The information about the boundary of all the cells can be condensed within two matrices: $M_1$ in Table I represents the boundary of the 1-cells, i.e., it relates the 1-cells with the 0-cells, while $M_2$ in Table II represents the boundary of the 2-cells that relates the 2-cells with the 1-cells.

Using these matrices, it is possible to compute the quotient group

$$ \mathcal{H}_k := \mathcal{Z}_k / \mathcal{B}_k, $$

called the k-th homology group of a complex K. For instance, $\mathcal{H}_1 := \mathcal{Z}_1 / \mathcal{B}_1$, where $\mathcal{B}_1$ can be computed as the linearly independent rows of the transpose of $M_2$, and $\mathcal{Z}_1$, the set of 1-cycles, is the null space of the transpose of $M_1$. If the elements of $\mathcal{B}_1$ are expressed as linear combinations of the elements of $\mathcal{Z}_1$, the result is the list of chains of 1-cycles that are homologous to zero. This yields the homology relations between the 1-cycles, which lead, in turn, to the 1-generators, namely the generators of $\mathcal{Z}_1$.

In the case of $K'$, there is only one 1-generator,

$$ h_1 = +\langle 0, 2 \rangle + \langle 2, 4 \rangle + \langle 4, 7 \rangle - \langle 0, 7 \rangle. $$

<table><caption>TABLE I. $M_1$ is the matrix of $\partial_1$, which represents the boundary of the 1-cells of $K'$.</caption><thead><tr><th></th><th>∂<sub>1</sub></th><th>⟨0⟩</th><th>⟨1⟩</th><th>⟨2⟩</th><th>⟨3⟩</th><th>⟨4⟩</th><th>⟨5⟩</th><th>⟨6⟩</th><th>⟨7⟩</th><th>⟨8⟩</th><th>⟨9⟩</th></tr></thead><tbody><tr><th>M<sub>1</sub><sup>T</sup> =</th><td>⟨0, 1⟩</td><td>-1</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨0, 2⟩</td><td>-1</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨0, 7⟩</td><td>-1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨0, 9⟩</td><td>-1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><th></th><td>⟨1, 3⟩</td><td>0</td><td>-1</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨1, 8⟩</td><td>0</td><td>-1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td></tr><tr><th></th><td>⟨2, 3⟩</td><td>0</td><td>0</td><td>-1</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨2, 4⟩</td><td>0</td><td>0</td><td>-1</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨3, 6⟩</td><td>0</td><td>0</td><td>0</td><td>-1</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨4, 5⟩</td><td>0</td><td>0</td><td>0</td><td>-1</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨4, 7⟩</td><td>0</td><td>0</td><td>0</td><td>-1</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨5, 6⟩</td><td>0</td><td>0</td><td>0</td><td>0</td><td>-1</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨5, 8⟩</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>-1</td><td>0</td><td>0</td><td>1</td><td>0</td></tr><tr><th></th><td>⟨6, 9⟩</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>-1</td><td>0</td><td>0</td><td>1</td></tr><tr><th></th><td>⟨7, 8⟩</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>-1</td><td>1</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨8, 9⟩</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>-1</td><td>1</td><td>0</td></tr></tbody></table>

Notice that this 1-hole is minimal, since it contours the hole corresponding to the focus-type hole in the attractor. In short, all 1-cycles in $K'$ are homologous to $h_1$. For $K'$, the lower level, $\mathcal{H}_0$, has one generator since the complex has only one connected component. The higher level, $\mathcal{H}_2$, equals the empty set since there are no enclosed cavities. Furthermore, in this complex, like in $K$ for L63, there are no orientability chains and therefore no torsion group.

Up to this point, we have calculated homologies as customary in algebraic topology textbooks.¹² But the cell complex contains information relevant to the characterization of a chaotic attractor that is neither included in the homology groups nor in the torsion groups per se. Locating the 1-cells which are shared by at least three 2-cells, we can also compute the joining 1-cells. In $K'$, $\langle 0, 1 \rangle$ is such

<table><caption>TABLE II. $M_2$ is the matrix of $\partial_2$, which represents the boundary of the 2-cells of $K'$.</caption><thead><tr><th></th><th>∂<sub>2</sub></th><th>γ<sub>1</sub></th><th>γ<sub>2</sub></th><th>γ<sub>3</sub></th><th>γ<sub>4</sub></th><th>γ<sub>5</sub></th><th>γ<sub>6</sub></th></tr></thead><tbody><tr><th>M<sub>2</sub> =</th><td>⟨0, 1⟩</td><td>1</td><td>0</td><td>0</td><td>0</td><td>-1</td><td>1</td></tr><tr><th></th><td>⟨0, 2⟩</td><td>-1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨0, 7⟩</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td></tr><tr><th></th><td>⟨0, 9⟩</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>-1</td></tr><tr><th></th><td>⟨1, 3⟩</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨1, 8⟩</td><td>0</td><td>0</td><td>0</td><td>0</td><td>-1</td><td>1</td></tr><tr><th></th><td>⟨2, 3⟩</td><td>-1</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨2, 4⟩</td><td>0</td><td>-1</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨3, 6⟩</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨4, 5⟩</td><td>0</td><td>-1</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨4, 7⟩</td><td>0</td><td>0</td><td>-1</td><td>0</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨5, 6⟩</td><td>0</td><td>-1</td><td>0</td><td>1</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨5, 8⟩</td><td>0</td><td>0</td><td>1</td><td>-1</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨6, 9⟩</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td></tr><tr><th></th><td>⟨7, 8⟩</td><td>0</td><td>0</td><td>-1</td><td>0</td><td>1</td><td>0</td></tr><tr><th></th><td>⟨8, 9⟩</td><td>0</td><td>0</td><td>0</td><td>-1</td><td>0</td><td>1</td></tr></tbody></table>
a joining 1-cell. The 2-cells that are glued to it can also be identified: they are $\gamma_1$, $\gamma_5$ and $\gamma_6$. Together they form what we call the joining 2-chain: $\gamma_1 - \gamma_5 + \gamma_6 = (1, 3, 2, 0) - (7, 8, 1, 0) + (1, 8, 9, 0)$.

Let us now endow the complex $K'$ shown in Fig. 12(a) with a directed graph $G'$, shown in Fig. 12(b), in order to form the complex $T'$ of the Rössler attractor in Fig. 12. The set of operations required to combine the complex $K'$ and the digraph $G'$ are detailed in Charó et al.<fnref n="3" /> The joining 1-cell has two ingoing 2-cells, $\gamma_5$ and $\gamma_6$, and one outgoing 2-cell, $\gamma_1$. This information can be viewed in the joining subgraph $\{\gamma_5 \to \gamma_1, \gamma_6 \to \gamma_1\}$ of Fig. 13(b).

The existence of an orientation change within a connected set of joining 1-cells reveals the existence of a splitting 0-cell. In this example, no splitting 0-cells are found for $T'$. The stripexes are given by

$$ \gamma_1 \to \gamma_2 \to \gamma_4 \to \gamma_6 \to \gamma_1, \quad (\text{D2a}) $$

$$ \gamma_1 \to \gamma_2 \to \gamma_3 \to \gamma_5 \to \gamma_1. \quad (\text{D2b}) $$

The Rössler attractor is thus homologically equivalent to a cylinder, but a cylinder as such does not have a joining locus. The complex $T'$ tells us much more about the attractor's structure: its joining locus has a single component—as expected for an attractor bounded by a genus-1 torus—and there are two stripexes, $S_1$ given by (D2a) and $S_2$ given by (D2b). The stripex $S_1$ is a Möbius band while $S_2$ is a cylinder or normal band (without torsion).

All the computations herein can be handled algorithmically with the Wolfram Mathematica code freely available at git.cima.fcen.uba.ar/sciamarella/templex-properties.git.

## REFERENCES

<footnote n="1">¹J. Birman and R. F. Williams, “Knotted periodic orbits in dynamical systems. I: Lorenz’s equations,” *Topology* **22**, 47–82 (1983).</footnote>

<footnote n="2">²R. F. Williams, “Expanding attractors,” *Publ. Math. Inst. Hautes Études Sci.* **43**, 169–203 (1974).</footnote>

<footnote n="3">³G. D. Charó, C. Letellier, and D. Sciamarella, “Templex: A bridge between homologies and templates for chaotic attractors,” *Chaos* **32**, 083108 (2022).</footnote>

<footnote n="4">⁴M. D. Chekroun, E. Simonnet, and M. Ghil, “Stochastic climate dynamics: Random attractors and time-dependent invariant measures,” *Physica D* **240**, 1685–1700 (2011).</footnote>

<footnote n="5">⁵D. Sciamarella and G. B. Mindlin, “Topological structure of chaotic flows from human speech data,” *Phys. Rev. Lett.* **82**, 1450 (1999).</footnote>

<footnote n="6">⁶D. Sciamarella and G. B. Mindlin, “Unveiling the topological structure of chaotic flows from data,” *Phys. Rev. E* **64**, 036209 (2001).</footnote>

<footnote n="7">⁷G. D. Charó, M. D. Chekroun, D. Sciamarella, and M. Ghil, “Noise-driven topological changes in chaotic dynamics,” *Chaos* **31**, 103115 (2021).</footnote>

<footnote n="8">⁸J. Bang-Jensen and G. Z. Gutin, *Digraphs: Theory, Algorithms and Applications* (Springer Science & Business Media, 2008).</footnote>

<footnote n="9">⁹T. M. Lenton, H. Held, E. Kriegler, J. W. Hall, W. Lucht, S. Rahmstorf, and H. J. Schellnhuber, “Tipping elements in the Earth’s climate system,” *Proc. Natl. Acad. Sci. U. S. A.* **105**, 1786–1793 (2008).</footnote>

<footnote n="10">¹⁰H. Poincaré, “Analysis situs,” J. Éc. Polythec. Mat. **1**, 1–121 (1895).</footnote>

<footnote n="11">¹¹E. N. Lorenz, “Deterministic nonperiodic flow,” *J. Atmos. Sci.* **20**, 130–141 (1963).</footnote>

<footnote n="12">¹²L. C. Kinsey, *Topology of Surfaces* (Springer Science & Business Media, 2012).</footnote>

<footnote n="13">¹³D. Siersma, “Poincare and analysis situs, the beginning of algebraic topology,” *Nieuw Arch. Wisk. Serie 5* **13**, 196–200 (2012).</footnote>

<footnote n="14">¹⁴C. Letellier, P. Dutertre, and B. Maheu, “Unstable periodic orbits and templates of the Rössler system: Toward a systematic topological characterization,” *Chaos* **5**, 271–282 (1995).</footnote>

<footnote n="15">¹⁵M. Ghil, M. D. Chekroun, and E. Simonnet, “Climate dynamics and fluid mechanics: Natural variability and related uncertainties,” *Physica D* **237**, 2111–2126 (2008).</footnote>

<footnote n="16">¹⁶T. Bódai and T. Tél, “Annual variability in a conceptual climate model: Snapshot attractors, hysteresis in extreme events, and climate sensitivity,” *Chaos* **22**, 023110 (2012).</footnote>

<footnote n="17">¹⁷T. Tél, T. Bódai, G. Drótos, T. Haszpra, M. Herein, B. Kaszás, and M. Vincze, “The theory of parallel climate realizations,” *J. Stat. Phys.* **179**, 1496–1530 (2020).</footnote>

<footnote n="18">¹⁸H. Crauel and F. Flandoli, “Attractors for random dynamical systems,” *Probab. Theory Relat. Fields* **100**, 365–393 (1994).</footnote>

<footnote n="19">¹⁹L. Arnold, *Random Dynamical Systems* (Springer, 1988).</footnote>

<footnote n="20">²⁰F. J. Romeiras, C. Grebogi, and E. Ott, “Multifractal properties of snapshot attractors of random maps,” *Phys. Rev. A* **41**, 784 (1990).</footnote>

<footnote n="21">²¹T. Caravallo, G. Łukaszewicz, and J. Real, “Pullback attractors for asymptotically compact non-autonomous dynamical systems,” *Nonlinear Anal.* **64**, 484–498 (2006).</footnote>

<footnote n="22">²²M. Ghil, M. D. Chekroun, and E. Simonnet, “Climate dynamics and fluid mechanics: Natural variability and related uncertainties,” *Physica D* **237**, 2111–2126 (2008).</footnote>

<footnote n="23">²³D. Kondrashov, M. D. Chekroun, and M. Ghil, “Data-driven non-Markovian closure models,” *Physica D* **297**, 33–55 (2015).</footnote>

<footnote n="24">²⁴M. Santos Gutiérrez, V. Lucarini, M. D. Chekroun, and M. Ghil, “Reduced-order models for coupled dynamical systems: Data-driven methods and the Koopman operator,” *Chaos* **31**(5), 053116 (2021).</footnote>

<footnote n="25">²⁵T. Caraballo and X. Han, *Applied Nonautonomous and Random Dynamical Systems: Applied Dynamical Systems* (Springer, 2017).</footnote>

<footnote n="26">²⁶P. Kloeden and M. Yang, *An Introduction to Nonautonomous Dynamical Systems and Their Attractors* (World Scientific, 2020), Vol. 21.</footnote>

<footnote n="27">²⁷M. Ghil and S. Childress, *Topics in Geophysical Fluid Dynamics: Atmospheric Dynamics, Dynamo Theory, and Climate Dynamics* (Springer, 1987, reissued, 2012).</footnote>

<footnote n="28">²⁸M. Ghil and V. Lucarini, “The physics of climate variability and climate change,” *Rev. Mod. Phys.* **92**, 035002 (2020).</footnote>

<footnote n="29">²⁹M. D. Chekroun, M. Ghil, and J. D. Neelin, “Pullback attractor crisis in a delay differential ENSO model,” in *Advances in Nonlinear Geosciences*, edited by A. A. Tsonis (Springer Science & Business Media, 2018), pp. 1–33.</footnote>

<footnote n="30">³⁰T. Tél, T. Bódai, G. Drótos, T. Haszpra, M. Herein, B. Kaszás, and M. Vincze, “The theory of parallel climate realizations: A new framework of ensemble methods in a changing climate: An overview,” *J. Stat. Phys.* **179**, 1496–1530 (2020).</footnote>

<footnote n="31">³¹K. Hasselmann, “Stochastic climate models. I: Theory,” *Tellus* **28**, 473–485 (1976).</footnote>

<footnote n="32">³²F. Ledrappier and L.-S. Young, “Entropy formula for random transformations,” *Probab. Theory Relat. Fields* **80**, 217–240 (1988).</footnote>

<footnote n="33">³³L.-S. Young, “What are SRB measures, and which dynamical systems have them?,” *J. Stat. Phys.* **108**, 733–754 (2002).</footnote>

<footnote n="34">³⁴J.-P. Eckmann and D. Ruelle, “Ergodic theory of chaos and strange attractors,” *Rev. Mod. Phys.* **57**, 617–656 (1985).</footnote>

<footnote n="35">³⁵L.-S. Young, “Generalizations of SRB measures to nonautonomous, random, and infinite dimensional systems,” *J. Stat. Phys.* **166**, 494–515 (2017).</footnote>

<footnote n="36">³⁶L. Hörmander, “Hypoelliptic second order differential equations,” *Acta Math.* **119**, 147–171 (1967).</footnote>

<footnote n="37">³⁷O. E. Rössler, “An equation for continuous chaos,” *Phys. Lett. A* **57**, 397–398 (1976).</footnote>