QuasiGeostrophic Equations

Quasi-Geostrophic Equations

CNRS

MEOM

In this section, we look at how we can solve the QG equations using elements from this package.

Generalized QG Equations¶

Note: This formational is largely based on [Amraoui et al., 2023]. See this paper for more information including a derivation about the energetics.

We can define two fields as via the potential vorticity (PV) and stream function (SF).

\begin{aligned} \text{Stream Function}: && \psi & =\boldsymbol{\psi}(\vec{\mathbf{x}},t) && && \vec{\mathbf{x}}\in\Omega\sub\mathbb{R}^{D_s} , && t\in\mathcal{T}\sub\mathbb{R}^{+} \\ \text{Potential Vorticity}: && q& =\boldsymbol{q}(\vec{\mathbf{x}},t) && && \vec{\mathbf{x}}\in\Omega\sub\mathbb{R}^{D_s}, && t\in\mathcal{T}\sub\mathbb{R}^{+} \\ \end{aligned}

(1)

Let's assume we have a stack of said fields, i.e. $q_k,\psi_k$ . So we have a state defined as:

\begin{aligned} \text{Stream Function}: && \vec{\boldsymbol{\psi}} = [\psi_1, \psi_2, \ldots, \psi_K]^\top \\ \text{Potential Vorticity}: && \vec{\boldsymbol{q}} = [q_1, q_2, \ldots, q_K]^\top \\ \end{aligned}

(2)

We can write the QG PDE to describe the spatiotemporal relationship between the PV and the SF which is defined as:

\begin{aligned} \partial_t q_k + \vec{\boldsymbol{u}}_k\cdot\boldsymbol{\nabla}q_k = F_k + D_k \end{aligned}

(3)

where at each layer $k$ , we have $\vec{\boldsymbol{u}}_k$ is the velocity vector, $q_k$ is the potential vorticity, $F_k$ are the forcing term(s), and $D_k$ are dissipation terms. The forcing term can be made up of any external forces we deem necessary for the PDE. For example, we could have wind stress that affects the top layer or bottom friction that affects the bottom layer. The dissipation term represents all of the diffusion terms. For example, we could have some lateral friction or hyper-viscosity terms.

The advection term in equation (3) includes the velocity vector, $\vec{\boldsymbol{u}}_k$ , which is defined in terms of the stream function, $\psi$ . This is given by:

\vec{\boldsymbol{u}} = [u, v]^\top = \left[ -\partial_y \psi, \partial_x \psi\right]^\top

(4)

The PV is the culmination of the potential energy contained within the dynamical forces, the thermal/stretching forces and the planetary forces.

\text{Potential Vorticity Forces} = \text{Dynamical} + \text{Thermal} + \text{Planetary} + \ldots

(5)

Concretely, we can define this the PV and its relationship to the SF and other forces as:

\begin{aligned} q_k = \boldsymbol{\Delta}\psi_k + (\mathbf{M}\psi)_k + f_k \end{aligned}

(6)

where at each layer, $k$ , we have the dynamical (relative) vorticity, $\boldsymbol{\Delta}\psi_k$ , the thermal/stretching vorticity, $(\mathbf{M}\psi)_k$ , and the planetary vorticity, $f_k$ .

Idealized QG Model¶

We can use the above formulation to describe an idealized QG model. This idealized model will be fairly abstract with parameters that represent real things but do not actually correspond to anything tangible. However, these types of flows are very well-studied and are very useful for having a controlled setting which removes all uncertainties. We will use the notation used in [Frezat et al. (2022)Srinivasan et al. (2023)] with some explanations taken from [Gupta & Brandstetter, 2022].

The domain is simplified to be periodic and defined on a $2\pi$ domain, i.e. $Lx = L_y = 2\pi$ . We will also assume a periodic domain. This domain is topologically equivalent to a two-dimensional torus.

We will define a final form for the idealized QG model of the vorticity-stream functions on this idealized, periodic domain as

\partial_t \omega + \det\boldsymbol{J}(\psi,\omega) = - \mu\omega + \nu\boldsymbol{\nabla}^2\omega - \beta\partial_x\psi + F

(7)

where $\omega$ is the vorticity and $\psi$ is the stream function. The $\beta\psi_x =\beta v$ -term is geophysical system approximation that captures the effect of the differential rotation. This force is experienced within the Earth system in a tangent plane approximation, i.e. $\beta = \partial_y f$ . This $\beta$ -term is important and it allows the flow to manifest different turbulent regimes. For example, Rossby waves are common within the planetary systems and can appear when $\beta=0$ . The determinant Jacobian term encapsulates the advection term seen in (3). It is defined as:

\begin{aligned} \text{Determinant Jacobian}: && && \det \boldsymbol{J}(\psi,\omega) &= \partial_x\psi\partial_y \omega - \partial_y\psi\partial_x \omega \end{aligned}

(8)

We see from that it is directly related to the advection term seen in (3) but written in a different way.

Determinant Jacobian

Like most things in physics, there are often many ways to express the same expression. Ultimately, they are all advection expressions. See the example on wikipedia for more details. The detemrinant Jacobian term can be written in many ways. Let's rewrite it as:

\text{Advection} = \det \boldsymbol{J}(\psi,\omega)

(9)

We can expand the full expression which gives us

\text{Advection} = \partial_x\psi\partial_y \omega - \partial_y\psi\partial_x \omega

(10)

We can plug int the velocity components of the stream function definition (4) into the above equation

\text{Advection} = v\partial_y \omega - u\partial_x \omega

(11)

The partial derivative operator is commutable so we can take out the operator of both terms

\text{Advection} = \partial_y (v\omega) - \partial_x (u\omega)

(12)

Alternatively, we can write the velocity and partial derivative operators in the vector format

\text{Advection} = \vec{\boldsymbol{u}} \cdot \boldsymbol{\nabla} \omega

(13)

and we see that we arrive at formulation in (7). I personally prefer this way of writing it as it is more general. Furthermore, it exposes the many ways we can express this term like the determinant Jacobian or simple partial derivatives.

Note: It is important to take note of the many ways to express this as it can be useful for numerical schemes. For example, an upwind scheme might benefit from advection term where the velocity components are multiplied with the partial derivatives. Alternatively the determinant Jacobian on an Arakawa C-grid is a well known formulation for dealing with this.

The forcing term is typically chosen to be a constant wind forcing

\boldsymbol{F}_\omega(\vec{\mathbf{x}}) = k_f \left[ \cos (k_f x) + \cos (k_f y)\right]

(14)

Relation to Navier-Stokes Equations

This derivation and explanation was largely taken from [Srinivasan et al., 2023] which has one of the best high-level explanation of the derivation without too much mathematical detail.

Taking a velocity field, $\vec{\boldsymbol{u}} = [u,v]^\top$ , we can write the non-dimensional form of the Navier-Stokes (N-S) equations as

\partial_t \vec{\boldsymbol{u}} + \vec{\boldsymbol{u}} \cdot \boldsymbol{\nabla}\vec{\boldsymbol{u}} + f (k \times \vec{\boldsymbol{u}}) = -\mu \vec{\boldsymbol{u}} + \frac{1}{Re}\boldsymbol{\nabla}^2\vec{\boldsymbol{u}} + \boldsymbol{F}_{\vec{\boldsymbol{u}}}(\vec{\mathbf{x}})

(15)

The, $f$ , is the Coriolis parameter which is the local rotation rate of the Earth and/or other planetary atmospheric forcings. The example in showcases an example with beta-plane forcing. $\boldsymbol{k}$ is the unit-vector normal to the $(x,y)$ -plane. The $\mu$ is the linear drag coefficient which represents the bottom friction. The $Re$ is the Reynolds number measuring the strength of the non-linear advection term, relative to the viscous term. In otherwords, the relationship give by:

\vec{\boldsymbol{u}} \cdot \boldsymbol{\nabla}\vec{\boldsymbol{u}}\propto \frac{1}{Re}\boldsymbol{\nabla}^2\vec{\boldsymbol{u}}

(16)

This Reynolds number is indirectly proportional to the viscosity and proportional to the absolute velocity [Gupta & Brandstetter, 2022]

[Srinivasan et al., 2023] choose a Reynolds number of $Re=2,500$ . The forcing is given by

\begin{aligned} \text{Forcing}: && \boldsymbol{F}_{\vec{\boldsymbol{u}}}(\vec{\mathbf{x}}) &= [-\sin(k_f y), \sin(k_f x)]^\top \end{aligned}

(17)

which is a sinusoidal time-invariant forcing field that continuously drives the flow. In [Guan et al., 2023], the wavenumber $k_f=4$ , was chosen.

The velocity field in equation (15) is required to satisfy the mass conservation principal given by the continuity equation

\begin{aligned} \text{Continuity Equation}: && \boldsymbol{\nabla}\cdot\vec{\boldsymbol{u}} &= \partial_yu + \partial_xv = 0 \end{aligned}

(18)

One can satisfy this by defining a stream function, $\psi(\vec{\mathbf{x}})$ , which is a scalar field that is defined as

\begin{aligned} u = - \psi_y, && && v =\psi_x \end{aligned}

(19)

The stream function and the continuity equations can be expressed as an evolution equation of a single scalar field, the vorticity. This scalar field is defined as the two-dimensional curl of the velocity field

\omega = \boldsymbol{\nabla}\times \vec{\boldsymbol{u}}=\partial_x v - \partial_y u = \boldsymbol{\nabla}^2\psi

(20)

This equation captures the local rotation of a fluid parcel. The final result is the incompressible 2D Navier-Stokes equations in the scalar vorticity stream function form. In other words, the Quasi-Geostrophic equations.

\partial_t \omega + \det\boldsymbol{J}(\psi,\omega) - \beta v= - \mu\omega + \frac{1}{Re}\boldsymbol{\nabla}^2\omega + F

(21)

Note: There are some small differences between this equation and . The first is the coefficient in front of the diffusion term, $\boldsymbol{\nabla}^2\omega$ . Here, we have the Reynolds number, $1/Re$ instead of the viscosity term, $\nu$ , as shown in (7). In addition, we have the $\beta$ term. In this formulation, it is $\beta v$ whereas in (7) it is expressed as $\beta \partial_x$ . However, these are equivalent because the first component of the stream function velocities (19) is defined as $v=\partial_x\psi$ . So we can plug this into the equation above.

NOTE: I am not sure about the sign issue of the $\beta$ -term in (21). I think it is a mistake and that it should be positive which would match the equation in (7) along with various other formulations [Frezat et al., 2022]

Parameter Configurations¶

Below are some experimental parameters found in [Frezat et al., 2022] which showcase 3 different flow regimes based on the parameter scheme.

Table 1:Table with idealized configuration

Name	Symbol	Units	Decay Flow	Forced Flow	$\beta$ -Plane Flow
Resolution	$N_x\times N_y$		$2,048\times 2,048$	$2,048\times 2,048$	$2,048\times 2,048$
Domain	$L_x\times L_y$	km	$10e3 \times 10e3$	$10e3 \times 10e3$	$10e3 \times 10e3$
Time Step	$\Delta_t$	s	$120$	$120$	$120$
Linear Drag Coefficient	$\mu$	m $^{-1}$	$0$	$1.25e-8$	$1.25e-8$
Viscosity	$\nu$	m $^2$ s $^{-1}$	$67.0$	$22.0$	$22.0$
Beta-Term	$\beta$	m $^{-1}$ s $^{-1}$	$0.0$	$0.0$	$1.14e-11$
Reynolds Number	$Re$		$32e3$	$22e4$	$34e4$

Toy QG¶

In this section, we will briefly outline an idealized QG model with terms that correspond to real quantities in nature. This was taken from the lab that's available online (PDF | Course | Notebook | Code). It features a great step-by-step introduction to a complete toy problem.

\partial_t q + \vec{\boldsymbol{u}}\cdot q = \frac{1}{\rho H}\boldsymbol{\nabla}_H\times\vec{\boldsymbol{\tau}} - \kappa\boldsymbol{\nabla}_H\psi + a_4\boldsymbol{\Delta}_H^2 - \beta\partial_x\psi

(22)

where $q=\boldsymbol{\nabla}^2\psi$ is the PV. This PDE describes the transport of the PV with respect to the following

\partial_t q + \text{Advection} = \text{Wind Stress} + \text{Bottom Friction} + \text{Dissipation} + \text{Planetary Forces}

(23)

which correspond to all of the terms in equation (22).

Horizontal Operators. Apart from the partial derivative operators, we have some horizontal operators which means that they only operate on the horizontal axis, i.e. $(x,y)$ . $\boldsymbol{\nabla}_H$ is the horizontal gradient operator defined as $\boldsymbol{\nabla}_H = [\partial_x, \partial_y]^\top$ . $\boldsymbol{\Delta}_H$ is the horizontal Laplacian operator defined as $\boldsymbol{\Delta}_H = [\partial^2_x, \partial^2_y]^\top$ . $\boldsymbol{\Delta}_H^2$ is the hyper-Laplacian operator which is normally applied sequentially, i.e. $\boldsymbol{\Delta}_H\circ\boldsymbol{\Delta}_H$ .

Advection. We have the advection term defined as the dot product velocity vector (equation (4)) defined via the stream function. As mentioned above, there are many ways to write this advection term, for example the determinant Jacobian.

Wind Stress. We have defined a wind stress term which is the cross product of the horizontal gradient and the wind stress, $\vec{\boldsymbol{\tau}}$ . As shown above, this term is a 2D vector field for the x and y directions.

\vec{\boldsymbol{\tau}} := \vec{\boldsymbol{\tau}}(x,y) = \left[ \tau_x(x,y), \tau_y(x,y)\right]^\top

(24)

An "ideal" example os to assume the ocean is a "box" that extends from the equator to a latitude of about 60 $^\circ$ . The winds are typically constant and easterly near the equation and turn westerly at middle latitudes. In ideal setting, this is prescribed as a constant cosine forcing $\vec{\boldsymbol{\tau}}=\vec\tau_{\text{max}}(-\cos y,0)$ which is similar to the idealized QG model listed above. This product is then scaled according to the inverse of the density of water, $\rho$ , and the height of the layer, $H$ .

Bottom Friction. The bottom friction parameter is a scaling constant that dissipates the energy. The constant, $\kappa$ , is defined as inversely proportional to the product of the height, $H$ , and the square-root of the vertical Eddy viscosity constant and the Coriolis parameter. This constant is multiplied with the Laplacian of the stream function.

Dissipation. There is a dissipation term which features the hyper-viscosity term. The horizontal hyper-Laplacian of the stream function is multiplied by some constant $\nu_H$ which represents the lateral Eddy viscosity that may occur in the system.

Planetary Forcing The planetary forcing within this toy model is given by the $\beta$ -plane approximation for the planetary vorticity. This is given by the $\beta$ term which is defined as

\beta \approx \partial_y f

(25)

where $f$ is the planetary forcing function given by equation . The only term in this function that depends upon $y$ is the second term so we approximate this as $\beta$ .

Parameter Details

There are many constants that are displayed in equation (22) and they each mean something. Below are some tables which outline each of the constants and parameters seen in the equation. The subsequent tables showcase some derived quantities. The last table showcases the nondimensional versions which might be useful for the non-dimensional version of the PDE.

The first table showcases the constants along with the range of values that it can take for the toy model.

Table 2:Table with constants

Name	Symbol	Units	Value
Earth's Radius	$R$	m	$6400e3$
Earth's Angular Frequency	$\Omega$	s $^{-1}$	$7.25e-5$
Depth of Active Layer	$H$	m	$100 \rightarrow 4,000$
Length of Ocean Basin	$L$	m	$1000e3 \rightarrow 5000e3$
Density of Water	$\rho$	kg m $^{-3}$	$1e3$
Lateral Eddy Viscosity	$\nu_H$	m $^2$ s $^{-1}$	$0$ or $1e1 \rightarrow 1e4$
Vertical Eddy Viscosity	$\nu_Z$	m $^2$ s $^{-1}$	$1e-4 \rightarrow 1e1$
Maximum Wind Stress	$\tau_{\text{max}}$	kg m $^{-1}$ s $^{-2}$	$1e-2\rightarrow 1$
Mean/Mid Latitude	$\theta_0$		$0\rightarrow \pi/3$

The second table lists some of the derived constants that we can calculate given the above constants.

Table 3:Table of derived quantities.

Name	Equation	Units	Value
$\beta$ -Plane	$\beta = 2\Omega\cos\theta_0/R$		$1.1 \rightarrow 2.3e11$
Coriolis Parameter	$f_0 = 2\Omega\sin\theta_0$	s $^{-1}$	$0.0 \rightarrow 1.3e-4$
Velocity Scale	$U_0 = \tau_{\text{max}}/(\beta\rho H L)$		$1e-5 \rightarrow 1e-1$
Bottom Friction	$\kappa = \frac{1}{H\sqrt{(\nu_Z f_0)/2}}$		$0.0 \rightarrow 1e-5$

The last table lists the non-dimensional derived quantities which are useful for the non-dimensional version of the toy QG PDE (22).

Table 4:Table of non-dimensional quantities.

Name	Equation	Value
	$\epsilon/\text{Vorticity Ratio} = U_0/(\beta L^2)$
	$\tau_{\text{max}}/(\epsilon\beta^2\rho HL^3)$	$1e-12 \rightarrow 1e-14$
	$\kappa/(\beta L)$	$4e-4 \rightarrow 6e1$
	$\nu_H/(\beta L^3)$	$1e-7 \rightarrow 1e-4$

QG for Sea Surface Height¶

In this section, we will briefly outline how the QG equations can be used to represent the

There is a known relationship between sea surface height (SSH) and the stream function. This is given by

\eta = \frac{f_0}{g}\psi

(26)

where $f_0$ and $g$ is the Coriolis and gravity constant respectively.

This relationship has been exploited in various works in the literature, especially within the data assimilation literature. For example [Guillou et al. (2021)Guillou et al. (2023)Amraoui et al. (2023)] used a 1 layer QG model to assimilate NADIR alongtrack SSH observations. A related example was in [Guillou et al., 2021] where they used a simple 1 layer QG model and a 1 layer linear SW model to model the balanced motions and the internal tides respectively. This was also used to assimilate NADIR and SWOT alongtrack SSH observations.

To write out the PDE, we will use the same formulation listed in equation (3). However, we will change this to reflect the relationship between SSH, the SF and now the PV.

q = (\boldsymbol{\nabla}^2 - \frac{f_0^2}{c^2})\psi + f

(27)

where $f$ is the planetary forcing given by equation , $f_0$ is the Coriolis parameter and $c$ is the phase speed. This formulation is very similar to the PV listed in (6) except for the coefficient listed in front of the $\psi$ component. This can be summarized as the barotrophic deformation wavenumber

k = \frac{f_0}{c}=\frac{1}{L_R}

(28)

where $L_R$ is the Rossby radius of deformation. In an idealized setting, this is equal to $L_R = \sqrt{gH}/f_0$ It is a helpful parameter that can be used for the QG model of choice. For example, in [Guillou et al., 2021], they used a Rossby radius of $30$ km for a $10^\circ\times 10^\circ$ domain.

Example Experimental Setup¶

In the [Amraoui et al., 2023] paper, the authors did a free-run of the QG model for an assimilation task of 21 days. Their target area was over the North Atlantic (NA) domain of $2,052$ km $\times$ $3,099$ km with a grid cell size, $dx$ , of $\sim18$ km. This resulted in a domain size of $N_x\times N_y=113\times 170$ . They used a leapfrog time integration method with a time step, $\Delta t$ , was $600$ s. They needed to do $\sim$ 10,000 spinup steps to get a good running simulation. A slight change to the above formulation is that they used a forcing term as the constant wind forcing term outlined above.

Stacked QG Model¶

In this section, we will look at the formulation from the stacked QG model. The majority of this section comes from papers by [Thiry et al. (2022)Thiry et al. (2023)] with some inspiration from the Q-GCM numerical model.

In this case, we are similar to the generalized QG model outlined in equation (3). We will change the notation slightly. Let's have a stacked set of PV and SF vectors denoted as $\vec{\boldsymbol{q}}$ and $\vec{\boldsymbol{\psi}}$ respectively. We consider the stream function and the potential vorticity to be $N_Z$ stacked isopycnal layers. This is already outlined in equation (2) Now, we will write the stacked model in vector format as

\partial_t \vec{\boldsymbol{q}} + \vec{\boldsymbol{u}}\cdot\boldsymbol{\nabla} \vec{\boldsymbol{q}} = \mathbf{BF} + \mathbf{D}

(29)

where $F$ is a vector of forcing terms to be applied at each layer, $k$ , $\mathbf{B}$ is a matrix of interlayer height scaling factors and $\mathbf{D}$ is a vector of dissipation terms to be applied at each layer, $k$ . This is analogous to equation (3) but in vector format instead of a component-wise.

We can also define the PV equation in vectorized format (analogous to equation (6)) as

\vec{\boldsymbol{q}} = \left(\boldsymbol{\nabla}^2 - f_0^2\mathbf{M}\right)\psi + \beta y

(30)

Below, we will go over some key components of these two formulations.

Inter-Connected Layers¶

We define $\mathbf{B}$ as the matrix which connects the inter-layer interactions. This is a bi-diagonal matrix which only maps the corresponding layer $k$ and $k-1$ .

\mathbf{B} = \begin{bmatrix} \frac{1}{H_1} & \frac{-1}{H_1} & 0 & \ldots & \ldots \\ 0 & \frac{1}{H_1} & \frac{-1}{H_2} & \ldots & \ldots \\ \ldots & \ldots & \ldots & \ldots & \ldots \\ \ldots & \ldots & \frac{-1}{H_{N-1} } & \frac{1}{H_{N-1}} & 0 \\ \ldots & \ldots& \ldots & \frac{-1}{H_N} & \frac{1}{H_N} \\ \end{bmatrix}

(31)

We also have $\mathbf{M}$ as the matrix connecting the layers of the stream function

\mathbf{M} = \begin{bmatrix} \frac{1}{H_1 g_1'} & \frac{-1}{H_1 g_2'} & \ldots & \ldots & \ldots \\ \frac{-1}{H_2 g_1'} & \frac{1}{H_1}\left(\frac{1}{g_1'} + \frac{1}{g_2'} \right) & \frac{-1}{H_2 g_2'} & \ldots & \ldots \\ \ldots & \ldots & \ldots & \ldots & \ldots \\ \ldots & \ldots & \frac{-1}{H_{n-1} g_{n-2}'} & \frac{1}{H_{n-1}}\left(\frac{1}{g_{n-2}'} + \frac{1}{g_{n-1}'} \right) & \frac{-1}{H_{n-1} g_{n-2}'} \\ \ldots & \ldots& \ldots & \frac{-1}{H_n g_{n-1}'} & \frac{1}{H_n g_{n-1}'} \\ \end{bmatrix}

(32)

This matrix is a tri-diagonal matrix as each layer, $k$ , within the stacked layers will influence or be influenced by the neighbours, $k-1,k+1$

Forcing¶

We have a forcing vector/matrix $\mathbf{F}$ which defines the forcing we apply at each layer. For example, the top layer could have wind forcing and/or atmospheric forcing. The bottom layer could have bottom friction due to the topography.

\mathbf{F} = \left[ F_0, 0, \ldots, 0, F_N\right]^\top

(33)

where we have the wind stress and the bottom friction. The wind forcing could be

\begin{aligned} F_0 =& \partial_x\tau_y - \partial_y\tau_x \\ \vec{\boldsymbol{\tau}} =& \frac{\tau_0}{\rho_0 H_1} \left[ -\cos(2\pi y/L_y), 0\right]^\top \end{aligned}

(34)

and the bottom friction could be described as

\begin{aligned} F_N =& \frac{\delta_{ek}}{2H_N} \boldsymbol{\nabla}\psi_N \end{aligned}

(35)

where $\delta_{ek}$ is the Ekman coefficient.

Dissipation¶

We have a very similar dissipation stategy as listed above in .... Here, we define this as

\begin{aligned} \text{Viscosity}: && && D_1 &= a_2 \boldsymbol{\Delta}^2\psi \\ \text{HyperViscosity}: && && D_2 &= - a_4 \boldsymbol{\Delta}^3\psi \end{aligned}

(36)

Parameter Configuration¶

In the paper [Thiry et al., 2023], they use the following formulation

\partial_t \vec{\boldsymbol{q}} + \vec{\boldsymbol{u}}\cdot\boldsymbol{\nabla} \vec{\boldsymbol{q}} = -a_2\boldsymbol{\Delta}_H^2\psi -a_4\boldsymbol{\Delta}_H^3\psi + \frac{\tau_0}{\rho_0H_1}\left[\partial_x\tau_y - \partial_y\tau_x, 0\cdots,0\right] - \frac{\delta_{ek}}{2H_{N_Z}} \left[0,\cdots,0,\Delta\psi_N\right]

(37)

Below is a table with the parameter configuration for their experiments.

Examples¶

Idealized QG (TODO)¶

We have a jupyter notebook showscasing how we can reconstruct the fields displayed in the above section using the idealized QG equation (7). This is a simple 1 Layer QG model on a periodic domain. We will use the experimental configuration found in Table 1 which will showcase the different flow regimes depending upon the parameters, i.e. decay flow, forced flow and $\beta$ -plane flow.

Realistic Idealized QG (TODO)¶

We have a jupyter notebook showcasing how we can reconstruct the fields displayed in the above section for the realistic idealized QG equation. This is another simple 1 Layer QG model but with parameters that mean something in real space.

Free-Run QG (TODO)¶

This is a simple 1.5 layer QG and the mapping problem for SSH. This comes from Guillou et al., 2021 paper where they showcase the relationship between SSH and the QG model.

Stacked QG (TODO)¶

This is an example that implements an example for the stacked QG model. The underlying mathematics is based on the Q-GCM numerical model. However, the configurations are based upon the papers [Thiry et al. (2022)Thiry et al. (2023)].

References

Amraoui, S., and Didier Auroux, Blum, J., & and, E. C. (2023). Back-and-forth nudging for the quasi-geostrophic ocean dynamics with altimetry: Theoretical convergence study and numerical experiments with the future SWOT observations. Discrete and Continuous Dynamical Systems - S, 16(2), 197–219. 10.3934/dcdss.2022058
Frezat, H., Sommer, J. L., Fablet, R., Balarac, G., & Lguensat, R. (2022). A posteriori learning for quasi-geostrophic turbulence parametrization. 10.48550/ARXIV.2204.03911
Srinivasan, K., Chekroun, M. D., & McWilliams, J. C. (2023). Turbulence closure with small, local neural networks: Forced two-dimensional and β-plane flows. arXiv. 10.48550/ARXIV.2304.05029
Gupta, J. K., & Brandstetter, J. (2022). Towards Multi-spatiotemporal-scale Generalized PDE Modeling. arXiv. 10.48550/ARXIV.2209.15616
Guan, Y., Subel, A., Chattopadhyay, A., & Hassanzadeh, P. (2023). Learning physics-constrained subgrid-scale closures in the small-data regime for stable and accurate LES. Physica D: Nonlinear Phenomena, 443, 133568. 10.1016/j.physd.2022.133568