Boundary Layer Parameterization

Size: px

Start display at page:

Download "Boundary Layer Parameterization"

Cody Townsend
5 years ago
Views:

1 Boundary Layer Parameterization Bob Plant With thanks to: R. Beare, S. Belcher, I. Boutle, O. Coceal, A. Grant, D. McNamara NWP Physics Lecture 3 Nanjing Summer School July 2014

2 Outline Motivation Some background on turbulence Explicit turbulence simulations Organization of the problem Surface layer parameterization K theory The boundary layer grey zone Boundary Layer parameterization p.1/48

3 Motivation Boundary Layer parameterization p.2/48

4 Effects of the boundary layer Control simulation, T+60 Simulation with no boundary layer turbulence, T mb 950mb Simulations with (left) and without (right) boundary layer, of storm of 30/10/00 Boundary Layer parameterization p.3/48

5 Role of friction: Ekman pumping tropopause ξ 2 free troposphere ξ 1 ξ1 > ξ 2 boundary layer top boundary layer Boundary layer convergence leads to ascent leads to spin-down of a barotropic vortex Barotropic vorticity equation, Dζ w =ζ, Dt z ζ = f +ξ surface Boundary Layer parameterization p.4/48

6 Effects on low-level stability Mid-level feature associated with dry intrusion Baroclinic frictional effects increase low-level stability over the low centre Reduces the strength of coupling between tropopause-level PV feature and surface temperature wave Boundary Layer parameterization p.5/48

7 Some background on turbulence Boundary Layer parameterization p.6/48

8 TKE KE = u = u + v + w + u + v + w +(uu + vv + ww ) inc. KE of mean flow, KE of turbulence and cross terms Reynolds average gives KE=mean flow KE + TKE 1 2 e= u + v 2 + w 2 2 Boundary Layer parameterization p.7/48

9 TKE Evolution De u g 1 w p w e ε = u w + wθ Dt z θ0 ρ z z Storage + Advection = Shear + Buoyancy + Pressure correlations + Dissipation A crucial distinction is between the convective BL (buoyancy generates TKE) and the stable BL (buoyancy destroys TKE) Boundary Layer parameterization p.8/48

10 Typical TKE Budgets Nocturnal, stable boundary layer boundary layer Morning, well-mixed Boundary Layer parameterization p.9/48

11 Split of eddy sizes Fourier transform of TKE energy into different k = 2π/λ A Energy-containing range. B Inertial subrange, k 5/3 C Viscous dissipation at Kolmogorov microscale η= 3 1/4 ν ε 1mm Boundary Layer parameterization p.10/48

12 Split of eddy sizes The turbulent energy cascade Boundary Layer parameterization p.11/48

13 Moisture transports Schematic based on boundary-layer moisture budget analysis: Divergence from high and convergence towards low within the boundary layer necessary to supply WCB with moisture Boundary Layer parameterization p.12/48

14 Explicit turbulence simulations Boundary Layer parameterization p.13/48

15 Direct Numerical Simulation Viscous dissipation important at 1mm, and BL height 500m Snapshot of DNS at Re= 500. yz slice through an array of cubes with large-scale flow out of page Fully resolving turbulence with DNS of Navier-Stokes needs 1018 grid points Inflate viscous scale and assume Reynolds number indpendence Boundary Layer parameterization p.14/48

16 Large Eddy Simulation Simulate only as far as the inertial subrange: capture the large eddies E= Z E(k)dk 0 Z kc =π/ x E(k)dk 0 Dissipation rate is ε = 2ν Z 0 k2 E(k)dk 6 2ν Z kc k2 E(k)dk 0 0 Without viscous scales there is no sink of TKE Boundary Layer parameterization p.15/48

17 Energy loss from LES Need a parameterization P of the energy drain ε 2 Z kc P(k)k2 E(k)dk 0 Ideally P acts close to kc only so well-resolved eddies are not affected Popular (and very simple) choice is a Smagorinsky scheme, which is effectively a diffusion with coefficient x Boundary Layer parameterization p.16/48

18 LES snapshots of w CBL simulated at x = y = 10m, z = 4m z = 80m (left) and 800m (right) Boundary Layer parameterization p.17/48

19 Organization of the problem Boundary Layer parameterization p.18/48

20 NWP and GCMs: Closure Problem On an NWP grid, no attempt to simulate turbulent eddies Parameterize full turbulent spectrum. u u u u 1 p u u v u w u +u +v +w f v+ = t x y z ρ x x y z Effects of turbulence are described by fluxes like u w Evolution equation for u w includes u w w etc Boundary Layer parameterization p.19/48

21 Mellor-Yamada hierachy Level 4 Carry (simplified) prognostic equations for all 2nd order moments with parameterization of 3rd order terms Level 3 Carry (simplified) prognostic equations for θ 2 and TKE with parameterization of 3rd order terms Level 2.5 Carry (simplified) prognostic equation for TKE with parameterization of 3rd order terms Level 2 Carry diagnostic equations for all 2nd order moments Level 1 Carry simplified diagnostic equations 2nd order moments, K theory Boundary Layer parameterization p.20/48

22 Level 3 e u g 1 = u w + wθ w p w e ε t z θ0 ρ z z θ 2 θ 2 = 2w θ w θ εθ t z z NWP/GCMs usually have diagnostic treatment but some use TKE approach Drop terms in red to get to the level 2.5, TKE approach Not well justified theoretically though: no fundamental reason to prefer kinetic to potential energy Boundary Layer parameterization p.21/48

23 Parameterization of terms After Kolmogorov, e e3/2 ε= τ l θ 2 e1/2 εθ l Terms involving pressure treated using return-to-isotropy ideas (Rotta) Many possibilities for 3rd order terms, from simple downgradient forms: eg, 2 u w u 2 z to much more sophisticated (ie, complicated!) methods Boundary Layer parameterization p.22/48

24 Surface layer parameterization Boundary Layer parameterization p.23/48

25 Surface Layer Similarity Similarity theory requires us to: 1. write down all of physically relevant quantities that we believe may control the strength and character of the turbulence 2. put these together into dimensionless combinations 3. any non-dimensional turbulent quantity must be a function of the dimensionless variables: we just need to measure that function Boundary Layer parameterization p.24/48

26 Monin-Obukhov theory: Step 1 Postulate that surface layer turbulence can be described by height z friction velocity (essentially the surface drag) u = q u w 0 a temperature scale that we can get from the surface heat flux w θ 0 T = u which is more conveniently expressed as a turbulent production of buoyancy (g/θ0 )T Boundary Layer parameterization p.25/48

27 Monin-Obukhov theory: Step 2 Construct dimensionless combinations from these three. Here z/l where L is the Obukhov length u2 θ0 L= kgt L measures relative strength of shear and buoyancy buoyancy becomes as important as shear at height z L +ve in stable conditions k is von Karman s constant, = 0.4 Boundary Layer parameterization p.26/48

28 Monin-Obukhov theory: Step 2 Express the variables we need in terms of dimensionless quantities u u = φm z kz ; θ T = φh z kz Now measure the dimensionless functions (obs or LES) Which must be universal if the scaling is correct Boundary Layer parameterization p.27/48

29 Monin-Obukhov Theory: Step 3 Functions φm (left) and φh (right) Boundary Layer parameterization p.28/48

30 Computing surface-layer fluxes Lowest model level typically at 10m Use Monin-Obukhov similarity theory u u = φm (z/l) z kz Integrate with height to get u for lowest level u z u= ln ψm k z0 ψm = Z z/l z0 /L (1 φm )(z/l) 1 d(z/l) Iterative calculations needed: u10 u L u10... Boundary Layer parameterization p.29/48

31 Connection to diagnostic equations We can rewrite MOST in the form w u Km = u = Km z 2 2 u (kz) φm z So that in the main u equation we have a diffusion structure u u = w u = Km t turbulence z z z Boundary Layer parameterization p.30/48

32 The outer boundary layer Boundary Layer parameterization p.31/48

33 Mixed Layer Similarity: Step 1 Surface heat flux w θ 0 drives convective eddies of the scale of h Scaling velocities for vertical velocity and temperature 1/3 g w θ 0h w = θ0 w θ 0 θ = w Dimensionless quantity z/h Boundary Layer parameterization p.32/48

34 Mixed Layer Similarity: Step 2 w 2 = w2 f1 (z/h) u 2 = w2 f2 (z/h) θ 2 = θ2 f3 (z/h) w θ = w θ f4 (z/h) w3 ε= f5 (z/h) h etc Boundary Layer parameterization p.33/48

35 Mixed layer similarity: Step 3 Left to right: w 2, u 2, θ 2, w θ, ε Boundary Layer parameterization p.34/48

36 K theory For the turbulent fluxes within the boundary layer, the K -theory formula is w u Km = u 2 λm fm (Ri) z u = Km z compare Km = 2 2 u (kz) φm z ie, generalizes MOST Typical eddy size kz or λm where l h = + λm kz l Boundary Layer parameterization p.35/48

37 Stability dependence z/l and Richardson number Ri play similar role in measuring relative importance of buoyancy and shear θg0 θ z Ri = u 2 z The flux Richardson number is sometimes used instead (g/θ0 )w θ Rf = w θ u z ratio of terms in TKE equation Boundary Layer parameterization p.36/48

38 Properties of K K theory should match to MOST for small z and to the relevant similarity theory in the outer layer, being based on the appropriate scaling parameters eg, for the well-mixed CBL a suitable choice is z 2 Km = w hg1 (z/h) = kw z 1 h eg, Holstlag and Bolville 1993; as used at ECMWF Boundary Layer parameterization p.37/48

39 A problem in the CBL w θ θ = Kh z = can diagnose Kh from LES/obs data Boundary Layer parameterization p.38/48

40 A solution for the CBL Simplest possibility is to introduce a non-local contribution w θ θ = Kh + Kh γ z where γ is simply a number Boundary Layer parameterization p.39/48

41 Parameterizations Based on Scalings Possible difficulties are: Developing good scalings for each boundary layer regime Good decison making needed for which regime to apply Handling transitions between regimes Boundary Layer parameterization p.40/48

42 Effects of radiation on buoyancy LW cooling at cloud top; SW heating within cloud produces source of buoyant motions Boundary Layer parameterization p.41/48

43 UM scheme Decision making about type is important Sc treated as BL cloud; shallow Cu separately Boundary Layer parameterization p.42/48

44 EDMF Eddy-diffusivity mass-flux treatment, w φ = Kdφ/dz + i Mi (φi φ) Mass flux component for largest, most energetic eddies of size h which produce the non-local transport Can be used as a treatment for Sc and (especially) shallow Cu Can be high sensitivity to bulk entrainment rate Boundary Layer parameterization p.43/48

45 The boundary layer grey zone Boundary Layer parameterization p.44/48

46 The grey zone NWP approaches to turbulence valid when all turbulence is parameterized LES approaches to turbulence valid Grey zone is difficult middle range when model grid is comparable to the size of energy-containing eddies Should we try to ignore or such eddies and use an NWP approach? Double counting? Or allow them and use an LES approach? Risks under counting? Boundary Layer parameterization p.45/48

47 A Perspective from LES Stochastic backscatter useful very near surface where l breaks down eg, improves profiles of dimensionless wind shear near surface Other LES models proposed for grey zone: dynamic model, tensorial model... Dry, neutral boundary layer, Weinbrecht 2006 Boundary Layer parameterization p.46/48

Perspective from NWP Small boundary layer fluctuations ( 0.1K) important for convective initiation Can easily shift the locations of precipitating

48 Perspective from NWP Small boundary layer fluctuations ( 0.1K) important for convective initiation Can easily shift the locations of precipitating cells e.g. Leoncini et al (2010) 1 Perturbation at 2000 UTC, 8 km Raining in both K Fraction Raining in control Raining in perturbed Time UTC Boundary Layer parameterization p.47/48

49 Conclusions Most NWP models have relatively simple turbulence modelling based on K -theory Simple methods are able to make use of similarity arguments These often work very well in the appropriate regimes, although transitions between regimes can be awkward and somewhat ad hoc...because the performance is not very much worse than using very much more complex and more expensive turbulence modelling approaches But as resolutions x h new approaches may be needed; ensemble-based modelling breaks down Boundary Layer parameterization p.48/48

VI Workshop Brasileiro de Micrometeorologia

VI Workshop Brasileiro de Micrometeorologia Validation of a statistic algorithm applied to LES model Eduardo Bárbaro, Amauri Oliveira, Jacyra Soares November 2009 Index Objective 1 Objective 2 3 Vertical Profiles Flow properties 4 Objective 1 The