The bioclimatology of thermal stress over Europe between 1979 and 2019 was analysed using the Universal Thermal Climate Index (UTCI) derived from ERA5‐HEAT reanalysis. The bioclimatology of different European regions was assessed using Köppen–Geiger climate classification. The annual number of hours with heat stress (UTCI > 32°C) increased significantly during the study period for all the analysed Köppen–Geiger climate subclasses, showing also a clear increase towards southern Europe. The highest percentage of hours (20% of all hours) with cold stress (UTCI < −13°C) occur over northern Europe. A significant increasing trend (>0.05 hr·year⁻¹) in the number of hours with heat stress was observed for 23 out of 32 analysed European cities representative for the Köppen–Geiger climate subclasses. For these cities not only the number of hours with heat stress has increased but also the heat stress is more persistent, while the number of cases and the persistence of the periods with cold stress have decreased over the last four decades. The UTCI values showed a statistically significant increase between 0.6 and 3.2°C for all the analysed cities over the study period reflecting the rising of global mean temperatures.

On the morning of 23 March 2018, an unusual phenomenon was observed over Romania where the southeastern part of the country was covered in a fresh-layer of orange snow. The event was extensively reported in mass-media and social-media and raised questions about the origin and the possible impact of the orange snow. Even if this type of events, intrusions of Saharan dust, have been reported before in Romania, and in Europe in general, their occurrence during negative temperature conditions is very rare. Saharan dust intrusion occurs over Europe mainly during spring and, in general, is not accompanied by snow at low altitudes. In this article, for the first time, the synoptic-scale conditions leading to the Saharan dust intrusion over Romania and the chemical and physical properties of the deposited dust particles in a snow layer were analyzed. The Saharan dust event affected a permanent atmospheric measurement research infrastructure located southwest of Bucharest, the capital city of Romania. In-situ and remote sensing measurements conducted at this research infrastructure allowed the identification of the dust source as the north Sahara. The source was confirmed by the elemental ratios of the main components (e.g., Al, Ca, Mg, Fe, K). For example, the (Ca+Mg)/Fe ratio of 1.39 was characteristic for the north Sahara. The dust morphology and the minerals were analyzed by scanning electron microscopy with energy disperse X-ray spectrometry (SEM/EDX). The size distribution of the particle geometric diameter showed that they are centred on 1 $\mathsf{\mu }$ m, but larger particles up to 40 $\mathsf{\mu }$ m are also present. To visualize the minerals, an approach was developed which emphasized the presence of the calcite, quartz or clay minerals. The optical parameters of dust were measured by re-suspending the particles. Values of the optical parameters (i.e., asymmetry parameter at 550 nm was 0.604, single scattering albedo was 0.84–0.89) were similar to those measured for Saharan dust intrusions over the Iberian Peninsula. Also, the non-refractory particles found in the dust-contaminated snow layer were analyzed, indicating the presence of HULIS-like compounds, most probably advected from the Mediterranean sea.

In this study, AERONET (Aerosol Robotic Network) and EARLINET (European Aerosol Research Lidar Network) data from 17 collocated lidar and sun photometer stations were used to characterize the optical properties of aerosol and their types for the 2008–2018 period in various regions of Europe. The analysis was done on six cluster domains defined using circulation types around each station and their common circulation features. As concluded from the lidar photometer measurements, the typical aerosol particles observed during 2008–2018 over Europe were medium-sized, medium absorbing particles with low spectral dependence. The highest mean values for the lidar ratio at 532 nm were recorded over Northeastern Europe and were associated with Smoke particles, while the lowest mean values for the Angstrom exponent were identified over the Southwest cluster and were associated with Dust and Marine particles. Smoke (37%) and Continental (25%) aerosol types were the predominant aerosol types in Europe, followed by Continental Polluted (17%), Dust (10%), and Marine/Cloud (10%) types. The seasonal variability was insignificant at the continental scale, showing a small increase in the percentage of Smoke during spring and a small increase of Dust during autumn. The aerosol optical depth (AOD) slightly decreased with time, while the Angstrom exponent oscillated between “hot and smoky” years (2011–2015) on the one hand and “dusty” years (2008–2010) and “wet” years (2017–2018) on the other hand. The high variability from year to year showed that aerosol transport in the troposphere became more and more important in the overall balance of the columnar aerosol load.

The changes in the characteristics of heatwaves over Romania have been analyzed using the excess heat factor calculated for two climate change scenarios (RCP4.5 and RCP8.5) from the EURO-CORDEX project. The changes were evaluated for the near future (2021–2050) using the historical period (1971–2000) as reference. The frequency of occurrence and the duration of heatwaves is projected to increase for both climate scenarios in particular over southern Romania. In this region, the percentage of change in the near future for the number of heatwaves is between 50 and 60% for the RCP4.5 scenario and 60–80% for the RCP8.5 scenario. Also for the same region, the duration of heatwaves will increase by 30–50% for the RCP4.5 scenario and 60–80% for the RCP8.5 scenario. These results indicate that the human exposure to heatwaves will increase in Romania in the near future. To increase awareness on heatwaves and their impact, we propose a series of immediate actions that include (1) improving the communication of the impact of heatwaves, (2) identification of the regions where the population is more vulnerable to heatwaves, and (3) better understanding of the mortality and morbidity associate with heatwaves in Romania.

A high-impact freezing rain event affected parts of southeastern Romania on 24–26 January 2019. The freezing rain caused extensive damages in Bucharest, the capital city of Romania. The meteorological analysis highlighted the presence of a particular synoptic pattern involving a high-pressure system advecting cold air mass at low levels, while at mid-levels a warm and humid intrusion was associated with a low-pressure system of Mediterranean origin. At Bucharest, the vertical profiles from ERA5 and radiosondes emphasized the presence of a thick warm layer between 1000–1400 m above the re-freezing layer close to the surface. A climatology of freezing rain events in Bucharest was built to understand the frequency and intensity of this phenomenon. On average, there were approximately 5 observations of freezing rain in Bucharest per year between 1980–2018. The number of consecutive freezing rain days was used as a proxy for the event severity. Moderate-duration events (2 consecutive days) represented 16 periods of all 59 non-overlapping freezing rain periods in Bucharest and long-duration events (3 consecutive days) represented 3 periods. The monthly distribution showed that freezing rain occurs more frequently between December–February with a maximum in December. The moderate and long-duration freezing rain events were associated with two main sub-synoptic patterns related to the Carpathians lee cyclogenesis.

The goal of this article is to investigate how the Romanian air space would be affected by a potential large eruption of Vesuvius volcano. In order to characterize the potential impact of such an event on the Romanian air space, we simulated the ash plume transport from a hypothetical explosive eruption with a significant magnitude, that is, Volcanic Explosivity Index (VEI 4). The volcanic ash cloud transport towards Romania was described by using the mass distribution from the HYbrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) for Volcanic Ash numerical model. To assess the dominant air circulation types, an adaptation of the WetterLagenKlassification (WLK) classification scheme was used. The results emphasized that a long‐wave trough and southwesterly flow in upper levels would favour the direct transport of volcanic ash cloud towards the Romanian air space. Non‐negative matrix factorization method has been used to identify the atmospheric layer with the maximum values of mass. The simulation results reveal that the maximum values of averaged mass were obtained at an altitude around 7 km.