Goal: Use understanding of physically-relevant scales to reduce the complexity of the governing equations

Similar documents
Governing Equations and Scaling in the Tropics

ρ x + fv f 'w + F x ρ y fu + F y Fundamental Equation in z coordinate p = ρrt or pα = RT Du uv tanφ Dv Dt + u2 tanφ + vw a a = 1 p Dw Dt u2 + v 2

Dust devils, water spouts, tornados

The Equations of Motion in a Rotating Coordinate System. Chapter 3

Dynamic Meteorology - Introduction

g (z) = 1 (1 + z/a) = 1 1 ( km/10 4 km) 2

Dynamics Rotating Tank

2 A: The Shallow Water Equations

ESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr. DeCaria

On side wall labeled A: we can express the pressure in a Taylor s series expansion: x 2. + higher order terms,

1/18/2011. Conservation of Momentum Conservation of Mass Conservation of Energy Scaling Analysis ESS227 Prof. Jin-Yi Yu

10 Shallow Water Models

The dynamics of high and low pressure systems

Radiative equilibrium Some thermodynamics review Radiative-convective equilibrium. Goal: Develop a 1D description of the [tropical] atmosphere

Conservation of Mass Conservation of Energy Scaling Analysis. ESS227 Prof. Jin-Yi Yu

g (z) = 1 (1 + z/a) = 1

Circulation and Vorticity

Models of ocean circulation are all based on the equations of motion.

Chapter 1. Governing Equations of GFD. 1.1 Mass continuity

The Hydrostatic Approximation. - Euler Equations in Spherical Coordinates. - The Approximation and the Equations

7 Balanced Motion. 7.1 Return of the...scale analysis for hydrostatic balance! CSU ATS601 Fall 2015

CHAPTER 4. THE HADLEY CIRCULATION 59 smaller than that in midlatitudes. This is illustrated in Fig. 4.2 which shows the departures from zonal symmetry

+ ω = 0, (1) (b) In geometric height coordinates in the rotating frame of the Earth, momentum balance for an inviscid fluid is given by

7 The General Circulation

2 Transport of heat, momentum and potential vorticity

Atmospheric Dynamics: lecture 2

Reynolds Averaging. We separate the dynamical fields into slowly varying mean fields and rapidly varying turbulent components.

Class exercises Chapter 3. Elementary Applications of the Basic Equations

4. Atmospheric transport. Daniel J. Jacob, Atmospheric Chemistry, Harvard University, Spring 2017

ATMO/OPTI 656b Spring 09. Physical properties of the atmosphere

Atmospheric Circulation

EART164: PLANETARY ATMOSPHERES

where p oo is a reference level constant pressure (often 10 5 Pa). Since θ is conserved for adiabatic motions, a prognostic temperature equation is:

Chapter 2. The continuous equations

APPENDIX B. The primitive equations

Atmosphere, Ocean and Climate Dynamics Answers to Chapter 8

Hurricanes are intense vortical (rotational) storms that develop over the tropical oceans in regions of very warm surface water.

Scale Analysis of the Equations of Motion

1/18/2011. From the hydrostatic equation, it is clear that a single. pressure and height in each vertical column of the atmosphere.

Balanced Flow Geostrophic, Inertial, Gradient, and Cyclostrophic Flow

The atmosphere: A general introduction Niels Woetmann Nielsen Danish Meteorological Institute

Sample Final Questions: Solutions Math 21B, Winter y ( y 1)(1 + y)) = A y + B

Chapter 4: Fundamental Forces

ATM 298, Spring 2013 Lecture 2 The Equa;ons of Fluid Mo;on April 3, Paul A. Ullrich (HH 251)

Parcel Model. Meteorology September 3, 2008

Dynamics and Kinematics

Part 4. Atmospheric Dynamics

Chapter 2. Quasi-Geostrophic Theory: Formulation (review) ε =U f o L <<1, β = 2Ω cosθ o R. 2.1 Introduction

Parcel Model. Atmospheric Sciences September 30, 2012

Thermodynamics Review [?] Entropy & thermodynamic potentials Hydrostatic equilibrium & buoyancy Stability [dry & moist adiabatic]

Chapter 5. Shallow Water Equations. 5.1 Derivation of shallow water equations

1 Introduction to Governing Equations 2 1a Methodology... 2

Today s Lecture: Atmosphere finish primitive equations, mostly thermodynamics

The dynamics of a simple troposphere-stratosphere model

Fundamentals of Weather and Climate

Geophysics Fluid Dynamics (ESS228)

Dynamics of the Extratropical Response to Tropical Heating

z g + F w (2.56) p(x, y, z, t) = p(z) + p (x, y, z, t) (2.120) ρ(x, y, z, t) = ρ(z) + ρ (x, y, z, t), (2.121)

Calculus Math 21B, Winter 2009 Final Exam: Solutions

Chapter 1. Introduction

Chapter 5. Fundamentals of Atmospheric Modeling

Thermospheric Winds. Astrid Maute. High Altitude Observatory (HAO) National Center for Atmospheric Science (NCAR) Boulder CO, USA

EATS Notes 1. Some course material will be online at

Q.1 The most abundant gas in the atmosphere among inert gases is (A) Helium (B) Argon (C) Neon (D) Krypton

Chapter 9. Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation

References: Parcel Theory. Vertical Force Balance. ESCI Cloud Physics and Precipitation Processes Lesson 3 - Stability and Buoyancy Dr.

Fundamental Meteo Concepts

1. The vertical structure of the atmosphere. Temperature profile.

d v 2 v = d v d t i n where "in" and "rot" denote the inertial (absolute) and rotating frames. Equation of motion F =

Dynamics of the Zonal-Mean, Time-Mean Tropical Circulation

2.5 Shallow water equations, quasigeostrophic filtering, and filtering of inertia-gravity waves

2 Equations of Motion

Meteorology 6150 Cloud System Modeling

Quick Recapitulation of Fluid Mechanics

Isentropic Analysis. Much of this presentation is due to Jim Moore, SLU

Eliassen-Palm Theory

OCN/ATM/ESS 587. The wind-driven ocean circulation. Friction and stress. The Ekman layer, top and bottom. Ekman pumping, Ekman suction

2. Conservation laws and basic equations

Synoptic Meteorology I: Other Force Balances

AE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.

1 The circulation of a zonally symmetric atmosphere. 1.1 Angular momentum conservation and its implications

Quasi-equilibrium Theory of Small Perturbations to Radiative- Convective Equilibrium States

Using simplified vorticity equation,* by assumption 1 above: *Metr 430 handout on Circulation and Vorticity. Equations (4) and (5) on that handout

1/27/2010. With this method, all filed variables are separated into. from the basic state: Assumptions 1: : the basic state variables must

Lecture 10a: The Hadley Cell

t tendency advection convergence twisting baroclinicity

A few examples Shallow water equation derivation and solutions. Goal: Develop the mathematical foundation of tropical waves

Objective: To introduce the equations of motion and describe the forces that act upon the Atmosphere

Boundary layer equilibrium [2005] over tropical oceans

Synoptic-Dynamic Meteorology in Midlatitudes

An Optimal Control Problem Formulation for. the Atmospheric Large-Scale Wave Dynamics

Lecture 1. Equations of motion - Newton s second law in three dimensions. Pressure gradient + force force

Modelling of dispersed, multicomponent, multiphase flows in resource industries. Section 3: Examples of analyses conducted for Newtonian fluids

Free-Fall Timescale of Sun

A semi-implicit non-hydrostatic covariant dynamical kernel using spectral representation in the horizontal and a height based vertical coordinate

Project 3 Convection and Atmospheric Thermodynamics

The perturbation pressure, p, can be represented as the sum of a hydrostatic pressure perturbation p h and a nonhydrostatic pressure perturbation p nh

196 7 atmospheric oscillations:

1/3/2011. This course discusses the physical laws that govern atmosphere/ocean motions.

Lecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2)

Transcription:

Scale analysis relevant to the tropics [large-scale synoptic systems]* Goal: Use understanding of physically-relevant scales to reduce the complexity of the governing equations *Reminder: Midlatitude scale analysis in dynamics

Momentum Equations in Spherical Coordinates du uv tan# " + uw = " 1 %p $ %x + 2&v sin# " 2&w cos# + F rx dv + u2 tan# + vw = " 1 %p $ %y " 2&usin# + F ry dw " u2 + v 2 = " 1 %p $ %z + 2&ucos# " g + F rz total derivative pressure gradient Coriolis gravity friction To what extent are these terms important?

Scale Analysis [following Holton 11.2] Goal: To determine relative importance of the terms in the basic equations for particular scales of motion. Approach: Estimate the following quantities 1) The magnitude of the field variables. 2) The amplitudes of fluctuations in the field variables. [To estimate derivatives.] 3) The characteristic length, depth and time scales on which these fluctuations occur. But first a coordinate change

Vertical coordinate transformation It is often convenient to replace height z by pressure p. Consider p=p(x,z), where p(x,z) is assumed to be monotonic in z and x represents any horizontal coordinate, and the system of level curves sketched below: z p 3 p 2 Along each level curve, p(x,z) = constant. So let s consider the differential of p: p 1 dp = 0 =! p!x z dx +! p!z x dz!!z!x p = "! p /!x z! p /!z x x Under hydrostatic equilibrium:!z =! p /!x z!x p "g For further convenience, we introduce the geopotential: " = g Thus:! h " p =! #1! h p z The notation h denotes horizontal components of the gradient operator. z # z 0 dz

Conservation of mass in pressure coordinate Consider a parcel of mass δm which is conserved following the motion of the flow: 1 "m d("m) = 0 Assuming a density ρ, δm = ρδv, and: 1!"V d(!"v ) = 0 = g ("x"y" p) d! # " "x"y" p g $ & % Pressure velocity! dp (" 1! 0 = * $ )#!x d!x + 1!y d!y + 1 d! p! p % + " '- =!u &,!x +!v!y +!" % $ ' #! p & "!!u!x +!v!y +!" % $ ' = 0 #! p & No density, hence no explicit time-dependence appearing in the mass conservation relationship expressed in pressure in the vertical.

log pressure coordinates For the purposes of scale analysis, we ll consider the natural logarithm of pressure rather than pressure itself: " z* =!H ln p % $ ' # p s & where p s is a standard reference pressure and H is a scale height defined as: In this coordinate: w* = dz * H " R d T s g =!H d ln " p % $ ' =! H! # p s & p d " # #t + v $ % + w * # h #z *!"! p =!w *!z *! w * H For an isothermal atmosphere at temperature Ts [~global surface temperature], z*=z. However, for realistic temperature profiles in the troposphere, the difference between z* and z is usually small.

Scaling quantities Scale Symbol Value Horizontal velocity U 10 m s -1 Horizontal length L 10 6 m Vertical depth H 10 4 m Time τ~l/u 10 5 s Vertical velocity W TBD Geopotential fluctuation Temperature fluctuation δφ δt TBD TBD

Physical Parameters g "10 m s #2 gravity "10 7 m $ 0 =? f = 2%sin$ radius of earth We ll consider both middle and tropical latitudes. Coriolis parameter & =10 #5 m 2 s #1 viscosity

Continuity + Horizontal Equations "u "x + "v "y + "w * "z * # w * H = 0 "u "x + "v "y # U L In the horizontal, for nondivergent flow, h v=0. W W " UH H L To perform the horizontal scaling, we ll consider the magnitudes of terms relative to the horizontal advective scale: U 2 du " uv tan# + uw * = " $% $x + 2&v sin# " 2&w *cos# + F rx dv + u2 tan# + vw * = " $% $y " 2&usin# + F ry L Re is the Reynolds Number 1 L ~ 0.1 LW U " H ~ 10 #3 "# U 2 fl U = LW"cos#! Ro"1 $ LU = Re!1 ~ 10!12 U 2 Ro is the Rossby Number "H U ~ 0.1

1 Continuity + Horizontal Equations du " uv tan# + uw * = " $% $x + 2&v sin# " 2&w *cos# + F rx dv + u2 tan# + vw * = " $% $y " 2&usin# + F ry L ~ 0.1 LW U " H ~ 10 #3 "# U 2 fl U = $! Ro"1 U 2 LU = Re!1 ~ 10!12 For mid-latitudes (λ~45 ), f~10-4, so Ro ~ 0.1. Thus, the Coriolis term (~10) must be balanced by the geopotential gradient term, implying: "# ~ ful [~1000 m 2 s -2 ] At low [tropical] latitudes, f 10-5, so Ro 1, so Coriolis may not balance the pressure gradient force. In fact, for Ro 10, we expect: "# ~ U 2 [~100 m 2 s -2 ] LW"cos# "H U ~ 0.1 Thus, for synoptic disturbances in the tropics, geopoential perturbations are an order of magnitude smaller than for similar-sized systems in mid latitudes, which has several important consequences

Vertical equation! 1 " p! "z =! 1 #"" %! $ " p & ( '!1 "" "z =! g #""& % (! $ " p '!1!!! p =!!!z *!z *! p = " H p!!!z *!" 1 " p! "z = g H p $ "# ' & )! %"z *( "1 = gr d T H $ "# ' & ) %"z *( "1 dw * " u2 + v 2 = gr dt H % #$ ( ' * &#z *) "1 + 2+ucos, " g + F rz WU L " U 2 H L 2 ~ 10 #4 U 2 ~ 10 "5 gr d T H $ "# ' & ) %"z *( The above scaling implies hydrostatic equilibrium: But it also provides a scaling for horizontal temperature perturbations*: *1 U"cos# $10 %3 10 "# "z * = R T d H!T = H "!!!! ~ ~ U 2 R d "z * R d "W H 2 # "U HL ~ 10$14 R d ~ 0.3 K It is assumed that the geopotential consists of a mean component, which is a function of the vertical coordinate only, and a synoptic perturbation; since we re evaluating synoptic scale motion, the scaling applies to the perturbation part.

(Dry) Thermodynamic equation % " "t +v #$ ( ' h h * T + w * HN 2 & ) For δt~0.3 K, the first term on the left-hand side is of order 0.3K x 10-5 s -1 ~ 0.3K day -1. When precipitation is absent, the diabatic heating term on the right-hand side comprises long-wave radiative emission, which is observed to be of order 1 K day -1. Thus, R d = J c P w * HN 2 R d " J c p # w* " J c p R d HN 2 For the tropical tropopshere, N -1 ~100 s, so: w* ~ W ~ 0.3 cms "1 In the absence of precipitation, synoptic scale tropical vertical velocities are much smaller than in similar-sized extratropical systems. Also, from the continuity equation, the divergence of horizontal wind is of order 10-7 ; thus the flow is essentially nondivergent.

(Brief overview of) precipitating tropical synoptic systems What modifications are anticipated for precipitating tropical synoptic systems? For precipitating tropical synoptic systems, rain rates of 20 mm day -1 are typical. For a unit area cross-section, this precipitation rate implies a daily condensation of 20 kg of H 2 0. [Think density x area, which is dimensionally mass/length.] Since L c ~2.5 x 10 6 J kg -1, the net condensational energy input per day into a unit area atmospheric column is: d ( m H2 OL c )! 5"10 7 Jm #2 day #1 Assuming this heating is spread uniformly over the entire tropospheric depth, or over a unit area column of tropospheric air mass, implies a daily heating rate per unit air mass, i.e., (p s! p t )g!1, of: Q cnv = J d cnv = ( m L H2 O c)g " 5 Kday!1 c p c p (p s! p t ) As we ll see later, the distribution of condensational heating via tropical convection is nonuniformly distributed in the vertical and tends to maximize in the mid- to uppertroposphere (~300 to 400 mb), with heating rates of order 10 K day -1. Recalling that the radiative heating rate is ~ 1 K day -1, the above implies a vertical velocity scale ~ 10x larger than in nonprecipitating regions. This in turn implies a relatively large component of the divergent flow.