Part IV: Physical Processes
Total Page：16
File Type：pdf, Size：1020Kb
Part IV: Physical Processes IFS DOCUMENTATION { Cy45r1 Operational implementation 5 June 2018 PART IV: PHYSICAL PROCESSES Table of contents Chapter 1 Overview Chapter 2 Radiation Chapter 3 Turbulent transport and interactions with the surface Chapter 4 Subgrid-scale orographic drag Chapter 5 Non-orographic gravity wave drag Chapter 6 Convection Chapter 7 Clouds and large-scale precipitation Chapter 8 Surface parametrization Chapter 9 Methane oxidation Chapter 10 Ozone chemistry parametrization Chapter 11 Climatological data Chapter 12 Basic physical constants and thermodynamic functions References c Copyright 2018 European Centre for Medium-Range Weather Forecasts Shinfield Park, Reading, RG2 9AX, England Literary and scientific copyrights belong to ECMWF and are reserved in all countries. This publication is not to be reprinted or translated in whole or in part without the written permission of the Director. Appropriate non-commercial use will normally be granted under the condition that reference is made to ECMWF. The information within this publication is given in good faith and considered to be true, but ECMWF accepts no liability for error, omission and for loss or damage arising from its use. IFS Documentation { Cy45r1 1 TABLE OF CONTENTS REVISION HISTORY Changes since CY43R3 • Chapter 6. Updates to parcel perturbations, phase of shallow clouds, positiveness of CAPE closure denominator, typos in previous versions. Added description of lightning parametrization. • Chapter 8. Updates for HRES now being coupled to NEMO. Changes since CY43R1 • Chapter 2. Aerosol climatological revision. • Chapter 6. Updates to sections on Large-scale budget equations, Freezing in convective updraughts, Generation of precipitation, Fallout of precipitation, Melting and freezing of precipitation, Link to the cloud scheme. • Chapter 7. New section on Cloud height. • Chapter 11.Climatological data plots have been updated to new cubic octahedral grid and text modified accordingly. Changes since CY41R2 • Convection chapter updates in sections: Detrainment rates; Deep convection; Freezing in convective updraughts; Generation of precipitation; Melting and freezing of precipitation; Temporal discretization; Diagnostics for postprocessing. • Climatological data chapter updates in section: Ozone • Surface parametrization chapter updates in sections: Sea Ice and Ocean Boundary Conditions. • Clouds chapter updates: New section on 'Description of output fields. Addition of 3 new figures. • Turbulent transport chapter updates: Removal of reference to ice in cloud description in section Outer Layer; section on Description of output fields changed to Diagnostic computations for post- processing and added to; corrections to equations in section Solution of the EDMF equations. 2 IFS Documentation { Cy45r1 Part IV: Physical Processes Chapter 1 Overview Table of contents 1.1 Introduction 1.2 Overview of the code 1.1 INTRODUCTION The physical processes associated with radiative transfer, turbulent mixing, convection, clouds, surface exchange, subgrid-scale orographic drag and non-orographic gravity wave drag have a strong impact on the large scale flow of the atmosphere. However, these mechanisms are often active at scales smaller than the resolved scales of the model grid. Parametrization schemes are then necessary in order to properly describe the impact of these subgrid-scale mechanisms on the large scale flow of the atmosphere. In other words the ensemble effect of the subgrid-scale processes has to be formulated in terms of the resolved grid- scale variables. Furthermore, forecast weather parameters, such as two-metre temperature, precipitation and cloud cover, are computed by the physical parametrization part of the model. This part (PartIV `Physical processes') of the IFS documentation describes the physical parametrization package. The radiation scheme (Chapter2) performs computations of the short-wave and long-wave radiative fluxes using the predicted values of temperature, humidity, cloud, and monthly-mean climatologies for aerosols and the main trace gases (CO2,O3, CH4,N2O, CFCl3 and CF2Cl2). The radiation code is based on the Rapid Radiation Transfer Model (RRTM, Mlawer et al., 1997; Iacono et al., 2008). Cloud- radiation interactions are taken into account in detail by using the values of cloud fraction and liquid, ice and snow water contents from the cloud scheme using the McICA (Monte Carlo Independent Column Approximation) method (McRad, Morcrette et al., 2008a). The solution of the radiative transfer equations Non-orographic wave drag O3 Chemistry Long-wave Short-wave CH4 Oxidation radiation radiation Cloud Cloud Subgrid-scale orographic drag Deep convection Shallow Turbulent diffusion convection Sensible Latent heat flux Long-wave Short-wave heat Wind Waves flux flux flux Surface Ocean model Figure 1.1 Schematic diagram of the different physical processes represented in the IFS model. IFS Documentation { Cy45r1 3 Chapter 1: Overview to obtain the fluxes is computationally very expensive, so depending on the model configuration, full radiation calculations are performed on a reduced (coarser) radiation grid and/or on a reduced time frequency. The results are then interpolated back to the original grid. However, the short-wave fluxes are updated at every grid point and timestep using values of the solar zenith angle. The turbulent diffusion scheme (Chapter3) represents the vertical exchange of heat, momentum and moisture through sub-grid scale turbulence. The vertical turbulent transport is treated differently in the surface layer and above. In the surface layer, the turbulence fluxes are computed using a first order K-diffusion closure based on the Monin{Obukhov (MO) similarity theory. Above the surface layer a K-diffusion turbulence closure is used everywhere, except for unstable boundary layers where an Eddy- Diffusivity Mass-Flux (EDMF) framework is applied, to represent the non-local boundary layer eddy fluxes (K¨ohler et al., 2011). The scheme is written in moist conserved variables (liquid static energy and total water) and predicts total water variance. A total water distribution function is used to convert from the moist conserved variables to the prognostic cloud variables (liquid/ice water content and cloud fraction), but only for the treatment of stratocumulus. Convective clouds are treated separately by the shallow convection scheme. The effects of unresolved orography on the atmospheric flow are parametrized as a sink of momentum (drag). The turbulent diffusion scheme includes a parametrization in the lower atmosphere to represent the turbulent orographic form drag induced by small scale (< 5 km) orography (Beljaars et al., 2004b). In addition, in stably stratified flow, the orographic drag parametrization (Chapter4) represents the effects of low-level blocking due to unresolved orography (blocked flow drag) and the absorption and/or reflection of vertically propagating gravity waves (gravity wave drag) on the momentum budget (Lott and Miller, 1997). The non-orographic gravity wave drag parametrization (Chapter5) accounts for the effects of unresolved non-orographic gravity waves. These waves are generated in nature by processes like deep convection, frontal disturbances, and shear zones. Propagating upward from the troposphere the waves break in the middle atmosphere, comprising the stratosphere and the mesosphere, where they exert a strong drag on the flow. The parametrization uses a globally uniform wave spectrum, and propagates it vertically through changing horizontal winds and air density, thereby representing the wave breaking effects due to critical level filtering and non-linear dissipation. A description of the scheme and its effects on the middle atmosphere circulation can be found in Orr et al. (2010). The moist convection scheme (Chapter6) is based on the mass-flux approach and represents deep (including congestus), shallow and mid-level (elevated moist layers) convection. The distinction between deep and shallow convection is made on the the basis of the cloud depth (≤ 200 hPa for shallow). For deep convection the mass-flux is determined by assuming that convection removes Convective Available Potential Energy (CAPE) over a given time scale. The intensity of shallow convection is based on the budget of the moist static energy, i.e. the convective flux at cloud base equals the contribution of all other physical processes when integrated over the subcloud layer. Finally, mid-level convection can occur for elevated moist layers, and its mass flux is set according to the large-scale vertical velocity. The scheme, originally described in Tiedtke(1989), has evolved over time and amongst many changes includes a modified entrainment formulation leading to an improved representation of tropical variability of convection (Bechtold et al., 2008), and a modified CAPE closure leading to a significantly improved diurnal cycle of convection (Bechtold et al., 2014). Clouds and large-scale precipitation (Chapter7) are parametrized with a number of prognostic equations for cloud liquid, cloud ice, rain and snow water contents and a sub-grid fractional cloud cover. The cloud scheme represents the sources and sinks of cloud and precipitation due to the major generation and destruction processes, including cloud formation by detrainment from cumulus convection, condensation, deposition, evaporation, collection, melting and freezing. The scheme is based on (Tiedtke, 1993) but with an enhanced representation of mixed-phase clouds and prognostic precipitation (Forbes and Tompkins,