Forty-five years of oceanographic and meteorological observations at a coastal station in the NW Mediterranean: a ground truth for satellite observations

Marine and atmospheric parameters, including temperature observations from surface to 80 m (at 6 depths) are measured since September 1973 on a higher-than-weekly frequency, at a coastal station 4 km offshore L’Estartit (Costa Brava; NW Mediterranean). This constitutes the longest available uninterrupted oceanographic time series in the Mediterranean Sea. The present contribution focuses on observed climatic trends in temperature (°C/year) of air (AT; 0.05), sea surface (SST; 0.03), sea at 80 m depth (S80T; 0.02) and sea level (SL; 3.1 mm/year) as well as comparison with trends estimated from coincident high-resolution satellite data. The trending evolution is not uniform across seasons, being significantly higher in spring for both AT and SST, while in autumn for S80T. Other climatological results are a stratification increase (0.02 °C/year in summer temperature difference between 20 m (S20T) and S80T), trends in summer conditions at sea (when S20T > 18 °C), estimated as 0.5 and 0.9 days/year for the starting day and period respectively, and a decreasing trend of nearly 2 days/year in the period of conditions favourable for marine evaporation (when AT < SST). This last trend may be related to the observed decrease of coastal precipitation in spring. The long-term consistency in the in situ SST measurements presents an opportunity to validate the multi-decadal trends. The good agreement for 2013–2018 (RMS 0.5–0.6, bias − 0.1 to − 0.2; trends of 0.09 °C/year in situ vs. 0.06 to 0.08 °C/year from satellite) allows considering this observational site as ground truth for satellite observations and a monitoring site for climate change.


Introduction
Long-term monitoring sites are scarce in the world ocean. Most of the financial sources have been traditionally addressed to achieve short-term oceanographic objectives. In comparison with weather stations, that were required to produce weather forecasts, there was absence of equivalent oceanographic forecasts. There was a general lack of acknowledgement that an oceanographic equivalent was necessary, possibly because the sea state depends essentially on winds, and marine traffic does not directly depend upon ocean temperature or salinity. Oceanographic information was known to be useful for fisheries, but does not have an immediate impact as to require high-frequency monitoring. Indeed, ocean motion is slow, and it was a widespread perception that ocean properties change smoothly. There was also clear evidence that ocean data acquisition costs were much higher than their atmospheric counterpart. Thus, early oceanographic efforts focused on producing maps of oceanic circulation and water masses. The underlying idea was that these properties could be considered as a permanent background, with some degree of variability (that is, a mean value with more or less variance).
Two of the first parameters that were regularly monitored were sea level and sea state, essentially because of their impact on shores, coastal equipment and harbours. Reconstructions of the oceanic history were derived from the comparison of data obtained from different oceanographic cruises. Although data was collected under diverse projects, comparisons evidenced that changes at sea were more rapid and much less smooth than foreseen. In most cases, such cruises had no common purpose nor belonged to any monitoring program.
By the end of the last century, the concept of Operational Oceanography emerged as an equivalent to the weather reports, with the associated requirement of regular monitoring. Other drivers fostering oceanographic monitoring have been the need to evaluate anthropogenic contributions to climate change. Although the relation between the ocean and climate is not new, it is recently vindicated in many papers with a classical introduction justifying the interest of ocean monitoring, satellite observations, modelling exercises or altogether. During the 1980s, when signs of a heating trend on climate started to be evident, international programs such as the World Ocean Circulation Experiment (WOCE 1988) started to foster long-term oceanographic studies. It was recognized "that successful simulation of the broad evolution of global climate beyond a few years will require the inclusion of the world oceans through their full depth". For the first time, a global pattern of oceanographic observations was established, with common methodologies that could be repeated-most of them have indeed been repeated several times-enabling interesting comparisons related to the impact ocean circulation on climate and vice-versa (e.g. Bryden et al. 2005;Ganachaud et al. 2000).
Nowadays, it is widely accepted that monitoring the ocean properties provides valuable information on climatological (long-term) mean values, variability ranges and long-term changes or trends. The recent growth of technological innovations fostered oceanic monitoring to unprecedented levels: from sea surface to ocean floor, from coast to open seas, by means of satellite, drifters, gliders, moorings, buoys, Argo profilers, ships of opportunity, as well as the classical oceanographic cruises. These networks are producing massive amounts of data that are then used, in near real-time or in data-assimilative models, to produce ocean state estimations and oceanic forecasts, and to give estimates of climate trends. This is, however, a very recent situation in terms of climate evolution, involving not more than the last 20-30 years. Forecasts and predictions face strong data constraints and challenges, including uncertainty quantification, due to the lack of long oceanic time series.
Indeed, long time series of ocean properties are rare in the surface ocean, and even rarer in the deep ocean. In days past, most of these initiatives took advantage of the presence of coastal laboratories that promoted fixed oceanographic stations easily reachable from their onshore installations. One example of fixed oceanographic stations is the California Cooperative Oceanic Fisheries Investigations (CalCOFI), formed in 1949 to study the ecological aspects of the sardine population collapse off California. Today, their focus has shifted to the study of the marine environment off the coast of California, the management of its living resources and monitoring the indicators of El Niño and climate change. CalCOFI conducts quarterly cruises off southern and central California, collecting a suite of hydrographic and biological data on station and while underway (Bograd et al. 2003). The North Atlantic region is an unusual example of long-term monitoring. It has a relatively large number of locations at which oceanographic data have been collected repeatedly for many years or decades; the longest records go back more than a century (e.g. ICES data sets). In the Mediterranean Sea, oceanographic data series during the last century are more scarce and unevenly distributed. One of the oldest fixed oceanographic stations is the observing point B in the Villefranche sur Mer bay (Ligurian Sea; France) that started weekly oceanographic data collection in May 1957 (Bougis and Carré 1960). Observations in this point have been maintained up to now although with several interruptions.
Likely due to the financial restrictions mentioned earlier, many of the atmosphere and oceanic long time series initiatives started (or were in part maintained) as a non-governmental, non-profit, voluntary initiative. One famous example is the oldest climate time series, the Mauna-Loa station, whose data was used as evidence for anthropogenic increase of the CO 2 concentration in the atmosphere. Similar origin had many initiatives collecting relevant oceanographic data. For example, the Oceanographic Society of Gipuzkoa has recorded sea surface temperature (SST) measurements in Donostia (southeastern Bay of Biscay; Spain), on a nearly daily basis since July 1946, providing more than 70 years of SST measurements to analyse local hydrographic trends and anomalies (Goikoetxea et al. 2009). The present paper refers to another of these voluntary initiatives: the Meteorological and Oceanographic coastal station at L'Estartit (Costa Brava, NW Mediterranean; Spain) that is carrying out weekly oceanographic observations uninterruptedly since September 1973. These data provide a unique opportunity for comparisons with satellite-derived SST in coastal regions. This paper aims to be not only a tribute to its origins and its scientific value, but also a wake-up call to start, today, similar initiatives and to insure their continuity through time.
The paper is outlined as follows. In section 2, we present the history of the L'Estartit station and the data sets produced, detailing the methods used to record each type of observation. In section 3, the methods to calculate climatic trends for the selected variables such as air, sea surface and sub-surface temperatures, changes in stratification, estimates of favourable conditions for evaporation (from the number of days during which sea surface temperature is higher than air temperature), sea level and precipitation. Methodology for Satellite SST products is also presented in this section. In section 4, we present the results obtained from the L'Estartit data. Section 5 places our results in a general context and are compared with previous works in the region and those based on L'Estartit data sets and compares the observed in situ SST with the available satellite products and their trends. Some general conclusions are outlined in section 6.

History and environmental considerations
As early as in 1969, Josep Pascual started a series of daily meteorological observations, measuring precipitation, air temperature and atmospheric pressure at L'Estartit, a small coastal town in the Catalan Costa Brava (NW Mediterranean). Later on, in September 1973, the inland meteorological observations were complemented with sea temperature observations at an ocean point located a few miles offshore (Fig. 1). L'Estartit is located at the north side of a small bay. It has a small harbour open to the SE and there are several small islands closing the bay. A hill, 228 m high, dominates the coastal cliff.
Ocean observations at a fixed point (AA) at about 4 km offshore L'Estartit, away from the bay and the islands, were obtained using calibrated reversing thermometers at 4 to 7 fixed levels, from surface to 80 m depth, with an average frequency of at least once per week, weather permitting, resulting in 60 to 80 visits per year, uninterrupted until today. The station, although quite close to shore, is located at the southern limit of the relatively narrow continental shelf, at the coastal side of the main along-slope current in the NW Mediterranean, the Northern Current, that flows in SW direction from the Ligurian Sea to near the Alboran Sea. Oceanographic conditions at this site are representative of those found upstream, at the southern limit of the continental shelf of the Gulf of Lions (Font et al. 1988).
Additional observations (mostly still on paper, not yet converted onto digital format) on marine environmental conditions at sight, such as water transparency (Secchi disc), clouds, sea-state, presence of birds, medusae or marine mammals and surface roughness and visual surface currents at several points along the track across the islands to reach the point AA were also recorded (Fig. 2). Inland meteorological observations of daily air temperature, precipitation, atmospheric pressure and humidity started in 1969. With time, these have been progressively complemented with potential daily evaporation (Piché) since 1976, and with sea state, wind speed and direction measurements since 1988. Since 1990, continuous analogic sea level data has been collected by a tide gauge installed inside the harbour (Fig. 3). Finally, since early 1990s, conductivity-temperature-depth (CTD) profiles along the track across the islands, and water samples for salinity analysis were also occasionally collected.
The zone near the station was declared marine-protected area (MPA) in 1990 (within the Parc Natural del Montgrí, les Illes Medes i el Baix Ter). The station is still maintained by Josep Pascual, as a voluntary observer, with the collaboration of the Institut de Ciències del Mar (CSIC) and the authority of the MPA. The series of observations is the longest uninterrupted within the Mediterranean Sea. Data (detailed in Table 1) are free and digitised data can be accessed upon demand (http://meteolestartit.cat/).

Data set
The main station (AA) is 90 m deep and located at~42°0 3′ N, 3°15′ E (Fig. 1). The average frequency of observations is higher than once per week (~60 to + 100 visits per year). Sea temperatures have been obtained using two Richter and Wiese protected reversing thermometers mounted on a Nansen bottle. Temperature readings, after correction according to the thermometer's calibration certificates, are averaged if the temperature difference among thermometers is lower than 0.2°C and repeated otherwise, so that the accuracy can be set to be lower than 0.2°C. Thermometers are periodically controlled at facilities of the Marine Research Institute (ICM) in Barcelona against a high-precision Pt 100 . which, in turn, is annually controlled using Ga and pure H 2 O melting points. Sea level gauge is a mechanical device attached to a continuous graphical recorder (Fig. 3). The record is digitized at 2-h interval and the originals preserved for detailed analyses if necessary. The position is periodically georeferenced by the Catalan Cartographic Institute (ICC) every 5 years and data is backwards linearly corrected for each period.
Atmospheric parameters (air temperature, pressure, relative humidity, precipitation and potential evaporation, Piché) have been measured in a meteorological box with standard instruments certified and regularly checked by the Spanish Meteorological Office (AEMET). These observations have been included in the official AEMET weather observation network. The meteorological box is located in front of the harbour at the sea level. Wind measurements at this location are not representative of the coastal zone because of the shelter from buildings and the coastal hills. For that reason, an anemometer is located on top of the hill, free of any sheltering. The anemometer is also controlled and certified by the Meteorological Office. Wind velocities and direction on top of the hill are very similar with those regularly measured offshore at the station AA (std error < 0.4 m/s). Wind velocities measured at the harbour station are Data sets of monthly-averaged data for air temperature, sea level and sea temperatures at surface, 20, 50 and 80 m, as well as monthly precipitation, up to December 2018, that have been used for the present paper, can be found and freely accessed through the PANGAEA repository (https://doi.org/10.1594/PANGAEA.902591) (Pascual and Salat 2019).

Trend estimations for in situ observations and derived parameters
For climate trend analyses, the primary measurements (all temperatures and sea level) have been monthly averaged. To remove the seasonal cycle, monthly anomalies are calculated by subtracting the monthly-averaged climatology (calculated over 45 years of data). The total number of data points for each trend estimation is thus 540. Trends are obtained by linear regression of the anomalies and converted into values-peryear (per annum). Significance of the trend estimates is evaluated using the "t test" compared against a "null hypothesis" of no trend (zero trend value). The confidence interval is computed based on the significance level α = 0.01 (99% confidence). The associated "p value", stating the probability of observing the given trend under the null hypothesis, quantifies how rare it is to observe such a trend only by chance.
Seasonal trends have been computed selecting January, February and March anomalies for the winter season, and so on for the rest of the seasons. Trends are always expressed in rates of variation per annum. The total number of data points for each seasonal trend is thus 135. Precipitation has been seasonally and annually accumulated so that the total number of precipitation data points is 45 either for annual and seasonal trends.
In temperate seas, the annual cycle of the water column show the formation and subsequent deepening of a seasonal thermocline from early spring to late summer. Vertical exchanges of water during the stratified period are almost restricted to the surface mixed layer above the thermocline. Storms and surface cooling in autumn and early winter enlarge surface mixed layer to deeper layers to reach the bottom by the end of the year. According to the careful analyses done by Coma et al. (2009) using the present data set (up to 2006), two different criteria were used to characterize the stratification conditions at the observing point: mean monthly temperature at 20 m (S20T) higher than 18°C and the difference  3 Levels marked in italics are those sampled at every visit, the rest are sampled with lower frequency between S20T and mean monthly temperatures at 80 m (S80T) is higher than 4.7°C. It is assumed that in spring, when the S20T reaches 18°C, stratification is consolidated and cannot be reversed until the autumn cooling. Similarly, limited vertical water exchanges will be reflected when higher than a threshold temperature difference is between 20 and 80 m. In the present paper, since there were almost no differences in the behaviour of these two criteria, we used only the first one for the duration of the stratified period, also referred to as summer conditions. The evolution of S20T-S80T will be used for climate trends. Evaporation over the sea cannot be easily measured and, in many cases, especially in modelling, bulk formulae have been used to get reasonable estimates of the evaporation rates. These estimates are essentially based on the capacity of an air mass to absorb water vapour assuming that air in contact with the sea surface will get the same temperature as the sea surface and reach the water saturation (Jacobs 1942). Although our data sets would provide information to have estimates of the sea evaporation, they would not be as consistent as other data reported here, so we decided not to include them in the results and consider only the evolution of the conditions favourable to an effective evaporation. As observed by Sverdrup et al. (1942, p 117-119), sea surface temperature higher than air temperature will favour evaporation while in the opposite case it will be greatly limited. Accordingly, within the annual cycle, there is a period favourable to sea evaporation, when SST is higher than air temperature (AT). In this paper, we used the monthly AT-SST to estimate the evaporation favourable conditions.
To estimate the duration of annual periods for a given condition such as that favourable to evaporation or the abovementioned stratification, linear interpolation between monthly data is used, assuming that the actual monthly mean value corresponded to the mid-month day (14 to 15.5 according to the month). Trends for these periods are obtained by linear regression with the same treatment as monthly or seasonal data but with only 45 data points, one per year.
Heat content within the upper 80 m of the water column (not shown) can also be calculated using the sea temperatures obtained at surface, 20 m, 50 m and 80 m, assuming these to be representative of temperatures at 0-10 m, 10-35 m, 35-65 m and 65-80 m, respectively.

Satellite products
Several satellite-based sea surface temperature (SST) analysis data sets were used for comparisons. These data sets are daily compositions of multiple satellite SST measurements of various horizontal resolutions ranging from 1 to 25 km, and they differ in selection of the input data sets as well as in how the input data are interpolated on to a grid whose resolution also varies from 1 to 25 km. They are the Multi-scale Ultra-high Resolution (MUR) analysis from NASA Jet Propulsion Laboratory (Chin et al. 2017), Italian Institute of Atmospheric Sciences and Climate G r u p p o d i O c e a n o g r a f i a d a S a t e l l i t e ( G O S ) Mediterranean analysis (Buongiorno Nardelli et al. 2013), the Canadian Meteorological Centre (CMC) analysis (Brasnett 2008), UK Met Office Operational Sea Surface Temperature and Sea Ice Analysis (OSTIA) (Donlon et al. 2012), and National Center for Environmental Information (NCEI) optimal interpolation analysis (Reynolds et al. 2007). All these data sets may be a c c e s s e d t h r o u g h r e g i s t r a t i o n a t t h e P h y s i c a l Oceanography Distributed Active Archive Center (PO.DAAC) and use of the PO:DAAC drive at https:// podaac-tools.jpl.nasa.gov/drive/. For comparison, an SST time series coincident with the in situ counterpart is obtained from each of these analysis data sets by bilinear interpolation to the L'Estartit sampling location on the sampling dates. Since the earliest satellite SST measurements are from 1980s and some of the satellite analyses go back in time much less, we make comparison only for the recent years, from 2013 to 2018.

Annual cycles
The average annual cycle of air and SST (Fig. 4a) shows that from late winter-early spring to mid-summer, SST is lower than air temperature while during the rest of the year, the situation is reversed. The initial and ending days for the period when SST are colder than air temperature is presented in Fig.  4b and the corresponding annual trend in Table 2. There is a successive advancement of the starting point of the period (p < 0.0001) while the delay in ending point is less significant (p < 0.01, i.e., still within the 99% of confidence). Such a behaviour has a relevant role in the air-sea interaction, in particular concerning sea water evaporation/condensation, and thus in the local water cycle.
The annual cycle of S20T is presented in Fig. 5a to indicate the initial and ending points of the stratified period. The evolution of these points is presented in Fig. 5b and the corresponding annual trends in Table 2. The evolution shows an enlargement of the period associated to both an advancement of the starting point and delay of the ending point. Although there are strong interannual differences in this process, from cold years when stratified period is short to hot waves with long periods of stratification, all these situations have significant biological implications that may affect local fisheries.

Monthly temperatures
Monthly AT, SST and S80T temperature anomalies, with the respective linear trends are presented in Fig. 6. All the three sequences show a high interannual variability and upward linear increasing annual rates (Table 2). AT a warming trend higher than SST and S80T is increasing at a lower rate than at surface. Seasonal trends at surface show their maxima in spring (p < 0.0001) and a less marked and also less significant value (p < 0.01) in summer and autumn (Table 2).
Monthly temperature differences between air and sea surface (AT-SST) are presented in Fig. 7a. According to the trends shown above, there is a significant upward annual trend. Seasonal trends are lower or similar than the annual trend, except in spring where it is nearly twice the annual and more than three times the winter trend (Table 2). Monthly temperature differences between sea temperatures at 20 and 80 m (S20T- Fig. 4 a Mean annual cycle (1974-2018) of monthly air and sea surface temperature. The period when air temperature is higher than sea surface temperature is marked. b Time evolution of the period when air temperature is higher than sea surface temperature. Blue points indicate the initial day of this period and pink points the end. Lines show the linear fitting S80T) are presented in Fig. 7b. In this case, the summer seasonal trend is about twice the annual trend-other seasons show lower or non-significant trends (Table 2).

Monthly sea level
The amplitude of the astronomical tide at L'Estartit is lower than 20 cm. Monthly averaged mean sea level as well as monthly maximum and minimum anomalies display a mean positive trend of more than 3 cm/decade ( Fig. 8; Table 2). Sea level records reflect some episodic seiches with high-frequency oscillations (1/600) and tide amplitudes reaching 100 cm, associated with shifts (> 2 mb/h) and rapid oscillations of the atmospheric pressure (Fig. 9). The positive trend in sea level would be at least partially compatible with the net heat gain shown above.

Local precipitation
Annual and Spring accumulated precipitation are presented in Fig. 10, and the corresponding trends appear in Table 2. Although anyone of these trends is not statistically significant (only in spring has a relatively low p value), they are included for its descriptive interest. In particular, it is worth to note that almost all the annual decay corresponds to the spring season. This trend is also confirmed by the decrease in spring precipitation of nearly a 30% by comparing the average from 1974 to 1983 and 2009-2018. It also can be related to the decay in sea evaporation in the season suggested by the trend to early dates of starting the period in which sea surface temperature is lower than air temperature ( Fig. 4; Table 2).

In situ data vs. current satellite products
The L'Estartit sea surface temperature (SST) time series agrees quite well with the SST time series derived from the   various satellite-based data sets (Fig. 11), with an RMS difference of consistently just over 0.5°C (Table 3 and Table 4). This is comparable with the typical minimum uncertainty in satellite SST measurement of 0.3-0.5°C. The satellite SST is cooler than the L'Estartit SST by 0.1 to 0.2°C on average. A potential cause for this bias is the diurnal warming of sea temperature: while the satellite data sets are aiming to report the pre-dawn temperature, the L'Estartit measurements are performed mostly in the late morning to noon. We note that the feature resolution (smoothness of the SST field) of the satellite data sets used here varies vastly from around 1 up to 100 km, with the OSTIA and NCEI data sets having a similar feature resolution of around 100 km despite having fairly different grid spacing (Reynolds and Chelton 2010;Chin et al. 2017). The comparison statistics (RMS and bias) are remarkably consistent across feature resolutions up to 25 km, and degrades only by 0.1 to 0.2°C for the OSTIA and NCEI data sets. This indicates that the temporal variability of the L'Estartit time series is representative of the areal average SST from the surrounding sea. The good statistical comparison of sea surface temperature values (Table 3 and Table 4) and trends (Table 5) obtained Fig. 10 Time series of annual (blue) and spring (red) precipitation and its corresponding linear trend (significant at 85%) Fig. 9 Example of high-frequency sea level oscillations (seixes) recorded on 2 June 1994, starting at 06:30 GMT. Enclosed is the atmospheric pressure record of the first 12 h of the same day from L'Estartit data set and from available satellite products confirms that the reported sea surface temperature trends are robust. Additionally, their good agreement for such a longterm (more than 5 consecutive years) is a first step towards the cross-validation of these products. That opportunity to validate the multi-decadal trends that can now be observed in satellite data gives more confidence to use these data to model and/or parameterize time changes in oceanic states (which have often assumed to be recurring in order to develop theories, e.g. for spatial patterns).

Comparison with previous and other estimates of climate trends
Most of the articles published within the two last decades show the Western Mediterranean has increased in temperature and salinity during the twentieth century. Most of these articles, however, focused on the intermediate and deep layers; works dealing with changes in the upper layers and shallow waters remain scarce. Vargas-Yáñez et al. (2005) analysed coastal data obtained at several fixed stations and detected strong warming trends in the near surface waters, usually ranging from 0.02 to 0.05°C/year, up to almost 0.1°C/year (in the Ligurian sea) during the 1990s. These results (of almost one order of magnitude higher than those found at intermediate and deep waters) contrasted with those from Krahmann and Schott (1998) and Sparnocchia et al. (1994) that found no significant changes in the temperature evolution of the upper layer in the Western Mediterranean since the mid-1970s. More recent works, that included data from satellite infrared data, evidenced the warming of the surface layer, both in the Western and Eastern Mediterranean since the mid-1980s to the first decade of the twenty-first century (Skliris et al. 2012;Nykjaer 2009). The trend estimates in these two papers are close to those obtained from our station. Similarly, Rixen et al. (2005), using data from the MEDATLAS database, found positive trends of about 0.02°C/year in the upper layer temperatures since mid-1980s, after an almost steady (zero trend) situation in the previous 30 years. Near the coast in   Vargas Yáñez et al. (2009 attempted to estimate the evolution of air and sea temperatures during the whole twentieth century, up to 2008, in the Spanish waters of the Western Mediterranean. The authors merged the information from MEDATLAS with some of the more recent time series at fixed stations (RADMED) in the region. The estimates for the whole period and the whole area showed a non-significant trend for the upper layer of sea temperature of 0.0038 ± 0.0032°C/year (partially due to the many gaps in the available series) and a significant trend of 0.0074 ± 0.0013°C for the air temperature. The evolution, however, showed an acceleration of the warming trends since around 1975 (from when data started to be more abundant although many locations still presented gaps). Trends similar to ours were only detected in the Balearic Sea, while much lower trends were found in the rest of NW Mediterranean areas. These analyses were updated in Vargas-Yáñez et al. (2017), and described the results for the period 1943-2015 in four separate areas, from the Alboran Sea to near the Gulf of Lions. Again, the upper layer temperatures did not show any significant trend, but while in the Alboran and Balearic areas, these trends were slightly negative, in the SE of Spain and in the north trends were slightly positive.
The Copernicus Marine Environment Monitoring Service (CMEMS) provided annual reports of the state of the global ocean and European regional seas based on all the available observations from 1993. In the 2016 report (von Schuckmann et al. 2016), the observed 1993-2015 SST trends for the whole ocean indicate a global increase of 0.016 ± 0.002°C/year (99% confidence interval), that was quite sensitive to the particular strong increase in SST during the 2015 El Niño event. The Arctic region also showed a high SST increase during 2015 while this signal was much weaker in the North Atlantic. Within the European regional seas, the Mediterranean (including the Black sea) is where trends are the highest (0.039 ± 0.009°C/year (95% confidence interval). However, the Western Mediterranean basin presented lower trends, in good agreement with other observations (e.g. Vargas-Yáñez et al. 2009 and. Again, the relative maxima of the NW basin appeared in open sea while at the coastal zone, including L'Estartit, remained within 0.02-0.03°C/year. The CMEMS paper also mentioned that trends for subsurface layers down to 100 m are dependent on SST. For that reason, there is no information about the evolution at 80 m. The update of the CMEMS report for 2018 (von Schuckmann et al. 2018) enlarged the length of the time series to 2016 and obtained a slightly higher trend for SST with reduced uncertainty 0.040 ± 0.004°C/year. This may be attributed to the increase of the time series and the similarly high temperatures in the last years of the series.
In terms of sea level rise, despite the global increasing trend, sea level in the Mediterranean sea has been decreasing from the 1960s to the 1990s, at an average rate of − 0.5 to − 1 mm/year, reaching rates as low as − 1.3 mm/year in some locations (Marcos and Tsimplis 2008). Such negative trend was attributed to an increase of the atmospheric pressure in the Mediterranean sea region, associated to the positive phase of the North Atlantic Oscillation (Gomis et al. 2008). Although Tsimplis and Baker (2000) suggested the increase would be associated to the increasing salinities of the sea during that period, Jordà and Gomis (2013) demonstrated that the influence of the reported changes in salinity was only residual. Since the mid-1990s, the atmospheric pressure recovered to its normal values for the Mediterranean region and the sea level started to rise at trends higher than 2 mm/ year, reaching up to 10 mm/year (Criado-Aldeanueva et al. 2008;Gomis et al. 2008). Marcos et al. (2016), using satellite altimetry data, estimated an average increase in absolute (geocentric) sea level of 2.6 ± 0.2 mm/year for the period 1993-2015.
The CMEMS report (von Schuckmann et al. 2016) estimated the sea level rising trend for 1993-2015 in the whole Mediterranean was 2.9 ± 0.9 mm/year (95% confidence interval), of which about one half corresponded to the thermosteric component. This figure is very similar to ours (Table 2). Interestingly, the whole-Mediterranean time series of sea level anomalies reproduced most of the peaks observed in our time series (Fig. 8)-namely in 1996, 1997, 2010 and 2013, but not  (2016) slightly reduced the previous estimate for the whole Mediterranean to 2.7 ± 0.9 mm/ year, as for other regional seas. The estimate for L'Estartit zone lies within 2.5 and 3.5 mm/year, in good agreement with CMEMS.
In terms of specific aspects of the annual cycle, such as the seasonal distribution of trends, the evolution of the difference between sea surface and air temperatures or stratification periods dealt in our paper, we hardly find any paper that specifically addressed these aspects in the Mediterranean, with the exception of those that used this set of data (e.g. Coma et al. 2009). Only scattered references to some of these subjects in the above-mentioned papers. In particular, it is worth to mention that in Skliris et al. (2012) seasonal trends were estimated for the Western, Eastern and the whole Mediterranean sea. In the Western Mediterranean, these trends were found to be much higher and significant in spring and summer than it was found in our series. The same paper also highlights the role of the latent heat exchanges as the dominant factor among the drivers of seasonality and long-term trends. In particular, they found that the temperature difference between air and sea surface is very well correlated with the latent heat exchanges, whose annual maximum falls in winter, as widely assumed. The relation of these temperature differences with latent heat exchanges, hence with evaporation, is a good support for our hypothesis. In particular, the decrease in spring evaporation may be related to the decrease in local precipitation during this season. Although we did not find any paper that highlighted this relation, many papers on Mediterranean climate (e.g. De Luis et al. 2009Luis et al. , 2010 show a decreasing trend in annual precipitation in the Western Mediterranean region that is more larger in spring. Specific events such as the heat waves of summer 2003 and 2006 have been clearly reflected in the series and, although it is a shallow coastal station, the mild to very mild 1990, 2014, 2016 and 2017 winters as well as the exceptional heat loss in the 2004-2005 winter. While this extreme heat loss conducted to an unprecedent volume of dense water formation in the Western Mediterranean (Canals et al. 2006;Font et al. 2007), the previously mentioned mild winters only produced small amounts of water not dense enough to get the deep layers (Mertens and Schott, 1998;Schoreder et al. 2017).
Discrepancies among results reported above, as well as with ours, can be attributed to different factors, including uncertainties resulting from merging data from wide areas, the use of different methodologies to calculate the trends and the collection of data from different seasons. Last but not least, the high variability in the upper ocean layer introduces large standard deviations in data. The good statistical comparison of sea surface temperature values and trends obtained from L'Estartit data set and from available satellite products is encouraging.

Conclusions
The Mediterranean Sea is a recognized hot-spot for future climate changes, due to its characteristics of a miniature ocean, since most oceanographic processes of the global ocean (with the exception of sea ice formation processes) occur in it. More importantly, its relatively small size makes it more accessible to monitoring than other regions of the world. According to its geographic position in mid latitudes, the Mediterranean is subjected to a marked seasonality that essentially affects heat exchanges with the atmosphere and ocean stratification. The set of observations presented in this contribution, not yet completely digitalized, constitute a unique reference document of the evolution of a Mediterranean coastal environment for about half a century. Among other uses, it permits one to derive trends for local climate evolution, from surface to stratification of the upper ocean. These trends can then be contrasted with those currently obtained through big data processes (satellite observations and numerical models). During the last years, the dataset has been used in many scientific papers related to climate research and to document extreme events such as cold or heat waves, storms and their consequences on neritic ecosystems. Robust estimates of annual trends of sea temperature from SST (0.03°C/year), S80T (0.02°C/year), AT (0.05°C/year) and sea level (3.1 mm/year) are presented. Seasonal trends and other properties related to the seasonal cycle are also presented, some of them for the first time. A careful inspection and comparisons among annual and seasonal trends highlight the relative importance of seasonal processes such the thermocline formation or destruction and conditions for air/sea heat and mass exchanges. For instance, a separation among impacts of climate evolution according to the seasons has been useful to explain several ecological impacts at a local scale. It can also be to assign climate change impacts to processes acting at different seasons that can be masked if only annual trends are taking into account. Finally, the results can be used as reference for other Mediterranean areas. Some relevant results associated to the annual cycles are the stretching of the period of positive difference among air and sea surface temperature, that shows an increasing trend of nearly 2 days/year. Such a trend may partially explain a 30% decrease (although not significant) of coastal precipitation in spring. The present contribution also highlights the potential of this coastal dataset for climate evolution studies and as ground truth for sea surface temperature (SST) satellite estimations. Although there have been comparisons with ARGO, etc., they have been over very limited periods of time. For instance, the comparison presented here extended along more than 16 consecutive years (Fig. 11) and it is up to our knowledge the longest ever done. The good statistical agreement between on-site SST measurements and SST obtained from satellite observations confirms the reported positive trends are robust. In views of the results presented here, we propose L'Estartit observational site should become a long-term ground truth for satellite observations and a regional monitoring site for climate change.