Next Article in Journal
Comparison of Air Pollutants during the Two COVID-19 Lockdown Periods in Winter 2019 and Spring 2022 in Shanghai, China
Next Article in Special Issue
Dust Transport from North Africa to the Middle East: Synoptic Patterns and Numerical Forecast
Previous Article in Journal
The Verification and Fusion Analysis of Passive Microwave Soil Moisture Products in the Three Northeastern Provinces of China
Previous Article in Special Issue
Expected Changes in Heating and Cooling Degree Days over Greece in the near Future Based on Climate Scenarios Projections
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Synoptic Analysis and Subseasonal Predictability of an Early Heatwave in the Eastern Mediterranean

by
Dimitris Mitropoulos
1,2,*,
Ioannis Pytharoulis
1,2,*,
Prodromos Zanis
1,2 and
Christina Anagnostopoulou
1
1
Department of Meteorology and Climatology, School of Geology, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
2
Center for Interdisciplinary Research and Innovation (CIRI-AUTH), Balkan Center, Buildings A & B, 10th km Thessaloniki-Thermi Rd, P.O. Box 8318, 57001 Thessaloniki, Greece
*
Authors to whom correspondence should be addressed.
Atmosphere 2024, 15(4), 442; https://doi.org/10.3390/atmos15040442
Submission received: 30 January 2024 / Revised: 21 March 2024 / Accepted: 26 March 2024 / Published: 2 April 2024
(This article belongs to the Special Issue Numerical Weather Prediction Models and Ensemble Prediction Systems)

Abstract

:
Greece and the surrounding areas experienced an early warm spell with characteristics of a typical summer Mediterranean heatwave in mid-May 2020. The maximum 2 m temperature at Kalamata (southern Greece) reached 40 °C on 16 May and at Aydin (Turkey), it was 42.6 °C on 17 May. There was a 10-standard deviation positive temperature anomaly (relative to the 1975–2005 climatology) at 850 hPa, with a southwesterly flow and warm advection over Greece and western Turkey from 11 to 20 May. At 500 hPa, a ridge was located over the Eastern Mediterranean, resulting in subsidence. The aims of this study were (a) to investigate the prevailing synoptic conditions during this event in order to document its occurrence and (b) to assess whether this out-of-season heatwave was predictable on subseasonal timescales. The subseasonal predictability is not a well-researched scientific topic in the Eastern Mediterranean Sea. The ensemble global forecasts from six international meteorological centres (European Centre for Medium-Range Weather Forecasts—ECMWF, United Kingdom Met Office—UKMO, China Meteorological Administration—CMA, Korea Meteorological Administration—KMA, National Centers for Environmental Prediction—NCEP and Hydrometeorological Centre of Russia—HMCR) and limited area forecasts using the Weather Research and Forecasting model with the Advanced Research dynamic solver (WRF) forced by Climate Forecast System version 2 (CFSv.2; NCEP) forecasts were evaluated for lead times ranging from two to six weeks using statistical scores. WRF was integrated using two telescoping nests covering Europe, the Mediterranean basin and large part of the Atlantic Ocean, with a grid spacing of 25 km, and Greece–western Turkey at 5 km. The results showed that there were some accurate forecasts initiated two weeks before the event’s onset. There was no systematic benefit from the increase of the WRF model’s resolution from 25 km to 5 km for forecasting the 850 hPa temperature, but regarding the prediction of maximum air temperature near the surface, the high resolution (5 km) nest of WRF produced a marginally better performance than the coarser resolution domain (25 km).

1. Introduction

Intense and prolonged heat is the leading cause of weather-related deaths around the world [1,2,3]. Heatwaves form when very warm air masses are advected over a region and high pressures on the surface strengthen and remain over it for several days up to weeks. Due to global warming, a change in heatwaves in terms of severity, frequency and duration is expected [4,5], causing economic losses [6] and increased human mortality [7,8]. Recent research has shown that there is rapidly increasing likelihood of temperatures over 50 °C in the Mediterranean and the Middle East due to climate change [9]. The Mediterranean is recognized as one of regions that is the most sensitive and vulnerable to anthropogenic climate change on Earth [10] with climate models projecting a consistent warming over the 21st century for various emission scenarios [11,12,13]. Therefore, it is necessary to extend the subseasonal heatwave predictability of existing model frameworks beyond two weeks to provide the state with sufficient time for preparations and mitigation strategies to be implemented.
In recent decades, there have been advances in weather prediction with the help of the continuous advances in computing power, weather modelling and a better understanding of the atmospheric processes [14]. However, the predictability horizon seemed to reach an apparent plateau for extreme weather events such as heatwaves and cold spells [15]. On the other hand, in the last decade, great progress has been made in the subseasonal forecasting of heatwaves [16]. The discovery of certain sources of subseasonal predictability associated with atmospheric, oceanic and land processes [17] has led closer to a seamless prediction across timescales. One of these sources is the memory of the soil moisture which can last several weeks and influence the atmosphere through changes in evaporation and the surface energy budget, affecting subseasonal forecasts of air temperature [18]. Another source of subseasonal predictability is the sea surface temperature (SST) anomalies which can lead to changes in air–sea heat fluxes [19,20]. Also, further improvements in numerical weather prediction (NWP) such as a finer model resolution, improved data assimilation techniques and advancements in ensemble forecasting have increased the predictability horizon beyond the two weeks (lead time) [21,22]. Finally, the demand from users across the world for subseasonal forecasts and the recent international effort in forecasting on this timescale (e.g., ECMWF) database—(https://apps.ecmwf.int/datasets/data/s2s-realtime-instantaneous-accum-ecmf; accessed on 20 January 2024) have drawn more attention from scientists for further research on more accurate forecasts and eventually achieving a connection between weather forecasting and climate prediction, in other words, a seamless prediction.
Over the past few years, there has been increased research interest in the subseasonal prediction of heatwaves. These studies evaluated Subseasonal-to-Seasonal (S2S) models [16,17,23], dynamical downscaling [24] and the predictability and the role of different factors such as atmospheric blocking, soil moisture deficits and SST anomalies [25,26,27,28,29,30]. Also, S2S predictability research has been conducted using deep learning (artificial intelligence) techniques to improve the predictability of heatwaves and droughts [31,32,33]. Additionally, health-related research has been performed, investigating how human mortality and morbidity is affected by heatwaves and how well these events are predicted in S2S timescales [34].
However, most of the research on heatwave predictability in S2S timescales refers to global-scale studies and case studies in North America, East Asia, western Europe and Russia [27,28,29]. The S2S research on the predictability of heatwaves for the Eastern Mediterranean region is very limited and needs improvement, as this highly populated region is affected by high temperatures during the summer months. During the period that extends from late spring to early autumn, the countries surrounding the eastern Mediterranean Sea attract many tourists who plan their vacations even months before. It is noted that the Mediterranean countries received one third of the global tourism which was 330 million tourists in 2016 [35]. The Eastern Mediterranean includes a few low-income countries which highly rely their economy on tourism, and they are more vulnerable to climate change. Also, the Eastern Mediterranean countries are prone to wildfires during heatwaves with the notable cases of the recent wildfires in Attica (Greece) and Rhodes (Greece) in 2021 and 2023, respectively. Another example is the case of the 2007 extensive and lethal wildfires in Greece, which were thoroughly investigated for a connection between extreme temperatures and wildfire creation [36,37]. Water shortages are also prevalent during heatwaves. According to the latest Intergovernmental Panel on Climate Change (IPCC) report, the southeastern Mediterranean countries will double their irrigation requirements by 2080–2090 [38].
For these reasons, it is necessary to evaluate the predictability of weather models in the subseasonal timescale (minimum two weeks ahead) to allow the state to take action to mitigate the outcomes of a heatwave. This becomes imperative due to the fact that in recent decades, the heatwave period is expanding into spring and autumn [39,40,41]. Specifically for the Eastern Mediterranean, there is evidence [42] that many regions are experiencing earlier onsets of heatwaves. Additionally, the city of Athens (Greece) has already experienced some out-of-season heatwaves during June (e.g., in 2007) [42,43,44]. The abnormally early appearance (11–20 May 2020) of the heatwave of this study that mainly affected the whole Greek region with maximum temperatures reaching 40 °C in Kalamata (southern Greece) and 42.6 °C in Aydin (southwest Turkey) shows that a subseasonal prediction of such an event is crucial and truly beneficial as the people and the state are not expecting a heatwave to occur earlier than summer and it is more difficult for the human body to adapt to the heat in mid-summer [45,46].
The aims of this study were (1) to examine the prevailing synoptic conditions during the heatwave in the central–eastern Mediterranean and more specifically in Greece, in order to document its occurrence and (2) to assess the predictability (with lead times of two to six weeks) of (a) global S2S forecasts from various leading meteorological centres and (b) limited area S2S forecasts with a higher resolution than the global models, through dynamical downscaling with the WRF numerical model. This paper includes a brief review of the definitions of a heatwave, followed by the data and methodology (Section 2). The observational and synoptic analyses of the event are presented in Section 3.1 and Section 3.2, respectively. Section 3.3 outlines the assessment of the prediction performance of the S2S forecasts, while Section 4 provides the discussion and conclusions.

2. Materials and Methods

2.1. Data

To achieve the goals of this study, different datasets were required. Their main characteristics are summarized in Table 1. The WRF forecasts are derived from the analyses and forecasts of the CFSv.2 [47] (https://www.ncei.noaa.gov/products/weather-climate-models/climate-forecast-system, accessed on 23 June 2023) from the National Centers for Environmental Prediction (NCEP). This dataset was selected because it provides all the required meteorological fields (i.e., air temperature, dew point temperature, geopotential heights, wind speeds, etc.) across multiple vertical levels to the stratopause (1 hPa) and more importantly, the forecast window of this dataset covers an adequate time range of lateral boundary forecast data for the subseasonal forecasts. Every WRF run was initialized with the control analyses of CFSv.2 from the forecast cycles of 00 and 12 UTC and uses the corresponding six-hourly CFSv.2 forecasts (at 00, 06, 12 and 18 UTC) as boundary conditions. The grid spacing of the CFSv.2 gridded fields (used by WRF) is 1° × 1°.
Furthermore, for the predictability estimation of the global forecasts, the S2S Prediction Project Database [30] was used (https://www.ecmwf.int/en/research/projects/s2s and s2s.cma.cn, accessed on 25 January 2024) which provides S2S forecasts and reforecasts (hindcasts) from various meteorological centres around the world (ECMWF, UKMO, CMA, KMA, NCEP, HMCR, etc.). The real-time data are available twice per week, once per week or daily (depending on the providing centre) (Table 1) with a grid spacing of 1.5° × 1.5°. This resolution is not the native resolution of each model, but for consistency and comparability reasons, every meteorological centre’s model has been archived and disseminated with 1.5° grid spacing in the S2S database. The main reason that these gridded outputs were not selected to provide the initial and boundary conditions of the WRF simulations is that the S2S database does not provide the full spectrum of data fields required for the initialization of the WRF model. Additionally, the reforecast data [47] are provided with the same resolution as the real-time global forecasts, but with different production frequencies from UKMO and KMA. UKMO and KMA provide reforecast data four times per month in contrast to their daily real-time forecasts. In this study, the UKMO and KMA reforecasts that were examined are from 9, 17 and 25 April. Reforecasts prior to 9 April were not utilized because their lead times do not cover the entire duration of the heatwave.
Table A1 presents the weekly frequencies at which the six selected meteorological centres produced their ensemble forecasts. The selection of these six centres for this study was based on the length of their forecasts (lead time > 40 days) and a sufficiently high native resolution (Table 2).
The analysis was based on datasets of reanalysis data and observations from weather stations in major Greek cities (i.e., Athens, Thessaloniki and Heraklion). The reanalysis data were obtained from the ERA5 [48] Copernicus Climate Change Service database (ECMWF—https://cds.climate.copernicus.eu; accessed on 25 January 2024). ERA5’s grid spacing is 0.25° × 0.25° and its gridded data are provided for the examination of the temperature, geopotential heights, winds and other meteorological parameters at specific levels (mainly at 850 hPa and 500 hPa), as well as the near surface air temperature (at 2 m). In addition to ERA5, ERA5-Land [49] reanalysis data (0.1° × 0.1°) were used for analysis of soil moisture. ERA5-Land data were not used for the investigation of the predictability of maximum temperatures of the WRF simulations because the Copernicus database does not provide maximum temperatures, only hourly data. Finally, weather station observations were utilized for the investigation of the predictability and the production of heat indices. The observations were obtained from the National Observatory of Athens (NOA) and the National Climatic Data Center of the National Oceanic and Atmospheric Administration (NOAA) (https://www.ncei.noaa.gov, accessed on 25 January 2024). In total, 25 weather stations, from numerous urban areas and airports in Greece and the surrounding countries, were taken into consideration for the study.

2.2. Methodology

2.2.1. Definition of a Heatwave

Heatwaves are prolonged periods with extremely high temperatures that exert a considerable impact not only on human health and ecosystems, but also on infrastructure and social services. There is no global definition of a heatwave because local acclimatization and adaptation influence the impact of excessive heat. Even at the local level, there can be multiple heatwave definitions based on varying temperature levels or time periods. In the literature, there is a plethora of heatwave definitions which are related to numerous factors such as human morbidity [50,51], mortality [52,53,54] and even cattle milk production [55]. Heatwaves are commonly characterized by temperature and duration thresholds [56], in addition to humidity and diurnal temperature cycle characteristics that are related to human morbidity and mortality [57].
Τhere are some extensively used heatwave indices like the THI (Temperature–Humidity Index) [58,59] and EHF (Excess Heat Factor index) [60], which are based on simple, readily available meteorological variables such as the maximum, minimum and mean temperature and the dew point temperature. These variables are easily found in various datasets (e.g., in ERA5 and CFSv.2). THI and EHF have been successfully used in heatwave studies in Greece [61,62]. Additionally, there are other more sophisticated heat indices which focus on specific aspects like human discomfort, but with additional requirements on data. One of these heat indices is the Universal Thermal Climate Index (UTCI) [63].
In this study, the EHF index was mainly used (in the analysis of the event) because it only requires temperature values, it is applicable to different climate regimes and it is supported by the recent literature [60,62,64,65,66,67,68]. Additionally, the THI, UTCI and the 95th percentile of maximum 2 m temperatures (T2max95) were also utilized for the assessment of the heatwave’s onset, duration and intensity.

2.2.2. Excessive Heat Factor, Temperature–Humidity Index and Universal Thermal Climate Index

The EHF [60] is composed of two factors. The first factor (Excessive Heat Index of significance—EHIsig) is a measure of how warm a three-day period is with respect to a climatic threshold, such as the T2max95, at each location. In the literature, for the calculation of EHIsig, the minimum, the mean or the maximum 2 m air temperature is used. In this study, in line with [27], only the T2max was utilized as the combination of mean, minimum and maximum air temperature is optional. Additionally, the 90th or 95th percentile of the temperature was used in the literature. Thus, the T2max95 of May (reference period: 1975–2005) was selected to assess the characteristics of the heatwave more rigorously. If the first factor is above a chosen temperature threshold for three consecutive days, then this period gives a hint of the presence of a heatwave. The second factor (Excessive Heat Index of acclimatization—EHIaccl) is a measure of how warm the same three-day period is relative to the previous 30-day period. This concept suggests that individuals may adapt to the typical temperature variations in their local climate, taking into account changes across latitudes and seasons. However, they may not be prepared for a sudden rise in temperature exceeding that experienced in the preceding 30 days. The formulas for the components of EHF are
E H I s i g = T i + T i + 1 + T i + 2 3 T 2 m a x 95
E H I a c c l = T i + T i + 1 + T i + 2 3 T i 1 + + T i 30 30
The product of these two indices results in the EHF, with the specificity that the EHF must have the same sign as EHIsig. The EHF formula is
EHF = EHIsig × max(1, EHIaccl)
The unit for the EHF is degrees Celsius squared (°C2) and a three consecutive day period of EHF > 0 °C2 is indicative of the presence of a heatwave since day i (Equations (1) and (2)).
The Temperature–Humidity Index (THI) represents the stress to the human body from the combination of humidity and heat [69,70]. Its formula is
THI = 0.4 (T + Tw) +15
where T is the near-surface (2 m) temperature and Tw is the wet-bulb temperature in degrees Fahrenheit. The threshold used in this study is THI = 70 °F. Above this value, half of the population of a particular area feels uncomfortable. The selection of this threshold is conservative mainly because the heatwave (of this study) occurred earlier than the hot period of summer and secondly due to the relatively low resolution of ERA5. Also, the THI does not need a three consecutive days period above its threshold to indicate a heatwave [71].
The Universal Thermal Climate Index (UTCI) serves as an additional tool for identifying the onset, end and severity of the heatwave using a threshold of 32 °C. The UTCI incorporates a complex set of factors (i.e., radiant temperature and thermophysiological factors like net body heat loss) [72] that are not readily available, easily acquired or implemented in the WRF model. For that reason, UTCI was only calculated using the ERA5 dataset without any modification to implement it in the WRF forecasts. Above the UTCI threshold of 32 °C, people can experience varying degrees of discomfort, ranging from “strong heat stress” to “extreme heat stress” [63,73].

2.2.3. Limited-Area WRF Simulations

In order to assess the effect of dynamical downscaling on the prediction ability of the model, limited area numerical forecasts were performed using version 4.1.2 of Advanced Research WRF model [74]. The outer domain (d01) employed a horizontal grid spacing of 25 × 25 km and covered Europe, northern Africa, western Asia and a large part of the northern Atlantic to represent the synoptic scale flow. The inner nested domain (d02) covered Greece and western Turkey with a finer grid spacing of 5 × 5 km in order to investigate the effect of fine resolution forecasts on the predictability of the heatwave. The NCEP/CFS analyses and forecast fields of the control runs were utilized as initial and boundary conditions (for d01, updated every six hours) because they provided all the necessary input fields for the WRF runs. In total, 42 simulations were performed with initialization dates from 9 April 2020 00 UTC to 29 April 12 UTC (two initializations starting at 00 and 12 UTC daily). The characteristics of the WRF forecasts are summarized in Table 3 and the telescoping domains are depicted in Figure 1. All the simulations had the same forecast horizon (until 31 May 2020). This forecast horizon was selected to identify any temporal shifts in the predictions of the heatwave as the heatwave occurred between 11 and 20 May. The Moderate Resolution Imaging Spectroradiometer 21-category dataset from the International Geosphere Biosphere Program (MODIS-IGBP) was used for the specification of the land-use categories and the United States Geological Survey (USGS) dataset was used for the topography.

2.2.4. Global S2S Reforecasts

The global reforecasts from the S2S database were utilized in order to take into account the systematic biases of each model. The reforecasts were compared with the real-time forecasts of temperature at 850 hPa at S2S timescales when the forecast coincides with the reforecast initialization day. This is justified because some meteorological centres (ECMWF, CMA, HMCR and NCEP) produce reforecasts along with forecasts, but other centres (UKMO and KMA) produce reforecasts only four days per month (on the 1st, 9th, 17th and 25th) in addition to their daily production of forecasts. Anomaly Correlation Coefficients (ACCs) were calculated for each global model using its reforecasts [81]. Moreover, the reforecasts were compared with the ERA5 climatology data.

2.2.5. Predictability Statistics and Scores

To produce the prediction statistics and performance scores, it was necessary to define common initialization periods between the limited-area S2S forecasts of WRF and the global ones of the S2S database. The initialization dates were separated into three seven-day periods (weeks): (1) 9–15 April, (2) 16–22 April and (3) 23–29 April. These three weeks were selected because all the models (global and limited-area models) have a certain number of ensemble members each week (Table A1).
The prediction performance scores were calculated for the temperature at 850 hPa (T850) and the maximum daily temperature at 2 m above the ground (T2max). The scores for T2max were only used in the WRF simulations (25 × 25 km and 5 × 5 km resolution) as the S2S global forecasts showed no prediction ability at the S2S timescale at the surface since they are archived at 1.5° × 1.5° and it was not possible to obtain their native resolution output. Additionally, the validation period was from 8 to 23 May 2020 which includes the heatwave (that occurred from 11 to 20 May, as it will be shown in Section 3.1 and Section 3.2) and non-heatwave days. In this way, the contingency table includes hits, misses, false alarms and true negatives [82]. If the validation period was the exact length of the heatwave, then false alarms and true negatives would be missing.
Before producing the statistics and the scores of the global and WRF S2S forecasts, all data (including reanalysis data) had to be brought to the same resolution. For the predictability of the 850 hPa temperature, a grid spacing of 1.5° was selected because it was the lowest resolution. Therefore, the ERA5 reanalysis and WRF forecasts were regridded to 1.5° × 1.5° using the “conserve” interpolation method. However, for the T2max, only the WRF forecasts and ERA5 were compared, so the WRF forecasts (25 km and 5 km) were regridded to the grid spacing of ERA5 (0.25°) and its grid geometry.
In the study, firstly, boxplots of the weekly ensembles were examined, accompanied by quantile–quantile (Q-Q) plots which revealed scale differences and identified the extremes of two different temperature distributions (observations and forecasts) [83]. Furthermore, receiver operating characteristic (ROC) curves [84] accompanied by reliability diagrams [85] were evaluated in terms of probabilistic forecasts using contingency tables with a yes/no decision threshold for the 95th percentile of the climatological temperature (T850 and T2max) at each grid point. ROC curves measure the ability of a forecast to discriminate between two alternative outcomes (resolution). Also, ROC curves are conditioned by the observations, answering the question “what was the corresponding forecast of a heatwave that has occurred” [86]. On the other hand, the reliability diagram is conditioned by the forecasts and shows how well the predicted probabilities of a heatwave’s occurrence correspond to their observed frequencies [87]. In other words, the reliability diagram indicates the reliability of a probabilistic forecast and whether conditional biases are encountered.
Finally, another statistical score that utilizes the hit rate (HR—Equation (5)) and false alarm rate (FAR—Equation (6)) is the extremal dependence index (EDI—Equation (7)) [88]. It is used to assess the association between the probabilistic forecasts and the observations. It stands out from the other scores, for example, the HR, because it is a non-degenerative measure of the predictability of rare binary events as the yes/no threshold for the heatwave for this study is the 95th percentile of air temperature. A perfect forecast has an EDI = 1 and an inaccurate forecast in S2S timescales has an EDI = 0.
H R = H i t s H i t s + M i s s e s
F A R = F a l s e   A l a r m s T r u e   N e g a t i v e s + F a l s e   A l a r m s
E D I = log F a l s e   A l a r m s T r u e   N e g a t i v e s + F a l s e   A l a r m s log H i t s H i t s + M i s s e s log F a l s e   A l a r m s T r u e   N e g a t i v e s + F a l s e   A l a r m s + log H i t s H i t s + M i s s e s

3. Results

3.1. Observational Analysis

This section aims to document the occurrence and the duration of the heatwave via the analysis of station observations and reanalyses. According to the EHF (weather station data), the early heatwave started on different dates (from 10 to 14 May) in the three main cities of Greece (Athens, Thessaloniki and Heraklion), its end occurred between 19 and 21 May and the highest maximum air temperatures occurred mostly on 16 May in Greece and the surrounding countries (Figure 2 and Table 4). The chronological differences are mostly justified by the latitudinal differences, because Thessaloniki (Northern Greece) and Heraklion (Crete, Southern Greece) are five degrees of latitude apart and on 21 May, the mainland of Greece was covered by a thick cloud layer of stratocumulus and convective clouds. Additionally, the EHF was calculated using ERA5 and ERA5-Land reanalyses to show that all three datasets (weather observations, ERA5 and ERA5-Land) indicated the presence of a heatwave because the EHF is above 0 °C2 for at least three consecutive days. ERA5-Land and ERA5 agreed on the maximum EHF values in Thessaloniki and Heraklion.
The maximum THI and UTCI values occurred on 16 and 17 May, except for the THI in Thessaloniki, which was marginally higher on 19 May compared to the previous three days (16 to 18 May). Both the THI and UTCI values were marginally higher than their thresholds, which were 70 °F (mild discomfort) and 32 °C (strong heat stress), respectively. It should be noted that the THI and UTCI thresholds are not equal in terms of heat stress severity; however, the factors that contribute to each index are quite different. Additionally, UTCI is much more modern, well tested, and more sophisticated and complex than THI. Thus, UTCI’s threshold is more appropriate for evaluating whether a heatwave has taken place.
As mentioned previously, the peak intensity of the heatwave was found on different dates, but the most robust results were provided by the EHF from the weather stations as the reanalysis data do not perfectly describe the state of the atmosphere [48,49]. The peak intensity of the EHF occurred on the 16th and 17th of May. The majority of the maximum temperatures, which were measured by the weather stations of NOA and the Hellenic National Meteorological Service (HΝMS) across Greece and the surrounding countries, were observed on 16 May 2020 (Figure 2). The highest temperature (40 °C) in Greece was recorded in Kalamata city (KAL in Figure 2). The KAL station is owned by HNMS and it does not have a complete temperature time series to derive its reference climatology 95th percentile maximum temperature (nan in KAL). However, it was included in the map because it exhibited the highest maximum air temperature recorded in Greece during the event. This was a record high (absolute maximum) 2 m temperature for Kalamata for the month of May. The same temperature (40 °C) was recorded again the next day (17 May). In agreement with [42], mesoscale features seemed to have affected the maximum temperature at the airport of Kalamata (where the maximum 2 m temperature of 40 °C was observed in Greece on 16 May 2020). Hot and dry light winds from the land, together with the insolation, appeared to have contributed to the synoptic scale warm advection at the lower troposphere, in order to reach a 2 m temperature of 40 °C at 12:20 p.m. local time (09:20 UTC) on 16 May. However, the onset of the sea breeze and the proximity to the coast (at a distance of 4 km) led to an abrupt decrease in the 2 m temperature by 6 °C in the next 30 min. Despite the clear sky and the insolation of the early afternoon hours, the 2 m temperature dropped further to 30 °C until 16:50 local time (13:50 UTC) under the influence of the sea breeze. The highest temperature (42.6 °C) in the whole region of interest appeared in Aydin, Turkey. This temperature was 11.7 °C higher than its T2max95 and it was also the highest maximum temperature recorded for the month of May in Aydin. Additionally, all the stations shown in Figure 2 exceeded their T2max95 of May by 2 to 11 °C (reference period: 1975–2005).
According to the synoptic systems in the vicinity of Greece, at the beginning of the heatwave, there were stable conditions, with a sea level pressure between 1012 and 1017 hPa; there were no apparent fronts and the surface winds were weakly blowing from SSW and SW. At the end of the heatwave (20 May 2020), a cold front approached from the west, bringing in cold air masses and lowering the mean daily temperature across the Eastern Mediterranean. Finally, the end of the heatwave in Greece was indicated by the EHF and THI indices and T2max95 (station data and ERA5) for Athens and Thessaloniki.

3.2. Synoptic Analysis

Extreme temperatures were also prominent in the lower troposphere. In Figure 3 at 850 hPa, the reanalysis data indicated temperature increases of ten to eleven standard deviations or around 14 °C above the mean value for the Aegean Sea, eastern mainland of Greece and western coast of Turkey. These maximum anomalies were observed on 16 May. The reason for this severe warming is the strong horizontal thermal advection (maximum value: 3 °C/h—Western Greece and Albania) from the southwest during the event. Using the HYSPLIT air parcel trajectory tool [89], we found that the warm air masses originating in the Libyan region and the Gulf of Sidra moved anticyclonically north and northeast, affecting Southern Italy, Greece (Figure 4), and afterwards western Turkey. Figure 5 illustrates the maximum extent of this warm intrusion in the European region at the maximum intensity of the heatwave, with temperatures up to 25 °C at 850 hPa. Additionally, a broad mid-tropospheric ridge prevailed at the eastern Mediterranean Sea and a deep trough was located over the Iberian Peninsula and northwest Africa. The ridge was associated with subsidence leading to adiabatic warming through thermal compression of the air.
Furthermore, in Figure 6, a strong positive anomaly in the geopotential heights was observed above Greece and western Turkey. Specifically, above the Aegean Sea and southwestern Turkey, the geopotential height anomalies at 500 hPa (Figure 6a) and 850 hPa (Figure 6b) were, respectively, at least ten and seven standard deviations higher than the climatological value of May. On the other hand, there was a negative anomaly in the geopotential heights at the trough above Morocco, which was intense, but not as extreme as the positive anomaly. This dipole of anomalies amplifies the poleward flow and the advection of warm air masses from northern Africa to the Mediterranean Sea and southern Europe [91,92].
At the level of the jet stream (Figure 7a) at 300 hPa, there is a meridional (southwesterly) flow of the jet stream above Algeria and the central Mediterranean Sea which turns anticyclonically and becomes zonal above the Balkans, where the jet streak is located. According to [93], a pattern of this kind, where the jet stream is north-northwest of Greece and anticyclonically curved, causes strong subsidence (Figure 7b) of the air masses leading to adiabatic heating of the lower troposphere above Greece. This contribution further enhances the heating from the warm air advection in the lower troposphere, which seemed to be the main driver of the heatwave.
The two factors that potentially affected the intensity of the heatwave were the soil moisture in the upper layer surface (0–7 cm) anomaly and the SST anomaly. The soil moisture anomaly (Figure 8a) showed a water deficit in the upper soil layer, which mostly interacts with the atmosphere. The maximum observed soil moisture deficit (negative anomaly) was −0.1 m3 in northwestern Greece. A persistent low water content has the potential to contribute to heating through the prevalence of the sensible heat fluxes [94,95,96,97]. Also, the SST anomaly (Figure 8b) showed a general warming of the surface waters at around two standard deviations above the climatological value. Specifically, around the island of Crete, the difference was four to five standard deviations higher with a maximum of six standard deviations at 23.6 °C. The warm SSTs were persistent throughout the duration of the heatwave, contributing to further heating of the atmosphere. The reasons for this strong heating might be the prolonged low soil moisture accompanied together with the prevailing ridge at 500 hPa. Finally, a wedge of warm air coming from the Sahara Desert raised the temperature of the lower troposphere and the position of the jet stream northwest of Greece resulted in a downward atmospheric motion that led to the amplification of the warming due to thermal compression of the air. These atmospheric and land factors have previously been studied for heatwaves in relation to soil moisture content [94,98,99,100,101], blocking highs [102,103,104], atmospheric heat transfer [105], jet stream position [93] and SSTs [106].

3.3. S2S Prediction Ability

Initially, the boxplots (or whisker plots) of the temperature at 850 hPa using the full distribution of the forecasts (Figure 9) were generated for the region 34.5° Ν–42° N, 19.5° E–29.5° E, hereafter called the area of interest. The boxplots were based on the deterministic and the perturbed forecasts and they include the variation between the ensemble members. The ERA5 boxplot (leftmost boxplot) indicated a median (Q2) value of 17 °C, a 75th percentile (Q3) of 22.3 °C, a 25th percentile (Q1) of 11.1 °C, a maximum (Q4) of 27.1 °C and a minimum (Q0) of 0.7 °C in the period of 8–23 May 2020. For the initialization period of 9–15 April (hereafter referred to as week 1), the median of all the S2S forecasts was close to the Q1 of ERA5 and their Q3 was below the ERA5 median. On the other hand, the extreme values (Q0 and Q4) of the S2S forecasts (global and WRF) were closer to the extremes of ERA5 with the exception of the HMCR model. Similar results were obtained in the 16–22 April initialization period (week 2) where the CMA model demonstrated the best resemblance to the reanalysis. Also, the CMA ensemble forecast extremes were higher than the extremes of ERA5. The WRF forecasts showed smaller compatibility with ERA5 temperatures in week 2. In the initialization period of 23–29 April (week 3), most of the models presented higher temperatures than in the previous weeks, except for HMCR and CMA. The WRF forecasts showed improvements with a higher resemblance to ERA5 presented by WRF 5 km. Its median was higher and closer to the median of reanalysis than WRF 25 km.
The boxplot diagrams of the WRF deterministic forecasts of the daily maximum 2 m temperature (T2max; Figure 10) had smaller temperature ranges than ERA5, spanning from 9.9 to 29.8 °C (WRF 5 km from week 3), while the range of the reanalysis was from 11.2 to 33 °C. Figure 10 shows that the WRF forecasts could not predict the extreme high temperatures in week 1 and week 2. Although WRF 5 km predicted higher temperatures in Q4 than WRF 25 km, its median was lower than that of WRF 25 km. Week 2 presented the same behaviour as week 1, but in week 3, both WRF 25 km and 5 km showed an improvement. More specifically, the Q3 and Q4 of WRF 25 km increased. The same occurred with WRF 5 km, which exhibited a larger increase in Q4 than WRF 25 km. Quantitatively, the difference between the high extremes was Q4-reanalysis − Q4-WRF-5km = 3.2 °C, while Q4-reanalysis − Q4-WRF-25km = 5.4 °C in the WRF forecasts from week 3.
Looking into the quantile–quantile (Q-Q) distributions of the predicted temperatures using the full distribution of the global forecasts, and the deterministic WRF forecasts and the reanalysis at 850 hPa from week 3 (Figure 11a), it is obvious that the highest temperature quantiles were underestimated by all the models. Regarding the upper tail (>90%), the highest resemblance with the reanalysis was achieved by KMA and UKMO. Additionally, UKMO had the closest distribution fit with ERA5 as its mean squared error (MSE = 16.57) was the smallest of all the models. Generally, the S2S models and reanalysis agreed (i.e., they are closer to the diagonal) mostly at the quantiles between 15% and 50%. WRF 5 km had higher temperatures than 25 km across all the distribution percentages at 850 hPa. In Figure 11b, WRF 25 km and 5 km had similar maximum daily temperature distributions at 2 m and both nests presented similar MSE values. In the quantiles above 70%, WRF 5 km presented slightly higher temperatures and was closer to ERA5 than WRF 25 km. However, the maximum temperatures were underestimated in the 850 hPa Q-Q profile. For example, the 90% quantile was underestimated by WRF by 5.2 °C. In higher quantiles, the gap between WRF and ERA5 increased.
The uncentered Anomaly Correlation Coefficients (ACCs) of the global models were calculated according to each model’s climatology using reforecast data. The ACC compares the predicted anomalies against the observed anomalies (reanalysis) and offers a way to quantify the prediction ability of each forecasting model and assess its prediction accuracy. The ACC is presented separately for the unperturbed (deterministic) forecasts and the forecasts with perturbed initial conditions for various initialization dates (9–29 April 2020 00 UTC—first row) (Table 5). The results of the perturbed forecasts are the ACCs averaged over all the corresponding ensemble members. Note that the ACC of each individual perturbed member was calculated before averaging. The blank cells indicate days when both real-time forecasts and reforecasts had not been produced as most of the centres (except NCEP) do not produce daily real-time forecasts and reforecasts. Also, although UKMO and KMA produce daily forecasts, their reforecasts are produced only four times per month (on the 1st, 9th, 17th and 25th). Red cells indicate the initialization dates with an ACC lower than 0 for the evaluation period of 8–23 May 2020. An ACC near zero indicates that the forecast is unreliable, while negative values show that the observed and the predicted anomalies are negatively correlated. According to Table 5, the highest and sustained (without large values fluctuations above the ACC = 0.6 threshold [107]—Figure A1) ACCs above 0.6 were found in the perturbed forecasts of ECMWF and NCEP after 20 and 25 April, respectively. The perturbed forecasts of ECMWF reached ACC values between 0.69 and 0.78 after 20 April 2020. Furthermore, the UKMO forecasts for 25 April exhibited high ACC values of 0.84 and 0.59 for the unperturbed and the perturbed members, respectively. The unperturbed forecasts of KMA from 25 April exhibited a high ACC of 0.72. The perturbed forecasts of NCEP from 25 April appeared to have the highest credibility because the daily ACC values were generally stable without a substantial decrease as seen in the NCEP unperturbed forecasts. If there were more frequent forecasts and reforecasts from ECMWF, UKMO and KMA, it would be possible to check the stability of their ACC values. ECMWF, UKMO and NCEP also demonstrated the lowest temperature bias (<0.3 °C) between ERA5 climatology and each model’s climatology (reforecasts) (Figure A2, Figure A3, Figure A4, Figure A5, Figure A6 and Figure A7), while the rest of the models overestimated the temperature at 850 hPa by 1–3 °C. Only HMCR underestimated the temperature by 1 °C.
In Figure 12, the timeseries of the hit rate (HR) and false alarm rate (FAR) is presented for each initialization week using the WRF forecasts and the perturbed members of the global S2S forecasts. The 850 hPa temperature forecasts (of all the ensemble members at each grid point of the area of interest during the validation period) were compared to their 95th percentile in the ERA5 climatology of 1975–2005. It was obvious that the HR of all the forecasts, except for the ones from HMCR, increased from week 1 to week 3, indicating that the models detected the heating during the heatwave when they were initialized in week 3. The highest HR (0.71) was obtained by the KMA model when it was initialized in week 3. However, this increase was not always monotonic as the initialization period approached the actual heatwave. On the other hand, the FAR generally increased at almost the same rate as the HR. This means that although the models detected the high temperature signal, their forecasts exhibited large spatiotemporal errors.
In Figure 13a, a ROC diagram is presented for week 3 using the perturbed members of the global forecasts. The vertical axis represents the hit rate (HR) and the horizontal axis the false alarm rate (FAR) for the different probability thresholds (PTs) of the S2S forecasts. The PT is the percentage of forecast ensembles (for each model) which exceeded the climatological 95th percentile of the ERA5 temperature at 850 hPa on every day of the validation period at 00 UTC at each grid point of the area of interest (34.5° Ν–42° N, 19.5° E–29.5° E at a grid spacing of 1.5° × 1.5°). The ERA5 climatology was used instead of each model’s climatology because the temperature biases from the ERA5 climatology, as mentioned previously, are relatively small and there are no WRF reforecasts to verify its model climatology. Each point in the diagram of Figure 13a represents an approximately 17% probability increase starting from the first point to the top-right corner of the plot. For instance, if a model has 100 ensembles members (like ECMWF), a 17% probability means that 17 out of the 100 ensemble members surpassed the 95th percentile temperature threshold derived from the climatology (ERA5). The first (second) point represents a PT of approximately 17% (34%), etc. The 17% increment was chosen to ensure that even the model with the smallest ensemble size (CMA with six members per week) had at least one member for each PT. Also, in the ROC diagram, the Area Under the Curve (AUC) score of each model is presented. An AUC higher than 0.5 indicates a prediction ability higher than random guessing. Therefore, a poor prediction line was plotted (dashed black diagonal). Theoretically, a perfect forecast has an AUC equal to 1 [108].
In week 1 and week 2 many models demonstrated an S2S prediction ability because their AUC was above 0.5. The three highest AUC scores, in descending order, were achieved by WRF 5 km, KMA and WRF 25 km for week 1 and CMA, WRF 25 km and WRF 5 km for week 2 (Figure 13b). In week 3 (Figure 13a), the highest AUC scores were exhibited by ECMWF, KMA and WRF 25 km. The remarkably high AUC score of CMA in week 2 (Figure 13b) was related to three very important ensemble members, which constituted a large part (50%) of the influence of all members and raised the PT of the 50% HR to FAR ratio well above the diagonal. However, it is possible that these members’ effects may have occurred due to random chance because the AUC of CMA dropped from 0.84 (week 2) to 0.5 (i.e., poor prediction in week 3), while the rest of the models showed an improvement in their forecasts or demonstrated similar AUC scores as the initialization time approaches the heatwave. The investigation of random successful S2S forecasts was not included in the aims of this study, but it is an interesting topic for further research. Also, a model with a sufficient amount of ensemble members (above ten) minimizes the effects of random hits on the prediction ability [109,110,111]. The WRF S2S forecasts presented AUC scores above 0.5 for all three weeks. The two nests of WRF did not exhibit noticeable AUC differences. A substantial difference was found when comparing WRF forecasts with the global forecasts, especially for weeks 1 and 2. The only model that exhibited similar AUC values to WRF was the KMA model. In week 3 the AUCs of NCEP and ECMWF increased significantly by 0.15 and 0.22, respectively. Lastly, the WRF AUC scores for all the initialization periods were better than those of NCEP, which provided the initial and boundary conditions for the WRF forecasts.
The ROC diagram in Figure 14a shows the performance of the deterministic forecasts of WRF 25 km (blue curve) and 5 km (red curve) for T2max from week 3. WRF 5 km showed a slightly better prediction performance (AUC = 0.65) compared to WRF 25 km (AUC = 0.64). In weeks 1 and 2, the models exhibited poor prediction abilities, as their AUC scores were 0.49 to 0.51 for both grid spacings (Figure 14b).
The reliability diagrams (Figure 15) for each global (perturbed members) and limited-area model (WRF; deterministic) indicate the reliability of the S2S forecasts, that is, how well the predicted probabilities (x-axis) of the validation period of the heatwave correspond to their observed frequencies (y-axis) [87,112]. The correspondence is perfect when the plot is identical to the diagonal line (dashed orange line). The dashed red trendline shows the forecast probability range where the S2S forecast has the greatest sharpness (see histogram bars with frequencies above median). This trendline allows for a comparison with the perfect reliability diagonal. If the reliability curve exhibits a positive slope, it signifies that as the predicted probability of the event rises, the confirmed likelihood of observing the event also increases. Hence, the forecasts demonstrate a degree of reliability. If the trendline has a flatter appearance, then the model has a low resolution, for example, in the UKMO forecast for week 3 (Figure 15, first row, second column). Resolution pertains to the ability of a forecast to distinguish between different categories of outcomes that can occur. When the slope is below the diagonal, it suggests an imperfect reliability. Lastly, the histogram bins were derived from the ensemble size of each model.
ECMWF (Figure 15) showed a more robust reliability than most S2S models for all three weeks (only week 3 is shown). Its trendline lies above the diagonal, which is a sign that the forecast frequency underestimated the observed frequency of the 850 hPa temperature and exceeding its 95th percentile (underforecasting). ECMWF in week 3 displayed a good reliability with a 40 to 55% probability of a correct forecast because the sharpness (observed frequencies in histogram) was high and the diagram was closer and nearly parallel to the diagonal. Similarly, the forecasts of KMA and NCEP from week 3 (Figure 15) exhibited good reliability and resolution, with the trendline of the former model closer to the diagonal. On the other hand, the other three global S2S models demonstrated a low resolution and low reliability.
The limited-area forecasts (WRF 25 km and 5 km; Figure 15) showed good reliability and both had a better reliability than NCEP in week 3, which is the initialization period in which NCEP demonstrated some S2S prediction ability. More specifically, in week 3, the trendlines of WRF 25 km and 5 km were closer and more parallel to the diagonal even though they had fewer ensemble members (WRF: 14 members; NCEP: 105 members). Furthermore, WRF 5 km was slightly closer to the diagonal than WRF 25 km. Finally, both WRF forecasts showed a bias in underforecasting (trendline above the diagonal) the temperatures at 850 hPa. Regarding the maximum daily 2 m temperatures (Figure 16) predicted from week 3, both WRF forecasts presented a better reliability than in the previous initialization weeks, but the WRF 5 km (blue dashed line) trendline was more parallel to the diagonal than WRF 25 km for all the forecast probabilities.
Figure 17 illustrates the EDI of the various global perturbed and WRF forecasts relative to the T95 of the ERA5 climatology of May 1975–2005 for the different probability thresholds (PTs: 17%, 34%, etc., on the x-axis). The maximum (minimum) value of the EDI is equal to 1 (−1) when the HR (FAR) is equal to 1. In comparison with the previously mentioned scores (such as HR, FAR and CSI), EDI complements them by focusing on the tail of the temperature distribution, which is crucial for rare or extreme events like the heatwave of this study [88]. In week 3 (Figure 17), KMA had the best profile (highest EDIs), followed by ECMWF and NCEP at the lower PTs and WRF forecasts at the higher PTs. KMA reached an EDI = 0 for the PT of 100%. UKMO reached the same point, but on lower probabilities, its EDI was around 0. ECMWF and KMA maintained an EDI = 1 (perfect score) for PTs up to 34%. The forecasts of week 3 tended to lose their accuracy in higher PTs (i.e., for larger numbers of ensemble members) than the ones of week 1 and week 2.
Μost of the models showed an increase in prediction ability in week 3 and the main contributors to the strong heating were the temperature advection at 850 hPa and the downward motion due to the presence of the high geopotential heights at 500 hPa. In Figure 18, the reanalysis showed positive geopotential height anomalies over the Eastern Mediterranean, which were associated with anomalous anticyclonic flow in the lower troposphere and subsidence over Greece and western Turkey for the period from 8 to 23 May 2020. Many models (Figure 19) identified the positive–negative dipoles of the anomalies at 850 hPa above the Eastern Mediterranean and northeastern Europe for week 3 (initialization period). These models were ECMWF, UKMO, KMA, NCEP and WRF 25 km. The NCEP and WRF anomalies had many similarities because the initial/boundary conditions of WRF were derived from the NCEP/CFSv.2 model. However, a better-defined anomalous pattern was predicted by WRF compared to NCEP. At 500 hPa, the same models identified the positive anomalies above Greece and Turkey. On the other hand, HMCR and CMA did not detect the southwest anomalous flow over the central Mediterranean at 850 hPa. Also, the CMA forecasts from week 2, which exhibited very high AUC values (0.84) relative to the other models, did not detect the dipole anomaly, which aligns with the argument that 50% of the CMA ensemble members detected high temperatures above Greece due to randomness.

4. Discussion and Conclusions

The heatwave of this case study occurred abnormally early before summer. According to the various heat indices, it started on 11 May and ended on 20 May 2020, affecting the whole Greek region and the west coast of Turkey. The intensity was also high given the fact that the heatwave took place in mid-May. The highest 2 m temperatures recorded on 16 May were 40 °C in southern Greece and 42.6 °C in southwestern Turkey on 17 May. Additionally, the tropospheric temperature (850 hPa) and geopotential height (850 and 500 hPa) anomalies were substantially high for this period. More specifically, the temperatures at 850 hPa were ten to eleven standard deviations warmer than the climatology above the Aegean Sea and the geopotential heights at 850 and 500 hPa were, respectively, six to eight and ten to eleven standard deviations higher than the climatology above the southeast Aegean Sea and southwest Turkey. The warm air masses that affected the area of interest originated from the central Sahara region and moved anticyclonically, reaching Greece and Turkey. The main causes of this strong heating were the thermal advection from the southwest and the subsidence (leading to thermal compression) caused by the upper air ridge and the position of the jet stream.
The S2S prediction abilities of six global models (ECMWF, UKMO, CMA, KMA, NCEP and HMCR) and the limited-area model WRF with nested forecasts at grid spacings of 25 km and 5 km were investigated using different statistical scores. The main conclusions of this analysis are as follows:
  • Most global forecasts demonstrated the ability to predict the temperature at 850 hPa, 14 to 21 days (S2S timescale) before the event, which is in agreement with [28]. As also presented in a previous study [17], the best models were ECMWF, KMA and NCEP which showed a good prediction ability even three to four weeks before the heatwave. The lowest prediction performance for this event was shown by the HMCR model, which is likely due to the small number of vertical levels (28) employed by this model and the fact that it was not coupled to an ocean model.
  • The WRF model (25 km and 5 km) has approximately the same performance as ECMWF, KMA and NCEP, but WRF, especially WRF 5 km, was more reliable (its trendline was closer and parallel to the perfect reliability line).
  • Similar to the study of [24], the WRF model, in terms of AUC and reliability, performs better than NCEP, which provides the boundary conditions to the downscaling forecasts of WRF.
  • There was no substantial difference in the S2S predictability of temperature at 850 hPa when WRF 25 km and WRF 5 km were compared. Therefore, there is no systematic benefit from the dynamical downscaling of the 850 hPa temperature forecasts to a higher spatial resolution (5 km × 5 km).
  • The S2S prediction ability for the maximum daily 2 m temperatures was slightly better for the WRF 5 km compared to WRF 25 km. This difference was only observed in week 3 (lead time of 15–20 days) in terms of reliability, AUC and higher (and closer to ERA5) maximum temperatures at 2 m that were forecasted by WRF 5 km.
  • All models underestimated the extreme high temperatures at 850 hPa and near the surface (as far as WRF is concerned), but the ACC scores for the global models showed indications of better prediction abilities for strong temperature anomalies two to three weeks before the heatwave.
Future work related to the predictability of heatwaves includes numerical weather prediction experiments to investigate the effects of soil moisture, SSTs and North Atlantic oscillations on extreme temperatures at 2 m and in the lower troposphere.

Author Contributions

Conceptualization, D.M. and I.P.; methodology, D.M. and I.P.; software, D.M.; validation, D.M.; formal analysis, D.M.; investigation, D.M., I.P., P.Z. and C.A.; resources, D.M. and I.P.; data curation, D.M.; writing—original draft preparation, D.M.; writing—review and editing, D.M., I.P., P.Z. and C.A.; visualization, D.M.; supervision, I.P.; project administration, I.P.; funding acquisition, D.M. All authors have read and agreed to the published version of the manuscript.

Funding

The research work was supported by the Hellenic Foundation for Research and Innovation (HFRI) under the 3rd Call for HFRI PhD Fellowships (Fellowship Number: 5413).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data of the S2S forecasts of the six international meteorological centers (ECMWF, UKMO, KMA, CMA, NCEP, HMCR) presented in this study are available at https://apps.ecmwf.int/datasets/data/s2s-realtime-instantaneous-accum-ecmf/levtype=sfc/type=cf/ (accessed on 20 March 2024), The ERA5 and ERA5-Land data are available at the Copernicus Climate Data Store in https://cds.climate.copernicus.eu/ and the weather stations data are available at https://www.ncei.noaa.gov/cdo-web/datasets (accessed on 20 March 2024) and https://www.meteo.gr. The HYSPLIT (Internet based) trajectory model is available at https://www.ready.noaa.gov/HYSPLIT_traj.php (accessed on 20 March 2024). The code of the WRF model is available at https://www2.mmm.ucar.edu/wrf/users/download/get_source.html (accessed on 20 March 2024). The initial and boundary conditions (CFSv.2) for the WRF forecasts are available at https://www.ncei.noaa.gov/products/weather-climate-models/climate-forecast-system (accessed on 20 March 2024). The S2S forecasts of WRF are available on request from the first author (corresponding author) because they are part of his ongoing PhD research.

Acknowledgments

This work was based on the S2S data hosted at ECMWF. S2S, a joint initiative of the World Weather Research Programme (WWRP) and the World Climate Research Programme (WCRP). The results presented in this work have been produced using the Aristotle University of Thessaloniki (AUTh) High Performance Computing Infrastructure and Resources. We thank the Hellenic National Meteorological Service, the National Observatory of Athens and the ECMWF for making the surface observations available. For the ERA5 and ERA5-Land reanalyses, we thank the Copernicus Climate Data Store. The authors gratefully acknowledge the NOAA Air Resources Laboratory (ARL) for the provision of the HYSPLIT transport and dispersion model and READY website (https://www.ready.noaa.gov; accessed on 15 November 2023) used in this publication. This work utilized the WRF model developed and supported by the National Center for Atmospheric Research (NCAR). We acknowledge the use of the Climate Forecast System Version 2 (CFSv2) data provided by the National Centers for Environmental Prediction (NCEP) in the WRF forecasts. We thank the three anonymous reviewers for their constructive reviews.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, interpretation of the data; in the writing of the manuscript; and in the decision to publish the results.

Appendix A

Table A1. S2S model forecast initialization frequency during each week. Each check mark indicates the initialization day. The last column (Total) indicates the total number of ensemble members of each model per week. The number in brackets next to the name of each centre shows the ensemble members produced on each day that has a tick mark.
Table A1. S2S model forecast initialization frequency during each week. Each check mark indicates the initialization day. The last column (Total) indicates the total number of ensemble members of each model per week. The number in brackets next to the name of each centre shows the ensemble members produced on each day that has a tick mark.
Meteorological CentreMonTueWedThuFriSatSunTotal
ECMWF (50) 100
UKMO (3)21
CMA (3) 6
KMA (3)21
NCEP (15)105
HMCR (19) 19
AUTH-WRF (2)14
Figure A1. Anomaly Correlation Coefficients for six global S2S ECMWF, UKMO, CMA, KMA, NCEP and HMCR forecasts across different initialization dates (9 to 29 April). Unperturbed forecast (red), mean of perturbed forecasts (blue) and ACC = 0.6 (black). The ACC value for a centre is only displayed when there is a forecast and reforecast for the same initialization date.
Figure A1. Anomaly Correlation Coefficients for six global S2S ECMWF, UKMO, CMA, KMA, NCEP and HMCR forecasts across different initialization dates (9 to 29 April). Unperturbed forecast (red), mean of perturbed forecasts (blue) and ACC = 0.6 (black). The ACC value for a centre is only displayed when there is a forecast and reforecast for the same initialization date.
Atmosphere 15 00442 g0a1
Figure A2. Temperature (K) at 850 hPa (solid lines) from reanalysis (ERA5; pink), climatology (ERA5; orange), ECMWF ensemble mean of 50 perturbed forecasts (blue), ECMWF unperturbed forecast (red) and ECMWF climatology from the corresponding reforecasts (green). Temperature anomalies (K) at 850 hPa (dashed lines) in ECMWF ensemble mean minus ECMWF model’s climatology (blue), ECMWF unperturbed forecast minus ECMWF model’s climatology (red), ERA5 reanalysis minus ERA5 climatology (pink). Evaluation period: 8–23 May 2020. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Figure A2. Temperature (K) at 850 hPa (solid lines) from reanalysis (ERA5; pink), climatology (ERA5; orange), ECMWF ensemble mean of 50 perturbed forecasts (blue), ECMWF unperturbed forecast (red) and ECMWF climatology from the corresponding reforecasts (green). Temperature anomalies (K) at 850 hPa (dashed lines) in ECMWF ensemble mean minus ECMWF model’s climatology (blue), ECMWF unperturbed forecast minus ECMWF model’s climatology (red), ERA5 reanalysis minus ERA5 climatology (pink). Evaluation period: 8–23 May 2020. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Atmosphere 15 00442 g0a2
Figure A3. Temperature (K) at 850 hPa (solid lines) from reanalysis (ERA5—pink), climatology (ERA5—orange), UKMO ensemble mean of three perturbed forecasts (blue), UKMO unperturbed forecast (red) and UKMO climatology from the corresponding reforecasts (green). Temperature anomalies (K) at 850 hPa (dashed lines) in UKMO ensemble mean minus UKMO model’s climatology (blue), UKMO unperturbed forecast minus UKMO model’s climatology (red) and ERA5 reanalysis minus ERA5 climatology (pink). Evaluation period: 8–23 May 2020. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Figure A3. Temperature (K) at 850 hPa (solid lines) from reanalysis (ERA5—pink), climatology (ERA5—orange), UKMO ensemble mean of three perturbed forecasts (blue), UKMO unperturbed forecast (red) and UKMO climatology from the corresponding reforecasts (green). Temperature anomalies (K) at 850 hPa (dashed lines) in UKMO ensemble mean minus UKMO model’s climatology (blue), UKMO unperturbed forecast minus UKMO model’s climatology (red) and ERA5 reanalysis minus ERA5 climatology (pink). Evaluation period: 8–23 May 2020. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Atmosphere 15 00442 g0a3
Figure A4. Temperature (K) at 850 hPa (solid lines) from reanalysis (ERA5—pink), climatology (ERA5—orange), CMA ensemble mean of three perturbed forecasts (blue), CMA unperturbed forecast (red) and CMA climatology from the corresponding reforecasts (green). Temperature anomalies (K) at 850 hPa (dashed lines) in CMA ensemble mean minus CMA model’s climatology (blue), CMA unperturbed forecast minus CMA model’s climatology (red) and ERA5 reanalysis minus ERA5 climatology (pink). Evaluation period: 8–23 May 2020. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Figure A4. Temperature (K) at 850 hPa (solid lines) from reanalysis (ERA5—pink), climatology (ERA5—orange), CMA ensemble mean of three perturbed forecasts (blue), CMA unperturbed forecast (red) and CMA climatology from the corresponding reforecasts (green). Temperature anomalies (K) at 850 hPa (dashed lines) in CMA ensemble mean minus CMA model’s climatology (blue), CMA unperturbed forecast minus CMA model’s climatology (red) and ERA5 reanalysis minus ERA5 climatology (pink). Evaluation period: 8–23 May 2020. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Atmosphere 15 00442 g0a4
Figure A5. Temperature (K) at 850 hPa (solid lines) from reanalysis (ERA5—pink), climatology (ERA5—orange), KMA ensemble mean of three perturbed forecasts (blue), KMA unperturbed forecast (red) and KMA climatology from the corresponding reforecasts (green). Temperature anomalies (K) at 850 hPa (dashed lines) in KMA ensemble mean minus KMA model’s climatology (blue), KMA unperturbed forecast minus KMA model’s climatology (red) and ERA5 reanalysis minus ERA5 climatology (pink). Evaluation period: 8–23 May 2020. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Figure A5. Temperature (K) at 850 hPa (solid lines) from reanalysis (ERA5—pink), climatology (ERA5—orange), KMA ensemble mean of three perturbed forecasts (blue), KMA unperturbed forecast (red) and KMA climatology from the corresponding reforecasts (green). Temperature anomalies (K) at 850 hPa (dashed lines) in KMA ensemble mean minus KMA model’s climatology (blue), KMA unperturbed forecast minus KMA model’s climatology (red) and ERA5 reanalysis minus ERA5 climatology (pink). Evaluation period: 8–23 May 2020. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Atmosphere 15 00442 g0a5
Figure A6. Temperature (K) at 850 hPa (solid lines) from reanalysis (ERA5—pink), climatology (ERA5—orange), NCEP ensemble mean of 15 perturbed forecasts (blue), NCEP unperturbed forecast (red) and NCEP climatology from the corresponding reforecasts (green). Temperature anomalies (K) at 850 hPa (dashed lines) in NCEP ensemble mean minus NCEP model’s climatology (blue), NCEP unperturbed forecast minus NCEP model’s climatology (red) and ERA5 reanalysis minus ERA5 climatology (pink). Evaluation period: 8–23 May 2020. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Figure A6. Temperature (K) at 850 hPa (solid lines) from reanalysis (ERA5—pink), climatology (ERA5—orange), NCEP ensemble mean of 15 perturbed forecasts (blue), NCEP unperturbed forecast (red) and NCEP climatology from the corresponding reforecasts (green). Temperature anomalies (K) at 850 hPa (dashed lines) in NCEP ensemble mean minus NCEP model’s climatology (blue), NCEP unperturbed forecast minus NCEP model’s climatology (red) and ERA5 reanalysis minus ERA5 climatology (pink). Evaluation period: 8–23 May 2020. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Atmosphere 15 00442 g0a6aAtmosphere 15 00442 g0a6bAtmosphere 15 00442 g0a6cAtmosphere 15 00442 g0a6d
Figure A7. Temperature (K) at 850 hPa (solid lines) from reanalysis (ERA5—pink), climatology (ERA5—orange), HMCR ensemble mean of 19 perturbed forecasts (blue), HMCR unperturbed forecast (red) and HMCR climatology from the corresponding reforecasts (green). Temperature anomalies (K) at 850 hPa (dashed lines) in HMCR ensemble mean minus HMCR model’s climatology (blue), HMCR unperturbed forecast minus HMCR model’s climatology (red) and ERA5 reanalysis minus ERA5 climatology (pink). Evaluation period: 8–23 May 2020. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Figure A7. Temperature (K) at 850 hPa (solid lines) from reanalysis (ERA5—pink), climatology (ERA5—orange), HMCR ensemble mean of 19 perturbed forecasts (blue), HMCR unperturbed forecast (red) and HMCR climatology from the corresponding reforecasts (green). Temperature anomalies (K) at 850 hPa (dashed lines) in HMCR ensemble mean minus HMCR model’s climatology (blue), HMCR unperturbed forecast minus HMCR model’s climatology (red) and ERA5 reanalysis minus ERA5 climatology (pink). Evaluation period: 8–23 May 2020. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Atmosphere 15 00442 g0a7

References

  1. Changnon, S.A.; Kunkel, K.E.; Reinke, B.C. Impacts and Responses to the 1995 Heat Wave: A Call to Action. Bull. Am. Meteorol. Soc. 1996, 77, 1497–1506. [Google Scholar] [CrossRef]
  2. Klinenberg, E. Heat Wave: A Social Autopsy of Disaster in Chicago; University of Chicago Press: Chicago, IL, USA, 2015; ISBN 978-0-226-27621-2. [Google Scholar]
  3. Lowe, R.; Ballester, J.; Creswick, J.; Robine, J.-M.; Herrmann, F.R.; Rodó, X. Evaluating the Performance of a Climate-Driven Mortality Model during Heat Waves and Cold Spells in Europe. Int. J. Environ. Res. Public Health 2015, 12, 1279–1294. [Google Scholar] [CrossRef] [PubMed]
  4. Marx, W.; Haunschild, R.; Bornmann, L. Heat Waves: A Hot Topic in Climate Change Research. Theor. Appl. Clim. 2021, 146, 781–800. [Google Scholar] [CrossRef] [PubMed]
  5. Perkins, S.E. A Review on the Scientific Understanding of Heatwaves—Their Measurement, Driving Mechanisms, and Changes at the Global Scale. Atmos. Res. 2015, 164–165, 242–267. [Google Scholar] [CrossRef]
  6. Miller, S.; Chua, K.; Coggins, J.; Mohtadi, H. Heat Waves, Climate Change, and Economic Output. J. Eur. Econ. Assoc. 2021, 19, 2658–2694. [Google Scholar] [CrossRef]
  7. Mitchell, D.; Heaviside, C.; Vardoulakis, S.; Huntingford, C.; Masato, G.; Guillod, B.P.; Frumhoff, P.; Bowery, A.; Wallom, D.; Allen, M. Attributing Human Mortality during Extreme Heat Waves to Anthropogenic Climate Change. Environ. Res. Lett. 2016, 11, 074006. [Google Scholar] [CrossRef]
  8. O’Neill, M.S.; Ebi, K.L. Temperature Extremes and Health: Impacts of Climate Variability and Change in the United States. J. Occup. Environ. Med. 2009, 51, 13–25. [Google Scholar] [CrossRef] [PubMed]
  9. Christidis, N.; Mitchell, D.; Stott, P.A. Rapidly Increasing Likelihood of Exceeding 50 °C in Parts of the Mediterranean and the Middle East Due to Human Influence. NPJ Clim. Atmos. Sci. 2023, 6, 45. [Google Scholar] [CrossRef]
  10. Giorgi, F. Climate Change Hot-Spots. Geophys. Res. Lett. 2006, 33, 1–4. [Google Scholar] [CrossRef]
  11. Zittis, G.; Hadjinicolaou, P.; Klangidou, M.; Proestos, Y.; Lelieveld, J. A Multi-Model, Multi-Scenario, and Multi-Domain Analysis of Regional Climate Projections for the Mediterranean. Reg. Environ. Chang. 2019, 19, 2621–2635. [Google Scholar] [CrossRef]
  12. Coppola, E.; Raffaele, F.; Giorgi, F.; Giuliani, G.; Xuejie, G.; Ciarlo, J.M.; Sines, T.R.; Torres-Alavez, J.A.; Das, S.; di Sante, F.; et al. Climate Hazard Indices Projections Based on CORDEX-CORE, CMIP5 and CMIP6 Ensemble. Clim. Dyn. 2021, 57, 1293–1383. [Google Scholar] [CrossRef]
  13. Georgoulias, A.K.; Akritidis, D.; Kalisoras, A.; Kapsomenakis, J.; Melas, D.; Zerefos, C.S.; Zanis, P. Climate Change Projections for Greece in the 21st Century from High-Resolution EURO-CORDEX RCM Simulations. Atmos. Res. 2022, 271, 106049. [Google Scholar] [CrossRef]
  14. Bauer, P.; Thorpe, A.; Brunet, G. The Quiet Revolution of Numerical Weather Prediction. Nature 2015, 525, 47–55. [Google Scholar] [CrossRef] [PubMed]
  15. Mariotti, A.; Ruti, P.M.; Rixen, M. Progress in Subseasonal to Seasonal Prediction through a Joint Weather and Climate Community Effort. NPJ Clim. Atmos. Sci. 2018, 1, 4. [Google Scholar] [CrossRef]
  16. Ford, T.W.; Dirmeyer, P.A.; Benson, D.O. Evaluation of Heat Wave Forecasts Seamlessly across Subseasonal Timescales. NPJ Clim. Atmos. Sci. 2018, 1, 20. [Google Scholar] [CrossRef]
  17. Vitart, F.; Robertson, A.W. The Sub-Seasonal to Seasonal Prediction Project (S2S) and the Prediction of Extreme Events. NPJ Clim. Atmos. Sci. 2018, 1, 3. [Google Scholar] [CrossRef]
  18. Koster, R.D.; Mahanama, S.P.P.; Yamada, T.J.; Balsamo, G.; Berg, A.A.; Boisserie, M.; Dirmeyer, P.A.; Doblas-Reyes, F.J.; Drewitt, G.; Gordon, C.T.; et al. The Second Phase of the Global Land–Atmosphere Coupling Experiment: Soil Moisture Contributions to Subseasonal Forecast Skill. J. Hydrometeorol. 2011, 12, 805–822. [Google Scholar] [CrossRef]
  19. Chen, M.; Wang, W.; Kumar, A. Prediction of Monthly-Mean Temperature: The Roles of Atmospheric and Land Initial Conditions and Sea Surface Temperature. J. Clim. 2010, 23, 717–725. [Google Scholar] [CrossRef]
  20. White, C.J.; Carlsen, H.; Robertson, A.W.; Klein, R.J.T.; Lazo, J.K.; Kumar, A.; Vitart, F.; Coughlan de Perez, E.; Ray, A.J.; Murray, V.; et al. Potential Applications of Subseasonal-to-Seasonal (S2S) Predictions. Meteorol. Appl. 2017, 24, 315–325. [Google Scholar] [CrossRef]
  21. Shapiro, M.; Shukla, J.; Brunet, G.; Nobre, C.; Béland, M.; Dole, R.; Trenberth, K.; Anthes, R.; Asrar, G.; Barrie, L.; et al. An Earth-System Prediction Initiative for the Twenty-First Century. Bull. Am. Meteorol. Soc. 2010, 91, 1377–1388. [Google Scholar] [CrossRef]
  22. Brunet, G.; Shapiro, M.; Hoskins, B.; Moncrieff, M.; Dole, R.; Kiladis, G.N.; Kirtman, B.; Lorenc, A.; Mills, B.; Morss, R.; et al. Collaboration of the Weather and Climate Communities to Advance Subseasonal-to-Seasonal Prediction. Bull. Am. Meteorol. Soc. 2010, 91, 1397–1406. [Google Scholar] [CrossRef]
  23. Ford, T.; Dirmeyer, P. Evaluation of the Ability of S2S and NMME Models to Predict Heat Waves Following Drought Events in the United States. AGU Fall Meet. Abstr. 2016, 2016, A43O-06. [Google Scholar]
  24. Oh, J.; Oh, J.; Huh, M. A Case Study of Bias Correction and the Dynamical Downscaling of CFSv2 S2S Forecasts Using a WRF Model: Heatwave in 2018 over South Korea. Atmosphere 2022, 13, 2086. [Google Scholar] [CrossRef]
  25. Thanigachalam, A.; AchutaRao, K.; Mitra, A.K.; Ashrit, R.; Gupta, A. Extended-Range Prediction of Heatwave Events over North India: Role of Atmospheric Blocking over North Atlantic. EGU Gen. Assem. Conf. Abstr. 2020, EGU2020-750. [Google Scholar] [CrossRef]
  26. White, R.H.; Kornhuber, K.; Martius, O.; Wirth, V. From Atmospheric Waves to Heatwaves: A Waveguide Perspective for Understanding and Predicting Concurrent, Persistent, and Extreme Extratropical Weather. Bull. Am. Meteorol. Soc. 2022, 103, E923–E935. [Google Scholar] [CrossRef]
  27. Kueh, M.-T.; Lin, C.-Y. The 2018 Summer Heatwaves over Northwestern Europe and Its Extended-Range Prediction. Sci. Rep. 2020, 10, 19283. [Google Scholar] [CrossRef] [PubMed]
  28. Lin, H.; Mo, R.; Vitart, F. The 2021 Western North American Heatwave and Its Subseasonal Predictions. Geophys. Res. Lett. 2022, 49, e2021GL097036. [Google Scholar] [CrossRef]
  29. Xie, J.; Yu, J.; Chen, H.; Hsu, P.-C. Sources of Subseasonal Prediction Skill for Heatwaves over the Yangtze River Basin Revealed from Three S2S Models. Adv. Atmos. Sci. 2020, 37, 1435–1450. [Google Scholar] [CrossRef]
  30. Vitart, F.; Ardilouze, C.; Bonet, A.; Brookshaw, A.; Chen, M.; Codorean, C.; Déqué, M.; Ferranti, L.; Fucile, E.; Fuentes, M.; et al. The Subseasonal to Seasonal (S2S) Prediction Project Database. Bull. Am. Meteorol. Soc. 2017, 98, 163–173. [Google Scholar] [CrossRef]
  31. Miloshevich, G.; Cozian, B.; Abry, P.; Borgnat, P.; Bouchet, F. Probabilistic Forecasts of Extreme Heatwaves Using Convolutional Neural Networks in a Regime of Lack of Data. Phys. Rev. Fluids 2023, 8, 040501. [Google Scholar] [CrossRef]
  32. Weyn, J.A.; Durran, D.R.; Caruana, R. Improving Data-Driven Global Weather Prediction Using Deep Convolutional Neural Networks on a Cubed Sphere. J. Adv. Model. Earth Syst. 2020, 12, e2020MS002109. [Google Scholar] [CrossRef]
  33. Drouard, M.; Pérez-Aracil, J.; Barriopedro, D.; Zaninelli, P.G.; Garrido-Perez, J.M.; Fister, D.; Salcedo-Sanz, S.; García-Herrera, R. S2S Prediction of Summer Heatwaves in the Iberian Peninsula Using Convolutional Networks. EGU Gen. Assem. Conf. Abstr. 2023, EGU23-5742. [Google Scholar]
  34. Lowe, R.; García-Díez, M.; Ballester, J.; Creswick, J.; Robine, J.-M.; Herrmann, F.R.; Rodó, X. Evaluation of an Early-Warning System for Heat Wave-Related Mortality in Europe: Implications for Sub-Seasonal to Seasonal Forecasting and Climate Services. Int. J. Environ. Res. Public Health 2016, 13, 206. [Google Scholar] [CrossRef] [PubMed]
  35. Tovar-Sánchez, A.; Sánchez-Quiles, D.; Rodríguez-Romero, A. Massive Coastal Tourism Influx to the Mediterranean Sea: The Environmental Risk of Sunscreens. Sci. Total Environ. 2019, 656, 316–321. [Google Scholar] [CrossRef] [PubMed]
  36. Knorr, W.; Pytharoulis, I.; Petropoulos, G.P.; Gobron, N. Combined Use of Weather Forecasting and Satellite Remote Sensing Information for Fire Risk, Fire and Fire Impact Monitoring. Comput. Ecol. Softw. 2011, 1, 112. [Google Scholar]
  37. Koutsias, N.; Arianoutsou, M.; Kallimanis, A.S.; Mallinis, G.; Halley, J.M.; Dimopoulos, P. Where Did the Fires Burn in Peloponnisos, Greece the Summer of 2007? Evidence for a Synergy of Fuel and Weather. Agric. For. Meteorol. 2012, 156, 41–53. [Google Scholar] [CrossRef]
  38. Fader, M.; Shi, S.; Von Bloh, W.; Bondeau, A.; Cramer, W. Mediterranean Irrigation under Climate Change: More Efficient Irrigation Needed to Compensate Increases in Irrigation Water Requirements. Hydrol. Earth Syst. Sci. Discuss. 2015, 12, 8459–8504. [Google Scholar] [CrossRef]
  39. Khodayar Pardo, S.; Paredes-Fortuny, L. Uneven Evolution of Regional European Summer Heatwaves under Climate Change. Weather. Clim. Extrem. 2024, 43, 100648. [Google Scholar] [CrossRef]
  40. Domeisen, D.I.V.; Eltahir, E.A.B.; Fischer, E.M.; Knutti, R.; Perkins-Kirkpatrick, S.E.; Schär, C.; Seneviratne, S.I.; Weisheimer, A.; Wernli, H. Prediction and Projection of Heatwaves. Nat. Rev. Earth Environ. 2023, 4, 36–50. [Google Scholar] [CrossRef]
  41. EPA. Seasonality and Climate Change. Available online: https://www.epa.gov/climate-indicators/seasonality-and-climate-change (accessed on 15 December 2023).
  42. Founda, D.; Varotsos, K.V.; Pierros, F.; Giannakopoulos, C. Observed and Projected Shifts in Hot Extremes’ Season in the Eastern Mediterranean. Glob. Planet. Chang. 2019, 175, 190–200. [Google Scholar] [CrossRef]
  43. Founda, D.; Giannakopoulos, C. The Exceptionally Hot Summer of 2007 in Athens, Greece—A Typical Summer in the Future Climate? Glob. Planet. Chang. 2009, 67, 227–236. [Google Scholar] [CrossRef]
  44. Galanaki, E.; Giannaros, C.; Kotroni, V.; Lagouvardos, K.; Papavasileiou, G. Spatio-Temporal Analysis of Heatwaves Characteristics in Greece from 1950 to 2020. Climate 2023, 11, 5. [Google Scholar] [CrossRef]
  45. Matzarakis, A.; Nastos, P.T. Human-Biometeorological Assessment of Heat Waves in Athens. Theor. Appl. Climatol. 2011, 105, 99–106. [Google Scholar] [CrossRef]
  46. Heat-Waves: Risks and Responses. Available online: https://www.who.int/publications-detail-redirect/9789289010948 (accessed on 9 August 2023).
  47. Saha, S.; Moorthi, S.; Wu, X.; Wang, J.; Nadiga, S.; Tripp, P.; Behringer, D.; Hou, Y.-T.; Chuang, H.; Iredell, M.; et al. The NCEP Climate Forecast System Version 2. J. Clim. 2014, 27, 2185–2208. [Google Scholar] [CrossRef]
  48. Hersbach, H.; Bell, B.; Berrisford, P.; Hirahara, S.; Horányi, A.; Muñoz-Sabater, J.; Nicolas, J.; Peubey, C.; Radu, R.; Schepers, D.; et al. The ERA5 Global Reanalysis. Q. J. R. Meteorol. Soc. 2020, 146, 1999–2049. [Google Scholar] [CrossRef]
  49. Muñoz-Sabater, J.; Dutra, E.; Agustí-Panareda, A.; Albergel, C.; Arduini, G.; Balsamo, G.; Boussetta, S.; Choulga, M.; Harrigan, S.; Hersbach, H.; et al. ERA5-Land: A State-of-the-Art Global Reanalysis Dataset for Land Applications. Earth Syst. Sci. Data 2021, 13, 4349–4383. [Google Scholar] [CrossRef]
  50. Dong, J.; Peng, J.; He, X.; Corcoran, J.; Qiu, S.; Wang, X. Heatwave-Induced Human Health Risk Assessment in Megacities Based on Heat Stress-Social Vulnerability-Human Exposure Framework. Landsc. Urban Plan. 2020, 203, 103907. [Google Scholar] [CrossRef]
  51. Napoli, C.D.; Pappenberger, F.; Cloke, H.L. Verification of Heat Stress Thresholds for a Health-Based Heat-Wave Definition. J. Appl. Meteorol. Climatol. 2019, 58, 1177–1194. [Google Scholar] [CrossRef]
  52. Hahn, G.L. Dynamic Responses of Cattle to Thermal Heat Loads. J. Anim. Sci. 1999, 77, 10–20. [Google Scholar] [CrossRef]
  53. Vitali, A.; Felici, A.; Esposito, S.; Bernabucci, U.; Bertocchi, L.; Maresca, C.; Nardone, A.; Lacetera, N. The Effect of Heat Waves on Dairy Cow Mortality. J. Dairy Sci. 2015, 98, 4572–4579. [Google Scholar] [CrossRef]
  54. Yılmaz, M.; Kara, Y.; Çulpan, H.C.; Can, G.; Toros, H. Detection and Regional Analysis of Heatwave Characteristics in İstanbul. Sustain. Cities Soc. 2023, 97, 104789. [Google Scholar] [CrossRef]
  55. Dunn, R.J.H.; Mead, N.E.; Willett, K.M.; Parker, D.E. Analysis of Heat Stress in UK Dairy Cattle and Impact on Milk Yields. Environ. Res. Lett. 2014, 9, 064006. [Google Scholar] [CrossRef]
  56. Russo, S.; Dosio, A.; Graversen, R.G.; Sillmann, J.; Carrao, H.; Dunbar, M.B.; Singleton, A.; Montagna, P.; Barbola, P.; Vogt, J.V. Magnitude of Extreme Heat Waves in Present Climate and Their Projection in a Warming World. J. Geophys. Res. Atmos. 2014, 119, 12500–12512. [Google Scholar] [CrossRef]
  57. Raymond, C.; Matthews, T.; Horton, R.M. The Emergence of Heat and Humidity Too Severe for Human Tolerance. Sci. Adv. 2020, 6, eaaw1838. [Google Scholar] [CrossRef] [PubMed]
  58. Schoen, C. A New Empirical Model of the Temperature–Humidity Index. J. Appl. Meteorol. Climatol. 2005, 44, 1413–1420. [Google Scholar] [CrossRef]
  59. Steadman, R.G. The Assessment of Sultriness. Part I: A Temperature-Humidity Index Based on Human Physiology and Clothing Science. J. Appl. Meteorol. Climatol. 1979, 18, 861–873. [Google Scholar] [CrossRef]
  60. Nairn, J.R.; Fawcett, R.J.B. The Excess Heat Factor: A Metric for Heatwave Intensity and Its Use in Classifying Heatwave Severity. Int. J. Environ. Res. Public Health 2015, 12, 227–253. [Google Scholar] [CrossRef]
  61. Karacostas, T.S.; Downing, T.E. Heat Wave Events in a Changing Climate. In Climate Change and Extreme Events: Altered Risk, Socio-Economic Impacts and Policy Responses; Institute for Environmental Studies: Amsterdam, The Netherlands, 1996; pp. 139–156. [Google Scholar]
  62. Tolika, K. Assessing Heat Waves over Greece Using the Excess Heat Factor (EHF). Climate 2019, 7, 9. [Google Scholar] [CrossRef]
  63. Jendritzky, G.; de Dear, R.; Havenith, G. UTCI—Why Another Thermal Index? Int. J. Biometeorol. 2012, 56, 421–428. [Google Scholar] [CrossRef]
  64. Hatvani-Kovacs, G.; Belusko, M.; Pockett, J.; Boland, J. Can the Excess Heat Factor Indicate Heatwave-Related Morbidity? A Case Study in Adelaide, South Australia. EcoHealth 2016, 13, 100–110. [Google Scholar] [CrossRef]
  65. Langlois, N.; Herbst, J.; Mason, K.; Nairn, J.; Byard, R.W. Using the Excess Heat Factor (EHF) to Predict the Risk of Heat Related Deaths. J. Forensic Leg. Med. 2013, 20, 408–411. [Google Scholar] [CrossRef] [PubMed]
  66. Varghese, B.M.; Hansen, A.; Nitschke, M.; Nairn, J.; Hanson-Easey, S.; Bi, P.; Pisaniello, D. Heatwave and Work-Related Injuries and Illnesses in Adelaide, Australia: A Case-Crossover Analysis Using the Excess Heat Factor (EHF) as a Universal Heatwave Index. Int. Arch. Occup. Environ. Health 2019, 92, 263–272. [Google Scholar] [CrossRef] [PubMed]
  67. Nairn, J.; Ostendorf, B.; Bi, P. Performance of Excess Heat Factor Severity as a Global Heatwave Health Impact Index. Int. J. Environ. Res. Public Health 2018, 15, 2494. [Google Scholar] [CrossRef]
  68. Piticar, A.; Croitoru, A.-E.; Ciupertea, F.-A.; Harpa, G.-V. Recent Changes in Heat Waves and Cold Waves Detected Based on Excess Heat Factor and Excess Cold Factor in Romania. Int. J. Climatol. 2018, 38, 1777–1793. [Google Scholar] [CrossRef]
  69. Schlatter, T.W. Temperature-Humidity Index. In Climatology; Encyclopedia of Earth Science; Springer: Boston, MA, USA, 1987; pp. 837–838. [Google Scholar]
  70. Thom, E.C. The Discomfort Index. Weatherwise 1959, 12, 57–61. [Google Scholar] [CrossRef]
  71. Gogoașe-Nistoran, D.-E.; Antohe, D.-M.; Opriș, I.; Ionescu, C.-S. A Method to Analyse Urban Heatwaves. Case Study—Bucharest City (Romania). E3S Web. Conf. 2021, 286, 04004. [Google Scholar] [CrossRef]
  72. Bröde, P.; Fiala, D.; Błażejczyk, K.; Holmér, I.; Jendritzky, G.; Kampmann, B.; Tinz, B.; Havenith, G. Deriving the Operational Procedure for the Universal Thermal Climate Index (UTCI). Int. J. Biometeorol. 2012, 56, 481–494. [Google Scholar] [CrossRef] [PubMed]
  73. Jendritzky, G. The Universal Thermal Climate Index for the Thermo-Physiologically Significant Assessment of the Atmospheric Environment. In Proceedings of the 3rd Symposium Urban Climatology, Davis, CA, USA, 15 August 2000; pp. 9–13. [Google Scholar]
  74. Skamarock, C.; Klemp, B.; Dudhia, J.; Gill, O.; Liu, Z.; Berner, J.; Wang, W.; Powers, G.; Duda, G.; Barker, D.; et al. A Description of the Advanced Research WRF Model Version 4.1. NCAR Tech. Note 2019. [Google Scholar] [CrossRef]
  75. Thompson, G.; Field, P.R.; Rasmussen, R.M.; Hall, W.D. Explicit Forecasts of Winter Precipitation Using an Improved Bulk Microphysics Scheme. Part II: Implementation of a New Snow Parameterization. Mon. Weather. Rev. 2008, 136, 5095–5115. [Google Scholar] [CrossRef]
  76. Grell, G.A.; Freitas, S.R. A Scale and Aerosol Aware Stochastic Convective Parameterization for Weather and Air Quality Modeling. Atmos. Chem. Phys. 2014, 14, 5233–5250. [Google Scholar] [CrossRef]
  77. Iacono, M.J.; Delamere, J.S.; Mlawer, E.J.; Shephard, M.W.; Clough, S.A.; Collins, W.D. Radiative Forcing by Long-Lived Greenhouse Gases: Calculations with the AER Radiative Transfer Models. J. Geophys. Res. Atmos. 2008, 113, 1–8. [Google Scholar] [CrossRef]
  78. Janjić, Z.I. The Step-Mountain Eta Coordinate Model: Further Developments of the Convection, Viscous Sublayer, and Turbulence Closure Schemes. Mon. Weather. Rev. 1994, 122, 927–945. [Google Scholar] [CrossRef]
  79. Tewari, M.; Chen, F.; Wang, W.; Dudhia, J.; LeMone, M.A.; Mitchell, K.; Ek, M.; Gayno, G.; Wegiel, J.; Cuenca, R.H. Implementation and Verification of the Unified NOAH Land Surface Model in the WRF Model. In Proceedings of the 20th Conference on Weather Analysis and Forecasting/16th Conference on Numerical Weather Prediction, Seattle, WA, USA, 14 January 2004; Volume 1115, pp. 2165–2170. [Google Scholar]
  80. Monin, A.S.; Obukhov, A.M. Basic Laws of Turbulent Mixing in the Surface Layer of the Atmosphere. Contrib. Geophys. Inst. Acad. Sci. USSR 1954, 151, e187. [Google Scholar]
  81. Miyakoda, K.; Hembree, G.D.; Strickler, R.F.; Shulman, I. Cumulative Results of Extended Forecast Experiments I. Model Performance for Winter Cases. Mon. Weather. Rev. 1972, 100, 836–855. [Google Scholar] [CrossRef]
  82. Wilks, D.S. Statistical Methods in the Atmospheric Sciences; Academic Press: Cantabrigian, MA, USA, 2011; ISBN 978-0-12-385022-5. [Google Scholar]
  83. Marden, J.I. Positions and QQ Plots. Stat. Sci. 2004, 19, 606–614. [Google Scholar] [CrossRef]
  84. Fan, J.; Upadhye, S.; Worster, A. Understanding Receiver Operating Characteristic (ROC) Curves. Can. J. Emerg. Med. 2006, 8, 19–20. [Google Scholar] [CrossRef] [PubMed]
  85. Hsu, W.; Murphy, A.H. The Attributes Diagram A Geometrical Framework for Assessing the Quality of Probability Forecasts. Int. J. Forecast. 1986, 2, 285–293. [Google Scholar] [CrossRef]
  86. Mason, I. A Model for Assessment of Weather Forecasts. Aust. Meteor. Mag. 1982, 30, 291–303. [Google Scholar]
  87. Jolliffe, I.T.; Stephenson, D.B. Forecast Verification: A Practitioner’s Guide in Atmospheric Science; John Wiley & Sons: Hoboken, NJ, USA, 2012; ISBN 978-1-119-96107-9. [Google Scholar]
  88. Ferro, C.A.T.; Stephenson, D.B. Extremal Dependence Indices: Improved Verification Measures for Deterministic Forecasts of Rare Binary Events. Weather. Forecast. 2011, 26, 699–713. [Google Scholar] [CrossRef]
  89. Stein, A.F.; Draxler, R.R.; Rolph, G.D.; Stunder, B.J.B.; Cohen, M.D.; Ngan, F. NOAA’s HYSPLIT Atmospheric Transport and Dispersion Modeling System. Bull. Am. Meteorol. Soc. 2015, 96, 2059–2077. [Google Scholar] [CrossRef]
  90. Mitropoulos, D.; Pytharoulis, I.; Zanis, P.; Anagnostopoulou, C. Subseasonal-to-Seasonal Predictability Assessment of an Early Heat Wave in the Eastern Mediterranean in May 2020. Environ. Sci. Proc. 2023, 26, 42. [Google Scholar] [CrossRef]
  91. Tolika, K.; Maheras, P.; Pytharoulis, I.; Anagnostopoulou, C. The Anomalous Low and High Temperatures of 2012 over Greece—An Explanation from a Meteorological and Climatological Perspective. Nat. Hazards Earth Syst. Sci. 2014, 14, 501–507. [Google Scholar] [CrossRef]
  92. Tolika, K.; Pytharoulis, I.; Maheras, P. The Anomalous High Temperatures of November 2010 over Greece: Meteorological and Climatological Aspects. Nat. Hazards Earth Syst. Sci. 2011, 11, 2705–2714. [Google Scholar] [CrossRef]
  93. Brikas, D.P.; Karacostas, T.S.; Pennas, P.J.; Flocas, A.A. The Role of the Subtropical Jet Stream during Heat Wave Events over North-Central Greece. Meteorol. Atmos. Phys. 2006, 94, 219–233. [Google Scholar] [CrossRef]
  94. Fischer, E.M.; Seneviratne, S.I.; Vidale, P.L.; Lüthi, D.; Schär, C. Soil Moisture–Atmosphere Interactions during the 2003 European Summer Heat Wave. J. Clim. 2007, 20, 5081–5099. [Google Scholar] [CrossRef]
  95. Herold, N.; Kala, J.; Alexander, L.V. The Influence of Soil Moisture Deficits on Australian Heatwaves. Environ. Res. Lett. 2016, 11, 064003. [Google Scholar] [CrossRef]
  96. Stéfanon, M.; Drobinski, P.; D’Andrea, F.; Lebeaupin-Brossier, C.; Bastin, S. Soil Moisture-Temperature Feedbacks at Meso-Scale during Summer Heat Waves over Western Europe. Clim. Dyn. 2014, 42, 1309–1324. [Google Scholar] [CrossRef]
  97. Zhang, J.; Liu, Z.; Chen, L. Reduced Soil Moisture Contributes to More Intense and More Frequent Heat Waves in Northern China. Adv. Atmos. Sci. 2015, 32, 1197–1207. [Google Scholar] [CrossRef]
  98. Lorenz, R.; Jaeger, E.B.; Seneviratne, S.I. Persistence of Heat Waves and Its Link to Soil Moisture Memory. Geophys. Res. Lett. 2010, 37, 1–5. [Google Scholar] [CrossRef]
  99. Hauser, M.; Orth, R.; Seneviratne, S.I. Role of Soil Moisture versus Recent Climate Change for the 2010 Heat Wave in Western Russia. Geophys. Res. Lett. 2016, 43, 2819–2826. [Google Scholar] [CrossRef]
  100. Hirschi, M.; Seneviratne, S.I.; Alexandrov, V.; Boberg, F.; Boroneant, C.; Christensen, O.B.; Formayer, H.; Orlowsky, B.; Stepanek, P. Observational Evidence for Soil-Moisture Impact on Hot Extremes in Southeastern Europe. Nature Geosci. 2011, 4, 17–21. [Google Scholar] [CrossRef]
  101. García-Herrera, R.; Díaz, J.; Trigo, R.M.; Luterbacher, J.; Fischer, E.M. A Review of the European Summer Heat Wave of 2003. Crit. Rev. Environ. Sci. Technol. 2010, 40, 267–306. [Google Scholar] [CrossRef]
  102. Dong, L.; Mitra, C.; Greer, S.; Burt, E. The Dynamical Linkage of Atmospheric Blocking to Drought, Heatwave and Urban Heat Island in Southeastern US: A Multi-Scale Case Study. Atmosphere 2018, 9, 33. [Google Scholar] [CrossRef]
  103. Dole, R.; Hoerling, M.; Perlwitz, J.; Eischeid, J.; Pegion, P.; Zhang, T.; Quan, X.-W.; Xu, T.; Murray, D. Was There a Basis for Anticipating the 2010 Russian Heat Wave? Geophys. Res. Lett. 2011, 38, 1–5. [Google Scholar] [CrossRef]
  104. Neal, E.; Huang, C.S.Y.; Nakamura, N. The 2021 Pacific Northwest Heat Wave and Associated Blocking: Meteorology and the Role of an Upstream Cyclone as a Diabatic Source of Wave Activity. Geophys. Res. Lett. 2022, 49, e2021GL097699. [Google Scholar] [CrossRef]
  105. Miralles, D.G.; Teuling, A.J.; van Heerwaarden, C.C.; Vilà-Guerau de Arellano, J. Mega-Heatwave Temperatures Due to Combined Soil Desiccation and Atmospheric Heat Accumulation. Nature Geosci. 2014, 7, 345–349. [Google Scholar] [CrossRef]
  106. Feudale, L.; Shukla, J. Role of Mediterranean SST in Enhancing the European Heat Wave of Summer 2003. Geophys. Res. Lett. 2007, 34, 1–4. [Google Scholar] [CrossRef]
  107. Mass, C. The Uncoordinated Giant II: Why U.S. Operational Numerical Weather Prediction Is Still Lagging and How to Fix It. Bull. Am. Meteorol. Soc. 2023, 104, E851–E871. [Google Scholar] [CrossRef]
  108. Marzban, C. The ROC Curve and the Area under It as Performance Measures. Weather. Forecast. 2004, 19, 1106–1114. [Google Scholar] [CrossRef]
  109. Murphy, J.M. The Impact of Ensemble Forecasts on Predictability. Q. J. R. Meteorol. Soc. 1988, 114, 463–493. [Google Scholar] [CrossRef]
  110. Kumar, A.; Chen, M. Inherent Predictability, Requirements on the Ensemble Size, and Complementarity. Mon. Weather. Rev. 2015, 143, 3192–3203. [Google Scholar] [CrossRef]
  111. Branković, Č.; Palmer, T.N. Atmospheric Seasonal Predictability and Estimates of Ensemble Size. Mon. Weather. Rev. 1997, 125, 859–874. [Google Scholar] [CrossRef]
  112. Hartmann, H.C.; Pagano, T.C.; Sorooshian, S.; Bales, R. CONFIDENCE BUILDERS: Evaluating Seasonal Climate Forecasts from User Perspectives. Bull. Am. Meteorol. Soc. 2002, 83, 683–698. [Google Scholar] [CrossRef]
Figure 1. The two nested domains of the WRF simulations, with horizontal grid spacings of 25 × 25 km in Domain 1 (d01) and 5 × 5 km in Domain 2 (d02).
Figure 1. The two nested domains of the WRF simulations, with horizontal grid spacings of 25 × 25 km in Domain 1 (d01) and 5 × 5 km in Domain 2 (d02).
Atmosphere 15 00442 g001
Figure 2. Maximum 2 m temperatures (blue) on 16 May 2020 and 95th percentile of maximum daily 2 m temperatures in May (ref. period: 1975–2005) (red). ATH, THE and HER stand for Athens, Thessaloniki and Heraklion, respectively. “nan” stands for “not available”.
Figure 2. Maximum 2 m temperatures (blue) on 16 May 2020 and 95th percentile of maximum daily 2 m temperatures in May (ref. period: 1975–2005) (red). ATH, THE and HER stand for Athens, Thessaloniki and Heraklion, respectively. “nan” stands for “not available”.
Atmosphere 15 00442 g002
Figure 3. Anomalies (°C; contours) and standardized anomalies (standard deviations; colour shading) in the 850 hPa temperature (ERA5) at 12 UTC 16 May 2020, relative to its climatology for May 1975–2005 at 12 UTC (after [90]).
Figure 3. Anomalies (°C; contours) and standardized anomalies (standard deviations; colour shading) in the 850 hPa temperature (ERA5) at 12 UTC 16 May 2020, relative to its climatology for May 1975–2005 at 12 UTC (after [90]).
Atmosphere 15 00442 g003
Figure 4. The 72 h airmass backward trajectories of the warm airmass above mainland Greece on 16 May 2020 at 12 UTC. Terminal point (star sign in the map) coordinates: 39.64° N, 22.42° E; altitude: 1500 m above mean sea level (data: CFS reanalysis) (https://www.arl.noaa.gov/hysplit, accessed on 25 February 2024). Map boundaries: 24° N–42° N, 2° E–30° E. The latitudes in ° N (longitudes in ° E) are displayed within the graph at 17.5° E (37.5° N).
Figure 4. The 72 h airmass backward trajectories of the warm airmass above mainland Greece on 16 May 2020 at 12 UTC. Terminal point (star sign in the map) coordinates: 39.64° N, 22.42° E; altitude: 1500 m above mean sea level (data: CFS reanalysis) (https://www.arl.noaa.gov/hysplit, accessed on 25 February 2024). Map boundaries: 24° N–42° N, 2° E–30° E. The latitudes in ° N (longitudes in ° E) are displayed within the graph at 17.5° E (37.5° N).
Atmosphere 15 00442 g004
Figure 5. Temperature (°C—colour shading) at 850 hPa and geopotential height (gpm—contours) at 500 hPa at 12 UTC on 16 May 2020 (ERA5).
Figure 5. Temperature (°C—colour shading) at 850 hPa and geopotential height (gpm—contours) at 500 hPa at 12 UTC on 16 May 2020 (ERA5).
Atmosphere 15 00442 g005
Figure 6. Anomalies (gpm; contours) and standardized anomalies (standard deviations; colour shading) in the (a) 500 hPa and (b) 850 hPa geopotential heights (ERA5) at 12 UTC on 16 May 2020, relative to the climatology for May 1975–2005 at 12 UTC.
Figure 6. Anomalies (gpm; contours) and standardized anomalies (standard deviations; colour shading) in the (a) 500 hPa and (b) 850 hPa geopotential heights (ERA5) at 12 UTC on 16 May 2020, relative to the climatology for May 1975–2005 at 12 UTC.
Atmosphere 15 00442 g006
Figure 7. (a) Wind speed (m/s; colour shading) at 300 hPa. (b) Cross-section at 23° E of omega (Pa/s—colour shading) and zonal wind speed (stronger than 10 m/s—contours) (ERA5). Green lines mark the minimum and maximum latitudes of the study region.
Figure 7. (a) Wind speed (m/s; colour shading) at 300 hPa. (b) Cross-section at 23° E of omega (Pa/s—colour shading) and zonal wind speed (stronger than 10 m/s—contours) (ERA5). Green lines mark the minimum and maximum latitudes of the study region.
Atmosphere 15 00442 g007
Figure 8. Anomalies ((a) m3 and (b) °C; contours) and standardized anomalies (standard deviations; colour shading) in the (a) soil moisture content (m3) in the upper soil layer (0–7 cm depth) and (b) SST (°C) at 12 UTC on 16 May 2020, relative to their climatology for May 1975–2005 at 12 UTC.
Figure 8. Anomalies ((a) m3 and (b) °C; contours) and standardized anomalies (standard deviations; colour shading) in the (a) soil moisture content (m3) in the upper soil layer (0–7 cm depth) and (b) SST (°C) at 12 UTC on 16 May 2020, relative to their climatology for May 1975–2005 at 12 UTC.
Atmosphere 15 00442 g008
Figure 9. Boxplots of the 850 hPa temperatures (°C) predicted by WRF (25 km and 5 km grid spacing) and the six global S2S forecasts for the initialization dates of 9–15 April 2020 (week 1), 16–22 April 2020 (week 2) and 23–29 April 2020 (week 3). The validation period was 8–23 May 2020 at 00 UTC of each day. All the forecasts (global models and WRF) and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Figure 9. Boxplots of the 850 hPa temperatures (°C) predicted by WRF (25 km and 5 km grid spacing) and the six global S2S forecasts for the initialization dates of 9–15 April 2020 (week 1), 16–22 April 2020 (week 2) and 23–29 April 2020 (week 3). The validation period was 8–23 May 2020 at 00 UTC of each day. All the forecasts (global models and WRF) and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Atmosphere 15 00442 g009
Figure 10. Boxplots of the predicted daily 2 m maximum temperatures (°C) of WRF (25 km and 5 km grid spacing) for the initialization dates of 9–15 April 2020 (week 1), 16–22 April 2020 (week 2) and 23–29 April 2020 (week 3). All the WRF forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 0.25° × 0.25°.
Figure 10. Boxplots of the predicted daily 2 m maximum temperatures (°C) of WRF (25 km and 5 km grid spacing) for the initialization dates of 9–15 April 2020 (week 1), 16–22 April 2020 (week 2) and 23–29 April 2020 (week 3). All the WRF forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 0.25° × 0.25°.
Atmosphere 15 00442 g010
Figure 11. Quantile–Quantile plots of the distribution of (a) predicted temperatures at 00 UTC from all the models at 850 hPa and (b) maximum daily temperatures at 2 m from WRF 25 km and 5 km against the distribution of reanalysis (ERA5) 850 hPa temperature and maximum daily 2 m temperature, respectively, from week 3 (23–29 April). The mean squared error (MSE in °C2) values are presented in the legend. All the forecasts and the ERA5 reanalyses were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of (a) 1.5° × 1.5° and (b) 0.25° × 0.25°.
Figure 11. Quantile–Quantile plots of the distribution of (a) predicted temperatures at 00 UTC from all the models at 850 hPa and (b) maximum daily temperatures at 2 m from WRF 25 km and 5 km against the distribution of reanalysis (ERA5) 850 hPa temperature and maximum daily 2 m temperature, respectively, from week 3 (23–29 April). The mean squared error (MSE in °C2) values are presented in the legend. All the forecasts and the ERA5 reanalyses were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of (a) 1.5° × 1.5° and (b) 0.25° × 0.25°.
Atmosphere 15 00442 g011
Figure 12. Hit rates (solid lines) and false alarm rates (dashed lines) of 850 hPa temperatures predicted by WRF (25 km and 5 km grid spacing) and the six meteorological centres for the initialization dates of 9–15 (week 1), 16–22 (week 2) and 23–29 (week 3) April 2020. Threshold used: 95th percentile of the 850 hPa temperature of ERA5 (1975–2005). The validation period was 8–23 May 2020 at 00 UTC of each day. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Figure 12. Hit rates (solid lines) and false alarm rates (dashed lines) of 850 hPa temperatures predicted by WRF (25 km and 5 km grid spacing) and the six meteorological centres for the initialization dates of 9–15 (week 1), 16–22 (week 2) and 23–29 (week 3) April 2020. Threshold used: 95th percentile of the 850 hPa temperature of ERA5 (1975–2005). The validation period was 8–23 May 2020 at 00 UTC of each day. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Atmosphere 15 00442 g012
Figure 13. (a) ROC diagram of 850 hPa temperatures predicted by WRF (25 km and 5 km grid spacing) and the six meteorological centres for initialization dates of 23–29 April 2020. (b) Timeseries of the AUC scores across the three initialization periods (9–15, 16–22 and 23–29 April 2020). Threshold used: 95th percentile of the 850 hPa temperature of ERA5. The validation period was 8–23 May 2020 at 00 UTC of each day. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Figure 13. (a) ROC diagram of 850 hPa temperatures predicted by WRF (25 km and 5 km grid spacing) and the six meteorological centres for initialization dates of 23–29 April 2020. (b) Timeseries of the AUC scores across the three initialization periods (9–15, 16–22 and 23–29 April 2020). Threshold used: 95th percentile of the 850 hPa temperature of ERA5. The validation period was 8–23 May 2020 at 00 UTC of each day. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Atmosphere 15 00442 g013
Figure 14. (a) ROC diagram for the initialization period of 23–29 April 2020 (week 3) and (b) AUC scores for all three initialization periods (9–15, 16–22 and 23–29 April 2020) of the predicted maximum daily 2 m temperatures (T2max) for WRF (25 km and 5 km grid spacings). Threshold used: 95th percentile of the maximum daily 2 m temperature of ERA5. The validation period was 8–23 May 2020. All WRF forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 0.25° × 0.25°.
Figure 14. (a) ROC diagram for the initialization period of 23–29 April 2020 (week 3) and (b) AUC scores for all three initialization periods (9–15, 16–22 and 23–29 April 2020) of the predicted maximum daily 2 m temperatures (T2max) for WRF (25 km and 5 km grid spacings). Threshold used: 95th percentile of the maximum daily 2 m temperature of ERA5. The validation period was 8–23 May 2020. All WRF forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 0.25° × 0.25°.
Atmosphere 15 00442 g014
Figure 15. Reliability diagrams of predicted temperatures at 850 hPa of WRF (25 km and 5 km grid spacings) and the six meteorological centres for initialization dates of 23–29 April 2020. Threshold used: 95th percentile of the 850 hPa temperature of ERA5. The evaluation period was 8–23 May 2020 at 00 UTC of each day. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Figure 15. Reliability diagrams of predicted temperatures at 850 hPa of WRF (25 km and 5 km grid spacings) and the six meteorological centres for initialization dates of 23–29 April 2020. Threshold used: 95th percentile of the 850 hPa temperature of ERA5. The evaluation period was 8–23 May 2020 at 00 UTC of each day. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Atmosphere 15 00442 g015
Figure 16. Reliability diagram of the predicted maximum daily 2 m temperatures from WRF 25 km and 5 km for initialization dates of 23–29 April 2020. Threshold used: 95th percentile of the maximum daily 2 m temperature of ERA5. The evaluation period was 8–23 May 2020. All the WRF forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 0.25° × 0.25°.
Figure 16. Reliability diagram of the predicted maximum daily 2 m temperatures from WRF 25 km and 5 km for initialization dates of 23–29 April 2020. Threshold used: 95th percentile of the maximum daily 2 m temperature of ERA5. The evaluation period was 8–23 May 2020. All the WRF forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 0.25° × 0.25°.
Atmosphere 15 00442 g016
Figure 17. Extreme dependence index (EDI) vs. forecast probability of temperature at 850 hPa of WRF (25 km and 5 km grid-spacing) and the six meteorological centres for initialization dates of 23–29 April 2020. Threshold used: 95th percentile of the 850 hPa temperature of ERA5. The evaluation period was 8–23 May 2020 at 00 UTC of each day. A perfect score is 1 and no prediction ability is 0. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Figure 17. Extreme dependence index (EDI) vs. forecast probability of temperature at 850 hPa of WRF (25 km and 5 km grid-spacing) and the six meteorological centres for initialization dates of 23–29 April 2020. Threshold used: 95th percentile of the 850 hPa temperature of ERA5. The evaluation period was 8–23 May 2020 at 00 UTC of each day. A perfect score is 1 and no prediction ability is 0. All the forecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Atmosphere 15 00442 g017
Figure 18. Anomalies (gpm) of the ERA5 mean geopotential height for 8–23 May 2020 at 00 UTC (relative to the ERA5 mean geopotential height of the same period during 1975–2005) at (a) 850 hPa and (b) 500 hPa. The ERA5 reanalysis was interpolated to a grid spacing of 1.5° × 1.5°.
Figure 18. Anomalies (gpm) of the ERA5 mean geopotential height for 8–23 May 2020 at 00 UTC (relative to the ERA5 mean geopotential height of the same period during 1975–2005) at (a) 850 hPa and (b) 500 hPa. The ERA5 reanalysis was interpolated to a grid spacing of 1.5° × 1.5°.
Atmosphere 15 00442 g018
Figure 19. Mean geopotential height anomalies (gpm) at 850 hPa (1st column) and 500 hPa (2nd column) during 8–23 May 2020 at 00 UTC (relative to the ERA5 mean geopotential height of the same period during 1975–2005) from the forecasts of ECMWF, UKMO, CMA, KMA, NCEP, HMCR and WRF 25 km that were initialized in week 3 (23–29 April 2020). All the forecasts were interpolated to the same grid (20° Ν–70° N, 30° W–55° E) with a grid spacing of 1.5° × 1.5°. Data void regions (dark blue colour) appear in the locations where the 850 hPa forecasts of WRF 25 km crossed the model topography.
Figure 19. Mean geopotential height anomalies (gpm) at 850 hPa (1st column) and 500 hPa (2nd column) during 8–23 May 2020 at 00 UTC (relative to the ERA5 mean geopotential height of the same period during 1975–2005) from the forecasts of ECMWF, UKMO, CMA, KMA, NCEP, HMCR and WRF 25 km that were initialized in week 3 (23–29 April 2020). All the forecasts were interpolated to the same grid (20° Ν–70° N, 30° W–55° E) with a grid spacing of 1.5° × 1.5°. Data void regions (dark blue colour) appear in the locations where the 850 hPa forecasts of WRF 25 km crossed the model topography.
Atmosphere 15 00442 g019aAtmosphere 15 00442 g019b
Table 1. The datasets used in this study and their characteristics. The temporal resolution of the WRF forecasts and the existing S2S forecasts refers to the initialization of each forecast. The CMA reforecasts for 20 April were not available in the S2S database.
Table 1. The datasets used in this study and their characteristics. The temporal resolution of the WRF forecasts and the existing S2S forecasts refers to the initialization of each forecast. The CMA reforecasts for 20 April were not available in the S2S database.
NameCentre(s)Grid SpacingTemporal Resolution
WRF
Forecasts
WRFAUTH * forecasts forced by NCEP CFSv.2
  • 25 km (Europe, N. Africa, N. Atlantic)
  • 5 km (Greece)
12 h (00 and 12 UTC)
S2S
Forecasts and
Reforecasts
S2S
Database
ECMWF *, UKMO, CMA, KMA, NCEP, HMCRDifferent native
resolutions—archived at 1.5° (S2S database)
  • 2/week (ECMWF, CMA)
  • 1/week (HMCR)
  • 24 h (UKMO ***, KMA, NCEP)
ReanalysisERA5
ERA5-Land
ECMWF0.25°
0.1°
Hourly
Climatic Parameters and
Heat Indices
Weather
Stations
NOA ** and NCDCPinpoint observationsDaily
* Aristotle University of Thessaloniki (AUTH), European Centre for Medium-Range Weather Forecasts (ECMWF), United Kingdom Met Office (UKMO), China Meteorological Administration (CMA), Korean Meteorological Administration (KMA), National Centers for Environmental Prediction (NCEP), Hydrometeorological Research Centre of Russia (HMCR). ** National Observatory of Athens, Greece (NOA); National Climatic Data Center (NCDC); *** UKMO and KMA reforecast frequency: four/month (1st, 9th, 17th and 25th day of each month).
Table 2. Global S2S model characteristics of real-time forecasts and reforecasts used in this study.
Table 2. Global S2S model characteristics of real-time forecasts and reforecasts used in this study.
Meteorological CentreMaximum Lead TimeNative Spatial ResolutionModel Top LevelReforecast PeriodReforecast Ensemble Size
ECMWF46 days32 km~0.29°0.01 hPaPast 20 years11
UKMO60 days60 km~0.54°85 km1993–20167
CMA60 days90 km~0.8°0.1 hPaPast 15 years4
KMA60 days60 km~0.54°85 km1991–20103
NCEP44 days100 km~0.9°0.02 hPa1999–20104
HMCR44 days100 km~0.9°5 hPa1985–201010
Table 3. WRF configuration characteristics and its parameters.
Table 3. WRF configuration characteristics and its parameters.
WRF
Input DataCFS (NCEP) 1° × 1°
Domain 1 (Europe–N. Africa–N. Atlantic)25 km × 25 km
Domain 2 (Greece)5 km × 5 km
Model Top Level2 hPa~39 km
Vertical Levels (Surface to model top)53
Frequency of Output DataThree hours
Update Frequency of Sea Surface TemperaturesSix hours from NCEP S2S forecast data
Forecast HorizonUntil 31 May 2020
Physical Parameters
MicrophysicsThompson [75] (option 8 of WRF)
Cumulus ConvectionGrell–Freitas [76] (option 3)
Shortwave/Longwave RadiationRapid Radiative Transfer Model for GCMs (RRTMG) [77] (option 4)
Planetary Boundary LayerMellor–Yamada–Janjic [78] (option 2)
Land SurfaceNoah Land Surface Model [79] (option 2)
Surface LayerEta [78,80] (option 2)
Table 4. The days from 9 May to 22 May 2020 with EHF > 0 °C2, THI > 70 °F, UTCI > 32 °C and T2max > 95th percentile (T2max95) (highlighted in orange) in the different datasets (stations, ERA5 and ERA5-land) for the three stations of Athens, Thessaloniki and Heraklion. Green blocks indicate the days when the EHF, THI, UTCI and T2max95 thresholds were not exceeded. Red blocks indicate the date when the maximum value of each index occurred in each city, along with its corresponding maximum value.
Table 4. The days from 9 May to 22 May 2020 with EHF > 0 °C2, THI > 70 °F, UTCI > 32 °C and T2max > 95th percentile (T2max95) (highlighted in orange) in the different datasets (stations, ERA5 and ERA5-land) for the three stations of Athens, Thessaloniki and Heraklion. Green blocks indicate the days when the EHF, THI, UTCI and T2max95 thresholds were not exceeded. Red blocks indicate the date when the maximum value of each index occurred in each city, along with its corresponding maximum value.
Athens *910111213141516171819202122
EHF (Station) > 0 °C2 136.7 °C2
EHF (ERA5) > 0 °C2 114.8 °C2
EHF (ERA5-Land) > 0 °C2 99 °C2
THI (ERA5) > 70 °F 75.4 °F
UTCI (ERA5) > 32 °C 41.7 °C
T2max > T2max95 (Station)Above Thres 37.8 °C
T2max > T2max95 (ERA5)Below Thres 39.7 °C
T2max > T2max95 (ERA5-L.)Max. Value 34.5 °C
Thessaloniki *910111213141516171819202122
EHF (Station) > 0 °C2 62.7 °C2
EHF (ERA5) > 0 °C2 62.3 °C2
EHF (ERA5-Land) > 0 °C2 64.1 °C2
THI (ERA5) > 70 °F 78.4 °F
UTCI (ERA5) > 32 °C 36.6 °C
T2max > T2max95 (Station) 36 °C
T2max > T2max95 (ERA5) 33.1 °C
T2max > T2max95 (ERA5-L.) 31.8 °C
Heraklion *910111213141516171819202122
EHF (Station) > 0 °C2 126.3 °C2
EHF (ERA5) > 0 °C2 146 °C2
EHF (ERA5-Land) > 0 °C2 147 °C2
THI (ERA5) > 70 °F 76.8 °F
UTCI (ERA5) > 32 °C 42 °C
T2max > T2max95 (Station) 38.5 °C
T2max > T2max95 (ERA5) 39.2 °C
T2max > T2max95 (ERA5-L.) 38.1 °C
* Athens Eleftherios Venizelos airport: 37.90° N 23.77° E; Thessaloniki Makedonia airport: 40.54° N 22.97° E; Heraklion Nikos Kazantzakis airport: 35.34° N 25.17° E.
Table 5. Uncentred Anomaly Correlation Coefficients (ACCs) for initialization dates from 9 April to 29 April 2020 (first row) from different global S2S forecasts (perturbed and unperturbed forecasts). The results of the perturbed forecasts are ACCs averaged over all the corresponding ensemble members. The ACC of each individual perturbed member was calculated before averaging. Evaluation period: 8–23 May 2020 at 00 UTC. ACC ≤ 0 (red), ACC > 0 (colour shading). All the global S2S forecasts, reforecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
Table 5. Uncentred Anomaly Correlation Coefficients (ACCs) for initialization dates from 9 April to 29 April 2020 (first row) from different global S2S forecasts (perturbed and unperturbed forecasts). The results of the perturbed forecasts are ACCs averaged over all the corresponding ensemble members. The ACC of each individual perturbed member was calculated before averaging. Evaluation period: 8–23 May 2020 at 00 UTC. ACC ≤ 0 (red), ACC > 0 (colour shading). All the global S2S forecasts, reforecasts and the ERA5 reanalysis were interpolated to the same grid (34.5° Ν–42° N, 19.5° E–29.5° E) with a grid spacing of 1.5° × 1.5°.
91011121314151617181920212223242526272829
ECMWF Unperturbed0.39 −0.15 0.44 −0.51 0.49 0.70
ECMWF Perturbed0.23 0.53 −0.30 0.72 0.69 0.78
UKMO Unperturbed0.33 0.30 0.84
UKMO Perturbed0.10 0.10 0.59
CMA Unperturbed0.78 −0.21 −0.85 −0.86 −0.88
CMA Perturbed−0.25 −0.80 −0.10 0.00 −0.30
KMA Unperturbed−0.71 −0.21 0.72
KMA Perturbed0.49 −0.30 0.32
NCEP Unperturbed0.40−0.480.300.56−0.32−0.450.320.57−0.510.31−0.66−0.40−0.17−0.10−0.400.10−0.260.820.810.320.14
NCEP Perturbed−0.210.40−0.480.330.830.52−0.850.23−0.760.800.32−0.21−0.75−0.21−0.040.580.810.830.600.700.68
HMCR Unperturbed−0.58 0.00 0.14
HMCR Perturbed−0.30 0.09 −0.17
Atmosphere 15 00442 i001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mitropoulos, D.; Pytharoulis, I.; Zanis, P.; Anagnostopoulou, C. Synoptic Analysis and Subseasonal Predictability of an Early Heatwave in the Eastern Mediterranean. Atmosphere 2024, 15, 442. https://doi.org/10.3390/atmos15040442

AMA Style

Mitropoulos D, Pytharoulis I, Zanis P, Anagnostopoulou C. Synoptic Analysis and Subseasonal Predictability of an Early Heatwave in the Eastern Mediterranean. Atmosphere. 2024; 15(4):442. https://doi.org/10.3390/atmos15040442

Chicago/Turabian Style

Mitropoulos, Dimitris, Ioannis Pytharoulis, Prodromos Zanis, and Christina Anagnostopoulou. 2024. "Synoptic Analysis and Subseasonal Predictability of an Early Heatwave in the Eastern Mediterranean" Atmosphere 15, no. 4: 442. https://doi.org/10.3390/atmos15040442

APA Style

Mitropoulos, D., Pytharoulis, I., Zanis, P., & Anagnostopoulou, C. (2024). Synoptic Analysis and Subseasonal Predictability of an Early Heatwave in the Eastern Mediterranean. Atmosphere, 15(4), 442. https://doi.org/10.3390/atmos15040442

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop