The Role of Sea Surface Temperature Forcing in the Life-Cycle of Mediterranean Cyclones

: Air–sea interface processes are highly associated with the evolution and intensity of marine-developed storms. Specifically, in the Mediterranean Sea, the air–ocean temperature deviations have a profound role during the several stages of Mediterranean cyclonic events. Subsequently, this enhances the need for better knowledge and representation of the sea surface temperature (SST). In this work, an analysis of the impact and uncertainty of the SST from different well-known datasets on the life-cycle of Mediterranean cyclones is attempted. Daily SST from the Real Time Global SST (RTG_SST) and hourly SST fields from the Operational SST and Sea Ice Ocean Analysis (OSTIA) and the NEMO ocean circulation model are implemented in the RAMS/ICLAMS-WAM coupled modeling system. For the needs of the study, the Mediterranean cyclones Trixi, Numa, and Zorbas were selected. Numerical experiments covered all stages of their life-cycles (five to seven days). Model results have been analyzed in terms of storm tracks and intensities, cyclonic structural characteristics, and derived heat fluxes. Remote sensing data from the Integrated Multi-satellitE Retrievals (IMERG) for Global Precipitation Measurements (GPM), Blended Sea Winds, and JASON altimetry missions were employed for a qualitative and quantitative comparison of modeled results in precipitation, maximum surface wind speed, and wave height. Spatiotemporal deviations in the SST forcing rather than significant differences in the maximum/minimum SST values, seem to mainly contribute to the differences between the model results. Considerable deviations emerged in the resulting heat fluxes, while the most important differences were found in precipitation exhibiting spatial and intensity variations reaching 100 mm. The employment of widely used products is shown to result in different outcomes and this point should be taken into consideration in forecasting and early warning systems.


Introduction
Air-ocean interface processes are highly associated with the development, evolution, and intensity of extreme weather phenomena. Heat, moisture, and momentum exchanges due to ocean drag and thermodynamic disequilibrium between the ocean surface and the upper air result in intense atmospheric perturbations.
The Mediterranean Sea is an area where the aforementioned air-sea processes very often contribute to intense cyclonic activity. The presence of cold cut-off lows in the middle and upper troposphere blended with the warmer water surfaces leads to cyclogenesis. However, despite the considerable number of cyclones over the Mediterranean basin, only a few of them on an annual basis demonstrate characteristics similar to tropical cyclones (TCs), known as Mediterranean tropical-like cyclones (TLCs) or medicanes. Some important features characterizing a Mediterranean TLC in the literature are a rounded shape with a cloudless core (cyclone eye), a drop in sea level pressure, a warm core structure at the mid-troposphere, heavy rainfall, and strong cyclonic winds [1][2][3].
Due to their formation in a marine environment, the air-ocean temperature deviations are one of the main triggering mechanisms of Mediterranean TLCs, equivalent to the tropical hurricanes [4][5][6]. The intrusions of cold air from the European mainland over the warmer sea surface enhance the role of the heat and moisture fluxes in the development and evolution of such systems.
The better knowledge and representation of sea surface temperatures (SSTs) during the several stages and the intensity of marine-developed storms have received significant attention in various associated studies. TCs have an impact on the upper-ocean layer causing alterations in the water surface temperature. Wind mixing and upwelling of colder water masses from deeper locations are some of the underlying processes taking place under TCs. This storm-induced SST cooling is expected to affect the air-sea enthalpy differences and the intensity of TCs [4], respectively.
An important limitation in numerical studies focusing on the sensitivity of the SST lies in the absence of large spatially distributed observational datasets. The sparsity of in-situ and satellite observations under extreme weather conditions [7] may lead to uncertainties both in the proper representation of the air-sea heat exchanges as well as the evaluation procedures. For example, in [4], a homogenous decrease of 3 °C in the SST fields during a Mediterranean cyclone led to a reduction in sea surface fluxes up to 150 W/m² and an elimination of the tropical characteristics. Moreover, in [8], the use of climatological SST in a modeling study of a Mediterranean cyclone decreased its lifetime. Towards the better description of the prevailing SST conditions, large data networks have been established in the last decade. This led to the construction of high-resolution SST gridded fields with the implementation of quality-controlled satellite, buoys and other ocean measurements [9,10].
Concurrently, the lack of information regarding the fast-evolving conditions in the air-water interface may conceal several features of these extreme events. To efficiently represent the on-going interplay between the atmospheric and oceanic environments, coupling methodologies have been developed [11]. Introducing the ocean effects in the operational European Centre Weather for Medium Range Weather Forecasting (ECMWF) high-resolution forecast improved the foresting intensity of hurricanes [12]. Likewise, the inclusion of air-ocean interaction methods in the Coupled Ocean/Atmosphere Mesoscale Prediction System for Tropical Cyclones (COAMPS-TC), enhanced the foresting capability in the track, the intensity, and the fine-scale structure in a number of hurricane cases [13].
In this context, this study aims at the understanding of how model results are affected in the simulation of Mediterranean cyclones using SSTs from different well-known datasets. More precisely, the current work examines the sensitivity of the atmospheric-ocean surface conditions on the formation and the evolution of Mediterranean TLCs. A coupled modeling system is used, consisting of atmospheric, wave, and ocean components. To consider the evolution of ocean temperature throughout the several stages of the cyclonic events, modeled and observational analysis SST fields enriched with various sources of remote sensing records are tested.
The manuscript is structured as follows. The modeling system used, the methodology followed, the description of the methodology, experimental cases, and the data used are described in Section 2. A brief analysis of the different input data used, the results of the performed experiments, and evaluation are discussed in Section 3. Further discussion and conclusions are presented in Section 4.

Atmospheric, Wave, and Ocean Components
A coupled modeling system was used to provide continuous feedback of information between the atmospheric-wave and ocean environments. The modeling system included an atmospheric and a wave model, online coupled. Moreover, the boundary conditions in the atmospheric model were continuously updated by a 2-D ocean component utilizing SST fields from several sources.

Atmospheric Component
The atmospheric component was an enhanced version of the Regional Atmospheric Modeling System (RAMS), the Integrated Community Limited Area Modeling System RAMS/ICLAMS [14][15][16][17]. Among its most notable characteristics is the online treatment of mineral dust and sea salt from wave breaking [18]. These natural aerosols contribute to model calculations through feedback mechanisms, including direct, semi-direct, and indirect effects in the radiation scheme [19] and the ice nuclei (IN) and cloud condensation nuclei (CCN) estimations [20].

Wave Component
The wave model used was the Wave Analysis Model (WAM) [21] version CY33R1 [22]. The model simulates the distribution of wave variance in different frequencies and propagation directions. The basic transport equation describes the evolution of the two-dimensional wave spectrum employed by the model. The solution of this equation leads to the calculation of different parameters, such as significant and swell wave height, peak frequency, and directional spread. The used version utilizes explicit source functions for the description of white-capping dissipation and bottom friction. Additional features, such as the calculation of depth induced wave breaking and shallow water effects, were also implemented.

Ocean Component
To consider the dynamic variation of the prevailing SST conditions, a 2-D model was constructed. This concerns an algorithm that inputs the SST from various sources and in the next step, interpolates it in a gridded domain corresponding to the spatial coverage of the atmospheric model. Unavoidably a small distortion of the initial SST fields is expected mainly in nearshore and complex areas due to the interpolation process.

Coupled Modeling System
For the model components to be able to operate synchronously in a coupled way, the OASIS-MCT [23] coupling module was used. The last allowed the exchange of parameters between the atmospheric-wave components in a two-way mode and the atmospheric-ocean in a one-way mode (Figure 1). The communication between the different components was performed in desired time intervals depending on its model time step. Surface wind speed components and air density were passed to the wave model, with wind speed to act as a forcing factor and density to be used in the computation of ocean surface stress. Simultaneously, the ocean surface roughness was transferred by the wave model to the atmospheric one. This allows a more realistic representation of the sea state conditions in the atmospheric model, affecting both the wind speed profile and the calculation of the surface heat fluxes. Precisely, the ocean surface roughness zo was derived in the atmospheric model parameterized, employing wave characteristics, such as wave significant height (hs) and wave slope (hs/Lp), following the Taylor and Yelland parameterization [24]: zo = 1200 · hs · (hs/Lp) 4.5 (1) with Lp the wavelength of the peak frequency. This relationship has been tested in a number of datasets and sea conditions, such as wind-sea, swell, and short/long fetch conditions [25][26][27].

Methodology and Data Used
For the needs of the study, three Mediterranean cyclones were simulated employing the coupled system and performing sensitivity tests with different SST fields. The analysis initially focused on the spatiotemporal evolution of the different ocean surface temperatures in terms of anomalies and comparison with data. Concerning the numerical experiments, an examination and comparison of storm tracks, maximum surface wind speed, and air-sea fluxes were performed, supported with available measurements. Furthermore, a phase analysis algorithm was applied for cyclone structure and characteristics based on the phase diagrams proposed by Hart [28]. A detailed description is given in Appendix A. Finally, the accumulated precipitation obtained from the different approaches was studied and assessed with the aid of remote sensing estimates.

Models Setup and Configuration
Concerning the atmospheric component, two nests were used, with the fine one to cover the entire Mediterranean and most parts of Europe and North Africa. This setup was followed to be able to capture the effects of the differential heating between Africa and Europe [29], the possible direct, indirect and semi-direct effects of dust particles and to include the whole Mediterranean Basin. The outer one was much larger to keep the lateral boundary conditions far from the area of interest. Subsequently, the wave model was running with a single domain covering the entire Mediterranean Sea ( Figure 2). This configuration was expected to minimize the lateral boundary conditions problems and was computationally affordable.  Table 1. The SST for the atmospheric model stem from three different sources briefly described below.

RTG-SST
The Real Time Global SST (RTG_SST) is a global operational analysis product available since 2001 [30]. An analysis algorithm operates once per day combining in-situ and satellite data of the last 24 h. In-situ data come from fixed buoys, drifting buoys, and ships. Satellite retrievals come from NOAA and METOP-A AVHRR data that are averaged in day and night fields following a physical algorithm. Currently, the RTG_SST data are provided on a daily basis on a resolution of 1/12 × 1/12 degrees.

OSTIA-SST
The Operational SST and Sea Ice Ocean Analysis (OSTIA), is a global coverage SST dataset produced at the United Kingdom Met Office [10] with a resolution of 1/4 × 1/4 degrees. The system blends in-situ data and satellite observations, a warm layer diurnal algorithm that assimilates the satellite data and a cool skin model. Remote sensing data of different coverage stem from eleven different sensors of Infrared and Microwave sensor type. Satellite retrievals come from SEVIRI, GOES-W, MTSAT 2, and METOP-A AVHRR data [31]. OSTIA SST analysis fields are used by the United Kingdom Met Office and European Centre Weather for Medium Range Weather Forecasting (ECMWF) models for operational forecasts.

NEMO-SST
Forecast SST data are available through the Operational Mercator global ocean analysis and forecast system. The system is based on the Nucleus for European Modeling of the Ocean (NEMO) ocean model, version 3.1 [32], with a horizontal resolution of 1/12 × 1/12 degrees. Starting from the sea surface, fifty vertical levels are used up to 5000 m depth, with the first 28 distributing within the first 450 m water depth. Bathymetry data come from the GEBCO8 [33] database for depths less than 300 m and from ETOPO1 [34] for deeper waters. Model physics and recent updates and improvements concerning the advection, assimilation, and sea/ice schemes, as well as evaluation studies, can be found in [35] and [9].

Case Studies
Three recent Mediterranean cyclones with tropical-like characteristics were analyzed. The storm cases were selected from the EUMESTAT database, identified by an operating network of satellites. The events were generated and grew in a mature stage over the Central and Eastern Mediterranean waters. These areas belong to a zone where such systems usually appear during late autumn and early winter due to relatively high SST [36]. The selection was based on the availability of the OSTIA SST and the data for evaluation.
The selected cases and periods are summarized in Table 2. For the convenience of the reader, the performed simulations throughout the manuscript are stated with the names of the origin of the SST as RTG-sst, NEMO-sst, and OSTIA-sst. The first Mediterranean storm took place during the period, 27 October to 1 November 2016. It was formed during the first hours of 27 October in the offshore area along the western shoreline of Italy. Moving southward, the system was located over relatively warm waters in the area east of Malta (Figure 3a), and from the first hours of 30 October (Figure 3b), an eye-like structure was developed and started moving eastward towards Greece. Afterward, the storm crossed the South Aegean Sea and rapidly lost its strength as it reached landfall on the southwest coast of Turkey during 1 November. The second Mediterranean cyclone under consideration occurred during November 2017. The cyclone was formed after the dissipation of an extra-tropical low that was located in the British Isles. The storm moved southeastward and passed France and Italy on 15 November. The last hours of the same day, the storm started obtaining sub-tropical characteristics. During 17 November, when it was located in the Ionian Sea (Figure 4a), the storm presented tropical characteristics [37] that maintained until the next day (Figure 4b). Moving eastward on 19 November, it dissipated over the Greek Peninsula. The third selected case was related to a storm developed over warm waters in the eastern Mediterranean Sea on 27 September 2018. Heading northeastward, the following day, the cyclone gradually intensified and obtained tropical-like characteristics (Figure 5a). The cyclone reached its peak as it approached the southwest part of Greece on 29 September (Figure 5b). Continuing northeastward, the storm passed along the Aegean Sea and dissipated over Turkey during the next day.

Data Used
For evaluation purposes, remote sensing, in-situ, and model reanalysis data were used. Satellite-based winds and precipitation were implemented from the Blended Sea Winds and the Integrated Multi-satellitE Retrievals for GPM (IMERG), respectively. Remotely sensed altimeter data from JASON2 and JASON3 satellite missions were employed for wave height comparison. Moreover, SST records from a buoy in the South Aegean Sea were used. To have an alternative estimate of the derived heat fluxes and rainfall, Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) reanalysis data were also considered.
The Blended Sea Winds is a gridded and high-resolution dataset [38,39]. It contains surface wind speed and direction data on a global 0.25° grid in multiple time resolutions. Wind speed vectors are generated from multiple satellite observations stemming from the Defense Meteorological Satellites Program (DMSP) SSM/I, the QuikSCAT, the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) and the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI). A spatiotemporal interpolation is made to fill possible temporal and spatial gaps of the individual satellite samplings. At the same time, correction methods for the reduction in the subsampling aliases and the imported random errors are used. Concerning the wind direction, two sources are employed depending on the products. For research delayed mode product, the source is the National Centers for Environmental Prediction (NCEP) Reanalysis 2, while for the near-real-time product, the source is the numerical weather prediction of ECMWF. Nevertheless, data gaps are not completely eliminated and while the representation of wind fields is sufficient, errors in both low and high wind speed categories were found [40,41]. In this work, the finest temporal resolution of six-hours was implemented to test the model behavior regarding the cyclonic formulation and the maximum wind speed recorded.
IMERG is a precipitation product based on the inter-calibration, merging, and interpolation of measurements from a network of satellites together with microwave-calibrated infrared (IR) satellite estimates, precipitation gauge analyses, and potentially other precipitation estimators. More specifically, it utilizes TRMM and GPM eras globally. It has a spatial resolution of 0.1° × 0.1° and a temporal of 30 min covering up to ±60° latitudes. The system runs several times providing different products depending on user requirements. The first is a quick estimate (IMERG "early run"), provided with a 6-h delay for real-time applications, such as hazard predictions. The second provides better estimates as more data are taken into consideration. The so-called IMERG "late run", has a delay of 18 h and could be employed in other applications, such as crop forecasting. The final step is the one also used in the present study (IMERG "final run"). It uses monthly gauge data to create research-level products and comes at a 4-month delay. More details can be found in the following articles and references within [42][43][44][45][46]. It should be mentioned that several studies [47][48][49] have indicated an overestimation of IMERG products in the estimated precipitation under heavy rainfall events.
JASON2 and JASON3 altimetry missions perform global measurements of sea surface height [50,51]. Both repeat 256 passes in each cycle with an approximately 10 day return period. Wave height, altimetry range and other relevant data are measured at C (5.3 GHz) and Ku (13.575 GHz) band frequencies. For the needs of the study, significant wave records from the Ku band frequency were used after a simple quality control procedure for the removal of erroneous data.
Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2), is a global atmospheric reanalysis with a resolution of 0.5° × 0.625° degrees [52]. Among its capabilities is the assimilation of satellite measurements and aerosol observations and new schemes for mass conservation and water balance in forecasting and assimilation procedures [53].

Results and Discussion
The study examined the uncertainty in model results using different SSTs in the simulation of Mediterranean cyclones. As a first step, a comparison between the time-series of the SST from the three selected datasets and an independent source that was data from a buoy with reference 61277 at a selected location was performed. The location corresponded to the coordinates (35.72° N, 25.13° E), , that was inside the affected zone of the selected events. Unfortunately, no other buoy covering all selected cases was available. Regarding the case of cyclone Trixi, SSTs from RTG and NEMO compared quite well, following the fluctuations of the observational data (Figure 6a). The SST from OSTIA tended to underestimate at this point, providing values around 21.0-21.5 °C, about one degree less than the other sources. The second case concerning cyclone Numa took place during November with lower temperatures as expected (Figure 6b). The NEMO and OSTIA products underestimated the values to a small degree, while RTG could be characterized by a slight overestimation. In any case, the first two datasets followed the trend of the measurements, while RTG did not. For the case of cyclone Zorbas, the sea was warmer than in the previous case since it was October (Figure 6c). NEMO and OSTIA datasets had a similar behavior as compared with the observations with values that gradually decreased. The RTG performance could be separated into two parts. During the first two days, RTG overestimated the water temperature, while for the rest of the duration of the event, it managed to better capture the SST. It should be noted that RTG is a daily product, something easily understood through the figures. The spatial variability of the SST alongside the different phases of the events under investigation is also discussed. Precisely, the differences between the SST in several stages of the simulations and the initial one were examined (Figure 7). For simplicity, the SST anomalies for the case of cyclone Zorbas are the ones illustrated. This event was also found to present the highest variations in the evolution of water surface temperatures between the various datasets. In all cases, sea surface tended to get colder as the system evolved. This is something expected as the cyclone gains energy from the warm water. Other processes contributing to SST reduction can be evaporative cooling processes, mixing of the top water layers due to wind and precipitation effect. It is well known that tropical cyclones can impact the upper ocean by causing SST cooling by as much as 6 °C [54]. In general, SST cooling can decrease the air-sea heat exchange and, therefore, reduce cyclone intensity [55][56][57]. Obviously, here, the cooling effects were smaller due to the scaling of the event and the local characteristics of the basin. Beginning with the first 24 h of the simulations, the largest differences were observed in the NEMO product with the SST dropping more than one degree in the area north of the cyclonic center. For the same time, the OSTIA dataset that is satellite oriented was characterized by a slight maximum drop around 0.5 °C. In contrast, the RTG daily differences seemed to be negligible. It should also be mentioned that the initial fields of RTG and NEMO showed a similar pattern and magnitude, while the OSTIA had slightly lower values in the initial fields. Continuing with the next 48 h, the maximum temperature drop of 2 °C was observed in NEMO-sst simulation, while both RTG and OSTIA had a maximum decrease of around 1 °C. Concurrently in RTG, this decrease was traced on the east side of the storm, while for the other two datasets, it was around the cyclone. Finally, the third day, 72 h from the beginning of the simulations, was characterized by even lower temperatures. The passage of the cyclone deepened the thermocline base, causing a further cooling of the ocean surface temperature. The maximum differences were also found in the NEMO-sst output, followed by RTG and OSTIA. Here the cyclone was in its dissipation phase, and the cooling was due to processes that took place earlier this day, such as wind mixing and upwelling. Model results in storm tracks and intensities based on surface pressure and particularly the mean sea level pressure (MSLP), are illustrated in Figure 8. Overall, small differences between the path and the translation speed appeared. This can also be explained due to the dependence of the cyclone tracks on large-scale processes. Spatially, divergences could be mainly traced in the initial stage of the storms, while differences during their lifetime emerged only in the 2016 storm ( Figure  8a). Concerning the evolution of MSLP in time, a consistency in the variation of the minimum values between the different numerical experiments was presented. Regarding the intensity, notable deviations appeared in cyclone Trixi and Zorbas cases. In cyclone Trixi, the numerical simulation of NEMO-sst provided the lowest MSLP values, while in cyclone Zorbas, the use of RTG SST distinctly increased the strength of the cyclone with a drop of 1 hPa as compared to the other two simulations. To understand the tropical development during the cyclonic evolution along with the level of the warm/cold core structure and symmetry, cyclone phase diagrams have been formed for each case. Precisely, the daily structural behavior for each cyclone is illustrated. Both parameters of thermal symmetry and thermal winds were calculated in a radius of 200 km around the location of MSLP for six-hour intervals.
Cyclone Trixi presented a symmetric formation for most of its lifetime, with a weak warm-core structure in the lower troposphere and cold-core in the upper, as displayed in Figure 9. A transition from extra-tropical to tropical behavior is illustrated for the 28 October at 12:00, where the minimum surface pressure appeared. The period before and after, the storm exhibited extra-tropical characteristics. Cyclone Numa also demonstrated extra-tropical type characteristics. The system underwent two different phases, as shown in Figure 10. In the first, starting from 16 November, a warm-core structure emerged in the lower atmospheric layers that lasted for three days. This period the cyclone also tended to be symmetric for less than twelve hours. During the second phase from 19 November, the cyclone obtained an asymmetric cold-core structure. Cyclone Zorbas rapidly evolved to a tropical cyclone twelve hours after its genesis on 27 September with a deep warm-core structure (Figure 11). A thermally direct circulation was evident in both the lower and upper troposphere, which lasted until 29 September at 18.00. Reaching the Greek Peninsula the next hours, the system weakened and obtained an asymmetric structure. Based on the phase diagram analysis, more intensified warm-core characteristics were illustrated in cyclone Trixi and Zorbas simulations using the RTG SST dataset (Figure 9b and Figure  11b). Apart from this point, no significant differences could be traced in the structural features from the implementation of different SSTs.
The coupling of ocean surface temperature in the atmospheric-wave environment can also be expressed in terms of air-sea fluxes. It should also be stated that air-sea temperature and humidity differences are not the only source of upward/downward energy transfer in the marine boundary layer. Wind stress and ocean surface roughness also have an impact on momentum, heat, and moisture fluxes. Both were also considered in the implemented atmospheric-wave/ocean modeling system. The maximum interfacial air-sea sensible and latent heat fluxes were estimated from the different SST used in the model simulations ( Figure 12). The heat fluxes from MERRA-2 reanalysis were also demonstrated to have an additional individual dataset for reference. Maximum values were based on the storm tracks of the surface pressure, where the sensible and latent heat maxima which emerged over water in a radius of 5° around the center of the cyclone were taken into account. Latent heat fluxes clearly overcame the sensible heat. The highest energy fluxes appeared in cyclone Zorbas's case due to warmer water temperatures prevailing in the area where the cyclone was generated. Overall, the intensity in fluxes was higher in RTG-sst simulations. This can also be related to the absence of SST cooling in RTG included in diurnal variability. Differences between the simulations were mainly traced in cyclone Trixi and secondly in cyclone Zorbas. In contrast, the resulting fluxes in cyclone Numa were characterized by smaller deviations between the different approaches. For cyclone Trixi, the peaks in the maximum sensible heat fluxes were almost the same for all the simulations. The highest values presented for the RTG-sst simulations exceeded the other two simulations in a range of 50 to 70 W/m². Concerning the maximum latent heat fluxes, again, the fluxes from RTG simulations in most of the cases were higher than those from the other two simulations in most of the cyclonic lifetimes. In this case, differences reached 150-200 W/m². The same tendency, but for smaller time scales, was illustrated in cyclone Zorbas. The heat fluxes stemming from RTG-sst simulations exceeded the ones from NEMO-sst and OSTIA-sst simulations. Precisely, the differences were up to 50 W/m² for the sensible and between 100 to 200 W/m² for the latent heat fluxes. For almost all cases, heat fluxes from MERRA-2 fluctuated in a similar manner to the ones from the performed simulations. Certainly, the variation of the maximum fluxes was traced in lower energy scales due to the coarser resolution of MERRA-2. Even though there were still some cases where heat fluxes values of MERRA-2 could be found close to the fluxes from the model simulations. To further analyze the differences between the applied methods, the temporal development of maximum surface wind speed was examined. In this case, wind speed magnitudes from blended sea winds were also used ( Figure 13). It should be noted that blended winds were employed as an indicator for the variation and tendency of the forecasted winds rather than a means for evaluation, due to the issues stated earlier in the description of the Blended Sea Winds product. Considering that wind speed and pressure are associated parameters through the pressure gradients, maximum wind speed also deviated similarly to the lowest pressure values for the different SST couplings. For all cases, wind speed compared well to the one presented by the Blended Sea Winds. Differences in the maximum values were evident in the range of 4 to 5 m/s and could be found in the cyclones Trixi and Zorbas. It should also be stated that maximum wind speed results in cyclone Numa are in good agreement with the ones presented in a similar study for the same event [37]. The resulting significant wave heights from the different runs were evaluated against the ones derived from JASON satellites. To perform the comparison, satellite data from the nearest hourly model time were employed. Figure 14 displays the cumulative scatter plots for each applied methodology covering the significant wave height (SWH) values from all the test cases among the results in statistical indicators of performance. Particularly, the bias (in m), the mean absolute error (MAE) (in m), the root mean square error (RMSE) and the Pierson correlation coefficient (Corr.Coef) were used [58]. For all cases, SWHs seemed to be gathered in the range of 0 to 2.5 m, while statistical results were almost identical. This was rather expected based on the aforementioned results in wind speed. Distances (errors) between observed and model values were increasing for higher SWHs with a tendency for over-prediction for waves over approximately 5 m. To examine the effects in precipitation, the spatial and temporal distributions of modeled and satellite estimated rainfalls from IMERG were compared. The purpose was to identify qualitatively whether the rain patterns and locations of the maxima between modeled and satellite estimates were in agreement. Therefore, the comparison concerned the accumulated precipitation fields in the modeling of each cyclone. Accumulated precipitations from MERRA-2 were also considered. To better compare the rain fields from the different simulations, RTG-sst simulations were used as a reference and subtracted from NEMO-sst and OSTIA-sst simulations.
For cyclone Trixi, both satellite and modeled accumulated rains were concentrated in the western region of Greece in the south Ionian Sea, as shown in Figure 15. The spatial differences between RTG-sst-OSTIA-sst and RTG-sst-NEMO-sst indicated that all modeled precipitation fields were positioned in the same area. The same pattern could be observed and for the MERRA-2 precipitation fields. A close agreement in the distribution, as well as the position of the maxima between RTG-sst modeled and IMERG satellite data, was evident. A notable underestimation with regard to the remote sensing rain data occurred in the peak values. This can be partially attributed to the resolution of the model, the forcing input, the cloud microphysical scheme used and to the fact that IMERG data tend to overestimate precipitation in extreme conditions. OSTIA-sst rainfall fields had a different spatial distribution, as compared to those of RTG-sst and could overall be characterized as slightly lower. This was due to the changes in the path of the cyclone and the differences in the distribution and magnitude of the SST fields. However, the same point was not presented in RTG-sst-NEMO-sst comparison, where a spatial shift of around 100 mm in the maximum values appeared. These were positively correlated to the deepening of the cyclone in this simulation. MERRA-2 precipitation data distributed similarly to the modeled and satellite rainfall fields but, presented significantly lower magnitudes due to the coarser resolution. Figure 15. (a1) Accumulated precipitation field in millimeters from IMERG satellite data, (a2) accumulated precipitation field in millimeters from MERRA-2 reanalysis (a3) accumulated precipitation in millimeters field from RTG-sst run (a4) differences in accumulated precipitation fields between NEMO-sst and RTG-sst, and (a5) differences in accumulated precipitation fields between OSTIA-sst and RTG-sst, for cyclone Trixi.
For cyclone Numa, maximum rainfall values were distinctly placed on the east coast of Italy (Gulf of Taranto) in a structure that also implied the position of the cyclone center ( Figure 16). Moreover, rain bands in a radius of approximately five degrees were also illustrated in the southern Mediterranean Sea. Modeled, satellite-derived and MERRA-2 reanalysis rain amounts followed similar patterns in most of the affected areas. Overall, an underprediction of model results with respect to satellite accumulated rain fields occurred. This also concerns the rain maxima that emerged in the satellite estimates that were not detectable in the simulations. However, these take place in limited spatial scales. Between the different simulations, minimal differences were presented, noting small deviations between the SST fields. MERRA-2 maximum values could be mainly found in the mainland of Greece and Italy in comparison with that of modeled and remote sensing that were also observed in the Adriatic Sea. The same result was evident in a relevant modeling study of cyclone Numa, also supported by radar data [37]. Figure 16. (a1) Accumulated precipitation fields in millimeters from IMERG satellite data, (a2) accumulated precipitation field in millimeters from MERRA-2 reanalysis, (a3) accumulated precipitation field in millimeters from RTG-sst run (a4) differences in accumulated precipitation fields between NEMO-sst and RTG-sst, and (a5) differences in accumulated precipitation fields between OSTIA-sst and RTG-sst, for cyclone Numa.
Concerning cyclone Zorbas, the accumulated precipitation fields derived from the satellite estimates tended to extreme rainfall values in the South Mediterranean ( Figure 17). Unfortunately, due to the area where the maxima appear, verification with ground observations was not feasible. Cumulatively, remote sensing, simulated and reanalysis fields from MERRA-2 correlated well spatially. Maximum values were also located in the same area. OSTIA-sst and NEMO-sst rain fields were considerably lower, as compared to the RTG. This was the result of the higher initial SST fields in RTG-sst simulation, as compared to the OSTIA-sst and NEMO-sst runs. The warmer sea surfaces provided more energy and water vapor to the system, intensifying it and creating conditions with increased rainfall rates. Figure 17. (a1) Accumulated precipitation fields in millimeter from IMERG satellite data, (a2) accumulated precipitation field in millimeters from MERRA-2 reanalysis, (a3) accumulated precipitation field in millimeters from RTG-sst run (a4) differences in accumulated precipitation fields between NEMO-sst and RTG-sst, and (a5) differences in accumulated precipitation fields between OSTIA-sst and RTG-sst, for cyclone Zorbas.

Summary and Conclusions
In this study, an analysis of the impact of the SST on the lifecycle of Mediterranean cyclones with tropical-like characteristics was attempted following different methodologies. A coupled atmospheric/wave-ocean modeling system was used to investigate the influence of the SST on the formation and evolution of cyclonic events.
The importance of the water surface temperature in the creation, evolution, and intensity of Mediterranean cyclones has received a lot of attention in previous studies. These were mainly focused on the implementation of uniform warm or cold SST anomalies [4,8,59]. The current work approached this issue with the use of static and dynamic SST fields from various sources in the simulation of three cyclonic events. These concern modeled and operational daily and hourly analysis SST results enriched with remote sensing data. Precisely, the RTG, OSTIA, and NEMO SST products with broad implementation were used. Therefore, it is arguable that spatiotemporal deviations in the SST forcing seem to mainly contribute to the differences between the model results as there are no significant differences in the maximum/minimum SST values. All numerical experiments performed with the different SST inputs seemed to be adequate for the modeling of the cyclonic systems. This was also supported by the comparison with remote sensing estimates of wind, wave, and precipitation. The utilization of different SSTs did not seem to affect the forecasting skill drastically. Maximum wind speed variations and spatial rain distributions were sufficiently described.
SST deviations did not lead to considerable deviations in model results regarding cyclone tracks, intensities, and wind speed. The most noticeable differences were found in the surface pressure in the simulations of cyclone Trixi and Zorbas. In contrast, the use of alternative SST products in the simulation of cyclone Numa did not seem to provide different model results. This could be attributed to the distinct dynamics and characteristics of this cyclone with a strong relation to large scale processes in the upper atmosphere and baroclinic instability rather than the influence of the SST dissimilarities.
The analysis of the different model results did not imply significant differences in the evolution of the structural characteristics of the cyclones. The RTG-sst simulations resulted in more intensified tropical-type properties, as shown in the peaks of cyclones Zorbas and Trixi. Both were characterized by a deep-warm core structure. However, cyclone Trixi for most of its lifetime, remained in an extra-tropical stage presenting tropical characteristics for a short period. Conversely, cyclone Numa was limited to an extra-tropical type storm without denoting significant differences between the different simulations.
Considerable deviations emerged in the resulted heat fluxes. RTG products, in general, corresponded to warmer seas as compared to OSTIA and NEMO. As a result, both sensible and latent heat fluxes from RTG-sst simulations exceeded the ones from OSTIA-sst and NEMO-sst ones. Apparently, the absence of diurnal variation, including cooling effects, in the RTG SST led to enhanced upward heat fluxes. Nevertheless, cold SST anomalies from the storm passage were evident in daily temporal intervals. This was indicatively shown in the case of cyclone Zorbas.
The most important differences were found in precipitation exhibiting spatial and intensity differences. The warmer sea surfaces presented in the RTG products provided energy and water vapor to the cyclonic systems and favored conditions with increased rainfall rates. On the other hand, model precipitation rates were lower as compared to satellite products.
Overall, SST conditions play a crucial role in the development and evolution of cyclones in the Mediterranean Sea. Even the employment of widely used products is shown to result in different outcomes and this point should be taken into consideration in forecasting and early warning systems.
Author Contributions: C.S. performed computations, analysis, and wrote the manuscript. P.P. aided in analysis and writing and C.T. contributed to suggestions. G.K. contributed to the final reading and suggestions. All authors have read and agreed to the published version of the manuscript.

Funding: This research was co-funded by the European Union (ESF) and
Greek national funds through the Operational Program "Human Resources Development, Education and Lifelong Learning", project title "Study of Extreme Environmental with Numerical Weather Prediction Models and Stochastic Processes" (project code: MIS 5007050).