Convection-permitting ensembles have led to improved forecasts of many atmospheric phenomena. However, to fully utilize these forecasts the dependence of predictability on synoptic conditions needs to be understood. In this study, convective regimes are diagnosed based on a convective time scale that identifies the degree to which convection is in equilibrium with the large-scale forcing. Six convective cases are examined in a convection-permitting ensemble constructed using the Met Office Unified Model. The ensemble members were generated using small-amplitude buoyancy perturbations added into the boundary layer, which can be considered to represent turbulent fluctuations close to the grid scale. Perturbation growth is shown to occur on different scales with an order of magnitude difference between the regimes [O(1) km for cases closer to nonequilibrium convection and O(10) km for cases closer to equilibrium convection]. This difference reflects the fact that cell locations are essentially random in the equilibrium events after the first 12 h of the forecast, indicating a more rapid upscale perturbation growth compared to the nonequilibrium events. Furthermore, large temporal variability is exhibited in all perturbation growth diagnostics for the nonequilibrium regime. Two boundary condition–driven cases are also considered and show similar characteristics to the nonequilibrium cases, implying that caution is needed to interpret the time scale when initiation is not within the domain. Further understanding of perturbation growth within the different regimes could lead to a better understanding of where ensemble design improvements can be made beyond increasing the model resolution and could improve interpretation of forecasts.

Recent surface-water and flash floods have caused millions of pounds worth of damage in the UK. These events form rapidly and are difficult to predict due to their short-lived and localised nature. The interdisciplinary Flooding From Intense Rainfall (FFIR) programme investigated the feasibility of enhancing the integration of an end-to-end forecasting system for flash and surface-water floods to help increase the lead time for warnings for these events. Here we propose developments to the integration of an operational end-to-end forecasting system based on the findings of the FFIR programme. The suggested developments include methods to improve radar-derived rainfall rates and understanding of the uncertainty in the position of intense rainfall in weather forecasts; the addition of hydraulic modelling components; and novel education techniques to help lead to effective dissemination of flood warnings. We make recommendations for future advances such as research into the propagation of uncertainty throughout the forecast chain. We further propose the creation of closer bonds to the end users to allow for an improved, integrated, end-to-end forecasting system that is easily accessible for users and end users alike, and will ultimately help mitigate the impacts of flooding from intense rainfall by informed and timely action.

The FRANC project (Forecasting Rainfall exploiting new data Assimilation techniques and Novel observations of Convection) has researched improvements in numerical weather prediction of convective rainfall via the reduction of initial condition uncertainty. This article provides an overview of the project’s achievements. We highlight new radar techniques: correcting for attenuation of the radar return; correction for beams that are over 90% blocked by trees or towers close to the radar; and direct assimilation of radar reflectivity and refractivity. We discuss the treatment of uncertainty in data assimilation: new methods for estimation of observation uncertainties with novel applications to Doppler radar winds, Atmospheric Motion Vectors, and satellite radiances; a new algorithm for implementation of spatially-correlated observation error statistics in operational data assimilation; and innovative treatment of moist processes in the background error covariance model. We present results indicating a link between the spatial predictability of convection and convective regimes, with potential to allow improved forecast interpretation. The research was carried out as a partnership between University researchers and the Met Office (UK). We discuss the benefits of this approach and the impact of our research, which has helped to improve operational forecasts for convective rainfall events.

We present a simple, physically consistent stochastic boundary layer scheme implemented in the Met Office’s Unified Model. It is expressed as temporally correlated multiplicative Poisson noise with a distribution that depends on physical scales. The distribution can be highly skewed at convection-permitting scales (horizontal grid lengths around 1 km) when temporal correlation is far more important than spatial. The scheme is evaluated using small ensemble forecasts of two case studies of severe convective storms over the UK. Perturbations are temporally correlated over an eddy-turnover timescale, and may be similar in magnitude to or larger than the mean boundary-layer forcing. However, their mean is zero and hence they, in practice, they have very little impact on the energetics of the forecast, so overall domain-averaged precipitation, for example, is essentially unchanged. Differences between ensemble members grow; after around 12 h they appear to be roughly saturated; this represents the time scale to achieve a balance between addition of new perturbations, perturbation growth and dissipation, not just saturation of initial perturbations. The scheme takes into account the area chosen to average over, and results are insensitive to this area at least where this remains within an order of magnitude of the grid scale.

Convection-permitting modelling has led to a step change in forecasting convective events. However, convection occurs within different regimes which exhibit different forecast behaviour. A convective adjustment timescale can be used to distinguish between these regimes and examine their associated predictability. The convective adjustment timescale is calculated from radiosonde ascents and found to be consistent with that derived from convection-permitting model forecasts. The model-derived convective adjustment timescale is then examined for three summers in the British Isles to determine characteristics of the convective regimes for this maritime region. Convection in the British Isles is predominantly in convective quasi-equilibrium, with 85% of convection having a timescale less than or equal to 3 h. This percentage varies spatially with more non-equilibrium events occurring in the south and southwest. The convective adjustment timescale exhibits a diurnal cycle over land. The non-equilibrium regime occurs more frequently at mid-range wind speeds and with winds from southerly to westerly sectors. Most non-equilibrium convective events in the British Isles are initiated near large coastal orographic gradients or on the European continent. Thus, the convective adjustment timescale is greatest when the location being examined is immediately downstream of large orographic gradients and decreases with distance from the convective initiation region. The dominance of convective quasi-equilibrium conditions over the British Isles argues for the use of large-member ensembles in probabilistic forecasts for this region.

Abstract. The dynamical and microphysical properties of a well-observed cyclone from the North Atlantic Waveguide and Downstream Impact Experiment (NAWDEX), called the Stalactite cyclone and corresponding to intensive observation period 6, is examined using two atmospheric components (ARPEGE-Climat 6.3 and LMDZ6A) of the global climate models CNRM-CM6-1 and IPSL-CM6A, respectively. The hindcasts are performed in “weather forecast mode”, run at approximately 150–200 km (low resolution, LR) and approximately 50 km (high resolution, HR) grid spacings, and initialised during the initiation stage of the cyclone. Cyclogenesis results from the merging of two relative vorticity maxima at low levels: one associated with a diabatic Rossby vortex (DRV) and the other initiated by baroclinic interaction with a pre-existing upper-level potential vorticity (PV) cut-off. All hindcasts produce (to some extent) a DRV. However, the second vorticity maximum is almost absent in LR hindcasts because of an underestimated upper-level PV cut-off. The evolution of the cyclone is examined via the quasi-geostrophic ω equation which separates the diabatic heating component from the dynamical one. In contrast to some previous studies, there is no change in the relative importance of diabatic heating with increased resolution. The analysis shows that LMDZ6A produces stronger diabatic heating compared to ARPEGE-Climat 6.3. Hindcasts initialised during the mature stage of the cyclone are compared with airborne remote-sensing measurements. There is an underestimation of the ice water content in the model compared to the one retrieved from radar-lidar measurements. Consistent with the increased heating rate in LMDZ6A compared to ARPEGE-Climat 6.3, the sum of liquid and ice water contents is higher in LMDZ6A than ARPEGE-Climat 6.3 and, in that sense, LMDZ6A is closer to the observations. However, LMDZ6A strongly overestimates the fraction of super-cooled liquid compared to the observations by a factor of approximately 50.

Case studies remain an important method for meteorological parameter sensitivity process studies. However, these types of study often use just a few case studies (typically up to three) and are not tested for statistical significance. This approach can be problematic at the convective scales, since uncertainty in the representation of an event increases, and the variability in the atmosphere arising from convective‐scale noise is not routinely taken into account. Here we propose a simple ensemble method for performing more robust sensitivity analysis without the need for an operational‐style ensemble prediction system and demonstrate it using a case study from the 2005 Convective Storm Initiation Project. Boundary‐layer stochastic potential temperature perturbations with Gaussian spatial structure are used to create small ensembles to examine the impact of increasing cloud droplet number concentration (CDNC) on precipitation. Whilst there is a systematic difference between the experiments, such that increasing the CDNC reduces the precipitation, there is also an overlap between the different ensembles implying that convective‐scale variability should be taken into account in case study process‐based sensitivity studies.

Convection-permitting forecasts have improved the forecasts of flooding from intense rainfall. However, probabilistic forecasts, generally based upon ensemble methods, are essential to quantify forecast uncertainty. This leads to a need to understand how different aspects of the model system affect forecast behaviour. We compare the uncertainty due to initial and boundary condition (IBC) perturbations and boundary-layer turbulence using a super ensemble (SE) created to determine the influence of 12 IBC perturbations vs. 12 stochastic boundary-layer (SBL) perturbations constructed using a physically-based SBL scheme. We consider two mesoscale extreme precipitation events. For each we run a 144–member SE. The SEs are analysed to consider the growth of differences between the simulations, and the spatial structure and scales of those differences. The SBL perturbations rapidly spin-up, typically within 12 h of precipitation commencing. The SBL perturbations eventually produce spread that is not statistically different from the spread produced by the IBC perturbations, though in one case there is initially increased spread from the IBC perturbations. Spatially, the growth from IBC occurs on larger scales than that produced by the SBL perturbations (typically by an order of magnitude). However, analysis across multiple scales shows that the SBL scheme produces a random relocation of precipitation up to the scale at which the ensemble members agree with each other. This implies that statistical post-processing can be used instead of running larger ensembles. Use of these statistical post-processing techniques could lead to more reliable probabilistic forecasts of convective events and their associated hazards.