Search

A generic approach to stochastic climate modelling is developed for the example of an idealized Atmosphere-Ocean model that rests upon Hasselmann's paradigm for stochastic climate models. Namely, stochasticity is incorporated into the fast moving atmospheric component of an idealized coupled model by means of stochastic Lie transport, while the slow moving ocean model remains deterministic. More specifically the stochastic model SALT (stochastic advection by Lie transport) is constructed by introducing stochastic transport into the material loop in Kelvin's circulation theorem. The resulting stochastic model preserves circulation, as does the underlying deterministic climate model. A variant of SALT called LA-SALT (Lagrangian-Averaged SALT) is introduced in this paper. In LA-SALT, we replace the drift velocity of the stochastic vector field by its expected value. The remarkable property of LA-SALT is that the evolution of its higher moments are governed by linear deterministic equations. Our modelling approach is substantiated by establishing local existence results, first, for the deterministic climate model that couples compressible atmospheric equations to incompressible ocean equation, and second, for the two stochastic SALT and LA-SALT models., Comment: 51 pages, 1 figure

We present a criterion for uniform in time convergence of the weak error of the Euler scheme for Stochastic Differential equations (SDEs). The criterion requires i) exponential decay in time of the space-derivatives of the semigroup associated with the SDE and ii) bounds on (some) moments of the Euler approximation. We show by means of examples (and counterexamples) how both i) and ii) are needed to obtain the desired result. If the weak error converges to zero uniformly in time, then convergence of ergodic averages follows as well. We also show that Lyapunov-type conditions are neither sufficient nor necessary in order for the weak error of the Euler approximation to converge uniformly in time and clarify relations between the validity of Lyapunov conditions, i) and ii). Conditions for ii) to hold are studied in the literature. Here we produce sufficient conditions for i) to hold. The study of derivative estimates has attracted a lot of attention, however not many results are known in order to guarantee exponentially fast decay of the derivatives. Exponential decay of derivatives typically follows from coercive-type conditions involving the vector fields appearing in the equation and their commutators; here we focus on the case in which such coercive-type conditions are non-uniform in space. To the best of our knowledge, this situation is unexplored in the literature, at least on a systematic level. To obtain results under such space-inhomogeneous conditions we initiate a pathwise approach to the study of derivative estimates for diffusion semigroups and combine this pathwise method with the use of Large Deviation Principles., Comment: 47 pages and 9 figures

We study the long time behaviour of a large class of diffusion processes on $R^N$, generated by second order differential operators of (possibly) degenerate type. The operators that we consider {\em need not} satisfy the H\"ormander condition. Instead, they satisfy the so-called UFG condition, introduced by Herman, Lobry and Sussman in the context of geometric control theory and later by Kusuoka and Stroock, this time with probabilistic motivations. In this paper we study UFG diffusions and demonstrate the importance of such a class of processes in several respects: roughly speaking i) we show that UFG processes constitute a family of SDEs which exhibit multiple invariant measures and for which one is able to describe a systematic procedure to determine the basin of attraction of each invariant measure (equilibrium state). ii) We use an explicit change of coordinates to prove that every UFG diffusion can be, at least locally, represented as a system consisting of an SDE coupled with an ODE, where the ODE evolves independently of the SDE part of the dynamics. iii) As a result, UFG diffusions are inherently "less smooth" than hypoelliptic SDEs; more precisely, we prove that UFG processes do not admit a density with respect to Lebesgue measure on the entire space, but only on suitable time-evolving submanifolds, which we describe. iv) We show that our results and techniques, which we devised for UFG processes, can be applied to the study of the long-time behaviour of non-autonomous hypoelliptic SDEs and therefore produce several results on this latter class of processes as well. v) Because processes that satisfy the (uniform) parabolic H\"ormander condition are UFG processes, our paper contains a wealth of results about the long time behaviour of (uniformly) hypoelliptic processes which are non-ergodic, in the sense that they exhibit multiple invariant measures., Comment: 66 pages

In the late seventies, Clark [In Communication Systems and Random Process Theory (Proc. 2nd NATO Advanced Study Inst., Darlington, 1977) (1978) 721-734, Sijthoff & Noordhoff] pointed out that it would be natural for $\pi_t$, the solution of the stochastic filtering problem, to depend continuously on the observed data $Y=\{Y_s,s\in[0,t]\}$. Indeed, if the signal and the observation noise are independent one can show that, for any suitably chosen test function $f$, there exists a continuous map $\theta^f_t$, defined on the space of continuous paths $C([0,t],\mathbb{R}^d)$ endowed with the uniform convergence topology such that $\pi_t(f)=\theta^f_t(Y)$, almost surely; see, for example, Clark [In Communication Systems and Random Process Theory (Proc. 2nd NATO Advanced Study Inst., Darlington, 1977) (1978) 721-734, Sijthoff & Noordhoff], Clark and Crisan [Probab. Theory Related Fields 133 (2005) 43-56], Davis [Z. Wahrsch. Verw. Gebiete 54 (1980) 125-139], Davis [Teor. Veroyatn. Primen. 27 (1982) 160-167], Kushner [Stochastics 3 (1979) 75-83]. As shown by Davis and Spathopoulos [SIAM J. Control Optim. 25 (1987) 260-278], Davis [In Stochastic Systems: The Mathematics of Filtering and Identification and Applications, Proc. NATO Adv. Study Inst. Les Arcs, Savoie, France 1980 505-528], [In The Oxford Handbook of Nonlinear Filtering (2011) 403-424 Oxford Univ. Press], this type of robust representation is also possible when the signal and the observation noise are correlated, provided the observation process is scalar. For a general correlated noise and multidimensional observations such a representation does not exist. By using the theory of rough paths we provide a solution to this deficiency: the observation process $Y$ is "lifted" to the process $\mathbf{Y}$ that consists of $Y$ and its corresponding L\'{e}vy area process, and we show that there exists a continuous map $\theta_t^f$, defined on a suitably chosen space of H\"{o}lder continuous paths such that $\pi_t(f)=\theta_t^f(\mathbf{Y})$, almost surely., Comment: Published in at http://dx.doi.org/10.1214/12-AAP896 the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org)