Source author record

Gary Froyland

Gary Froyland appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2022arXiv

Perturbation formulae for quenched random dynamics with applications to open systems and extreme value theory

We consider quasi-compact linear operator cocycles $\mathcal{L}^{n}_ω:=\mathcal{L}_{σ^{n-1}ω}\circ\cdots\circ\mathcal{L}_{σω}\circ \mathcal{L}_ω$ driven by an invertible ergodic process $σ:Ω\toΩ$, and their small perturbations $\mathcal{L}_{ω,ε}^{n}$. We prove an abstract $ω$-wise first-order formula for the leading Lyapunov multipliers. We then consider the situation where $\mathcal{L}_ω^{n}$ is a transfer operator cocycle for a random map cocycle $T_ω^{n}:=T_{σ^{n-1}ω}\circ\cdots\circ T_{σω}\circ T_ω$ and the perturbed transfer operators $\mathcal{L}_{ω,ε}$ are defined by the introduction of small random holes $H_{ω,ε}$ in $[0,1]$, creating a random open dynamical system. We obtain a first-order perturbation formula in this setting, which reads $λ_{ω,ε}=λ_ω-θ_ωμ_ω(H_{ω,ε})+o(μ_ω(H_{ω,ε})),$ where $μ_ω$ is the unique equivariant random measure (and equilibrium state) for the original closed random dynamics. Our new machinery is then deployed to create a spectral approach for a quenched extreme value theory that considers random dynamics with general ergodic invertible driving, and random observations. An extreme value law is derived using the first-order terms $θ_ω$. Further, in the setting of random piecewise expanding interval maps, we establish the existence of random equilibrium states and conditionally invariant measures for random open systems via a random perturbative approach. Finally we prove quenched statistical limit theorems for random equilibrium states arising from contracting potentials. We illustrate the theory with a variety of explicit examples.

preprint2021arXiv

Deep Lagrangian connectivity in the global ocean inferred from Argo floats

We describe the application of a new technique from nonlinear dynamical systems to infer the Lagrangian connectivity of the deep global ocean. We approximate the dynamic Laplacian using Argo trajectories from January 2011 to January 2017 and extract the eight dominant coherent (or dynamically self-connected) regions at 1500m depth. Our approach overcomes issues such as sparsity of observed data, and floats continually leaving and entering the dataset; only 10\% of floats record for the full six years. The identified coherent regions maximally trap water within them over the six-year time frame, providing a distinct analysis of the deep global ocean, and relevant information for planning future float deployment. While our study is concerned with ocean circulation at a multi-year, global scale, the dynamic Laplacian approach may be applied at any temporal or spatial scale to identify coherent structures in ocean flow from positional time series information arising from observations or models.

preprint2021arXiv

Spectral analysis of climate dynamics with operator-theoretic approaches

The Earth's climate system is a classical example of a multiscale, multiphysics dynamical system with an extremely large number of active degrees of freedom, exhibiting variability on scales ranging from micrometers and seconds in cloud microphysics, to thousands of kilometers and centuries in ocean dynamics. Yet, despite this dynamical complexity, climate dynamics is known to exhibit coherent modes of variability. A primary example is the El Niño Southern Oscillation (ENSO), the dominant mode of interannual (3-5 yr) variability in the climate system. The objective and robust characterization of this and other important phenomena presents a long-standing challenge in Earth system science, the resolution of which would lead to improved scientific understanding and prediction of climate dynamics, as well as assessment of their impacts on human and natural systems. Here, we show that the spectral theory of dynamical systems, combined with techniques from data science, provides an effective means for extracting coherent modes of climate variability from high-dimensional model and observational data, requiring no frequency prefiltering, but recovering multiple timescales and their interactions. Lifecycle composites of ENSO are shown to improve upon results from conventional indices in terms of dynamical consistency and physical interpretability. In addition, the role of combination modes between ENSO and the annual cycle in ENSO diversity is elucidated.

preprint2021arXiv

Thermodynamic Formalism for Random Interval Maps with Holes

We develop a quenched thermodynamic formalism for open random dynamical systems generated by finitely branched, piecewise-monotone mappings of the interval. The openness refers to the presence of holes in the interval, which terminate trajectories once they enter; the holes may also be random. Our random driving is generated by an invertible, ergodic, measure-preserving transformation $σ$ on a probability space $(Ω,\mathscr{F},m)$. For each $ω\inΩ$ we associate a piecewise-monotone, surjective map $T_ω:I\to I$, and a hole $H_ω\subset I$; the map $T_ω$, the random potential $φ_ω$, and the hole $H_ω$ generate the corresponding open transfer operator $\mathcal{L}_ω$. For a contracting potential, under a condition on the open random dynamics in the spirit of Liverani--Maume-Deschamps, we prove there exists a unique random probability measure $ν_ω$ supported on the survivor set ${X}_{ω,\infty}$ satisfying $ν_{σ(ω)}(\mathcal{L}_ωf)=λ_ων_ω(f)$. We also prove the existence of a unique random family of functions $q_ω$ that satisfy $\mathcal{L}_ωq_ω=λ_ωq_{σ(ω)}$. These yield an ergodic random invariant measure $μ=νq$ supported on the global survivor set, while $q$ combined with the random closed conformal measure yields a unique random absolutely continuous conditional invariant measure (RACCIM) $η$ supported on $I$. We prove quasi-compactness of the transfer operator cocycle and exponential decay of correlations for $μ$. Finally, the escape rates of the random closed conformal measure and the RACCIM $η$ coincide, and are given in terms of the expected pressure, as is the Hausdorff dimension of the surviving set $X_{ω,\infty}$. We provide examples of our general theory, including random $β$-transformations and random Lasota-Yorke maps.

preprint2020arXiv

Computation and optimal perturbation of finite-time coherent sets for aperiodic flows without trajectory integration

Understanding the macroscopic behavior of dynamical systems is an important tool to unravel transport mechanisms in complex flows. A decomposition of the state space into coherent sets is a popular way to reveal this essential macroscopic evolution. To compute coherent sets from an aperiodic time-dependent dynamical system we consider the relevant transfer operators and their infinitesimal generators on an augmented space-time manifold. This space-time generator approach avoids trajectory integration, and creates a convenient linearization of the aperiodic evolution. This linearization can be further exploited to create a simple and effective spectral optimization methodology for diminishing or enhancing coherence. We obtain explicit solutions for these optimization problems using Lagrange multipliers and illustrate this technique by increasing and decreasing mixing of spatial regions through small velocity field perturbations.

preprint2016arXiv

Hölder continuity of Oseledets splittings for semi-invertible operator cocycles

For Hölder continuous cocycles over an invertible, Lipschitz base, we establish the Hölder continuity of Oseledets subspaces on compact sets of arbitrarily large measure. This extends a result of Araújo, Bufetov, and Filip by considering possibly noninvertible cocycles, which in addition may take values in the space of compact operators on a Hilbert space. As a by-product of our work, we also show that a noninvertible cocycle with nonvanishing Lyapunov exponents exhibits nonuniformly hyperbolic behaviour (in the sense of Pesin) on a set of full measure.

preprint2015arXiv

A rough-and-ready cluster-based approach for extracting finite-time coherent sets from sparse and incomplete trajectory data

We present a numerical method to identify regions of phase space that are approximately retained in a mobile compact neighbourhood over a finite time duration. Our approach is based on spatio-temporal clustering of trajectory data. The main advantages of the approach are the ability to produce useful results (i) when there are relatively few trajectories and (ii) when there are gaps in observation of the trajectories as can occur with real data. The method is easy to implement, works in any dimension, and is fast to run.

preprint2015arXiv

Dynamic isoperimetry and the geometry of Lagrangian coherent structures

The study of transport and mixing processes in dynamical systems is particularly important for the analysis of mathematical models of physical systems. We propose a novel, direct geometric method to identify subsets of phase space that remain strongly coherent over a finite time duration. This new method is based on a dynamic extension of classical (static) isoperimetric problems; the latter are concerned with identifying submanifolds with the smallest boundary size relative to their volume. The present work introduces \emph{dynamic} isoperimetric problems; the study of sets with small boundary size relative to volume \emph{as they are evolved by a general dynamical system}. We formulate and prove dynamic versions of the fundamental (static) isoperimetric (in)equalities; a dynamic Federer-Fleming theorem and a dynamic Cheeger inequality. We introduce a new dynamic Laplacian operator and describe a computational method to identify coherent sets based on eigenfunctions of the dynamic Laplacian. Our results include formal mathematical statements concerning geometric properties of finite-time coherent sets, whose boundaries can be regarded as Lagrangian coherent structures. The computational advantages of our new approach are a well-separated spectrum for the dynamic Laplacian, and flexibility in appropriate numerical approximation methods. Finally, we demonstrate that the dynamic Laplacian operator can be realised as a zero-diffusion limit of a newly advanced probabilistic transfer operator method (Froyland, 2013) for finding coherent sets, which is based on small diffusion. Thus, the present approach sits naturally alongside the probabilistic approach (Froyland, 2013), and adds a formal geometric interpretation.

preprint2015arXiv

On fast computation of finite-time coherent sets using radial basis functions

Finite-time coherent sets inhibit mixing over finite times. The most expensive part of the transfer operator approach to detecting coherent sets is the construction of the operator itself. We present a numerical method based on radial basis function collocation and apply it to a recent transfer operator construction that has been designed specifically for purely advective dynamics. The construction is based on a "dynamic" Laplacian operator and minimises the boundary size of the coherent sets relative to their volume. The main advantage of our new approach is a substantial reduction in the number of Lagrangian trajectories that need to be computed, leading to large speedups in the transfer operator analysis when this computation is costly.

preprint2014arXiv

A computational method to extract macroscopic variables and their dynamics in multiscale systems

This paper introduces coordinate-independent methods for analysing multiscale dynamical systems using numerical techniques based on the transfer operator and its adjoint. In particular, we present a method for testing whether an arbitrary dynamical system exhibits multiscale behaviour and for estimating the time-scale separation. For systems with such behaviour, we establish techniques for analysing the fast dynamics in isolation, extracting slow variables for the system, and accurately simulating these slow variables at a large time step. We illustrate our method with numerical examples and show how the reduced slow dynamics faithfully represents statistical features of the full dynamics which are not coordinate dependent.

preprint2014arXiv

A trajectory-free framework for analysing multiscale systems

We develop algorithms built around properties of the transfer operator and Koopman operator which 1) test for possible multiscale dynamics in a given dynamical system, 2) estimate the magnitude of the time-scale separation, and finally 3) distill the reduced slow dynamics on a suitably designed subspace. By avoiding trajectory integration, the developed techniques are highly computationally efficient. We corroborate our findings with numerical simulations of a test problem.

preprint2013arXiv

Stochastic stability of Lyapunov exponents and Oseledets splittings for semi-invertible matrix cocycles

We establish (i) stability of Lyapunov exponents and (ii) convergence in probability of Oseledets spaces for semi-invertible matrix cocycles, subjected to small random perturbations. The first part extends results of Ledrappier and Young to the semi-invertible setting. The second part relies on the study of evolution of subspaces in the Grassmannian.

preprint2013arXiv

Ulam's method for Lasota-Yorke maps with holes

Ulam's method is a rigorous numerical scheme for approximating invariant densities of dynamical systems. The phase space is partitioned into connected sets and an inter-set transition matrix is computed from the dynamics; an approximate invariant density is read off as the leading left eigenvector of this matrix. When a hole in phase space is introduced, one instead searches for \emph{conditional} invariant densities and their associated escape rates. For Lasota-Yorke maps with holes we prove that a simple adaptation of the standard Ulam scheme provides convergent sequences of escape rates (from the leading eigenvalue), conditional invariant densities (from the corresponding left eigenvector), and quasi-conformal measures (from the corresponding right eigenvector). We also immediately obtain a convergent sequence for the invariant measure supported on the survivor set. Our approach allows us to consider relatively large holes. We illustrate the approach with several families of examples, including a class of Lorenz maps.

preprint2012arXiv

An analytic framework for identifying finite-time coherent sets in time-dependent dynamical systems

The study of transport and mixing processes in dynamical systems is particularly important for the analysis of mathematical models of physical systems. Barriers to transport, which mitigate mixing, are currently the subject of intense study. In the autonomous setting, the use of transfer operators (Perron-Frobenius operators) to identify invariant and almost-invariant sets has been particularly successful. In the nonautonomous (time-dependent) setting, \emph{coherent sets}, a time-parameterised family of minimally dispersive sets, are a natural extension of almost-invariant sets. The present work introduces a new analytic transfer operator construction that enables the calculation of \emph{finite-time} coherent sets (sets are that minimally dispersive over a finite time interval). This new construction also elucidates the role of diffusion in the calculation and we show how properties such as the spectral gap and the regularity of singular vectors scale with noise amplitude. The construction can also be applied to general Markov processes on continuous state space.

preprint2012arXiv

Computing covariant vectors, Lyapunov vectors, Oseledets vectors, and dichotomy projectors: a comparative numerical study

Covariant vectors, Lyapunov vectors, or Oseledets vectors are increasingly being used for a variety of model analyses in areas such as partial differential equations, nonautonomous differentiable dynamical systems, and random dynamical systems. These vectors identify spatially varying directions of specific asymptotic growth rates and obey equivariance principles. In recent years new computational methods for approximating Oseledets vectors have been developed, motivated by increasing model complexity and greater demands for accuracy. In this numerical study we introduce two new approaches based on singular value decomposition and exponential dichotomies and comparatively review and improve two recent popular approaches of Ginelli et al. (2007) and Wolfe and Samelson (2007). We compare the performance of the four approaches via three case studies with very different dynamics in terms of symmetry, spectral separation, and dimension. We also investigate which methods perform well with limited data.

preprint2012arXiv

Metastability, Lyapunov exponents, escape rates, and topological entropy in random dynamical systems

We explore the concept of metastability in random dynamical systems, focussing on connections between random Perron-Frobenius operator cocycles and escape rates of random maps, and on topological entropy of random shifts of finite type. The Lyapunov spectrum of the random Perron-Frobenius cocycle and the random adjacency matrix cocycle is used to decompose the random system into two disjoint random systems with rigorous upper and lower bounds on (i) the escape rate in the setting of random maps, and (ii) topological entropy in the setting of random shifts of finite type, respectively.

preprint2012arXiv

Stability and approximation of random invariant densities for Lasota-Yorke map cocycles

We establish stability of random absolutely continuous invariant measures (acims) for cocycles of random Lasota-Yorke maps under a variety of perturbations. Our family of random maps need not be close to a fixed map; thus, our results can handle very general driving mechanisms. We consider (i) perturbations via convolutions, (ii) perturbations arising from finite-rank transfer operator approximation schemes and (iii) static perturbations, perturbing to a nearby cocycle of Lasota-Yorke maps. The former two results provide a rigorous framework for the numerical approximation of random acims using a Fourier-based approach and Ulam's method, respectively; we also demonstrate the efficacy of these schemes.

preprint2011arXiv

Estimating long term behavior of flows without trajectory integration: the infinitesimal generator approach

The long-term distributions of trajectories of a flow are described by invariant densities, i.e. fixed points of an associated transfer operator. In addition, global slowly mixing structures, such as almost-invariant sets, which partition phase space into regions that are almost dynamically disconnected, can also be identified by certain eigenfunctions of this operator. Indeed, these structures are often hard to obtain by brute-force trajectory-based analyses. In a wide variety of applications, transfer operators have proven to be very efficient tools for an analysis of the global behavior of a dynamical system. The computationally most expensive step in the construction of an approximate transfer operator is the numerical integration of many short term trajectories. In this paper, we propose to directly work with the infinitesimal generator instead of the operator, completely avoiding trajectory integration. We propose two different discretization schemes; a cell based discretization and a spectral collocation approach. Convergence can be shown in certain circumstances. We demonstrate numerically that our approach is much more efficient than the operator approach, sometimes by several orders of magnitude.

preprint2011arXiv

Spectral degeneracy and escape dynamics for intermittent maps with a hole

We study intermittent maps from the point of view of metastability. Small neighbourhoods of an intermittent fixed point and their complements form pairs of almost-invariant sets. Treating the small neighbourhood as a hole, we first show that the absolutely continuous conditional invariant measures (ACCIMs) converge to the ACIM as the length of the small neighbourhood shrinks to zero. We then quantify how the escape dynamics from these almost-invariant sets are connected with the second eigenfunctions of Perron-Frobenius (transfer) operators when a small perturbation is applied near the intermittent fixed point. In particular, we describe precisely the scaling of the second eigenvalue with the perturbation size, provide upper and lower bounds, and demonstrate $L^1$ convergence of the positive part of the second eigenfunction to the ACIM as the perturbation goes to zero. This perturbation and associated eigenvalue scalings and convergence results are all compatible with Ulam's method and provide a formal explanation for the numerical behaviour of Ulam's method in this nonuniformly hyperbolic setting. The main results of the paper are illustrated with numerical computations.

preprint2010arXiv

A semi-invertible Oseledets Theorem with applications to transfer operator cocycles

Oseledets' celebrated Multiplicative Ergodic Theorem (MET) is concerned with the exponential growth rates of vectors under the action of a linear cocycle on R^d. When the linear actions are invertible, the MET guarantees an almost-everywhere pointwise splitting of R^d into subspaces of distinct exponential growth rates (called Lyapunov exponents). When the linear actions are non-invertible, Oseledets' MET only yields the existence of a filtration of subspaces, the elements of which contain all vectors that grow no faster than exponential rates given by the Lyapunov exponents. The authors recently demonstrated that a splitting over R^d is guaranteed even without the invertibility assumption on the linear actions. Motivated by applications of the MET to cocycles of (non-invertible) transfer operators arising from random dynamical systems, we demonstrate the existence of an Oseledets splitting for cocycles of quasi-compact non-invertible linear operators on Banach spaces.

preprint2010arXiv

Optimally coherent sets in geophysical flows: A new approach to delimiting the stratospheric polar vortex

The "edge" of the Antarctic polar vortex is known to behave as a barrier to the meridional (poleward) transport of ozone during the austral winter. This chemical isolation of the polar vortex from the middle and low latitudes produces an ozone minimum in the vortex region, intensifying the ozone hole relative to that which would be produced by photochemical processes alone. Observational determination of the vortex edge remains an active field of research. In this letter, we obtain objective estimates of the structure of the polar vortex by introducing a new technique based on transfer operators that aims to find regions with minimal external transport. Applying this new technique to European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-40 three-dimensional velocity data we produce an improved three-dimensional estimate of the vortex location in the upper stratosphere where the vortex is most pronounced. This novel computational approach has wide potential application in detecting and analysing mixing structures in a variety of atmospheric, oceanographic, and general fluid dynamical settings.

preprint2010arXiv

Transport in time-dependent dynamical systems: Finite-time coherent sets

We study the transport properties of nonautonomous chaotic dynamical systems over a finite time duration. We are particularly interested in those regions that remain coherent and relatively non-dispersive over finite periods of time, despite the chaotic nature of the system. We develop a novel probabilistic methodology based upon transfer operators that automatically detects maximally coherent sets. The approach is very simple to implement, requiring only singular vector computations of a matrix of transitions induced by the dynamics. We illustrate our new methodology on an idealized stratospheric flow and in two and three dimensional analyses of European Centre for Medium Range Weather Forecasting (ECMWF) reanalysis data.