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.

Aerosols and organic source apportionment were characterized using data collected during two measurement campaigns. These campaigns were conducted during the summer and winter seasons at Măgurele, a site located southwest of Bucharest, the capital of Romania and one of the largest cities in southeastern Europe (raking seven in Europe based on population). The summer campaign was conducted between 7 June–18 July 2012, and the winter campaign from 14 January–6 February 2013. Approximately 50% of the organic fraction contribution to the total submicron particulate matter sampled by aerosol mass spectrometer was evidenced during both seasons. Submicronic organic aerosol sources were quantified using the positive matrix factorization approach. For warm (summer) and cold (winter) seasons, more than 50% from total organics was represented by oxidized factors. For the summer season, separate analyses were conducted on data influenced by urban and non-urban sources. The influence of pollution from Bucharest on the measurement site was observed in aerosol concentration and composition. The primary organic aerosols have different contribution percentage during summer, depending on their main origin. The influence of Bucharest, during summer, included cooking contribution of 13%. The periods with more regional influence were characterized by lower contribution from traffic and biomass burning in a total proportion of 28%. In winter, the influence of local non-traffic sources was dominant. For more than 99% of the measurements, the biomass burning indicator, $f60$ , exceeded the background value, with residential heating being an important source in this area. Fossil fuel contribution was confirmed for one week during the winter campaign, when ${}^{14}$ C analysis of total and elemental carbon revealed the presence of 17% fossil contributions to total carbon. Mass spectrometry, ${}^{14}$ C and absorption data suggest biomass burning as the predominant primary source of organic aerosols for the winter season.

An intensive winter campaign was organized for measuring the surface air pollutants in southeastern Europe. For a three months period, the gas concentrations of NO ${}_x$ , SO ${}_2$ , CO, O ${}_3$ , and CH ${}_4$ as well as meteorological parameters were simultaneously sampled to evaluate the variations and characteristic reactions between the gases during winter at the measuring site. The photochemical production of the ozone was observed through the diurnal variation of ozone and the solar radiation, the maximum concentration for ozone being reached one hour after the maximum value for solar radiation. A non-parametric wind regression method was used to highlight the sources of the air pollutants. The long-range transport of SO ${}_2$ and two hotspots for CO from traffic and from residential heating emissions were emphasized. The traffic hotspot situated north of the measuring site, close to the city ring road, is also a hotspot for NO ${}_x$ . The air quality during the cold season was evaluated by comparing the measured gas concentration with the European limits. During the measuring period, the values for NO ${}_2$ , CO, and SO ${}_2$ concentration were at least two times lower than the European Union pollution limits. Only twice during the study period was the concentration of O ${}_3$ higher than the established limits.

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.