Prior to IFS cycle 47r1, the EFI (and SOT) for both CAPE and CAPE-shear were based on instantaneous CAPE (ShortName=cape; parameter ID=59) and CAPE-shear (ShortName=capes; parameter ID=228044) fields, in 6-hour steps, from which 24-hour maximum values were computed (more details can be found in ECMWF Newsletter 144). For example for the nominal T+24-48h EFI forecast, four CAPE values, at T+30h, T+36h, T+42h and T+48h, were used to compute the maximum for this forecast period. In IFS cycle 45r1 two new model output parameters, maximum CAPE in the last 6 hours (ShortName=mxcape6; parameter ID=228035) and maximum CAPE-shear in the last 6 hours (ShortName=mxcapes6; parameter ID=228036), are introduced (see details here). In IFS cycle 47r1, for the respective EFI and SOT computations, mxcape6 and mxcapes6 will be used instead of cape and capes respectively. This change is aiming for a better sampling in the computation of the 24-hour maxima needed for the EFI. In effect with the change we extract the maximum within 24 hourly values, instead of using 4 6-hourly values. The impacts of this change are described in detail below. 

Firstly this changes has a neutral impact in terms of EFI skill, as measured by the area under the Relative Operating Characteristic (ROC) curve.

The EFI fields using mxcape6 and mxcapes6 generally look smoother and in the case of a tornadic outbreak in the USA earlier in March 2020 shown in Fig. 1 they even fit better to the severe weather reports. Please note that the EFI shows a wide area where the environment favours convective hazards whilst the actual severe thunderstorms develop along relatively narrow bands where sufficient lift is present to initiate deep, moist convection, e.g. along cold fronts or dry lines. 

 

Fig. 1. T+24-48h CAPE-shear EFI/SOT forecasts for the tornadic outbreak on 3rd March 2020 over Tennessee, USA.

As a result of taking instantaneous cape and capes at discrete time steps the EFI can exhibit a stripy structure in the case of fast moving squall lines or very convective fronts as in the example in Fig. 2. This is an issue because the cold front is not jumping in reality, but is moving continuously from west to east. Replacing cape and capes with mxcape6 and mxcapes6 removes the stripy behaviour and provides a smoother and more realistic forecast field due to better sampling. Smaller scale discretization apparent on the top centre panel for HRES (attributable to the use of values at 1h intervals) is not manifested on the EFI because it will effectively be smoothed out by frontal-timing spread within the ensemble. The plot also provides one nice example of why we might expect to see more areally-extended regions of high EFI/SOT with the new formulation (and why these should be more realistic).

Fig. 2. EFI (and SOT) for CAPE-shear as well as the CAPE-shear high-resolution forecast for a case of a fast moving cold front. Stripy fields of instantaneous CAPE-shear and EFI/SOT (see black lines) are due to taking instantaneous CAPE-shear values at discrete time steps. The EFI forecast looks smoother when using the smoother fields of maximum CAPE-shear in 6-hour periods. The sequence of Air mass RGB imagery shows snapshots of the cold front approaching the south-western parts of the Iberian Peninsular.

The EFI is a climate-related product which provides information on how different are the cumulative distribution functions (CDFs) of the ensemble forecast and the model climate. So an anomalous forecast does not necessarily translate into a high-impact one. In the case of convection this means that anomalous values of CAPE do not necessarily suggest that convective hazards are likely. For example, climatological values of CAPE in the cold season over the Arctic are so low that severe convective hazards are impossible even if the forecast shows up "extreme" CAPE. To filter out anomalous but insignificant signals in the EFI, CAPE values less than 10 J.kg-1 are set to 0 before computing both the model climate and the ensemble forecast CDFs. With cape replaced by mxcape6 and capes with mxcapes6, CAPE over 10 J.kg-1 will happen more often and therefore this may show up on the EFI (Fig. 3). To correctly interpret the EFI one should always account for the model climate as well.

  

Fig. 3. EFI/SOT for CAPE and max CAPE* in Siberia, illustrating that the latter may show up a bit more often.

Finally, as a result of this change, the model climate also changes slightly across the world (Fig. 4a). In the example in Fig. 4b the EFI suggests a bit more extreme CAPE values over NW Spain. The CDFs for a representative location there show what this signal corresponds to in terms of model climate and ensemble forecast distributions.

a) b)

Fig. 4 a) Model climate 99th percentile for cape and mxcape6* valid for the end of April and b) an example of the EFI/SOT for cape and mxcape* and CDFs of corresponding model climate and real-time ensemble forecasts.
 

*Naming Conventions: We refer on this web page to "max CAPE" (or mxcape6) and "max CAPE-shear" (or mxcapes6) for distinguishing the EFI/SOT of these in cycles from 47r1, from the equivalent EFI/SOT fields used up to and including cycle 46r1, but in fact when the model changes the former literally replaces the latter. Therefore users do not see any change in the naming - e.g. in ecCharts - and there will not be any additional EFI/SOT; instead the EFI and SOT for CAPE and CAPE-shear will continue to exist but they will use instead mxcape6 and mxcapes6.

5 Comments

  1. I would like to know why 925hPA an not 10 m was chosen in CAPES ID 228044 and I think also
    in the computation mxcapes6. In literature the DLS 0-6 km is most used.


    1. Dear Alessio, 10m wind magnitude is strongly affected by the roughness of the underlying surface. But what we want for CAPES is a metric of the innate shear within the low-to-mid troposphere. So whilst over relatively flat grasslands (say) 10m winds can be representative of the lower troposphere, and could be OK for a CAPES calculation, in other areas this is not the case. This is particularly true over mountain where the IFS 10m wind speeds tend to be very light, and very probably too light. So if we were to use 10m winds in CAPES one would see disproportionately large values, on average, over mountains, because of the artificial inflation of the shear component there by friction. Which in turn would have meant misleading guidance for forecasters. So we needed an alternative. We chose the standard atmospheric level of 925hPa. This has the advantages of: (i) being a lower tropospheric level; (ii) not having the roughness-related disadvantages of 10m winds; (iii) being a level archived from all the ensemble members, which helps with reproducibility and any post-event analysis and (iv) also being a standard level reported from radiosonde ascents. There are some minor disadvantages too but on balance this was deemed the best option available. 

      1. Dear Timothy,
        thank you for your kind reply, but I don't agree with your arguments.
        Organized storms, supercell or tornadic supercells occur only in the plains and not in the mountains, because roughness/obstacles disturb upflow and downflow. Forecaster don't care about high values.
        925hPA corresponds approx 800m and this is to high for plains, and this level is also disturbed over mountains.
        CAPES is also like an index, it would be better to consider the two separate quantities: CAPE and shear.

        1. Dear Alessio,

          Thank you very much for your comments. We started computing CAPE-shear composite parameter (CAPES) few years ago in order to use this quantity for the Extreme Forecast Index (EFI) and provide forecasters with suitable ensemble-based guidance on organised severe convection.  We used 925-500 hPa wind difference as one to represent deep-layer wind shear. On the other hand, I agree that forecasters use 0-6 km (or 10 m to 500 hPa) wind vector difference to assess the deep-layer shear. In the ecCharts we provide the wind shear computed for various layers that one can select themselves. I noticed that we have 100 m level as the lowest one so far – we can add 10m as well. In the ecCharts one can select the wind shear for any layer from the option menu and add CAPE for example as a separate layer; therefore one can look how CAPE and shear superimpose when forecasting deep, moist convection. In the EFI we use CAPES in a relative sense with respect to model climatology and there is no evidence that 925-500 shear changes considerably the EFI signal compared to 0-6 km shear. With European Severe Storms Laboratory (ESSL) we are currently preparing an overview of different shear parameters for forecasting severe convective storms. The choice of 925 hPa pressure level is somewhat arbitrary and it can always be revised if necessary but I haven't read any study to point out that 10 m is better or worse than any other  level in the lower troposphere when considering deep-layer shear.

          1. Thank you for adding the layer 10 m in the ecCharts updates - November 2020.