BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Georg A. Gottwald

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

ResearcherUnclaimed source record
math.DSnlin.CDphysics.ao-phphysics.data-annlin.AOmath.NAnlin.PSBiological Physicsmath.PROther Quantitative Biologyq-fin.PM

Catalog footprint

What is connected

24works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

24 published item(s)

preprint2022arXiv

On financial market correlation structures and diversification benefits across and within equity sectors

We study how to assess the potential benefit of diversifying an equity portfolio by investing within and across equity sectors. We analyse 20 years of US stock price data, which includes the global financial crisis (GFC) and the COVID-19 market crash, as well as periods of financial stability, to determine the `all weather' nature of equity portfolios. We establish that one may use the leading eigenvalue of the cross-correlation matrix of log returns as well as graph-theoretic diagnostics such as modularity to quantify the collective behaviour of the market or a subset of it. We confirm that financial crises are characterised by a high degree of collective behaviour of equities, whereas periods of financial stability exhibit less collective behaviour. We argue that during times of increased collective behaviour, risk reduction via sector-based portfolio diversification is ineffective. Using the degree of collectivity as a proxy for the benefit of diversification, we perform an extensive sampling of equity portfolios to confirm the old financial adage that 30-40 stocks provide sufficient diversification. Using hierarchical clustering, we discover a `best value' equity portfolio for diversification consisting of 36 equities sampled uniformly from 9 sectors. We further show that it is typically more beneficial to diversify across sectors rather than within. Our findings have implications for cost-conscious retail investors seeking broad diversification across equity markets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Model Reduction for the Kuramoto-Sakaguchi Model: The Importance of Non-entrained Rogue Oscillators

The Kuramoto-Sakaguchi model for coupled phase oscillators with phase-frustration is often studied in the thermodynamic limit of infinitely many oscillators. Here we extend a model reduction method based on collective coordinates to capture the collective dynamics of finite size Kuramoto-Sakaguchi models. We find that the inclusion of the effects of rogue oscillators is essential to obtain an accurate description, in contrast to the original Kuramoto model where we show that their effects can be ignored. We further introduce a more accurate ansatz function to describe the shape of synchronized oscillators. Our results from this extended collective coordinate approach reduce in the thermodynamic limit to the well-known mean-field consistency relations. For finite networks we show that our model reduction describes the collective behavior accurately, reproducing the order parameter, the mean frequency of the synchronized cluster, and the size of the cluster at given coupling strength, as well as the critical coupling strength for partial and for global synchronization.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stochastically perturbed bred vectors

The breeding method is a computationally cheap procedure to generate initial conditions for ensemble forecasting which project onto relevant synoptic growing modes. Ensembles of bred vectors, however, often lack diversity and align with the leading Lyapunov vector, which severely impacts their statistical reliability. In previous work we developed stochastically perturbed bred vectors (SPBVs) and random draw bred vectors (RDBVs) in the context of multi-scale systems. Here we explore when this method can be extended to systems without scale separation, and examine the performance of the stochastically modified bred vectors in the single scale Lorenz 96 model. In particular, we show that the performance of SPBVs crucially depends on the degree of localisation of the bred vectors. It is found that, contrary to the case of multi-scale systems, localisation is detrimental for applications of SPBVs in systems without scale-separation when initialised from assimilated data. In the case of weakly localised bred vectors, however, ensembles of SPBVs constitute a reliable ensemble with improved ensemble forecasting skills compared to classical bred vectors, while still preserving the low computational cost of the breeding method. RDBVs are shown to have superior forecast skill and form a reliable ensemble in weakly localised situations, but in situations when they are strongly localised they do not constitute a reliable ensemble and are over-dispersive.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Linear response for macroscopic observables in high-dimensional systems

The long-term average response of observables of chaotic systems to dynamical perturbations can often be predicted using linear response theory, but not all chaotic systems possess a linear response. Macroscopic observables of complex dissipative chaotic systems, however, are widely assumed to have a linear response even if the microscopic variables do not, but the mechanism for this is not well-understood. We present a comprehensive picture for the linear response of macroscopic observables in high-dimensional coupled deterministic dynamical systems, where the coupling is via a mean field and the microscopic subsystems may or may not obey linear response theory. We derive stochastic reductions of the dynamics of these observables from statistics of the microscopic system, and provide conditions for linear response theory to hold in finite dimensional systems and in the thermodynamic limit. In particular, we show that for large systems of finite size, linear response is induced via self-generated noise. We present examples in the thermodynamic limit where the macroscopic observable satisfies LRT, although the microscopic subsystems individually violate LRT, as well a converse example where the macroscopic observable does not satisfy LRT despite all microscopic subsystems satisfying LRT when uncoupled. This latter, maybe surprising, example is associated with emergent non-trivial dynamics of the macroscopic observable. We provide numerical evidence for our results on linear response as well as some analytical intuition.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Broadband nature of power spectra for intermittent Maps with summable and nonsummable decay of correlations

We present results on the broadband nature of the power spectrum $S(ω)$, $ω\in(0,2π)$, for a large class of nonuniformly expanding maps with summable and nonsummable decay of correlations. In particular, we consider a class of intermittent maps $f:[0,1]\to[0,1]$ with $f(x)\approx x^{1+γ}$ for $x\approx 0$, where $γ\in(0,1)$. Such maps have summable decay of correlations when $γ\in(0,\frac12)$, and $S(ω)$ extends to a continuous function on $[0,2π]$ by the classical Wiener-Khintchine Theorem. We show that $S(ω)$ is typically bounded away from zero for Hölder observables. Moreover, in the nonsummable case $γ\in[\frac12,1)$, we show that $S(ω)$ is defined almost everywhere with a continuous extension $\tilde S(ω)$ defined on $(0,2π)$, and $\tilde S(ω)$ is typically nonvanishing.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the detection of superdiffusive behaviour in time series

We present a new method for detecting superdiffusive behaviour and for determining rates of superdiffusion in time series data. Our method applies equally to stochastic and deterministic time series data (with no prior knowledge required of the nature of the data) and relies on one realisation (ie one sample path) of the process. Linear drift effects are automatically removed without any preprocessing. We show numerical results for time series constructed from i.i.d. $α$-stable random variables and from deterministic weakly chaotic maps. We compare our method with the standard method of estimating the growth rate of the mean-square displacement as well as the $p$-variation method, maximum likelihood, quantile matching and linear regression of the empirical characteristic function.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stochastic Climate Theory

In this chapter we review stochastic modelling methods in climate science. First we provide a conceptual framework for stochastic modelling of deterministic dynamical systems based on the Mori-Zwanzig formalism. The Mori-Zwanzig equations contain a Markov term, a memory term and a term suggestive of stochastic noise. Within this framework we express standard model reduction methods such as averaging and homogenization which eliminate the memory term. We further discuss ways to deal with the memory term and how the type of noise depends on the underlying deterministic chaotic system. Secondly, we review current approaches in stochastic data-driven models. We discuss how the drift and diffusion coefficients of models in the form of stochastic differential equations can be estimated from observational data. We pay attention to situations where the data stems from multi scale systems, a relevant topic in the context of data from the climate system. Furthermore, we discuss the use of discrete stochastic processes (Markov chains) for e.g. stochastic subgrid-scale modeling and other topics in climate science.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A data-driven method for the stochastic parametrisation of subgrid-scale tropical convective area fraction

Observations of tropical convection from precipitation radar and the concurring large-scale atmospheric state at two locations (Darwin and Kwajalein) are used to establish effective stochastic models to parameterise subgrid-scale tropical convective activity. Two approaches are presented which rely on the assumption that tropical convection induces a stationary equilibrium distribution. In the first approach we parameterise convection variables such as convective area fraction as an instantaneous random realisation conditioned on the large-scale vertical velocities according to a probability density function estimated from the observations. In the second approach convection variables are generated in a Markov process conditioned on the large-scale vertical velocity, allowing for non-trivial temporal correlations. Despite the different prevalent atmospheric and oceanic regimes at the two locations, with Kwajalein being exposed to a purely oceanic weather regime and Darwin exhibiting land-sea interaction, we establish that the empirical measure for the convective variables conditioned on large-scale mid-level vertical velocities for the two locations are close. This allows us to train the stochastic models at one location and then generate time series of convective activity at the other location. The proposed stochastic subgrid-scale models adequately reproduce the statistics of the observed convective variables and we discuss how they may be used in future scale-independent mass-flux convection parameterisations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Central limit theorems and suppression of anomalous diffusion for systems with symmetry

We give general conditions for the central limit theorem and weak convergence to Brownian motion (the weak invariance principle / functional central limit theorem) to hold for observables of compact group extensions of nonuniformly expanding maps. In particular, our results include situations where the central limit theorem would fail, and anomalous behaviour would prevail, if the compact group were not present. This has important consequences for systems with noncompact Euclidean symmetry and provides the rigorous proof for a conjecture made in our paper: A Huygens principle for diffusion and anomalous diffusion in spatially extended systems. Proc. Natl. Acad. Sci. USA 110 (2013) 8411-8416.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Homogenization for Deterministic Maps and Multiplicative Noise

A recent paper of Melbourne & Stuart, A note on diffusion limits of chaotic skew product flows, Nonlinearity 24 (2011) 1361-1367, gives a rigorous proof of convergence of a fast-slow deterministic system to a stochastic differential equation with additive noise. In contrast to other approaches, the assumptions on the fast flow are very mild. In this paper, we extend this result from continuous time to discrete time. Moreover we show how to deal with one-dimensional multiplicative noise. This raises the issue of how to interpret certain stochastic integrals; it is proved that the integrals are of Stratonovich type for continuous time and neither Stratonovich nor Ito for discrete time. We also provide a rigorous derivation of superdiffusive limit where the stochastic differential equation is driven by a stable Levy process. In the case of one-dimensional multiplicative noise, the stochastic integrals are of Marcus type both in the discrete and continuous time contexts.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Model reduction for networks of coupled oscillators

We present a collective coordinate approach to describe coupled phase oscillators. We apply the method to study synchronisation in a Kuramoto model. In our approach an N-dimensional Kuramoto model is reduced to an n-dimensional ordinary differential equation with n<<N, constituting an immense reduction in complexity. The onset of both local and global synchronisation is reproduced to good numerical accuracy, and we are able to describe both soft and hard transitions. By introducing 2 collective coordinates the approach is able to describe the interaction of two partially synchronised clusters in the case of bimodally distributed native frequencies. Furthermore, our approach allows us to accurately describe finite size scalings of the critical coupling strength. We corroborate our analytical results by comparing with numerical simulations of the Kuramoto model with all-to-all coupling networks for several distributions of the native frequencies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On convergence of higher order schemes for the projective integration method for stiff ordinary differential equations

We present a convergence proof for higher order implementations of the projective integration method (PI) for a class of deterministic multi-scale systems in which fast variables quickly settle on a slow manifold. The error is shown to contain contributions associated with the length of the microsolver, the numerical accuracy of the macrosolver and the distance from the slow manifold caused by the combined effect of micro- and macrosolvers, respectively. We also provide stability conditions for the PI methods under which the fast variables will not diverge from the slow manifold. We corroborate our results by numerical simulations.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A test for a conjecture on the nature of attractors for smooth dynamical systems

Dynamics arising persistently in smooth dynamical systems ranges from regular dynamics (periodic, quasiperiodic) to strongly chaotic dynamics (Anosov, uniformly hyperbolic, nonuniformly hyperbolic modelled by Young towers). The latter include many classical examples such as Lorenz and Hénon-like attractors and enjoy strong statistical properties. It is natural to conjecture (or at least hope) that most dynamical systems fall into these two extreme situations. We describe a numerical test for such a conjecture/hope and apply this to the logistic map where the conjecture holds by a theorem of Lyubich, and to the Lorenz-96 system in 40 dimensions where there is no rigorous theory. The numerical outcome is almost identical for both (except for the amount of data required) and provides evidence for the validity of the conjecture.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Huygens principle for diffusion and anomalous diffusion in spatially extended systems

We present a universal view on diffusive behaviour in chaotic spatially extended systems for anisotropic and isotropic media. For anisotropic systems, strong chaos leads to diffusive behaviour (Brownian motion with drift) and weak chaos leads to superdiffusive behaviour (Lévy processes with drift). For isotropic systems, the drift term vanishes and strong chaos again leads to Brownian motion. We establish the existence of a nonlinear Huygens principle for weakly chaotic systems in isotropic media whereby the dynamics behaves diffusively in even space dimension and exhibits superdiffusive behaviour in odd space dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On convergence of the projective integration method for stiff ordinary differential equations

We present a convergence proof of the projective integration method for a class of deterministic multi-dimensional multi-scale systems which are amenable to centre manifold theory. The error is shown to contain contributions associated with the numerical accuracy of the microsolver, the numerical accuracy of the macrosolver and the distance from the centre manifold caused by the combined effect of micro- and macrosolvers, respectively. We corroborate our results by numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On finite-size Lyapunov exponents in multiscale systems

We study the effect of regime switches on finite size Lyapunov exponents (FSLEs) in determining the error growth rates and predictability of multiscale systems. We consider a dynamical system involving slow and fast regimes and switches between them. The surprising result is that due to the presence of regimes the error growth rate can be a non-monotonic function of initial error amplitude. In particular, troughs in the large scales of FSLE spectra is shown to be a signature of slow regimes, whereas fast regimes are shown to cause large peaks in the spectra where error growth rates far exceed those estimated from the maximal Lyapunov exponent. We present analytical results explaining these signatures and corroborate them with numerical simulations. We show further that these peaks disappear in stochastic parametrizations of the fast chaotic processes, and the associated FSLE spectra reveal that large scale predictability properties of the full deterministic model are well approximated whereas small scale features are not properly resolved.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Controlling overestimation of error covariance in ensemble Kalman filters with sparse observations: A variance limiting Kalman filter

We consider the problem of an ensemble Kalman filter when only partial observations are available. In particular we consider the situation where the observational space consists of variables which are directly observable with known observational error, and of variables of which only their climatic variance and mean are given. To limit the variance of the latter poorly resolved variables we derive a variance limiting Kalman filter (VLKF) in a variational setting. We analyze the variance limiting Kalman filter for a simple linear toy model and determine its range of optimal performance. We explore the variance limiting Kalman filter in an ensemble transform setting for the Lorenz-96 system, and show that incorporating the information of the variance of some un-observable variables can improve the skill and also increase the stability of the data assimilation procedure.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Data assimilation in slow-fast systems using homogenized climate models

A deterministic multiscale toy model is studied in which a chaotic fast subsystem triggers rare transitions between slow regimes, akin to weather or climate regimes. Using homogenization techniques, a reduced stochastic parametrization model is derived for the slow dynamics. The reliability of this reduced climate model in reproducing the statistics of the slow dynamics of the full deterministic model for finite values of the time scale separation is numerically established. The statistics however is sensitive to uncertainties in the parameters of the stochastic model. It is investigated whether the stochastic climate model can be beneficial as a forecast model in an ensemble data assimilation setting, in particular in the realistic setting when observations are only available for the slow variables. The main result is that reduced stochastic models can indeed improve the analysis skill, when used as forecast models instead of the perfect full deterministic model. The stochastic climate model is far superior at detecting transitions between regimes. The observation intervals for which skill improvement can be obtained are related to the characteristic time scales involved. The reason why stochastic climate models are capable of producing superior skill in an ensemble setting is due to the finite ensemble size; ensembles obtained from the perfect deterministic forecast model lacks sufficient spread even for moderate ensemble sizes. Stochastic climate models provide a natural way to provide sufficient ensemble spread to detect transitions between regimes. This is corroborated with numerical simulations. The conclusion is that stochastic parametrizations are attractive for data assimilation despite their sensitivity to uncertainties in the parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The large core limit of spiral waves in excitable media: A numerical approach

We modify the freezing method introduced by Beyn & Thuemmler, 2004, for analyzing rigidly rotating spiral waves in excitable media. The proposed method is designed to stably determine the rotation frequency and the core radius of rotating spirals, as well as the approximate shape of spiral waves in unbounded domains. In particular, we introduce spiral wave boundary conditions based on geometric approximations of spiral wave solutions by Archimedean spirals and by involutes of circles. We further propose a simple implementation of boundary conditions for the case when the inhibitor is non-diffusive, a case which had previously caused spurious oscillations. We then utilize the method to numerically analyze the large core limit. The proposed method allows us to investigate the case close to criticality where spiral waves acquire infinite core radius and zero rotation frequency, before they begin to develop into retracting fingers. We confirm the linear scaling regime of a drift bifurcation for the rotation frequency and the core radius of spiral wave solutions close to criticality. This regime is unattainable with conventional numerical methods.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the Implementation of the 0-1 Test for Chaos

In this paper we address practical aspects of the implementation of the 0-1 test for chaos in deterministic systems. In addition, we present a new formulation of the test which significantly increases its sensitivity. The test can be viewed as a method to distill a binary quantity from the power spectrum. The implementation is guided by recent results from the theoretical justification of the test as well as by exploring better statistical methods to determine the binary quantities. We give several examples to illustrate the improvement.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Testing for Chaos in Deterministic Systems with Noise

Recently, we introduced a new test for distinguishing regular from chaotic dynamics in deterministic dynamical systems and argued that the test had certain advantages over the traditional test for chaos using the maximal Lyapunov exponent. In this paper, we investigate the capability of the test to cope with moderate amounts of noisy data. Comparisons are made between an improved version of our test and both the ``tangent space'' and ``direct method'' for computing the maximal Lyapunov exponent. The evidence of numerical experiments, ranging from the logistic map to an eight-dimensional Lorenz system of differential equations (the Lorenz 96 system), suggests that our method is superior to tangent space methods and that it compares very favourably with direct methods.

0 citations0 reviewsNo code linkedOpen access