An Optimal Control Problem Formulation for. the Atmospheric Large-Scale Wave Dynamics
|
|
- Laureen O’Brien’
- 5 years ago
- Views:
Transcription
1 pplied Mathematical Sciences, Vol. 9, 5, no. 8, HIKRI Ltd, n Optimal Control Problem Formulation for the tmospheric Large-Scale Wave Dynamics Sergei Soldatenko Center for ustralian Weather and Climate Research 7 Collins Street, VIC 38, Melbourne, ustralia Rafael Yusupov St. Petersburg State Polytechnical University and SPIIRS 9 Polytechnical, St. Petersburg 955, Russia Copyright 5 Sergei Soldatenko and Rafael Yusupov. This is an open access article distributed under the Creative Commons ttribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. bstract In this paper, an optimal control problem is formulated for the large-scale atmospheric dynamics, which is mainly characterized by slow-moving planetary-scale waves with zonal wavenumbers to 6, also known as Rossby waves. n equivalent barotropic model of the atmosphere is used to describe the evolution of Rossby waves. Necessary conditions for optimality are obtained. Keywords: Optimal control, Rossby waves, distributed parameter system Introduction Weather modification is the man-made activity for intentionally manipulating or altering the weather and atmospheric processes. This activity has so far been considered outside the scope of control theory [, 3, 7]. For this reason the boundaries of weather modification, its goals and methods of achieving the objectives are commonly formulated in general terms and expected results are fairly vague. However, weather modification can be designed and implemented on the basis of control theory. This paper aims to formulate the optimal control problem for large-scale atmospheric dynamics, which is mainly characterized by slow-moving planetary-scale waves, also known as Rossby waves [4]. On weather maps these waves are observed as large-scale meanders of the pressure field in the
2 876 Sergei Soldatenko and Rafael Yusupov middle troposphere. Since Rossby waves are quasi-vertical and quasi-barotropic, in this study an equivalent barotropic model (EBM) of the atmosphere [5] is used as a tool to formulate the optimal control problem and to derive the necessary conditions for optimality. Equivalent barotropic model of the atmosphere tmospheric models, used for numerical weather prediction (NWP), are formulated as initial value problems for systems of nonlinear partial differential equations called primitive equations (PEs), which are comprised of the horizontal momentum equations, the equation of hydrostatic equilibrium, the thermodynamic equation and the continuity equation. We consider the inviscid adiabatic PEs written in the normalized isobaric coordinates ( xy,, ) [4] u u uu fv, () t P x v v uv fu, () t P y RT, u, (3) P where u, T u T S, (4) t uv is the horizontal velocity vector; is the two-dimensional gradient operator on a surface of constant ξ; is the geopotential; f is the Coriolis parameter; T is the temperature; dp dt is pressure vertical velocity, where p is pressure; R is the gas constant for dry air; S is the static stability measure, pp is the normalized pressure, P hpa is a standard pressure. The equation of the EBM is derived from the vorticity equation. Taking xof () and subtracting yof () gives the equation for the vertical component of vorticity v x u y, which represents the most interest in meteorology. Retaining the terms of order - s - yields the vorticity equation that is valid for large-scale atmospheric dynamics [4, 5]: f f f u v t x y P, (5) where we have used the continuity equation (4).
3 n optimal control problem formulation 877 The following basic assumptions are made in the derivation of EBM: the wind speeds vary little with height, and the variations of wind velocity with height are similar at all vertical levels. With these assumptions, at a certain vertical level the motion of the baroclinic atmosphere corresponds to the motion of a barotropic atmosphere and the vorticity equation can be written as [5]: where f f f u v, (6) t x y P is an empirical function used to describe the vertical variation of the horizontal wind speed, and is a pressure vertical velocity at the level ( z ). The level is known as the equivalent barotropic level, and the equation (6) is called as the equivalent barotropic vorticity equation. Climatological data shows that there exist two equivalent barotropic levels in the atmosphere: one around the 5 hpa (~5.5 km) level, and the other around hpa (~6 km) level. The mid-tropospheric level (5 hpa) is of particular interest for the NWP. Since the large-scale atmospheric motions are quasi-geostrophic [4], the vorticity equation (6) can be recast in terms of the geostrophic streamfunction f : f, f, (7) t P where a, b a xb y a yb x. By the definition [4, 5]: dp p p p u v gw. dt t x y Here w dz dt is the vertical velocity at the level z, is the air density and g is the gravity acceleration. The subscript refers to the level z. On the lower boundary the kinematic boundary condition requires that w. For the geostrophic motion it follows that up. Consequently, the value of reduces to p t. Substitution of this expression into (7) gives: f p, f. (8) t P t Suppose that p t and t are correlated with each other, where is the geopotential at the level. Since f p t R f t, then where R c is a correlation coefficient, and the equation (8) becomes Rc f, f. (9) t P t The coefficient Rc f P has the dimension of an inverse square length. Thus, from (9) we can obtain the following inhomogeneous Helmholtz equation: c
4 878 Sergei Soldatenko and Rafael Yusupov where q t, F, f q q L q F, () and L = Rc f P. We shall consider the equation () in a closed domain D of the xy - plane with a piecewise continuous boundary D. Let [, T ] be a time interval on which the solution of equation () is defined. Let the boundary curve be represented in parametric form: x, y, here and are piecewise continuous functions of the parameter. ssuming the coordinate origin coincides with the North Pole, we shall obtain: x r cos, y r sin, () where r, is a polar coordinate of a generic point on the circle. The right-hand side of () can be calculated if the geopotential field is given at the initial time x, y, x, y. Specifying const along the boundary D, t : and introducing the new dependent variable z t t, where t is the integration time step, we can obtain the following Dirichlet problem: where z L z G in D, () z on D, G F t. The problem () is solved numerically using an appropriate z t t and, therefore, the iteration method to obtain the function streamfunction: t z t Consequently, we can obtain the forecast of geopotential field f t, (3), which characterizes the Rossby wave dynamics. The forward difference (3) is used only on the first time step and then a central difference is applied. 3 Statement of the optimal control problem When the equation (6) is equal to zero on the right-hand side, we refer to it as a barotropic vorticity equation. Thus, (6) is a forced barotropic equation in which the forcing term is determined by the vertical motion at the lower boundary. The control objective is to manipulate the phase velocity of Rossby waves. It is known that the natural forcing of Rossby waves in the atmosphere is mainly of orographic or thermal origin. These two effects are formally described by the right-hand term in the vorticity equation (6), i.e. by the vertical velocity. Therefore, the variable can be chosen as a control variable. Let us consider the linearized vorticity equation on a periodic mid-latitude beta-plane [4]:
5 n optimal control problem formulation 879 ' ' ' u u, (4) t x L t x where f y is the latitudinal gradient of Coriolis parameter, u is a zonal unperturbed flow, u d u dy, ' is the perturbation of streamfunction. We seek a solution of the form ikxct ' x, y, t Re e cos ly. (5) Here is an amplitude of perturbation, k Lx and l Ly are the zonal and meridional wavenumbers respectively, and c is a perturbation phase speed. fter substituting (5) into (4) we can obtain the dispersion relationship for the Rossby waves: u k l u c k l L. (6) This formula shows that the atmospheric vertical structure through the factor L contributes to the dispersion properties of Rossby waves. However, the factor L depends on the correlation coefficient R, which in turn relates the vertical velocity c with the streamfunction ψ at the mid-tropospheric vertical level. Hence, the phase speed of Rossby waves is also a function of the vertical velocity. Let us now consider the following control system with the state variable, U U x, y : z z x y and the control variable z L z U G in D, z on D. (7) The term G is calculated using the geostrophic streamfunction at the initial time t : G G, f t, where f. Let us assume that the control variable U U, where U is the set of all permissible controls. The control U, from a physical standpoint, is a measure of additional vertical velocity on the bottom of the atmosphere near the ground. It is very important to note that the set U should be defined on the basis of physical and technical feasibility, taking into account the properties of the atmosphere as a physical object. Let us introduce the performance index as z z U dd, (8) D
6 88 Sergei Soldatenko and Rafael Yusupov where the pair zu, satisfies the equation (7) and z is the desired spatial distribution of z. The term U ( ) is proportional to the consumed energy. The optimal control problem is defined as follows: find the control U U generating the system state z X such that the performance index (8) is minimized. Here X is a set of additional constraints on the state variable z. The control problem, from a physical viewpoint, reflects the ability to manipulate phase speeds of large-scale atmospheric waves by changing a vertical velocity on the lower boundary of the atmosphere. Let us transform the equations (7) to a normal form by introducing new dependent and parametric variables [, 6]: z z z z z, z3,,, x y x y (9) z3 z3, U z G. x y L djoin the equations (9) to the performance index (8) with spatial-varying Lagrange multipliers,, 3,, and 3: z z z z z U z z3 x y x D z z3 z U z G dd. () y x L y The Hamiltonian function H associated with the optimal control problem is introduced as H z z U z 3 z3 3 U z G. L () The augmented performance index is rewritten then as: z z z3 z z z3 H 3 3 dd, x x x y y y D or in the compact form z z H dd, () x y D
7 n optimal control problem formulation 88 where z z, z, z 3,,, 3 and,, 3 subtract the term z x y. If we add and to the integrand of (), we shall obtain: H zdd z z dd. (3) D x y D x y By the use of Green's theorem, the second integral in (3) becomes a line integral, and then takes the following form: H zdd zd. D x y D The first variation in due to variations in the control and parametric variables is: H H z U dd zd, (4) D z x y U D where,, U U. We choose the functions and to cause the first term in parentheses in the double integral in (4) to vanish: H. (5) x y z Since the function z does not change along the boundary D ( z ), its variation z on D. For z and z 3, which are not defined on the boundary D, we require the relations to hold along D. Then, taking into account (), we can obtain the following boundary conditions for the system (9): cos sin, cos sin. (6) 3 3 The first variation in then becomes: H U dd. (7) U D For a minimum, it is necessary that.this can only be achieved if
8 88 Sergei Soldatenko and Rafael Yusupov H U. (8) Equations (5) and (8) with the boundary conditions (7) and (6) represent the necessary conditions for optimum, which can be written as 3 z z,, x y L x y (9) 3 3, U 3, 3, 3. x y Therefore, to find the control U that minimize the performance index (8) we must solve the system of partial differential equations (9) and (9) with unknown two-dimensional variables and with given boundary conditions: z,. (3) D D D Let us consider a particular case of the optimal control problem in which the constant is equals to zero. From (9) we can obtain that:, 3 and the necessary conditions (9) for an extremum of takes the form: z z,, x y y (3) 3, 3. x Since 3, the second and third equations of the system (3) are similar to the first two equations of the system (9) and we might make the following assumptions [, 6]: 3 z, z, 3 z, (3) where is an arbitrary constant. These assumptions are compatible with the boundary conditions (3). If we substitute (3) into the first equation of the system (3) then we get: from which z z z z, xy xy z z in the whole domain D. It follows that the performance index can achieve its absolute minimum under the constraint (7).
9 n optimal control problem formulation 883 To find the optimal control U that minimizes, we need to solve the following problem U z L z G, z D in the domain D. Rossby waves propagate westward relative to the mean quasi-zonal flow [4], so that it is possible for them to be stationary (with respect to the surface) in a westerly atmospheric flow. If the optimal control objective is to achieve the stationarity of Rossby waves by manipulating the lower boundary vertical velocity, i.e. z t, then the control U in the domain D can be find by t solving the equation U G. t first glance the results obtained seem almost obvious. However, in this paper we rigorously prove the existence of absolute extremum of the performance index (8) under the constraint (7). If various additional constraints are imposed on the control variable U, then the performance index (8) may not achieve its absolute extremum, and then the expression for calculating the control variable U is not so obvious. 4 Conclusion Based on the optimal control theory of distributed parameter systems, an optimal control problem for Rossby wave dynamics in the atmosphere has been theoretically formulated. The equivalent barotropic atmospheric model was used as an instrument to obtain the necessary conditions for optimality. The application of the optimal control theory to manipulate weather and climate (geoengineering) represents a new multidisciplinary research area. However, this problem is extremely complex due to the uniqueness of the control object, the atmosphere, which is very sophisticated to model and predict. This paper provides a basis for further research in the field of meteorological cybernetics that is a new research area of a self-regulating cybernetic system in which the atmosphere represents the control object and the human society plays the role of controller. There is no doubt that this problem requires also the consideration of physical, technical, ethical and legal aspects and limitations. cknowledgments. This work was supported in part by the Government of the Russian Federation under Grant 74-U. References [].-L. rmand, pplications of the theory of optimal control of distributed parameter systems to structural optimization, NS CR-44, Washington, D.C., 97, 54 p.
10 884 Sergei Soldatenko and Rafael Yusupov [] Critical Issues in Weather Modification Research. Committee on the Status and Future Directions in U.S Weather Modification Research and Operations, National Research Council. National Research Council. National cademic Press, Washington, D.C., 3, 44 p. [3]. R. Fleming, The pathological history of weather and climate modification: Three cycles of promise and hype, Historical Studies in the Physical and Biological Sciences, 37 (6), [4]. R. Holton, n introduction to dynamic meteorology, London, Elsevier, 4, 535 p. [5] I.. Kibel, n introduction to the hydro-dynamical methods of short period weather forecasting. Pergamon Press, New York, 963, 383 p. [6] K.. Lurie, pplied optimal control theory of distributed systems, New York, Springer, 993, 499 p. [7] D. L. Mitchell and W. Finnegan, Modification of Cirrus clouds to reduce global warming, Environmental Research Letters, 4 (9), 45, Received: anuary 6, 5; Published: anuary 9, 5
Chapter 9. Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation
Chapter 9 Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation 9.1 Geostrophy and scaling. We examined in the last chapter some consequences of the dynamical balances for low frequency, nearly
More information1/27/2010. With this method, all filed variables are separated into. from the basic state: Assumptions 1: : the basic state variables must
Lecture 5: Waves in Atmosphere Perturbation Method With this method, all filed variables are separated into two parts: (a) a basic state part and (b) a deviation from the basic state: Perturbation Method
More information2.5 Shallow water equations, quasigeostrophic filtering, and filtering of inertia-gravity waves
Chapter. The continuous equations φ=gh Φ=gH φ s =gh s Fig..5: Schematic of the shallow water model, a hydrostatic, incompressible fluid with a rigid bottom h s (x,y), a free surface h(x,y,t), and horizontal
More informationEliassen-Palm Theory
Eliassen-Palm Theory David Painemal MPO611 April 2007 I. Introduction The separation of the flow into its zonal average and the deviations therefrom has been a dominant paradigm for analyses of the general
More informationSt. Petersburg Institute for Informatics and Automation, The Russian Academy of Sciences, No. 39, 14th Line, St. Petersburg , Russia
Journal of Control Science and Engineering Volume 216, Article ID 1547462, 9 pages http://d.doi.org/1.1155/216/1547462 Research Article The Determination of Feasible Control Variables for Geoengineering
More informationIn two-dimensional barotropic flow, there is an exact relationship between mass
19. Baroclinic Instability In two-dimensional barotropic flow, there is an exact relationship between mass streamfunction ψ and the conserved quantity, vorticity (η) given by η = 2 ψ.the evolution of the
More informationESCI 343 Atmospheric Dynamics II Lesson 11 - Rossby Waves
ESCI 343 Atmospheric Dynamics II Lesson 11 - Rossby Waves Reference: An Introduction to Dynamic Meteorology (4 rd edition), J.R. Holton Atmosphere-Ocean Dynamics, A.E. Gill Fundamentals of Atmospheric
More informationCirculation and Vorticity
Circulation and Vorticity Example: Rotation in the atmosphere water vapor satellite animation Circulation a macroscopic measure of rotation for a finite area of a fluid Vorticity a microscopic measure
More information( ) (9.1.1) Chapter 9. Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation. 9.1 Geostrophy and scaling.
Chapter 9 Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation 9.1 Geostrophy and scaling. We examined in the last chapter some consequences of the dynamical balances for low frequency, nearly
More informationChapter 2. Quasi-Geostrophic Theory: Formulation (review) ε =U f o L <<1, β = 2Ω cosθ o R. 2.1 Introduction
Chapter 2. Quasi-Geostrophic Theory: Formulation (review) 2.1 Introduction For most of the course we will be concerned with instabilities that an be analyzed by the quasi-geostrophic equations. These are
More informationThe dynamics of a simple troposphere-stratosphere model
The dynamics of a simple troposphere-stratosphere model by Jessica Duin under supervision of Prof. Dr. A. Doelman, UvA Dr. W.T.M. Verkley, KNMI August 31, 25 Universiteit van Amsterdam Korteweg-de vries
More informationDynamics and Kinematics
Geophysics Fluid Dynamics () Syllabus Course Time Lectures: Tu, Th 09:30-10:50 Discussion: 3315 Croul Hall Text Book J. R. Holton, "An introduction to Dynamic Meteorology", Academic Press (Ch. 1, 2, 3,
More informationGeophysics Fluid Dynamics (ESS228)
Geophysics Fluid Dynamics (ESS228) Course Time Lectures: Tu, Th 09:30-10:50 Discussion: 3315 Croul Hall Text Book J. R. Holton, "An introduction to Dynamic Meteorology", Academic Press (Ch. 1, 2, 3, 4,
More informationThe Eady problem of baroclinic instability described in section 19a was shown to
0. The Charney-Stern Theorem The Eady problem of baroclinic instability described in section 19a was shown to be remarkably similar to the Rayleigh instability of barotropic flow described in Chapter 18.
More informationLecture #2 Planetary Wave Models. Charles McLandress (Banff Summer School 7-13 May 2005)
Lecture #2 Planetary Wave Models Charles McLandress (Banff Summer School 7-13 May 2005) 1 Outline of Lecture 1. Observational motivation 2. Forced planetary waves in the stratosphere 3. Traveling planetary
More informationGravity Waves. Lecture 5: Waves in Atmosphere. Waves in the Atmosphere and Oceans. Internal Gravity (Buoyancy) Waves 2/9/2017
Lecture 5: Waves in Atmosphere Perturbation Method Properties of Wave Shallow Water Model Gravity Waves Rossby Waves Waves in the Atmosphere and Oceans Restoring Force Conservation of potential temperature
More informationA few examples Shallow water equation derivation and solutions. Goal: Develop the mathematical foundation of tropical waves
A few examples Shallow water equation derivation and solutions Goal: Develop the mathematical foundation of tropical waves Previously: MCS Hovmoller Propagating w/ wave velocity From Chen et al (1996)
More informationρ 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
Fundamental Equation in z coordinate p = ρrt or pα = RT Du uv tanφ + uw Dt a a = 1 p ρ x + fv f 'w + F x Dv Dt + u2 tanφ + vw a a = 1 p ρ y fu + F y Dw Dt u2 + v 2 = 1 p a ρ z g + f 'u + F z Dρ Dt + ρ
More informationDynamics Rotating Tank
Institute for Atmospheric and Climate Science - IACETH Atmospheric Physics Lab Work Dynamics Rotating Tank Large scale flows on different latitudes of the rotating Earth Abstract The large scale atmospheric
More informationOn the effect of forward shear and reversed shear baroclinic flows for polar low developments. Thor Erik Nordeng Norwegian Meteorological Institute
On the effect of forward shear and reversed shear baroclinic flows for polar low developments Thor Erik Nordeng Norwegian Meteorological Institute Outline Baroclinic growth a) Normal mode solution b) Initial
More informationt tendency advection convergence twisting baroclinicity
RELATIVE VORTICITY EQUATION Newton s law in a rotating frame in z-coordinate (frictionless): U + U U = 2Ω U Φ α p U + U U 2 + ( U) U = 2Ω U Φ α p Applying to both sides, and noting ω U and using identities
More informationChapter 2. The continuous equations
Chapter. The continuous equations Fig. 1.: Schematic of a forecast with slowly varying weather-related variations and superimposed high frequency Lamb waves. Note that even though the forecast of the slow
More information2. Baroclinic Instability and Midlatitude Dynamics
2. Baroclinic Instability and Midlatitude Dynamics Midlatitude Jet Stream Climatology (Atlantic and Pacific) Copyright 26 Emily Shuckburgh, University of Cambridge. Not to be quoted or reproduced without
More information1/18/2011. Conservation of Momentum Conservation of Mass Conservation of Energy Scaling Analysis ESS227 Prof. Jin-Yi Yu
Lecture 2: Basic Conservation Laws Conservation Law of Momentum Newton s 2 nd Law of Momentum = absolute velocity viewed in an inertial system = rate of change of Ua following the motion in an inertial
More informationChapter 3. Stability theory for zonal flows :formulation
Chapter 3. Stability theory for zonal flows :formulation 3.1 Introduction Although flows in the atmosphere and ocean are never strictly zonal major currents are nearly so and the simplifications springing
More informationQuasi-geostrophic ocean models
Quasi-geostrophic ocean models March 19, 2002 1 Introduction The starting point for theoretical and numerical study of the three dimensional large-scale circulation of the atmosphere and ocean is a vorticity
More informationThe General Circulation of the Atmosphere: A Numerical Experiment
The General Circulation of the Atmosphere: A Numerical Experiment Norman A. Phillips (1956) Presentation by Lukas Strebel and Fabian Thüring Goal of the Model Numerically predict the mean state of the
More information4. Atmospheric transport. Daniel J. Jacob, Atmospheric Chemistry, Harvard University, Spring 2017
4. Atmospheric transport Daniel J. Jacob, Atmospheric Chemistry, Harvard University, Spring 2017 Forces in the atmosphere: Gravity g Pressure-gradient ap = ( 1/ ρ ) dp / dx for x-direction (also y, z directions)
More informationNonlinear baroclinic dynamics of surface cyclones crossing a zonal jet
Nonlinear baroclinic dynamics of surface cyclones crossing a zonal jet Jean-Baptiste GILET, Matthieu Plu and Gwendal Rivière CNRM/GAME (Météo-France, CNRS) 3rd THORPEX International Science Symposium Monterey,
More informationDust devils, water spouts, tornados
Balanced flow Things we know Primitive equations are very comprehensive, but there may be a number of vast simplifications that may be relevant (e.g., geostrophic balance). Seems that there are things
More informationConservation of Mass Conservation of Energy Scaling Analysis. ESS227 Prof. Jin-Yi Yu
Lecture 2: Basic Conservation Laws Conservation of Momentum Conservation of Mass Conservation of Energy Scaling Analysis Conservation Law of Momentum Newton s 2 nd Law of Momentum = absolute velocity viewed
More information[1]{Izaña Atmospheric Research Centre (AEMET), Santa Cruz de Tenerife, Spain}
Supplement of Pivotal role of the North African Dipole Intensity (NAFDI) on alternate Saharan dust export over the North Atlantic and the Mediterranean, and relationship with the Saharan Heat Low and mid-latitude
More informationMeasurement of Rotation. Circulation. Example. Lecture 4: Circulation and Vorticity 1/31/2017
Lecture 4: Circulation and Vorticity Measurement of Rotation Circulation Bjerknes Circulation Theorem Vorticity Potential Vorticity Conservation of Potential Vorticity Circulation and vorticity are the
More information8 Mechanisms for tropical rainfall responses to equatorial
8 Mechanisms for tropical rainfall responses to equatorial heating More reading: 1. Hamouda, M. and Kucharski, F. (2019) Ekman pumping Mechanism driving Precipitation anomalies in Response to Equatorial
More informationNWP Equations (Adapted from UCAR/COMET Online Modules)
NWP Equations (Adapted from UCAR/COMET Online Modules) Certain physical laws of motion and conservation of energy (for example, Newton's Second Law of Motion and the First Law of Thermodynamics) govern
More informationDynamics of the Extratropical Response to Tropical Heating
Regional and Local Climate Modeling and Analysis Research Group R e L o C l i m Dynamics of the Extratropical Response to Tropical Heating (1) Wegener Center for Climate and Global Change (WegCenter) and
More informationNonlinear Balance on an Equatorial Beta Plane
Nonlinear Balance on an Equatorial Beta Plane David J. Raymond Physics Department and Geophysical Research Center New Mexico Tech Socorro, NM 87801 April 26, 2009 Summary Extension of the nonlinear balance
More information9 Rossby Waves. 9.1 Non-divergent barotropic vorticity equation. CSU ATS601 Fall (Holton Chapter 7, Vallis Chapter 5)
9 Rossby Waves (Holton Chapter 7, Vallis Chapter 5) 9.1 Non-divergent barotropic vorticity equation We are now at a point that we can discuss our first fundamental application of the equations of motion:
More information2. Outline of the MRI-EPS
2. Outline of the MRI-EPS The MRI-EPS includes BGM cycle system running on the MRI supercomputer system, which is developed by using the operational one-month forecasting system by the Climate Prediction
More informationThe Shallow Water Equations
If you have not already done so, you are strongly encouraged to read the companion file on the non-divergent barotropic vorticity equation, before proceeding to this shallow water case. We do not repeat
More informationHEIGHT-LATITUDE STRUCTURE OF PLANETARY WAVES IN THE STRATOSPHERE AND TROPOSPHERE. V. Guryanov, A. Fahrutdinova, S. Yurtaeva
HEIGHT-LATITUDE STRUCTURE OF PLANETARY WAVES IN THE STRATOSPHERE AND TROPOSPHERE INTRODUCTION V. Guryanov, A. Fahrutdinova, S. Yurtaeva Kazan State University, Kazan, Russia When constructing empirical
More informationAtmospheric Fronts. The material in this section is based largely on. Lectures on Dynamical Meteorology by Roger Smith.
Atmospheric Fronts The material in this section is based largely on Lectures on Dynamical Meteorology by Roger Smith. Atmospheric Fronts 2 Atmospheric Fronts A front is the sloping interfacial region of
More informationFundamentals of Weather and Climate
Fundamentals of Weather and Climate ROBIN McILVEEN Environmental Science Division Institute of Environmental and Biological Sciences Lancaster University CHAPMAN & HALL London Glasgow Weinheim New York
More informationDynamic Meteorology - Introduction
Dynamic Meteorology - Introduction Atmospheric dynamics the study of atmospheric motions that are associated with weather and climate We will consider the atmosphere to be a continuous fluid medium, or
More informationEATS Notes 1. Some course material will be online at
EATS 3040-2015 Notes 1 14 Aug 2015 Some course material will be online at http://www.yorku.ca/pat/esse3040/ HH = Holton and Hakim. An Introduction to Dynamic Meteorology, 5th Edition. Most of the images
More informationGoal: Use understanding of physically-relevant scales to reduce the complexity of the governing equations
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
More informationEART164: PLANETARY ATMOSPHERES
EART164: PLANETARY ATMOSPHERES Francis Nimmo Last Week Radiative Transfer Black body radiation, Planck function, Wien s law Absorption, emission, opacity, optical depth Intensity, flux Radiative diffusion,
More informationIntroduction to Isentropic Coordinates: a new view of mean meridional & eddy circulations. Cristiana Stan
Introduction to Isentropic Coordinates: a new view of mean meridional & eddy circulations Cristiana Stan School and Conference on the General Circulation of the Atmosphere and Oceans: a Modern Perspective
More information7 The General Circulation
7 The General Circulation 7.1 The axisymmetric state At the beginning of the class, we discussed the nonlinear, inviscid, axisymmetric theory of the meridional structure of the atmosphere. The important
More informationCan a Simple Two-Layer Model Capture the Structure of Easterly Waves?
Can a Simple Two-Layer Model Capture the Structure of Easterly Waves? Cheryl L. Lacotta 1 Introduction Most tropical storms in the Atlantic, and even many in the eastern Pacific, are due to disturbances
More informationIntroduction to Isentropic Coordinates:! a new view of mean meridional & eddy circulations" Cristiana Stan
Introduction to Isentropic Coordinates:! a new view of mean meridional & eddy circulations" Cristiana Stan School and Conference on the General Circulation of the Atmosphere and Oceans: a Modern Perspective!
More information6 Two-layer shallow water theory.
6 Two-layer shallow water theory. Wewillnowgoontolookatashallowwatersystemthathastwolayersofdifferent density. This is the next level of complexity and a simple starting point for understanding the behaviour
More informationChapter 13 Instability on non-parallel flow Introduction and formulation
Chapter 13 Instability on non-parallel flow. 13.1 Introduction and formulation We have concentrated our discussion on the instabilities of parallel, zonal flows. There is the largest amount of literature
More informationNote that Rossby waves are tranverse waves, that is the particles move perpendicular to the direction of propagation. f up, down (clockwise)
Ocean 423 Rossby waves 1 Rossby waves: Restoring force is the north-south gradient of background potential vorticity (f/h). That gradient can be due to either the variation in f with latitude, or to a
More informationQuasi-geostrophic system
Quasi-eostrophic system (or, why we love elliptic equations for QGPV) Charney s QG the motion of lare-scale atmospheric disturbances is overned by Laws of conservation of potential temperature and potential
More information) 2 ψ +β ψ. x = 0. (71) ν = uk βk/k 2, (74) c x u = β/k 2. (75)
3 Rossby Waves 3.1 Free Barotropic Rossby Waves The dispersion relation for free barotropic Rossby waves can be derived by linearizing the barotropic vortiticy equation in the form (21). This equation
More informationQuasi-equilibrium Theory of Small Perturbations to Radiative- Convective Equilibrium States
Quasi-equilibrium Theory of Small Perturbations to Radiative- Convective Equilibrium States See CalTech 2005 paper on course web site Free troposphere assumed to have moist adiabatic lapse rate (s* does
More information+ ω = 0, (1) (b) In geometric height coordinates in the rotating frame of the Earth, momentum balance for an inviscid fluid is given by
Problem Sheet 1: Due Thurs 3rd Feb 1. Primitive equations in different coordinate systems (a) Using Lagrangian considerations and starting from an infinitesimal mass element in cartesian coordinates (x,y,z)
More informationChapter 7: Circulation and Vorticity
Chapter 7: Circulation and Vorticity Circulation C = u ds Integration is performed in a counterclockwise direction C is positive for counterclockwise flow!!! Kelvin s Circulation Theorem The rate of change
More informationAtmospheric dynamics and meteorology
Atmospheric dynamics and meteorology B. Legras, http://www.lmd.ens.fr/legras III Frontogenesis (pre requisite: quasi-geostrophic equation, baroclinic instability in the Eady and Phillips models ) Recommended
More informationweak mean flow R. M. Samelson
An effective-β vector for linear planetary waves on a weak mean flow R. M. Samelson College of Oceanic and Atmospheric Sciences 14 COAS Admin Bldg Oregon State University Corvallis, OR 97331-553 USA rsamelson@coas.oregonstate.edu
More informationA mechanistic model study of quasi-stationary wave reflection. D.A. Ortland T.J. Dunkerton NorthWest Research Associates Bellevue WA
A mechanistic model study of quasi-stationary wave reflection D.A. Ortland T.J. Dunkerton ortland@nwra.com NorthWest Research Associates Bellevue WA Quasi-stationary flow Describe in terms of zonal mean
More information7 Balanced Motion. 7.1 Return of the...scale analysis for hydrostatic balance! CSU ATS601 Fall 2015
7 Balanced Motion We previously discussed the concept of balance earlier, in the context of hydrostatic balance. Recall that the balanced condition means no accelerations (balance of forces). That is,
More informationLecture 12: Angular Momentum and the Hadley Circulation
Lecture 12: Angular Momentum and the Hadley Circulation September 30, 2003 We learnt last time that there is a planetary radiative drive net warming in the tropics, cooling over the pole which induces
More informationAtmosphere, Ocean and Climate Dynamics Answers to Chapter 8
Atmosphere, Ocean and Climate Dynamics Answers to Chapter 8 1. Consider a zonally symmetric circulation (i.e., one with no longitudinal variations) in the atmosphere. In the inviscid upper troposphere,
More informationPlanetary Atmospheres. Structure Composition Clouds Photochemistry Meteorology Atmospheric Escape
Planetary Atmospheres Structure Composition Clouds Photochemistry Meteorology Atmospheric Escape Photochemistry We can characterize chemical reactions in the atmosphere in the following way: 1. Photolysis:
More informationInfluence of Gravity Waves on the Atmospheric Climate
Influence of Gravity Waves on the Atmospheric Climate François Lott, LMD/CNRS, Ecole Normale Supérieure, Paris flott@lmd.ens.fr 1)Dynamical impact of mountains on atmospheric flows 3)Non-orographic gravity
More informationFour ways of inferring the MMC. 1. direct measurement of [v] 2. vorticity balance. 3. total energy balance
Four ways of inferring the MMC 1. direct measurement of [v] 2. vorticity balance 3. total energy balance 4. eliminating time derivatives in governing equations Four ways of inferring the MMC 1. direct
More informationThe Euler Equation of Gas-Dynamics
The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u
More informationwhere p oo is a reference level constant pressure (often 10 5 Pa). Since θ is conserved for adiabatic motions, a prognostic temperature equation is:
1 Appendix C Useful Equations Purposes: Provide foundation equations and sketch some derivations. These equations are used as starting places for discussions in various parts of the book. C.1. Thermodynamic
More informationTraveling planetary-scale Rossby waves in the winter stratosphere: The role of tropospheric baroclinic instability
GEOPHYSICAL RESEARCH LETTERS, VOL.???, XXXX, DOI:.29/, 1 2 Traveling planetary-scale Rossby waves in the winter stratosphere: The role of tropospheric baroclinic instability Daniela I.V. Domeisen, 1 R.
More informationAFRICAN EASTERLY WAVES IN CURRENT AND FUTURE CLIMATES
AFRICAN EASTERLY WAVES IN CURRENT AND FUTURE CLIMATES Victoria Dollar RTG Seminar Research - Spring 2018 April 16, 2018 Victoria Dollar ASU April 16, 2018 1 / 26 Overview Introduction Rossby waves and
More information3. Midlatitude Storm Tracks and the North Atlantic Oscillation
3. Midlatitude Storm Tracks and the North Atlantic Oscillation Copyright 2006 Emily Shuckburgh, University of Cambridge. Not to be quoted or reproduced without permission. EFS 3/1 Review of key results
More informationROSSBY WAVE PROPAGATION
ROSSBY WAVE PROPAGATION (PHH lecture 4) The presence of a gradient of PV (or q.-g. p.v.) allows slow wave motions generally called Rossby waves These waves arise through the Rossby restoration mechanism,
More informationExamples of Pressure Gradient. Pressure Gradient Force. Chapter 7: Forces and Force Balances. Forces that Affect Atmospheric Motion 2/2/2015
Chapter 7: Forces and Force Balances Forces that Affect Atmospheric Motion Fundamental force - Apparent force - Pressure gradient force Gravitational force Frictional force Centrifugal force Forces that
More informationLecture 1 ATS 601. Thomas Birner, CSU. ATS 601 Lecture 1
Lecture 1 ATS 601 Thomas Birner, CSU About your Instructor: Thomas Birner Assistant Professor, joined CSU 10/2008 M.Sc. Physics (Condensed Matter Theory), U of Leipzig, Germany Ph.D. Atmospheric Science
More informationPart-8c Circulation (Cont)
Part-8c Circulation (Cont) Global Circulation Means of Transfering Heat Easterlies /Westerlies Polar Front Planetary Waves Gravity Waves Mars Circulation Giant Planet Atmospheres Zones and Belts Global
More informationThe dynamics of high and low pressure systems
The dynamics of high and low pressure systems Newton s second law for a parcel of air in an inertial coordinate system (a coordinate system in which the coordinate axes do not change direction and are
More informationATMOSPHERIC MOTION I (ATM S 441/503 )
http://earth.nullschool.net/ ATMOSPHERIC MOTION I (ATM S 441/503 ) INSTRUCTOR Daehyun Kim Born in 1980 B.S. 2003 Ph.D. 2010 2010-2013 2014- Assistant Professor at Dept. of Atmospheric Sciences Office:
More informationMT Dynamic Meteorology II MWF Spring :00 a.m. - 8:50 a.m. Dr. Jim Koermer Boyd 306 DYNAMIC METEOROLOGY II SYLLABUS
MT5320.01 - Dynamic Meteorology II MWF Spring 2011 8:00 a.m. - 8:50 a.m. Dr. Jim Koermer Boyd 306 DYNAMIC METEOROLOGY II SYLLABUS COURSE DESCRIPTION: Dynamic Meteorology II involves the application of
More informationThe stratospheric response to extratropical torques and its relationship with the annular mode
The stratospheric response to extratropical torques and its relationship with the annular mode Peter Watson 1, Lesley Gray 1,2 1. Atmospheric, Oceanic and Planetary Physics, Oxford University 2. National
More informationCHAPTER 4. THE HADLEY CIRCULATION 59 smaller than that in midlatitudes. This is illustrated in Fig. 4.2 which shows the departures from zonal symmetry
Chapter 4 THE HADLEY CIRCULATION The early work on the mean meridional circulation of the tropics was motivated by observations of the trade winds. Halley (1686) and Hadley (1735) concluded that the trade
More informationBaroclinic Rossby waves in the ocean: normal modes, phase speeds and instability
Baroclinic Rossby waves in the ocean: normal modes, phase speeds and instability J. H. LaCasce, University of Oslo J. Pedlosky, Woods Hole Oceanographic Institution P. E. Isachsen, Norwegian Meteorological
More informationModelling the atmosphere. Hennie Kelder University of Technology Eindhoven
Modelling the atmosphere Hennie Kelder University of Technology Eindhoven Content Basics of the atmosphere Atmospheric dynamics Large scale circulation Planetary waves Brewer-Dobson circulation Some Key
More informationEliassen-Palm Cross Sections Edmon et al. (1980)
Eliassen-Palm Cross Sections Edmon et al. (1980) Cecily Keppel November 14 2014 Eliassen-Palm Flux For β-plane Coordinates (y, p) in northward, vertical directions Zonal means F = v u f (y) v θ θ p F will
More information196 7 atmospheric oscillations:
196 7 atmospheric oscillations: 7.4 INTERNAL GRAVITY (BUOYANCY) WAVES We now consider the nature of gravity wave propagation in the atmosphere. Atmospheric gravity waves can only exist when the atmosphere
More informationLinear model for investigation of nonlinear NWP model accuracy. Marko Zirk, University of Tartu
Linear model for investigation of nonlinear NWP model accuracy Marko Zirk, University of Tartu Introduction A method for finding numerical solution of non-hydrostatic linear equations of atmospheric dynamics
More informationUsing simplified vorticity equation,* by assumption 1 above: *Metr 430 handout on Circulation and Vorticity. Equations (4) and (5) on that handout
Rossby Wave Equation A. Assumptions 1. Non-divergence 2. Initially, zonal flow, or nearly zonal flow in which u>>v>>w. 3. Initial westerly wind is geostrophic and does not vary along the x-axis and equations
More information1/25/2010. Circulation and vorticity are the two primary
Lecture 4: Circulation and Vorticity Measurement of Rotation Circulation Bjerknes Circulation Theorem Vorticity Potential Vorticity Conservation of Potential Vorticity Circulation and vorticity are the
More informationFrictional Damping of Baroclinic Waves
Frictional Damping of Baroclinic Waves HHH, 26th May 2006 Bob Plant With thanks to Stephen Belcher, Brian Hoskins, Dan Adamson Motivation Control simulation, T+60 Simulation with no boundary layer turbulence,
More informationDynamics of the Atmosphere. Large-scale flow with rotation and stratification
12.810 Dynamics of the Atmosphere Large-scale flow with rotation and stratification Visualization of meandering jet stream Upper level winds from June 10th to July 8th 1988 from MERRA Red shows faster
More information1/18/2011. From the hydrostatic equation, it is clear that a single. pressure and height in each vertical column of the atmosphere.
Lecture 3: Applications of Basic Equations Pressure as Vertical Coordinate From the hydrostatic equation, it is clear that a single valued monotonic relationship exists between pressure and height in each
More informationAtmospheric Dynamics: lecture 11
Atmospheric Dynamics: lecture 11 (http://www.staff.science.uu.nl/~delde102/) Chapter 9 Baroclinic waves and cyclogenesis What is a baroclinic wave? Quasi-geostrophic equations Omega equation Original articles:
More informationReynolds Averaging. We separate the dynamical fields into slowly varying mean fields and rapidly varying turbulent components.
Reynolds Averaging Reynolds Averaging We separate the dynamical fields into sloly varying mean fields and rapidly varying turbulent components. Reynolds Averaging We separate the dynamical fields into
More informationMechanisms for influence of the stratosphere on the troposphere
Mechanisms for influence of the stratosphere on the troposphere Radiative Mass transfer/chemical Dynamical Alan Plumb M. I. T. Apr 2003 What dynamical effect does the stratosphere have on the troposphere?
More informationwarmest (coldest) temperatures at summer heat dispersed upward by vertical motion Prof. Jin-Yi Yu ESS200A heated by solar radiation at the base
Pole Eq Lecture 3: ATMOSPHERE (Outline) JS JP Hadley Cell Ferrel Cell Polar Cell (driven by eddies) L H L H Basic Structures and Dynamics General Circulation in the Troposphere General Circulation in the
More informationUnderstanding inertial instability on the f-plane with complete Coriolis force
Understanding inertial instability on the f-plane with complete Coriolis force Abstract Vladimir Zeitlin Laboratory of Dynamical Meteorology, University P. and M. Curie and Ecole Normale Supérieure, Paris,
More informationInfluence of Gravity Waves on the Atmospheric Climate
Influence of Gravity Waves on the Atmospheric Climate François Lott, LMD/CNRS, Ecole Normale Supérieure, Paris flott@lmd.ens.fr 1)Dynamical impact of mountains on atmospheric flows 2)Representations of
More informationQuick Recapitulation of Fluid Mechanics
Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change
More informationChapter 9. Barotropic Instability. 9.1 Linearized governing equations
Chapter 9 Barotropic Instability The ossby wave is the building block of low ossby number geophysical fluid dynamics. In this chapter we learn how ossby waves can interact with each other to produce a
More information