Stochastic geometric models with non-stationary spatial correlations in Lagrangian fluid flows

TL;DR

This paper develops data-driven stochastic models for oceanic surface flows that incorporate non-stationary spatial correlations, inspired by satellite observations, and introduces two new models with symmetry-breaking mechanisms.

Contribution

It introduces two novel stochastic models with non-stationary spatial correlations for geophysical fluid dynamics, extending previous models by incorporating flow-advection of correlations.

Findings

Models capture non-stationary spatial correlations in ocean flows

New models derived using stochastic variational principles and Hamiltonian structures

Enhanced understanding of oceanic current dynamics through advanced stochastic modeling

Abstract

Inspired by spatiotemporal observations from satellites of the trajectories of objects drifting near the surface of the ocean in the National Oceanic and Atmospheric Administration's `Global Drifter Program', this paper develops data-driven stochastic models of geophysical fluid dynamics (GFD) with non-stationary spatial correlations representing the dynamical behaviour of oceanic currents. Three models are considered. Model 1 from Holm [2015] is reviewed, in which the spatial correlations are time independent. Two new models, called Model 2 and Model 3, introduce two different symmetry breaking mechanisms by which the spatial correlations may be advected by the flow. These models are derived using reduction by symmetry of stochastic variational principles, leading to stochastic Hamiltonian systems, whose momentum maps, conservation laws and Lie-Poisson bracket structures are used in developing the new stochastic Hamiltonian models of GFD.

Equations

$$d{g}_t(X)=u_t(g_t(X))dt+\xi (g_t(X))\circ dW(t),with{g}_0(X)=X\in \mathcal{D},$$

$$d{g}_t{g}_t^{-1}(x)=u_t(x)dt+\xi (x)\circ dW(t),with{g}_0(X)=X\in \mathcal{D}.$$

$$L(g,v)={\int }_{\mathcal{D}}\frac{1}{2}|v(X){|}^2{d}^{n}X,$$

$$\ell (u)={\int }_{\mathcal{D}}\frac{1}{2}|u(x){|}^2{d}^{n}x,$$

$$q\in V\mapsto g^{\ast }q\in V,$$

$$q\in V\mapsto {\pounds }_{u}q:={\frac{d}{d\epsilon }|}_{\epsilon =0}{g}_{\epsilon }^{\ast }q\in V,$$

$${⟨p,q⟩}_V={\int }_{\mathcal{D}}⟨p(x),q(x)⟩d^{n}x.$$

$$\delta\int_{0}^{T}\Big{[}\ell(u)dt+\left\langle p,dq+\pounds_{dx_{t}}q\right\rangle_{V}\Big{]}=0,$$

$$d{x}_t:=u(t,x)dt+\sum _{i=1}^N{\xi }_i(x)\circ d{W}_i(t),$$

$$\frac{\delta \ell }{\delta u}=p\diamond q,dq+{\pounds }_{d{x}_t}q=0,dp-{\pounds }_{d{x}_t}^{\mathsf{T}}p=0,$$

$${⟨p\diamond q,u⟩}_{\mathfrak{g}}={⟨p,{\pounds }_{u}q⟩}_V,{⟨\frac{\delta \ell }{\delta u},\delta u⟩}_{\mathfrak{g}}:={\frac{d}{d\epsilon }|}_{\epsilon =0}\ell (u+\epsilon \delta u),{⟨{\pounds }_u^{\mathsf{T}}p,q⟩}_V={⟨p,{\pounds }_{u}q⟩}_V,$$

$$d\frac{\delta \ell }{\delta u}+{\mathrm{ad}}_{d{x}_t}^{\ast }\frac{\delta \ell }{\delta u}=0,$$

$${⟨m,u⟩}_{\mathfrak{g}}={\int }_{\mathcal{D}}m(x)\cdot u(x)d^{n}x,$$

$$du+\mathbb{P}(u\cdot \nabla u)dt+\sum _{i=1}^N\mathbb{P}(\mathrm{curl}u\times {\xi }_i)\circ d{W}_i(t)=0.$$

$$d\omega +(d{x}_t\cdot \nabla )\omega -(\omega \cdot \nabla )d{x}_t=0,$$

$${⟨\varpi ,\psi ⟩}_{\mathfrak{g}}={\int }_{\mathcal{D}}\varpi (x)\psi (x)d^{2}x.$$

$$d\omega +\{\omega ,\psi \}dt+\sum _{i=1}^N\{\omega ,{\psi }_i\}\circ d{W}_i(t)=0,$$

$$\delta {\int }_0^{T}L(g,\dot{g})dt=0,$$

$$\delta\int_{0}^{T}\Big{[}L(g,v)dt+\big{\langle}\pi,dg-vdt-\sum_{i=1}^{N}\xi_{i}g\circ dW_{i}(t)\big{\rangle}\Big{]}=0\,,$$

$$\delta\int_{0}^{T}\Big{[}\ell(u)dt+\big{\langle}m,dgg^{-1}-udt-\sum_{i=1}^{N}\xi_{i}\circ dW_{i}(t)\big{\rangle}_{\mathfrak{g}}\Big{]}=0\,,$$

$$\delta\int_{0}^{T}\Big{[}L(g,v)+\big{\langle}\pi,dg-vdt\big{\rangle}-\sum_{i=1}^{N}H_{i}(g,\pi;\xi_{i})\circ dW_{i}(t)\Big{]}=0\,,$$

$$H_i(g,\pi ;{\xi }_i):=⟨\pi ,{\xi }_{i}g⟩={⟨\pi g^{-1},{\xi }_i⟩}_{\mathfrak{g}},i=1,...,N.$$

$$d\frac{\partial L}{\partial v}=\frac{\partial L}{\partial g}dt-\sum _{i=1}^N\frac{\partial H_i}{\partial g}\circ d{W}_i(t)=0,dg=vdt+\frac{\partial H_i}{\partial \pi }\circ d{W}_i(t),\pi =\frac{\partial L}{\partial v}.$$

$$dg=\frac{\partial H}{\partial \pi }dt+\sum _{i=1}^N\frac{\partial H_i}{\partial \pi }\circ d{W}_i(t),d\pi =-\frac{\partial H}{\partial g}dt-\sum _{i=1}^N\frac{\partial H_i}{\partial g}\circ d{W}_i(t).$$

$$dF=\{F,H{\}}_{\mathrm{can}}dt+\sum _{i=1}^N\{F,H_i{\}}_{\mathrm{can}}\circ d{W}_i(t),$$

$$dm+{\mathrm{ad}}_{\frac{\delta h}{\delta m}}^{\ast }mdt+\sum _{i=1}^N{\mathrm{ad}}_{\frac{\delta h_i}{\delta m}}^{\ast }m\circ d{W}_i(t)=0,$$

$$h(m)={\int }_{\mathcal{D}}\frac{1}{2}|m(x){|}^2{d}^{n}x\text{and}{h}_i(m)={⟨m,{\xi }_i⟩}_{\mathfrak{g}}={\int }_{\mathcal{D}}m(x)\cdot {\xi }_i(x)d^{n}x.$$

$$\{f,h{\}}_{\mathrm{LP}}(m)={⟨m,\left[\frac{\delta f}{\delta m},\frac{\delta h}{\delta m}\right]⟩}_{\mathfrak{g}},$$

$$df=\{f,h{\}}_{\mathrm{LP}}dt+\sum _{i=1}^N\{f,h_i{\}}_{\mathrm{LP}}\circ d{W}_i(t),$$

$$\{f,g{\}}_{\mathrm{LP}}(\omega )={\int }_{\mathcal{D}}\omega \left\{\frac{\delta f}{\delta \omega },\frac{\delta g}{\delta \omega }\right\}{d}^{2}x$$

Full text

1. 1 Introduction
2. 2 Model 1: Review of stochastic variational principles for fluids
3. 3 Model 2: Frozen-in correlations for non-stationary statistics
4. 4 Model 3: Eigenvectors depending on advected quantities
5. 5 Conclusions

1 Introduction

This paper develops data-driven stochastic models of fluid dynamics, inspired by spatiotemporal observations from satellites of the spatial paths of objects drifting near the surface of the ocean in the National Oceanic and Atmospheric Administration’s “Global Drifter Program”. The Lagrangian paths of these freely drifting instruments track the ocean currents. That is, the satellite readings of their positions approximate the motion of a fluid parcel as a curve parameterised by time.

Figure 1.1 Lilly [2017] displays the global array of surface drifter trajectories from the National Oceanic and Atmospheric Administration’s “Global Drifter Program” (www.aoml.noaa.gov/phod/dac). In total, more than 10,000 drifters have been deployed since 1979, representing nearly 30 million data points of positions along the Lagrangian paths of the drifters at six-hour intervals. This large spatiotemporal data set is a major source of information regarding ocean circulation, which in turn is an important component of the global climate system; see for example Lumpkin and Pazos [2007]; Griffa et al. [2007]; Sykulski et al. [2016]. An important feature of this data is that the ocean currents show up as time-varying spatial correlations, easily recognised visually by the concentrations of colours representing individual paths. These spatial correlations exhibit a variety of spatial scales for the trajectories of the drifters, corresponding to the variety of spatiotemporal scales in the evolution of the ocean currents which transport the drifters.

Figure 1.2 shows a sample of the Lagrangian trajectories of drifters released in the vicinity of Cape Cod. Here one sees concentrations of recirculating drifters following Western boundary currents, and splitting into three main streams just off the coast of Cape Cod. One of these main streams forms a large-scale re-circulation some distance away from the boundary and to the East of Cape Cod. Small-scale, erratic deviations of the main streams of the drifter paths are also visible in this Figure. These erratic trajectories will be represented stochastically in the paper, and the larger scale spatial correlations they follow will be modelled as spatiotemporal modulations of the stochasticity, which follow the resolved drift currents.

Figure 1.3 displays the Lagrangian path taken by Drifter 36256, deployed on Sep 18, 2005, and successfully recovered on Feb 21, 2007, at Brest, France, after a journey of 521 days across the North Atlantic Ocean. Along this path one sees the effects of the interactions of the drifter with a variety of space and time scales of of evolving fluid motion, strongly suggesting a need for modelling the non-stationary statistics of Lagrangian paths Sykulski et al. [2016].

Why introduce stochasticity into fluid dynamics?

In developing parameterisations in weather and climate prediction, the term “unresolved scales” refers to those fluid motions and thermodynamic properties which are either not computable in a given set of numerical simulations, or are not measurable in a given set of observations. The effects of the smaller, faster, unresolved scales on the larger, slower, resolved variables are often modelled deterministically, in terms of the resolved scales and their gradients. See, e.g., Foias et al. [2001, 2002]; Holm [2002] and references therein. Although this deterministic approach may be pragmatic and even advisable, due to the limitations of simulations and measurements, it is clear from the variety of deterministic models which have been proposed that these models are by no means unique. For example, the plethora of turbulence models and averaging methods, as well as the stark differences in solution behaviour as simulations of the fluid equations are achieved at finer and finer resolutions, show that the resolved states may be associated with many possible unresolved states. For example, phenomena such as vortex tubes of a certain radius, as seen in simulations at one scale, or resolution, may no longer even be present when simulated at smaller scales, achieved with better resolution. This is not the usual situation in the development of science. Usually, the advent of better, more accurate, measurements lead to new effects and new laws governing them. Instead, in the case of weather and climate prediction, the fundamental deterministic equations and the laws of thermodynamics are completely known.

No new deterministic laws of fluid motion should be expected. However, there can be new statistical approaches to the physical and mathematical descriptions of weather and climate. For example, the recognition of stochastic vector fields as the basic paradigm in fluid dynamics has long been the province of turbulence modelling Monin and Yaglom [1971]. However, recently the use of stochastic vector fields has also been recognized in estimating statistical model uncertainty in numerical weather prediction, Franzke et al. [2015]. From this viewpoint, the uncertainty and variability of the predictions are crucial aspects of the solution. In the statistical science of numerical weather and climate prediction, stochastic methods offer systematic approaches toward quantitative estimates of uncertainties due to model error and inaccuracy of data assimilation, as well as improved estimates of long-term climate variability, including estimates of the probability of extreme events. Following Berner et al. [2013] we agree that “stochasticity must be incorporated at a very basic level within the design of physical process parameterizations and improvements to the dynamical core.” The purpose of this paper is to offer new approaches at this basic level. However, it is beyond the scope of the present paper to analyse the large dataset of drifter trajectories that has inspired our investigation. The work is in progress to fully analyse the drifter dataset and will be discussed elsewhere. Here we will discuss both 2D and 3D stochastic fluid dynamics models.

How to do it?

How does one use stochasticity to improve the physical and mathematical basis for designing statistical model uncertainty schemes? Recently, this question was answered by using Hamilton’s principle to derive a new class of mathematical models of stochastic transport in fluid dynamics Holm [2015]. In this class of models, the effects of the small, fast, unresolved fine scales of motion on the coarser ones are modelled by introducing stochastic uncertainty into the transport velocity of fluid parcels in the dynamics at the resolvable coarse scale. This stochastic transport velocity decomposition encompasses both Newtonian and variational perspectives of mechanics and also leads to proper Kelvin circulation dynamics. For a rigorous derivation of the same decomposition into mean and fluctuating velocities using multi-time homogenisation methods, see Cotter et al. [2017]. Applying this decomposition to create stochastic fluid dynamics preserves the fundamental mathematical properties of their deterministic counterparts Crisan et al. [2017]. It also enables new approaches to sub-grid scale parameterization, expressed both in terms of fluctuation distributions, and spatial/temporal correlations. As such, it introduces stochastic corrections that are amenable to statistical inference from high-resolution data (either observed, or numerical). Moreover, this new class of models forms the ideal approach for the development of a novel data assimilation technology based on particle filters Beskos et al. [2017]. Filtering and ensemble techniques require de facto a stochastic representation of the dynamics. This randomization is most often achieved through random perturbations of the initial conditions. However, this approach tends to yield insufficient spreading of the ensemble and produce a poor representation of the error dynamics Berner et al. [2013]. The stochastic transport class of models establishes the much-needed randomization via its rigorous derivation at the fundamental level, rather than via ad hoc empiricism.

In this paper, we will introduce two different stochastic extensions of Holm [2015] in applying the geometric mechanics framework to estimate the contribution of stochastic transport to statistical uncertainty (error) in fluid dynamics models.

Stochastic transport means that the Lagrangian fluid parcel motion has a stochastic component. In this context, we ask the following question in the context of numerical weather and climate prediction: What fundamental properties of the deterministic fluid equations would persist in a stochastic vector field representation of continuum fluid motions? First, even if the fluid parcel velocity were stochastic, the fluid continuum motion would still be describable as a spatially smooth but now temporally stochastic flow, $g_{t}$, depending on time $t$. For the sake of brevity and simplicity of the presentation, we will take the domain of flow $\mathcal{D}$ to lie in $\mathbb{R}^{2}$ or $\mathbb{R}^{3}$, and we will neglect considerations of boundary conditions in the examples discussed in this paper. Even in the presence of stochasticity, the stochastic path of the fluid parcel which is initially at position $X$ in the domain of flow can still be represented by the formula for the Lagrange-to-Euler flow map, $x_{t}=g_{t}(X)$, so that $g_{0}$ at time $t=0$ is the identity map, $g_{0}(X)=X$. Since the stochasticity is Markovian, the flow map still corresponds to a stochastic time-dependent curve $g_{t}$ on the group of compositions of smooth invertible maps, i.e., the diffeomorphisms, acting on flow domain $\mathcal{D}$, see Crisan et al. [2017].

Thus, following Kraichnan [1994], Mikulevicius and Rozovskii [2004], and Holm [2015], we may begin by assuming that the stochastic paths $x_{t}=g_{t}(X)$ solve a Lagrangian stochastic differential equation (SDE) with prescribed spatially dependent function $\xi_{t}=\xi(g_{t}(X))$

[TABLE]

where $g_{t}:\mathcal{D}\to\mathcal{D}$ is a spatially smooth invertible map depending on time, $t\in\mathbb{R}$. The corresponding Eulerian stochastic velocity decomposition is given in terms of cylindrical noise, introduced in Schaumlöffel [1988], as

[TABLE]

In this approach, the particle-and-field duality of the Lagrangian and Eulerian descriptions of continuum fluid motion is still available, even through the fluid parcel motions are stochastic. This means the symmetry of the Eulerian description under relabelling of Lagrangian fluid parcels persists, even when the Lagrangian fluid paths are stochastic, so that the Lagrange-to-Euler fluid map is also stochastic. This relabelling symmetry implies the Kelvin circulation theorems for fluid dynamics models Holm et al. [1998]. Of course, these remarks generalise to any number of dimensions.

What does this paper do?

Given the stochasticity in the flow map, how does one derive the corresponding Eulerian equations of continuum motion? The paper lays out a geometric framework for deriving stochastic Eulerian motion equations for geophysical fluid dynamics (GFD) using the method of symmetry reduction for a modified Hamilton’s principle for fluids with advected quantities.

Although the paper is based on previous developments in Kraichnan [1994], Mikulevicius and Rozovskii [2004], and Holm [2015], it introduces two new approaches for incorporating non-stationary statistics due to flow dependence, as seen in the NOAA drifter data shown in Figures 1.1, 1.2 and 1.3, and analysed in Sykulski et al. [2016].

In particular, the paper allows for flow dependence of the eigenvectors $\xi_{t}$ of the spatial correlations in the stochastic process in (1.2). To obtain this flow dependence, we postulate two different approaches for allowing the eigenvectors for the stochastic process to evolve along with the advected quantities. These two models are each modifications of the approach in Holm [2015]. The two models provide alternative systematic avenues for forecasting with evolving time-dependent statistics, such as those seen in the NOAA drifter data analysed in Sykulski et al. [2016], rather than using spatially dependent but steady statistics. While steady statistics may be appropriate for climate science, analysis of flows at shorter time scales in forecasting weather variability, for example, may require flow-dependent, evolving statistics. Our purpose here is to present a systematic framework for modelling non-stationary statistics in stochastic fluid flows.

Plan of the paper.

In section 2, we quickly review the approach of Holm [2015] as it applies to the Euler equations of a perfect fluid. The methodology of this approach has two primary features: (1) stochastic variational principles; and (2) stochastic Hamiltonian formulations. These two primary features will allow us to introduce two new stochastic extensions of Geophysical Fluid Dynamics (GFD), one with advection by the drift velocity of the eigenvectors $\xi_{t}$ discussed in section 3, and the other with eigenvectors $\xi_{t}$ depending on advected fluid quantities discussed in section 4. Both Stratonovich and Itô forms of the equations are provided. We conclude by comparing the three stochastic models in terms of their Kelvin circulation theorems in section 5.

2 Model 1: Review of stochastic variational principles for fluids

Following the geometric approach of Arnold [1966], we consider the group $G$ of volume preserving diffeomorphisms of the fluid domain $\mathcal{D}$, as the configuration manifold for incompressible fluids. Curves $g_{t}\in G$ in this group describe Lagrangian trajectories $x_{t}=g_{t}(X)$ of the fluid motion. To simplify our discussion, we will take the domain $\mathcal{D}$ to lie in $\mathbb{R}^{2}$ or $\mathbb{R}^{3}$ and neglect considerations of boundary conditions. Our developments extend easily to the case where $\mathcal{D}$ is a manifold with smooth boundary.

The Lagrangian of the incompressible fluid is defined on the tangent bundle $TG$ of the group $G$ and is given by the kinetic energy, i.e.,

[TABLE]

for $n=2,3$. By a change of variables, we note that $L$ is right invariant: $L(gh,vh)=L(g,v)$, for all $h$ in $G$. We can thus write $L(g,v)=\ell(vg^{-1})$, where $\ell:\mathfrak{g}\rightarrow\mathbb{R}$. That is, we write

[TABLE]

as the reduced Lagrangian defined on the Lie algebra $\mathfrak{g}$ of $G$, given by the space of divergence free vector fields, denoted $u:=vg^{-1}\in\mathfrak{g}$.

Hamilton’s principle $\delta\int_{0}^{T}L(g,\dot{g})dt=0$ yields the Euler equations in Lagrangian description, i.e., geodesics on $G$, whereas the variational principle induced on $\ell$, the Euler-Poincaré principle, yields the Euler equations in their standard spatial description.

In particular, the stochastic variational principle in Holm [2015] is obtained by selecting a space $V$ of tensor fields on $\mathcal{D}$, denoted $q(t,x)$, on which the group $G$ acts linearly by the pull-back operation

[TABLE]

which is the natural transport operation of tensor fields by the fluid motion. The associated Lie algebra action of divergence free vector fields $u\in\mathfrak{g}$ is given by the Lie derivative as

[TABLE]

where $g_{\varepsilon}$ is the flow of $u$. We fix a space $V^{\ast}$ of tensor fields in nondegenerate duality with $V$ relative to the $L^{2}$ pairing

[TABLE]

Remark 2.1** (Determining the correlation eigenvectors $\xi_{i}(x)$).**

How may the eigenvectors $\xi_{i}(x)$ be found in practice? One approach for determining them would be based on running an ensemble of computer simulations of the flow started with observed initial velocity data measured in the NOAA data base at each grid point $X$ on a computational mesh erected in the flow domain. This computational mesh would be significantly coarser that the resolution of the data. An ensemble of these computer runs would produce a distribution of simulated end points $X^{r}=U^{r}\,dt$ on the coarse mesh, where $r$ labels the ensemble member. The distribution of differences $\Delta X^{r}$ of the simulations from the observed data around each coarse grid point would be modelled (up to some specified tolerance) as a local stochastic process given by $dX=\sum_{j=1}^{N}\xi_{j}(X)\circ dW_{j}(t)$, where $\{\xi_{j}(X)\}$ is a set of orthogonal functions whose squared amplitudes are given by $\lambda_{j}$ (kinetic energy) for the spectrum of the correlation tensor $C(X,Y)=\xi(\{X\})\otimes\xi^{T}(\{Y\})$ where vectors $\{X\}$ and $\{Y\}$ represent all of the points on the coarse grid. Finally, the choice of tolerance determines the choice of $N$. See Cotter et al. [2017] for an application of this stochastic data assimilation approach, in which the “observed data”, or “truth”, is given by a highly resolved numerical simulation. In practice, the noise correlation eigenvectors must be inferred from observed data, on a case by case basis which will depend on the quality and completeness of the data. However, the inference of velocity-velocity correlation tensors is standard practice in fluid turbulence experiments, cf. Berkooz et al. [1993].

Given $N$ time independent divergence free vector fields $\xi_{i}(x)$, $i=1,...,N$, the stochastic variational principle in Holm [2015] is formally written as

[TABLE]

with respect to variations $\delta u$, $\delta q$, $\delta p$ and where $dx_{t}$ is defined as

[TABLE]

in which the vector fields $\xi_{i}(x)$ represent spatial correlations of the stochasticity and $dW_{i}(t)$ are independent Brownian motions, introduced in the Stratonovich sense. The stationarity conditions are computed by fixing a space $\mathfrak{g}^{\ast}$ in nondegenerate duality with $\mathfrak{g}$, relative to a pairing $\left\langle m,u\right\rangle_{\mathfrak{g}}$, $u\in\mathfrak{g}$, $m\in\mathfrak{g}^{\ast}$. The variations of (2.3) in $\delta u$, $\delta p$ and $\delta q$ yield, respectively, the conditions

[TABLE]

where $p\diamond q\in\mathfrak{g}^{\ast},\frac{\delta\ell}{\delta u}\in\mathfrak{g}^{\ast}$, and $\pounds^{\mathsf{T}}_{u}p\in V^{\ast}$ are defined as

[TABLE]

for $q\in V$, $p\in V^{\ast}$, and $u,\delta u\in\mathfrak{g}$. The conditions (2.5) imply the following stochastic equation:

[TABLE]

where $\operatorname{ad}^{*}_{u}:\mathfrak{g}^{\ast}\rightarrow\mathfrak{g}^{\ast}$ denotes the coadjoint operator defined by $\left\langle\operatorname{ad}^{*}_{u}m,v\right\rangle_{\mathfrak{g}}=\left\langle m,[u,v]\right\rangle_{\mathfrak{g}}$, with $[u,v]=v\cdot\nabla u-u\cdot\nabla v$, and $dx_{t}$ is given in (2.4).

The notations used in (2.7) are general enough to make this equation valid for any Lie group $G$ and Lagrangian $\ell:\mathfrak{g}\rightarrow\mathbb{R}$. See Arnaudon et al. [2017a] for a parallel treatment of Model 1 for the rigid body and the group $SO(3)$, as well as for the heavy top, which involves advected quantities arising from symmetry breaking from $SO(3)$ to $SO(2)$.

Upon choosing for $\mathfrak{g}^{\ast}$ the space of divergence free vector fields on $\mathcal{D}$, i.e., $\mathfrak{g}^{\ast}=\mathfrak{g}$, and the duality pairing

[TABLE]

the coadjoint operator is $\operatorname{ad}^{*}_{u}m=\mathbb{P}(u\cdot\nabla m+\nabla u^{\mathsf{T}}\cdot m)$, where $\mathbb{P}$ is the Hodge projection onto divergence free vector fields. With the Lagrangian (2.1), the stochastic Euler equation (2.7) becomes, in 3D,

[TABLE]

Equation (2.8) can be written equivalently in vorticity form as

[TABLE]

where $\omega={\rm curl}\,u$ is the vorticity and the stochastic vector field ${dx_{t}}$ is given in equation (2.4).

In 2D we identify the space of divergence free vector fields with the space of functions on $\mathcal{D}$ modulo constants, the stream functions $\psi$. As explained in Marsden and Weinstein [1983], the dual space $\mathfrak{g}^{\ast}$ is identified with functions on $\mathcal{D}$ with zero integrals, the absolute vorticities $\varpi$, via the duality pairing

[TABLE]

The absolute vorticity $\varpi$ is related to the total fluid momentum $m$ as $\varpi=\operatorname{curl}m\cdot\mathbf{z}$. For instance, for the Euler equation, the absolute vorticity coincide with the vorticity $\omega=\operatorname{curl}u\cdot\mathbf{z}$, whereas for the rotating Euler equation, we have $\varpi=\operatorname{curl}u\cdot\mathbf{z}+f=\omega+f$, where $f$ is Earth’s frequency rotation.

In 2D, the stochastic Euler equation (2.7) becomes

[TABLE]

where, for two functions $f,g$ on $\mathcal{D}$, the function $\{f,g\}$ is the Jacobian defined by $\{f,g\}:=\partial_{x_{1}}f\partial_{x_{2}}g-\partial_{x_{2}}f\partial_{x_{1}}g$, with $x=(x_{1},x_{2})$. In (2.10), $\psi(t,x)$ is the stream function of the fluid velocity $u(t,x)$, $\omega(t,x)=-\Delta\psi(t,x)$ is its vorticity, and the functions $\psi_{i}(x)$ are the stream functions of the divergence free vector fields $\xi_{i}(x)$. The deterministic Euler equations are recovered in (2.8) and (2.10) when $\xi_{i}=0$, for all $=1,...,N$.

We shall now make two crucial observations about this approach, that will allow us later to extend the approach further to develop the other two stochastic models treated in the paper.

Observation 1: Stochastic variational principles.

Knowing that the deterministic Euler equations in the Lagrangian fluid description arise from the Hamilton principle

[TABLE]

for the right invariant Lagrangian $L:TG\rightarrow\mathbb{R}$ given by the kinetic energy, we expect the equations (2.7) to arise, in the Lagrangian description, via a stochastic extension of Hamilton’s principle (2.11). This is indeed the case if one proceeds formally here and below by considering the stochastic Hamilton-Pontryagin (SHP) principle444This is a variational principle on the Pontryagin bundle $TG\oplus T^{*}G\rightarrow G$, defined as the vector bundle over $G$ with vector fiber at $g\in G$ given by $T_{g}G\oplus T^{*}_{g}G$.

[TABLE]

for variations $\delta g$, $\delta v$, $\delta\pi$. The variables $v$ and $\pi$ are, respectively, the material fluid velocity and material fluid momentum. In (2.12), $\left\langle\,\cdot\,,\,\cdot\,\right\rangle$ denotes the pairing between elements in $T_{g}^{\ast}G$, and $T_{g}G$, the cotangent and tangent space to $G$ at $g$.

The notation $\xi_{i}g$ indicates the composition of the vector field $\xi_{i}$ on the right by the diffeomorphism $g$.

Stochastic Hamilton-Pontryagin principles (SHP) have been considered for finite dimensions in Bou-Rabee and Owhadi [2009]. The present paper considers SHP in infinite dimensions for the first time. The SHP affords a systematic derivation of the stochastic equations that preserves their deterministic mathematical properties, both geometrical and analytical.

In the present paper, we shall consider stochastic variational principles in infinite dimensions only in a formal sense for the purpose of modelling time-dependent spatial correlations. The corresponding questions in analysis, for example, the questions of local in time existence and uniqueness of solutions answered in Crisan et al. [2017] for the stochastic 3D Euler fluid model, will all be left open for the two new stochastic geometric fluid models that are introduced in this paper.

Note that (2.12) imposes the stochastic process (1.1) as a constraint on the variations by using the Lagrange multiplier $\pi$. From the $G$-invariance of both the Lagrangian and the constraint, this principle can be equivalently written formally in the reduced Eulerian description as

[TABLE]

with respect to variations $\delta u,\delta g,\delta m$, and where $u=vg^{-1}\in\mathfrak{g}$, $m=\pi g^{-1}\in\mathfrak{g}^{\ast}$. This is the reduced stochastic Hamilton-Pontryagin (RSHP) principle.

One then directly checks that the stochastic variational principle (2.13) also yields the stochastic equation (2.7). Thus, the two variational principles (2.3) and (2.13) both yield the same stochastic equations. Moreover, in absence of stochasticity, equation (2.12) recovers the Hamilton-Pontryagin principle for Lagrangian mechanics, see Yoshimura and Marsden [2006].

Remark 2.2**.**

The RSHP principle in (2.13) has several interesting properties: (i) it allows a formulation of reduction by symmetry in the stochastic context; (ii) it does not need the introduction of the extra advected quantities $q,p$; and (iii) it does not restrict the values of the Eulerian fluid momentum $m\in\mathfrak{g}^{\ast}$ to be of the form, $m=p\,\diamond\,q$. In addition, as we will show later, the unreduced SHP principle (2.12) allows us to consistently implement the new Model 2 and Model 3, in which the spatial correlation eigenvectors $\xi_{i}(x)$ which are fixed functions of the spatial coordinates in Model 1 become time dependent through their flow dependence in Model 2 and Model 3.**

Observation 2: Stochastic Hamiltonian formulations.

We note that the SHP principle (2.12) can be equivalently written as

[TABLE]

for the $G$-invariant functions $H_{i}(\_\,,\_\,;\xi_{i}):T^{*}G\rightarrow\mathbb{R}$ defined by

[TABLE]

The principle in (2.14) yields the following stochastic extension of the Euler-Lagrange equations with Lagrangian $L$:

[TABLE]

This is the Lagrangian description of the stochastic equations (2.7). Denoting by $H:T^{*}G\rightarrow\mathbb{R}$, the Hamiltonian associated to $L$ by the Legendre transform, we can rewrite these equations in stochastic Hamiltonian form

[TABLE]

Consequently, we can call the functions $H_{i}$ the stochastic Hamiltonians. Stochastic Hamiltonian systems of the form (2.17) have been developed in Bismut [1982]. The intrinsic form of equations (2.16) and (2.17) on the Lie group $G$ would require introducing a covariant derivative. The formulation above is only valid locally.

These equations can be written in terms of the canonical Poisson bracket $\{\,\cdot\,,\,\cdot\,\}_{\rm can}$ on $T^{*}G$ as

[TABLE]

for arbitrary functionals $F=F(g,\pi):T^{*}G\rightarrow\mathbb{R}$.

Consistently with this observation, we note that the stochastic equation (2.7) can also be written in Hamiltonian form as

[TABLE]

where $h:\mathfrak{g}^{\ast}\rightarrow\mathbb{R}$ and $h_{i}:\mathfrak{g}^{\ast}\rightarrow\mathbb{R}$ are the reduced Hamiltonians associated to $H$ and $H_{i}$ in (2.15), i.e., $H(g,\pi)=h(\pi g^{-1})$ and $H_{i}(g,\pi;\xi_{i})=h_{i}(\pi g^{-1})$. We find

[TABLE]

The expression (2.19) is the reduced (or Eulerian) formulation of the Hamiltonian formulation (2.17).

In terms of the Lie-Poisson bracket $\{\,,\}_{\rm LP}$ on $\mathfrak{g}^{\ast}$, given by

[TABLE]

equation (2.7) and hence (2.19) can be formulated in the Stratonovich-Lie-Poisson form

[TABLE]

for arbitrary functions $f:\mathfrak{g}^{\ast}\rightarrow\mathbb{R}$, which is just the reduced form of (2.18).

For example, the stochastic 2D Euler equations (2.10), can be written in the Stratonovich-Lie-Poisson form (2.21) with the Lie-Poisson bracket written on the space of vorticities as Marsden and Weinstein [1983]

[TABLE]

and with the stochastic Hamiltonians

[TABLE]

Aim of the paper.

By exploiting the two observations discussed above, this paper will develop two new stochastic models by appropriate modifications of the stochastic Hamiltonians $H_{i}$ and of their symmetries.

3 Model 2: Frozen-in correlations for non-stationary statistics

As discussed in the Introduction, we consider a model in which the eigenvectors $\xi_{i}(X)$ are advected by the flow map $g_{t}$, giving the time dependent vector fields $\zeta_{i}(t,x)$. Mathematically, the advection of a vector field by a smooth invertible map $g_{t}$ corresponds to the push-forward operation,

[TABLE]

meaning that

[TABLE]

From a Lie group point of view, the push-forward is given by the adjoint action Ad$:G\times\mathfrak{g}\to\mathfrak{g}$ of the Lie group $G$ on its Lie algebra $\mathfrak{g}$, that is,

[TABLE]

The notation used in this section is general enough to make our developments valid for any Lie group $G$. In our examples, we will take $G$ to be the either the group of diffeomorphisms or the group of volume preserving diffeomorphisms. The corresponding coadjoint operators are

[TABLE]

The stochastic model considered here for advection of the eigenvectors is obtained by modifying the expression and the symmetries of the stochastic Hamiltonians $H_{i}(g,\pi;\xi_{i})$ in (2.14). Namely, given $N$ vector fields $\xi_{i}$, $i=1,...,N$, we consider the principle

[TABLE]

where each of the stochastic Hamiltonians $H_{i}(\_\,,\_\,;\xi_{i}):T^{*}G\rightarrow\mathbb{R}$ is right invariant only under the action of the isotropy subgroup of the eigenvector $\xi_{i}$ with respect to the adjoint action (i.e., the push-forward action), namely,

[TABLE]

That is, we have $H_{i}(gh,\pi h;\xi_{i})=H_{i}(g,\pi;\xi_{i})$, for all $h\in G_{\xi_{i}}$. The SHP principle (3.3) yields equations in the same general form as (2.17), but with stochastic Hamiltonians which are not $G$-invariant. As we will show, this difference due to symmetry breaking from $G$ to $G_{\xi_{i}}$ induces significant changes in the reduced Eulerian representation.

Physically, the symmetry breaking means that the initial conditions for the correlation eigenvectors are “frozen” into the subsequent flow, as a property carried along with individual Lagrangian fluid parcels, and which is not exchanged with other fluid parcels.

Being only $G_{\xi_{i}}$-invariant, the stochastic Hamiltonian function $H_{i}$ induces, in the Eulerian description, the reduced stochastic Hamiltonians

[TABLE]

defined by

[TABLE]

where $\mathcal{O}_{\xi_{i}}:=\{\operatorname{Ad}_{g}\xi_{i}\mid g\in G\}=\{g_{\ast}\xi_{i}\mid g\in G\}\subset\mathfrak{g}$ is the adjoint orbit of $\xi_{i}$. The SHP principle (3.3) can thus be written in the reduced Eulerian form as

[TABLE]

The stationarity conditions with respect to the variations $\delta u$, $\delta m$, $\delta g$, yield

[TABLE]

where the advected eigenvector $\zeta_{i}:=\operatorname{Ad}_{g}\xi_{i}$ obeys the auxiliary equation $d\zeta_{i}=[dgg^{-1},\zeta_{i}]$, obtained from its definition. The expression of the coadjoint operator $\operatorname{ad}^{\ast}$ is given in (3.2). Upon using the second equation in (3.7), this auxiliary equation becomes

[TABLE]

Remark 3.1**.**

In our applications of this model, we will always assume that the stochastic Hamiltonians $h_{i}$ in (3.8) do not depend on $m$. That is, $\delta h_{j}/\delta m=0$, so that the eigenvectors $\zeta_{i}$ are advected only by the drift velocity, $u$. That is,

[TABLE]

This model using correlation eigenvectors frozen into the drift velocity is quite different from the model in Holm [2015], which chose $h_{i}(m)=\langle m\,,\,\xi_{i}\rangle_{\mathfrak{g}}$ and all fluid properties were advected by a velocity vector field comprising the sum of both the drift component and the stochastic component.

Hamiltonian structure.

Denoting by $h:\mathfrak{g}^{\ast}\rightarrow\mathbb{R}$ the Hamiltonian associated to $\ell$, the above system can be equivalently written as

[TABLE]

One may check that this system is Hamiltonian with respect to the Poisson bracket

[TABLE]

on $\mathfrak{g}^{\ast}\times\mathcal{O}_{\xi}\ni(m,\zeta_{1},...,\zeta_{N})$, where $\mathcal{O}_{\xi}\subset\mathfrak{g}^{N}$ is the orbit of $\xi:=(\xi_{1},...,\xi_{N})$ under the adjoint action of $G$. This reduced Poisson bracket is inherited from Poisson reduction of the canonical Poisson bracket $\{\,\cdot\,,\,\cdot\,\}_{\rm can}$ on $T^{*}G$, by the isotropy subgroup $G_{\xi}\subset G$ of $\xi:=(\xi_{1},....,\xi_{N})$. Namely, the map

[TABLE]

is Poisson with respect to the Poisson brackets $\{\,\cdot\,,\,\cdot\,\}_{\rm can}$ and $\{\,\cdot\,,\,\cdot\,\}_{\rm red}$, see the Appendix. The system (3.10) admits the Stratonovich-Poisson formulation,

[TABLE]

for arbitrary functions $f:\mathfrak{g}^{\ast}\times\mathcal{O}\rightarrow\mathbb{R}$. Note that in (3.12) the Hamiltonians $h$ and $h_{i}$, $i=1,...,N$, depend a priori on all the variables $(m,\zeta_{1},...,\zeta_{N})$. The system (3.10) is recovered when $h$ depends only on $m$, while the $h_{i}$ depend only on $m$ and $\zeta_{i}$ (not on $\zeta_{j}$, for $j\neq i$). The Poisson tensor at $(m,\zeta_{1},...,\zeta_{N})$ reads

[TABLE]

Remark 3.2** (Itô form).**

In the special case that ${\delta h_{i}}/{\delta m}=0$, as assumed in Remark 3.1, the Itô form of equations (3.10) does not introduce any additional drift terms. That is, in this special case, the Itô form of (3.10) is obtained by simply removing the Stratonovich symbol $(\,\circ\,)$. However, in the general case, if ${\delta h_{i}}/{\delta m}\neq 0$, the Itô form does contain additional drift terms. These additional drift terms in the Itô form for the general case can be computed in a standard way, but the equations then may take a more complicated form. By making use of the Poisson formulation in terms of the bracket $\{\,\cdot\,,\,\cdot\,\}_{\rm red}$ in (3.12), we can write the additional drift terms in the Itô form for the general case in a concise way as

[TABLE]

Example 3.3** (Incompressible 2D models).**

In the 2D incompressible case, we can identify the Lie algebra $\mathfrak{g}$ with the space of differentiable functions on $\mathcal{D}$, modulo constants. These are the stream functions, denoted $\psi$. We use the $L^{2}$ duality pairing $\left\langle\varpi,\psi\right\rangle_{\mathfrak{g}}=\int_{\mathcal{D}}\varpi(x)\psi(x)d^{2}x$ and identify $\mathfrak{g}^{\ast}$ with the space of functions on $\mathcal{D}$ with zero integral. These are the absolute vorticities, denoted $\varpi$, as explained in Marsden and Weinstein [1983].

Let $\psi_{i}^{0}(X)$ be the stream function associated to the eigenvector $\xi_{i}(X)$, $i=1,...,N$. The stream function $\psi_{i}(t,x)$ of the advected eigenvector $\zeta_{i}(t)=(g_{t})_{\ast}\xi_{i}$ is found to be

[TABLE]

The stochastic model (3.10) applied to 2D incompressible fluid dynamics with Hamiltonian $h(\varpi)$ and stochastic Hamiltonians $h_{i}(\varpi,\psi_{i})=h_{i}(\psi_{i})$ is given by

[TABLE]

where we have used the formula $\operatorname{ad}^{*}_{\psi}\varpi=\{\varpi,\psi\}$ for the coadjoint operator for 2D incompressible fluids. For example, with the appropriate choice of the Hamiltonian, we can write the stochastic model (3.15) for the following cases:

• (a)

2D perfect fluid: $\psi=\frac{\delta h}{\delta\varpi}=-\Delta^{-1}\varpi$;

• (b)

2D rotating perfect fluid: $\psi=\frac{\delta h}{\delta\varpi}=-\Delta^{-1}(\varpi-f)$;

• (c)

2D rotating quasigeostrophy (QG): $\psi=\frac{\delta h}{\delta\varpi}=-(\Delta-\mathcal{F})^{-1}(\varpi-f)$;

where $\mathcal{F}$ and $f$ denote, respectively, the square of the inverse Rossby radius and the rotation frequency. For the stochastic Hamiltonians $h_{i}$, $i=1,...,N$, one may choose

[TABLE]

in which case $\frac{\delta h_{i}}{\delta\psi_{i}}=-\Delta\psi_{i}$.

Similarly to the discussion in Remark 3.2, the Itô form of equations (3.15) takes the same expression, since we have chosen stochastic Hamiltonians for which $\delta h_{i}/\delta\varpi=0$. In this case, one may obtain the Itô forms by simply replacing the Stratonovich noise $\circ\,dW_{i}(t)$ by the Ito noise $dW_{i}(t)$ without modifying the drift terms.

Remark 3.4** (Conserved correlation enstrophies).**

The stochastic model (3.15) does not preserve the well-known vorticity enstrophies $\Lambda(\omega)=\int_{\mathcal{D}}\Phi(\omega)d^{2}x\,,$ which are preserved in the deterministic case. However, the stochastic equation for $\psi_{i}$ in (3.15) does preserve the following correlation enstrophy functionals

[TABLE]

when the stochastic model (3.10) is applied to 2D incompressible fluid dynamics. In general, for the case $h=h(m)$ and $h_{i}=h_{i}(\zeta_{i})$, by equation (3.10), the corresponding functionals $\Lambda_{i}(\zeta_{i})$ verify

[TABLE]

which vanishes for 2D incompressible fluids after integration by parts and imposition of homogeneous boundary conditions.

Example 3.5** (3D incompressible Euler).**

We consider the stochastic Hamitonians

[TABLE]

where $F_{i}$ are smooth functions. Upon using the formula from (3.2) for incompressible flows,

[TABLE]

where $\mathbb{P}$ is the Hodge projector onto divergence free vector fields, the stochastic model (3.10) reads

[TABLE]

The stochastic terms can be written equivalently with the help of the stress tensors

[TABLE]

and pressures $p_{i}$, $i=1,...,N$. With these definitions, equation (3.18) becomes

[TABLE]

where the divergence is defined as $(\operatorname{div}\sigma_{i})_{b}=\partial_{a}({\sigma_{i}})^{a}_{b}$ for all $i$, and the individual $p_{i}$ are each found by solving a Poisson equation, with boundary conditions given by $\hat{n}_{a}(\operatorname{div}{\sigma_{i}})^{a}=0$. Recall that the vector fields $\zeta_{i}(t,x)$ are obtained from the given eigenvectors field $\xi_{i}(X)$ by the push-forward operation (3.1). As mentioned earlier in Remark 3.2, the Itô forms of the equations have the same expression.

Remark 3.6**.**

The three dimensional system (3.18) is reminiscent of incompressible magnetohydrodynamics (MHD), except it has a stochastic “$J\times B\,$” force depending on the sum over all of the $\zeta_{i}$. In following this analogy with MHD, we may introduce vector potentials $\alpha_{i}$ by writing $\zeta_{i}=:{\rm curl\,}\alpha_{i}$. Having done so, one notices that evolution under the stochastic system (3.18) preserves the integrals

[TABLE]

Proof.

By the second equation in the equation set (3.18), we have

[TABLE]

∎

The $\Lambda_{i}$ integrals are topological quantities known as correlation helicities which measure the number of linkages of the lines of each vector field $\zeta_{i}$ with itself. Conservation of the correlation helicity $\Lambda_{i}$ means that evolution by the stochastic system (3.18) cannot unlink the linkages of each divergence free vector field $\zeta_{i}$ with itself. This conclusion is the analogue of conservation of magnetic helicity in MHD.

Inclusion of additional advected tensor fields.

More generally, suppose that the fluid model involves a tensor field $q(t,x)$ advected by the fluid flow as in (2.2). The evolution of this advected field is this given by $q(t)=(g_{t})_{\ast}q_{0}$, where $q_{0}(X)$ is the initial value and $g_{t}\in G$ is the fluid flow. In this case, the variational principle is written in reduced Eulerian form as

[TABLE]

for variations $\delta u$, $\delta m$, and $\delta g$. The stationarity conditions yield the same first two equations of (3.7), whereas the third one becomes

[TABLE]

From its definition, the quantity $q(t)=(g_{t})_{\ast}q_{0}$ verifies $dq+\pounds_{dgg^{-1}}q=0$. For the case that $\delta h_{i}/\delta m=0$, this quantity is governed by the ordinary advection equation,

[TABLE]

Remark 3.7** (Itô form).**

For the case that $\delta h_{i}/\delta m=0$, passing to the Ito formulation does not introduce any change in the drift terms.

Example 3.8** (Rotating shallow water).**

Equations (3.22) and (3.23) for the inclusion of such advected quantities into Model 2 may be illustrated with the example of the rotating shallow water equation. In this case $\mathcal{D}$ is a two dimensional domain and $G$ is the group of diffeomorphisms of $\mathcal{D}$. Let us denote by $\eta(t,x)$ the water depth, by $B(x)$ the bottom topography, by $R(x)$ the Coriolis vector field, and $x=(x_{1},x_{2})\in\mathcal{D}$. The Lagrangian of the rotating shallow water system is

[TABLE]

where $g$ is the gravity acceleration. Taking the stochastic Hamiltonians $h_{i}(\zeta_{i})$ in (3.16), the stochastic variational principle in (3.21) produces the equations

[TABLE]

where the stochastic stress $\sigma_{i}$ is defined as above in (3.19). The effect of $\sigma_{i}$ can be seen by writing the Kelvin circulation theorem obtained by integrating the first equation in (3.25) around a loop $c(u)$ moving with the drift velocity $u(t,x)$, to find

[TABLE]

Thus, the total stochastic stress generates circulation of the total velocity $(u+R)$ around any loop moving with the relative fluid velocity $u$ in the rotating frame.

Example 3.9** (Rotating compressible barotropic fluid).**

The above developments easily extend to other fluid models such as compressible barotropic fluid flow in a rotating frame, whose Lagrangian is

[TABLE]

where $z=x_{3}$ is the vertical coordinate, $\rho$ is the mass density and $e$ is the specific internal energy. With the choice (3.16), one gets from (3.22), the stochastic balance of momentum

[TABLE]

with advection equations

[TABLE]

4 Model 3: Eigenvectors depending on advected quantities

The stochastic model detailed in this section applies to fluids with advected quantities. As for the model described in §3, the model in this section also introduces a modification of the symmetries of the stochastic Hamiltonians $H_{i}$ in the SHP principle (2.14). Namely, in the presence of an advected tensor field $q(t,x)$, evolving as by the push-forward $q(t)=(g_{t})_{\ast}q_{0}$, we assume that $H_{i}$ depends on $q_{0}$. This dependence breaks the $G$-symmetry of $H_{i}$. In particular, we shall take $H_{i}(\_\,,\_\,;q_{0}):T^{*}G\rightarrow\mathbb{R}$ to be a $G_{q_{0}}$-invariant function, i.e.,

[TABLE]

where $G_{q_{0}}=\{g\in G\mid g^{\ast}q_{0}=q_{0}\}\subset G$ is the isotropy subgroup of the tensor field $q_{0}$ under the pull-back action. This is similar to the symmetry assumed on the Lagrangian $L(\_\,,\_\,;q_{0}):TG\rightarrow\mathbb{R}$ for such fluids, studied in Holm et al. [1998]. With this assumption, the SHP principle (2.14) becomes

[TABLE]

If we also assume that $H_{i}$ is linear in the material fluid momentum $\pi$, then it necessarily takes the form

[TABLE]

for a function $\xi_{i}:\mathcal{O}_{q_{0}}\rightarrow\mathfrak{g}$ defined on the $G$-orbit of $q_{0}$, $\mathcal{O}_{q_{0}}=\{g^{\ast}q_{0}\mid g\in G\}$. In the simplest case, the vector fields $\xi_{i}$ only depend on the pointwise values of the tensor fields, denoted as $\xi_{i}(q(x),x)$, but more general dependencies are possible, in which for example, the vector fields $\xi_{i}$ would depend on the spatial gradients of the tensor fields.

The reduced Eulerian version of this SHP principle reads

[TABLE]

The stochastic Hamilton-Pontryagin principle (4.1) and its reduced form (4.3) extend to the stochastic case the Hamilton-Pontryagin principles with advection developed in Gay-Balmaz and Yoshimura [2015].

By comparing (4.3) with (2.13), we observe that after the introduction of symmetry breaking in this SHP principle (4.1) in allowing the stochastic Hamiltonians in (4.2) to depend functionally on the initial advected quantities $q_{0}$ the reduction process has allowed the vector fields $\xi_{i}$, $i=1,...,N$ in the RSHP (4.3) to depend on the advected tensor field $q(t)=(g_{t})_{\ast}q_{0}$, as we have sought.

Applying Hamilton’s principle to the reduced action integral in (4.3) now results in the stochastic motion equation for the model that we will discuss in this section,

[TABLE]

where $dx_{t}:=u\,dt+\sum_{i=1}^{N}\xi_{i}(q)\circ dW_{i}(t)$ and the advected tensor field $q(t)=(g_{t})_{\ast}q_{0}$ verifies

[TABLE]

The expression of the coadjoint operators $\operatorname{ad}^{*}$ in the compressible and incompressible cases are recalled in equation (3.2).

In (4.4), the dot product $(\frac{\delta\ell}{\delta u}\cdot\frac{\partial\xi_{i}}{\partial q})$ is the composition of two linear maps. The contraction in the dot product is taken on the vector indices of the variational derivative $\frac{\delta\ell}{\delta u}$ with the derivative of the vector field $\frac{\partial\xi_{i}}{\partial q}$ of $\xi_{i}(q)$ with respect to $q$, so that, using (2.6),

[TABLE]

for $v\in\mathfrak{g}$. For a density, $q=\rho$, we have $\langle\phi\diamond\rho\,,\,v\rangle=\langle\phi\,,\,{\rm div}(\rho v)\rangle$, so $\langle\phi\diamond\rho\,,\,v\rangle=-\langle\rho\nabla\phi\,,\,v\rangle$ for a scalar function $\phi$. Hence, in this example, we have

[TABLE]

Hamiltonian structure.

From right-invariance of the Hamiltonian (4.2) under the isotropy subgroup $G_{q_{0}}$, we obtain the reduced stochastic Hamiltonians

[TABLE]

see the Appendix for more discussion. Upon denoting by $h:\mathfrak{g}^{\ast}\times\mathcal{O}_{q_{0}}\rightarrow\mathbb{R}$, the Hamiltonian associated to $\ell$, the above system can be written equivalently as

[TABLE]

One may check that system (4.8) is Hamiltonian with respect to the Poisson bracket

[TABLE]

on $\mathfrak{g}^{\ast}\times\mathcal{O}_{q_{0}}\ni(m,q)$. This Poisson bracket is inherited by Poisson reduction of the canonical Poisson bracket on $T^{*}G$, by the isotropy subgroup $G_{q_{0}}\subset G$ of $q_{0}$, namely, the map

[TABLE]

is Poisson with respect to these brackets, as discussed in the Appendix. Thus, the system (4.8) admits the Stratonovich-Poisson formulation

[TABLE]

where both $h$ and $h_{i}$ depend on $m$ and $q$.

Remark 4.1** (Itô forms).**

The Itô forms of the equations in (4.8) are rather involved. In particular, the Itô correction in the drift term of the momentum equation is

[TABLE]

Likewise, the Itô correction in the drift term of the advection equation is

[TABLE]

The drift terms $M_{\rm Ito}dt$ and $A_{\rm Ito}dt$ above follow from the double bracket terms in equation (3.14).

Example 4.2** (Rotating shallow water).**

Upon choosing the Lagrangian (3.24) for the rotating shallow water equation, we obtain from (4.4)

[TABLE]

where $dx_{t}:=udt+\sum_{i=1}^{N}\xi_{i}(\eta)\circ dW_{i}(t)$. This equation may be written more explicitly as

[TABLE]

The advection equation for the surface elevation in this model is

[TABLE]

Taking the curl of the momentum equation yields the corresponding RSW vorticity equation

[TABLE]

This, together with the advection equation (4.14) yields the potential vorticity (PV) equation,

[TABLE]

which expresses conservation of potential vorticity along the stochastic fluid particle path $dx_{t}$. **

Remark 4.3**.**

We shall now consider a specific expression for $\xi_{i}(h)$, that yields a simplified expression for the stochastic PDE (4.13). First, we note that by assembling the stochastic terms in (4.13), we get, in slightly abbreviated notation,

[TABLE]

Upon making the choice $\xi_{i}(\eta):=-\eta^{-1}X_{i}$, for a set of fixed vector fields $X_{i}$, $i=1,...,N$, we have $\frac{\partial\xi_{i}}{\partial\eta}=\eta^{-2}X_{i}$ and $\frac{\partial\xi_{i}}{\partial\eta}^{\ast}\cdot m=\eta^{-2}(m\cdot X_{i})$, so that

[TABLE]

in the expression above. In this case, the stochastic RSW equations (4.13) take the simpler form

[TABLE]

Having chosen the simpler form above with $\xi_{i}(\eta):=-\eta^{-1}X_{i}$, we may write the stochastic term of (4.15) in terms of the potential vorticity $Q$ as

[TABLE]

Example 4.4** (Rotating compressible barotropic fluid).**

Upon choosing the tensor field $q(t,x)$ to be the mass density $\rho(t,x)$, the stochastic model (4.4) yields

[TABLE]

where $dx_{t}:=udt+\sum_{i=1}^{N}\xi_{i}(\rho)\circ dW_{i}(t)$. Upon using the advection equation,

[TABLE]

the previous equation can be written equivalently as

[TABLE]

Upon taking the Lagrangian (3.26) for the compressible barotropic fluid, we obtain the stochastic motion equation,

[TABLE]

in which the thermodynamic pressure is defined by $p=\rho^{2}\frac{\partial e}{\partial\rho}$.

Taking the curl of the momentum equation and using the advection equation $d\rho+\operatorname{div}(\rho\,dx_{t})=0$ now implies PV conservation as

[TABLE]

Similarly to the case of the RSW equations, with the choice $\xi_{i}(\rho)=-\,X_{i}/\rho$ for a set of fixed vector fields $X_{i}$, $i=1,...,N$, the previous motion equation simplifies to

[TABLE]

Remark 4.5**.**

The last two examples preserve the integrals

[TABLE]

for smooth functions $\Phi$, as may be seen either from the Hamiltonian formulation (4.10) with Poisson bracket (4.9) with $q=\rho$, or by direct verification. These conserved quantities are Casimirs for the Poisson bracket. They are the same as the conserved potential vorticity integrals in the deterministic case, because the Poisson bracket persists in passing to the stochastic case.

5 Conclusions

This paper has used standard methods from symmetry breaking in geometric continuum mechanics, as discussed for example in Holm et al. [1998], Gay-Bamaz and Tronci [2010], to set out three different approaches for incorporating stochastic transport into ideal fluid dynamics. As mentioned in the Introduction, our approach preserves the relabelling symmetry of fluid dynamics, up to isotropy of the advected quantities. Therefore, we may summarise the results of these approaches, simply by comparing their advection laws and the symmetry breaking terms of their respective Kelvin-Noether circulation theorems.

Model 1 – Time independent spatial correlation eigenfunctions.

In Model 1, reviewed in section 2, the time-independent spatial statistical correlation eigenvectors $\xi_{i}(x)$ are obtained as eigenvectors of an appropriate correlation function, which is assumed to be time-independent. From equations (2.19) and (2.20), possibly extended to include advected quantities $q$, the corresponding Kelvin circulation theorem is given by

[TABLE]

where $\rho$ is the mass density obeying the continuity equation, $d\rho+\operatorname{div}(\rho dx_{t})=0$, and the advected tensor field $q=(g_{t})_{*}q_{0}$ satisfies

[TABLE]

where $\pounds_{v}q$ denotes Lie derivative of advected quantity $a$ by the stochastic vector field

[TABLE]

This is the model introduced in Holm [2015].

Model 2 – Time-dependent advected statistical correlation eigenvectors.

From equation (3.22) for Model 2, discussed in section 3 we have the circulation dynamics

[TABLE]

and with ${\delta h_{i}}/{\delta m}=0$ the advected tensor fields $q(t)=(g_{t})_{\ast}q_{0}$ and $\zeta_{i}(t)={\rm Ad}_{g_{t}}\xi_{i}$ verify

[TABLE]

In particular, the mass density verifies the ordinary continuity equation $d\rho+\operatorname{div}(\rho u)dt=0$.

This means that in Model 2, treated in section 3, the choice of stochastic Hamiltonians satisfying $\delta h_{i}/\delta m=0$ implies that the transport velocity of the circulation loop in (5.3) reduces to the drift velocity, $u(x,t)dt$, alone, rather than the entire stochastic velocity $dx_{t}$, as in Model 1. Moreover, for $\delta h_{i}/\delta m=0$, the Stratonovich and Itô representations of Model 2 take the same form.

Remark 5.1**.**

In deriving the Kelvin theorem in equation (5.3) from (3.22), we have used the relation

[TABLE]

which arises from equivalence of the definitions of ${\rm ad}^{\ast}$ and the diamond operator $(\diamond)$ when the advected quantity is a vector field. This notation makes it clear that the drift term and the stochastic term on the right-hand side of equation (5.3) both arise from the same source; namely, symmetry breaking of the Lie group of diffeomorphisms to the isotropy subgroup of the initial condition for an advected quantity.**

Model 3 – Time-dependent spatial correlation eigenvectors depending on advected quantities.

From equation (4.4) for Model 3, discussed in section 4, we have the circulation dynamics

[TABLE]

where $dx_{t}:=udt+\sum_{i=1}^{N}\xi_{i}(q)\circ dW_{i}(t)$ and the advected tensor fields $q=(g_{t})_{\ast}q_{0}$ verify

[TABLE]

where $\pounds_{dx_{t}}q$ denotes Lie derivative of advected quantity $q$ by the stochastic vector field

[TABLE]

In particular the mass density satosfies $d\rho+\operatorname{div}(\rho dx_{t})=0$.

In Model 3, treated in section 4, the velocity of the circulation loop is $dx_{t}$. As in Model 1 this is the full stochastic velocity, except now its spatial correlation eigenvectors $\xi_{i}(q)$ respond to the motion of the advected quantities which also contribute to the thermodynamic and potential energy properties of the fluid.

Remark 5.2**.**

The Itô forms of these three models have been computed in the body of the paper, although the additional Itô drift terms may take a more complicated form than in the Stratonovich case. An example is the additional Itô drift term $M_{\rm Ito}dt$ in (4.11) for Model 3. Nonetheless, the equivalent Itô forms of their circulation theorems may be computed by the standard methods from their Stratonovich forms; namely, by transforming both the integrand and the velocity of the loop into Itô form, using the corresponding motion equations.

Remark 5.3**.**

Besides the initial and boundary conditions, the key to applying the methods of this paper efficiently begins with the choice of the spatial correlation eigenvectors. Once a data assimilation step for determining the eigenvectors is completed, such as the one discussed in Remark 2.1 based on Cotter et al. [2017], then the stochastic Hamilton’s principle provides the dynamics of the drift velocity. Thus, the present approach may be regarded as a data driven approach to dynamically and self-consistently separating the Lagrangian paths into their drift dynamics and stochastic parts. This question needs to be confronted in any application to any data set. Model 1 chooses to model the stochastic parts by using optimal eigenvectors of an equilibrium spatially-constant time-independent correlation function. Models 2 and 3 provide two different approaches for making these correlation eigenvectors time-dependent. For a detailed discussion of methods for modelling the separation between drift and stochasticity in comparing numerical simulations of fluid flows at fine and coarse resolutions, see Cotter et al. [2017].

In applications, these three models of the stochastic fluid velocity vector field could be used either separately, or in combinations, by taking combinations of the correlation eigenvectors to evolve in any way that might be needed for tuning the stochastic transport, e.g., in applying them for data assimilation.

In conclusion, we mention that the three models of stochastic fluid dynamics treated here may also be transferred into data assimilation methods in biomedical image analysis for computational anatomy. In particular, the approach of Arnaudon et al. [2017b] for stochastic image analysis based on Model 1 could be extended to create new approaches to computational anatomy based on the ideas underlying Model 2 and Model 3.

Appendix

The Appendix provides additional background for some of the standard notations from geometric mechanics used in the paper and gives some details on the structure of the reduced Poisson brackets $\{\,\cdot\,,\cdot\,\}_{\rm red}$ arising in Model 2 and Model 3. For further geometric mechanics background, see Marsden and Ratiu [1994]; Holm [2011].

Lie group notation.

Given a Lie group $G$, and two elements $g,h\in G$, we denote by $gh$ the Lie group multiplication. The spaces $T_{g}G$ and $T^{*}_{g}G$ refer to the tangent and cotangent spaces of $G$ at $g$. For $v\in T_{g}G$ and $\pi\in T^{*}_{g}G$, we denote by $vh\in T_{gh}G$ and $\pi h\in T^{*}_{gh}G$, the right translation of $v$ and $\pi$ by $h\in G$. These right translations are defined as follows:

[TABLE]

where $g_{\varepsilon}\in G$ is a path with $g_{\varepsilon=0}=g$ and $\left.\frac{d}{d\varepsilon}\right|_{\varepsilon=0}g_{\varepsilon}=v$, and $\left\langle\,\cdot\,,\cdot\,\right\rangle$ is the duality pairing between the tangent and cotangent spaces of $G$. We have $\mathfrak{g}=T_{e}G$ and $\mathfrak{g}^{\ast}=T_{e}^{\ast}G$, where $e$ is the neutral element in $G$. The adjoint action of $g\in G$ on $u\in\mathfrak{g}$ is defined by

[TABLE]

where $h_{\varepsilon}\in G$ is a path with $h_{\varepsilon=0}=e$ and $\left.\frac{d}{d\varepsilon}\right|_{\varepsilon=0}h_{\varepsilon}=u$.

When $G$ is the group of diffeomorphism of the fluid domain $\mathcal{D}$, the multiplication $gh$ corresponds to the composition of diffeomorphisms. Elements $v\in T_{g}G$ are vector field on $\mathcal{D}$ covering the diffeomorphism $g$ and $vh$ is the composition on the right by the diffeomorphism $h$. In particular, elements in the Lie algebra $\mathfrak{g}$ are vector fields on $\mathcal{D}$. In the paper, we have extensively used the fact that on groups of diffeomorphisms, $\operatorname{Ad}_{g}u=g_{\ast}u$, the push-forward of vector fields.

Reduced Poisson brackets.

The Poisson brackets (3.11) and (4.9) of models 2 and 3 both have the same structure, which we now explain. Both brackets are inherited from the canonical Poisson bracket $\{\,\cdot\,,\cdot\,\}_{\rm can}$ on $T^{*}G$, by a reduction by symmetry associated to a subgroup of $G$.

Consider a right action of $G$ on a manifold $P$, denoted as follows $(n,g)\in P\times G\mapsto\phi_{g}(n)\in P$. Given $n_{0}\in P$, we denote by $G_{n_{0}}=\{g\in G\mid\phi_{g}(n_{0})=n_{0}\}\subset G$ the isotropy group of $n_{0}$ and by $\mathcal{O}_{n_{0}}=\{\phi_{g}(n)\mid g\in G\}\subset P$ the orbit of $n_{0}$. The Lie algebra action of $u\in\mathfrak{g}$ on $n\in N$ is written as

[TABLE]

were $g_{\varepsilon}\in G$ is a path with $g_{\varepsilon=0}=e$ and $\left.\frac{d}{d\varepsilon}\right|_{\varepsilon=0}g_{\varepsilon}=u$.

The quotient space of $T^{*}G$ by $G_{n_{0}}$ relative to the action given by $(g,\pi)\in T^{*}G\mapsto(gh,\pi h)\in T^{*}G$, $h\in G_{n_{0}}$, is denoted as $(T^{*}G)/G_{n_{0}}$. One observes that we have the identification $(T^{*}G)/G_{n_{0}}\simeq\mathfrak{g}^{\ast}\times\mathcal{O}_{n_{0}}$ given by

[TABLE]

A $G_{n_{0}}$-invariant Hamiltonian $H:T^{*}G\rightarrow\mathbb{R}$ thus induces the reduced Hamiltonian $h:\mathfrak{g}^{\ast}\times\mathcal{O}_{n_{0}}\rightarrow\mathbb{R}$ defined by

[TABLE]

This formula has been used several times in the paper, for instance in (3.5) and (4.7) for the stochastic Hamiltonians. The reduced Poisson bracket on $\mathfrak{g}^{\ast}\times\mathcal{O}_{n_{0}}$ induced by the canonical Poisson bracket $\{\,\cdot\,,\cdot\,\}_{\rm can}$ on $T^{*}G$ is given by

[TABLE]

for functions $f,g:\mathfrak{g}^{\ast}\times\mathcal{O}_{n_{0}}\rightarrow\mathbb{R}$, see, e.g, Gay-Bamaz and Tronci [2010] (Theorem 2.7) for a proof of this fact. The last two terms involve the pairing $\left\langle\,\cdot\,,\cdot\,\right\rangle$ between $T^{\ast}\mathcal{O}_{n_{0}}$ and $T\mathcal{O}_{n_{0}}$. We now explain how this Poisson reduction $T^{*}G\rightarrow\mathfrak{g}^{\ast}\times\mathcal{O}_{n_{0}}$ applies to Models 2 and 3.

For Model 2, we have $P=\mathfrak{g}\times...\times\mathfrak{g}\ni(\zeta_{1},...,\zeta_{N})$, $n_{0}=(\xi_{1},...,\xi_{N})=\xi$ and the action $\phi_{g}$ given by the push-forward of vector fields $\operatorname{Ad}_{g}\xi=g_{\ast}\xi$ on each term. With these choices, the reduced Poisson bracket (5.9) recovers the expression (3.11). For the system (3.22), we take $P=V\times\mathfrak{g}\times...\times\mathfrak{g}\ni(q,\zeta_{1},...,\zeta_{N})$, $n_{0}=(q_{0},\xi_{1},...,\xi_{N})$ and the action $\phi_{g}$ given by the push-forward action on each terms.

For Model 3, we have $P=V\ni q$, $n_{0}=q_{0}$ and the action $\phi_{g}$ given by the push-forward of tensor fields, $\phi_{g}(q)=g^{\ast}q$. With these choices, the reduced Poisson bracket (5.9) recovers the expression (4.9).

Acknowledgements.

We are very grateful to C. J. Cotter, D. Crisan, E. Mémin, S. Olhede, W. Pan, V. Resseguier, I. Shevshenko, and A. Sykulski for valuable discussions and presentations during the course of this work. During this work, FGB was partially supported by the ANR project GEOMFLUID 14-CE23-0002-01 and DDH was partially supported by the European Research Council Advanced Grant 267382 FCCA and EPSRC Standard Grant EP/N023781/1.