Numerical Weather Prediction Parametrization of diabatic processes Convection I

1. Numerical Weather Prediction Parametrization of diabatic processesConvection I Peter Bechtold and Christian Jakob

2. Convection • Lectures: • The nature of convection • Parametrisation of convection • The ECMWF mass-flux parametrisation and Tracer transport • Forecasting of Convection • Wavelets, Neural Networks and EOFs (in preparation) • Exercises • The big secret !!!

3. Convection • Aim of Lectures: The aim of the lecture is only to give a rough overview of convective phenomena and parameterisation concepts in numerical models. The student is not expected to be able to directly write a new convection code- the development and full validation of a new convection scheme takes time (years). There are many details in a parameterisation, and the best exercise is to start with an existing code, run some offline examples on Soundings and dig in line by line ….. there is already a trend toward explicit representation of convection in limited area NWP (no need for parameterization) but for global we are not there yet, and still will need parameterizations for the next decade • Offline convection Code: Can be obtained from peter.bechtold@ecmwf.int

4. Convection Parametrisation and Dynamics - Text Books • Emanuel, 1994: Atmospheric convection, OUP • Houze R., 1993: Coud dynamics, AP • Holton, 2004: An introduction to Dynamic Meteorology, AP • Bluestein, 1993: Synoptic-Dynamic meteorology in midlatitudes, Vol II. OUP • Peixoto and Ort, 1992: The physics of climate. American Institute of Physics • Emanuel and Raymond, 1993: The representation of cumulus convection in numerical models. AMS Meteor. Monogr. • Smith, 1997: The physics and parametrization of moist atmospheric convection. Kluwer • Dufour et v. Mieghem: Thermodynamique de l’Atmosphère, 1975: Institut Royal météorologique de Belgique • Bohren and Albrecht, 1998: Atmospheric Thermodynamics.OUP AP=Academic Press; OUP=Oxford University Press

5. How does it look like ?

6. Intense deep Deep and shallow convection ITCZ at 10ºN Sc convection SPCZ African Squall lines Moist convection : Global IR GOES METEOSAT 7/04/2003

7. Convection and role of water vapor Interaction of Tropics and midlatitudes: Dry air intrusions modulate convection (Rossby wave breaking)

8. Convection and upper-level Divergence(determine divergence from variation of cold cloud top areas) Mc is the convective mass flux (see later)

9. Outline General: • Convection and tropical circulations • Midlatitude Convection • Shallow Convection Useful concepts and tools: • Buoyancy • Convective Available Potential Energy • Soundings and thermodynamic diagrams • Convective quasi-equilibrium • Large-scale observational budgets

10. Convection and tropical circulations (1)It’s raining again… 2000/2001 rainfall rate as simulated by IFS CY30R2 and compared to GPCP obs about 3 mm/day is falling globally, but most i.e. 5-7 mm/day in the Tropics

11. Model Tendencies – Tropical Equilibria Nevertheless, the driving force for atmospheric dynamics and convection is the radiation Above the boundary layer, there is an equilibrium Radiation-Clouds-Dynamics-Convection for Temperature, whereas for moisture there is roughly an equilibrium between dynamical transport (moistening) and convective drying. - Global Budgets are very similar

12. Convection and tropical circulations (2)proxy distribution of deep and shallow convective clouds as obtained from IFS

13. Shallow Convection (1)Field of tropical oceanic Cumuli

14. Shallow (Boundary-Layer) Convection (5) Basic physics of the trade-cumulus boundary layer Emanuel, 1994

15. A third convective mode Recent studies indicate, that there is a third important mode of convection (besides deep and shallow) in the tropics consisting of mainly cumulus congestus clouds terminating near the melting level at around 5 km. Johnson et al., 1999, JCL

16. Convection and tropical circulations (3)ITCZ and the Hadley meridional circulation: the role of trade-wind cumuli and deep tropical towers

17. Convection and tropical circulations (4)The Walker zonal Circulation From Salby (1996)

18. Convection and tropical circulations (5)Tropical waves: Rossby, Kelvin, Gravity, African easterly waves a Squall line

19. Convection and tropical circulations (7) The KELVIN wave 50/50 rotational/divergent 50/50 KE/PE Strong zonal wind along the Equator Symmetric around the Equator Eastward moving ~18 m/s

20. Convection and tropical circulations (8)The Kelvin wave, OLR composite

21. Convection and tropical circulations (9) The (n=1) Equatorial Rossby wave Symmetric KE>PE KE max at Equator,PE max off the Equator Westward moving ~ 5 m/s

22. Convection and tropical circulations (10)The Equatorial Rossby wave, OLR composite

23. African easterly waves • African easterly waves have periods of 2-6 days, typical wavelengths of about 2500 km and propagation speeds around -8 m/s. • These waves are thought to originate from barotropic and baroclinic instability, but the effects of diabatic (Cumulus) heating, the diurnal cycle and orography also modulate the waves. For instability to exist, the quasi-potential vorticity gradient must change sign in the domain, i.e for tropical North Africa it must become negative. • Although the shear instability associated with the jet is present throughout the rainy season, the waves appear to contribute to the development of rainfall systems only during late summer as only then the can access the necessary moisture in the low-level monsoon flow. • The exit region of the Tropcial easterly Jet – which is the consequence of the outflow=divergence due to the Asian Monsoon - might also have an effect on convection over tropical North Africa • The waves are generally confined to a latitudinal zone close and south of the Jet Further reading: Diedhiou et al. (1998, GRL), Nicholson and Grist (2003,J. Clim), Hsieh and Cook (2004, MWR), Grist (2002, MWR)

24. West-African meteorology – easterly waves Mid-level dry “Harmattan” Low-levelMonsoonflow Upper-level easterlies Monsoon flow ,Easterly waves, and midlatitude-tropical mixing

25. Hovmoeller plots as an easy way to plot wave (propagation) Analysis 10.8-9.9 2005 Comparison Analysis – Forecast for African easterly waves. No filtering required from a series of 2-day Forecasts

26. Convection and tropical circulations Summary oftropical motions and scales • There are still uncertainties concerning our knowledge about the interaction between convective and synoptic scales in the Tropics. • Horizontal temperature fluctuations in the Tropics are small <1K/1000 km; and in the absence of precipitation the vertical motions(subsidence) tend to balance the cooling through IR radiation loss: w dθ/dz = dθ/dt_rad = -1-2 K/day => w ~ -.5 cm/s • In the absence of condensation heating, tropical motions must be barotropic and cannot convert PE in KE. Therefore they must be driven by precipitating disturbances or lateral coupling with midlatitude systems. • When precipitation takes place, heating rates are strong; e.g. 100 mm/day precip ~ energy flux of 2900 W/m2 or an average 30 K/day heating of the atmospheric column => w ~ 8.6 cm/s. However, this positive mean motion is composed of strong ascent of order w ~ 1 m/s in the Cumulus updrafts and slow descending motion around (“compensating subsidence”) • when analysing the vorticity equation it appears that in precipitating disturbances the vertical transport of vorticity (momentum) through Cumulus is important to balance the divergence term

27. Midlatitude Convection (1) Convection associated to synoptic forcing, orographic uplift, and/or strong surface fluxes A Supercell over Central US, Mai 1998, flight level 11800 m

28. Midlatitude Convection (2)It’s raining again… Europe climatology (Frei and Schär, 1998) In Europe most intense precipitation is associated with orography, especially around the Mediterranean, associated with strong large-scale forcing and mesoscale convective systems

29. Midlatitude Convection (3)European MCSs (Morel and Sénési, 2001)Density Map of Triggering ….. over Orography

30. Midlatitude Convection (4)European MCSs (Morel and Sénési, 2001)Time of Trigger and mean propagation European (midlatitude) MCSs essentially form over orography (convective inhibition –see later- offset by uplift) and then propagate with the midtropospheric flow (from SW to NE)

31. Midlatitude Convection (5) along the main cold frontal band and in the cold core of the main depression – 17/02/97 during FASTEX A Supercell over Central US, Mai 1998, flight level 11800 m

32. Midlatitude Convection (6) Forcing of ageostrophic circulations/convection in the right entrance and left exit side of upper-level Jet Acceleration/deceleration of Jet Thermally indirect circulation Total energy is conserved: e.g. at the exit region where the Jet decelerates kinetic energy is converted in potential energy Thermally direct circulation

33. Midlatitude Convection (7) Conceptual model of a Squall line system with a trailing stratiform area (from Houze et al. 1989) • Evaporation of precipitation creates negatively buoyant air parcels. This can lead to the generation of convective-scale penetrative downdraughts. • In the stratiform part there is heating/cooling couple with an upper-level mesoscale ascent, and a lower-level mesoscale downdraught, due to the inflow of dry environmental air and the evaporation of stratiform rain.

34. Forward Flank downdraft induced by evaporation of precipitation Rear Flank Downdraft induced by dynamic pressure perturbation: Interaction of updraft with shear vector of environment: The linear part of the dynamic pressure perturbation is proportional to the horizontal gradient of the vertical velocity perturbation (updraft) times the environmental shear vector Midlatitude Convection (8a) Conceptual model of a rotating mesoscale convective system – tornadic thunderstorm (from Lemon and Doswell, 1979)

35. Conversion of horizontal vorticity at surface frontal boundary in vertical vorticity by tilting in updraft Midlatitude Convection (8b) Origin and mechanism of generation of vertical vorticity A useful quantity in estimating the storm intensity is the “bulk” Richardson number R=CAPE/S2, where CAPE is the convective available energy (see later) and S is the difference between the mean wind vector at 500 and 925 hPa

36. Summary: What is convection doing, where does it occur • Convection transports heat, water vapor, momentum … and chemical constituents upwards …. Water vapor then condenses and falls out -> net convective heating/drying • Deep Convection (precipitating convection) stabilizes the environment, an approximate not necessarily complete picture is to consider it as reacting to the large-scale environment (e.g. tropical waves, mid-latitude frontal systems) =“quasi-equilibrium”; shallow convection redistributes the surface fluxes • The tropical atmosphere is in radiative(cooling) / convective(heating) equilibrium 2K/day cooling in lowest 15 km corresponds to about 5 mm/day precipitation. • The effect of convection (local heat source) is fundamentally different in the midlatitudes and the Tropics. In the Tropics the Rossby radius of deformation R=N H/f (N=Brunt Vaisala Freq, f=Coriolis parameter, H=tropopause height) is infinite, and therefore the effects are not locally bounded, but spread globally via gravity waves – “throwing a stone in a lake”

37. What we have not talked about • Organization of convection: Squall lines, Mesoscale convective systems, tropical superclusters, and the influence of vertical wind shear • The diurnal cycle of convection over land (see lecture Notes) Follow some Tools and Concepts !

38. Assume fluid to be in hydrostatic equlibrium Top Bottom Gravity Net Force: Acceleration: Buoyancy - physics of Archimedes (1) Body in a fluid Forces: Emanuel, 1994

39. Buoyancy (2) Vertical momentum equation: Neglect second order terms

40. = = = B - buoyancy acceleration Buoyancy (3)

41. Buoyancy acceleration: Often (but not always): Then Hence Buoyancy (4)Contributions Dry air:

42. Buoyancy (5)Contributions Cloudy air: effects of humidity and condensate need to be taken into account In general all 3 terms are important. 1 K perturbation in T is equivalent to 5 g/kg perturbation in water vapor or 3 g/kg in condensate

43. 15 P B 10 Z (km) 5 0.02 -0.04 -0.02 0 0.04 (ms-2) Non-hydrostat. Pressure gradient effects Physics: Vector field of the buoyancy pressure-gradient force for a uniformly buoyant parcel of finite dimensions in the x-z-plane. (Houze, 1993, Textbook) CRM analysis of the terms Guichard and Gregory

44. Example: Much larger than observed - what’s going on ? Convective Available Potential Energy (CAPE) Definition: CAPE represents the amount of potential energy of a parcel lifted to its level of neutral buoyancy. This energy can potentially be released as kinetic energy in convection.

45. Constant temperature Moist adiabat Dry adiabat Constant mixing ratio Constant pressure Thermodynamic diagrams Tephigram

46. LNB CAPE CIN LFC LCL Convection in thermodynamic diagrams (1)using Tephigram/Emagram Idealised Profile

47. Convection in thermodynamic diagrams (2)using equivalent Potential Temperature and saturated equivalent Potential Temperature GATE Sounding θ Θeis conserved during moist adiabatic ascent CAPE Θesat(T) Θe(T,q) Note that no CAPE is available for parcels ascending above 900 hPa and that the tropical atmosphere is stable above 600 hPa (θe increases) – downdrafts often originate at the minimum level of θe in the mid-troposphere.

48. Importance of choice of moist adiabat in CAPE calculations Reversible moist adiabat: Condensate remains in parcel at all time. Consequences: Water loading (gravity acting on condensate) Condensate needs to be heated - different heat capacity than dry air Phase transition from water to ice leads to extra heating Irreversible moist adiabat (Pseudo-adiabat): Condensate is removed from parcel instantly

49. Importance of choice of moist adiabat in CAPE calculations CAPE - reversible adiabat without freezing vs. irreversible adiabat Reversible CAPE much smaller, typically by a factor of 2 with respect to irreversible Emanuel, 1994

50. From high-resolution LES simulation (dx=dy=50 m) Vaillancourt, You, Grabowski, JAS 1997 Mixing and 3D flowsubcloud and cloud-layer Circulations