The continuous warming of SST in the Eastern Mediterranean Sea has been well-studied in the literature (Samuel-Rhoads et al. 2013; Shaltout and Omstedt 2014; Ozer et al. 2017; Pastor et al. 2018, 2020; Mohamed et al. 2019; Pisano et al. 2020; El-Geziry 2021; García-Monteiro et al. 2022). Most of these studies are conducted spatially on a basin scale and temporally on an annual scale. The effects of warming SST in the Cilician Basin, particularly in the southern Cilician Basin (northern coast of Cyprus), are not well-known, especially on monthly, seasonal, and regional scales. This study investigated the temporal evolution (i.e. annual, seasonal and monthly) of SST in the southern Cilician Basin between 1983 and 2023 with two datasets (i.e. NOAA and ERA5) and between 1995 and 2023 with three datasets (i.e. NOAA, ERA5 and In-Situ) to analyse warming rates and trends and also compared the consistency and differences between the datasets (i.e. satellite, reanalyses and in-situ) to determine the optimal methods to be used in the region for the future studies.
Remarkably, the mean SST and the warming rate of the NOAA (22.14 oC and 0.043 oC/year) and In-Situ (22.23 oC and 0.042 oC/year) datasets are very similar and consistent on the annual scale between 1995 and 2023. On the other hand, the annual mean SST and the warming rate of the ERA5 revealed a higher warming rate (0.054 oC /year) and a lower mean SST (21.85 oC) in the same period (Fig. 2 and Table 1). This difference could be because of either the spatial resolution (gridded points) of the ERA5 in the region, where there is only one gridded point in the studied area or the estimation issues (over or under-estimation) of the ERA5, especially in the coastal regions, as concluded in the literature (Romanou et al. 2010; Jordà et al. 2017). The increasing trends of the annual mean SST on all datasets in this study are consistent and within the same range as the results of several other studies in the literature, which reported an SST increase of about 0.025–0.065 ℃/year in the Levantine Basin in the last few decades (Skliris et al. 2012; Macias et al. 2013; Samuel-Rhoads et al. 2013; Shaltout and Omstedt 2014; Ozer et al. 2017, 2022a; Pastor et al. 2018; Bengil and Mavruk 2019; Pisano et al. 2020; El-Geziry 2021; Saraçoğlu et al. 2021).
The annual mean SST is increased by a total of 1.48 ℃ in the NOAA and 2.12 ℃ in the ERA5 datasets from 1983 to 2023. The increase of the annual mean SST from 1995 to 2023 is a total of 1.20 ℃, 1.51 ℃ and 1.18 ℃, respectively, in the NOAA, ERA5 and In-Situ datasets. However, analyses of the annual mean SST over ten-year periods (1983–1993, 1993–2003, 2003–2013 and 2013–2023) revealed that warming in the region is not linear (Fig. 3 and Table 1). The highest increasing trends in all datasets are recorded between 1993 and 2003 (NOAA 0.106 ℃/year, ERA5 0.142 ℃/year and In-Situ 0.088 ℃/year).
On the other hand, decreasing (negative) trends of the annual mean SST of about − 0.045 ℃/year in the NOAA and − 0.054 ℃/year in the ERA5 datasets are recorded in the region between 1983 and 1993 (Fig. 3 and Table 1). Although no data is available in the In-Situ dataset before 1995, the highest increasing trend recorded between 1995 and 2003 may indicate a decreasing trend in the previous period, as in the NOAA and ERA5 datasets.
The decreasing trend of the annual mean SST between 1983 and 1993 is consistent with the study (El-Geziry 2021), which revealed a decreasing trend of about − 0.06 ℃/year in the southern Levantine Basin between 1975 and 1991. The cooling of the SST in the 1980s and early 1990s may be related to the external oceanic processes such as the AMO, which is highly correlated with the results (Table 3) and was in a negative phase during the period (Skliris et al. 2012; Pisano et al. 2020), or the regime shift that occurred in the 1980s due to volcanic eruptions, which had a global impact on climate (Reid et al. 2016), or internal dynamics such as intense atmospheric forcing and oceanic processes which affect the thermohaline circulation of the Eastern Mediterranean Sea (Schroeder et al. 2017). The highest increasing trend observed between 1993 and 2003, following the cooling of the SST in the region, shall be explained by the recovery of the climate system, enhanced by the anthropogenic forcing.
On the annual scale, all datasets (i.e. NOAA, ERA5 and In-Situ) generally revealed consistent and satisfactory results, especially in temporal evolution, trends and anomalies observed in the annual mean SST (Fig. 2 and Table 1). It is important to note that the highest three annual mean SSTs of the last forty years are recorded in the last six years (Fig. 2), and there has been a continuous and significant positive annual mean SST anomaly in the last ten years (Fig. 4). In all datasets, the warmest annual mean SST was recorded in 2018, followed by 2023, 2021 and 2010, respectively (Fig. 2). Despite 2023 being affected by a strong El Nino, which leads an unprecedented increase in the SST all around the globe (Li et al. 2024), 2018 is still the year that the highest annual mean SST and the highest positive annual mean SST anomalies are recorded in all datasets (Fig. 2 and Fig. 4).
The record-breaking annual mean SST of 2018 is also reported by the study of (Ozer et al. 2022b) conducted in the southern Levantine Basin, which associated the warming in 2018 with the BIOS mechanism (reversal of the North Ionian Gyre). The BIOS mechanism affects the thermohaline circulation of the Eastern Mediterranean Sea by limiting the exchange of water masses (i.e. Atlantic Water and the LIW) at the Strait of Sicily (Ozer et al. 2017; Schroeder et al. 2017; Estournel et al. 2021).
This study does not aim to discuss the thermohaline circulation or circulation time of surface waters; however, it shall be noted that the effects of the BIOS mechanism may not be similar simultaneously in different spaces, and it is not the first time that the Northern Ionian Gyre has reversed, but such warming, as experienced in 2018, has never been evidenced. The warming in 2018 may be related to the BIOS mechanism, as concluded by (Ozer et al. 2022b), or along with the BIOS mechanism, it may be related to a combination of regional atmospheric and oceanic conditions (i.e. changes in wind regimes, atmospheric pressure and circulation), which requires further investigation.
The record-breaking annual mean SST and anomalies recorded in 2018 are mainly caused by the unusual warming observed in the winter and especially in the spring of 2018 (Fig. 6), which may also be related to and influenced by the high SSTs recorded in the summer and fall of 2017. The summer of 2017 was one of the hottest seasons of the last forty years in all datasets, and the fall of 2017 was recorded as one of the hottest in the In-Situ dataset (not observed in the NOAA and ERA5 datasets) (Fig. 6). It should be noted that with the lack of the local Poyraz wind (not shown here), a cold and dry north-northeasterly wind blowing in mid-fall and winter, the SST is not cooling down efficiently in the Cilician Basin.
Contrary to the annual mean SST, significant differences between the datasets were observed in monthly and seasonal SST trends and warming rates. At this point, it is worth mentioning that the climatological seasons in the Eastern Mediterranean Sea are defined as long, hot and dry summers with mild and wet winters. Early spring and mid-fall shall be characterised as the transition period between summer and winter (Özsoy et al. 1989). In the Levantine Basin, the time lag of the effect of atmospheric forcing (i.e. air temperature) on the SST is about fifteen to thirty days (Ozer et al. 2022a). Therefore, in this study, the seasons (oceanographical seasons) are defined according to the monthly mean SSTs, in which January, February, and March are classified as Winter; April, May, and June are classified as Spring; July, August, and September are classified as Summer; and October, November, and December are classified as Fall.
The oceanographic spring season, particularly June, has the lowest increasing trend in all datasets between 1995 and 2023. Moreover, in June, even decreasing trends are recorded in the NOAA (-0.002 oC/year) and In-Situ (-0.017 oC/year) datasets in the same period. The ERA5 dataset revealed an increasing trend with 0.030 oC/year in June (Fig. 5 and Fig. 6). This result contradicts the results of (Skliris et al. 2012; Pastor et al. 2018; El-Geziry 2021), which concluded that a high warming rate of the SST in the region occurred in the spring.
On the other hand, variabilities are observed on the highest warming trend between all datasets on monthly and seasonal scales. In the NOAA dataset, the highest warming rate is recorded in early to mid-fall (September to November). In the ERA5 dataset, the highest warming occurred in mid-summer (August, September) and November. In the In-Situ dataset, the highest warming trend is recorded in the fall, especially in November and December (Fig. 5 and Fig. 6).
One of the most significant differences between the datasets on a monthly/seasonal scale is observed in November and December of 2017 and 2018, with excessive warming recorded in the In-Situ dataset but not in the NOAA and ERA5 datasets. This difference between the datasets shall be explained by local atmospheric and oceanographic processes, which cannot monitored or measured by satellite and reanalyses, or it should be an observational error. The latter is unlikely due to the observation methodology, and it was observed consecutively in 2017 and 2018.
The results of the monthly/seasonal ERA5 dataset in this study revealed consistent and similar results to the studies that conclude the highest increasing trend on the SST is observed in the summer season (Samuel-Rhoads et al. 2013; Shaltout and Omstedt 2014). However, contrary to the ERA5 dataset of this study and the studies mentioned above, the highest increasing trend of the SST in the NOAA and In-Situ datasets from 1995 to 2023 is recorded in the fall, which supports the results of (Saraçoğlu et al. 2021) conducted in the northern Cilician Basin.
Probably the most important outcome of this study is the warming of the SST recorded in the fall. In the Levantine Basin, the SST in late fall (November and December) is crucial for vertical mixing processes and water mass formations. Cooling of the SST in November and December is the precondition of these processes (The POEM Group 1992; The LIWEX Group 2003; Kucuksezgin and Pazi 2006; Salihoglu et al. 2019; Deliceirmak and Salihoglu 2020; Fach et al. 2021). Therefore, excessive warming recorded in these months may affect the vertical mixing and stratification of the water column and the formation of the LIW, hence the circulation and the coastal ecosystem.
Furthermore, as concluded in the past (Hurrell and Deser 2009; Trenberth 2009; Deser et al. 2010; Gastineau and Frankignoul 2015; Pisano et al. 2020; Rencurrel and Rose 2020), the results of this study revealed that the rising trend of SST is non-linear on the temporal scale. Hence, the effect of warming SST on the seawater’s physical properties may be observed in long-term (i.e. decadal) and short-term (i.e. annual, seasonal and monthly) scales. The short-term effect of warming SST shall be reflected as salinity shifts and changes in physical and biochemical oceanographic characteristics observed on the upper thermocline in the Levantine Basin (Ozer et al. 2017; Fach et al. 2021; Deliceirmak et al. 2024, unpublished document) and may be correlated with the increasing SST in the region, especially in the last ten years.
As a final point, methods used widely in ocean sciences, such as reanalysis, remote sensing and model simulation, are continually improving, and they are powerful and beneficial tools, primarily due to temporal and spatial coverage and resolution. However, limitations in their estimations, especially in coastal regions, are subject to discussions (Romanou et al. 2010; Jordà et al. 2017; O’Carroll et al. 2019; Pisano et al. 2020; Meng et al. 2023). Keeping this in mind, the encouraging results of this study may shed light on improving the accuracy of the methods mentioned above for coastal regions and, if consistently verified, may encourage the use of in-situ data sets of coastal meteorological stations as a supplementary or validation method in coastal ocean studies.