Abstract. In September 2015 one of the severest and most unusual dust events on record occurred in the Eastern Mediterranean. Surprisingly, operational dust transport models were unable to forecast the event. This study details the reasons for this failure and presents simulations of the event at convection-permitting resolution using the modelling system ICON-ART. The results allow for an in-depth analysis of the influence of the synoptic situation, the complex interaction of multiple driving atmospheric systems and the mineral dust radiative effect on the dust event. A comparison of the results with observations reveals the quality of the simulation results with respect to structure and timing of the dust transport. The forecast of the dust event is improved decisively. The event is triggered by the unusually early occurrence of an active Red Sea trough situation with an easterly axis over Mesopotamia. The connected sustained organized mesoscale convection produces multiple cold-pool outflows responsible for intense dust emissions. Complexity is added by the interaction with an intense heat low, the inland-penetrating Eastern Mediterranean sea breeze and the widespread occurrence of supercritical flow conditions and subsequent hydraulic jumps in the vicinity of the Dead Sea Rift Valley. The newly implemented mineral dust radiation interaction leads to systematically more intense and faster propagating coldpool outflows.

No abstract

Shallow convection is a subgrid process in cloud-resolving models for which their grid box is larger than the size of small cumulus clouds (Cu). At the same time such Cu substantially affect radiation properties and thermodynamic parameters of the low atmosphere. The main microphysical parameters used for calculation of radiative properties of Cu in cloud-resolving models are liquid water content (LWC), effective droplet radius, and cloud fraction (CF). In this study, these parameters of fields of small, warm Cu are calculated using large-eddy simulations (LESs) performed using the System for Atmospheric Modeling (SAM) with spectral bin microphysics. Despite the complexity of microphysical processes, several fundamental properties of Cu were found. First, despite the high variability of LWC and droplet concentration within clouds and between different clouds, the volume mean and effective radii per specific level vary only slightly. Second, the values of effective radius are close to those forming during adiabatic ascent of air parcels from cloud base. These findings allow for characterization of a cloud field by specific vertical profiles of effective radius and of mean liquid water content, which can be calculated using the theoretical profile of adiabatic liquid water content and the droplet concentration at cloud base. Using the results of these LESs, a simple parameterization of cloud-field-averaged vertical profiles of effective radius and of liquid water content is proposed for different aerosol and thermodynamic conditions. These profiles can be used for calculation of radiation properties of Cu fields in large-scale models. The role of adiabatic processes in the formation of microstructure of Cu is discussed.

Abstract. The significance of planetary boundary layer (PBL) height detection is apparent in various fields, especially in air pollution dispersion assessments. Numerical weather models produce a high spatial and temporal resolution of PBL heights; however, their performance requires validation. This necessity is addressed here by an array of eight ceilometers; a radiosonde; and two models – the Integrated Forecast System (IFS) global model and COnsortium for Small-scale MOdeling (COSMO) regional model. The ceilometers were analyzed with the wavelet covariance transform method, and the radiosonde and models with the parcel method and the bulk Richardson method. Good agreement for PBL height was found between the ceilometer and the adjacent Bet Dagan radiosonde (33 m a.s.l.) at 11:00 UTC launching time (N=91 d, ME =4 m, RMSE =143 m, R=0.83). The models' estimations were then compared to the ceilometers' results in an additional five diverse regions where only ceilometers operate. A correction tool was established based on the altitude (h) and distance from shoreline (d) of eight ceilometer sites in various climate regions, from the shoreline of Tel Aviv (h=5 m a.s.l., d=0.05 km) to eastern elevated Jerusalem (h=830 m a.s.l., d=53 km) and southern arid Hazerim (h=200 m a.s.l., d=44 km). The tool examined the COSMO PBL height approximations based on the parcel method. Results from a 14 August 2015 case study, between 09:00 and 14:00 UTC, showed the tool decreased the PBL height at the shoreline and in the inner strip of Israel by ∼100 m and increased the elevated sites of Jerusalem and Hazerim up to ∼400 m, and ∼600 m, respectively. Cross-validation revealed good results without Bet Dagan. However, without measurements from Jerusalem, the tool underestimated Jerusalem's PBL height by up to ∼600 m.

Abstract. Flash floods are among the most devastating and lethal natural hazards. In 2018, three flash-flood episodes resulted in 46 casualties in the deserts of Israel and Jordan alone. This paper presents the hydrometeorological analysis and forecasting of a substantial storm (25–27 April 2018) that hit an arid desert basin (Zin, ∼1400 km2, southern Israel) claiming 12 human lives. This paper aims to (a) spatially assess the severity of the storm, (b) quantify the timescale of the hydrological response, and (c) evaluate the available operational precipitation forecasting. Return periods of the storm's maximal rain intensities were derived locally at 1 km2 resolution using weather radar data and a novel statistical methodology. A high-resolution grid-based hydrological model was used to study the intra-basin flash-flood magnitudes which were consistent with direct information from witnesses. The model was further used to examine the hydrological response to different forecast scenarios. A small portion of the basin (1 %–20 %) experienced extreme precipitation intensities (75- to 100-year return period), resulting in a local hydrological response of a high magnitude (10- to 50-year return period). Hillslope runoff, initiated minutes after the intense rainfall occurred, reached the streams and resulted in peak discharge within tens of minutes. Available deterministic operational precipitation forecasts poorly predicted the hydrological response in the studied basins (tens to hundreds of square kilometers) mostly due to location inaccuracy. There was no gain from assimilating radar estimates in the numerical weather prediction model. Therefore, we suggest using deterministic forecasts with caution as it might lead to fatal decision making. To cope with such errors, a novel cost-effective methodology is applied by spatially shifting the forecasted precipitation fields. In this way, flash-flood occurrences were captured in most of the subbasins, resulting in few false alarms.