You are viewing an old version of this page. View the current version.

Compare with Current View Page History

« Previous Version 49 Next »

Introduction

This page describes two studies of different convective cases; one over N. America associated with formation of severe tornadoes, the other over central Africa. The N.America case has strong large scale forcing whereas the central African case is driven by the diurnal cycle.

Both cases are studied starting the forecast from the same date/time (initial conditions).

In these case studies, you will carry out a control forecast, followed by any number of suggested sensitivities experiments.

US Tornado convection case (Arkansas)

On the 27 April 7pm local time (00UTC 28 April), tornadoes hit towns north and west of Little Rock, Arkansas killing approx 17 people.  (http://edition.cnn.com/2014/04/28/us/severe-weather/index.html?hpt=hp_c2). On the evening on the 28 April fatal tornadoes occurred over Mississippi (http://www.bbc.co.uk/news/world-us-canada-27199071).

This case study will look at the role of convection and the large scale in these events.

More information can also be found on the ECMWF Severe Event Catalogue 201404 - Convection - Arkansas U.S.

Tornado damage

African diurnal deep convection (Central Africa)

Over tropical land masses, incoming radiation strongly heats the surface leading to the development of deep convection and precipitation. Observations show that convective activity and precipitation peak in the late afternoon or early evening. Until very recently numerical weather prediction models struggled to reproduce this diurnal cycle, often predicting convection to peak too early in the day. In this case study, aspects of the convective parameterization scheme can be altered to see how the intensity and the diurnal cycle of convection responses.

On this page...

 

 

Initial conditions

Case study initial conditions for this case study are provided on the OpenIFS ftp site.

The initial conditions are available at a range of different resolutions and start dates for a 30hr forecast. The experiment ids are created at ECMWF and used for identifying the model forecasts on the ECMWF archive system (for those with access).

Note that ERA-Interim has a resolution of T255.

As OpenIFS is a spectral model, the 'T' number refers to the triangular truncation in spectral space. Equivalent grid resolutions are:
T159 ~ 125km resolution, T255 ~ 80km, T511 ~ 40km, T799 ~ 25km, T1279 ~ 16km.

TODO: table of data to download

Download instructions

Example using T159
% mkdir -p runs/convection/t255
% cd runs
% ftp ftp.ecmwf.int
ftp> cd case_studies/lothar_storm
ftp> binary
ftp> get 1999122412_T159_fqar.tgz
ftp> quit
% tar zxf 1999122412_T159_fqar.tgz
% ls
1999122412_T159.tgz  ICMCLfqarINIT  ICMGGfqarINIT  ICMGGfqarINIUA  ICMSHfqarINIT  ecmwf
% ls ecmwf
NODE.001_01  namelistfc

The 'ecmwf' directory contains the files produced at ECMWF when this experiment was run:

  • namelistfc : copy this file to 'fort.4' to run the experiment (modify as required)
  • NODE.001_01 : this is the model output file as run at ECMWF. If your run fails, it may be useful to compare with this file.

Run the control forecast

The first step is to run the control forecast. Both cases can be studied with a single forecast.

See below for tasks and key questions to address for the control forecast before moving on to the sensitivity experiments.

Case study: N.America deep convection

On 27 April 2014 7pm local time (00UTC 28 April), tornadoes hit towns north and west of Little Rock, Arkansas.



Key questions and tasks using the control forecast


  1. Understand the weather situation resulting in tornadoes
  2. Evaluate the control forecast and compare it to the ECMWF reanalysis and observations
  3. What is the area of threat according to the control forecast? Area of threat = the area where severe weather can expected. This can be identified by considering parameters such as CAPE, CIN, 850-hPa equivalent potential temperature. 
  4. How does the convective adjustment process takes place and and what is the role of large scale forcing (why and where it happens)?

 

Case study: African deep convection

 

 

Key questions and tasks using the control forecast

 

  1. Understand the weather situation over Africa.
  2. What is the role of large scale in this case (compare with N.America case).
  3. Look at the diurnal variation of key parameters (2m temperature, surface heat fluxes, precipitation, outgoing-longwave-radiation) for location 0N,25E.
  4. Compare differences in convection profiles between Central Africa and (i) open ocean, and (ii) Amazonia.

 

Sensitivity experiments

The IFS is highly tuned to give the best forecast over a range of initial conditions. However, it is instructive to try some sensitivity experiments to understand the role of various physical and dynamical processes.

Not all of the suggested experiments are applicable to both cases, indicated in brackets.

  • What's the impact of the different 'lead times' on the forecast of the convection (i.e. starting from different dates)? (N.America only)
  • What's the impact of resolution on the forecast of the convection? (both)
  • Does reducing the model timestep improve or worsen the forecast? (both)
  • Turn off deep convection (both)

    Do this by editing fort.4, find the namelist block NAMCUMF and add a line:

    LMFPEN=false,            ! disable deep convection

  • Impact of the improved diurnal cycle of convection. (Africa only)
    In this sensitivity experiment, look at the timing of convective and precipitation events by changing how the model parametrizes the diurnal cycle.

    OpenIFS has 3 options for the controlling the diurnal cycle. To change between them:

    - Edit the fort.4 file

    - Find the namelist NAMCUMF and change the parameter RCAPDCYCL accordingly:

    RCAPDCYCL = 2 (default) activates the diurnal cycle using sub-cloud CAPE,

    RCAPDCYCL = 1 diurnal cycle using surface sensible heat flux,

    RCAPDCYCL = 0 reverts the code to a setting before the diurnal cycle for convection was implemented.

  • Increase the precipitation auto conversion rate - what impact does this have? (both)

    Edit the source code to increase the auto conversion rate by 20%

    File: ifs/phys_ec/sucldp.F90, change:

    line 123: RKCONV=1._JPRB/6000._JPRB   ! 1/autoconversion time scale (s)

    to:

    line 123: ! RKCONV=1._JPRB/6000._JPRB   ! 1/autoconversion time scale (s)
    line 124: RKCONV=1.2_JPRB/6000._JPRB    ! default scaled by 20%: 1/autoconversion time scale (s)
  • Impact of the convective time scale adjustment (both)
    An optimization factor in the parametrization is used for tuning the diurnal cycle. This can be altered by changing a value in the model namelist.

    To change the timescale:

    - Edit the fort.4 file

    - Find the namelist  NAMCUMF, parameter RTAUA.

    - The default value is RTAUA=1.

    - Run two sensitivity experiments with values of RTAUA = 0.33 and 3.

    The ratio between the actual cloud base mass flux and the unit (initial) cloud base mass flux:

    \[ \frac{M_{base}}{M^*_{base}} = \frac{PCAPE - PCAPE_{bl}}{\tau} \]

    Look at the amplitude of precipitation.

  • Sensitivity to entrainment rate (both)

    To change the entrainment rate:

    - Edit the fort.4 file

    - Find the namelist block NAMCUMF, parameter ENTRORG

    - The default value is ENTRORG= 1.75E-3

    ENTRORG= 5.8E-4 reduced by factor 3 (mostly shallow convection regime)

    ENTRORG= 5.25E-3 multiplied by factor 3 (mostly deep convection regime)

    Look at the cloud top height, precipitation and eventually changes in temperature and moisture fields with respect to the reference. Note also this is having less impact with the diurnal cycle activated.

Additional questions

  • How important is the correct diurnal cycle of precipitation and radiation for 2m temperature and dewpoint forecast?

Further reading

 

Comments

The forecasting system at ECMWF makes use of "ensembles" of forecasts to account for errors in the initial state. In reality, the forecast depends on the initial state in a much more complex way than just the model resolution or starting date.  At ECMWF many initial states are created for the same starting time by use of "singular vectors" and "ensemble data assimilation" techniques which change the vertical structure of the initial perturbations.

As further reading and an extension of this case study, research how these methods work.

Acknowledgements

We are grateful to: Peter Bechtold, Filip Vana, Sandor Kertesz in preparing the material for the OpenIFS user workshop in Stockholm 2014, from which most of the material on this page is derived. We also thank the forecast department for their material on the ECMWF Severe Event Catalogue that was used in preparing these cases.

 

 


 

  • No labels