Session 2 | Numerical modelling of storms, storm-scale data assimilation

Session 2

Numerical modelling of storms, storm-scale data assimilation
Orals
| Thu, 11 May, 17:00–18:00 (EEST)|Main Conference Room
Posters
| Attendance Thu, 11 May, 14:30–16:00 (EEST) | Display Wed, 10 May, 09:00–Thu, 11 May, 18:30|Exhibition area
Orals |
Thu, 17:00
Thu, 14:30

Orals: Thu, 11 May | Main Conference Room

17:00–17:15
|
ECSS2023-33
Killian P. Brennan, Michael Sprenger, Marco Arpagaus, André Walser, and Heini Wernli
On the 28th of June 2021, a supercell traversed the complex topography of Switzerland. While propagating over a distance of 350 km and lasting 8 hours, the cell produced severe hail with diameters up to 9 cm. A case study is performed, aimed at understanding the development of the supercell environment, and the physical processes associated with the cell’s initiation, intensification, propagation, and dissipation. To this end, ensemble hindcast simulations are conducted with a convection-resolving numerical weather prediction model (COSMO-1E) coupled to a one-dimensional hail growth parameterization (HAILCAST). The simulation setup is identical to the operational ensemble forecasts performed at MeteoSwiss. Object-based analysis of the simulated hail cell is facilitated by a tracking algorithm, developed specifically to address the many challenges associated with tracking hail storms in output from numerical simulations. The hail cell tracks enable a storm-relative analysis of the storm environment and the temporal evolution of the storm structure during its phases of development. The inflow of air feeding the storm updraft is investigated with the aid of Lagrangian air parcel trajectories initiated along the hail track, giving novel insights into the low-level flow and moisture sources of the storm. All ensemble members produce thunderstorms comparable to the observed storm in terms of track position and intensity. However, the observations and the simulations differ with respect to initiation location and the lifetime of the cell. The analysis focuses on the intensification and dissipation stages of the storm, and especially how and why the cell tracks in the ensemble members differ from each other.

How to cite: Brennan, K. P., Sprenger, M., Arpagaus, M., Walser, A., and Wernli, H.: A modeling case study of a severe hail storm in complex topography, 11th European Conference on Severe Storms, Bucharest, Romania, 8–12 May 2023, ECSS2023-33, https://doi.org/10.5194/ecss2023-33, 2023.

17:15–17:30
|
ECSS2023-142
Cloé David, Clotilde Augros, Benoît Vié, and François Bouttier

Thunderstorms are among the most devastating weather phenomena, mostly because of hail, heavy rainfall and strong winds. Although significant progress in convective scale numerical weather prediction models have been done these last years, the forecast of thunderstorms from a few hours to a few days remains difficult. It is partly because of the low predictability of these small scale phenomena but also due to deficiencies in the models, in particular in their microphysical schemes that describe the evolution of hydrometeor concentrations and sizes. Dual-polarization radar observations are very useful for microphysical schemes evaluation (Ryzhkov et al., 2020) as they provide information related to the size, shape, composition and orientation of the hydrometeors.

Recent studies (Kumjian and Ryzhkov, 2008; Kumjian et al., 2014; Johnson et al., 2016; Snyder et al., 2017a, b) have shown that supercell thunderstorms are associated with recurrent distinctive polarimetric signatures, such as the ZDR columns that consist of quasi-vertical continuous columns of enhanced ZDR extending above the environmental 0°C level and provide information about the location and intensity of a storm’s updraft (Snyder et al., 2015).

This work aims to statistically evaluate the ability of the AROME convective scale model (Brousseau et al., 2016) to replicate this signature when associated with different microphysics schemes : ICE3 (Pinty and Jabouille, 1998, one-moment) or LIMA (Vié, 2016; two-moment for rain, cloud water and ice cristals) which both can explicitely predict hail.

To do so, an automated detection of ZDR columns will be first implemented following Snyder et al. (2015) and Starzec et al. (2017). The algorithm will be applied on a few convective cases that occurred in France in 2022. Simulations with the AROME model will be performed on the same cases with different microphysics options. Simulated ZDR columns will be obtained thanks to the Augros et al. (2016) polarimetric radar forward operator applied to the AROME simulations. Finally, the characteristics (width, height) of the observed and simulated ZDR columns will be compared to assess which microphysics scheme is able to best reproduce the observed signatures. Besides evaluating the model microphysics, this study will enable to better understand the discrepancies between simulated and observed ZDR columns, a necessary first step towards using theses signatures in an assimilation context to improve storm-scale analysis and forecast.

How to cite: David, C., Augros, C., Vié, B., and Bouttier, F.: Statistical comparisons between observed and simulated ZDR columns using AROME model with different microphysics schemes, 11th European Conference on Severe Storms, Bucharest, Romania, 8–12 May 2023, ECSS2023-142, https://doi.org/10.5194/ecss2023-142, 2023.

17:30–17:45
|
ECSS2023-90
Valerian Jewtoukoff, Gustavo Abade, and Marta Waclawczyk

Turbulent clouds are challenging to model and simulate due to uncertainties in processes occurring at unresolved subgrid scales (SGS). These processes include the transport of cloud particles, supersaturation fluctuations, turbulent mixing, and the resulting stochastic droplet activation and growth by condensation. In this work, we use a stochastic Filtered Density Function (FDF) particle method to describe both temperature and the amount of water vapor distribution at unresolved scales in rising thermals. The particle method is solved in combination with an Eulerian scheme (here we use the CM1 model, Bryan and Fritsch 2002). In this hybrid implementation, the particle method plays two important roles: (a) it models the SGS variability of the environment where cloud particles would activate and grow by condensation, and (b) it provides the SGS unresolved fluxes that close the filtered equations underlying the Eulerian scheme, devoid of any approximation. The hybrid algorithm is tested in 2D and 3D simulations of dry and moist thermals against schemes that use standard subgrid models (such as Deardorff’s parametrization of the eddy viscosity). Our simulations suggest that FDF models may expand the ability of Large Eddy Simulations to represent cloud entrainment and associated microphysical details at the edge of cumulus clouds. Furthermore, the intrinsic character of this particle-based method is also useful to study generation and reorientation of vorticity around and within thunderstorm clouds.

How to cite: Jewtoukoff, V., Abade, G., and Waclawczyk, M.: Simulations of Rising Thermals Using a Lagrangian Stochastic Closure for the Subgrid Turbulent Fluxes, 11th European Conference on Severe Storms, Bucharest, Romania, 8–12 May 2023, ECSS2023-90, https://doi.org/10.5194/ecss2023-90, 2023.

17:45–18:00
|
ECSS2023-28
Mateusz Taszarek, Bartosz Czernecki, and Piotr Szuster

ThundeR is a freeware R language package for sounding and hodograph visualization, and rapid computation of convective parameters commonly used in the research and operational prediction of severe convective storms. Core algorithm is based on C++ code seamlessly integrated into the R language within the RCPP library. This solution allows to compute a large number of thermodynamic and kinematic parameters within hundredths of a second per atmospheric profile. Such performance enables to process large numerical datasets such as reanalyses or weather prediction models for the research and operational purposes. ThundeR package has been developed since 2017 and is constantly updated with new features and parameters following requests from the community and the most up to date research findings on severe storm environments. The most recent version of the package (v1.1) allows users to calculate 201 parameters, manually specify mixing and altitude of a convective parcel, input a manual storm motion vector, and control plotting of CAPE, CIN, DCAPE and SRH polygons. An online tool available at www.rawinsonde.com allows users to use thundeR package in visualizing rawinsonde measurements and historical atmospheric profiles from ERA5 reanalysis since 1950. As of January 2023, thundeR package has been used in more than 20 peer-reviewed studies.

How to cite: Taszarek, M., Czernecki, B., and Szuster, P.: thundeR - a rawinsonde package for processing convective parameters and visualizing atmospheric profiles, 11th European Conference on Severe Storms, Bucharest, Romania, 8–12 May 2023, ECSS2023-28, https://doi.org/10.5194/ecss2023-28, 2023.

Posters: Thu, 11 May, 14:30–16:00 | Exhibition area

Display time: Wed, 10 May, 09:00–Thu, 11 May, 18:30
P4
|
ECSS2023-7
Flavio T. Couto, Jean-Baptiste Filippi, Roberta Baggio, and Rui Salgado

The development of PyroCumulonimbus clouds during mega fire events has high impact in the evolution of the fire fronts and their development is frequently associated with strong convective processes due the heat and moisture released by the combustion. In 2017, Portugal was affected by several episodes of extreme wildfires with such a cloud system. The “Pedrógão Grande” mega fire caused more than 60 fatalities and burned a total area of almost 29,000 ha. In general, the atmospheric models do not account for fire-atmosphere interactions. Aiming to investigate the pyro-convective activity during the Pedrógão Grande mega fire, a numerical simulation was run with the Meso-NH model coupled to the ForeFire model. The Meso-NH model was configured into three nested domains. The horizontal resolution is 2000 m for the outer domain (600 km × 600 km). The inner computational grids have grid increments of respectively 400 m (120 km × 120 km) and 80 m (24 km × 24 km) for the innermost model. Initial and lateral boundary conditions for the outer domain are provided by ECMWF analysis, with updates every 6 h. The simulation with the coarsest resolution began on 17th June 2017 at 0600 UTC, with a progressive downscaling up to the finest resolution beginning at 1300 UTC. The vertical resolution is the same for all the nested domains, with 50 levels and a first level above the ground at 30 m height. The study used the reference fire propagation deduced from the official investigation (forced fire) and the emission of heat and vapour into the atmosphere was made using the ForeFire model. The results highlight the importance of the use of cloud resolving models configured with very-high spatial and temporal resolutions (80m, 10s) for representing the development of phenomena associated to pyro-convective activity, namely those occurring in the micro-scale from the cloud microphysics processes, like very-localised microbursts. This study was funded by national funds through FCT-Foundation for Science and Technology, I.P. under the PyroC.pt project (Ref. PCIF/MPG/0175/2019).

How to cite: Couto, F. T., Filippi, J.-B., Baggio, R., and Salgado, R.: Modelling fire-generated thunderstorms: Pedrógão Grande case study, 11th European Conference on Severe Storms, Bucharest, Romania, 8–12 May 2023, ECSS2023-7, https://doi.org/10.5194/ecss2023-7, 2023.

P5
|
ECSS2023-8
Vojtech Bliznak and Petr Zacharov

Atmospheric reanalysis are powerful tool to obtain information about the past state of the atmosphere by combination of the observations and numerical weather prediction (NWP) modelling through various data assimilation schemes. Major developments in NWP modelling have recently enabled to increase model resolution by nesting the regional reanalysis into the global reanalysis using a limited-area model. However, the relatively small area limits the use of those reanalyses for further meteorological, climatological and/or hydrological applications. A new atmospheric reanalysis ALADIN/PERUN was recently calculated to simulate various meteorological variables at high horizontal (2.3 km) and temporal (1 h) resolution over most of Europe in the 1989–2020 time span. Due to the high-resolution of the reanalysed data, precipitation fields can well capture the small-scale processes and thus well reproduce heavy convective precipitation. The contribution will evaluate this ability using adjusted radar-derived precipitation estimates as ground truth in the warm parts of the years 2002-2019 over the Czech Republic. Heavy convective precipitation will be selected from the dataset using defined criteria and their location and total amount will be evaluated for two different NWP model runs. The first one (A) includes the full assimilation of the observed data every 6 hours using the 4D-VAR assimilation scheme. The second one (B) uses only the boundary conditions from the ERA-5 global reanalysis and the prognostic data are not modified in any way. Comparing the two runs will provide us with information about the level of physical description in the NWP model as well as the effect of assimilation on the resulting precipitation fields.

How to cite: Bliznak, V. and Zacharov, P.: Ability of the new high-resolution atmospheric reanalysis ALADIN/PERUN to simulate heavy convective precipitation, 11th European Conference on Severe Storms, Bucharest, Romania, 8–12 May 2023, ECSS2023-8, https://doi.org/10.5194/ecss2023-8, 2023.

P6
|
ECSS2023-31
Natalia Pilguj, Mateusz Taszarek, Maciej Kryza, and Harold Brooks

Around 2–3 violent tornadoes hit Europe per decade. Although they are less frequent compared to the Great Plains of the United States, Europe has a population density higher by the order of the magnitude, thus extreme convective events typically produce large societal impacts and considerable damage to infrastructure. As operational prediction of such events is of high priority to emergency managers and forecasters, in this study we reconstruct environmental conditions for selected 12 violent tornadoes from the period between 1950 and 2021. Because numerical modeling of severe weather events requires convection-permitting resolutions, we used the Weather Research and Forecasting model for dynamical downscaling of ERA5 reanalysis. Results show that out of 12 high-resolution simulations (3 km), in 10 cases the modelled updraft helicity (UH) tracks cross the neighborhood of tornado reports, however, in 2 cases the timing of convection was more than 3 hours away from the time of the reports. Analysis of vertical profiles show that all simulations represented favorable tornado environments with cyclonically curved hodographs, comparable to significant tornadoes on the Great Plains of the United States. The combination of UH tracks with convective environments offers promising results for operational forecasting in Europe, and should be explored in future studies.

How to cite: Pilguj, N., Taszarek, M., Kryza, M., and Brooks, H.: Reconstruction of violent tornado environments in Europe: High-resolution dynamical downscaling of ERA5, 11th European Conference on Severe Storms, Bucharest, Romania, 8–12 May 2023, ECSS2023-31, https://doi.org/10.5194/ecss2023-31, 2023.

P7
|
ECSS2023-105
Robert Kvak, Petr Zacharov, Luboslav Okon, Vojtech Bliznak, Ladislav Meri, and Marek Kaspar

For decades, the role of mountains in influencing both a supercell environment during initiation stage and later storm dynamics was poorly understood. Since the consequences of terrain-induced effects on supercells seem to play a crucial part in the distribution of convective precipitation at least in mountainous regions, several topic-focused papers have already been published particularly in the last fifteen years. Our contribution to the issue is based on the examination of two-dimensional radar data of 62 supercells that initiated over the western part of the Carpathian Mountains in Central Europe in the years 2015–2019. To investigate the basic spatiotemporal characteristics and estimated precipitation intensity of storms, mainly Doppler and pseudo-CAPPI products from the Slovak weather radar network were employed in the study. Potential impact of terrain asymmetry on supercells precipitation intensity through their lifetime was evaluated by the relationship between morphometric variables of the terrain and radar reflectivity factor within the proximity (5, 10, 15-km radius) of detected mesocyclones in 5-min radar scans. Subsequently, atmospheric conditions along with convective parameters during the events were simulated in the COSMO model. Initial environments were compared with actual radar reflectivity properties (e.g., area of reflecting radar pixels >50 dBZ) of storms and their propagation qualities such as translation speed, duration, and track length. Although no general connection of radar reflectivity to underlying terrain even with a time shift up to 30 min was found, increasing reflectivity of some supercells showed preferences to more unstable environments which clearly were related to orographic modification. Besides, the poster highlights several other intriguing findings of supercells behavior over complex terrain.

How to cite: Kvak, R., Zacharov, P., Okon, L., Bliznak, V., Meri, L., and Kaspar, M.: Impact of mountains on environments and estimated precipitation intensity of supercells, 11th European Conference on Severe Storms, Bucharest, Romania, 8–12 May 2023, ECSS2023-105, https://doi.org/10.5194/ecss2023-105, 2023.

P8
|
ECSS2023-138
Sven Ulbrich, Philipp Zschenderlein, Martin Rempel, Alberto de Lozar, Tomas Pucik, and Ulrich Blahak

The German Meteorological Service (DWD) aims to enhance its ability to provide warnings on shorter notice (e.g. of extreme convective events) and therefore implements a rapid update cycle (RUC). This system produces more frequent and faster available numerical ensemble forecasts compared to the standard ICON-D2 short range numerical weather prediction (SRNWP) forecast procedure. The RUC’s focus is on lead times of up to 12 hours.

The RUC’s capability in predicting extreme convective events is tested by applying its setup to re-forecast the violent tornado case near Hodonin in the Czech Republic in June 2021. As this area is outside of the ICON-D2 domain, the model domain was shifted. While several German radar stations could not contribute to the new domain, 3D radar reflectivities and radial winds from several Czech and Slovakian radars were used for this case study. The assimilation cycle is started from the ICON-EU state the day before the event and hourly ensemble forecasts are started before the actual event. In addition, we included a nested domain, which uses a resolution of ~1km to evaluate the impact on the atmospheric parameters and dynamics.

We present results of the simulation of the atmospheric condition as well as tracks and dynamics of the developing convective systems with respect to different lead times and horizontal resolutions. As expected, there is an increase in forecast accuracy of the convective event with decreasing lead time. With respect to the different horizontal resolution, amplitudes of wind speed and updraft helicity are larger for the 1-km domain. Further, the rotation signature of the simulated tornadic supercell (e.g. radial wind) is comparable to the actually observed one.

How to cite: Ulbrich, S., Zschenderlein, P., Rempel, M., de Lozar, A., Pucik, T., and Blahak, U.: Tornado in Czech Republic – An ICON-RUC case study, 11th European Conference on Severe Storms, Bucharest, Romania, 8–12 May 2023, ECSS2023-138, https://doi.org/10.5194/ecss2023-138, 2023.