Hydrostatic, Quasihydrostatic, and Nonhydrostatic Ocean Modeling

Total Page:16

File Type:pdf, Size:1020Kb

Hydrostatic, Quasihydrostatic, and Nonhydrostatic Ocean Modeling JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 102, NO. C3, PAGES 5733-5752, MARCH 15, 1997 Hydrostatic, quasi-hydrostatic, and nonhydrostatic ocean modeling John Marshall, Chris Hill, Lev Perelman, and Alistair Adcroft Center for Meteorologyand PhysicalOceanography, Department of Earth, Atmosphericand Planetary Sciences,Massachusetts Institute of Technology,Cambridge Abstract. Ocean modelsbased on consistenthydrostatic, quasi-hydrostatic, and nonhydrostaticequation setsare formulated and discussed.The quasi-hydrostaticand nonhydrostaticsets are more accuratethan the widely used hydrostaticprimitive equations.Quasi-hydrostatic models relax the precisebalance between gravity and pressuregradient forcesby includingin a consistentmanner cosine-of-latitudeCoriolis termswhich are neglectedin primitive equation models.Nonhydrostatic models employ the full incompressibleNavier Stokesequations; they are required in the studyof small- scalephenomena in the oceanwhich are not in hydrostaticbalance. We outline a solution strategyfor the Navier Stokesmodel on the spherethat performs efficientlyacross the whole range of scalesin the ocean,from the convectivescale to the global scale,and so leadsto a model of great versatility.In the hydrostaticlimit the Navier Stokesmodel involvesno more computationaleffort than those modelswhich assumestrict hydrostatic balanceon all scales.The strategyis illustrated in simulationsof laboratory experimentsin rotating convectionon scalesof a few centimeters,simulations of convectiveand baroclinicinstability of the mixed layer on the 1- to 10-km scale,and simulationsof the global circulation of the ocean. 1. Introduction Indeed, there are many important phenomena in the ocean, for example,wind- and buoyancy-driventurbulence in the sur- The ocean is a stratified fluid on a rotating Earth driven face mixedlayers of the ocean(on scalesL < 1 km), which are from its upper surfaceby patternsof momentum and buoyancy fundamentallynonhydrostatic and so cannot be studied using fluxes.The detailed dynamicsare very accuratelydescribed by hydrostaticmodels. the Navicr Stokes equations.These equations admit, and the In the present study wc outline, discuss,and illustrate the ocean contains,a wide variety of phenomena on a plethora of use of models based on equation sets that are more accurate spacescales and timescales.Modeling of the ocean is a formi- than the HPEs: quasi-hydrostaticmodels (QH), in which the dable challenge;it is a turbulent fluid containingenergetically precisebalance betweengravity and pressuregradient forces is active scalesranging from the global down to order 1-10 km relaxed, and fully nonhydrostaticmodels (NH), in which thc horizontallyand some tens of meters vertically; see Figure 1. incompressibleNavier Stokesequations are employed.Quasi- Important scale interactionsoccur over the entire spectrum. hydrostaticmodels treat the Coriolis force exactly,by including Numerical models of the ocean circulation, and the ocean in a consistentmanner cos (latitude) Coriolis terms that are models used in climate research, are rooted in the Navier conventionallyneglected in the HPEs. These cos (latitude) Stokes equations but employ approximate forms. Most are terms can becomc significant,particularly as the equator is based on the "hydrostatic primitive equations" (HPEs) in approached,and their inclusionendows the model with a com- which the vertical momentum equation is reduced to a state- plete angularmomentum principle. Nonhydrostatic models arc ment of hydrostaticbalance and the •'traditional approxima- required in the study of the smallestscales in the ocean. In tion" is made in which the Coriolis force is treated approxi- principle,of course,NH is also applicableon the largestscales; mately and the shallow atmosphereapproximation is made; we will demonstratethat modelsbased on algorithmsrooted in see section2. On the •'large scale" the terms omitted in the the Navier Stokes equations can be made etficient and used HPEs are generallythought to be small, but on "small scales" with economy even in the hydrostaticregime, leading to a the scalingassumptions implicit in them become increasingly singlealgorithm that can be employedacross the whole range problematical. of scalesdepicted in Figure 1. The global circulationof the oceanand its great wind-driven Our considerationsin this paper are independentof a par- gyres(on scalesL •--1000 km) are very accuratelydescribed by ticular numerical rendition or discretization.The strategy set the HPEs, as are the geostrophiceddies and rings associated out here could be employed in any model. From a single algorithmic base rooted in the Navier Stokes equations, NH, with its hydrodynamicalinstability (L --• 10-100 km). The QH, and HPE models are outlined. In a companion paper HPEs presumablybegin to break down somewherebetween 10 and I km, as the horizontal scale of the motion becomes [Marshallet al., this issue]we describethe details of a finite- volume, incompressibleNavier Stokes model which imple- comparablewith its vertical scale,the "grey area" in Figure 1. ments the ideas set out here. Copyright 1997 by the American GeophysicalUnion. In section 2 we critique the HPEs and review the assump- Paper number 96JC02776. tions made in their derivation, assessingtheir validity across 0148-0227/97/96J C-02776509.00 the range of scalesin the ocean. In section3 we write down the 5733 5734 MARSHALL ET AL.' HYDROSTATIC AND NONHYDROSTATIC OCEAN MODELING L ~ 1m I km 10 km 1,000 km 10,000 km which perform efficientlyacross the whole range of scalesin the ocean, from the convectiveto the global scale. Navier Stokesmodels are specificallydesigned for the studyof small- scalephenomenon such as convection.But when deployedto study hydrostaticscales, they need be no more demanding computationallythan hydrostaticmodels. Equation setsbased on more approximateforms, QH and HPE, are readily imple- gravity convectiongeostrophic gyres global waves 'eddies' circulation mentedas special cases of NH. A comparisonof integrationsof HPE, QH, and NH at large scale(1 ø horizontalresolution) givesessentially the samenumerical solutions. The neglectof cosqb Coriolis terms is the mostquestionable assumption made by the HPEs, but we find that their inclusion(in QH and NH) HYDROSTATIC yieldsdifferences in horizontalcurrents in oceangyres of only a fewmillimeters per second.Thus it is clearthat solutions based on the HPEs are not grosslyin error. Nevertheless, modelsbased on QH (orNH) shouldbe preferred in studiesof large-scalephenomenon. NON-HYDRQ..............................................................................................Finally, we have attemptedhere, as far as is possible,to present an accountwhich does not make strongassumptions Figure 1. A schematicdiagram showingthe range of scales about particular numerical choices. Details of the data- in the ocean. Global/basin-scalephenomenon are fundamen- parallel, finite-volume,incompressible Navier Stokes model tallyhydrostatic; convective processes on the kilometerscale usedhere to illustrateour ideasare givenby Marshall et al. [this are fundamentally nonhydrostatic.Somewhere between the issue].The mappingof the algorithmonto parallel machines, ge0strophicand convectivescales (10 km to 1 km) the hydro- in data-parallel FORTRAN on the Correction Machine staticapproximation breaks down: the "greyarea" in the fig- (CM5) and in the implicitlyparallel languageId on the data- ure. Modelswhich make the hydrostaticapproximation are flow machineMONSOON, is describedby Arvind et al. (A designedfor studyof large-scaleprocesses but are commonly usedat resolutionsthat encroach on thisgrey area (the left- comparisonof implicitly parallel multi-threaded and data- pointingarrow). Nonhydrostaticmodels based on the incom- parallel implementationsof an ocean model based on the pressibleNavier Stokesequations are valid acrossthe whole Navier Stokesequations, submitted to Journalof Parallel and range of scalesbut in oceanographyhave been hitherto used Distributed Computing, 1996) (hereinafter referred to as for processstudies on the convectivescale. We show in this Arvind et al., submittedmanuscript, 1996). paper that the computationaloverhead incurred by employing the unapproximatedequations is slightin the hydrostaticlimit. Thus modelsrooted in the Navier Stokesequations can alsobe 2. Critique of the Hydrostatic Primitive used with economyfor study of the large scale (the right- Equations pointingarrow in thefigure). The HPEs, which assumea precise balance between the pressureand densityfields, are almostaxiomatic to many me- teorologistsand oceanographers.They are widely used in nu- Navier Stokesequations on the sphereand discusshydrostatic, merical weather forecastingand climate simulationsof both quasi-hydrostatic,and nonhydrostaticregimes. Section 4 is the atmosphereand ocean. The terms omitted from the full concernedwith the diagnosisof the pressurefield required to Navier Stokesequations are customarilythought to be smallon ensurethat the evolvingvelocity field remainsnondivergent. In large scales(see Lorenz [1967] and Phillips[1973] for good HPE and QH a two-dimensional(2-D) elliptic equationmust discussions).However, the HPEs precludethe studyof non- be inverted for the surface pressure; in NH a three- hydrostaticsmall-scale phenomena, such as deep convection, dimensional(3-D) ellipticequation must be invertedsubject to the understandingof whichis of greatimportance to climate. Neumannboundary conditions. A strategyis developedfor the The HPEs can alsobe
Recommended publications
  • Arxiv:1902.03186V3 [Math.AP]
    PRIMITIVE EQUATIONS WITH HORIZONTAL VISCOSITY: THE INITIAL VALUE AND THE TIME-PERIODIC PROBLEM FOR PHYSICAL BOUNDARY CONDITIONS AMRU HUSSEIN, MARTIN SAAL, AND MARC WRONA Abstract. The 3D-primitive equations with only horizontal viscosity are con- sidered on a cylindrical domain Ω = (−h,h) × G, G ⊂ R2 smooth, with the physical Dirichlet boundary conditions on the sides. Instead of considering a vanishing vertical viscosity limit, we apply a direct approach which in particu- lar avoids unnecessary boundary conditions on top and bottom. For the initial value problem, we obtain existence and uniqueness of local z-weak solutions for initial data in H1((−h,h),L2(G)) and local strong solutions for initial data 1 1 2 q in H (Ω). If v0 ∈ H ((−h,h),L (G)), ∂zv0 ∈ L (Ω) for q > 2, then the z- weak solution regularizes instantaneously and thus extends to a global strong solution. This goes beyond the global well-posedness result by Cao, Li and Titi (J. Func. Anal. 272(11): 4606-4641, 2017) for initial data near H1 in the periodic setting. For the time-periodic problem, existence and uniqueness of z-weak and strong time periodic solutions is proven for small forces. Since this is a model with hyperbolic and parabolic features for which classical results are not directly applicable, such results for the time-periodic problem even for small forces are not self-evident. 1. Introduction and main results The 3D-primitive equations are one of the fundamental models for geophysical flows, and they are used for describing oceanic and atmospheric dynamics. They are derived from the Navier-Stokes equations assuming a hydrostatic balance.
    [Show full text]
  • Lecture 18 Ocean General Circulation Modeling
    Lecture 18 Ocean General Circulation Modeling 9.1 The equations of motion: Navier-Stokes The governing equations for a real fluid are the Navier-Stokes equations (con­ servation of linear momentum and mass mass) along with conservation of salt, conservation of heat (the first law of thermodynamics) and an equation of state. However, these equations support fast acoustic modes and involve nonlinearities in many terms that makes solving them both difficult and ex­ pensive and particularly ill suited for long time scale calculations. Instead we make a series of approximations to simplify the Navier-Stokes equations to yield the “primitive equations” which are the basis of most general circu­ lations models. In a rotating frame of reference and in the absence of sources and sinks of mass or salt the Navier-Stokes equations are @ �~v + �~v~v + 2�~ �~v + g�kˆ + p = ~ρ (9.1) t r · ^ r r · @ � + �~v = 0 (9.2) t r · @ �S + �S~v = 0 (9.3) t r · 1 @t �ζ + �ζ~v = ω (9.4) r · cpS r · F � = �(ζ; S; p) (9.5) Where � is the fluid density, ~v is the velocity, p is the pressure, S is the salinity and ζ is the potential temperature which add up to seven dependent variables. 115 12.950 Atmospheric and Oceanic Modeling, Spring '04 116 The constants are �~ the rotation vector of the sphere, g the gravitational acceleration and cp the specific heat capacity at constant pressure. ~ρ is the stress tensor and ω are non-advective heat fluxes (such as heat exchange across the sea-surface).F 9.2 Acoustic modes Notice that there is no prognostic equation for pressure, p, but there are two equations for density, �; one prognostic and one diagnostic.
    [Show full text]
  • Uniqueness of Solution for the 2D Primitive Equations with Friction Condition on the Bottom∗ 1 Introduction and Motivation
    Monograf´ıas del Semin. Matem. Garc´ıa de Galdeano. 27: 135{143, (2003). Uniqueness of solution for the 2D Primitive Equations with friction condition on the bottom∗ D. Breschy, F. Guill´en-Gonz´alezz, N. Masmoudix, M. A. Rodr´ıguez-Bellido{ Abstract Uniqueness of solution for the Primitive Equations with Dirichlet conditions on the bottom is an open problem even in 2D domains. In this work we prove a result of additional regularity for a weak solution v for the Primitive Equations when we replace Dirichlet boundary conditions by friction conditions. This allows to obtain uniqueness of weak solution global in time, for such a system [3]. Indeed, we show weak regularity for the vertical derivative of the solution, @zv for all time. This is because this derivative verifies a linear pde of convection-diffusion type with convection velocity v, and the pressure belongs to a L2-space in time with values in a weighted space. Keywords: Boundary conditions of type Navier, 2D Primitive Equations, uniqueness AMS Classification: 35Q30, 35B40, 76D05 1 Introduction and motivation. Primitive Equations are one of the models used to forecast the fluid velocity and pressure in the ocean. Such equations are obtained from the dimensionless Navier-Stokes equations, letting the aspect ratio (quotient between vertical dimension and horizontal dimensions) go to zero. The first results about existence of solution (weak, in the sense of the Navier- Stokes equations) are proved for boundary conditions of Dirichlet type on the bottom of the domain and with wind traction on the surface, in the works by Lions-Temam-Wang, ∗The second and fourth authors have been financed by the C.I.C.Y.T project MAR98-0486.
    [Show full text]
  • Shallow-Water Equations and Related Topics
    CHAPTER 1 Shallow-Water Equations and Related Topics Didier Bresch UMR 5127 CNRS, LAMA, Universite´ de Savoie, 73376 Le Bourget-du-Lac, France Contents 1. Preface .................................................... 3 2. Introduction ................................................. 4 3. A friction shallow-water system ...................................... 5 3.1. Conservation of potential vorticity .................................. 5 3.2. The inviscid shallow-water equations ................................. 7 3.3. LERAY solutions ........................................... 11 3.3.1. A new mathematical entropy: The BD entropy ........................ 12 3.3.2. Weak solutions with drag terms ................................ 16 3.3.3. Forgetting drag terms – Stability ............................... 19 3.3.4. Bounded domains ....................................... 21 3.4. Strong solutions ............................................ 23 3.5. Other viscous terms in the literature ................................. 23 3.6. Low Froude number limits ...................................... 25 3.6.1. The quasi-geostrophic model ................................. 25 3.6.2. The lake equations ...................................... 34 3.7. An interesting open problem: Open sea boundary conditions .................... 41 3.8. Multi-level and multi-layers models ................................. 43 3.9. Friction shallow-water equations derivation ............................. 44 3.9.1. Formal derivation ....................................... 44 3.10.
    [Show full text]
  • Vortex Stretching in Incompressible and Compressible Fluids
    Vortex stretching in incompressible and compressible fluids Esteban G. Tabak, Fluid Dynamics II, Spring 2002 1 Introduction The primitive form of the incompressible Euler equations is given by du P = ut + u · ∇u = −∇ + gz (1) dt ρ ∇ · u = 0 (2) representing conservation of momentum and mass respectively. Here u is the velocity vector, P the pressure, ρ the constant density, g the acceleration of gravity and z the vertical coordinate. In this form, the physical meaning of the equations is very clear and intuitive. An alternative formulation may be written in terms of the vorticity vector ω = ∇× u , (3) namely dω = ω + u · ∇ω = (ω · ∇)u (4) dt t where u is determined from ω nonlocally, through the solution of the elliptic system given by (2) and (3). A similar formulation applies to smooth isentropic compressible flows, if one replaces the vorticity ω in (4) by ω/ρ. This formulation is very convenient for many theoretical purposes, as well as for better understanding a variety of fluid phenomena. At an intuitive level, it reflects the fact that rotation is a fundamental element of fluid flow, as exem- plified by hurricanes, tornados and the swirling of water near a bathtub sink. Its derivation from the primitive form of the equations, however, often relies on complex vector identities, which render obscure the intuitive meaning of (4). My purpose here is to present a more intuitive derivation, which follows the tra- ditional physical wisdom of looking for integral principles first, and only then deriving their corresponding differential form. The integral principles associ- ated to (4) are conservation of mass, circulation–angular momentum (Kelvin’s theorem) and vortex filaments (Helmholtz’ theorem).
    [Show full text]
  • Studies on Blood Viscosity During Extracorporeal Circulation
    Nagoya ]. med. Sci. 31: 25-50, 1967. STUDIES ON BLOOD VISCOSITY DURING EXTRACORPOREAL CIRCULATION HrsASHI NAGASHIMA 1st Department of Surgery, Nagoya University School of Medicine (Director: Prof Yoshio Hashimoto) ABSTRACT Blood viscosity was studied during extracorporeal circulation by means of a cone in cone viscometer. This rotational viscometer provides sixteen kinds of shear rate, ranging from 0.05 to 250.2 sec-1 • Merits and demerits of the equipment were described comparing with the capillary viscometer. In experimental study it was well demonstrated that the whole blood viscosity showed the shear rate dependency at all hematocrit levels. The influence of the temperature change on the whole blood viscosity was clearly seen when hemato· crit was high and shear rate low. The plasma viscosity was 1.8 cp. at 37°C and Newtonianlike behavior was observed at the same temperature. ·In the hypothermic condition plasma showed shear rate dependency. Viscosity of 10% LMWD was over 2.2 fold as high as that of plasma. The increase of COz content resulted in a constant increase in whole blood viscosity and the increased whole blood viscosity after C02 insufflation, rapidly decreased returning to a little higher level than that in untreated group within 3 minutes. Clinical data were obtained from 27 patients who underwent hypothermic hemodilution perfusion. In the cyanotic group, w hole blood was much more viscous than in the non­ cyanotic. During bypass, hemodilution had greater influence upon whole blood viscosity than hypothermia, but it went inversely upon the plasma viscosity. The whole blood viscosity was more dependent on dilution rate than amount of diluent in mljkg.
    [Show full text]
  • Global Well-Posedness of the Three-Dimensional Viscous Primitive Equations of Large Scale Ocean and Atmosphere Dynamics
    Annals of Mathematics, 166 (2007), 245–267 Global well-posedness of the three-dimensional viscous primitive equations of large scale ocean and atmosphere dynamics By Chongsheng Cao and Edriss S. Titi Abstract In this paper we prove the global existence and uniqueness (regularity) of strong solutions to the three-dimensional viscous primitive equations, which model large scale ocean and atmosphere dynamics. 1. Introduction Large scale dynamics of oceans and atmosphere is governed by the primi- tive equations which are derived from the Navier-Stokes equations, with rota- tion, coupled to thermodynamics and salinity diffusion-transport equations, which account for the buoyancy forces and stratification effects under the Boussinesq approximation. Moreover, and due to the shallowness of the oceans and the atmosphere, i.e., the depth of the fluid layer is very small in compar- ison to the radius of the earth, the vertical large scale motion in the oceans and the atmosphere is much smaller than the horizontal one, which in turn leads to modeling the vertical motion by the hydrostatic balance. As a result one obtains the system (1)–(4), which is known as the primitive equations for ocean and atmosphere dynamics (see, e.g., [20], [21], [22], [23], [24], [33] and references therein). We observe that in the case of ocean dynamics one has to add the diffusion-transport equation of the salinity to the system (1)–(4). We omitted it here in order to simplify our mathematical presentation. However, we emphasize that our results are equally valid when the salinity effects are taking into account. Note that the horizontal motion can be further approximated by the geostrophic balance when the Rossby number (the ratio of the horizontal ac- celeration to the Coriolis force) is very small.
    [Show full text]
  • The Influence of Negative Viscosity on Wind-Driven, Barotropic Ocean
    916 JOURNAL OF PHYSICAL OCEANOGRAPHY VOLUME 30 The In¯uence of Negative Viscosity on Wind-Driven, Barotropic Ocean Circulation in a Subtropical Basin DEHAI LUO AND YAN LU Department of Atmospheric and Oceanic Sciences, Ocean University of Qingdao, Key Laboratory of Marine Science and Numerical Modeling, State Oceanic Administration, Qingdao, China (Manuscript received 3 April 1998, in ®nal form 12 May 1999) ABSTRACT In this paper, a new barotropic wind-driven circulation model is proposed to explain the enhanced transport of the western boundary current and the establishment of the recirculation gyre. In this model the harmonic Laplacian viscosity including the second-order term (the negative viscosity term) in the Prandtl mixing length theory is regarded as a tentative subgrid-scale parameterization of the boundary layer near the wall. First, for the linear Munk model with weak negative viscosity its analytical solution can be obtained with the help of a perturbation expansion method. It can be shown that the negative viscosity can result in the intensi®cation of the western boundary current. Second, the fully nonlinear model is solved numerically in some parameter space. It is found that for the ®nite Reynolds numbers the negative viscosity can strengthen both the western boundary current and the recirculation gyre in the northwest corner of the basin. The drastic increase of the mean and eddy kinetic energies can be observed in this case. In addition, the analysis of the time-mean potential vorticity indicates that the negative relative vorticity advection, the negative planetary vorticity advection, and the negative viscosity term become rather important in the western boundary layer in the presence of the negative viscosity.
    [Show full text]
  • Circulation, Lecture 4
    4. Circulation A key integral conservation property of fluids is the circulation. For convenience, we derive Kelvin’s circulation theorem in an inertial coordinate system and later transform back to Earth coordinates. In inertial coordinates, the vector form of the momentum equation may be written dV = −α∇p − gkˆ + F, (4.1) dt where kˆ is the unit vector in the z direction and F represents the vector frictional acceleration. Now define a material curve that, however, is constrained to lie at all times on a surface along which some state variable, which we shall refer to for now as s,is a constant. (A state variable is a variable that can be expressed as a function of temperature and pressure.) The picture we have in mind is shown in Figure 3.1. 10 The arrows represent the projection of the velocity vector onto the s surface, and the curve is a material curve in the specific sense that each point on the curve moves with the vector velocity, projected onto the s surface, at that point. Now the circulation is defined as C ≡ V · dl, (4.2) where dl is an incremental length along the curve and V is the vector velocity. The integral is a closed integral around the curve. Differentiation of (4.2) with respect to this gives dC dV dl = · dl + V · . (4.3) dt dt dt Since the curve is a material curve, dl = dV, dt so the integrand of the last term in (4.3) can be written as a perfect derivative and so the term itself vanishes.
    [Show full text]
  • 4. Circulation and Vorticity This Chapter Is Mainly Concerned With
    4. Circulation and vorticity This chapter is mainly concerned with vorticity. This particular flow property is hard to overestimate as an aid to understanding fluid dynamics, essentially because it is difficult to mod- ify. To provide some context, the chapter begins by classifying all different kinds of motion in a two-dimensional velocity field. Equations are then developed for the evolution of vorticity in three dimensions. The prize at the end of the chapter is a fluid property that is related to vorticity but is even more conservative and therefore more powerful as a theoretical tool. 4.1 Two-dimensional flows According to a theorem of Helmholtz, any two-dimensional flowV()x, y can be decom- posed as V= zˆ × ∇ψ+ ∇χ , (4.1) where the “streamfunction”ψ()x, y and “velocity potential”χ()x, y are scalar functions. The first part of the decomposition is non-divergent and the second part is irrotational. If we writeV = uxˆ + v yˆ for the flow in the horizontal plane, then 4.1 says that u = –∂ψ⁄ ∂y + ∂χ ∂⁄ x andv = ∂ψ⁄ ∂x + ∂χ⁄ ∂y . We define the vertical vorticity ζ and the divergence D as follows: ∂v ∂u 2 ζ =zˆ ⋅ ()∇ × V =------ – ------ = ∇ ψ , (4.2) ∂x ∂y ∂u ∂v 2 D =∇ ⋅ V =------ + ----- = ∇ χ . (4.3) ∂x ∂y Determination of the velocity field givenζ and D requires boundary conditions onψ andχ . Contours ofψ are called “streamlines”. The vorticity is proportional to the local angular velocity aboutzˆ . It is known entirely fromψ()x, y or the first term on the rhs of 4.1.
    [Show full text]
  • Nonlinear Balance and Potential Vorticity Thinking at Large Rossby Number
    Nonlinear Balance and Potential Vorticity Thinking at Large Rossby Number D. J. Raymond Physics Department and Geophysical Research Center, New Mexico Institute of Mining and Technology, Socorro, NM 87801 SUMMARY Two rational approximations are made to the divergence, potential vorticity, and potential temperature equations resulting in two different nonlinear balance models. Semibalance is similar (but not identical to) the nonlinear balance model of Lorenz. Quasibalance is a simpler model that is equivalent to quasigeostrophy at low Rossby number and the barotropic model at high Rossby number. Practical solution methods that work for all Rossby number are outlined for both models. A variety of simple initial value problems are then solved with the aim of fortifying our insight into the behaviour of large Rossby number flows. The flows associated with a potential vorticity anomaly at very large Rossby number differ in significant ways from the corresponding low Rossby number results. In particular, an isolated anomaly has zero vertical radius of influence at infinite Rossby number, while the induced tangential velocity in a horizontal plane containing the anomaly decreases as the inverse of radius rather than the inverse square. The effects of heating and frictional forces are approached from a point of view that is somewhat different from that recently expressed by Haynes and McIntyre, though the physical content is the same. 1 Introduction In a highly influential paper, Hoskins, McIntyre, and Robertson (1985) (hereafter HMR – see also Hoskins, McIntyre, and Robertson, 1987) presented a conceptual framework for thinking about large scale flow patterns that are nearly balanced. In this framework potential vorticity plays a central role.
    [Show full text]
  • Different Forms of the Governing Equations for Atmospheric Motions
    Different Forms of the Governing Equations for Atmospheric Motions Dale Durran Dept. of Atmospheric Sciences University of Washington Outline • Conservation relations and approximate equations for motion on the sphere • Nonhydrostatic “deep” equations • Hydrostatic Primitive equations • Filtering acoustic waves in nonhydrostatic models • Boussinesq approximation • Anelastic and Pseudo-incompressible systems • Application to global models Approximating the Equations on the Sphere Spherical Coordinate System Nonhydrostatic Deep Momentum Equations Material derivative is: • r is the distance to the center of the earth • All Coriolis terms retained • All metric terms retained Conservation Relations for NDE Axial angular momentum Ertel potential vorticity Total energy per unit mass Hydrostatic Primitive Equations Material derivative is: • hydrostatic approximation • z replaces r as the vertical coordinate • a replaces r as the distance to the center of the earth • cosφ Coriolis terms neglected • only two of six metric terms retained Conservation Relations for HPE Axial angular momentum Potential vorticity Total energy per unit mass Neglect of cosφ Coriolis terms • Required to obtain consistent conservation relations in the primitive equations. • “Subject of quiet controversy in meteorology and oceanography for many years.” (White, et al., 2005) • Not important according to linearized adiabatic analyses. • Various scaling analyses suggest they may nevertheless be non-negligible, particularly in the tropics. Simple Argument for Including the cosφ
    [Show full text]