1.4.1. Hydrostatic Primitive Equations for the Atmosphere in Pressure Coordinates

The hydrostatic primitive equations (HPE’s) in \(p-\)coordinates are:

(1.45)\[\frac{D\vec{\mathbf{v}}_{h}}{Dt}+f\hat{\boldsymbol{k}}\times \vec{\mathbf{v}}_{h}+ \nabla _{p}\phi = \vec{\boldsymbol{\mathcal{F}}}\]
(1.46)\[\frac{\partial \phi }{\partial p}+\alpha = 0\]
(1.47)\[\nabla _{p}\cdot \vec{\mathbf{v}}_{h}+\frac{\partial \omega }{\partial p} = 0\]
(1.48)\[p\alpha = RT\]
(1.49)\[c_{v}\frac{DT}{Dt}+p\frac{D\alpha }{Dt} = \mathcal{Q}\]

where \(\vec{\mathbf{v}}_{h}=(u,v,0)\) is the ‘horizontal’ (on pressure surfaces) component of velocity, \(\frac{D}{Dt}=\frac{\partial}{\partial t}+\vec{\mathbf{v}}_{h}\cdot \nabla _{p}+\omega \frac{\partial }{\partial p}\) is the total derivative, \(f=2\Omega \sin \varphi\) is the Coriolis parameter, \(\phi =gz\) is the geopotential, \(\alpha =1/\rho\) is the specific volume, \(\omega =\frac{Dp }{Dt}\) is the vertical velocity in the \(p-\)coordinate. Equation (1.49) is the first law of thermodynamics where internal energy \(e=c_{v}T\), \(T\) is temperature, \(Q\) is the rate of heating per unit mass and \(p\frac{D\alpha }{Dt}\) is the work done by the fluid in compressing.

It is convenient to cast the heat equation in terms of potential temperature \(\theta\) so that it looks more like a generic conservation law. Differentiating (1.48) we get:

\[p\frac{D\alpha }{Dt}+\alpha \frac{Dp}{Dt}=R\frac{DT}{Dt}\]

which, when added to the heat equation (1.49) and using \(c_{p}=c_{v}+R\), gives:

(1.50)\[c_{p}\frac{DT}{Dt}-\alpha \frac{Dp}{Dt}=\mathcal{Q}\]

Potential temperature is defined:

(1.51)\[\theta =T(\frac{p_{c}}{p})^{\kappa }\]

where \(p_{c}\) is a reference pressure and \(\kappa =R/c_{p}\). For convenience we will make use of the Exner function \(\Pi (p)\) which is defined by:

(1.52)\[\Pi (p)=c_{p}(\frac{p}{p_{c}})^{\kappa }\]

The following relations will be useful and are easily expressed in terms of the Exner function:

\[c_{p}T=\Pi \theta \;\;;\;\;\frac{\partial \Pi }{\partial p}=\frac{\kappa \Pi }{p}\;\;;\;\;\alpha =\frac{\kappa \Pi \theta }{p}=\frac{\partial \ \Pi }{ \partial p}\theta \;\;;\;\;\frac{D\Pi }{Dt}=\frac{\partial \Pi }{\partial p} \frac{Dp}{Dt}\]

where \(b=\frac{\partial \ \Pi }{\partial p}\theta\) is the buoyancy.

The heat equation is obtained by noting that

\[c_{p}\frac{DT}{Dt}=\frac{D(\Pi \theta )}{Dt}=\Pi \frac{D\theta }{Dt}+\theta \frac{D\Pi }{Dt}=\Pi \frac{D\theta }{Dt}+\alpha \frac{Dp}{Dt}\]

and on substituting into (1.50) gives:

(1.53)\[\Pi \frac{D\theta }{Dt}=\mathcal{Q}\]

which is in conservative form.

For convenience in the model we prefer to step forward (1.53) rather than (1.49).

1.4.1.1. Boundary conditions

The upper and lower boundary conditions are:

(1.54)\[\begin{aligned}\mbox{at the top:}\;\;p=0 &\text{, }\omega =\frac{Dp}{Dt}=0\end{aligned}\]
(1.55)\[\begin{aligned}\mbox{at the surface:}\;\;p=p_{s} &\text{, }\phi =\phi _{\rm topo}=g~Z_{\rm topo}\end{aligned}\]

In \(p-\)coordinates, the upper boundary acts like a solid boundary (\(\omega=0\) ); in \(z-\)coordinates the lower boundary is analogous to a free surface (\(\phi\) is imposed and \(\omega \neq 0\)).

1.4.1.2. Splitting the geopotential

For the purposes of initialization and reducing round-off errors, the model deals with perturbations from reference (or ‘standard’) profiles. For example, the hydrostatic geopotential associated with the resting atmosphere is not dynamically relevant and can therefore be subtracted from the equations. The equations written in terms of perturbations are obtained by substituting the following definitions into the previous model equations:

(1.56)\[\theta = \theta _{o}+\theta ^{\prime }\]
(1.57)\[\alpha = \alpha _{o}+\alpha ^{\prime }\]
(1.58)\[\phi = \phi _{o}+\phi ^{\prime }\]

The reference state (indicated by subscript ‘o’) corresponds to horizontally homogeneous atmosphere at rest (\(\theta _{o},\alpha _{o},\phi_{o}\)) with surface pressure \(p_{o}(x,y)\) that satisfies \(\phi_{o}(p_{o})=g~Z_{\rm topo}\), defined:

\[\begin{split}\theta _{o}(p) = f^{n}(p) \\\end{split}\]
\[\begin{split}\alpha _{o}(p) = \Pi _{p}\theta _{o} \\\end{split}\]
\[\phi _{o}(p) = \phi _{\rm topo}-\int_{p_{0}}^{p}\alpha _{o}dp\]

The final form of the HPE’s in \(p-\)coordinates is then:

(1.59)\[\frac{D\vec{\mathbf{v}}_{h}}{Dt}+f\hat{\boldsymbol{k}}\times \vec{\mathbf{v}} _{h}+ \nabla _{p}\phi ^{\prime } = \vec{\boldsymbol{\mathcal{F}}}\]
(1.60)\[\frac{\partial \phi ^{\prime }}{\partial p}+\alpha ^{\prime } = 0\]
(1.61)\[\nabla _{p}\cdot \vec{\mathbf{v}}_{h}+\frac{\partial \omega }{ \partial p} = 0\]
(1.62)\[\frac{\partial \Pi }{\partial p}\theta ^{\prime } = \alpha ^{\prime }\]
(1.63)\[\frac{D\theta }{Dt} = \frac{\mathcal{Q}}{\Pi }\]