Increasing Trends in Air and Sea Surface Temperature in the Central Adriatic Sea (Croatia)

: The Adriatic Sea and its coastal region have experienced signiﬁcant environmental changes in recent decades, aggravated by climate change. The most prominent effects of climate change (namely, an increase in sea surface and air temperature together with changes in the precipitation regime) could have an adverse effect on social and environmental processes. In this study, we analyzed the time series of sea surface temperature and air temperature measured at three meteorological stations in the Croatian part of the Adriatic Sea. To assess the trends and variations in the time series of sea surface and air temperature, different statistical methods were employed, i.e., linear and quadratic regressions, Mann–Kendall test, Rescaled Adjusted Partial Sums method, and autocorrelation. The results evidenced increasing trends in the mean annual sea surface temperature and air temperature; furthermore, sudden variations in values were observed in 1998 and 1992, respectively. Increasing trends in the mean monthly sea surface temperature and air temperature occurred in the warmer parts of the year (from March to August). The results of this study could provide a foundation for stakeholders, decision–makers, and other scientists for developing effective measures to mitigate the negative effects of climate change in the scattered environment of the Adriatic islands and coastal region.


Introduction
The Adriatic Sea is the northernmost part of the Mediterranean Sea and represents a vast bay that is indented far into the European mainland. The surface of the Adriatic Sea is 138,595 km 2 , which constitutes 4.6% of the total surface of the Mediterranean Sea. The Adriatic Sea stretches in NW-SE direction for 870 km, with an average width of approximately 200 km. It is connected to the Mediterranean Sea by the 70 km-wide Strait of Otranto. The Adriatic Sea is heterogeneous in terms of physical, chemical, and biological properties [1]. Based on the average depth, the Adriatic is divided into (i) very shallow Northern Adriatic (<35 m), (ii) Central Adriatic Sea (~140 m), and (iii) Southern Adriatic with an average depth of >200 m and maximum depth up to 1228 m [2]. The majority of the islands in the Adriatic Sea are within the Croatian territory. The total length of the Croatian coast is 5835 km, i.e., the lengths of the island and mainland coast are 4058 km and 1777 km, respectively.
One of the key environmental questions is how climate change, particularly global warming, will affect the environmental and social processes in different regions of the planet. The investigation of the effects of climate change in the coastal environment is a particularly complex task because of the simultaneous influence of vast land and water masses, which react differently to climate change. Since climate change continues to intensify, both the Adriatic and the Mediterranean regions are considered hot spots of climate change. In particular, the Adriatic ecosystems suffer the combined effect of climate change in the Northern hemisphere and local climate variations [3]. Regional atmospheric oscillations, as a result of the air pressure variations in the Northern hemisphere, and the intensification of water inflow from the Eastern Mediterranean affected the Adriatic Sea temperature and salinity. Further local factors impacting the Adriatic marine ecosystems are: (i) increasing sea surface temperature, (ii) negative trends in the precipitation patterns, particularly in the winter season, and (iii) reduced inflow of fresh water and nutrients as a result of the decreasing inflow from the Po river [4]. This complex situation is fostered by the variable geometry and topography of the Adriatic region where each bay, island, or channel in the Adriatic Sea is site-specific in terms of oceanographic properties [5] and reactions to climate change [6,7]. The devastation of the environment, the increasing demand for fresh water for potable and irrigation purposes, and the changes in the already unfavorable water balances of the Mediterranean semi-arid climate could lead to adverse environmental and social consequences in the near future. Commonly, the karst environments are scarce in surface hydrography and groundwater is the main freshwater resource. Hence, the Mediterranean water resources are under significant stress due to over-abstraction, climate change, and the high possibility of seawater intrusions into karst aquifers in the coastal areas or the islands.
To understand the trends and the behavior of the air temperature on small islands and in the coastal areas, it is necessary to understand the interaction between the sea surface temperature and the air temperature. Vlahakis and Pollatou [8] emphasized that the sea surface temperature is the key factor in the assessment of the climate and climaterelated processes on all scales, especially on the islands and in the coastal areas. Sea surface temperature has a direct influence on global energy transfer, atmospheric processes, precipitation, evapotranspiration, moisture in air and soil, wind development, hydrological cycle, and other ecological or social processes [9]. However, sea surface temperature is not considered a standard meteorological parameter even though its effect on the atmosphere is immense. Furthermore, the time series of the sea surface temperature are considerably shorter when compared to the air temperature. The sea has a higher thermal capacity than the land, resulting in slower heating and cooling. During winter and summer, the sea acts as a buffer that moderates the air temperature, resulting in a lower range of the sea surface temperature than the range of air temperature over land [10]. The flow of the seawater and other turbulent processes cause the deeper layers of the sea to mix with the surficial ones, influencing the air temperature. The air temperature reacts differently to changes in the sea surface temperature, and this is particularly pronounced on the small islands and in the coastal areas. Therefore, it is necessary to analyze each location based on a detailed and reliable time series of measured data.
Time series of the sea surface temperature and the air temperature measured at three meteorological stations located on Central Adriatic islands (Hvar and Vis) and on the mainland (Split) were analyzed in this study. Firstly, we will describe the study area and the available dataset. The dataset will be analyzed using different statistical approaches and the obtained results will be used to assess the occurrence of trends within the time series of the sea surface temperature and the air temperature. These data will provide a foundation for stakeholders, decision-makers, and other scientists for developing effective measures to mitigate the negative effects of climate change in the Adriatic region and beyond.

Study Area
The study area is located in the central part of the Adriatic Sea and consists of the two meteorological stations on the islands of Hvar and Vis, and the meteorological station in the city of Split on the mainland. The map of the study area showing the analyzed meteorological stations is in Figure 1. The main characteristics of the analyzed stations used in this study are summarized in Table 1. Note: H-the altitude of the meteorological station; SST-sea surface temperature; TA-air temperature; L-the distance between the SST and the air temperature measuring point; Ls-minimum distance between meteorological stations and the seashore.
The meteorological station Split is located in the eastern part of Split, on the Marjan peninsula. Marjan is a forest park used for recreational purposes. The meteorological station is in the vicinity of the second tallest peak on Marjan (178 m a.s.l.). In 2011, the population of Split was 178,192. Despite rapid urbanization of the Split agglomeration, the meteorological station Split is not affected by urbanization due to its position within a protected forest area. The sea surface temperature (SST hereafter) is measured at the endpoint of the pier at the easternmost point of the Marjan peninsula.
The small island of Hvar belongs to the Middle-Dalmatian island group. With a surface area of 297.32 km 2 and a coastline of 270 km, Hvar is the fourth biggest island in the Croatian part of the Adriatic Sea [11]. In 2011, the city of Hvar had a population of 4251. The meteorological station Hvar is situated in a small forest grove, away from the The main characteristics of the analyzed stations used in this study are summarized in Table 1. Note: H-the altitude of the meteorological station; SST-sea surface temperature; T A -air temperature; L-the distance between the SST and the air temperature measuring point; L s -minimum distance between meteorological stations and the seashore.
The meteorological station Split is located in the eastern part of Split, on the Marjan peninsula. Marjan is a forest park used for recreational purposes. The meteorological station is in the vicinity of the second tallest peak on Marjan (178 m a.s.l.). In 2011, the population of Split was 178,192. Despite rapid urbanization of the Split agglomeration, the meteorological station Split is not affected by urbanization due to its position within a protected forest area. The sea surface temperature (SST hereafter) is measured at the endpoint of the pier at the easternmost point of the Marjan peninsula.
The small island of Hvar belongs to the Middle-Dalmatian island group. With a surface area of 297.32 km 2 and a coastline of 270 km, Hvar is the fourth biggest island in the Croatian part of the Adriatic Sea [11]. In 2011, the city of Hvar had a population of 4251. The meteorological station Hvar is situated in a small forest grove, away from the urbanized city center. The SST is measured at the endpoint of a small pier 100 m from the meteorological station.
Vis is a small remote island in the Adriatic Sea. Its distance from the mainland is 43 km and the island is exposed to very strong winds. With a surface area of 89.72 km 2 and coastline of 85 km, Vis is the ninth biggest island in the Croatian part of the Adriatic Sea [11]. The meteorological station is located in the city of Komiža, on the northern edge of the urbanized area. In 2011, Komiža had a population of 1526. The SST is measured at the endpoint of a pier in Komiža, 190 m from the meteorological station. Vis Island has favorable geological and hydrogeological conditions that have enabled the formation of high-quality karst aquifers from which the fresh water is abstracted. Hence, Vis Island and its groundwater resources are considered highly vulnerable to climate change [7].
According to Köppen-Geiger climate classification, the study area is characterized by the Csa climate type [12], sometimes called the "olive" climate. It is a semiarid variety of Mediterranean climate characterized by mild and humid winters, and dry and hot summers.

Data Collection
The SST and the air temperature data used in this study were provided by the Croatian Meteorological and Hydrological Service (DHMZ). The conventional method of measuring the SST is performed using a thermometer immersed in the sea at a depth of 30 cm, given that the sea is not shallower than 180 cm at the measuring point. The thermometer is immersed for three minutes, and after that, it is taken out and read quickly [13]. The SST is measured three times a day at 7, 14, and 21 h local time. The measurement of the SST is performed by personnel from the nearby meteorological station. Values of the mean monthly and the mean annual SST were analyzed in this study. Furthermore, the time series of air temperature were measured at main meteorological stations. An interesting fact is that the meteorological station Hvar started to operate in 1858, followed by the Split station a year after, and Komiža in 1959. The World Meteorological Organization (WMO) recognized the quality of the station in Hvar and awarded it a Centennial Observing Station title. The air temperature is measured hourly at a standard height of 2 m above the ground. In this study, the mean monthly and the mean annual air temperature for the period 1960-2019 were analyzed. Despite several approaches (e.g., [14][15][16]), the method employed by the DHMZ [13], which is the most common in Europe, defines the mean daily air temperature as: where t 7 , t 14 , and t 21 are the air temperature values measured at 7, 14, and 21 h (local time), respectively. The same procedure is applied for the SST.

Statistical Methods
Linear and quadratic regressions were performed on the time series of mean monthly and mean annual SST and air temperature from three stations analyzed in this study. The linear regression equation is given as: and the quadratic regression equation as: where T is the mean monthly or mean annual SST or air temperature in year t, a and b are linear regression coefficients, and c, d, and e are quadratic regression coefficients. All five coefficients are calculated by the least-squares method. The coefficient a represents the slope of the regression line whose dimension is • C/year, and it is the indicator of the average intensity of the increasing or decreasing trends in the values of the analyzed time series. The correlation coefficients r 2 and R 2 were calculated for the linear and the quadratic regressions, respectively. Both coefficients show the strength and the direction of linear and quadratic correlation between variables x (time) and y (the mean monthly or the mean annual SST or air temperature).
To assess whether the time series have monotonic increasing or decreasing trends, the Mann-Kendall (M-K hereafter) non-parametric test was used [17]. The null hypothesis for this test is that there is no monotonic trend within the analyzed time series, while the alternate hypothesis is that the trend exists. As a criterion to accept the alternate hypothesis (i.e., the presence of an increasing or decreasing trend), p-values < 0.05 were used in this study.
Furthermore, the Rescaled Adjusted Partial Sums method (RAPS hereafter) was used to detect statistically significant peaks or declines in values (i.e., the trends variations) within the analyzed time series [18,19]. This method allows overcoming random and irregular fluctuations as well as rough errors in values within the time series, which may be hidden from the common plots of values of the time series. Based on the RAPS results, sub-periods with similar characteristics or a larger number of trends within the time series could be distinguished. The formula for the calculation of RAPS is: where Y t is the value of the observed parameter at time t, Y m is the mean value of observed time series, S y is the standard deviation of the observed time series, and k is the number of observations. The breakpoints between the sub-periods were established when the trend of RAPS results showed a significant variation. The differences in statistical parameters of the two neighboring sub-periods defined by the RAPS were evaluated by the F-test and the t-test [20]. In particular, the F-test was used to assess the equality of variances between the two normally distributed populations (i.e., sub-periods). The t-test was used to determine whether there is a statistical difference between the mean values of the two sub-periods. In this study, both tests accept the null hypothesis for p-values < 0.05.
Furthermore, the autocorrelation of the time series was determined. Autocorrelation is a mathematical function representing the degree of similarity between the specific time series and a lagged version of the same time series over successive time intervals. Autocorrelation coefficient r ranges from −1 to 1 and it measures the strength of a relationship between the current value of the variable with its shifted value. In this study, the interval of the shifting variable was set to 1 year. For r < 0.2, the time series is not autocorrelated meaning that the behavior of the values does not depend on the previous values [21].

Results and Discussion
In this chapter, temporal changes in the mean annual and the mean monthly SST and air temperature were analyzed. The analyzed time series are not identical in duration, which will affect the reliable comparison of the results obtained at three analyzed stations to a minor extent. Despite the differences in the duration of the analyzed time series, contemporaneous time series were available for the last 29-year period (i.e., 1991-2019), when the most significant increases in the SST and the air temperature were observed. This fact will allow a more reliable conclusion about the recent and future behavior of the SST and the air temperature in the study area. Table 2 shows a summary of the exploratory statistical analysis (minimum, average, maximum, and range) of the mean annual SST, the air temperature (T A ), and their difference (∆T = SST-T A ) at the analyzed stations. The longest time series analyzed in this study were recorded in Split (from 1960 to 2019), followed by Hvar (from 1964 to 2019), and Komiža with the shortest time series (from 1991 to 2019). Table 2. Statistics (minimum, average, maximum) of the mean annual SST, air temperature (T A ), and their difference (∆T = SST-T A ). The minimum, average, and maximum values of the mean annual SST measured in Split were 16 The ∆T values (SST-T A ) showed the same distribution as the SST and the T A . The ∆T of the minimum values were 1.2 • C, 1.5 • C, and 1.7 • C, in Split, Hvar, and Komiža, respectively. The ∆T average values were 1.1 • C, 1.4 • C, and 1.7 • C, while the ∆T maximum values were 0.9 • C, 1.1 • C, and 1.5 • C at the same stations, respectively. The positive values of ∆T indicated that the mean annual SST was always higher than the mean annual air temperature. It should be noted that ∆T of the minimum values were higher than the ∆T of the maximum values due to the smaller amplitude of the SST than the air temperature. Figure 2 shows the time series of the mean annual SST measured at Split (shown in blue), Hvar (shown in green), and Komiža (shown in red) stations. The linear regressions evidenced increasing trends in the SST at all analyzed stations, which were corroborated by the results of the M-K test (i.e., p < 0.01). However, to achieve better fitting to the data, quadratic regressions were performed on the time series of Split and Hvar stations. The results showed slightly higher R 2 values than the coefficient of linear regressions r 2 . Quadratic regression was not performed on the time series from Komiža due to the missing data before 1991.

Analyses of the Mean Annual Sea Surface Temperature and Air Temperature
The RAPS method has been used on the time series of SST from all analyzed stations. The results are in the Supplementary Material (S1) and they evidenced the presence of two sub-periods: (i) from the beginning of the respective time series until 1997, and (ii) from 1998 to 2019 ( Figure 3).
Despite the differences in the duration of the time series, increasing (Hvar and Komiža) or decreasing (Split) trends in SST were not statistically significant in the first sub-period defined by the RAPS (p-values > 0.05; Figure 3). Statistically significant increasing trends in the SST occurred in the second sub-period at stations in Hvar and Split (p-values < 0.05), and in Komiža there was not a statistically significant trend in the second sub-period (p-values > 0.05).
The average values of the mean annual SST within the sub-periods defined by the RAPS method are shown in Figure 3 and Table 3. The statistical analyses evidenced statistically significant differences between the average values of the mean annual SST in the two sub-periods, with the rejection of the null hypothesis of the t-test (low p-values, i.e., p < 0.01), and similar variances of the sub-periods reflecting the failure to reject the null hypothesis of the F-test (high p-value). The linear regressions evidenced increasing trends in the SST at all analyzed stations, which were corroborated by the results of the M-K test (i.e., p < 0.01). However, to achieve better fitting to the data, quadratic regressions were performed on the time series of Split and Hvar stations. The results showed slightly higher R 2 values than the coefficient of linear regressions r 2 . Quadratic regression was not performed on the time series from Komiža due to the missing data before 1991.
The RAPS method has been used on the time series of SST from all analyzed stations. The results are in the Supplementary Material (S1) and they evidenced the presence of two sub-periods: (i) from the beginning of the respective time series until 1997, and (ii) from 1998 to 2019 (Figure 3).
Despite the differences in the duration of the time series, increasing (Hvar and Komiža) or decreasing (Split) trends in SST were not statistically significant in the first sub-period defined by the RAPS (p-values > 0.05; Figure 3). Statistically significant increasing trends in the SST occurred in the second sub-period at stations in Hvar and Split (p-values < 0.05), and in Komiža there was not a statistically significant trend in the second sub-period (p-values > 0.05).
The average values of the mean annual SST within the sub-periods defined by the RAPS method are shown in Figure 3 and Table 3. The statistical analyses evidenced statistically significant differences between the average values of the mean annual SST in the two sub-periods, with the rejection of the null hypothesis of the t-test (low p-values, i.e., p < 0.01), and similar variances of the sub-periods reflecting the failure to reject the null hypothesis of the F-test (high p-value).   Figure 4 shows the time series of the mean annual air temperatures measured at Split, Hvar, and Komiža stations. The average values of the mean annual air temperature during the analyzed periods were identical at Hvar and Komiža (16.7 °C), while in Split they were slightly lower (16.4 °C). All three stations showed statistically significant increasing trends in the mean annual air temperature (i.e., low p-values < 0.01). To achieve better fitting to the data, quadratic regressions were preferred over simple linear regression for the time series from all analyzed stations.
The RAPS method evidenced the presence of two sub-periods in the mean annual air temperature time series: (i) 1960-1991, and (ii) 1992-2019 ( Figure 5). The M-K test evidenced that the trends within the first sub-period were statistically insignificant (pvalues > 0.05) at all analyzed stations. However, within the second sub-period statistically significant increasing trends were observed at Hvar and Split stations (p < 0.01).   Figure 4 shows the time series of the mean annual air temperatures measured at Split, Hvar, and Komiža stations. The average values of the mean annual air temperature during the analyzed periods were identical at Hvar and Komiža (16.7 • C), while in Split they were slightly lower (16.4 • C). All three stations showed statistically significant increasing trends in the mean annual air temperature (i.e., low p-values < 0.01). To achieve better fitting to the data, quadratic regressions were preferred over simple linear regression for the time series from all analyzed stations.   The average values of the mean annual air temperature within a sub-period defined by the RAPS method are shown in Figure 5 and Table 4. The statistical analyses evidenced statistically significant differences between the average values the mean annual air temperature in the two sub-periods, with the rejection of the null hypothesis of the t-test (low p-values, i.e., p < 0.01), and similar variances of the sub-periods reflecting the failure to reject the null hypothesis of the F-test (high p-value). The results indicated that the air temperature had started to increase considerably at the beginning of the 1990s at all analyzed stations. These results fit the regional warming patterns observed in Croatia and the western Balkans [22][23][24][25]. Furthermore, it can be concluded that the rapid increase in air temperature had occurred 6 years before the increase in SST at all analyzed stations.
The correlation between the mean annual SST and the mean annual air temperature time series was the highest in Hvar, with the values of r 2 = 0.796 in the period from 1964 to 2019, followed by Split with r 2 = 0.688 in the period from 1960 to 2019, and the lowest was in Komiža with r 2 = 0.6183 in the period from 1991 to 2019. Table 5 shows the r-squared values of the linear correlation coefficients (i) between the pairs of the time series of the mean annual SST during periods of contemporaneous measurements at all three stations (from 1991 to 2019), and (ii) between the pairs of the time series of the mean annual air temperature (from 1960 to 2019).    The high values of r 2 indicated the similarity of the SST and the air temperature regimes at analyzed stations. The highest r 2 from the SST time series was observed between the closest stations, Komiža and Hvar, r 2 = 0.87, and the lowest between Komiža and Split, r 2 = 0.75. The highest r 2 value from the time series of the mean annual air temperature was observed between Split and Hvar stations, r 2 = 0.95, and the lowest between Hvar and Komiža, r 2 = 0.80. Furthermore, the autocorrelation method was performed on the time series of the mean annual SST and air temperature measured at Split and Hvar stations, for the period from 1960 to 2019, and from 1964 to 2019, respectively. The time series from Komiža did not qualify for autocorrelation due to the insufficient duration of the time series of SST (from 1991 to 2019). The results indicated the similar behavior of the SST and the air temperature in Hvar and the air temperature in Split having a long-term autocorrelation ( Figure 6). However, the autocorrelogram of the SST from the Split station is significantly different. The values of r were steady at approximately 0.5 until 6 years when a significant drop occurred. After 8 years, the values of the autocorrelation coefficient were lower than the significance threshold (0.2) meaning that the "memory of the system" was lost. Plausible causes that decreased the correlation of the time series of the mean annual SST include higher variability of SST in Split, pronounced coastal effect and local variability of climate, and bay-like topography. Considering spatially and temporally limited data measured at the station in Split, a more detailed study should be conducted to evaluate the driving force of this different behavior.
J. Mar. Sci. Eng. 2021, 9, x FOR PEER REVIEW 11 of 17 temperature in Hvar and the air temperature in Split having a long-term autocorrelation ( Figure 6). However, the autocorrelogram of the SST from the Split station is significantly different. The values of r were steady at approximately 0.5 until 6 years when a significant drop occurred. After 8 years, the values of the autocorrelation coefficient were lower than the significance threshold (0.2) meaning that the "memory of the system" was lost. Plausible causes that decreased the correlation of the time series of the mean annual SST include higher variability of SST in Split, pronounced coastal effect and local variability of climate, and bay-like topography. Considering spatially and temporally limited data measured at the station in Split, a more detailed study should be conducted to evaluate the driving force of this different behavior.

Analyses of the Mean Monthly Sea Surface Temperature and Air Temperature
The summary of the statistical analysis (minimum, average, maximum, and range) of the mean monthly SST and air temperature (TA) time series, as well as their differences (ΔT = SST-TA), is shown in Table S2 of the supplementary material. The statistical analyses evidenced that the ΔT of the minimum values coincided at all three stations during the warmer period of the year (i.e., from May to August), and ranged from -1.8 to 6.4 °C in Split, from -2.6 to 7 °C in Hvar, and from -1.8 to 7 °C in Komiža. The amplitude was slightly higher at island stations (i.e., Hvar and Komiža) than at the Split station. The ΔT of the maximum values were the lowest during the warmer period of the year (from April to September) and the highest during December at all three stations, and ranged from -4.6 to 5.3 °C in Split, from -3.1 to 5 °C in Hvar, and from -2.8 to 4.7 °C in Komiža. The distribution of the ΔT of the average values showed a similar pattern as the ΔT of minimum and maximum, and the values were the lowest in July and the highest in December at all analyzed stations. The ΔT of the average values ranged from -2.6 to 5 °C in Split, from -2.2 to 5.2 °C in Hvar, and from -1.6 to 5.1 °C in Komiža. The ranges of the mean monthly ΔT values were the lowest (i.e., highest negative values) during the winter period (from January to March), and the highest during warmer periods of the year at all

Analyses of the Mean Monthly Sea Surface Temperature and Air Temperature
The summary of the statistical analysis (minimum, average, maximum, and range) of the mean monthly SST and air temperature (T A ) time series, as well as their differences (∆T = SST-T A ), is shown in Table S2 of the supplementary material. The statistical analyses evidenced that the ∆T of the minimum values coincided at all three stations during the warmer period of the year (i.e., from May to August), and ranged from −1.8 to 6.4 • C in Split, from −2.6 to 7 • C in Hvar, and from −1.8 to 7 • C in Komiža. The amplitude was slightly higher at island stations (i.e., Hvar and Komiža) than at the Split station. The ∆T of the maximum values were the lowest during the warmer period of the year (from April to September) and the highest during December at all three stations, and ranged from −4.6 to 5.3 • C in Split, from −3.1 to 5 • C in Hvar, and from −2.8 to 4.7 • C in Komiža. The distribution of the ∆T of the average values showed a similar pattern as the ∆T of minimum and maximum, and the values were the lowest in July and the highest in December at all analyzed stations. The ∆T of the average values ranged from −2.6 to 5 • C in Split, from −2.2 to 5.2 • C in Hvar, and from −1.6 to 5.1 • C in Komiža. The ranges of the mean monthly ∆T values were the lowest (i.e., highest negative values) during the winter period (from January to March), and the highest during warmer periods of the year at all analyzed stations. The ranges of the ∆T values were from −4.4 to −0.8 • C in Split, from −4.1 to -0.2 • C in Hvar, and from −3.7 to −0.8 • C in Komiža. At all analyzed stations, the ranges of the mean monthly air temperature were significantly higher than the ranges of the SST in each month of the year.
The average values of the mean monthly SST, air temperature, as well as their difference (∆T = SST-T A ), are shown in Figure 7. Figure 7 shows the comparison of the average values of the mean monthly SST and air temperatures at Split, Hvar, and Komiža stations. In the warmer part of the year (from May to August), the air temperature was higher than the SST at all analyzed stations. Furthermore, the air temperature and the SST were nearly identical in April at stations in Split and Hvar. The most significant difference occurred in July when their difference was 2.64 • C in Split, 2.25 • C in Hvar, and 1.63 • C in Komiža. In the colder parts of the year, the SST was higher than the air temperature at all analyzed stations, with the highest difference in December, when the mean monthly SST was on average 5 • C higher than the mean monthly air temperature. The results evidenced a very similar behavior of temperature at all analyzed stations, despite the differences in the duration of the time series. The smallest ∆T values were observed at the station in Komiža, and the highest at the Split station. This fact could be partly explained by the differences in the duration of the time series, but also by the local effect of the position of the meteorological station and its distance from the SST measuring point. Furthermore, the position of the meteorological station in terms of the distance from the landmass also plays an important role. Table 6 shows the slope of the linear equation, a, squared values of the linear correlation coefficient, r 2 , and M-K probability values, p, for the analyzed time series of the mean monthly SST. The results indicated that the statistically significant increasing trends in SST were observed in the warmer parts of the year (March-August) in Split, throughout the year in Hvar, while the statistically more complex situation was observed in Komiža, where increasing trends occurred in March, June, July, September, and December. The average values of the mean monthly SST, air temperature, as well as their difference (ΔT = SST-TA), are shown in Figure 7.    Table 7 shows the slope of the linear equation, a, squared values of the linear correlation coefficient, r 2 , and M-K probability values, p, calculated from the time series of the mean monthly air temperature. The results indicated the nearly identical behavior of the air temperature at all analyzed stations. Statistically significant increasing trends in the mean monthly air temperature occurred from March to August at all analyzed stations, but also in December at the Split and Komiža stations. Statistical analyses performed on the time series of the mean monthly SST and air temperature undoubtedly evidenced that the most significant increasing trends in the SST and the air temperature occurred during warmer parts of the year, i.e., during spring and summer. Similar warming trends could be present over the entire Adriatic Sea and its coast, but further detailed studies on the time series from the other meteorological stations, coupled with studies using gridded datasets from remote sensing or numerical simulation, are needed to assess whether these trends are present on a regional scale. Furthermore, the results of this study are in concordance with findings from Bartolini et al. [23], who conducted a regional climatological study where they analyzed air temperatures at 21 stations in the Mediterranean region, i.e., in Tuscany, Italy, and they have concluded that the most rapid and intensive warming trends occurred from March to August. Figure 8 shows the ratio of the average values of the mean monthly SST and air temperature in the period of contemporaneous measurements at all three stations (from 1991 to 2019). The data loops from all three stations are similar in shape but are slightly shifted in values. The analyses of the time series from the period of contemporaneous measurements confirmed the previous conclusion which was based on the divergent time series. Figure 8 shows the ratio of the average values of the mean monthly SST and air temperature in the period of contemporaneous measurements at all three stations (from 1991 to 2019). The data loops from all three stations are similar in shape but are slightly shifted in values. The analyses of the time series from the period of contemporaneous measurements confirmed the previous conclusion which was based on the divergent time series. Figure 8. The ratio of the mean monthly SST and air temperature from the period of contemporaneous measurements at all three analyzed stations. Figure 8. The ratio of the mean monthly SST and air temperature from the period of contemporaneous measurements at all three analyzed stations.

Conclusions
In the past 40 years, the Adriatic Sea and adjacent coastal areas have faced an increase in sea surface temperature, air temperature, and changes in the precipitation regime. Statistical analyses conducted within this study evidenced increasing trends in both the investigated temperature time series (i.e., SST and T A ). The results of RAPS indicated that the increase has been sharper since the 1990s but it occurred with a significant temporal shift (6 years) between the mean annual air temperature and SST. The observed lag in the warming of the Adriatic Sea is most likely a result of the slower response of the sea to the warming process, due to the inherent ability of the sea to absorb vast amounts of energy. Furthermore, the most significant increasing trends in the mean monthly air temperature and SST occurred during warmer parts of the year, i.e., during spring and summer. These results are in accordance with regional climate models [25,26] for the Adriatic Sea.
The climate changes described in this and other works (e.g., [22][23][24]) have a strong impact on the environment, the marine species, and the population living in the Adriatic region. For example, the habitats of many thermophilic species had migrated horizontally and vertically towards the deeper and the colder parts of the sea. If the sea surface temperature will continue to increase, the geographical distribution of these species will continue to decrease and will eventually cause extirpation and possibly even extinction. In addition, changes in the composition and quantity of zooplankton were observed, particularly in coastal areas of the Adriatic. Moreover, frequent blooms of marine phytoplankton and the spread of bacteria and thermophilic species of tropical algae were also observed [4].
Regional climate models for the Mediterranean region showed the continuation or even an increase in warming trends, and therefore, it is realistic to assume that the negative changes in the Adriatic Sea and its coast will be even more pronounced in the near future. In conjunction with increasing anthropogenic pressures, such as overfishing, urban and industrial pollution, the devastation of habitats, seasonal tourism pressures, and hydrocarbon exploitation, the negative consequences could be even more drastic. The lack of reliable indicators of climate change or variability, i.e., the sea surface and the air temperature, measured over a dense network of meteorological stations, has a significant influence on the development of effective measures for mitigation of negative effects of climate change [27].
Furthermore, as a result of the variable distance from the mainland, and local or regional topography, the effects of climate change manifest differently in specific islands or coastal regions. Limited and vulnerable groundwater resources along the coast of the Croatian part of the Adriatic Sea and in the related islands, in combination with unsustainable anthropic activities (e.g., mass tourism, land-use changes, groundwater over-abstraction, and urbanization), significantly reduce the options for adaption to current and future climate change. Due to the vast cultural, historical, social, geographical, and biological diversity, the Mediterranean region requires urgent and effective measures that will foster its sustainable development. The fundamental problem is that the Croatian part of the Adriatic Sea, similar to the other countries in the Mediterranean region, does not have a sufficiently dense network of meteorological stations and sufficiently long time series of measured data, especially sea surface temperature data.
The availability of high-resolution data on climate change and variability could enable island communities to enhance their resilience and design site-specific measures to mitigate possible negative effects on water availability and natural ecosystems. Besides the structural modifications (e.g., re-use of purified domestic and industrial wastewater, reduction in losses from water supply systems, desalinization plants, managed aquifer recharge), a holistic approach could also be fostered by increasing awareness and education of the local population on correct utilization of the water resource (e.g., promotion of water savings during dry months, rainwater harvesting, planting of crops that require little or no irrigation, reduction in carbon footprint, and preservation of ecosystems and their services to mitigate floods and droughts).
The authors hope that this study will contribute to a better understanding of this topic and that it will initiate interdisciplinary cooperation and discussion on the more intensive and coordinated investigation of this complex and exceedingly important issue for Croatia, the Adriatic, and the Mediterranean region.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/jmse9040358/s1. The following Supplementary Materials are submitted alongside the manuscript: Figure S1a. The RAPS visualization of the mean annual SST time series; Figure S1b. The RAPS visualization of the mean annual air temperature time series; Table S2. Statistics (minimum, average, maximum, and range) of the mean monthly SST, air temperature (TA), and their difference (∆T = SST-T A ) from the time series measured at Split, Hvar, and Komiža stations. The negative values of ∆T were highlighted in red.

Data Availability Statement:
The data used within this study are the property of the Croatian Meteorological and Hydrological Service (DHMZ). Terms of use, data availability, and contact can be found at: https://klima.hr/razno/katalog_i_cjenikDHMZ.pdf.