Feedback and Sensitivity in Climate

Size: px
Start display at page:

Download "Feedback and Sensitivity in Climate"

Transcription

1 Feedback and Sensitivity in Climate Frédéric Hourdin June 19, 2009

2 2

3 Boundary layer processes in climate 2

4 Boundary layer processes in climate ps2pdf OK. Figure: ClimSI WV feedback model responses log(t y), f(x)=365*10 x -0.4 Oc Bk Tau 28j Tau prec 3m. Rep Eff "sigma.data" u (log10($1/365)): *(exp( e-02*f(x))-1) *(exp( e-03*f(x))-1) e-03*(exp( *f(x))-1) *(exp( e-02*f(x))-1) "/data/al1/pablo/cyclsta/sigma.datacysta" u (log10($1/365)):58 "/data/al1/pablo/ocblok/sigma.data" u (log10($1/365)): ClimSI : Comparaison entre Reponses efficaces

5 Boundary layer processes in climate 2

6 Boundary layer processes in climate 2

7 2

8 of the conservation equation Conservation equation v : wind field c : conserved quantity (dc/dt = 0), Advective form : Flux form : X : "average" or "large scale" variable X = X X : turbulent fluctuation vc = vc + v c c + vgradc = 0 t (1) ρc + div (ρvc) = 0 t (2) q t + V.grad q + 1 ρ div ( ρv c ) = 0 (3)

9 Under boundary layer approximations ( / x << / z) : c t + v.grad c = S c + 1 ρ z w c (4) 3D Dynamical core 20 km Physical parametrizations One grid mesh or atmospheric column.? 200 km c : θ, u, v, water (vapor and others), chemical compounds...

10 2

11 Diffusive or local formulations for the PBL w c = K z c z c t = ( ) c K z z z (5) Analogy with molecular viscosity. Down-gradient fluxes. Turbulence acts as a "mixing"

12 Turbulent diffusivity K z Prandlt (1925) mixing length : K z = l w or K z = l 2 v z Accounting for static stability (Ex. Louis 1979) K z = f (Ri)l 2 v z, Turbulent kinetic energy w 2 e = 1 2 with Ri = g θ θ z ( v z [u 2 + v 2 + w 2 ] ) 2 (6) e t = w u u z w v v z g θ w θ 1 w p w e ɛ (7) ρ z z Ex : Mellor and Yamada w φ φ = K φ z with K φ = l 2eS φ (Ri) Note : Imposing e t gives a coefficient of the form Eq. 6

13 Feedback Limitations of turbulent diffusion Diffusive approach : Random process 1 Small scale turbulence of size l << h with h = 1c c z Real turbulence Long range vertical transport (from the bottom to PBL top) Organized structures 15km Cloud streets on North of France (March 2009, MSG) Radar echoes in a dry convective boundary layer

14 Limitations of turbulent diffusion Idealized view of the dry convective boundary layer. Potential temperature initial final z h z i 0 Inversion layer Neutral (slightly stable) mixed layer Unstable surface layer θ Heat flux w θ In the mixed layer θ a Thermal plume wa Diffusive formulation w θ θ = K z = 0 or slightly < 0 α z (8) Uniform heating by the surface θ t w θ 0 z i (Cste > 0) (9) w θ z z i w z θ 0 > 0 (10) i

15 2

16 Extension of diffusive formulations Introduction of a countergradient term [ w θ = K z Γ θ ] = 0 with Γ 1K/km (11) z Imposed countergradient Deardorf, 1966 Revisited by Troen & Mart, 1986, Holtzlag & Boville, 1993, based on a similarity approach. Higher order closures Abdella & Mc Farlane, 1997, Randall et al., 1992, Lapen & Ransall, 2002 (12)

17 Boundary layer, clouds and convection 2

18 Boundary layer, clouds and convection 2

19 Boundary layer, clouds and convection 2

Outline. university-logo. Boundary layer parameterization Introduction. university-logo. Boundary layer parameterization Introduction.

Outline. university-logo. Boundary layer parameterization Introduction. university-logo. Boundary layer parameterization Introduction. and climate s in climate models Frédéric Hourdin June, 009 Boundary layer in te climate system Boundary layer in te "Eart System" Driven by te Global Cange studies, climate models are more and more complex

More information

Vertical resolution of numerical models. Atm S 547 Lecture 8, Slide 1

Vertical resolution of numerical models. Atm S 547 Lecture 8, Slide 1 Vertical resolution of numerical models Atm S 547 Lecture 8, Slide 1 M-O and Galperin stability factors Atm S 547 Lecture 8, Slide 2 Profile vs. forcing-driven turbulence parameterization Mellor-Yamada

More information

Model description of AGCM5 of GFD-Dennou-Club edition. SWAMP project, GFD-Dennou-Club

Model description of AGCM5 of GFD-Dennou-Club edition. SWAMP project, GFD-Dennou-Club Model description of AGCM5 of GFD-Dennou-Club edition SWAMP project, GFD-Dennou-Club Mar 01, 2006 AGCM5 of the GFD-DENNOU CLUB edition is a three-dimensional primitive system on a sphere (Swamp Project,

More information

A Simple Turbulence Closure Model

A Simple Turbulence Closure Model A Simple Turbulence Closure Model Atmospheric Sciences 6150 1 Cartesian Tensor Notation Reynolds decomposition of velocity: Mean velocity: Turbulent velocity: Gradient operator: Advection operator: V =

More information

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

1 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 information

Atmospheric Boundary Layers

Atmospheric Boundary Layers Lecture for International Summer School on the Atmospheric Boundary Layer, Les Houches, France, June 17, 2008 Atmospheric Boundary Layers Bert Holtslag Introducing the latest developments in theoretical

More information

Incorporation 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 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 information

Atm S 547 Boundary Layer Meteorology

Atm S 547 Boundary Layer Meteorology Lecture 9. Nonlocal BL parameterizations for clear unstable boundary layers In tis lecture Nonlocal K-profile parameterization (e. g. WRF-YSU) for dry convective BLs EDMF parameterizations (e. g. ECMWF)

More information

Prototype Instabilities

Prototype Instabilities Prototype Instabilities David Randall Introduction Broadly speaking, a growing atmospheric disturbance can draw its kinetic energy from two possible sources: the kinetic and available potential energies

More information

The atmospheric boundary layer: Where the atmosphere meets the surface. The atmospheric boundary layer:

The 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 information

A Discussion on The Effect of Mesh Resolution on Convective Boundary Layer Statistics and Structures Generated by Large-Eddy Simulation by Sullivan

A 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 information

Convective Fluxes: Sensible and Latent Heat Convective Fluxes Convective fluxes require Vertical gradient of temperature / water AND Turbulence ( mixing ) Vertical gradient, but no turbulence: only very

More information

τ xz = τ measured close to the the surface (often at z=5m) these three scales represent inner unit or near wall normalization

τ xz = τ measured close to the the surface (often at z=5m) these three scales represent inner unit or near wall normalization τ xz = τ measured close to the the surface (often at z=5m) these three scales represent inner unit or near wall normalization Note that w *3 /z i is used to normalized the TKE equation in case of free

More information

Unified Cloud and Mixing Parameterizations of the Marine Boundary Layer: EDMF and PDF-based cloud approaches

Unified 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 information

A Simple Turbulence Closure Model. Atmospheric Sciences 6150

A Simple Turbulence Closure Model. Atmospheric Sciences 6150 A Simple Turbulence Closure Model Atmospheric Sciences 6150 1 Cartesian Tensor Notation Reynolds decomposition of velocity: V = V + v V = U i + u i Mean velocity: V = Ui + V j + W k =(U, V, W ) U i =(U

More information

Atm S 547 Boundary Layer Meteorology

Atm 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 information

Mellor-Yamada Level 2.5 Turbulence Closure in RAMS. Nick Parazoo AT 730 April 26, 2006

Mellor-Yamada Level 2.5 Turbulence Closure in RAMS. Nick Parazoo AT 730 April 26, 2006 Mellor-Yamada Level 2.5 Turbulence Closure in RAMS Nick Parazoo AT 730 April 26, 2006 Overview Derive level 2.5 model from basic equations Review modifications of model for RAMS Assess sensitivity of vertical

More information

6A.3 Stably stratified boundary layer simulations with a non-local closure model

6A.3 Stably stratified boundary layer simulations with a non-local closure model 6A.3 Stably stratified boundary layer simulations with a non-local closure model N. M. Colonna, E. Ferrero*, Dipartimento di Scienze e Tecnologie Avanzate, University of Piemonte Orientale, Alessandria,

More information

Interaction of boundary-layer turbulence and convection in the AROME model: Preliminary results

Interaction of boundary-layer turbulence and convection in the AROME model: Preliminary results Interaction of boundary-layer turbulence and convection in the AROME model: Preliminary results Balázs Szintai Hungarian Meteorological Service COST ES0905 Workshop Cambridge, England 23 March 2010 Outline

More information

Introduction to Heat and Mass Transfer. Week 12

Introduction to Heat and Mass Transfer. Week 12 Introduction to Heat and Mass Transfer Week 12 Next Topic Convective Heat Transfer» Heat and Mass Transfer Analogy» Evaporative Cooling» Types of Flows Heat and Mass Transfer Analogy Equations governing

More information

Bulk Boundary-Layer Model

Bulk 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 information

3-Fold Decomposition EFB Closure for Convective Turbulence and Organized Structures

3-Fold Decomposition EFB Closure for Convective Turbulence and Organized Structures 3-Fold Decomposition EFB Closure for Convective Turbulence and Organized Structures Igor ROGACHEVSKII and Nathan KLEEORIN Ben-Gurion University of the Negev, Beer-Sheva, Israel N.I. Lobachevsky State University

More information

Surface layer parameterization in WRF

Surface 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 information

Parameterization of the Dry Convective Boundary Layer Based on a Mass Flux Representation of Thermals

Parameterization of the Dry Convective Boundary Layer Based on a Mass Flux Representation of Thermals 15 MARCH 00 HOURDIN ET AL. 1105 Parameterization of the Dry Convective Boundary Layer Based on a Mass Flux Representation of Thermals FRÉDÉRIC HOURDIN, FLEUR COUVREUX, AND LAURENT MENUT Laboratoire de

More information

A Note on the Estimation of Eddy Diffusivity and Dissipation Length in Low Winds over a Tropical Urban Terrain

A Note on the Estimation of Eddy Diffusivity and Dissipation Length in Low Winds over a Tropical Urban Terrain Pure appl. geophys. 160 (2003) 395 404 0033 4553/03/020395 10 Ó Birkhäuser Verlag, Basel, 2003 Pure and Applied Geophysics A Note on the Estimation of Eddy Diffusivity and Dissipation Length in Low Winds

More information

PBL and precipitation interaction and grey-zone issues

PBL and precipitation interaction and grey-zone issues PBL and precipitation interaction and grey-zone issues Song-You Hong, Hyun-Joo Choi, Ji-Young Han, and Young-Cheol Kwon (Korea Institute of Atmospheric Prediction Systems: KIAPS) Overview of KIAPS (seminar

More information

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

Lecture 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 information

Logarithmic velocity profile in the atmospheric (rough wall) boundary layer

Logarithmic velocity profile in the atmospheric (rough wall) boundary layer Logarithmic velocity profile in the atmospheric (rough wall) boundary layer P =< u w > U z = u 2 U z ~ ε = u 3 /kz Mean velocity profile in the Atmospheric Boundary layer Experimentally it was found that

More information

Wind Flow Modeling The Basis for Resource Assessment and Wind Power Forecasting

Wind Flow Modeling The Basis for Resource Assessment and Wind Power Forecasting Wind Flow Modeling The Basis for Resource Assessment and Wind Power Forecasting Detlev Heinemann ForWind Center for Wind Energy Research Energy Meteorology Unit, Oldenburg University Contents Model Physics

More information

The parametrization of the planetary boundary layer May 1992

The parametrization of the planetary boundary layer May 1992 The parametrization of the planetary boundary layer May 99 By Anton Beljaars European Centre for Medium-Range Weather Forecasts Table of contents. Introduction. The planetary boundary layer. Importance

More information

2σ 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

2σ 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 information

Lecture 7. Science A-30 February 21, 2008 Air may be forced to move up or down in the atmosphere by mechanical forces (wind blowing over an obstacle,

Lecture 7. Science A-30 February 21, 2008 Air may be forced to move up or down in the atmosphere by mechanical forces (wind blowing over an obstacle, Lecture 7. Science A-30 February 21, 2008 Air may be forced to move up or down in the atmosphere by mechanical forces (wind blowing over an obstacle, like a mountain) or by buoyancy forces. Air that is

More information

A Combined Local and Nonlocal Closure Model for the Atmospheric Boundary Layer. Part I: Model Description and Testing

A Combined Local and Nonlocal Closure Model for the Atmospheric Boundary Layer. Part I: Model Description and Testing SEPTEMBER 2007 P L E I M 1383 A Combined Local and Nonlocal Closure Model for the Atmospheric Boundary Layer. Part I: Model Description and Testing JONATHAN E. PLEIM Atmospheric Sciences Modeling Division,*

More information

Testing and Improving Pacific NW PBL forecasts

Testing and Improving Pacific NW PBL forecasts Testing and Improving Pacific NW PBL forecasts Chris Bretherton and Matt Wyant University of Washington Eric Grimit 3Tier NASA MODIS Image Testing and Improving Pacific NW PBL forecasts PBL-related forecast

More information

WQMAP (Water Quality Mapping and Analysis Program) is a proprietary. modeling system developed by Applied Science Associates, Inc.

WQMAP (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 information

A B C D PROBLEMS Dilution of power plant plumes. z z z z

A 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 information

PALM - Cloud Physics. Contents. PALM group. last update: Monday 21 st September, 2015

PALM - 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 information

Assessment of Physical Parameterizations Using a Global Climate Model with Stretchable Grid and Nudging

Assessment of Physical Parameterizations Using a Global Climate Model with Stretchable Grid and Nudging 1474 M O N T H L Y W E A T H E R R E V I E W VOLUME 135 Assessment of Physical Parameterizations Using a Global Climate Model with Stretchable Grid and Nudging O. COINDREAU Commisariat à l Energie Atomique,

More information

IMPROVEMENT OF THE K-PROFILE MODEL FOR THE PLANETARY BOUNDARY LAYER BASED ON LARGE EDDY SIMULATION DATA

IMPROVEMENT OF THE K-PROFILE MODEL FOR THE PLANETARY BOUNDARY LAYER BASED ON LARGE EDDY SIMULATION DATA IMPROVEMENT OF THE K-PROFILE MODEL FOR THE PLANETARY BOUNDARY LAYER BASED ON LARGE EDDY SIMULATION DATA Y. NOH,W.G.CHEONandS.Y.HONG Department of Atmospheric Sciences, Yonsei University, 134 Shinchon-dong,

More information

Use of Satellite Observations to Measure Air-Sea Coupling and to Validate Its Estimates from Numerical Atmospheric Models

Use of Satellite Observations to Measure Air-Sea Coupling and to Validate Its Estimates from Numerical Atmospheric Models Use of Satellite Observations to Measure Air-Sea Coupling and to Validate Its Estimates from Numerical Atmospheric Models Natalie Perlin, Dudley Chelton, Simon de Szoeke College of Earth, Ocean, and Atmospheric

More information

Chapter (3) TURBULENCE KINETIC ENERGY

Chapter (3) TURBULENCE KINETIC ENERGY Chapter (3) TURBULENCE KINETIC ENERGY 3.1 The TKE budget Derivation : The definition of TKE presented is TKE/m= e = 0.5 ( u 2 + v 2 + w 2 ). we recognize immediately that TKE/m is nothing more than the

More information

Development and validation of an operational model for air quality predictions SIRANE Salizzoni P., Garbero V., Soulhac L.,

Development and validation of an operational model for air quality predictions SIRANE Salizzoni P., Garbero V., Soulhac L., Ecole Centrale de Lyon Laboratoire de Mécanique des Fluides et d'acoustique Development and validation of an operational model for air quality predictions SIRANE Salizzoni P., Garbero V., Soulhac L., Flow

More information

Anomalous solar heating dependence of Venus s cloud-level convection

Anomalous solar heating dependence of Venus s cloud-level convection Anomalous solar heating dependence of Venus s cloud-level convection T. Higuchi (Univ. Tokyo), T. Imamura (JAXA), Y. Maejima (MRI, JMA), M. Takagi (Kyoto Sangyo Univ.), N. Sugimoto (Keio Univ.), K. Ikeda

More information

Turbulence (January 7, 2005)

Turbulence (January 7, 2005) http://www.tfd.chalmers.se/gr-kurs/mtf071 70 Turbulence (January 7, 2005) The literature for this lecture (denoted by LD) and the following on turbulence models is: L. Davidson. An Introduction to Turbulence

More information

Météo-France Physics (AROME & ARPEGE) Eric Bazile, Yves Bouteloup, Rachel Honnert,

Météo-France Physics (AROME & ARPEGE) Eric Bazile, Yves Bouteloup, Rachel Honnert, Recent developments in Météo-France Physics (AROME & ARPEGE Eric Bazile, Yves Bouteloup, Rachel Honnert, Sébastien Riette, Yann Seity 34th Ewglam & 19th SRNWP joined meetings Helsinki 8-11 october 2012

More information

Taking into account the gustiness due to free and deep convection for the representation of air-sea fluxes. - in the LMDZ model -

Taking into account the gustiness due to free and deep convection for the representation of air-sea fluxes. - in the LMDZ model - Taking into account the gustiness due to free and deep convection for the representation of air-sea fluxes - in the LMDZ model - Alina Găinuşă-Bogdan, Frédéric Hourdin, Abdoul Traore, Catherine Rio, Jean-Yves

More information

Importance of minor treatments in parameterizations in GCMs for the cloud representations and the cloud feedbacks

Importance of minor treatments in parameterizations in GCMs for the cloud representations and the cloud feedbacks The WCRP Grand Challenge Workshop, 25 March. 2014 Importance of minor treatments in parameterizations in GCMs for the cloud representations and the cloud feedbacks Hideaki Kawai Meteorological Research

More information

Radiation, Sensible Heat Flux and Evapotranspiration

Radiation, 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 information

Lecture 12. The diurnal cycle and the nocturnal BL

Lecture 12. The diurnal cycle and the nocturnal BL Lecture 12. The diurnal cycle and the nocturnal BL Over flat land, under clear skies and with weak thermal advection, the atmospheric boundary layer undergoes a pronounced diurnal cycle. A schematic and

More information

Arctic Boundary Layer

Arctic Boundary Layer Annual Seminar 2015 Physical processes in present and future large-scale models Arctic Boundary Layer Gunilla Svensson Department of Meteorology and Bolin Centre for Climate Research Stockholm University,

More information

A study of the moist Entropy in Meteorology.

A study of the moist Entropy in Meteorology. Workshop on the Atmospheric Modelling 8-10 February, 2011.Toulouse, France. A study of the moist Entropy in Meteorology. Pascal MARQUET (Météo-France. DPrévi / LABO) Contents 1) Motivations : to study

More information

OCN/ATM/ESS 587. Ocean circulation, dynamics and thermodynamics.

OCN/ATM/ESS 587. Ocean circulation, dynamics and thermodynamics. OCN/ATM/ESS 587 Ocean circulation, dynamics and thermodynamics. Equation of state for seawater General T/S properties of the upper ocean Heat balance of the upper ocean Upper ocean circulation Deep circulation

More information

NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS

NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS June - July, 5 Melbourne, Australia 9 7B- NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS Werner M.J. Lazeroms () Linné FLOW Centre, Department of Mechanics SE-44

More information

Implementation of the Quasi-Normal Scale Elimination (QNSE) Model of Stably Stratified Turbulence in WRF

Implementation of the Quasi-Normal Scale Elimination (QNSE) Model of Stably Stratified Turbulence in WRF Implementation of the Quasi-ormal Scale Elimination (QSE) odel of Stably Stratified Turbulence in WRF Semion Sukoriansky (Ben-Gurion University of the egev Beer-Sheva, Israel) Implementation of the Quasi-ormal

More information

Diffusive Transport Enhanced by Thermal Velocity Fluctuations

Diffusive Transport Enhanced by Thermal Velocity Fluctuations Diffusive Transport Enhanced by Thermal Velocity Fluctuations Aleksandar Donev 1 Courant Institute, New York University & Alejandro L. Garcia, San Jose State University John B. Bell, Lawrence Berkeley

More information

The Atmospheric Boundary Layer. The Surface Energy Balance (9.2)

The Atmospheric Boundary Layer. The Surface Energy Balance (9.2) The Atmospheric Boundary Layer Turbulence (9.1) The Surface Energy Balance (9.2) Vertical Structure (9.3) Evolution (9.4) Special Effects (9.5) The Boundary Layer in Context (9.6) Fair Weather over Land

More information

TURBULENT KINETIC ENERGY

TURBULENT KINETIC ENERGY TURBULENT KINETIC ENERGY THE CLOSURE PROBLEM Prognostic Moment Equation Number Number of Ea. fg[i Q! Ilial.!.IokoQlI!!ol Ui au. First = at au.'u.' '_J_ ax j 3 6 ui'u/ au.'u.' a u.'u.'u k ' Second ' J =

More information

Part I: Dry Convection

Part I: Dry Convection Turbulent dispersion and chemical transformation in the atmospheric boundary layer: Part I: Dry Convection DISPERSION Thanks: Alessandro Dosio Jordi Vilà-Guerau de Arellano WA G E N I N G E N U N I VE

More information

Planetary boundary layer schemes in RegCM: evaluation and impact on the climate change signal over Europe and Mediterranean region

Planetary boundary layer schemes in RegCM: evaluation and impact on the climate change signal over Europe and Mediterranean region Seventh ICTP Workshop on the Theory and Use of Regional Climate Models Planetary boundary layer schemes in RegCM: evaluation and impact on the climate change signal over Europe and Mediterranean region

More information

TABLE OF CONTENTS CHAPTER TITLE PAGE

TABLE OF CONTENTS CHAPTER TITLE PAGE v TABLE OF CONTENTS CHAPTER TITLE PAGE TABLE OF CONTENTS LIST OF TABLES LIST OF FIGURES LIST OF SYMBOLS LIST OF APPENDICES v viii ix xii xiv CHAPTER 1 INTRODUCTION 1.1 Introduction 1 1.2 Literature Review

More information

The Ocean-Atmosphere System II: Oceanic Heat Budget

The 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 information

Atmospheric 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 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 information

M.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. 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 information

Presentation A simple model of multiple climate regimes

Presentation A simple model of multiple climate regimes A simple model of multiple climate regimes Kerry Emanuel March 21, 2012 Overview 1. Introduction 2. Essential Climate Feedback Processes Ocean s Thermohaline Circulation, Large-Scale Circulation of the

More information

Frictional boundary layers

Frictional boundary layers Chapter 12 Frictional boundary layers Turbulent convective boundary layers were an example of boundary layers generated by buoyancy forcing at a boundary. Another example of anisotropic, inhomogeneous

More information

A parameterization of third-order moments for the dry convective boundary layer

A parameterization of third-order moments for the dry convective boundary layer Boundary-Layer Meteorol (2006) 120:437 454 DOI 10.1007/s10546-006-9071-7 ORIGINAL ARTICLE A parameterization of third-order moments for the dry convective boundary layer Séverine Tomas Valéry Masson Received:

More information

CFD calculations of the test 2-4 experiments. Author: G. de With

CFD calculations of the test 2-4 experiments. Author: G. de With CFD calculations of the test 2-4 experiments Author: G. de With 34. Model setup and boundary conditions Dimensions CFD model: x=1000m / y=100m / z=2000m. CFD Model: Transient simulation, with steady-state

More information

Improving non-local parameterization of the convective boundary layer

Improving non-local parameterization of the convective boundary layer Boundary-Layer Meteorol DOI ORIGINAL PAPER Improving non-local parameterization of the convective boundary layer Zbigniew Sorbjan 1, 2, 3 Submitted on 21 April 2008 and revised on October 17, 2008 Abstract

More information

Sungsu Park, Chris Bretherton, and Phil Rasch

Sungsu Park, Chris Bretherton, and Phil Rasch Improvements in CAM5 : Moist Turbulence, Shallow Convection, and Cloud Macrophysics AMWG Meeting Feb. 10. 2010 Sungsu Park, Chris Bretherton, and Phil Rasch CGD.NCAR University of Washington, Seattle,

More information

The applicability of Monin Obukhov scaling for sloped cooled flows in the context of Boundary Layer parameterization

The applicability of Monin Obukhov scaling for sloped cooled flows in the context of Boundary Layer parameterization Julia Palamarchuk Odessa State Environmental University, Ukraine The applicability of Monin Obukhov scaling for sloped cooled flows in the context of Boundary Layer parameterization The low-level katabatic

More information

Introducing LMDZ physical schemes in WRF

Introducing LMDZ physical schemes in WRF Introducing LMDZ physical schemes in WRF L. Fita 1, F. Hourdin, L. Farihead 1 Laboratorie de Météorologie Dynamique, IPSL, Paris-6, UPMC, Tr 44-45 2nd fl., Jussieu, 75005, Paris, France LMD: Reunion Climat

More information

Modelling Convection Over Arctic Leads With LES and a Non-Eddy-Resolving Microscale Model

Modelling Convection Over Arctic Leads With LES and a Non-Eddy-Resolving Microscale Model JOURNAL OF GEOPHYSICAL RESEARCH, VOL.???, XXXX, DOI:1.129/, 1 2 Modelling Convection Over Arctic Leads With LES and a Non-Eddy-Resolving Microscale Model C. Lüpkes, 1 V.M. Gryanik, 1,2 B. Witha, 3 M. Gryschka,

More information

Multi-Scale Modeling of Turbulence and Microphysics in Clouds. Steven K. Krueger University of Utah

Multi-Scale Modeling of Turbulence and Microphysics in Clouds. Steven K. Krueger University of Utah Multi-Scale Modeling of Turbulence and Microphysics in Clouds Steven K. Krueger University of Utah 10,000 km Scales of Atmospheric Motion 1000 km 100 km 10 km 1 km 100 m 10 m 1 m 100 mm 10 mm 1 mm Planetary

More information

LECTURE 28. The Planetary Boundary Layer

LECTURE 28. The Planetary Boundary Layer LECTURE 28 The Planetary Boundary Layer The planetary boundary layer (PBL) [also known as atmospheric boundary layer (ABL)] is the lower part of the atmosphere in which the flow is strongly influenced

More information

Development of a High-Resolution Coupled Atmosphere-Ocean-Land General Circulation Model for Climate System Studies

Development of a High-Resolution Coupled Atmosphere-Ocean-Land General Circulation Model for Climate System Studies Chapter 1 Earth Science Development of a High-Resolution Coupled Atmosphere-Ocean-Land General Circulation Model for Climate System Studies Project Representative Tatsushi Tokioka Frontier Research Center

More information

A new theory for moist convection in statistical equilibrium

A new theory for moist convection in statistical equilibrium A new theory for moist convection in statistical equilibrium A. Parodi(1), K. Emanuel(2) (2) CIMA Research Foundation,Savona, Italy (3) EAPS, MIT, Boston, USA True dynamics: turbulent, moist, non-boussinesq,

More information

3D experiments with a stochastic convective parameterisation scheme

3D experiments with a stochastic convective parameterisation scheme 3D experiments with a stochastic convective parameterisation scheme R. J. Keane and R. S. Plant 3D experiments with a stochastic convective parameterisation scheme p.1/17 Outline Introduction to the Plant-Craig

More information

Parcel Model. Atmospheric Sciences September 30, 2012

Parcel Model. Atmospheric Sciences September 30, 2012 Parcel Model Atmospheric Sciences 6150 September 30, 2012 1 Governing Equations for Precipitating Convection For precipitating convection, we have the following set of equations for potential temperature,

More information

Hydrogen Vehicle Leak Modelling in Indoor Ventilated Environments.

Hydrogen Vehicle Leak Modelling in Indoor Ventilated Environments. Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Hydrogen Vehicle Leak Modelling in Indoor Ventilated Environments. A. Hallgarth *1, A. Zayer 1, A. Gatward 2 and J. Davies 2 1 Hazard Research

More information

Needs work : define boundary conditions and fluxes before, change slides Useful definitions and conservation equations

Needs work : define boundary conditions and fluxes before, change slides Useful definitions and conservation equations Needs work : define boundary conditions and fluxes before, change slides 1-2-3 Useful definitions and conservation equations Turbulent Kinetic energy The fluxes are crucial to define our boundary conditions,

More information

From small-scale turbulence to large-scale convection: a unified scale-adaptive EDMF parameterization

From small-scale turbulence to large-scale convection: a unified scale-adaptive EDMF parameterization From small-scale turbulence to large-scale convection: a unified scale-adaptive EDMF parameterization Kay Sušelj 1, Joao Teixeira 1 and Marcin Kurowski 1,2 1 JET PROPULSION LABORATORY/CALIFORNIA INSTITUTE

More information

Thermodynamics of Atmospheres and Oceans

Thermodynamics 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 information

3D Mesocale Modeling of the Venus Atmosphere

3D Mesocale Modeling of the Venus Atmosphere 3D Mesocale Modeling of the Venus Atmosphere Maxence Lefèvre, Sébastien Lebonnois and Aymeric Spiga maxence.lefevre@lmd.jussieu.fr Laboratoire de Météorologie Dynamique, Paris, FRANCE CPS, 29th March 2018

More information

Modelling ice shelf basal melt with Glimmer-CISM coupled to a meltwater plume model

Modelling ice shelf basal melt with Glimmer-CISM coupled to a meltwater plume model Modelling ice shelf basal melt with Glimmer-CISM coupled to a meltwater plume model Carl Gladish NYU CIMS February 17, 2010 Carl Gladish (NYU CIMS) Glimmer-CISM + Plume February 17, 2010 1 / 24 Acknowledgements

More information

Pollutant dispersion in urban geometries

Pollutant dispersion in urban geometries Pollutant dispersion in urban geometries V. Garbero 1, P. Salizzoni 2, L. Soulhac 2 1 Politecnico di Torino - Department of Mathematics 2 Ecole Centrale de Lyon - Laboratoire de Méchaniques des Fluides

More information

Today s Lecture: Atmosphere finish primitive equations, mostly thermodynamics

Today 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 information

Radiative-Convective Instability. Kerry Emanuel Massachusetts Institute of Technology

Radiative-Convective Instability. Kerry Emanuel Massachusetts Institute of Technology Radiative-Convective Instability Kerry Emanuel Massachusetts Institute of Technology Program Basic radiative-convective equilibrium Macro-instability of the RC state Some consequences Radiative Equilibrium

More information

The Simulation of Wraparound Fins Aerodynamic Characteristics

The Simulation of Wraparound Fins Aerodynamic Characteristics The Simulation of Wraparound Fins Aerodynamic Characteristics Institute of Launch Dynamics Nanjing University of Science and Technology Nanjing Xiaolingwei 00 P. R. China laithabbass@yahoo.com Abstract:

More information

HYCOM Overview. By George Halliwell, 13 February 2002

HYCOM Overview. By George Halliwell, 13 February 2002 HYCOM Overview By George Halliwell, 13 February 2002 General Remarks The HYbrid Coordinate Ocean Model (HYCOM; (Halliwell et al., 1998; 2000; Bleck, 2001) is a primitive equation ocean general circulation

More information

and Meteorology Amir A. Aliabadi July 16, 2017 The 2017 Connaught Summer Institute in Arctic Science: Atmosphere, Cryosphere, and Climate

and Meteorology Amir A. Aliabadi July 16, 2017 The 2017 Connaught Summer Institute in Arctic Science: Atmosphere, Cryosphere, and Climate The 2017 Connaught Summer Institute in Science: Atmosphere, Cryosphere, and Climate July 16, 2017 Modelling Outline Modelling Modelling A zonal-mean flow field is said to be hydro-dynamically unstable

More information

Meteorology 6150 Cloud System Modeling

Meteorology 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 information

A TWO-SCALE MIXING FORMULATION FOR THE ATMOSPHERIC BOUNDARY LAYER

A TWO-SCALE MIXING FORMULATION FOR THE ATMOSPHERIC BOUNDARY LAYER A TWO-SCALE MIXING FORMULATION FOR THE ATMOSPHERIC BOUNDARY LAYER M. FRECH* and L. MAHRT Oceanic and Atmospheric Sciences, Oregon State University, Corvallis, OR 97331, USA (Received in final form 10 August,

More information

Using Cloud-Resolving Models for Parameterization Development

Using 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 information

Numerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling

Numerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline

More information

Turbulent transport of the energy in the entrainment interface layer

Turbulent transport of the energy in the entrainment interface layer Turbulent transport of the energy in the entrainment interface layer Marta K. Kopec 1, S. P. Malinowski1, Z. P. Piotrowski2 1 nstitute 2 nstitute of Geophysics, Faculty of Physics, University of Warsaw,

More information

A Reexamination of the Emergy Input to a System from the Wind

A Reexamination of the Emergy Input to a System from the Wind Emergy Synthesis 9, Proceedings of the 9 th Biennial Emergy Conference (2017) 7 A Reexamination of the Emergy Input to a System from the Wind Daniel E. Campbell, Laura E. Erban ABSTRACT The wind energy

More information

Analysis of the role of the planetary boundary layer schemes during a severe convective storm

Analysis of the role of the planetary boundary layer schemes during a severe convective storm Annales Geophysicae (2004) 22: 1861 1874 SRef-ID: 1432-0576/ag/2004-22-1861 European Geosciences Union 2004 Annales Geophysicae Analysis of the role of the planetary boundary layer schemes during a severe

More information

Why do GCMs have trouble with the MJO?

Why do GCMs have trouble with the MJO? Why do GCMs have trouble with the MJO? The Madden-Julian Oscillation West East 200 [hpa] 500 Cool & dry Cool & dry p 700 850 SST Lag Day +20 +15 +10 +5 0-5 -10-15 -20 ~20 days ~10 days ~10-15 days

More information

The 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 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 information

LECTURE NOTES ON THE Planetary Boundary Layer

LECTURE NOTES ON THE Planetary Boundary Layer LECTURE NOTES ON THE Planetary Boundary Layer Chin-Hoh Moeng prepared for lectures given at the Department of Atmospheric Science, CSU in 1994 & 1998 and at the Department of Atmospheric and Oceanic Sciences,

More information