Model description of AGCM5 of GFD-Dennou-Club edition. SWAMP project, GFD-Dennou-Club
|
|
- Paulina Cameron
- 5 years ago
- Views:
Transcription
1 Model description of AGCM5 of GFD-Dennou-Club edition SWAMP project, GFD-Dennou-Club Mar 01, 2006
2 AGCM5 of the GFD-DENNOU CLUB edition is a three-dimensional primitive system on a sphere (Swamp Project, 1998). The model is basically the same as that used in Numaguti (1993), and its details are described in Numaguti(1992). In the following subsections, we briefly describe the model system and its implementation. 1 Basic equations The dynamical part of the model consists of the vorticity, divergence, hydrostatic, humidity and thermodynamic equations as follows. dζ dt = 1 u { (ζ + f)u σ a cos ϕ λ σ RT π a cos ϕ λ } 1 u [{(ζ + f)v σ a cos ϕ ϕ σ RT π diff } cos ϕ] Fζ, (1) a ϕ dd dt = 1 u {(ζ + f)v σ a cos ϕ λ σ xrt π a ϕ } + 1 u [{(ζ + f)v σ a cos ϕ ϕ σ RT π } cos ϕ] a ϕ 2 (Φ + RT π + E) F diff D, (2) dπ dt = 1 u a cos ϕ λ 1 σ (v cos ϕ) a cos ϕ ϕ σ, (3) Φ σ = RT σ, (4) dq dt = g Fq vdf p s σ + F q vdf + Sq cond, (5) dt dt = RT { π c p t + u π a cos ϕ λ + v π a ϕ + σ } σ + 1 ( g F vdf T c p p s σ + g F vdf ) rad + F diff T + LSq cond. (6) p s σ In these equations, d dt ζ t + u a cos ϕ λ + v a ϕ + σ σ, (7) 1 v a cos ϕ λ 1 (u cos ϕ), (8) a cos ϕ ϕ 1
3 D 1 u a cos ϕ λ + 1 (v cos ϕ), a cos ϕ ϕ (9) Φ gz, (10) π ln p s, (11) where, (λ, ϕ) are the latitudinal and longitudinal coordinates; σ is the vertical coordinate; u, v are the horizontal wind components; σ is the vertical wind in σ coordinate; T is temperature; q is specific humidity; p s is surface pressure; T g is surface temperature; Φ is geopotential; f is the Coriolis parameter; a is the radius of the earth; g is the acceleration of gravity; R is the gas constant of dry air; c p is the atmospheric specific heat at constant pressure; L is the Fq diff are the horizontal diffusion terms; Fu vdf, Fv vdf, F vdf T, Fq vdf are the vertical diffusion terms; and Sq cond is the condensation source (loss) of the specific humidity. The horizontal diffusion terms are represented by a hyper-viscosity type formula, F diff = ( 1) n K 2n (12) latent heat of water; F diff ζ, F diff D, F diff T for which n = 8 is adopted. The value of K is determined such that the e-folding time of the component with wavenumber 42 is four hours. 2 Discretization The equations are discretized horizontally by a non-alias latitude-longitude grid and are integrated by the pseudo-spectral method with a triangular truncation. The scheme of Arakawa and Suarez (1983) is utilized for discretizing in the vertical direction. The time integration is performed by the use of a leap-frog scheme for the terms of dynamical processes other than dissipative processes, a backward scheme for the terms of physical processes, and an Euler scheme for the horizontal diffusion terms. The time filter of Asselin (1972) is used to remove the computational mode caused by the leap-frog scheme. 3 Physical processes The model contains the following simplified physical processes. 2
4 3.1 Cumulus parameterization The cumulus parameterization schemes implemented are the adjustment scheme (Manabe et al., 1965) and Kuo scheme (Kuo, 1974). In the followings, the description of adjustment schmes is given. Condensed water is removed from the system immediately after condensation and hence there are neither rain drops nor clouds. Moist convective adjustment occurs when the following conditions are satisfied: In two neighboring levels, the atmosphere becomes moist unstable, that is, the temperature lapse rate exceeds the moist adiabatic lapse rate. Both of the neighboring levels are saturated or oversaturated. In the process of adjustment, the temperature and specific humidity of the two neighboring levels are instantaneously adjusted to values so that the levels are saturated and the temperature lapse rate between them coincides with the moist adiabatic lapse rate. In due course, moist static energy has to be conserved: k l=k 1 (C p T l + Lq l ) p l = k l=k 1 (C p ˆTl + Lˆq l ) p l, (13) where k is the index of the vertical level. Here, ˆT, ˆq denote the values before adjustment, while T, q denote the values after the adjustment. The moist unstable condition for layer (k, k 1) is evaluated by ST k 1/2 ˆT k ˆT k 1 + L ( q ( C ˆT k, p k ) q ( ˆT k 1, p k 1 ) ) R p k 1/2 p C p p k 1/2 > 0. ˆTk + ˆT k 1 2 (14) The saturation condition of two neighboring levels is simply expressed by ˆq k 1 > q ( ˆT k 1, p k 1 ), ˆq k > q ( ˆT k, p k ). (15) The adjusted values (T k, q k ) are determined as T k 1 = ˆT k 1 + T k 1, (16) T k = ˆT k + T k, (17) 3
5 q k 1 = q ( ˆT k 1, p k 1 ) + q T T k 1, (18) k 1 q k = q ( ˆT k, p k ) + q T T k, (19) k T k = where (1 + γ k 1 ) p k 1 S k 1/2 + [1 + γ k 1 κ p k 1/2 2p k 1/2 ] L C p ˆQ (1 + γ k 1 )(1 + γ k )( p k 1 + p k ) + κ p k 1/2 2p k 1/2 [(1 + γ k 1 ) p k 1 (1 + γ k ) p k ], (20) T k 1 = (1 + γ k) p k 1 L T k + T k, (21) (1 + γ k 1 ) p k 1 (1 + γ k 1 ) p k 1 C p γ k = L C p q T, γ k 1 = L q, (22) k C p T k 1 ˆQ = (ˆq k 1 q ( ˆT k 1, p k 1 )) p k 1 + (ˆq k q ( ˆT k, p k )) p k. (23) Precipitation of one column is calculated as K P = (ˆq k q k ) p k k=1 g. (24) These adjustment procedures are applied to all levels until the vertical temperature and moisture profiles converge. 3.2 Large scale condensation The large scale condensation scheme is applied at each grid point when the value of the specific humidity exceeds its saturation value. The large scale condensation scheme demands that the value of the specific humidity be adjusted to the saturation value. Condensed water is assumed to fall out as rain and is removed from the system immediately. There are no rain drops or clouds, as in the case of the convective adjustment scheme. In the time integration procedure, the large scale condensation scheme is applied just after the convective adjustment scheme. The condition for large scale condensation can be denoted by ˆq > q ( ˆT, p), (25) 4
6 where ˆT, ˆq are the values of temperature and specific humidity before the adjustment. After the large scale condensation adjustment, the values of the temperature and specific humidity become and q = q (T, p) (26) C p T + Lq = C p ˆT + Lˆq. (27) The solutions to these are obtained by Newton method; T = ˆT + L ˆq q ( ˆT, p) C p 1 + L q, (28) C p T q = ˆq L C p q T The amount of precipitation is calculated by 3.3 Radiation ˆq q ( ˆT, p) 1 + L q. (29) C p T P = p (ˆq q). (30) g Radiative flux is calclated by band model. The effects of scattering of solar radiation is not considered. It is assumed that there are no cloud nor rain drops. As for longwave radiation, we include only the absorption and emission of water and dry air; the effect of scattering is not included. The absorption coefficients of the longwave radiation are chosen so that the cooling profile of the atmosphere roughly resembles the observed one. Radiation flux F R,i is F R (z) = (πb(t g ) πb(t s )) T f (z, 0)+πB(T (z T ))T f zt (z, z T ) 0 dπb dξ T f (z, ξ)dξ, (31) where T f (z 1, z 2 ) is the flux transfer function between z = z 1, z 2 and πb σ SB T 4 is the source function. The flux transfer function T f (z 1, z 2 ) is assumed to be given as N R T f (z 1, z 2 ) = b i exp ( δ R τ R,i (z 1 ) τ R,i (z 2 ) ). (32) i=1 5
7 Here, τ i (z) is the optical depth whose value at the top of the atmosphere is 0 and τ R,i (z) = k R,i ρqdz + k R,i ρdz, (33) z where k R,i and k R,i are the absorption coefficients of water vapor and dry air for wavenumber band i, respectively. The values of these absorption coefficients are assumed to be constant. b i is the ratio of the energy for wavenumber band i to the total radiative energy and is also assumed to be constant. The default value for δ R is 1.5. In most cases, N R is set to be 4, three representing water vapor bands and one representing dry air band. The default values of k R,i, k R,i and b i are summarized in Table 1. z Table 1: The values of constants used in the longwave radiation scheme. Band Number k R,i k R,i b i Vertical mixing Vertical turbulent mixing is represented by the level II scheme of Mellor and Yamada (1974). In this scheme of Mellor and Yamada (1974), diffusion coefficients K M, K T, and K q depends on stability and are given as functions of u, v, and T : K M = l 2 v z S M, (34) K T = K q = l 2 v z S H. (35) Following Blakadar(1962), the mixing length l is evaluated by l = where l 0 = 300 m and κ = 0.4 is the Kárman constant. κz 1 + κz/l 0, (36) 6
8 In these expressions, S M and S H are given as functions of the Richardson number. S M = B 1/2 1 (1 R f ) 1/2 S M 3/2, (37) S H = B 1/2 1 (1 R f ) 1/2 S M 1/2 S H, (38) where S H = α 1 α 2 R f 1 R f, (39) R f is the flux Richardson number, given as S M = β 1 β 2 R f β 3 β 4 R f S H. (40) R f = 1 [ ] β 1 + β 4 R ib (β 1 + β 4 R ib ) 2β 2 4β 2 β 3 R ib, (41) 2 where R ib is the bulk Richardson number and R ib is defined as α 1 = 3A 2 γ 1, (42) α 2 = 3A 2 (γ 1 + γ 2 ), (43) β 1 = A 1 B 1 (γ 1 C 1 ), (44) β 2 = A 1 [B 1 (γ 1 C 1 ) + 6A 1 + 3A 2 ], (45) β 3 = A 2 B 1, (46) β 4 = A 2 [B 1 (γ 1 + γ 2 ) 3A 1 ], (47) (A 1, B 1, A 2, B 2, C 1 ) = (0.92, 16.6, 0.74, 10.1, 0.08), (48) γ 1 = 1 3 2A 1 B 1, γ 2 = B 2 B A 1 B 1. (49) R ib = ( u z q θ θ s z ) 2 + ( v z ) 2. (50) Here, the value of ( ) at the vertical level k 1/2 is evaluated as k k 1. 7
9 3.5 Surface flux The surface fluxes of momentum, heat and water vapor are evaluated by the bulk formulae (F λ ) 1/2 = ρ s C M v 1 u 1, (51) (F ϕ ) 1/2 = ρ s C M v 1 u 1, (52) (F T ) 1/2 = ρ s C p C H v 1 (T g ( T Π 1)), (53) (F q ) 1/2 = ρ s βc E v 1 (q (T g, p s ) q 1 ), (54) where the suffix 1/2 indicates that these are the value at the vertical level 1/2. Following Louis et al. (1982), the bulk coefficients C M, C H, and C E are determined by ( C M RiB / ) R ib R ib 0, C M = ( C M0 [1 10R ib 1 + z1 75CM0 z M0 R ib ) 1 ] R ib < 0, (55) ( C H RiB / ) R ib R ib 0, C H = C E = ( C H0 [1 15R ib 1 + z1 75CH0 z M0 R ib ) 1 ] R ib < 0, (56) where R ib is the bulk Richardson number evaluated by R ib = g θ s (θ 1 θ g )/g (u 1 /z 1 ) 2 = g θ s T 1 (p s /p 1 ) κ T g u 2 1/z 1. (57) C M0 and C H0 are the bulk coefficients for neutral case: κ C M0 = ln ( ), (58) z 1 z M0 κ C H0 = ln ( ), (59) z 1 z H0 where z M0 and z H0 are the roughness parameters of the surface. The values are 10 4 m. 8
10 References Arakawa, A. and M.J. Suarez, 1983: Vertical differencing of the primitive equations in sigma coordinates. J. Atmos. Sci., 111, Asselin, R. 1972: Frequency filter for time integrations. Mon. Wet. Rev., 100, Blakadar, A.K., 1962: The vertical distribution of wind and turbulent exchange in neutral atmosphere. J. Geophys. Res., 67, Kuo, H. L., 1974: Further studies of the parameterization of the influence of cumulus convection on large-scale flow. J. Atmos. Sci., 31, Louis, J.-F. and M. Tiedtke and J.-F. Geleyn, 1982: A short history of the PBL parameterization at ECMWF. Workshop on Planetary Boundary layer Parameterization, 59-80, ECMWF, Reading U.K. Manabe, S., J. Smagorinsky and R.F. Strickler, 1965: Simulated climatology of a general circulation model with a hydrologic cycle. Mon. Wea. Rev., 93, Mellor, G.L. and T. Yamada, 1974: A hierarchy of turbulence closure models for the planetary boundary layers. J. Atmos. Sci., 31, Numaguti, A., 1992: Numerical experiments on the large scale structure of cumulus activity in the tropics. Ph.D. Thesis, University of Tokyo. (in Japanese). Numaguti, A., 1993: Dynamics and energy balance of the Hadley circulation and the tropical precipitation zones: Significance of the distribution of evaporation. J. Atmos. Sci., 50, Swamp Project (1998) AGCM5. GFD Dennou Club. 9
5. General Circulation Models
5. General Circulation Models I. 3-D Climate Models (General Circulation Models) To include the full three-dimensional aspect of climate, including the calculation of the dynamical transports, requires
More information2.1 Effects of a cumulus ensemble upon the large scale temperature and moisture fields by induced subsidence and detrainment
Atmospheric Sciences 6150 Cloud System Modeling 2.1 Effects of a cumulus ensemble upon the large scale temperature and moisture fields by induced subsidence and detrainment Arakawa (1969, 1972), W. Gray
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 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 informationRadiative equilibrium Some thermodynamics review Radiative-convective equilibrium. Goal: Develop a 1D description of the [tropical] atmosphere
Radiative equilibrium Some thermodynamics review Radiative-convective equilibrium Goal: Develop a 1D description of the [tropical] atmosphere Vertical temperature profile Total atmospheric mass: ~5.15x10
More informationRadiation in climate models.
Lecture. Radiation in climate models. Objectives:. A hierarchy of the climate models.. Radiative and radiative-convective equilibrium.. Examples of simple energy balance models.. Radiation in the atmospheric
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 informationPALM - Cloud Physics. Contents. PALM group. last update: Monday 21 st September, 2015
PALM - Cloud Physics PALM group Institute of Meteorology and Climatology, Leibniz Universität Hannover last update: Monday 21 st September, 2015 PALM group PALM Seminar 1 / 16 Contents Motivation Approach
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 informationA Numerical Study on Appearance of the Runaway Greenhouse State of a Three-Dimensional Gray Atmosphere
15 NOVEMBER 2002 ISHIWATARI ET AL. 3223 A Numerical Study on Appearance of the Runaway Greenhouse State of a Three-Dimensional Gray Atmosphere MASAKI ISHIWATARI Division of Ocean and Atmospheric Sciences,
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 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 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 informationESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr. DeCaria
ESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr DeCaria References: An Introduction to Dynamic Meteorology, Holton MOMENTUM EQUATIONS The momentum equations governing the ocean or atmosphere
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 informationMeteorology 6150 Cloud System Modeling
Meteorology 6150 Cloud System Modeling Steve Krueger Spring 2009 1 Fundamental Equations 1.1 The Basic Equations 1.1.1 Equation of motion The movement of air in the atmosphere is governed by Newton s Second
More informationAtmosphere, Ocean and Climate Dynamics Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 12.003 Atmosphere, Ocean and Climate Dynamics Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Contents
More informationAn Introduction to Climate Modeling
An Introduction to Climate Modeling A. Gettelman & J. J. Hack National Center for Atmospheric Research Boulder, Colorado USA Outline What is Climate & why do we care Hierarchy of atmospheric modeling strategies
More informationECMWF Overview. The European Centre for Medium-Range Weather Forecasts is an international. organisation supported by 23 European States.
ECMWF Overview The European Centre for Medium-Range Weather Forecasts is an international organisation supported by 3 European States. The center was established in 1973 by a Convention and the real-time
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 informationAtmospheric Sciences 321. Science of Climate. Lecture 13: Surface Energy Balance Chapter 4
Atmospheric Sciences 321 Science of Climate Lecture 13: Surface Energy Balance Chapter 4 Community Business Check the assignments HW #4 due Wednesday Quiz #2 Wednesday Mid Term is Wednesday May 6 Practice
More informationIMPACT OF SOIL FREEZING ON THE CONTINENTAL-SCALE SEASONAL CYCLE SIMULATED BY A GENERAL CIRCULATION MODEL
IMPACT OF SOIL FREEZING ON THE CONTINENTAL-SCALE SEASONAL CYCLE SIMULATED BY A GENERAL CIRCULATION MODEL Kumiko Takata 1, Masahide Kimoto 2 1. Domestic Research Fellow, National Institute of Environmental
More informationThermodynamics Review [?] Entropy & thermodynamic potentials Hydrostatic equilibrium & buoyancy Stability [dry & moist adiabatic]
Thermodynamics Review [?] Entropy & thermodynamic potentials Hydrostatic equilibrium & buoyancy Stability [dry & moist adiabatic] Entropy 1. (Thermodynamics) a thermodynamic quantity that changes in a
More informationONE DIMENSIONAL CLIMATE MODEL
JORGE A. RAMÍREZ Associate Professor Water Resources, Hydrologic and Environmental Sciences Civil Wngineering Department Fort Collins, CO 80523-1372 Phone: (970 491-7621 FAX: (970 491-7727 e-mail: Jorge.Ramirez@ColoState.edu
More informationAPPENDIX B. The primitive equations
APPENDIX B The primitive equations The physical and mathematical basis of all methods of dynamical atmospheric prediction lies in the principles of conservation of momentum, mass, and energy. Applied to
More informationRadiative-Convective Models. The Hydrological Cycle Hadley Circulation. Manabe and Strickler (1964) Course Notes chapter 5.1
Climate Modeling Lecture 8 Radiative-Convective Models Manabe and Strickler (1964) Course Notes chapter 5.1 The Hydrological Cycle Hadley Circulation Prepare for Mid-Term (Friday 9 am) Review Course Notes
More informationWaVaCS summerschool Autumn 2009 Cargese, Corsica
Introduction Part I WaVaCS summerschool Autumn 2009 Cargese, Corsica Holger Tost Max Planck Institute for Chemistry, Mainz, Germany Introduction Overview What is a parameterisation and why using it? Fundamentals
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 informationToday s Lecture: Atmosphere finish primitive equations, mostly thermodynamics
Today s Lecture: Atmosphere finish primitive equations, mostly thermodynamics Reference Peixoto and Oort, Sec. 3.1, 3.2, 3.4, 3.5 (but skip the discussion of oceans until next week); Ch. 10 Thermodynamic
More informationBoundary layer equilibrium [2005] over tropical oceans
Boundary layer equilibrium [2005] over tropical oceans Alan K. Betts [akbetts@aol.com] Based on: Betts, A.K., 1997: Trade Cumulus: Observations and Modeling. Chapter 4 (pp 99-126) in The Physics and Parameterization
More informationSome remarks on climate modeling
Some remarks on climate modeling A. Gettelman & J. J. Hack National Center for Atmospheric Research Boulder, Colorado USA Selected overheads by Doug Nychka Outline Hierarchy of atmospheric modeling strategies
More informationWeather Research and Forecasting Model. Melissa Goering Glen Sampson ATMO 595E November 18, 2004
Weather Research and Forecasting Model Melissa Goering Glen Sampson ATMO 595E November 18, 2004 Outline What does WRF model do? WRF Standard Initialization WRF Dynamics Conservation Equations Grid staggering
More information1. The vertical structure of the atmosphere. Temperature profile.
Lecture 4. The structure of the atmosphere. Air in motion. Objectives: 1. The vertical structure of the atmosphere. Temperature profile. 2. Temperature in the lower atmosphere: dry adiabatic lapse rate.
More informationMEA 716 Exercise, BMJ CP Scheme With acknowledgements to B. Rozumalski, M. Baldwin, and J. Kain Optional Review Assignment, distributed Th 2/18/2016
MEA 716 Exercise, BMJ CP Scheme With acknowledgements to B. Rozumalski, M. Baldwin, and J. Kain Optional Review Assignment, distributed Th 2/18/2016 We have reviewed the reasons why NWP models need to
More informationIncorporation of 3D Shortwave Radiative Effects within the Weather Research and Forecasting Model
Incorporation of 3D Shortwave Radiative Effects within the Weather Research and Forecasting Model W. O Hirok and P. Ricchiazzi Institute for Computational Earth System Science University of California
More informationThermodynamics of Atmospheres and Oceans
Thermodynamics of Atmospheres and Oceans Judith A. Curry and Peter J. Webster PROGRAM IN ATMOSPHERIC AND OCEANIC SCIENCES DEPARTMENT OF AEROSPACE ENGINEERING UNIVERSITY OF COLORADO BOULDER, COLORADO USA
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 informationMore on Diabatic Processes
More on Diabatic Processes Chapter 3 Write Qtotal = Qrad + Qcond + Qsen total heating radiative heating condensationa l heating sensible heating While diabatic processes drive atmospheric motions, the
More informationA brief overview of the scheme is given below, taken from the whole description available in Lopez (2002).
Towards an operational implementation of Lopez s prognostic large scale cloud and precipitation scheme in ARPEGE/ALADIN NWP models F.Bouyssel, Y.Bouteloup, P. Marquet Météo-France, CNRM/GMAP, 42 av. G.
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 information11 days (00, 12 UTC) 132 hours (06, 18 UTC) One unperturbed control forecast and 26 perturbed ensemble members. --
APPENDIX 2.2.6. CHARACTERISTICS OF GLOBAL EPS 1. Ensemble system Ensemble (version) Global EPS (GEPS1701) Date of implementation 19 January 2017 2. EPS configuration Model (version) Global Spectral Model
More informationWQMAP (Water Quality Mapping and Analysis Program) is a proprietary. modeling system developed by Applied Science Associates, Inc.
Appendix A. ASA s WQMAP WQMAP (Water Quality Mapping and Analysis Program) is a proprietary modeling system developed by Applied Science Associates, Inc. and the University of Rhode Island for water quality
More informationSensitivity of Precipitation in Aqua-Planet Experiments with an AGCM
ATMOSPHERIC AND OCEANIC SCIENCE LETTERS, 2014, VOL. 7, NO. 1, 1 6 Sensitivity of Precipitation in Aqua-Planet Experiments with an AGCM YU Hai-Yang 1,2, BAO Qing 1, ZHOU Lin-Jiong 1,2, WANG Xiao-Cong 1,
More informationAOSC 614. Chapter 4 Junjie Liu Oct. 26, 2006
AOSC 614 Chapter 4 Junjie Liu Oct. 26, 2006 The role of parameterization in the numerical model Why do we need parameterization process? Chapter 2 discussed the governing equations. Chapter 3 discussed
More informationThe Effect of Sea Spray on Tropical Cyclone Intensity
The Effect of Sea Spray on Tropical Cyclone Intensity Jeffrey S. Gall, Young Kwon, and William Frank The Pennsylvania State University University Park, Pennsylvania 16802 1. Introduction Under high-wind
More informationClouds and turbulent moist convection
Clouds and turbulent moist convection Lecture 2: Cloud formation and Physics Caroline Muller Les Houches summer school Lectures Outline : Cloud fundamentals - global distribution, types, visualization
More informationImpact of different cumulus parameterizations on the numerical simulation of rain over southern China
Impact of different cumulus parameterizations on the numerical simulation of rain over southern China P.W. Chan * Hong Kong Observatory, Hong Kong, China 1. INTRODUCTION Convective rain occurs over southern
More informationThe atmospheric boundary layer: Where the atmosphere meets the surface. The atmospheric boundary layer:
The atmospheric boundary layer: Utrecht Summer School on Physics of the Climate System Carleen Tijm-Reijmer IMAU The atmospheric boundary layer: Where the atmosphere meets the surface Photo: Mark Wolvenne:
More informationAdiabatic expansion Isothermal compression Adiabatic compression
Tropical Cyclones: Steady State Physics 1 Energy Production 2 Carnot Theorem: Maximum efficiency results from a pa rticular energy e cycle: Isothermal expansion Adiabatic expansion Isothermal compression
More informationLecture 7: The Monash Simple Climate
Climate of the Ocean Lecture 7: The Monash Simple Climate Model Dr. Claudia Frauen Leibniz Institute for Baltic Sea Research Warnemünde (IOW) claudia.frauen@io-warnemuende.de Outline: Motivation The GREB
More informationAnalysis of Cloud-Radiation Interactions Using ARM Observations and a Single-Column Model
Analysis of Cloud-Radiation Interactions Using ARM Observations and a Single-Column Model S. F. Iacobellis, R. C. J. Somerville, D. E. Lane, and J. Berque Scripps Institution of Oceanography University
More informationThe Tropical Atmosphere: Hurricane Incubator
The Tropical Atmosphere: Hurricane Incubator Images from journals published by the American Meteorological Society are copyright AMS and used with permission. A One-Dimensional Description of the Tropical
More informationProject 3 Convection and Atmospheric Thermodynamics
12.818 Project 3 Convection and Atmospheric Thermodynamics Lodovica Illari 1 Background The Earth is bathed in radiation from the Sun whose intensity peaks in the visible. In order to maintain energy balance
More informationGoverning Equations and Scaling in the Tropics
Governing Equations and Scaling in the Tropics M 1 ( ) e R ε er Tropical v Midlatitude Meteorology Why is the general circulation and synoptic weather systems in the tropics different to the those in the
More informationREGIONAL CLIMATE SIMULATION OVER ASIA UNDER THE GLOBAL WARMING NESTED IN THE CCSR/NIES AGCM
REGIONAL CLIMATE SIMULATION OVER ASIA UNDER THE GLOBAL WARMING NESTED IN THE CCSR/NIES AGCM KOJI DAIRAKU Atmospheric Environment Division, National Institute for Environmental Studies, 16-2 Onogawa Tsukuba,
More informationLecture 10a: The Hadley Cell
Lecture 10a: The Hadley Cell Geoff Vallis; notes by Jim Thomas and Geoff J. Stanley June 27 In this short lecture we take a look at the general circulation of the atmosphere, and in particular the Hadley
More information1 Introduction to Governing Equations 2 1a Methodology... 2
Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................
More informationSimple Climate Models
Simple Climate Models Lecture 3 One-dimensional (vertical) radiative-convective models Vertical atmospheric processes The vertical is the second important dimension, because there are... Stron radients
More informationRadiation, Sensible Heat Flux and Evapotranspiration
Radiation, Sensible Heat Flux and Evapotranspiration Climatological and hydrological field work Figure 1: Estimate of the Earth s annual and global mean energy balance. Over the long term, the incoming
More informationLecture 9: Climate Sensitivity and Feedback Mechanisms
Lecture 9: Climate Sensitivity and Feedback Mechanisms Basic radiative feedbacks (Plank, Water Vapor, Lapse-Rate Feedbacks) Ice albedo & Vegetation-Climate feedback Cloud feedback Biogeochemical feedbacks
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 information282 Journal of the Meteorological Society of Japan Vol. 60, No. 1. A Theory and Method of Long-Range Numerical
282 Journal of the Meteorological Society of Japan Vol. 60, No. 1 A Theory and Method of Long-Range Numerical Weather Forecasts By Chao Jih-Ping, Guo Yu-Fu and Xin Ru-Nan Institute of Atmospheric Physics,
More informationDevelopment of a Coupled Atmosphere-Ocean-Land General Circulation Model (GCM) at the Frontier Research Center for Global Change
Chapter 1 Atmospheric and Oceanic Simulation Development of a Coupled Atmosphere-Ocean-Land General Circulation Model (GCM) at the Frontier Research Center for Global Change Project Representative Tatsushi
More informationUnified Cloud and Mixing Parameterizations of the Marine Boundary Layer: EDMF and PDF-based cloud approaches
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Unified Cloud and Mixing Parameterizations of the Marine Boundary Layer: EDMF and PDF-based cloud approaches LONG-TERM
More informationA B C D PROBLEMS Dilution of power plant plumes. z z z z
69 PROBLEMS 4. Dilution of power plant plumes Match each power plant plume (-4) to the corresponding atmospheric lapse rate (A-D, solid lines; the dashed line is the adiabatic lapse rate Γ). Briefly comment
More informationA Gray-Radiation Aquaplanet Moist GCM. Part II: Energy Transports in Altered Climates
1680 J O U R N A L O F T H E A T M O S P H E R I C S C I E N C E S VOLUME 64 A Gray-Radiation Aquaplanet Moist GCM. Part II: Energy Transports in Altered Climates DARGAN M. W. FRIERSON University of Chicago,
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 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 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 informationCourse , General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 4: Water Vapor Budget
Course 12.812, General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 4: Water Vapor Budget Water Vapor Distribution First let us look at the distribution of specific humidity, q. The
More informationHurricanes are intense vortical (rotational) storms that develop over the tropical oceans in regions of very warm surface water.
Hurricanes: Observations and Dynamics Houze Section 10.1. Holton Section 9.7. Emanuel, K. A., 1988: Toward a general theory of hurricanes. American Scientist, 76, 371-379 (web link). http://ww2010.atmos.uiuc.edu/(gh)/guides/mtr/hurr/home.rxml
More informationHabitable Planets. Much of it stolen from. Yutaka ABE University of Tokyo
Habitable Planets Much of it stolen from Yutaka ABE University of Tokyo 1. Habitability and Water Why water? Importance of Liquid Gas: highly mobile, but low material density. Solid: high density but very
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More information2σ e s (r,t) = e s (T)exp( rr v ρ l T ) = exp( ) 2σ R v ρ l Tln(e/e s (T)) e s (f H2 O,r,T) = f H2 O
Formulas/Constants, Physics/Oceanography 4510/5510 B Atmospheric Physics II N A = 6.02 10 23 molecules/mole (Avogadro s number) 1 mb = 100 Pa 1 Pa = 1 N/m 2 Γ d = 9.8 o C/km (dry adiabatic lapse rate)
More informationChapter 4 Water Vapor
Chapter 4 Water Vapor Chapter overview: Phases of water Vapor pressure at saturation Moisture variables o Mixing ratio, specific humidity, relative humidity, dew point temperature o Absolute vs. relative
More informationEST 5101 Climate Change Science Chapter 5. Climate Models. Rezaul Karim Environmental Science & Technology Jessore University of science & Technology
EST 5101 Climate Change Science Chapter 5. Climate Models Rezaul Karim Environmental Science & Technology Jessore University of science & Technology Chapter outline overview of climate modelling; How Are
More informationThe Ocean-Atmosphere System II: Oceanic Heat Budget
The Ocean-Atmosphere System II: Oceanic Heat Budget C. Chen General Physical Oceanography MAR 555 School for Marine Sciences and Technology Umass-Dartmouth MAR 555 Lecture 2: The Oceanic Heat Budget Q
More informationEarth s Energy Balance and the Atmosphere
Earth s Energy Balance and the Atmosphere Topics we ll cover: Atmospheric composition greenhouse gases Vertical structure and radiative balance pressure, temperature Global circulation and horizontal energy
More informationDepartment of Meteorology, School of Ocean and Earth Science and Technology, University of Hawaii at Manoa, Honolulu, Hawaii
478 J O U R N A L O F C L I M A T E VOLUME 0 Horizontal and Vertical Structures of the Northward-Propagating Intraseasonal Oscillation in the South Asian Monsoon Region Simulated by an Intermediate Model*
More informationBrief Overview of the Global Atmosphere
Brief Overview of the Global Atmosphere Images on pages 3, 5 through 9, 19, 77 through 80 and 90 are copyrighted by Academic Press, NY, 1999. The images are in a chapter entitled "Quasi-equilibrium thinking"
More informationLecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2)
Lecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2) The ABL, though turbulent, is not homogeneous, and a critical role of turbulence is transport and mixing of air properties, especially in the
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 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 informationJRA-55 Product Users' Handbook
Climate Prediction Division Global Environment and Marine Department Japan Meteorological Agency September 2013 Change record Version Date Remarks 1.0 30 September 2013 First version 2.0 3 March 2014 Corrected
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 informationLecture 10: Climate Sensitivity and Feedback
Lecture 10: Climate Sensitivity and Feedback Human Activities Climate Sensitivity Climate Feedback 1 Climate Sensitivity and Feedback (from Earth s Climate: Past and Future) 2 Definition and Mathematic
More informationP1M.4 COUPLED ATMOSPHERE, LAND-SURFACE, HYDROLOGY, OCEAN-WAVE, AND OCEAN-CURRENT MODELS FOR MESOSCALE WATER AND ENERGY CIRCULATIONS
P1M.4 COUPLED ATMOSPHERE, LAND-SURFACE, HYDROLOGY, OCEAN-WAVE, AND OCEAN-CURRENT MODELS FOR MESOSCALE WATER AND ENERGY CIRCULATIONS Haruyasu NAGAI *, Takuya KOBAYASHI, Katsunori TSUDUKI, and Kyeongok KIM
More informationParameterization of effects of unresolved clouds and precipitation
Parameterization of effects of unresolved clouds and precipitation eas471_cumparam.odp JDW, EAS, U. Alberta Last modified: 29 Mar. 2016 (from Physical Parameterizations in Canadian Operational Models,
More informationOutline. Aim. Gas law. Pressure. Scale height Mixing Column density. Temperature Lapse rate Stability. Condensation Humidity.
Institute of Applied Physics University of Bern Outline A planetary atmosphere consists of different gases hold to the planet by gravity The laws of thermodynamics hold structure as vertical coordinate
More informationReferences: Parcel Theory. Vertical Force Balance. ESCI Cloud Physics and Precipitation Processes Lesson 3 - Stability and Buoyancy Dr.
References: ESCI 340 - Cloud Physics and Precipitation Processes Lesson 3 - Stability and Buoyancy Dr. DeCaria Glossary of Meteorology, 2nd ed., American Meteorological Society A Short Course in Cloud
More information2. Meridional atmospheric structure; heat and water transport. Recall that the most primitive equilibrium climate model can be written
2. Meridional atmospheric structure; heat and water transport The equator-to-pole temperature difference DT was stronger during the last glacial maximum, with polar temperatures down by at least twice
More informationBulk Boundary-Layer Model
Bulk Boundary-Layer Model David Randall Ball (1960) was the first to propose a model in which the interior of the planetary boundary layer (PBL) is well-mixed in the conservative variables, while the PBL
More informationAtmospheric modeling in the Climate System. Joe Tribbia NCAR.ESSL.CGD.AMP
Atmospheric modeling in the Climate System Joe Tribbia NCAR.ESSL.CGD.AMP The climate represents a coupled system consisting of an atmosphere, hydrosphere, biosphere, and cryosphere What is CCSM? Bio Geochemistry
More informationAtmospheric Circulation
Atmospheric Circulation (WAPE: General Circulation of the Atmosphere and Variability) François Lott, flott@lmd.ens.fr http://web.lmd.jussieu.fr/~flott 1) Mean climatologies and equations of motion a)thermal,
More informationA Gray-Radiation Aquaplanet Moist GCM. Part I: Static Stability and Eddy Scale
2548 J O U R N A L O F T H E A T M O S P H E R I C S C I E N C E S VOLUME 63 A Gray-Radiation Aquaplanet Moist GCM. Part I: Static Stability and Eddy Scale DARGAN M. W. FRIERSON Program in Applied and
More informationCONTENTS 1 MEASURES OF ATMOSPHERIC COMPOSITION
i CONTENTS 1 MEASURES OF ATMOSPHERIC COMPOSITION 1 1.1 MIXING RATIO 1 1.2 NUMBER DENSITY 2 1.3 PARTIAL PRESSURE 6 PROBLEMS 10 1.1 Fog formation 10 1.2 Phase partitioning of water in cloud 10 1.3 The ozone
More informationMay 3, :41 AOGS - AS 9in x 6in b951-v16-ch13 LAND SURFACE ENERGY BUDGET OVER THE TIBETAN PLATEAU BASED ON SATELLITE REMOTE SENSING DATA
Advances in Geosciences Vol. 16: Atmospheric Science (2008) Eds. Jai Ho Oh et al. c World Scientific Publishing Company LAND SURFACE ENERGY BUDGET OVER THE TIBETAN PLATEAU BASED ON SATELLITE REMOTE SENSING
More informationSurface layer parameterization in WRF
Surface layer parameteriation in WRF Laura Bianco ATOC 7500: Mesoscale Meteorological Modeling Spring 008 Surface Boundary Layer: The atmospheric surface layer is the lowest part of the atmospheric boundary
More informationA. Parodi 1, (1) CIMA Research Foundation, Italy. in cooperation with: K. A. Emanuel 2, and E. Foufoula-Georgiou 3 (2) EAPS, MIT, USA
Spatial and temporal evolution of deep moist convective processes: the role of microphysics A. Parodi 1, (1) CIMA Research Foundation, Italy in cooperation with: K. A. Emanuel 2, and E. Foufoula-Georgiou
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 information