A near-surface sea temperature time series from Trieste, northern Adriatic Sea (1899–2015)

A time series of near-surface sea temperature was built from observations performed in the harbour of Trieste from 14 July 1899 to 31 December 2015. The description of the observation sites and instruments was possible thanks to historical documents. The measurements consist of two data sets: the first consists of analogue data obtained by means of thermometer and thermograph measurements, recorded one or two times per day, in the periods 1899–1923 and 1934–2008; the second consists of digital records obtained by thermistors on an hourly basis in the period 1986–2015. A quasi-homogeneous time series of daily temperatures at 2 m depth is formed from direct observations at that depth and from temperatures estimated from observations at shallower depths. From this time series a mean temperature rise rate of 1.1±0.3 ^∘C per century was estimated, while in 1946–2015 it is 1.3±0.5 ^∘C per century. The data are available through SEANOE (https://doi.org/10.17882/58728; Raicich and Colucci, 2019).

1 Introduction

Knowledge of the processes and evolution of the earth system depends on observations. As the oceans cover most of the Earth's surface, marine data are particularly valuable and, therefore, are in high demand by public authorities, private enterprises, the scientific community and civil society, in a wide range of research fields and applications. The public release of marine data and related products has been fostered in the framework of several European initiatives, such as, for instance, the European Marine Observation and Data Network (EMODNet) (http://www.emodnet.eu, last access: 19 May 2019) and SeaDataNet (https://www.seadatanet.org/, last access: 19 May 2019) (Shepherd, 2018).

Long time series of observations represent a key element in climate studies, based on both data analysis and model reconstructions, and thus the sustainability of suitable observing systems is critical, as emphasised by IFSOO (2012). Although sea temperature measurements date back to over 150 years ago (International Comprehensive Ocean-Atmosphere Data Set, Woodruff et al., 2011), time series longer than a few decades are uncommon, even in coastal areas.

Thanks to data retrieved from the archives of the Institute of Marine Sciences of CNR in Trieste, a quasi-homogeneous time series of daily temperature at 2 m depth was formed from near-surface sea temperatures measured in Trieste harbour, at the northernmost end of the Adriatic Sea, from 1899 to 2015. Even though these data are coastal, and therefore sensitive to local natural and anthropogenic processes, the considerable length of the time series makes it not only useful for local studies but also a long-term indicator of climate change in the northern Mediterranean.

The observation site is located in the Gulf of Trieste, which is a bay of about 20 × 25 km² with a maximum depth of 25 m (Fig. 1). The shallowness and semi-enclosed character of the basin makes the near-surface temperature strongly dependent on the atmospheric forcing. As a result, a marked seasonal cycle is observed, ranging from 9 ^∘C in February to 25 ^∘C in July (period 1991–2003; Malačič et al., 2006); moreover, fast temperature changes often occur on synoptic timescales, due to bursts of north-easterly Bora that causes coastal upwelling (e.g. Crisciani and Raicich, 2004; Crise et al., 2006) and intense air–sea heat fluxes (Stravisi and Crisciani, 1986; Picco, 1991; Supić and Orlić, 1999; Malačič and Petelin, 2001; Raicich et al., 2013). No significant influence on sea temperature is expected from the freshwater discharge in the Gulf of Trieste. The river runoff is concentrated in the northern part of the gulf, about 20 km from Trieste, where the mean annual discharge is 114 m³ s⁻¹, i.e. over 95 % of the total, while only 5 m³ s⁻¹ discharge occurs from small streams along the Slovenian coast, at about 10 km from Trieste (Cozzi et al., 2012).

Figure 1 Aerial image of the area of Trieste harbour with the observation sites: Molo Sartorio (SAR), Molo Santa Teresa (TER), Bagno Savoia (SAV), Porto Lido (LID) and Molo Bandiera (BAN). (Images extracted from Google Earth; ©2018 Landsat/Copernicus, ©2018 CNES/Airbus, ©2018 Digital Globe, ©2018 TerraMetric.)

In the next section the data used in this work and their sources will be described. The methods used to derive the 2 m temperature time series and to estimate the related errors will be outlined in Sect. 3. Section 4 will include basic information on the data availability. Concluding remarks will be presented in Sect. 5.

2 Data description

2.1 Measurement sites and instruments

All the sites are located in Trieste harbour (Fig. 1) within an area of approximately 0.1 km²; their geographical positions are summarised in Table 1. Molo (Pier) Sartorio (SAR) and Molo Santa Teresa (TER) face a semi-enclosed bay, while Bagno (bathing establishment) Savoia (SAV), Porto (harbour) Lido (LID) and Molo Bandiera (BAN) face the open sea. For each site, Table 1 summarises the depths of measurement, the instruments used and the periods of observation.

Table 1 Summary of observation sites, their geographical position and instruments (b.m.s.l.: below mean sea level).

The earliest regular observations of sea temperature date back to summer 1899. The first instrument used was a Pinsel-Thermometer, literally translated as brush thermometer, manufactured by the H. Kappeller firm in Vienna, with a scale pitch of 0.2 ^∘C. In order to measure temperature, the instrument was manually deployed to the prescribed depth and kept there until thermal equilibrium was reached. Its bulb was wrapped in a sort of brush which could absorb water and keep it insulated from the surrounding environment while the thermometer was being extracted from the water and read (Attlmayr et al., 1883).

Figure 2 Weekly thermograph chart from Bagno Savoia for 1–8 July 1935. Values written above the curve are temperatures at 12:00 noon. A remarkable cooling from 25 to 18 ^∘C, due to Bora-induced upwelling, is evident on 6–7 July.

On 8 January 1934 a thermograph, manufactured by Paul Wegener GmbH in Ballenstedt, Germany, was put into operation. The instrument was connected by a 5 m long capillary to a bulb placed at 2 m below the mean sea level at about 70 m from the shore, with chart (fixed on a rotating drum) speed of 45 mm d⁻¹ (∼1.9 mm h⁻¹) and non-linear vertical scale corresponding to 2.1 mm ^∘C⁻¹ at 15 ^∘C; the chart was changed weekly. Figure 2 shows the chart for 1–8 July 1935. On 19 June 1943 the observations ended because the connection between bulb and capillary was broken and the bulb was lost; attempts to obtain a replacement failed probably because of the ongoing Second World War (Geophysical Institute, 1943). Another thermograph was employed from 27 April 1948 to 1 December 1952 (Picotti, 1955). It was manufactured by the SIAP firm in Bologna, Italy, with chart speed of 40 mm d⁻¹ (∼1.7 mm h⁻¹) and linear vertical scale corresponding to 1.5 mm ^∘C⁻¹.

While manual measurements were performed at prescribed depths below the sea surface, the thermograph bulbs were at fixed heights therefore their depths would change according to sea level variations. In a few cases a very low sea level caused the bulb to be close to the sea surface, thus affecting the measurements. Trieste harbour is subject to an astronomical tide range of 1.2 m on syzygies and a mean seasonal sea level range of approximately 15 cm, mainly modulated by the atmospheric forcing (pressure and dominant winds) and the steric effect. On synoptic timescales the atmospheric forcing has produced positive sea level surges (storm surges) up to 1.7 m above the mean sea level and negative surges up to 1 m below (Beretta et al., 2005; Tsimplis et al., 2009; Raicich, 2010).

The effect on temperature would consist of an anomalous increase due to solar heating of the surface layer, observed at low tide and in weak wind conditions, particularly during late spring and summer sunny days. We know that in the 1934–1943 period the curve was usually interpolated or smoothed manually to correct the anomalous temperature values (Polli, 1946). Although no explicit information is available for the 1948–1952 period, the same procedure was also likely adopted.

Well thermometers have been usually employed for direct measurements, generally to fill gaps in the thermograph records. The instruments were manufactured by SIAP and their scale pitch was 0.1 ^∘C.

In 1986, automatic temperature recording started at Molo Bandiera using PT100 thermistors having a nominal accuracy of approximately 0.1 ^∘C. Usually, ocean temperature is reported to 0.01 ^∘C precision, i.e. 1 order of magnitude better; nevertheless, the PT100 accuracy is considered sufficient to measure the highly variable near-surface sea temperature in Trieste harbour, which is strongly influenced by atmospheric forcing.

We will distinguish two main data sets. Set A consists of manual and analogue automatic measurements performed from 14 July 1899 to 31 December 2008. Set B consists of digital records from 18 November 1986 to 31 December 2015. Several gaps exist due to instrumental malfunctions and failures, particularly in the earlier years of each data set, and some missing data may have simply been lost.

2.2 Data set A (1899–2008)

The Maritime Observatory (Osservatorio Marittimo) of Trieste managed the earliest sea temperature observations from 14 July 1899 to 31 December 1920. They were either reported as additional columns on tables of hourly air temperature data, obtained from the digitisation of thermograph charts, or on separate sheets. The sheet for 1916 is shown in Fig. 3. The data for 1905, 1909–1914 and 1918–1919 are missing. According to archive documents, the observations were also performed at least in the second half of 1918 and the first half of 1919 (Maritime Observatory, 1919), but the data could not be found. No information is available for the other data gaps.

Figure 3 Sea temperature summary sheet for 1916. On top, the indication of site and observation time: the temperature was measured at Trieste lighthouse in degrees Celsius at 08:00.

The data sheet for July 1899 reports that the measurements were performed at the head of Molo Sartorio at 0.75 m depth at 09:00 (Central European Time). This was the situation until 1904, while in 1906–1908 the time changed to 10:30. In 1915–1917 the site was moved by approximately 300 m to Molo Santa Teresa near the lighthouse and the observation time was 08:00. In 1920, the temperature was measured at 09:00, again at Molo Sartorio.

We can only assume that the nominal depth remained 0.75 m, as after the very first data sheet of 1899 no further information was provided. This assumption seems reasonable because the instrument was only suitable for near-surface water measurements. In 1915–1917 the data were also published in the daily telegraphic weather bulletin (Bollettino Meteorologico Telegrafico) issued by the Maritime Observatory (Maritime Observatory, 1915–1917).

No data are available for the 1921–1933 period except for a few sparse measurements in 1922–1923.

In 1934 the observations were resumed by the Geophysical Institute (Istituto Geofisico), which was renamed the Thalassographic Institute (Istituto Talassografico) in 1941. The temperature at 12:00 noon was usually extracted from the thermograph records for publication. Manual measurements with a well thermometer were performed in case of thermograph failure; nevertheless, several gaps still affect the record, particularly during 1941–1944.

Manual observations continued until 2008 and daily values were reported in monthly summaries alongside the meteorological data. From 1944 to 1979 they were performed at Bagno Savoia at 08:00 and 19:00; in 1952–1957 an additional observation was added at 12:00 noon. From 1964 to 1968 temperature was also measured at Molo Sartorio in the morning (times were reported). On an unknown date between 1971 and 1975 the site was moved to Porto Lido, 180 m from Bagno Savoia (Fig. 1). From 1980 to 2008 the observation was nominally made at 10:00.

The measurement sites can be divided in two groups. The first consists of Molo Sartorio and Molo Santa Teresa, located inside the Sacchetta Basin (Fig. 1); the second is composed of Bagno Savoia, Porto Lido and Molo Bandiera, outside the Sacchetta, facing the open sea. Almost simultaneous (within half an hour) observations allow the comparisons of Molo Sartorio with Bagno Savoia (using 1295 “surface” data from 1964 to 1968) and with Porto Lido (using 299 data at 2 m depth from 1983 to 1986). The average temperature differences are negligible in both cases (0.0±0.6 ^∘C) and it is therefore possible to merge the different data records into a quasi-homogeneous time series.

Regarding the times of observation, it is certain that permanent staff were on duty in offices close to the observation sites at Molo Sartorio (in 1899–1908) and at Molo Santa Teresa (in 1915–1917 and probably in 1920), and thus deviations from nominal times were probably a few minutes. By contrast, from 1934 onwards direct measurements required the observer to move from his office to the observation site; therefore, the observation time was variable and sometimes weather conditions were so bad that they prevented the observations from being performed. The actual time was reported in many cases, but not always.

These analogue records were reported to 0.1 ^∘C precision and occasionally to 0.05 ^∘C.

2.3 Data set B (1986–2015)

From 1986 to 1989 a PT100 probe was used at 0.4 m depth, which can be deemed equivalent to the “surface” depth of data set A (Table 1). Observations were interrupted until 1992, when the depth was changed to 2 m and in 1994 temperature began being measured at both 0.4 and 2 m depths. In 1999 three probes were installed at 0.4, 2 and 6 m (sea bottom) depths; this configuration has since been kept. These probes were attached to a float in order to keep their depth constant.

The original data record consists of hourly temperatures, represented by 10 min averages at the end of the hour; from mid-2008 onwards, average, minimum and maximum temperatures were recorded every 10 min. The record is affected by several interruptions due to system failures; sometimes the probes were damaged or torn away by waves during storms. The measurements were originally reported to 0.01 ^∘C precision.

3 The time series at 2 m depth

A composite long-term time series of mean daily temperatures at 2 m depth was built merging the two data sets described in the previous section. When not available from direct measurements, the 2 m temperature was estimated from the observations at other depths, as explained in Sect. 3.1.

3.1 The estimate of mean daily temperature

A provisional daily mean was first obtained by averaging the original observations for each calendar day, even when only one observation was available. Although the mean daily temperature range is generally less than 0.5 ^∘C, a daily mean computed with few observations may significantly differ from a mean accounting for the full daily cycle, which is well represented by 24-hourly observations. To account for the bias, mean corrections were estimated using data set B (see Sect. 2.3).

For each depth (0.4, 2 and 6 m) and each calendar day (1 January–31 December), climatological values for the 1999–2015 period were obtained by averaging hourly (0–23) temperatures and mean daily temperatures of days when all the 24-hourly observations were available, in order to account for the complete daily cycle.

If h represents the hour, d the day, m the month, y the year and z the depth, let $T_{\mathrm{o}}(h,d,m,y,z)$ be the observed temperature, $T_{\mathrm{c}}(h,d,m,z)$ the climatological hourly temperature and $T_{\mathrm{24}\mathrm{c}}(d,m,z)$ the climatological mean daily temperature:

$$\begin{array}{ll}{\displaystyle }& {\displaystyle }{T}_{\mathrm{c}}\left(h,d,m,z\right)=\sum _{y=\mathrm{1999}}^{\mathrm{2015}}[T_{\mathrm{o}}\left(h,d,m,y,z\right)\\ \text{(1)}& {\displaystyle }& {\displaystyle }\cdot p\left(h,d,m,y,z\right)]/\sum _{y=\mathrm{1999}}^{\mathrm{2015}}p\left(h,d,m,y,z\right),\text{(2)}& {\displaystyle }& {\displaystyle }{T}_{\mathrm{24}\mathrm{c}}\left(d,m,z\right)=\sum _{h=\mathrm{0}}^{\mathrm{23}}{T}_{\mathrm{c}}\left(h,d,m,z\right)/\mathrm{24},\end{array}$$

p being a weighing factor, equal to 1 if T_o is available and 0 if it is not. For a given $(h,d,m,z)$, from 1999 to 2015, there are generally between 13 and 17 values of T_o. A 31-day running mean is subsequently applied to T_c and T_24c in order to smooth out the effect of outliers. The values for 0.75 m depth are obtained by linear interpolation between those of 0.4 and 2 m depths, and the values for 29 February are interpolated using those of 28 February and 1 March.

The mean daily temperature $T(d,m,y,z)$ is estimated as follows:

$$\begin{array}{ll}{\displaystyle }& {\displaystyle }T\left(d,m,y,z\right)=T_{\mathrm{24}\mathrm{c}}\left(d,m,z\right)+{\displaystyle \frac{\mathrm{1}}{N}}\sum _{k=\mathrm{1}}^N[T_{\mathrm{o}}\left(t_k,d,m,y,z\right)\\ \text{(3)}& {\displaystyle }& {\displaystyle }-T_{\mathrm{c}}\left(t_k,d,m,z\right)],\end{array}$$

where $N=N(d,m,y,z)$ is the number of available observations on the relevant day and t_k is the time of observation; when T_c is not a full hour, it is interpolated at the proper time using the nearest hourly values. The term in square brackets of Eq. (3) represents the departure of the observed value from climatology. Note that if T_o is available at full hour from 0 to 23, then T is the arithmetic average of the 24 observations.

In other words, the difference between an observation and the corresponding climatological value provides a constant that is used to rescale the climatological daily cycle; then, the average of the rescaled daily cycle represents the estimated daily temperature. If more observations are available in a given day, the final estimate is the average of the individual estimated values.

When the 2 m observations are unavailable, the daily means are estimated from data at other depths when possible. If Z=2 m, then

$$\begin{array}{ll}{\displaystyle }& {\displaystyle }T\left(d,m,y,Z\right)=T\left(d,m,y,z\right)-[T_{\mathrm{24}\mathrm{c}}\left(d,m,z\right)\\ \text{(4)}& {\displaystyle }& {\displaystyle }-T_{\mathrm{24}\mathrm{c}}\left(d,m,Z\right)],\end{array}$$

and therefore

$$\begin{array}{ll}{\displaystyle }& {\displaystyle }T\left(d,m,y,Z\right)=T_{\mathrm{24}\mathrm{c}}\left(d,m,Z\right)+{\displaystyle \frac{\mathrm{1}}{N}}\sum _{k=\mathrm{1}}^N[T_{\mathrm{o}}\left(t_k,d,m,y,z\right)\\ \text{(5)}& {\displaystyle }& {\displaystyle }-T_{\mathrm{c}}\left(t_k,d,m,z\right)].\end{array}$$

The term in square brackets on the right-hand side of Eq. (4) represents the correction that accounts for the incomplete representation of the daily cycle. Considering that, overall, the observations were made in the 08:00–19:00 interval (Sect. 2.2, above), the largest corrections, in absolute values, occur in June and July when values between −0.3 and 0.4 ^∘C are expected, depending on the time of observation, with an average of +0.1 ^∘C. The smallest corrections are found between mid-October and mid-February with values between −0.1 and +0.1 ^∘C and average around 0.0 ^∘C. These corrections are associated with estimates from single observations, as in data set A, while they are negligible in data set B where the complete daily cycle is generally available.

The term in square brackets on the right-hand side of Eq. (5) represents the correction related to the daily temperature estimate from data at depths other that 2 m. In this case, the corrections range between 0.02 ^∘C in late November to early December and 0.68 ^∘C in late June to early July. Larger corrections are a consequence of the water column stratification which typically occurs in summer.

3.2 The error in mean daily temperature

The error in T, namely σ, is computed from the errors in the observation, σ_o, and in the climatologies, σ_24c and σ_c, respectively. These errors are assessed semi-empirically as follows.

An observation is affected by an instrumental error and an environmental error. The error in temperature measured with manual thermometers can be estimated in half the scale pitch, namely 0.05 or 0.1 ^∘C, depending on the instrument. In the case of thermographs, the distance of 1.5–2 mm between the chart markings and the curve thickness determines a reading error of approximately 0.2 ^∘C. The environmental error is caused by small temperature fluctuations occurring in the water body at sub-hourly frequencies, due to turbulence and circulation. Therefore, an instantaneous measurement may not be representative of the hourly value. In order to estimate this error, average hourly temperature ranges are computed using the temperature extremes that are available for 2008–2015. The mean 1 h temperature ranges turn out to be 0.24, 0.30 and 0.29 ^∘C at 0.4, 2 and 6 m, respectively. A ±0.15 ^∘C error (0.3 ^∘C band) can be adopted for all depths.

An additional error often affects the observations of data set A because only nominal times and not the actual observation times are reported. This time uncertainty affects temperature because of the involvement of T_c, which is a function of h, in Eq. (1). As discussed in Sect. 2.2, the time uncertainty until 1917 can be considered negligible since the observation site and the observers' office were nearby. An estimate of the time uncertainty in manual observations from 1920 onwards can be made from the known observation times from 1964 to 1968 and from 1983 to 1986, obtaining a value of approximately ±2 h. Temperature data can be reliably obtained from the thermograph curves every half an hour; therefore, the time error can be estimated at about ±15 min. The time-related error is estimated in the same way used for the environmental error, in this case by analysing 4 and 0.5 h intervals. The mean 4 h temperature ranges are 0.50, 0.54 and 0.52 ^∘C, and the mean 0.5 h ranges are 0.17, 0.23 and 0.21 ^∘C at 0.4, 2 and 6 m, respectively. Representative errors of ±0.25 ^∘C (0.5 ^∘C band) and ±0.1 ^∘C (0.2 ^∘C band) can be adopted in these cases. Time uncertainty does not affect data set B significantly, because time is correct within a few minutes (at worst) depending on clock stability.

Table 2 Errors in estimated mean daily temperature. σ_I is the instrumental error, σ_E is the environmental error, σ_T is the time-related error, σ_o is the observational error, No. of obs. is the number of observations per day, σ is the error of the estimated daily mean. Labels (a) and (b) correspond to zero and non-zero uncertainties in time, respectively. When No. of obs. is characterised by a range, the corresponding σ range is shown. See Sect. 3.2 for detailed explanations.

Table 2 summarises the errors that can be expected in the different cases, as a function of σ_T, which is related to a known or nominal observation time, and the number of daily data available. Provided that T_c and T_24c are computed from Eqs. (1) and (2), from the observational error σ_o=0.18 ^∘C we obtain σ_c=0.05 ^∘C and σ_24c=0.01 ^∘C. Finally, the overall error of the daily temperature estimate (Eq. 3) can be written as follows:

$$\begin{array}{}\text{(6)}& \mathit{\sigma }(d,m,y,Z)={\left[{\mathit{\sigma }}_{\mathrm{24}\mathrm{c}}^{\mathrm{2}}+{\displaystyle \frac{\mathrm{1}}{N^{\mathrm{2}}}}\sum _{k=\mathrm{1}}^N\left({\mathit{\sigma }}_{\mathrm{o}}^{\mathrm{2}}+{\mathit{\sigma }}_{\mathrm{c}}^{\mathrm{2}}\right)\right]}^{\frac{\mathrm{1}}{\mathrm{2}}}.\end{array}$$

Again, $N=N(d,m,y,z)$ is the number of available observations. To the precision of 0.1 ^∘C, the daily means at 2 m depth are characterised by overall errors between 0.2 and 0.4 ^∘C for data set A and between 0.1 and 0.2 ^∘C for data set B.

3.3 The composite time series

Before merging data sets A and B, a consistency check was done by comparing the data sets in the 1986–2008 period when temperature was measured both at Porto Lido (LID) and Molo Bandiera (BAN) (Fig. 1). We computed mean temperature differences ($\mathrm{\Delta }T=T_{\mathrm{A}}-T_{\mathrm{B}}$) when simultaneous measurements were available. Two periods can be distinguished, namely 18 November 1986–30 November 1989, in which $\mathrm{\Delta }T=\mathrm{0.5}\pm \mathrm{0.5}$ ^∘C, and 1 March 1993–31 December 2008, with $\mathrm{\Delta }T=-\mathrm{0.1}\pm \mathrm{0.4}$ ^∘C. As the procedures and instruments used for direct measurements (data set A) did not significantly change before and after 1986, the temperature difference in the first period is likely related to the early digital records (data set B), perhaps due to instrumental drifts or inaccurate calibration. In the second period ΔT is negligible.

The composite time series was built using data set A as a basis and data set B to fill gaps and extend the time series after 2008. For a few gaps in 1987–1989 the temperature from data set B was increased by 0.5 ^∘C on the basis of the above-mentioned temperature difference.

4 Data availability

Figure 4 illustrates the monthly data availability as the percentage of the number of days with temperature estimate to the number of days per month. Mean annual temperatures were computed for complete years, i.e. having no missing months; a missing month has less than 80 % available days (Fig. 5). The threshold was chosen arbitrarily.

The 1899–2015 composite time series is characterised by fluctuations on decadal timescales and a linear trend of 1.1±0.3 ^∘C per century (at the 95 % confidence level). The effect of the biases induced by poor time samplings, discussed in Sect. 3.1, can be estimated by analysing the time series of the original observations; in this case the trend is 1.0±0.3 ^∘C per century. Because of the many missing years in the first part of the time series, this result should be taken with caution. The trend of the continuous 1946–2015 period is 1.3±0.5 ^∘C per century, using both the normalised and the original data.

The data used in this work are available from SEANOE as “Trieste 1899–2015 near-surface sea temperature”, which includes the composite time series, data sets A and B and the monthly and annual mean temperatures (https://doi.org/10.17882/58728; Raicich and Colucci, 2019).

Figure 4 Data availability on a monthly basis for data set A (1899–2008) and data set B (1986–2015).

Figure 5 Annual mean sea temperature at 2 m depth from 1899 to 2015 (blue dots and dashed lines) and linear trends for 1899 to 2015 (red line) and 1946 to 2015 (green line).

5 Conclusive remarks

Near-surface sea temperature data measured at different sites in the harbour of Trieste were collected, forming two data sets. The first consists of analogue data from 1899 to 2008, while the second consists of digital records from 1986 to 2015. Their merger allowed the building of a 117-year-long quasi-homogeneous time series of mean daily temperature at 2 m depth. Although the data accuracy is lower than the modern standards for ocean temperature measurements, these data sets represent a valuable tool to study sea temperature variability on timescales from the synoptic (connected to the meteorological forcing) to decadal and century timescales (related to climate changes).

The search for undiscovered data will continue, in order to possibly fill the existing gaps.