OZO V.1.0: Software for Solving a Generalised Omega Equation and the Zwack–Okossi Height Tendency Equation Using WRF Model Output

Total Page:16

File Type:pdf, Size:1020Kb

OZO V.1.0: Software for Solving a Generalised Omega Equation and the Zwack–Okossi Height Tendency Equation Using WRF Model Output Geosci. Model Dev., 10, 827–841, 2017 www.geosci-model-dev.net/10/827/2017/ doi:10.5194/gmd-10-827-2017 © Author(s) 2017. CC Attribution 3.0 License. OZO v.1.0: software for solving a generalised omega equation and the Zwack–Okossi height tendency equation using WRF model output Mika Rantanen1, Jouni Räisänen1, Juha Lento2, Oleg Stepanyuk1, Olle Räty1, Victoria A. Sinclair1, and Heikki Järvinen1 1Department of Physics, University Of Helsinki, Helsinki, Finland 2CSC – IT Center for Science, Espoo, Finland Correspondence to: Mika Rantanen (mika.p.rantanen@helsinki.fi) Received: 19 August 2016 – Discussion started: 29 September 2016 Revised: 3 February 2017 – Accepted: 6 February 2017 – Published: 21 February 2017 Abstract. A software package (OZO, Omega–Zwack– ucation, and the distribution of the software will be supported Okossi) was developed to diagnose the processes that af- by the authors. fect vertical motions and geopotential height tendencies in weather systems simulated by the Weather Research and Forecasting (WRF) model. First, this software solves a gen- eralised omega equation to calculate the vertical motions as- 1 Introduction sociated with different physical forcings: vorticity advection, thermal advection, friction, diabatic heating, and an imbal- Today, high-resolution atmospheric reanalyses provide a ance term between vorticity and temperature tendencies. Af- three-dimensional (3-D) view on the evolution of synoptic- ter this, the corresponding height tendencies are calculated scale weather systems (Dee et al., 2011; Rienecker et al., with the Zwack–Okossi tendency equation. The resulting 2011). On the other hand, simulations by atmospheric models height tendency components thus contain both the direct ef- allow for exploring the sensitivity of both real-world and ide- fect from the forcing itself and the indirect effects (related alised weather systems to factors such as the initial state (e.g. to the vertical motion induced by the same forcing) of each Leutbecher et al., 2002; Hoskins and Coutinho, 2005), lower physical mechanism. This approach has an advantage com- boundary conditions (e.g. Elguindi et al., 2005; Hirata et al., pared with previous studies with the Zwack–Okossi equa- 2016), and representation of sub-grid scale processes (e.g. tion, in which vertical motions were used as an independent Wernli et al., 2002; Liu et al., 2004; Beare, 2007). Neverthe- forcing but were typically found to compensate the effects of less, the complexity of atmospheric dynamics often makes other forcings. the physical interpretation of reanalysis data and model out- The software is currently tailored to use input from WRF put far from simple. Therefore, there is also a need for diag- simulations with Cartesian geometry. As an illustration, re- nostic methods that help to separate the effects of individual sults for an idealised 10-day baroclinic wave simulation are dynamical and physical processes on the structure and evo- presented. An excellent agreement is found between OZO lution of weather systems. and the direct WRF output for both the vertical motion and Two variables that are of special interest in the study of the height tendency fields. The individual vertical motion and synoptic-scale weather systems are the geopotential height height tendency components demonstrate the importance of tendency and vertical motion (Holton and Hakim, 2012). both adiabatic and diabatic processes for the simulated cy- Height tendencies are directly related to the movement clone. OZO is an open-source tool for both research and ed- and intensification or decay of low- and high-pressure sys- tems. Vertical motions affect atmospheric humidity, cloudi- ness, and precipitation. They also play a crucial role in Published by Copernicus Publications on behalf of the European Geosciences Union. 828 M. Rantanen et al.: OZO software atmospheric dynamics by inducing adiabatic temperature bining the QG vorticity and thermodynamic equations, infers changes, by generating cyclonic or anticyclonic vorticity, and vertical motion from geostrophic advection of absolute vor- by converting available potential energy to kinetic energy ticity and temperature (Holton, 1992): (Lorenz, 1955; Holton and Hakim, 2012). For the need of diagnostic tools, some software pack- LQG.!/ D FV (QG) C FT (QG); (1) ages have been developed to separate the contributions of each forcing to the vertical motion and height tendency. where DIONYSOS (Caron et al., 2006), a tool for analysing numer- @2! ically simulated weather systems, provides currently online L .!/ D σ .p/r2! C f 2 (2) QG 0 @p2 daily diagnostics for the output of numerical weather predic- tion models. RIP4 (Stoelinga, 2009) can calculate Q-vectors and the two right-hand side (RHS) terms are (Hoskins et al., 1978; Holton and Hakim, 2012) and a quasi- geostrophic (QG) vertical motion (Holton and Hakim, 2012) @ F D f V · r ζ C f ; (3) from Weather Forecast and Research (WRF) model output, V (QG) @p g g but the division of ! into contributions from various atmo- R F D r2 V · rT : (4) spheric processes is not possible in RIP4. Furthermore, many T (QG) p g research groups have developed tools for their own needs but do not have resources to distribute the software. (the notation is conventional, see Table1 for an explanation Here, we introduce a software package Omega–Zwack– of the symbols.) Okossi (OZO) that can be used for diagnosing the contribu- Qualitatively, the QG omega equation indicates that cy- tions of different dynamical and physical processes to atmo- clonic (anticyclonic) vorticity advection increasing with spheric vertical motions and height tendencies. OZO calcu- height and a maximum of warm (cold) advection should in- lates vertical motion from a quasi-geostrophic and a gener- duce rising (sinking) motion in the atmosphere. However, alised omega equation (Räisänen, 1995) while height tenden- when deriving this equation, ageostrophic winds, diabatic cies are calculated using the Zwack–Okossi tendency equa- heating, and friction are neglected. In addition, hydrostatic tion (Zwack and Okossi, 1986; Lupo and Smith, 1998). The stability is treated as a constant and several terms in the vor- current, first version of OZO has been tailored to use out- ticity equation are omitted. Although Eq. (1) often provides a put from the WRF model (Wang et al., 2007; Shamarock reasonable estimate of synoptic-scale vertical motions at ex- et al., 2008) simulations run with idealised Cartesian geome- tratropical latitudes (Räisänen, 1995), these approximations try. Due to the wide use of the WRF model, we expect OZO inevitably deteriorate the accuracy of the QG omega equation to be a useful open-source tool for both research and educa- solution. tion. The omega equation can be generalised by relaxing the QG In the following, we first introduce the equations solved approximations (e.g. Krishnamurti, 1968; Pauley and Nie- by OZO: the two forms of the omega equation in Sect. 2.1 man, 1992; Räisänen, 1995). Here we use the formulation and the Zwack–Okossi height tendency equation in Sect. 2.2. The numerical techniques used in solving these equations L.!/ D FV C FT C FV C FQ C FA; (5) are described in Sect. 3. We have tested the software us- ing output from an idealised 25 km resolution WRF simu- where lation described in Sect. 4. Section 5 provides some compu- @2! @2ζ tational aspects of the software. The next two sections give L.!/ Dr2(σ !/ C f (ζ C f / − f ! @p2 @p2 an overview of the vertical motion (Sect. 6) and height ten- @ @V dency (Sect. 7) calculations for this simulation. Software lim- C f k · × r! (6) itations and plans for future development are presented in @p @p Sect. 8, and the conclusions are given in Sect. 9. Finally, information about the data and code availability is given in is a generalised form of Eq. (2) and the RHS terms have the Sect. 10. expressions @ F D f V · r .ζ C f /; (7) V @p 2 Equations R 2 FT D r .V · rT /; (8) 2.1 Omega equation p @ FF D −f [k · .r × F /]; (9) The omega equation is a diagnostic tool for estimating atmo- @p spheric vertical motions and studying their physical and dy- R 2 FQ D − r Q; (10) namical causes. Its well-known QG form, obtained by com- cpp Geosci. Model Dev., 10, 827–841, 2017 www.geosci-model-dev.net/10/827/2017/ M. Rantanen et al.: OZO software 829 Table 1. List of mathematical symbols. −1 cp = 1004 J kg specific heat of dry air at constant volume f Coriolis parameter F forcing in the omega equation F friction force per unit mass g = 9.81 m s−2 gravitational acceleration k unit vector along the vertical axis L linear operator in the left-hand side of the omega equation p pressure Q diabatic heating rate per mass R = 287 J kg−1 gas constant of dry air D − @lnθ S T @p stability parameter in pressure coordinates t time T temperature V horizontal wind vector V g geostrophic wind vector α relaxation coefficient D − RT @θ σ pθ @p hydrostatic stability σ0 isobaric mean of hydrostatic stability ζ vertical component of relative vorticity ζg relative vorticity of geostrophic wind ζag relative vorticity of ageostrophic wind D dp ! dt isobaric vertical motion r horizontal nabla operator r2 horizontal laplacian operator @ @ζ R 2 @T FA D f C r : (11) from the geostrophic vorticity tendency. Neglecting the vari- @p @t p @t ation of the Coriolis parameter, Apart from the reorganisation of the terms in Eq. (6), this g ζ D r2Z (12) generalised omega equation is identical with the one used by g f Räisänen(1995). It follows directly from the isobaric prim- and hence itive equations, which assume hydrostatic balance but omit @Z f @ζg the other approximations in the QG theory.
Recommended publications
  • REVIEW the Quasigeostrophic Omega Equation: Reappraisal
    VOLUME 143 MONTHLY WEATHER REVIEW JANUARY 2015 REVIEW The Quasigeostrophic Omega Equation: Reappraisal, Refinements, and Relevance HUW C. DAVIES Institute for Atmospheric and Climate Science, ETH Zurich, Zurich, Switzerland (Manuscript received 29 March 2014, in final form 1 September 2014) ABSTRACT A two-component study is undertaken of the classical quasigeostrophic (QG) omega equation. First, a reappraisal is undertaken of extant formulations of the equation’s so-called forcing function. It pinpoints shortcomings of various formulations and prompts consideration of alternative forms. Particular consider- ation is given to the contribution of flow deformation to the forcing function, and to the role of the advection of the geostrophic flow by the thermal wind (the R vector). The latter is closely related to the Q vector, the horizontal component of the ageostrophic vorticity, and the forcing function itself. The reexamination pro- motes further examination of the physical interpretation and diagnostic use of the omega equation particu- larly for assessing richly structured subsynoptic flow features. Second, consideration is given to the dynamics associated with the equation and its more general utility. It is shown that the R vector is intrinsic to a quasigeostrophic cascade to finer-scaled flow, and that a fundamental feature of the QG omega equation—the in-phase relationship between cloud-diabatic heating and the at- tendant vertical velocity—has important potential ramifications for the assimilation of data in numerical weather prediction (NWP) models. Finally, it is shown that, in the context of considering NWP model output, mild generalizations of the quasigeostrophic R vector retain interpretative value for flow settings beyond geostrophy and warrant consideration when addressing some contemporary NWP challenges.
    [Show full text]
  • Derivation and Solution of the Omega Equation Associated with A
    PUBLICATIONS Journal of Advances in Modeling Earth Systems RESEARCH ARTICLE Derivation and Solution of the Omega Equation Associated 10.1002/2017MS000992 With a Balance Theory on the Sphere Key Points: John F. Dostalek1 , Wayne H. Schubert2 , and Mark DeMaria3 A global balance theory, similar to quasi-geostrophic theory, can be 1Cooperative Institute for Research in the Atmosphere, Colorado State University, Fort Collins, CO, USA, 2Department of developed by decomposing the wind Atmospheric Science, Colorado State University, Fort Collins, CO, USA, 3National Hurricane Center, Miami, FL, USA into nondivergent and irrotational components A global omega equation, similar to the quasi-geostrophic omega Abstract The quasi-geostrophic omega equation has been used extensively to examine the large-scale equation, can be derived from the vertical velocity patterns of atmospheric systems. It is derived from the quasi-geostrophic equations, a global balance system The global omega equation may be balanced set of equations based on the partitioning of the horizontal wind into a geostrophic and an solved over the sphere using vertical ageostrophic component. Its use is limited to higher latitudes, however, as the geostrophic balance is normal modes and spherical undefined at the equator. In order to derive an omega equation which can be used at all latitudes, a new harmonics balanced set of equations is developed. Three key steps are used in the formulation. First, the horizontal wind is decomposed into a nondivergent and an irrotational component. Second, the Coriolis parameter is Correspondence to: J. F. Dostalek, assumed to be slowly varying, such that it may be moved in and out of horizontal derivative operators as [email protected] necessary to simplify the derivation.
    [Show full text]
  • A New Method of Calculating Vertical Motion in Isentropic Space
    A NEW METHOD OF CALCULATING VERTICAL MOTION IN ISENTROPIC SPACE _______________________________________ A Thesis presented to the Faculty of the Graduate School at the University of Missouri-Columbia _______________________________________________________ In Partial Fulfillment of the Requirements for the Degree Soil, Environmental, and Atmospheric Sciences _____________________________________________________ by MICHEAL JOSEPH SIMPSON Dr. Patrick Market, Thesis Supervisor MAY 2014 The undersigned, appointed by the dean of the Graduate School, have examined the thesis entitled A New Method of Calculating Vertical Motion in Isentropic Space presented by Micheal Simpson, a candidate for the degree of master of atmospheric science, and hereby certify that, in their opinion, it is worthy of acceptance. Professor Patrick Market Professor Anthony Lupo Professor Jason Hubbart Professor Scott Rochette This thesis is dedicated to my parents, who have never stopped believing in me. Without their love and support, the opportunity to write this paper would not have been possible. ACKNOWLEDGEMENTS The Weather Forecast Office in Springfield, Missouri provided the impetus for the writing of this paper; without the individuals at this office, this paper would not have been possible. I am eternally grateful to everyone at the Springfield WFO. Dr. Scott Rochette and Dr. Patrick Market must be thanked for serving as my advisors and helping me get to where I am today. My committee members: Dr. Patrick Market, Dr. Scott Rochette, Dr. Anthony Lupo, and Dr. Jason Hubbart for taking the time and effort to be on my committee; their guidance as professionals is inspiring and appreciated beyond words. Dr. Neil Fox and Andrew Jensen have helped with the writing and mathematics of this paper, and their help is greatly appreciated.
    [Show full text]
  • Applying Trenberth's Approximation of the Omega Equation
    Western Region Technical Attachment No. 91-48 December 10, 1991 APPLYING TRENBERTH'S APPROXIMATION OF THE OMEGA EQUATION The quasi-geostrophic omega equation is a valuable tool for diagnosing areas of synoptic scale vertical motion from common weather maps. The most traditional way of estimating vertical motion is to identify areas of positive (or cyclonic) vorticity advection (PVA) at 500mb. However, the quasi-geostrophic vertical motion in the omega equation is actually related to the PVA increasing with height and the Laplacian of temperature advection. Many studies have shown that the temperature advection term can be as important to the overall vertical motion field as the PVA term and, in many cases, may cancel the contribution from the PVA term. Furthermore, the 500mb PVA pattern doesn't always look the same as the pattern of where PVA is increasing with height (for example, Durran and Snellman, 1987). Thus, several ways of interpreting the complete quasi-geostrophic omega equation that do not rely on a two-term approach have been developed. Although the Q-vector method is accurate (see WRTA 90-07), it is difficult to apply with current forecast maps available on AFOS. Perhaps the easiest to use with forecast maps that are currently available on AFOS is the Trenberth approximation (Trenberth, 1978). Trenberth shows that the advection of positive (cyclonic) vorticity by the thermal wind is a relatively good approximation to the sum of the PVA and temperature advection terms of the full omega equation. This approximation is most valid in the middle atmosphere between about 600mb and 400mb.
    [Show full text]
  • The Q-Vector Form of the Quasi-Geostrophic Omega Equation
    Synoptic Meteorology II: The Q-Vector Form of the Omega Equation Readings: Section 2.3 of Midlatitude Synoptic Meteorology. Motivation: Why Do We Need Another Omega Equation? The quasi-geostrophic omega equation is an excellent diagnostic tool and has served as the one of the foundations for introductory synoptic-scale weather analysis for over 50 years. However, it is not infallible. Rather, it has two primary shortcomings: • The two primary forcing terms in the omega equation, the differential geostrophic-vorticity advection and Laplacian of the potential-temperature advection terms, often have different signs from one another. Without computing each term’s actual magnitude, it is impossible to assess which one is larger (and is the primary control on synoptic-scale vertical motions). • The terms in the omega equation are sensitive to the reference frame – stationary or moving with the flow – in which they are computed. In other words, different results are obtained if the reference frame is changed, even if the meteorological features are the same! This is not a primary concern for us, however, as we are primarily considering stationary reference frames (e.g., as depicted on standard weather charts). These shortcomings motivate a desire to obtain a new equation for synoptic-scale vertical motions that does not have these problems. This equation, known as the Q-vector form of the quasi- geostrophic omega equation, is derived below. Obtaining The Q-Vector Form of the Quasi-Geostrophic Omega Equation To obtain the Q-vector form of the quasi-geostrophic omega equation, we will start with the quasi- geostrophic horizontal momentum and thermodynamic equations.
    [Show full text]
  • Synoptic Meteorology II: Quasi-Geostrophic Omega
    Synoptic Meteorology II: Quasi-Geostrophic Omega Equation Examples The below images, obtained from the University of Arizona Real-Time QG Diagnostics webpage (https://tgalarneau.faculty.arizona.edu/content/real-time-qg-diagnostics), provide data that we can use to evaluate the forcing terms to the quasi-geostrophic omega equation. Two notes: • When evaluating the differential vorticity advection term, vorticity advection on the lower of the two isobaric surfaces (e.g., at 850 hPa if evaluating the omega equation at 700 hPa) is often assumed to be negligible compared to that on the upper of the two isobaric surfaces (e.g., at 500 hPa). Thus, the analyses below – which emphasize forcing for vertical motion at 700 hPa – use 500 hPa geostrophic advection of geostrophic relative vorticity as a proxy for the differential advection term. On the website above, the relevant plots are listed under the quasi-geostrophic height tendency equation section. • As before, we neglect the friction and diabatic heating terms to the full omega equation. We first focus on a shortwave trough in the central United States at 1200 UTC 20 February 2019. We see that 500 hPa geostrophic absolute vorticity is maximized in the base of this trough along 100°W (Fig. 1). Given the first note above, this implies that there is differential cyclonic vorticity advection to the east and northeast and differential anticyclonic vorticity advection to the west and southwest, which is confirmed by Fig. 2. With respect to the Laplacian of the temperature advection term, we see that 700 hPa temperature advection is a local maximum east of the trough and local minimum southwest of the trough (Fig.
    [Show full text]
  • An Introduction to Dynamic Meteorology.Pdf
    January 24, 2004 12:0 Elsevier/AID aid An Introduction to Dynamic Meteorology FOURTH EDITION January 24, 2004 12:0 Elsevier/AID aid This is Volume 88 in the INTERNATIONAL GEOPHYSICS SERIES A series of monographs and textbooks Edited by RENATADMOWSKA, JAMES R. HOLTON and H.THOMAS ROSSBY A complete list of books in this series appears at the end of this volume. January 24, 2004 12:0 Elsevier/AID aid AN INTRODUCTION TO DYNAMIC METEOROLOGY Fourth Edition JAMES R. HOLTON Department of Atmospheric Sciences University of Washington Seattle,Washington Amsterdam Boston Heidelberg London New York Oxford Paris San Diego San Francisco Singapore Sydney Tokyo January 24, 2004 12:0 Elsevier/AID aid Senior Editor, Earth Sciences Frank Cynar Editorial Coordinator Jennifer Helé Senior Marketing Manager Linda Beattie Cover Design Eric DeCicco Composition Integra Software Services Printer/Binder The Maple-Vail Manufacturing Group Cover printer Phoenix Color Corp Elsevier Academic Press 200 Wheeler Road, Burlington, MA 01803, USA 525 B Street, Suite 1900, San Diego, California 92101-4495, USA 84 Theobald’s Road, London WC1X 8RR, UK This book is printed on acid-free paper. Copyright c 2004, Elsevier Inc. All rights reserved. No part of this publication may be reproduced or transmitted in any form or by any means, electronic or mechanical, including photocopy, recording, or any information storage and retrieval system, without permission in writing from the publisher. Permissions may be sought directly from Elsevier’s Science & Technology Rights Department in
    [Show full text]
  • 1 QG Omega Equation 1. Traditional Form a Diagnostic Expression For
    QG Omega Equation 1. Traditional form A diagnostic expression for vertical motion is derived by manipulating the QG vorticity and thermodynamic equations. The resulting equation for frictionless, adiabatic flow is known as the QG omega equation [eq. 5.6.11 in Bluestein (1992), p. 329] $ �' � �' � �; !∇$ + + � = − /−� ∙ � 5� + �89 − ∇$ 5−� ∙ � �8 # � ��$ � �� � � 7 �� # � � $ where � is the vertical velocity, �� is the gradient operator on a pressure surface, ∇# is the Laplacian operator on a pressure surface, �� is the geostrophic wind, �7 is the geostrophic relative vorticity, �; is the gas constant for dry air, � is the pressure, � is the Coriolis parameter, -4 -1 �' is the Coriolis constant (10 s ), � is the static stability parameter (assumed constant at 2.0 × 10BC �$��B$�B$), and � is the temperature. The right-hand-side (RHS) is calculated to determine the forcing for vertical motion at the 700 hPa level. Positive values indicate forcing for ascent. Negative values indicate forcing for descent. Note that this diagnostic calculation results in the forcing for vertical motion not the actual vertical motion because we do not solve the Laplacian on the left-hand-side. The first term (A) on the RHS is differential advection of geostrophic absolute vorticity by the geostrophic wind. This term is evaluated using vorticity advection at 900 and 500 hPa. Finite differencing is used to evaluate the vertical derivative of vorticity advection. Cyclonic vorticity advection increasing with height and maximized at 500 hPa (term A > 0) is associated with forcing for ascent at 700 hPa. Anticyclonic vorticity advection increasing with height and maximized at 500 hPa (term A < 0) is associated with forcing for descent at 700 hPa.
    [Show full text]
  • Lectures on Dynamical Meteorology
    LECTURES ON DYNAMICAL METEOROLOGY Roger K. Smith Version: June 16, 2014 Contents 1 INTRODUCTION 5 1.1 Scales ................................... 6 2 EQUILIBRIUM AND STABILITY 9 3 THE EQUATIONS OF MOTION 16 3.1 Effectivegravity.............................. 16 3.2 TheCoriolisforce............................. 16 3.3 Euler’s equation in a rotating coordinate system . 18 3.4 Centripetalacceleration .. .. .. 19 3.5 Themomentumequation.. .. .. 20 3.6 TheCoriolisforce............................. 20 3.7 Perturbationpressure. .. .. .. 21 3.8 Scale analysis of the equation of motion . 22 3.9 Coordinate systems and the earth’s sphericity . 23 3.10 Scale analysis of the equations for middle latitude synoptic systems . 25 4 GEOSTROPHIC FLOWS 28 4.1 TheTaylor-ProudmanTheorem . 30 4.2 Blocking.................................. 34 4.3 Analogy between blocking and axial Taylor columns . 35 4.4 Stability of a rotating fluid . 38 4.5 Vortexflows: thegradientwindequation . 38 4.6 Theeffectsofstratification. 41 4.7 Thermaladvection ............................ 45 4.8 Thethermodynamicequation . 46 4.9 Pressurecoordinates ........................... 47 4.10 Thicknessadvection. .. .. .. 48 4.11 Generalizedthermalwindequation . 49 5 FRONTS, EKMAN BOUNDARY LAYERS AND VORTEX FLOWS 54 5.1 Fronts ................................... 54 5.2 Margules’model.............................. 54 5.3 Viscous boundary layers: Ekman’s solution . 59 2 CONTENTS 3 5.4 Vortexboundarylayers. .. .. .. 63 6 THE VORTICITY EQUATION FOR A HOMOGENEOUS FLUID 67 6.1 Planetary,orRossbyWaves . 68 6.2 Largescaleflowoveramountainbarrier . 74 6.3 Winddrivenoceancurrents . 75 6.4 Topographicwaves ............................ 79 6.5 Continentalshelfwaves. .. .. .. 81 7 THE VORTICITY EQUATION IN A ROTATING STRATIFIED FLUID 83 7.1 The vorticity equation for synoptic-scale atmospheric motions . ... 85 8 QUASI-GEOSTROPHIC MOTION 89 8.1 More on the approximated thermodynamic equation . 92 8.2 The quasi-geostrophic equation for a compressible atmosphere ...
    [Show full text]
  • Omega Equation Atmos 5110 Synoptic–Dynamic Meteorology I Instructor: Jim Steenburgh [email protected] 801-581-8727 Suite 480/Office 488 INSCC
    QG Theory and Applications: Omega Equation Atmos 5110 Synoptic–Dynamic Meteorology I Instructor: Jim Steenburgh [email protected] 801-581-8727 Suite 480/Office 488 INSCC Suggested reading: Lackmann (2011), Section 2.3 Usage The omega equation is used to diagnose large-scale vertical motion. Derivation Derived from the QG thermodynamic and vorticity equations. See Lakcmann (2011) for details. The Omega Equation �# �# � � 1 1 �Φ # & & 0⃗ # # 0⃗ !∇ + #+ � = .�2 ∙ ∇ 4 ∇ Φ + f89 + ∇ [�2 ∙ ∇ 4− 8] � �� � �� �& � �� Yeah, it’s daunting, but easier if you think of it as A=B+C where: �# �# � = !∇# + & + � � ��# ?@ B F # � = D�0⃗2 ∙ ∇ E ∇ Φ + fGH = Differential vorticity advection term A BC ?@ F BK � = �#[�0⃗ ∙ � E− G] = Temperature advection term A 2 BC Term A Essentially the 3-D Laplacian acting on ω. For sinusoidal (wave-like) patterns, the Laplacian can be approximated by a minus sign. 1 �# �# � = !∇# + & + � ~ − � ∝ � � ��# Therefore, �& � 1 # 1 # �Φ � ∝ −� ∝ .�0⃗2 ∙ ∇ 4 ∇ Φ + f89 + ∇ [�0⃗2 ∙ ∇ 4− 8] � �� �& � �� Term B Differential vorticity advection term � 0 � 00⃗ 1 2 � 00⃗ � ∝ !�� ∙ ∇ Q ∇ Φ + fS+ ∝ D−�� ∙ ∇U�2 + fWH � �� �0 �� Therefore • Cyclonic vorticity advection (positive in NH) increasing with height indicates ascent • Anticyclonic vorticity advection (negative in NH) increasing with height indicates descent Typically we assume that the vorticity advection at low levels is weak so that • Cyclonic vorticity advection (CVA) @ 500 mb indicates ascent • Anticyclonic vorticity advection (AVA) @ 500 mb indicates descent 2 0⃗ 0⃗ In
    [Show full text]
  • Combining the Two Forcing Terms in the QG Omega Equation
    Combining the Two Forcing Terms in the QG­ Omega Equation In this example, it appears that the vorticity advection and thermal advection terms are opposing one another. How can one decide which is more significant? In this case, it "appears" that the warm advection term is liable to be stronger than the differential negative vorticity advection. This is very difficult to estimate, however, from a visual inspection of patterns. Click on chart to see full size version. Click on chart to see full size version. The two terms can be combined mathematically. One method yields an expression that states that quasigeostrophic omega is proportional to the ADVECTION OF THE GEOSTROPHIC VORTICITY BY THE THERMAL WIND. There are other terms in the equation that drop out on an order of magnitude basis. This is called the TRENBERTH APPROXIMATION. To visualize this you need only two charts. One showing the absolute (geostrophic) vorticity at a given level (we are cheating and using the real absolute vorticity) (above right) and the other showing the thickness field where the level at which you would like to estimate omega is somewhere in the middle of the layer. In this case, we are using 1000­500 mb thickness (above left). The "thermal wind" blows parallel to the thickness contours, clockwise around highs and counterclockwise around lows. (It is a kind of geostrophic wind). Where the thermal wind blows from high to low values of vorticity (positive isothermal vorticity advection­­piva or positive or cyclonic voriticty advection by the thermal wind), there is quasigeostrophic forcing for upwards omega and vice versa.
    [Show full text]
  • Subtropical Floods and Tropical Cyclones Daniel A. Shaevitz
    Extreme weather: subtropical floods and tropical cyclones Daniel A. Shaevitz Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in the Graduate School of Arts and Sciences COLUMBIA UNIVERSITY 2016 c 2016 Daniel A. Shaevitz All Rights Reserved ABSTRACT Extreme weather: subtropical floods and tropical cyclones Daniel A. Shaevitz Extreme weather events have a large effect on society. As such, it is important to understand these events and to project how they may change in a future, warmer climate. The aim of this thesis is to develop a deeper understanding of two types of extreme weather events: subtropical floods and tropical cyclones (TCs). In the subtropics, the latitude is high enough that quasi-geostrophic dynamics are at least qualitatively relevant, while low enough that moisture may be abun- dant and convection strong. Extratropical extreme precipitation events are usu- ally associated with large-scale flow disturbances, strong ascent, and large latent heat release. In the first part of this thesis, I examine the possible triggering of convection by the large-scale dynamics and investigate the coupling between the two. Specifically two examples of extreme precipitation events in the subtropics are analyzed, the 2010 and 2014 floods of India and Pakistan and the 2015 flood of Texas and Oklahoma. I invert the quasi-geostrophic omega equation to decom- pose the large-scale vertical motion profile to components due to synoptic forcing and diabatic heating. Additionally, I present model results from within the Column Quasi-Geostrophic framework. A single column model and cloud-revolving model are forced with the large-scale forcings (other than large-scale vertical motion) com- puted from the quasi-geostrophic omega equation with input data from a reanalysis data set, and the large-scale vertical motion is diagnosed interactively with the simu- lated convection.
    [Show full text]