Quasi-3D Multiscale Modeling Framework as a Physics Option in CAM-SE
|
|
- Augustus Tyrone Hudson
- 5 years ago
- Views:
Transcription
1 Quasi-3D Multiscale Modeling Framework as a Physics Option in CAM-SE Joon-Hee Jung, Celal S. Konor, Don Dazlich, David A. Randall, Department of Atmospheric Sciences, Colorado State University, USA Peter Lauritzen, and Steve Goldhaber National Center for Atmospheric Research, Boulder, Colorado, USA This research has been supported by NSF AGS (and partly by DOE ACME DE-SC and DOE CMDV DE-SC ). CESM: Atmosphere Model Working Group Meeting NCAR, Mesa Lab, Boulder, CO, February 12-14, 2018
2 Multiscale Modeling Framework (aka Superparameterization) GCM Cloud Resolving Model (CRM) Traditional subgrid-scale parameterizations in the GCMs are replaced with explicit simulations by a 2D CRM embedded in each GCM grid column. MMF community: CSU, State Univ. of New York at Stony Brook, Univ. of Washington, Univ. of Chicago, UC at Irvine, Univ. of Oxford, Harvard Univ., George Mason Univ., NCAR, NASA GSFC & LRC, ECMWF, LBNL, SIO, ESRL, PNNL, MIT, IITM, DOE (E3SM)
3 Q3D MMF Jung and Arakawa (2010, 2014), Jung (2016) CRMs in GCM grid columns are seamlessly connected; Two perpendicular sets of CRMs are used; CRMs are three-dimensional, although covering only channel-like domains. Q3D MMF MMF y y z z x x GCM grid cell GCM grid CRM grid CRM ghost grid
4 Q3D MMF Jung and Arakawa (2010, 2014), Jung (2016) CRMs in GCM grid columns are seamlessly connected; Two perpendicular sets of CRMs are used; CRMs are three-dimensional, although covering only channel-like domains. Convections are allowed to propagate across the GCM cell boundaries; The surface orography can be resolved by the CRMs; Vertical momentum transport by convection and waves can be simulated. Q3D MMF MMF y y z z x x GCM grid cell GCM grid CRM grid CRM ghost grid
5 Q3D MMF inherits the structure of the conventional GCMs in Coupling the Dynamical Core and Physics FORCING GCM Effects on CRM GCM Dynamical Core FEEDBACK CRM Effects on GCM CRM Parameterized Physics This is like the first-generation MMF.
6 Why does the Q3D MMF still need a separate GCM dynamical core?
7 Why does the Q3D MMF still need a separate GCM dynamical core? The network of CRM channels cannot replace a GCM. The CRM channels do not and cannot directly interact with each other; Averaging the solutions of the CRM channels does not represent fully interacting 3D large-scale dynamical processes; To simulate well-behaving large-scale features, a reasonable 3D GCM dynamics is necessary; The computational cost of a GCM dynamics is not a burden comparing to that of CRM channels; It is good to maintain compatibility with conventional GCMs.
8 There is no double counting or spurious competition between the GCM and CRM solutions.
9 There is no double counting or spurious competition between the GCM and CRM solutions. GCM solution at t=t0 + dtgcm CRM-mean within a GCM cell Eddy transport and diabatic effects added to GCM s dynamical effects t = t0 + dtgcm t ep s GCM solution at t=t0 t0 x GCM grid GCM solution interpolated to CRM grids R M st dtgcm C G C M st ep t t = t0 Decomposition of a CRM variable: background-field CRM grid x (GCM grid-scale) + eddy GCM gathers only eddy transport and diabatic effects (i.e., subgrid-scale effects) from the CRMs. CRM maintains its GCM grid-scale solution close to GCM s solution through a relaxation.
10 Scale-awareness of the Q3D MMF All subgrid-scale processes are represented by the CRMs with the same grid size, regardless of the resolution of the GCM; CRMs can figure out the GCM s resolution; For example, with a fine-grid, the vertical velocities of the GCM can be much larger (e.g., in a mesoscale system) and also the horizontal gradients can be much larger (e.g., near a front). The feedbacks of CRMs are sensitive to the GCM s resolution, because the collective effects of subgrid-scale processes will correspond to the smallest scale (i.e., grid size) of the GCM.
11 Development of a Global Q3D MMF A North-South CRM channel CAM-SE dynamical core + Vector Vorticity Model (CRM developed at CSU/UCLA) An East-West CRM channel
12 Development of a Global Q3D MMF A North-South CRM channel CAM-SE dynamical core + Vector Vorticity Model (CRM developed at CSU/UCLA) An East-West CRM channel The VVM has been developed and applied as a limited-area model. To be applied globally (i.e., wrapped around the planet), it has been extended to the curvilinear coordinates. To couple the VVM with the CAM-SE dynamical core, a computational infrastructure has been designed and is now under construction.
13 The VVM predicts the Vorticity Vector. Large-scale motions are controlled by the vertical component of vorticity Small-scale motions are controlled by the horizontal component of vorticity
14 The VVM predicts the Vorticity Vector. Large-scale motions are controlled by the vertical component of vorticity Small-scale motions are controlled by the horizontal component of vorticity 3D nonhydrostatic anelastic model; Pressure gradient force is eliminated; Vertical velocity is a solution of a 3D elliptic equation; (Lower boundary condition over steep topography is easily formulated.) CRM-type physics parameterizations are included. Suitable for the Q3D-type application
15 Governing Equations of the VVM in the Curvilinear Coordinates Vorticity Equation ξ 1 = t G x η 1 = t G x ( G uξ + y ) ( u u 1 θ u G vξ ( wξ ) + ξ + η + (ζ + f ) + g y x z G y θ 0 z ) ( G uη + y ) ( v v 1 θ v G vη ( wη ) + ξ + η + (ζ + f ) g y x z G x θ 0 z ) w w ζ 1 w = Gu ζ + f + Gv ζ + f w ζ + f + ξ + η + ζ + f ( ) ( ) ( ) ( ) z t x y z x y G where x = q1 2 y = q u = u1 2 v = u u! = u1 v! = u2 G : Jacobian of transformation 1 ξ = (ω H ) = 2 η = (ω H ) = ζ = ω z = 1 w v! y z G 1 u! w z x G 1 v! u! x y G
16 Governing Equations of the VVM (Cont.) Nondivergence of 3D vorticity ζ = z zt 1 G x ( ) Gξ + y ( ) Gη dz + (ζ )z=z T Continuity Equation 1 G x ( ) Gρ0u + y ( (ρ0 w ) Gρ0 v + =0 z ) w-equation 1 G 11 w 21 w 12 w 22 w +g + G g +g G g x y y x y x! ξ! 1 (ρ0 w ) 1 η + = z ρ0 z G x y
17 Latitude (degrees) EQ - VVM: Barotropic Instability Test SPHERICAL LOW RESOLUTION t = 96 hr dx=dy ~ 100 km EQ - (dx=dy ~ 100 km) CUBED-SPHERE dx max =dy max ~ 100 km dx min =dy min ~ 84 km t = 96 hr Latitude (degrees) EQ - t = 144 hr Longitude (degrees) EQ - t = 144 hr Longitude (degrees) (s -1 )
18 Latitude (degrees) EQ - VVM: Barotropic Instability Test SPHERICAL LOW RESOLUTION t = 96 hr dx=dy ~ 100 km EQ - (dx=dy ~ 100 km) CUBED-SPHERE dx max =dy max ~ 100 km dx min =dy min ~ 84 km t = 96 hr Latitude (degrees) Latitude (degrees) Latitude (degrees) EQ Longitude (degrees) EQ - EQ - SPHERICAL Longitude (degrees) t = 144 hr EQ - HIGH RESOLUTION t = 96 hr t = 144 hr dx=dy ~ 5 km t = 144 hr Longitude (degrees) EQ EQ - (dx=dy ~ 5 km) CUBED-SPHERE Longitude (degrees) dx max =dy max ~ 5 km dx min =dy min ~ 4 km t = 96 hr t = 144 hr (s -1 )
19 VVM: Baroclinic Instability Test LOW RESOLUTION (dx = dy ~ 100 km) Zonally Uniform Initial State ( ) Idealized setting: f = 2Ωsin ϕ + π 4 Height (km) Latitude (deg) Initial Random Perturbation Longitude (deg) Latitude (deg) (K) Latitude (deg) (m/s) Potential Temp. at the lowest model layer (Day 12) SPHERICAL CUBED-SPHERE Longitude (deg) Latitude (deg) Longitude (deg)
20 Computational Infrastructure for coupling the VVM with CAM-SE Dycore Collaboration with Peter Lauritzen and Steve NCAR
21 Computational Infrastructure for coupling the VVM with CAM-SE Dycore Collaboration with Peter Lauritzen and Steve NCAR Using the finite-volume physics grid for CAM-SE, including the coupling to the spectral element dynamics grid (CAM-SE-physgrid). This physics grid forms the virtual GCM (vgcm) grid for the Q3D system. dynamics grid y physics grid y y x vgcm grid x x CRM grid CRM ghost grid
22 Computational Infrastructure (Cont.) Computing the (cube) face-edge ghost cells in the dycore, which have never existed in CAM before (these ghost cells are not part of the original physics grid). Background Field (GCM grid scale) of the VVM is interpolated from the vgcm fields. y x Moving the vgcm data from the dynamics decomposition to the physics decomposition.
23 Computational Infrastructure (Cont.) (In the physics decomposition) F5 Organizing all vgcm data (including the ghost cells) in a way that they can be properly assigned to the CRM channels. F4 F4 F1 F2 F6 Group 1 Group 2 Group 3 F5 F5 F5 F1 F6 F2 F3 F4 F1 F6 F2 F3 F3 F4 F1 F6 F2 F3
24 Computational Infrastructure (Cont.) Moving the vgcm data from the physics decomposition to the channel decomposition: Each vgcm s data is placed to the correct location of CRM channel. (Also, moving the mean CRM data from the channel decomposition to the physics decomposition.) channel decomposition is where each task has all the data for one or more CRM channels. (No MPI communications are necessary for the CRM calculations.)
25 POST-CRM Procedures DYCORE COMPUTE Distribution of CRM feedbacks to Dycore points vgcm communications (moving data to physics decomposition) CRM PREDICTION Calculation of mean CRM feedbacks (physics tendencies) CRM COMPUTE PRE-CRM Procedures CRM COMPUTE CRM COMPUTE CRM COMPUTE Interpolation of vgcm data to CRM points vgcm communications (moving data to channel decomposition) Generation of data on physics grids Q3D MMF Computation Algorithm Similar to the coupling approach of CAM-SE Q3D CAM-SE PREDICTION (Dycore + CRM effects) physics t0 t
26 Summary Our goal is to make global versions of the Q3D MMF: Q3D MMF as a Physics Option in CAM-SE (also in E3SM) Finished the preparation of the VVM with the standard cubedsphere grid, which will be coupled with the CAM-SE Dycore to form a global Q3D MMF. Created the Q3D physics package in a CAM branch. Making progress in designing and constructing a computational infrastructure (started implementing the codes and the communication layout).
MODEL UNIFICATION my latest research excitement Akio Arakawa
MODEL UNIFICATION my latest research excitement Akio Arakawa Department of Atmospheric and Oceanic Sciences, UCLA CMMAP, January 7, 24 Wayne Schubert ` 7 Cumulus/ L-S interaction David Randall Wayne Schubert
More informationNCAR Global Atmospheric Core Workshop, Boulder, June 2008
NCAR Global Atmospheric Core Workshop, Boulder, June 2008 D. Majewski based on Christiane Jablonowski (email: cjablono@umich.edu) University of Michigan Goals of the Test Suite NASA/GFDL Test cases should
More informationCESM2 release of CAM-SE (& CAM-SE-CSLAM)
CESM2 release of CAM-SE (& CAM-SE-CSLAM) Peter Hjort Lauritzen National Center for Atmospheric Research Boulder, Colorado, USA Collaborators: S. Goldhaber (NCAR), J. Bacmeister (NCAR), R.D.Nair (NCAR),
More informationMODEL TYPE (Adapted from COMET online NWP modules) 1. Introduction
MODEL TYPE (Adapted from COMET online NWP modules) 1. Introduction Grid point and spectral models are based on the same set of primitive equations. However, each type formulates and solves the equations
More informationAfrican Easterly Waves and Convection: A Potential Vorticity Perspective. James Russell & Anantha Aiyyer North Carolina State University
African Easterly Waves and Convection: A Potential Vorticity Perspective James Russell & Anantha Aiyyer North Carolina State University Background 1. Introduction Early studies described AEWs as dry adiabatic
More informationWojciech Grabowski. National Center for Atmospheric Research, USA
Numerical modeling of multiscale atmospheric flows: From cloud microscale to climate Wojciech Grabowski National Center for Atmospheric Research, USA This presentation includes results from collaborative
More informationUsing Cloud-Resolving Models for Parameterization Development
Using Cloud-Resolving Models for Parameterization Development Steven K. Krueger University of Utah! 16th CMMAP Team Meeting January 7-9, 2014 What is are CRMs and why do we need them? Range of scales diagram
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 informationReduced Complexity Frameworks for Exploring Physics Dynamics Coupling Sensitivities
Reduced Complexity Frameworks for Exploring Physics Dynamics Coupling Sensitivities Kevin A. Reed & Adam R. Herrington School of Marine and Atmospheric Sciences Stony Brook University, Stony Brook, New
More informationIdealized Nonhydrostatic Supercell Simulations in the Global MPAS
Idealized Nonhydrostatic Supercell Simulations in the Global Joe Klemp, Bill Skamarock, and Sang-Hun Park National Center for Atmospheric Research Boulder, Colorado Typical characteristics: Supercell Thunderstorms
More informationA Truncated Model for Finite Amplitude Baroclinic Waves in a Channel
A Truncated Model for Finite Amplitude Baroclinic Waves in a Channel Zhiming Kuang 1 Introduction To date, studies of finite amplitude baroclinic waves have been mostly numerical. The numerical models,
More informationA moist static energy budget analysis of the MJO in the Superparameterized Community Atmosphere Model
A moist static energy budget analysis of the MJO in the Superparameterized Community Atmosphere Model Mike Pritchard University of Washington in collaboration w. Chris Bretherton & The Center for Multiscale
More informationParameterizing large-scale dynamics using the weak temperature gradient approximation
Parameterizing large-scale dynamics using the weak temperature gradient approximation Adam Sobel Columbia University NCAR IMAGe Workshop, Nov. 3 2005 In the tropics, our picture of the dynamics should
More informationModeling multiscale interactions in the climate system
Modeling multiscale interactions in the climate system Christopher S. Bretherton Atmospheric Sciences and Applied Mathematics University of Washington 08.09.2017 Aqua Worldview Motivation Weather and climate
More informationA Proposed Test Suite for Atmospheric Model Dynamical Cores
A Proposed Test Suite for Atmospheric Model Dynamical Cores Christiane Jablonowski (cjablono@umich.edu) University of Michigan, Ann Arbor PDEs on the Sphere Workshop Monterey, CA, June/26-29/2006 Motivation
More informationATS 421/521. Climate Modeling. Spring 2015
ATS 421/521 Climate Modeling Spring 2015 Lecture 9 Hadley Circulation (Held and Hou, 1980) General Circulation Models (tetbook chapter 3.2.3; course notes chapter 5.3) The Primitive Equations (tetbook
More informationPUBLICATIONS. Journal of Advances in Modeling Earth Systems. Simulation of subgrid orographic precipitation with an embedded 2-D cloud-resolving model
PUBLICATIONS Journal of Advances in Modeling Earth Systems RESEARCH ARTICLE 10.1002/2015MS000539 Key Points: Subgrid orographic precipitation is simulated with an embedded 2-D CRM The simulated mean precipitation
More informationMathematical Concepts & Notation
Mathematical Concepts & Notation Appendix A: Notation x, δx: a small change in x t : the partial derivative with respect to t holding the other variables fixed d : the time derivative of a quantity that
More informationModeling the atmosphere of Jupiter
Modeling the atmosphere of Jupiter Bruce Turkington UMass Amherst Collaborators: Richard S. Ellis (UMass Professor) Andrew Majda (NYU Professor) Mark DiBattista (NYU Postdoc) Kyle Haven (UMass PhD Student)
More informationTHE WEATHER RESEARCH AND FORECAST MODEL VERSION 2.0
THE WEATHER RESEARCH AND FORECAST MODEL VERSION 2.0 J. MICHALAKES, J. DUDHIA, D. GILL J. KLEMP, W. SKAMAROCK, W. WANG Mesoscale and Microscale Meteorology National Center for Atmospheric Research Boulder,
More informationUsing hybrid isentropic-sigma vertical coordinates in nonhydrostatic atmospheric modeling
Using hybrid isentropic-sigma vertical coordinates in nonhydrostatic atmospheric modeling Michael D. Toy NCAR Earth System Laboratory/ Advanced Study Program Boulder, Colorado, USA 9th International SRNWP
More informationABSTRACT 2 DATA 1 INTRODUCTION
16B.7 MODEL STUDY OF INTERMEDIATE-SCALE TROPICAL INERTIA GRAVITY WAVES AND COMPARISON TO TWP-ICE CAM- PAIGN OBSERVATIONS. S. Evan 1, M. J. Alexander 2 and J. Dudhia 3. 1 University of Colorado, Boulder,
More informationDynamical System Approach to Organized Convection Parameterization for GCMs. Mitchell W. Moncrieff
Dynamical System Approach to Organized Convection Parameterization for GCMs Mitchell W. Moncrieff Atmospheric Modeling & Predictability Section Climate & Global Dynamics Laboratory NCAR Year of Tropical
More informationSynoptic-Dynamic Meteorology in Midlatitudes
Synoptic-Dynamic Meteorology in Midlatitudes VOLUME II Observations and Theory of Weather Systems HOWARD B. BLUESTEIN New York Oxford OXFORD UNIVERSITY PRESS 1993 Contents 1. THE BEHAVIOR OF SYNOPTIC-SCALE,
More informationModel equations for planetary and synoptic scale atmospheric motions associated with different background stratification
Model equations for planetary and synoptic scale atmospheric motions associated with different background stratification Stamen Dolaptchiev & Rupert Klein Potsdam Institute for Climate Impact Research
More informationPV Generation in the Boundary Layer
1 PV Generation in the Boundary Layer Robert Plant 18th February 2003 (With thanks to S. Belcher) 2 Introduction How does the boundary layer modify the behaviour of weather systems? Often regarded as a
More informationM.Sc. in Meteorology. Physical Meteorology Prof Peter Lynch. Mathematical Computation Laboratory Dept. of Maths. Physics, UCD, Belfield.
M.Sc. in Meteorology Physical Meteorology Prof Peter Lynch Mathematical Computation Laboratory Dept. of Maths. Physics, UCD, Belfield. Climate Change???????????????? Tourists run through a swarm of pink
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 informationMacroturbulent cascades of energy and enstrophy in models and observations of planetary atmospheres
Macroturbulent cascades of energy and enstrophy in models and observations of planetary atmospheres Peter Read + Roland Young + Fachreddin Tabataba-Vakili + Yixiong Wang [Dept. of Physics, University of
More informationLatest thoughts on stochastic kinetic energy backscatter - good and bad
Latest thoughts on stochastic kinetic energy backscatter - good and bad by Glenn Shutts DARC Reading University May 15 2013 Acknowledgments ECMWF for supporting this work Martin Leutbecher Martin Steinheimer
More informationA Nodal High-Order Discontinuous Galerkin Dynamical Core for Climate Simulations
A Nodal High-Order Discontinuous Galerkin Dynamical Core for Climate Simulations Institute for Mathematics Applied to the Geosciences (IMAGe) National Center for Atmospheric Research (NCAR) Boulder CO
More information8/21/08. Modeling the General Circulation of the Atmosphere. Topic 4: Equatorial Wave Dynamics. Moisture and Equatorial Waves
Modeling the General Circulation of the Atmosphere. Topic 4: Equatorial Wave Dynamics D A R G A N M. W. F R I E R S O N U N I V E R S I T Y O F W A S H I N G T O N, D E P A R T M E N T O F A T M O S P
More informationToward unification of the multiscale modeling of the atmosphere
Atmos. Chem. Phys., 11, 3731 3742, 2011 doi:10.5194/acp-11-3731-2011 Author(s) 2011. CC Attribution 3.0 License. Atmospheric Chemistry and Physics Toward unification of the multiscale modeling of the atmosphere
More informationWRF MODEL STUDY OF TROPICAL INERTIA GRAVITY WAVES WITH COMPARISONS TO OBSERVATIONS. Stephanie Evan, Joan Alexander and Jimy Dudhia.
WRF MODEL STUDY OF TROPICAL INERTIA GRAVITY WAVES WITH COMPARISONS TO OBSERVATIONS. Stephanie Evan, Joan Alexander and Jimy Dudhia. Background Small-scale Gravity wave Inertia Gravity wave Mixed RossbyGravity
More informationOCEAN MODELING II. Parameterizations
OCEAN MODELING II Parameterizations Gokhan Danabasoglu Oceanography Section Climate and Global Dynamics Division National Center for Atmospheric Research NCAR is sponsored by the National Science Foundation
More informationFigure 1 - Resources trade-off. Image of Jim Kinter (COLA)
CLIMATE CHANGE RESEARCH AT THE EXASCALE Giovanni Aloisio *,, Italo Epicoco *,, Silvia Mocavero and Mark Taylor^ (*) University of Salento, Lecce, Italy ( ) Euro-Mediterranean Centre for Climate Change
More informationVariable-Resoluiton Global Atmospheric Modeling Spanning Convective to Planetary Scales
Variable-Resoluiton Global Atmospheric Modeling Spanning Convective to Planetary Scales Bill Skamarock, NCAR/MMM MPAS consists of geophysical fluid-flow solvers based on unstructured centroidal Voronoi
More informationCAM-SE: Lecture I. Peter Hjort Lauritzen
CAM-SE: Lecture I Peter Hjort Lauritzen Atmospheric Modeling and Predictability Section Climate and Global Dynamics Laboratory National Center for Atmospheric Research 2nd WCRP Summer School on Climate
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 informationAppendix B. Coordinate Transformation
B Appendix B. Coordinate Transformation B.1. Coordinate Transformation Relations In this appendix we derive and discuss the transformation relations between curvilinear and orthogonal coordinate systems,
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 informationTesting of a New Dynamical Core at CSU
Testing of a New Dynamical Core at CSU David A. Randall Department of Atmospheric Science, Colorado State University Fort Collins, Colorado, 80526 USA randall@atmos.colostate.edu 1 Introduction Since 2001,
More informationDiabatic processes and the structure of extratropical cyclones
Geophysical and Nonlinear Fluid Dynamics Seminar AOPP, Oxford, 23 October 2012 Diabatic processes and the structure of extratropical cyclones Oscar Martínez-Alvarado R. Plant, J. Chagnon, S. Gray, J. Methven
More informationDynamics of the Zonal-Mean, Time-Mean Tropical Circulation
Dynamics of the Zonal-Mean, Time-Mean Tropical Circulation First consider a hypothetical planet like Earth, but with no continents and no seasons and for which the only friction acting on the atmosphere
More informationSusan Bates Ocean Model Working Group Science Liaison
Susan Bates Ocean Model Working Group Science Liaison Climate Simulation Laboratory (CSL) Accelerated Scientific Discovery (ASD) NCAR Strategic Capability (NSC) Climate Process Teams (CPTs) NSF Earth System
More informationTheories on the Optimal Conditions of Long-Lived Squall Lines
Theories on the Optimal Conditions of Long-Lived Squall Lines References: Thorpe, A. J., M. J. Miller, and M. W. Moncrieff, 1982: Two-dimensional convection in nonconstant shear: A model of midlatitude
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 informationTheories on the Optimal Conditions of Long-Lived Squall Lines
Theories on the Optimal Conditions of Long-Lived Squall Lines References: Thorpe, A. J., M. J. Miller, and M. W. Moncrieff, 1982: Two -dimensional convection in nonconstant shear: A model of midlatitude
More informationTurbulence and Energy Transfer in Strongly-Stratified Flows
Turbulence and Energy Transfer in Strongly-Stratified Flows James J. Riley University of Washington Collaborators: Steve debruynkops (UMass) Kraig Winters (Scripps IO) Erik Lindborg (KTH) First IMS Turbulence
More informationEnforcing conservation of atmospheric axial angular momentum in CAM FV: Method and results in CESM2 and NorESM2 simulations
Enforcing conservation of atmospheric axial angular momentum in CAM FV: Method and results in CESM2 and NorESM2 simulations Thomas Toniazzo AMWG meeting, Boulder 28 February 2017 PH EW! Enforcing conservation
More informationFIM: A vertically flow-following, finite-volume icosahedral model
FIM: A vertically flow-following, finite-volume icosahedral model J.L. Lee 1, R. Bleck 2, A.E. MacDonald 1, J.-W. Bao 1, S. Benjamin 1, J. Middlecoff 3, N. Wang 3, and J. Brown 1 1 NOAA/Earth System Research
More informationSIMULATION OF GEOPHYSICAL TURBULENCE: DYNAMIC GRID DEFORMATION
SIMULATION OF GEOPHYSICAL TURBULENCE: DYNAMIC GRID DEFORMATION Piotr K Smolarkiewicz, and Joseph M Prusa, National Center for Atmospheric Research, Boulder, Colorado, U.S.A. Iowa State University, Ames,
More informationNonhydrostatic Icosahedral Model (NIM) A 3-D finite-volume NIM Jin Lee (+ other contributors)
Nonhydrostatic Icosahedral Model (NIM) A 3-D finite-volume NIM Jin Lee (+ other contributors) Earth System Research Laboratory (ESRL) NOAA/OAR GFDL,NSSL,ARL,AOML,GLERL,PMEL Aeronomy Lab. Climate Diagnostic
More informationInternal Dynamic Control of Hurricane Intensity Change: The Dual Nature of Potential Vorticity Mixing
Internal Dynamic Control of Hurricane Intensity Change: The Dual Nature of Potential Vorticity Mixing James Kossin University of Wisconsin Madison Space Science and Engineering Center Cooperative Institute
More informationAn Introduction to Physical Parameterization Techniques Used in Atmospheric Models
An Introduction to Physical Parameterization Techniques Used in Atmospheric Models J. J. Hack National Center for Atmospheric Research Boulder, Colorado USA Outline Frame broader scientific problem Hierarchy
More informationDevelopment of Yin-Yang Grid Global Model Using a New Dynamical Core ASUCA.
Development of Yin-Yang Grid Global Model Using a New Dynamical Core ASUCA. M. Sakamoto, J. Ishida, K. Kawano, K. Matsubayashi, K. Aranami, T. Hara, H. Kusabiraki, C. Muroi, Y. Kitamura Japan Meteorological
More information( u,v). For simplicity, the density is considered to be a constant, denoted by ρ 0
! Revised Friday, April 19, 2013! 1 Inertial Stability and Instability David Randall Introduction Inertial stability and instability are relevant to the atmosphere and ocean, and also in other contexts
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 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 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 informationChapter 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 informationClimate of an Earth- like Aquaplanet: the high- obliquity case and the <dally- locked case
Climate of an Earth- like Aquaplanet: the high- obliquity case and the
More informationClimate Change and Variability in the Southern Hemisphere: An Atmospheric Dynamics Perspective
Climate Change and Variability in the Southern Hemisphere: An Atmospheric Dynamics Perspective Edwin P. Gerber Center for Atmosphere Ocean Science Courant Institute of Mathematical Sciences New York University
More informationDevelopment and Validation of WACCM-X Thermosphere and Ionosphere
Development and Validation of WACCM-X Thermosphere and Ionosphere Han-Li Liu and WACCM-X Team NCAR/HAO: Ben Foster, Jing Liu, Gang Lu, Astrid Maute, Joe McInerney, Nick Pedatella, Liying Qian, Art Richmond,
More informationNon-local Momentum Transport Parameterizations. Joe Tribbia NCAR.ESSL.CGD.AMP
Non-local Momentum Transport Parameterizations Joe Tribbia NCAR.ESSL.CGD.AMP Outline Historical view: gravity wave drag (GWD) and convective momentum transport (CMT) GWD development -semi-linear theory
More informationChallenges and Opportunities in Modeling of the Global Atmosphere
Challenges and Opportunities in Modeling of the Global Atmosphere Zavisa Janjic, Vladimir Djurdjevic, Ratko Vasic, and Tom Black NOAA/NCEP/EMC, College Park, MD 20740, U.S.A. (Email: zavisa.janjic@noaa.gov)
More informationA Global Atmospheric Model. Joe Tribbia NCAR Turbulence Summer School July 2008
A Global Atmospheric Model Joe Tribbia NCAR Turbulence Summer School July 2008 Outline Broad overview of what is in a global climate/weather model of the atmosphere Spectral dynamical core Some results-climate
More informationEarth System Modeling Domain decomposition
Earth System Modeling Domain decomposition Graziano Giuliani International Centre for Theorethical Physics Earth System Physics Section Advanced School on Regional Climate Modeling over South America February
More informationCirculation and Vorticity. The tangential linear velocity of a parcel on a rotating body is related to angular velocity of the body by the relation
Circulation and Vorticity 1. Conservation of Absolute Angular Momentum The tangential linear velocity of a parcel on a rotating body is related to angular velocity of the body by the relation V = ωr (1)
More informationExploring Gravity Wave Dynamics, Sources, and Predictability in DeepWave
Exploring Gravity Wave Dynamics, Sources, and Predictability in DeepWave James D. Doyle 1, Stephen D. Eckermann 2, Eric Hendricks 1, Qingfang Jiang 1, P. Alex Reinecke 1, Carolyn A. Reynolds 1, David C.
More informationUncertainty in Ocean General Circulation Model Mixing Tensors
Uncertainty in Ocean General Circulation Model Mixing Tensors Baylor Fox-Kemper, University of Colorado Boulder, Brown U. (>Jan. 2013) Cooperative Insititute for Research in Environmental Sciences and
More informationSuper-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models PI Isaac Ginis Graduate School of Oceanography
More informationKinetic energy backscatter for NWP models and its calibration
Kinetic energy backscatter for NWP models and its calibration Glenn Shutts Met Office, Fitzroy Road EX1 3PB, United Kingdom glenn.shutts@metoffice.gov.uk ABSTRACT A form of kinetic energy backscatter (CASBS)
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 informationContents. Parti Fundamentals. 1. Introduction. 2. The Coriolis Force. Preface Preface of the First Edition
Foreword Preface Preface of the First Edition xiii xv xvii Parti Fundamentals 1. Introduction 1.1 Objective 3 1.2 Importance of Geophysical Fluid Dynamics 4 1.3 Distinguishing Attributes of Geophysical
More informationAtmospheric Dynamics
Atmospheric Dynamics the CAM-FV dynamical core Peter Hjort Lauritzen Atmospheric Modeling and Predictability Section Climate and Global Dynamics Division National Center for Atmospheric Research Community
More informationThe next-generation supercomputer and NWP system of the JMA
The next-generation supercomputer and NWP system of the JMA Masami NARITA m_narita@naps.kishou.go.jp Numerical Prediction Division (NPD), Japan Meteorological Agency (JMA) Purpose of supercomputer & NWP
More informationEddy-resolving Simulation of the World Ocean Circulation by using MOM3-based OGCM Code (OFES) Optimized for the Earth Simulator
Chapter 1 Atmospheric and Oceanic Simulation Eddy-resolving Simulation of the World Ocean Circulation by using MOM3-based OGCM Code (OFES) Optimized for the Earth Simulator Group Representative Hideharu
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 informationNext Genera*on Compu*ng: Needs and Opportuni*es for Weather, Climate, and Atmospheric Sciences. David Randall
Next Genera*on Compu*ng: Needs and Opportuni*es for Weather, Climate, and Atmospheric Sciences David Randall Way back I first modified, ran, and analyzed results from an atmospheric GCM in 1972. The model
More informationAssessing and Improving the Numerical Solution of Atmospheric Physics in an Earth System Model
Funding for the studies presented here was provided by US DOE s Office of Biological and Environmental Research and Office of Advanced Scientific Computing Research Assessing and Improving the Numerical
More informationThe momentum budget of the QBO in WACCM. Rolando Garcia and Yaga Richter National Center for Atmospheric Research Boulder, CO
The momentum budget of the QBO in WACCM Rolando Garcia and Yaga Richter National Center for Atmospheric Research Boulder, CO a few model details Simulations use WACCM 5 at 1 horizontal resolution and 110
More informationBALANCED FLOW: EXAMPLES (PHH lecture 3) Potential Vorticity in the real atmosphere. Potential temperature θ. Rossby Ertel potential vorticity
BALANCED FLOW: EXAMPLES (PHH lecture 3) Potential Vorticity in the real atmosphere Need to introduce a new measure of the buoyancy Potential temperature θ In a compressible fluid, the relevant measure
More informationA Discussion on The Effect of Mesh Resolution on Convective Boundary Layer Statistics and Structures Generated by Large-Eddy Simulation by Sullivan
耶鲁 - 南京信息工程大学大气环境中心 Yale-NUIST Center on Atmospheric Environment A Discussion on The Effect of Mesh Resolution on Convective Boundary Layer Statistics and Structures Generated by Large-Eddy Simulation
More informationDynamical coupling between the middle atmosphere and lower thermosphere
Dynamical coupling between the middle atmosphere and lower thermosphere Anne Smith, Dan Marsh, Nick Pedatella NCAR* Tomoko Matsuo CIRES/NOAA NCAR is sponsored by the National Science Foundation Model runs
More informationThe effect of varying forcing on the transport of heat by transient eddies.
The effect of varying forcing on the transport of heat by transient eddies. LINDA MUDONI Department of Atmospheric Sciences, Iowa State University May 02 20 1. Introduction The major transport of heat
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 informationMass-conserving and positive-definite semi-lagrangian advection in NCEP GFS: Decomposition for massively parallel computing without halo
Mass-conserving and positive-definite semi-lagrangian advection in NCEP GFS: Decomposition for massively parallel computing without halo Hann-Ming Henry Juang Environmental Modeling Center, NCEP NWS,NOAA,DOC,USA
More informationSuper-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models PI Isaac Ginis Graduate School of Oceanography
More informationDiabatic processes and the structure of the warm conveyor belt
2 nd European Windstorm Workshop Leeds, 3-4 September 2012 Diabatic processes and the structure of the warm conveyor belt Oscar Martínez-Alvarado J. Chagnon, S. Gray, R. Plant, J. Methven Department of
More informationAtm S 547 Boundary Layer Meteorology
Lecture 8. Parameterization of BL Turbulence I In this lecture Fundamental challenges and grid resolution constraints for BL parameterization Turbulence closure (e. g. first-order closure and TKE) parameterizations
More informationImperial College London
Solar Influence on Stratosphere-Troposphere Dynamical Coupling Isla Simpson, Joanna D. Haigh, Space and Atmospheric Physics, Imperial College London Mike Blackburn, Department of Meteorology, University
More information2. Conservation laws and basic equations
2. Conservation laws and basic equations Equatorial region is mapped well by cylindrical (Mercator) projection: eastward, northward, upward (local Cartesian) coordinates:,, velocity vector:,,,, material
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 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 informationOverview of the Numerics of the ECMWF. Atmospheric Forecast Model
Overview of the Numerics of the Atmospheric Forecast Model M. Hortal Seminar 6 Sept 2004 Slide 1 Characteristics of the model Hydrostatic shallow-atmosphere approimation Pressure-based hybrid vertical
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 informationThe Fifth-Generation NCAR / Penn State Mesoscale Model (MM5) Mark Decker Feiqin Xie ATMO 595E November 23, 2004 Department of Atmospheric Science
The Fifth-Generation NCAR / Penn State Mesoscale Model (MM5) Mark Decker Feiqin Xie ATMO 595E November 23, 2004 Department of Atmospheric Science Outline Basic Dynamical Equations Numerical Methods Initialization
More informationGeneralizing the Boussinesq Approximation to Stratified Compressible Flow
Generalizing the Boussinesq Approximation to Stratified Compressible Flow Dale R. Durran a Akio Arakawa b a University of Washington, Seattle, USA b University of California, Los Angeles, USA Abstract
More informationTorres del Paine, Patagonia, Chile
Torres del Paine, Patagonia, Chile Alan Rhoades 1, Paul Ullrich 1,3, Colin Zarzycki 2, Hans Johansen 1, Zexuan Xu 1, and William Collins 1 Lawrence Berkeley National Laboratory 1, National Center for Atmospheric
More information