On the non-hydrostatic equations in pressure and sigma coordinates

The non-hydrostatic equations governing the inviscid, adiabatic motion of a perfect gas are formulated using pressure as vertical coordinate; and M. J. Miller's 1974 approximate quasi-non-hydrostatic pressure coordinate equations are derived by applying a systematic scale analysis and power series expansion. The derivation makes clear that these equations are the pressure coordinate counterparts of the anelastic height coordinate equations obtained by Y. Ogura and N. A. Phillips in 1962. The two sets cannot be interconverted by coordinate transformation and so they are not physically equivalent; but the differences are small at the order of validity of both sets. Consideration of a quasi-hydrostatic approximation emphasizes the non-hydrostatic character of Miller's equations. Sigma coordinate quasi-non-hydrostatic equations are obtained by direct transformation of the pressure coordinate forms, and consistent energy equations are derived for both sets. Convenient diagnostic partial differential equations for the geopotential field are obtained for both pressure and sigma coordinate forms. As shown by Miller, the quasi-non-hydrostatic formulation does not permit vertically propagating acoustic waves. Horizontally propagating acoustic waves (Lamb waves) are in general allowed, but can be removed from the pressure coordinate system by applying suitable boundary conditions. Some aspects of the treatment of the Lamb wave problem are corrected in this study. The quasi-non-hydrostatic sigma coordinate system permits Lamb waves, but it may still be considered suitable for convective and (especially) mesoscale modelling with or without orography. The possible use of the quasi-non-hydrostatic system in large-scale theory and modelling is also discussed.

[1]  On the use of pressure as vertical co‐ordinate in modelling convection , 1974 .

[2]  M. Moncrieff,et al.  The dynamics and simulation of tropical cumulonimbus and squall lines , 1976 .

[3]  T. Clark A small-scale dynamic model using a terrain-following coordinate transformation , 1977 .

[4]  R. P. Pearce,et al.  A three‐dimensional primitive equation model of cumulonimbus convection , 1974 .

[5]  N. A. Phillips,et al.  A COORDINATE SYSTEM HAVING SOME SPECIAL ADVANTAGES FOR NUMERICAL FORECASTING , 1957 .

[6]  George J. Haltiner,et al.  Numerical weather prediction , 1971 .

[7]  N. Phillips,et al.  Scale Analysis of Deep and Shallow Convection in the Atmosphere , 1962 .

[8]  M. Miller,et al.  The Hampstead storm: A numerical simulation of a quasi-stationary cumulonimbus system , 1978 .

[9]  Leslie Farnell Solution of poisson equations on a nonuniform grid , 1980 .

[10]  A. Gilchrist,et al.  A general circulation model of the atmosphere suitable for long period integrations , 1972 .

[11]  M. Miller,et al.  Numerical simulations showing the role of the downdraught in cumulonimbus motion and splitting , 1978 .

[12]  A. W. Nielsen A study of the baroclinic instability problem and the lower boundary condition , 1971 .

[13]  Gareth P. Williams Numerical integration of the three-dimensional Navier-Stokes equations for incompressible flow , 1969, Journal of Fluid Mechanics.

[14]  M. Moncrieff,et al.  Dynamical models of two‐dimensional downdraughts , 1980 .

[15]  Richard C. J. Somerville,et al.  Numerical solution of the navier-stokes equations with topography , 1975 .

[16]  I. Orlanski The Quasi-Hydrostatic Approximation , 1981 .

[17]  Richard C. J. Somerville,et al.  On the use of a coordinate transformation for the solution of the Navier-Stokes equations , 1975 .

[18]  A. Kasahara Various Vertical Coordinate Systems Used for Numerical Weather Prediction , 1974 .

[19]  F. Bretherton Lamb waves in a nearly isothermal atmosphere , 1969 .

[20]  Mitchell W. Moncrieff,et al.  Two‐dimensional convection in non‐constant shear: A model of mid‐latitude squall lines , 1982 .

[21]  Traveling Convective Storms over Venezuela , 1977 .