Showing 1–24 of 24 results for author: Gottwald, A

Data assimilation for networks of coupled oscillators: Inferring unknown model parameters from partial observations

Authors: Lauren D. Smith, Georg A. Gottwald

Abstract: We develop a data-driven method combining a forecast model with unknown parameters and partial noisy observations for networks of coupled oscillators. Employing network-specific localization of the forecast covariance, an Ensemble Kalman Filter with state space augmentation is shown to yield highly accurate estimates of both the oscillator phases and unknown model parameters in the case where only… ▽ More We develop a data-driven method combining a forecast model with unknown parameters and partial noisy observations for networks of coupled oscillators. Employing network-specific localization of the forecast covariance, an Ensemble Kalman Filter with state space augmentation is shown to yield highly accurate estimates of both the oscillator phases and unknown model parameters in the case where only a subset of oscillator phases are observed. In contrast, standard data assimilation methods yield poor results. We demonstrate the effectiveness of our approach for Kuramoto oscillators and for networks of theta neurons. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Heterostructure films of SiO$_2$ and HfO$_2$ for high power laser optics prepared by plasma-enhanced atomic layer deposition

Authors: Shawon Alam, Pallabi Paul, Vivek Beladiya, Olaf Stenzel, Marcus Trost, Steffen Wilbrandt, Sven Schroeder, Gabor Matthaeus, Stefan Nolte, Sebastian Riese, Felix Otto, Torsten Fritz, Alexander Gottwald, Adriana Szeghalmi

Abstract: Absorption losses and laser-induced damage threshold (LIDT) are considered as the major constraint for the development of optical coatings for high-power laser optics. Such coatings require paramount properties like low losses due to optical absorption, high mechanical stability, and enhanced damage resistance to withstand high-intensity laser pulses. In this work, heterostructure films were devel… ▽ More Absorption losses and laser-induced damage threshold (LIDT) are considered as the major constraint for the development of optical coatings for high-power laser optics. Such coatings require paramount properties like low losses due to optical absorption, high mechanical stability, and enhanced damage resistance to withstand high-intensity laser pulses. In this work, heterostructure films were developed by the intermixing of SiO$_2$ and HfO$_2$ using plasma-enhanced atomic layer deposition (PEALD) technique. Thin film characterization techniques such as spectroscopic ellipsometry, spectrophotometry, substrate curvature measurements, x-ray reflectivity, and Fourier transform infrared spectroscopy were employed for extracting optical constants, spectral inter-pretation, residual stress, layer formation, and functional groups present in the heterostructures, respectively. These heterostructures demonstrate tunable refractive index, bandgap, and improved optical losses and LIDT properties. The films were incorporated into antireflection coatings (multilayer stacks and graded index coatings) and the LIDT was determined at 355 nm wavelength by the R-on-1 method. Optical absorptions at the reported wavelengths were characterized using photothermal common-path interferometry and laser-induced deflection techniques. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Benchmarking theoretical electronic structure methods with photoemission orbital tomography

Authors: Anja Haags, Xiaosheng Yang, Larissa Egger, Dominik Brandstetter, Hans Kirschner, Alexander Gottwald, Mathias Richter, Georg Koller, Michael G. Ramsey, François C. Bocquet, Serguei Soubatch, F. Stefan Tautz, Peter Puschnig

Abstract: In the past decade, photoemission orbital tomography (POT) has evolved into a powerful tool to investigate the electronic structure of organic molecules adsorbed on (metallic) surfaces. By measuring the angular distribution of photoelectrons as a function of binding energy and making use of the momentum-space signature of molecular orbitals, POT leads to an orbital-resolved picture of the electron… ▽ More In the past decade, photoemission orbital tomography (POT) has evolved into a powerful tool to investigate the electronic structure of organic molecules adsorbed on (metallic) surfaces. By measuring the angular distribution of photoelectrons as a function of binding energy and making use of the momentum-space signature of molecular orbitals, POT leads to an orbital-resolved picture of the electronic density of states at the organic/metal interface. In this combined experimental and theoretical work, we apply POT to the prototypical organic $π$-conjugated molecule bisanthene (C$_{28}$H$_{14}$) which forms a highly oriented monolayer on a Cu(110) surface. Experimentally, we identify an unprecedented number of 13 $π$ and 12 $σ$ orbitals of bisanthene and measure their respective binding energies and spectral lineshapes at the bisanthene/Cu(110) interface. Theoretically, we perform density functional calculations for this interface employing four widely used exchange-correlation functionals from the families of the generalized gradient approximations as well as global and range-separated hybrid functionals. By analyzing the electronic structure in terms of orbital-projected density of states, we arrive at a detailed orbital-by-orbital assessment of theory vs. experiment. This allows us to benchmark the performance of the investigated functionals with regards to their capability of accounting for the orbital energy alignment at organic/metal interfaces. △ Less

Submitted 27 September, 2022; v1 submitted 23 September, 2022; originally announced September 2022.

Lévy flights as an emergent phenomenon in a spatially extended system

Authors: Chunxi Jiao, Georg A. Gottwald

Abstract: Anomalous diffusion and Lévy flights, which are characterized by the occurrence of random discrete jumps of all scales, have been observed in a plethora of natural and engineered systems, ranging from the motion of molecules to climate signals. Mathematicians have recently unveiled mechanisms to generate anomalous diffusion, both stochastically and deterministically. However, there exists to the b… ▽ More Anomalous diffusion and Lévy flights, which are characterized by the occurrence of random discrete jumps of all scales, have been observed in a plethora of natural and engineered systems, ranging from the motion of molecules to climate signals. Mathematicians have recently unveiled mechanisms to generate anomalous diffusion, both stochastically and deterministically. However, there exists to the best of our knowledge no explicit example of a spatially extended system which exhibits anomalous diffusion without being explicitly driven by Lévy noise. We show here that the Landau-Lifshitz-Gilbert equation, a stochastic partial differential equation (SPDE), despite only driven by Gaussian white noise, exhibits superdiffusive behaviour. The anomalous diffusion is an entirely emergent behaviour and manifests itself in jumps in the location of its travelling front solution. Using a collective coordinate approach we reduce the SPDE to a set of stochastic differential equations (SDEs) driven by Gaussian white noise. This allows us to identify the mechanism giving rise to the anomalous diffusion as random widening events of the front interface. △ Less

Submitted 1 September, 2023; v1 submitted 3 September, 2022; originally announced September 2022.

Momentum-space imaging of σ-orbitals for chemical analysis

Authors: Anja Haags, Xiaosheng Yang, Larissa Egger, Dominik Brandstetter, Hans Kirschner, François C. Bocquet, Georg Koller, Alexander Gottwald, Mathias Richter, J. Michael Gottfried, Michael G. Ramsey, Peter Puschnig, Serguei Soubatch, F. Stefan Tautz

Abstract: Tracing the modifications of molecules in surface chemical reactions benefits from the possibility to image their orbitals. While delocalized frontier orbitals with π-character are imaged routinely with photoemission orbital tomography, they are not always sensitive to local chemical modifications, particularly the making and breaking of bonds at the molecular periphery. For such bonds, σ-orbitals… ▽ More Tracing the modifications of molecules in surface chemical reactions benefits from the possibility to image their orbitals. While delocalized frontier orbitals with π-character are imaged routinely with photoemission orbital tomography, they are not always sensitive to local chemical modifications, particularly the making and breaking of bonds at the molecular periphery. For such bonds, σ-orbitals would be far more revealing. Here, we show that these orbitals can indeed be imaged in a remarkably broad energy range, and that the plane wave approximation, an important ingredient of photoemission orbital tomography, is also well fulfilled for these orbitals. This makes photoemission orbital tomography a unique tool for the detailed analysis of surface chemical reactions. We demonstrate this by identifying the reaction product of a dehalogenation and cyclodehydrogenation reaction. △ Less

Submitted 30 October, 2021; originally announced November 2021.

Journal ref: Sci. Adv. 8, eabn0819 (2022)

Combining machine learning and data assimilation to forecast dynamical systems from noisy partial observations

Authors: Georg A. Gottwald, Sebastian Reich

Abstract: We present a supervised learning method to learn the propagator map of a dynamical system from partial and noisy observations. In our computationally cheap and easy-to-implement framework a neural network consisting of random feature maps is trained sequentially by incoming observations within a data assimilation procedure. By employing Takens' embedding theorem, the network is trained on delay co… ▽ More We present a supervised learning method to learn the propagator map of a dynamical system from partial and noisy observations. In our computationally cheap and easy-to-implement framework a neural network consisting of random feature maps is trained sequentially by incoming observations within a data assimilation procedure. By employing Takens' embedding theorem, the network is trained on delay coordinates. We show that the combination of random feature maps and data assimilation, called RAFDA, outperforms standard random feature maps for which the dynamics is learned using batch data. △ Less

Submitted 2 September, 2021; v1 submitted 7 August, 2021; originally announced August 2021.

Supervised learning from noisy observations: Combining machine-learning techniques with data assimilation

Authors: Georg A. Gottwald, Sebastian Reich

Abstract: Data-driven prediction and physics-agnostic machine-learning methods have attracted increased interest in recent years achieving forecast horizons going well beyond those to be expected for chaotic dynamical systems. In a separate strand of research data-assimilation has been successfully used to optimally combine forecast models and their inherent uncertainty with incoming noisy observations. The… ▽ More Data-driven prediction and physics-agnostic machine-learning methods have attracted increased interest in recent years achieving forecast horizons going well beyond those to be expected for chaotic dynamical systems. In a separate strand of research data-assimilation has been successfully used to optimally combine forecast models and their inherent uncertainty with incoming noisy observations. The key idea in our work here is to achieve increased forecast capabilities by judiciously combining machine-learning algorithms and data assimilation. We combine the physics-agnostic data-driven approach of random feature maps as a forecast model within an ensemble Kalman filter data assimilation procedure. The machine-learning model is learned sequentially by incorporating incoming noisy observations. We show that the obtained forecast model has remarkably good forecast skill while being computationally cheap once trained. Going beyond the task of forecasting, we show that our method can be used to generate reliable ensembles for probabilistic forecasting as well as to learn effective model closure in multi-scale systems. △ Less

Submitted 8 March, 2021; v1 submitted 14 July, 2020; originally announced July 2020.

A model for Dansgaard-Oeschger events and millennial-scale abrupt climate change without external forcing

Authors: Georg A. Gottwald

Abstract: We propose a conceptual model which generates abrupt climate changes akin to Dansgaard-Oeschger events. In the model these abrupt climate changes are not triggered by external perturbations but rather emerge in a dynamic self-consistent model through complex interactions of the ocean, the atmosphere and an intermittent process. The abrupt climate changes are caused in our model by intermittencies… ▽ More We propose a conceptual model which generates abrupt climate changes akin to Dansgaard-Oeschger events. In the model these abrupt climate changes are not triggered by external perturbations but rather emerge in a dynamic self-consistent model through complex interactions of the ocean, the atmosphere and an intermittent process. The abrupt climate changes are caused in our model by intermittencies in the sea-ice cover. The ocean is represented by a Stommel two-box model, the atmosphere by a Lorenz-84 model and the sea-ice cover by a deterministic approximation of correlated additive and multiplicative noise (CAM) process. The key dynamical ingredients of the model are given by stochastic limits of deterministic multi-scale systems and recent results in deterministic homogenisation theory. The deterministic model reproduces statistical features of actual ice-core data such as non-Gaussian $α$-stable behaviour. The proposed mechanism for abrupt millenial-scale climate change only relies on the existence of a quantity, which exhibits intermittent dynamics on an intermediate time scale. We consider as a particular mechanism intermittent sea-ice cover where the intermittency is generated by emergent atmospheric noise. However, other mechanisms such as freshwater influxes may also be formulated within the proposed framework. △ Less

Submitted 23 September, 2020; v1 submitted 30 May, 2020; originally announced June 2020.

Stochastically perturbed bred vectors

Authors: Brent Giggins, Georg A. Gottwald

Abstract: 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… ▽ More 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. △ Less

Submitted 9 August, 2020; v1 submitted 10 March, 2020; originally announced March 2020.

Detecting regime transitions in time series using dynamic mode decomposition

Authors: Georg A. Gottwald, Federica Gugole

Abstract: We employ the framework of the Koopman operator and dynamic mode decomposition to devise a computationally cheap and easily implementable method to detect transient dynamics and regime changes in time series. We argue that typically transient dynamics experiences the full state space dimension with subsequent fast relaxation towards the attractor. In equilibrium, on the other hand, the dynamics ev… ▽ More We employ the framework of the Koopman operator and dynamic mode decomposition to devise a computationally cheap and easily implementable method to detect transient dynamics and regime changes in time series. We argue that typically transient dynamics experiences the full state space dimension with subsequent fast relaxation towards the attractor. In equilibrium, on the other hand, the dynamics evolves on a slower time scale on a lower dimensional attractor. The reconstruction error of a dynamic mode decomposition is used to monitor the inability of the time series to resolve the fast relaxation towards the attractor as well as the effective dimension of the dynamics. We illustrate our method by detecting transient dynamics in the Kuramoto-Sivashinsky equation. We further apply our method to atmospheric reanalysis data; our diagnostics detects the transition from a predominantly negative North Atlantic Oscillation (NAO) to a predominantly positive NAO around 1970, as well as the recently found regime change in the Southern Hemisphere atmospheric circulation around 1970. △ Less

Submitted 3 September, 2019; v1 submitted 19 April, 2019; originally announced April 2019.

On the impossibility of solitary Rossby waves in meridionally unbounded domains

Authors: Georg A. Gottwald, Dmitry E. Pelinovsky

Abstract: Evolution of weakly nonlinear and slowly varying Rossby waves in planetary atmospheres and oceans is considered within the quasi-geostrophic equation on unbounded domains. When the mean flow profile has a jump in the ambient potential vorticity, localized eigenmodes are trapped by the mean flow with a non-resonant speed of propagation. We address amplitude equations for these modes. Whereas the li… ▽ More Evolution of weakly nonlinear and slowly varying Rossby waves in planetary atmospheres and oceans is considered within the quasi-geostrophic equation on unbounded domains. When the mean flow profile has a jump in the ambient potential vorticity, localized eigenmodes are trapped by the mean flow with a non-resonant speed of propagation. We address amplitude equations for these modes. Whereas the linear problem is suggestive of a two-dimensional Zakharov-Kuznetsov equation, we found that the dynamics of Rossby waves is effectively linear and moreover confined to zonal waveguides of the mean flow. This eliminates even the ubiquitous Korteweg-de Vries equations as underlying models for spatially localized coherent structures in these geophysical flows. △ Less

Submitted 7 August, 2018; originally announced August 2018.

Transverse resonance island buckets for synchrotron-radiation based electron time-of-flight spectroscopy

Authors: T. Arion, W. Eberhardt, J. Feikes, A. Gottwald, P. Goslawski, A. Hoehl, H. Kaser, M. Kolbe, J. Li, C. Lupulescu, M. Richter, M. Ries, F. Roth, M. Ruprecht, T. Tydecks, G. Wüstefeld

Abstract: At the Metrology Light Source (MLS), the compact electron storage ring of the Physikalisch-Technische Bundesanstalt (PTB) with a circumference of 48\,m, a specific operation mode with two stable closed orbits for stored electrons was realized by transverse resonance island buckets. One of these orbits is closing only after three turns. In combination with single-bunch operation, the new mode was a… ▽ More At the Metrology Light Source (MLS), the compact electron storage ring of the Physikalisch-Technische Bundesanstalt (PTB) with a circumference of 48\,m, a specific operation mode with two stable closed orbits for stored electrons was realized by transverse resonance island buckets. One of these orbits is closing only after three turns. In combination with single-bunch operation, the new mode was applied for electron time-of-flight spectroscopy with an interval of the synchrotron radiation pulses which is three times the revolution period at the MLS of 160\,ns. The achievement is of significant importance for PTB's future programs of angular-resolved electron spectroscopy with synchrotron radiation and similar projects at other compact electron storage rings. Moreover, the applied scheme for orbit and source spot selection via optical imaging at the insertion device beamline of the MLS and may be relevant for the BESSY VSR project of the Helmholtz-Zentrum Berlin. △ Less

Submitted 26 June, 2018; originally announced June 2018.

Report number: Rev. Sci. Instrum. 89, 103114 (2018)

Stochastically perturbed bred vectors in multi-scale systems

Authors: Brent Giggins, Georg A. Gottwald

Abstract: The breeding method is a computationally cheap way to generate flow-adapted ensembles to be used in probabilistic forecasts. Its main disadvantage is that the ensemble may lack diversity and collapse to a low-dimensional subspace. To still benefit from the breeding method's simplicity and its low computational cost, approaches are needed to increase the diversity of these bred vector (BV) ensemble… ▽ More The breeding method is a computationally cheap way to generate flow-adapted ensembles to be used in probabilistic forecasts. Its main disadvantage is that the ensemble may lack diversity and collapse to a low-dimensional subspace. To still benefit from the breeding method's simplicity and its low computational cost, approaches are needed to increase the diversity of these bred vector (BV) ensembles. We present here such a method tailored for multi-scale systems. We describe how to judiciously introduce stochastic perturbations to the standard bred vectors leading to stochastically perturbed bred vectors. The increased diversity leads to a better forecast skill as measured by the RMS error, as well as to more reliable ensembles quantified by the error-spread relationship, the continuous ranked probability score and reliability diagrams. Our approach is dynamically informed and in effect generates random draws from the fast equilibrium measure conditioned on the slow variables. We illustrate the advantage of stochastically perturbed bred vectors over standard BVs in numerical simulations of a multi-scale Lorenz 96 model. △ Less

Submitted 9 December, 2018; v1 submitted 24 April, 2018; originally announced April 2018.

Comments: accepted for publication in Q.J.R. Meteorolog. Soc

Optimal balance via adiabatic invariance of approximate slow manifolds

Authors: Georg A. Gottwald, Haidar Mohamad, Marcel Oliver

Abstract: We analyze the method of optimal balance which was introduced by Viudez and Dritschel (J. Fluid Mech. 521, 2004, pp. 343-352) to provide balanced initializations for two-dimensional and three-dimensional geophysical flows, here in the simpler context of a finite dimensional Hamiltonian two-scale system with strong gyroscopic forces. It is well known that when the potential is analytic, such system… ▽ More We analyze the method of optimal balance which was introduced by Viudez and Dritschel (J. Fluid Mech. 521, 2004, pp. 343-352) to provide balanced initializations for two-dimensional and three-dimensional geophysical flows, here in the simpler context of a finite dimensional Hamiltonian two-scale system with strong gyroscopic forces. It is well known that when the potential is analytic, such systems have an approximate slow manifold that is defined up to terms that are exponentially small with respect to the scale separation parameter. The method of optimal balance relies on the observation that the approximate slow manifold remains an adiabatic invariant under slow deformations of the nonlinear interactions. The method is formulated as a boundary value problem for a homotopic deformation of the system from a linear regime, where the slow-fast splitting is known exactly, to the full nonlinear regime. We show that, providing the ramp function which defines the homotopy is of Gevrey class 2 and satisfies vanishing conditions to all orders at the temporal end points, the solution of the optimal balance boundary value problem yields a point on the approximate slow manifold that is exponentially close to the approximation to the slow manifold via exponential asymptotics, albeit with a smaller power of the small parameter in the exponent. In general, the order of accuracy of optimal balance is limited by the order of vanishing derivatives of the ramp function at the temporal end points. We also give a numerical demonstration of the efficacy of optimal balance, showing the dependence of accuracy on the ramp time and the ramp function. △ Less

Submitted 28 August, 2017; originally announced August 2017.

Comments: accepted for publication in SIAM Multiscale Modeling and Simulation

Stochastic partial differential fluid equations as a diffusive limit of deterministic Lagrangian multi-time dynamics

Authors: Colin J Cotter, Georg A Gottwald, Darryl D Holm

Abstract: In {\em{Holm}, Proc. Roy. Soc. A 471 (2015)} stochastic fluid equations were derived by employing a variational principle with an assumed stochastic Lagrangian particle dynamics. Here we show that the same stochastic Lagrangian dynamics naturally arises in a multi-scale decomposition of the deterministic Lagrangian flow map into a slow large-scale mean and a rapidly fluctuating small scale map. We… ▽ More In {\em{Holm}, Proc. Roy. Soc. A 471 (2015)} stochastic fluid equations were derived by employing a variational principle with an assumed stochastic Lagrangian particle dynamics. Here we show that the same stochastic Lagrangian dynamics naturally arises in a multi-scale decomposition of the deterministic Lagrangian flow map into a slow large-scale mean and a rapidly fluctuating small scale map. We employ homogenization theory to derive effective slow stochastic particle dynamics for the resolved mean part, thereby justifying stochastic fluid partial equations in the Eulerian formulation. To justify the application of rigorous homogenization theory, we assume mildly chaotic fast small-scale dynamics, as well as a centering condition. The latter requires that the mean of the fluctuating deviations is small, when pulled back to the mean flow. △ Less

Submitted 1 June, 2017; originally announced June 2017.

Comments: 11 pages, geometric mechanics, stochastic fluid models, stochastic processes, multi-scale fluid dynamics, symmetry reduced variational principles, homogenisation

Comparison of variational balance models for the rotating shallow water equations

Authors: David G. Dritschel, Georg A. Gottwald, Marcel Oliver

Abstract: We present an extensive numerical comparison of a family of balance models appropriate to the semi-geostrophic limit of the rotating shallow water equations, and derived by variational asymptotics in Oliver (2006) for small Rossby numbers ${\mathrm{Ro}}$. This family of generalized large-scale semi-geostrophic (GLSG) models contains the $L_1$-model introduced by Salmon (1983) as a special case. We… ▽ More We present an extensive numerical comparison of a family of balance models appropriate to the semi-geostrophic limit of the rotating shallow water equations, and derived by variational asymptotics in Oliver (2006) for small Rossby numbers ${\mathrm{Ro}}$. This family of generalized large-scale semi-geostrophic (GLSG) models contains the $L_1$-model introduced by Salmon (1983) as a special case. We use these models to produce balanced initial states for the full shallow water equations. We then numerically investigate how well these models capture the dynamics of an initially balanced shallow water flow. It is shown that, whereas the $L_1$-member of the GLSG family is able to reproduce the balanced dynamics of the full shallow water equations on time scales of ${\mathcal{O}}(1/{\mathrm{Ro}})$ very well, all other members develop significant unphysical high wavenumber contributions in the ageostrophic vorticity which spoil the dynamics. △ Less

Submitted 16 May, 2017; originally announced May 2017.

Comments: accepted for publication in Journal of Fluid Mechanics

Stochastic Climate Theory

Authors: Georg A. Gottwald, Daan T. Crommelin, Christian L. E. Franzke

Abstract: 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 averagi… ▽ More 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. △ Less

Submitted 22 December, 2016; originally announced December 2016.

Comments: book chapter in Nonlinear and Stochastic Climate Dynamics, Cambridge University Press (C. Franzke and T. O'Kane, eds.) (2016)

On the detection of superdiffusive behaviour in time series

Authors: Georg A. Gottwald, Ian Melbourne

Abstract: 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 preproces… ▽ More 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. △ Less

Submitted 21 December, 2016; v1 submitted 25 July, 2016; originally announced July 2016.

Journal ref: J Stat Mech (12) (2016) 123205

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

Authors: Georg A. Gottwald, Karsten Peters, Laura Davies

Abstract: 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 fir… ▽ More 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. △ Less

Submitted 16 August, 2015; originally announced August 2015.

Comments: accepted for publication in Quarterly Journal of the Royal Meteorological Society (QJRMS) (2015)

On finite-size Lyapunov exponents in multiscale systems

Authors: Lewis Mitchell, Georg A. Gottwald

Abstract: 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 par… ▽ More 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. △ Less

Submitted 9 April, 2012; originally announced April 2012.

Comments: Accepted for publication in Chaos

Journal ref: Chaos, 22 (2), 023115, 2012

Data assimilation in slow-fast systems using homogenized climate models

Authors: Lewis Mitchell, Georg A. Gottwald

Abstract: 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 m… ▽ More 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. △ Less

Submitted 30 October, 2011; originally announced October 2011.

Comments: Accepted for publication in Journal of the Atmospheric Sciences

ACM Class: J.2

Journal ref: Journal of the Atmospheric Sciences, 69(4), pp. 1359-1377, 2012

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

Authors: Georg A. Gottwald, Lewis Mitchell, Sebastian Reich

Abstract: 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… ▽ More 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. △ Less

Submitted 29 August, 2011; originally announced August 2011.

Comments: 32 pages, 11 figures

Journal ref: Monthly Weather Review, 139(8), pp. 2650-2667, 2011

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

Authors: Sebastian Hermann, Georg A. Gottwald

Abstract: 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 ap… ▽ More 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. △ Less

Submitted 30 March, 2010; originally announced March 2010.

Comments: 32 pages, 17 figures, as accepted by SIAM Journal on Applied Dynamical Systems on 20/03/10

Cuspons, peakons and regular gap solitons between three dispersion curves

Authors: Roger Grimshaw, Boris A. Malomed, Georg A. Gottwald

Abstract: A general wave model with the cubic nonlinearity is introduced to describe a situation when the linear dispersion relation has three branches, which would intersect in the absence of linear couplings between the three waves. Actually, the system contains two waves with a strong linear coupling between them, to which a third wave is then coupled. This model has two gaps in its linear spectrum. Re… ▽ More A general wave model with the cubic nonlinearity is introduced to describe a situation when the linear dispersion relation has three branches, which would intersect in the absence of linear couplings between the three waves. Actually, the system contains two waves with a strong linear coupling between them, to which a third wave is then coupled. This model has two gaps in its linear spectrum. Realizations of this model can be made in terms of temporal or spatial evolution of optical fields in, respectively, a planar waveguide or a bulk-layered medium resembling a photonic-crystal fiber. Another physical system described by the same model is a set of three internal wave modes in a density-stratified fluid. A nonlinear analysis is performed for solitons which have zero velocity in the reference frame in which the group velocity of the third wave vanishes. Disregarding the self-phase modulation (SPM) term in the equation for the third wave, we find two coexisting families of solitons: regular ones, which may be regarded as a smooth deformation of the usual gap solitons in a two-wave system, and cuspons with a singularity in the first derivative at their center. Even in the limit when the linear coupling of the third wave to the first two vanishes, the soliton family remains drastically different from that in the linearly uncoupled system; in this limit, regular solitons whose amplitude exceeds a certain critical value are replaced by peakons. While the regular solitons, cuspons, and peakons are found in an exact analytical form, their stability is tested numerically, which shows that they all may be stable. If the SPM terms are retained, we find that there again coexist two different families of generic stable soliton solutions, namely, regular ones and peakons. △ Less

Submitted 26 March, 2002; originally announced March 2002.

Comments: a latex file with the text and 10 pdf files with figures. Physical Review E, in press