# Seawater physics and chemistry along the Med-SHIP transects in the Mediterranean Sea in 2016

Katrin Schroeder, Vedrana Kovačević, Giuseppe Civitarese, Dimitris Velaoras, Marta Álvarez, Toste Tanhua, Loïc Jullion, Laurent Coppola, Manuel Bensi, Laura Ursella, et al.

<figure><img src="image_1.png" /></figure>

To cite this version:

Katrin Schroeder, Vedrana Kovačević, Giuseppe Civitarese, Dimitris Velaoras, Marta Álvarez, et al. "Seawater physics and chemistry along the Med-SHIP transects in the Mediterranean Sea in 2016." Scientific Data, 2024, 11 (1), pp.52. 10.1038/s41597-023-02835-3. hal-04887472

HAL Id: hal-04887472

https://hal.science/hal-04887472v1

Submitted on 16 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
scientific **data**

<figure><img src="image_3.png" /></figure>

Check for updates

DATA DESCRIPTOR

OPEN

# Seawater physics and chemistry along the Med-SHIP transects in the Mediterranean Sea in 2016

Katrin Schroeder<sup>1, ∗</sup>, Vedrana Kovačević<sup>1, †</sup>, Giuseppe Civitarese<sup>2</sup>, Dimitris Velaoras<sup>1, ‡</sup>

Marta Álvarez<sup>1, ‡</sup>, Toste Tanhua<sup>5</sup>, Loïc Jullion, Laurent Coppola<sup>6,7</sup>, Manuel Bensi<sup>1, ‡</sup>

Laura Ursella<sup>2</sup>, Chiara Santinelli<sup>8</sup>, Michele Giani<sup>2</sup>, Jacopo Chiggiato<sup>1</sup>, Mohamed

Aly-Eldeen<sup>1, ‡</sup>, Georgia Assimakopoulou<sup>3</sup>, Giancarlo Bachi<sup>8</sup>, Boie Bogner<sup>5</sup>, Mireno Borghini<sup>1</sup>

Vanessa Cardin<sup>2</sup>, Marin Cornec<sup>6</sup>, Antonia Giannakourou<sup>3</sup>, Louisa Giannoudi<sup>3</sup>,

Alexandra Gogou<sup>3</sup>, Melek Golbol<sup>6</sup>, Or Hazan<sup>10</sup>, Clarissa Karthäuser<sup>5</sup>, Martina Kralj<sup>2</sup>,

Evangelia Krasakopoulou<sup>11</sup>, Frano Matić<sup>12</sup>, Hrvoje Mihanović<sup>13</sup>, Stipe Muslim<sup>12</sup>,

Vassilis P. Papadopoulos<sup>3</sup>, Constantine Parinos<sup>3</sup>, Anne Paulitschke<sup>5</sup>, Alexandra Pavlidou<sup>3</sup>,

Elli Pitta<sup>3</sup>, Maria Protopapa<sup>3</sup>, Eyal Rahav<sup>10</sup>, Ofrat Raveh<sup>10</sup>, Panagiotis Renieris<sup>3</sup>,

Nydia C. Reyes-Suarez<sup>2</sup>, Eleni Rousselaki<sup>3</sup>, Jacop Silverman<sup>10</sup>, Ekaterini Souvermezoglou<sup>3</sup>,

Lidia Urbini<sup>2</sup>, Christina Zeri<sup>3</sup> & Soultana Zervoudaki<sup>3</sup>

The Mediterranean Sea has been sampled irregularly by research vessels in the past, mostly by national expeditions in regional waters. To monitor the hydrographic, biogeochemical and circulation changes in the Mediterranean Sea, a systematic repeat oceanographic survey programme called Med-SHIP was recommended by the Mediterranean Science Commission (CIESM) in 2011, as part of the Global Ocean Ship-based Hydrographic Investigations Program (GO-SHIP). Med-SHIP consists of zonal and meridional surveys with different frequencies, where comprehensive physical and biogeochemical properties are measured with the highest international standards. The first zonal survey was done in 2011 and repeated in 2018. In addition, a network of meridional (and other key) hydrographic sections were designed: the first cycle of these sections was completed in 2016, with three cruises funded by the EU project EUROFLEETS2. This paper presents the physical and chemical data of the meridional and key transects in the Western and Eastern Mediterranean Sea collected during those cruises.

## Background & Summary

**Introduction and scientific objectives.** The ocean is variable and plays a significant role in the Earth's energy and freshwater balance<fnref n="1" />, but it is poorly sampled, particularly below 2000 m. The Mediterranean Sea is a marginal sea that plays a significant role in the regional freshwater, heat and carbon budget. It is often described as a miniature ocean, where typical oceanic processes (such as dense water formation, mesoscale dynamics, thermaline circulation) occur faster and on smaller scales than in the open ocean [e.g. ²]. The Mediterranean Sea has also experienced remarkable changes in its physical and biogeochemical properties over the past decades, such as the Eastern Mediterranean Transient event (see next section "Geographical and oceanographic settings"), the

<footnote n="1">¹Consiglio Nazionale delle Ricerche-Istituto di Scienze Marine (CNR-ISMAR), Venice and La Spezia, Italy. ²National Institute of Oceanography and Applied Geophysics - OGS, Trieste, Italy. ³Hellenic Centre for Marine Research, HCMR, Institute of Oceanography, Anavyssos, Greece. ⁴Instituto Español de Oceanografía, IEO-CSIC, A Coruña, Spain. ⁵GEOMAR Helmholtz Centre for Ocean Research Kiel, Kiel, Germany. ⁶Sorbonne Université, CNRS, Laboratoire d'Océanographie de Villefranche, Villefranche-sur-Mer, France. ⁷Sorbonne Université, CNRS, OSU STAMAR, UAR2017, 4 Place Jussieu, 75252, Paris, cedex 05, France. ⁸Consiglio Nazionale delle Ricerche-Istituto di Biofisica (CNR-IBF), Pisa, Italy. ⁹National Institute of Oceanography and Fisheries, NIOF, Egypt, EG. ¹⁰Israel Oceanographic and Limnological Research, IOLR, Haifa, Israel. ¹¹University of the Aegean, Department of Marine Sciences, Mytilene, Greece. ¹²University Department of Marine Studies, University of Split, Split, Croatia. ¹³Institute of Oceanography and Fisheries, Split, Croatia. e-mail: katrin.schroeder@ismar.cnr.it</footnote>
<figure><img src="image_4.png" /></figure>

**Fig. 1** Station map of the Med-SHIP cruises: TALPro2016 (grey circles), CRELEV2016 (black circles), ESAW2015 + ESAW2016 common stations in red (ESAW2015 only in orange and ESAW2016 only in yellow). The ports of these cruises are indicated in red letters (Pa = Palermo, Ba = Barcelona, Pi = Piraeus, Sp = Split). For completeness, MED01 GO-SHIP stations (repeated in 2011 and 2018, not presented in this paper) are shown by light grey stars. All cruises are ideally performed coast-to-coast, however permissions to work in southern and easternmost shore waters are rarely given. For schematics of main currents refer to e.g.¹¹¹.

increase of salinity and temperature in the deep layers, and the decrease of dissolved oxygen concentrations in the intermediate waters. These changes have important implications for the climate, ocean biogeochemistry, and ecosystem functioning of the region.

To observe and understand these changes, regular surveys of the Mediterranean circulation with comprehensive physical and biogeochemical measurements are needed. However, in the past, the Mediterranean has been sampled sporadically in time and space, mostly by national expeditions in regional waters. The first basin wide zonal transect has been described by³, in his early studies of Mediterranean water masses. During the 1980s and 1990s four large international programs, i.e. the Gibraltar Experiment (GIBEX⁴), the Physical Oceanography of the Eastern Mediterranean (POEM⁵), the Western Mediterranean Circulation Experiment (WMCE⁶), and the Programme de Recherche International en Mediterranée Occidentale (PRIMO⁷), set the baseline to a thorough definition of the main Mediterranean Sea characteristics. We particularly note the pioneering studies of Wolfgang Roether⁸⁻¹², who used observations of chemical transient tracers (e.g., chlorofluorocarbons, tritium), that provided insights on ventilation rates, water mass formation processes of the Mediterranean Sea and the Eastern Mediterranean Transient event. Since the end of these programs, uncoordinated ship-based research efforts, driven mainly by national interests, provided fragmented and sporadic observations.

Marginal seas such as the Mediterranean Sea and others were not originally considered in the international GO-SHIP (Global Ocean Ship-based Hydrographic Investigations Program) or in global data synthesis efforts such as GLODAP¹³,¹⁴. In 2011, an outcome of the 43rd CIESM (Mediterranean Science Commission) Workshop in Supetar (Croatia) was the recommendation for repeat oceanographic surveys within the Mediterranean Sea, in a programme called Med-SHIP¹⁵. Similarly to the international GO-SHIP programme for the global ocean, Med-SHIP aims to detect the impact of climate change in the Mediterranean Sea by comprehensive observations of physical and biogeochemical variables. These variables should be measured to the highest international standards on a regular basis along a number of zonal and meridional hydrographic sections following the GO-SHIP guidelines¹⁶. There are two primary objectives for the Med-SHIP repeat hydrography cruises: (i) to observe long-term changes in physical and biogeochemical properties and (ii) to observe changes in the thermohaline circulation, both in the Mediterranean Sea as a whole, but also on sub-basin scales. Therefore, the Med-SHIP program consists of a zonal section from the Strait of Gibraltar to the easternmost Mediterranean and a number of meridional or key sections in the Eastern and Western Mediterranean basins (Fig. 1). The zonal survey is already officially part of the GO-SHIP reference sections, named as MED01 line (shown with light grey stars in Fig. 1), with the CTD and bottle data available at CCHDO¹⁷ and biogeochemical data included in the GLODAPv2 (Global Ocean Data Analysis Project) data product¹⁸. The western and eastern sections of Med-SHIP, the subject of this paper, are in the process of being proposed to GO-SHIP as "associated lines".

This paper is the culmination of efforts in 2016 to obtain the first quasi-synoptic meridional surveys of the western and eastern Mediterranean Sea as part of MED-SHIP. Scientists from eight countries worked on three national research vessels to complete these cruises. The data reported here include temperature, salinity, pressure, inorganic nutrients, dissolved oxygen, total alkalinity, dissolved inorganic carbon, chlorophyll-a, pH, and chlorofluorocarbons (details on measured variables are given in Tables 1, 2). The results from these cruises adhere to the data policy of GO-SHIP, i.e., all data are published in open and free data repositories. The data are curated by specialized data centres, which are listed in Table 3 of section "Data Records". The Med-SHIP programme also strives to achieve the accuracy and precision requirements articulated by GO-SHIP¹⁹.
<table><thead><tr><td>Parameter name</td><td>Short name</td><td>Unit</td><td>Method</td><td>Cruise</td></tr></thead><tbody><tr><td>date/time</td><td>Date/time</td><td>UTC</td><td></td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>latitude/longitude</td><td>Latitude/ Longitude</td><td>degrees North/degrees<br>East</td><td></td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>pressure, water</td><td>Press</td><td>dbar</td><td>CTD, SBE 911+</td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>temperature, water</td><td>Temp</td><td>°C</td><td>CTD, SBE 911+</td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>salinity</td><td>Sal</td><td></td><td>CTD, SBE 911+</td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>dissolved oxygen</td><td>O<sub>2</sub></td><td>µmol kg<sup>-1</sup></td><td>CTD with a SBE 43 calibrated using<br>Winkler titration (except for ESAW)</td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>fluorescence, chlorophyll</td><td>F chl</td><td>µg L<sup>-1</sup><br>mg m<sup>-3</sup></td><td>WET Labs ECO-AFL/FL fluorometer<br>Chelsea AquaTracka III<br>WET Labs ECO-AFL/FL fluorometer</td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>turbidity beam transmission</td><td>Turbidity<br>BT</td><td>NTU%</td><td>WET Labs turbidimeter<br>Chelsea AlphaTracka Mk II<br>WET Labs transmissometer</td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>current velocity east-west</td><td>U</td><td>m s<sup>-1</sup></td><td>lowered acoustic doppler current profiler<br>(LADCP)</td><td>TalPro2016<br>CRELEV2016</td></tr><tr><td>current velocity north-south</td><td>V</td><td>m s<sup>-1</sup></td><td>lowered acoustic doppler current profiler<br>(LADCP)</td><td>TalPro2016<br>CRELEV2016</td></tr><tr><td>current velocity east-west</td><td>U</td><td>m s<sup>-1</sup></td><td>vessel-mounted acoustic doppler current profiler<br>(VM-ADCP)</td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>current velocity north-south</td><td>V</td><td>m s<sup>-1</sup></td><td>vessel-mounted acoustic doppler current profiler<br>(VM-ADCP)</td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>temperature, water</td><td>Temp</td><td>°C</td><td>SBE 21 Seacat Thermosalinograph</td><td>TalPro2016,<br>CRELEV2016</td></tr><tr><td>salinity, water</td><td>Sal</td><td></td><td>SBE 21 Seacat Thermosalinograph</td><td>TalPro2016<br>CRELEV2016</td></tr></tbody></table>

**Table 1.** Details (parameter name, short name, unit, method, cruise) on the parameters measured by the CTD multiparametric probe and underway during each cruise.

In 2016, the meridional sections in the Western Mediterranean Sea (WMED) were carried out during the TalPro2016 cruise (Tyrrhenian and Algero-Provençal cruise; stations are shown by grey dots in Fig. 1²⁰). The meridional and key sections in the Eastern Mediterranean Sea (EMED) were carried out during the CRELEV2016 cruise (Cretan Sea and Levantine cruise; stations are shown by black dots in Fig. 1²¹), and during the ESAW cruises (Evolution and Spreading of the Southern Adriatic Waters cruises; stations are shown by yellow-orange-red dots in Fig. 1²²), respectively. For the sake of clarity, it is necessary to note that the ESAW2016 cruise (called “ESAW leg 2” in the submitted datasets) was the second leg of a broader programme that included a first leg in December 2015, ESAW2015 (ESAW leg 1 in the submitted datasets). The aim was to conduct repeated measurements in the same area during two different periods (pre- and post- winter convection phases) and this is the reason why also one 2015 cruise is included in this paper.

In addition to the cruises we report on here, it is worth noting that other Med-SHIP sections have been repeated or are planned for future years: the trans-Mediterranean zonal Med-SHIP section (MED01 line) was repeated in 2011 and 2018¹⁷,²³ and is planned again for 2025 or 2026; while the western meridional sections were recently repeated in spring 2022 (cruise TalPro2022²⁴).

The Med-SHIP program also aims to foster and facilitate regional collaboration between nations along the northern shores and countries in the Middle East and North Africa. This commitment is underscored by the rich diversity within our cruise teams and the authors of this paper, reflecting a wide geographical spectrum. Central to the Med-SHIP mission is capacity building, collaborative efforts, and knowledge exchange, which serve as vital mechanisms for bridging the research and expertise gap that exists between the northern shore countries and those in the Middle East and North Africa.

**Geographical and oceanographic settings.** The Algero-Provençal basin is the main site of Western Mediterranean Deep Water (WMDW) formation driven by intense heat loss from the sea in the northern part of the basin. Since 2005 it became evident that the gradual temperature and salinity trends in the WMDW were interrupted by an abrupt shift towards higher temperature and salinity²⁵. The onset of this shift, called the Western Mediterranean Transition, has been a major dense water formation event in winter 2004/2005 in the northern part of the WMED. Gradually this anomaly spread from its formation region and filled up the basin interior. The Tyrrhenian Sea acts as a blender of various water masses originating from the WMED and EMED¹²,²⁶,²⁷ and therefore plays an important role in the preconditioning the water
<table><caption>**Table 2.** Details (parameter name, short name, unit, method, cruise) on the parameters measured on discrete samples during each cruise.</caption><thead><tr><th>Parameter name</th><th>Short name</th><th>Unit</th><th>Method</th><th>Cruise</th></tr></thead><tbody><tr><td>dissolved oxygen</td><td>O<sub>2</sub></td><td>µmol kg<sup>-1</sup></td><td>Winkler titration</td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>salinity</td><td>Sal</td><td></td><td>Salinometer</td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>nitrate</td><td>NO<sub>3</sub><sup>-</sup></td><td>µmol kg<sup>-1</sup></td><td>Autoanalyzer</td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>nitrite</td><td>NO<sub>2</sub><sup>-</sup></td><td>µmol kg<sup>-1</sup></td><td>Autoanalyzer</td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>ammonium</td><td>NH<sub>4</sub><sup>+</sup></td><td>µmol kg<sup>-1</sup></td><td>Autoanalyzer</td><td>ESAW</td></tr><tr><td>phosphate</td><td>PO<sub>4</sub><sup>3-</sup></td><td>µmol kg<sup>-1</sup></td><td>Autoanalyzer</td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>silicate</td><td>Si(OH)<sub>4</sub></td><td>µmol kg<sup>-1</sup></td><td>Autoanalyzer</td><td>TalPro2016<br>CRELEV2016</td></tr><tr><td>dissolved inorganic carbon</td><td>DIC</td><td>µmol kg<sup>-1</sup></td><td>Coulometric titration</td><td>TalPro2016<br>CRELEV2016</td></tr><tr><td>total alkalinity</td><td>TA</td><td>µmol kg<sup>-1</sup></td><td>Potentiometric titration in an open cell</td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>pH</td><td>pH</td><td>total scale<br>at 25 °C</td><td>Single beam spectrophotometer, un-purified dye<br>Single beam Spectrophotometer, purified dye</td><td>TalPro2016<br>CRELEV2016<br>ESAW</td></tr><tr><td>carbonate</td><td>CO<sub>3</sub><sup>2-</sup></td><td>µmol kg<sup>-1</sup></td><td>Single beam spectrophotometer, Pb(ClO<sub>4</sub>)<sub>2</sub></td><td>TalPro2016</td></tr><tr><td>chlorofluorocarbon (CCl<sub>3</sub>F<sub>2</sub>)</td><td>CFC-12</td><td>pmol kg<sup>-1</sup></td><td>Purge and trap Gas-Chromatograph with an Electron Capture<br>Detector system</td><td>TalPro2016<br>CRELEV2016<br>ESAW2016</td></tr><tr><td>sulfur hexafluorid</td><td>SF<sub>6</sub></td><td>fmol kg<sup>-1</sup></td><td>Purge and trap Gas-Chromatograph with an Electron Capture<br>DetectorGC-ECD system</td><td>TalPro2016<br>CRELEV2016<br>ESAW2016</td></tr><tr><td>dissolved organic carbon</td><td>DOC</td><td>µmol kg<sup>-1</sup></td><td>High Temperature Catalytic Oxidation</td><td>CRELEV2016<br>ESAW</td></tr><tr><td>total dissolved phosphorus<br>and nitrogen</td><td>TDP and<br>TDN</td><td>µmol kg<sup>-1</sup></td><td>Wet chemical oxidation and automated segmented flow<br>colorimetric analysis</td><td>CRELEV2016<br>ESAW</td></tr></tbody></table>

column leading to the formation of WMDW<sup>28,29</sup>. Fresh and warm Atlantic Water (AW) enters the WMED via Gibraltar as an eastward boundary current along the African coast referred to as the Algerian Current<sup>30</sup>. This current, subject to strong mesoscale variability leading to the generation of long-lived anticyclonic eddies<sup>27,31</sup>, separates in a weak branch entering the Tyrrhenian Sea and a stronger branch entering the EMED through the Sicily Channel. The surface circulation along the northern slope of the Algero-Provençal basin is dominated by the Northern Current<sup>27</sup> characterised by a strong seasonal variability in its mesoscale activity.

The return flow from the EMED towards the WMED is made of warm and salty Eastern Intermediate Water (EIW): this term refers generally to a mixture of intermediate water masses formed in the Levantine Sea (also known as Levantine Intermediate Water, LIW), and in the Cretan Sea (also known as Cretan Intermediate Water, CIW). The EIW crosses the Sicily Channel and is topographically constrained to veer north along the western coast of Sicily into the Tyrrhenian Sea. The Adriatic Sea used to be the main source of dense water for the whole EMED (i.e.<sup>32,33</sup>), and is also strongly influenced by climate change affecting the Mediterranean region<sup>34</sup>. The increase in temperature and salinity observed in the Adriatic Sea is due to processes that are not only local and have strong effects on the renewal and oxygenation of the deep waters throughout the Mediterranean Sea<sup>35,36</sup>. In the early 1990s, the Eastern Mediterranean Transient event took place, marking a shift in the dominant deep water source from the Adriatic Sea to the Aegean and Cretan Sea, characterized by massive overflows of dense (salty, but warm) waters<sup>8,37</sup>. After 2002 the Adriatic deep water formation restarted, and the recent Adriatic dense waters are warmer and saltier than before the Eastern Mediterranean Transient<sup>38–42</sup>. The Cretan Sea, however, remains a crucial region, which participates in Mediterranean thermohaline circulation with episodic and limited dense water formation after the Eastern Mediterranean Transient. One of the drivers for dense water formation is attributed to the decadal reversals of the upper layer Ionian Sea circulation from cyclonic to anti-cyclonic, which decreases the AW flux into the Levantine Basin, increasing its overall salinity (<sup>43</sup> and references therein).

From the biogeochemical point of view, the deep Mediterranean is characterized by a low dissolved inorganic carbon/total alkalinity ratio (DIC/TA) and a low Revelle factor, which make it more prone to anthropogenic carbon uptake than any other basin<sup>44,45</sup>. It is therefore important to quantify the rate and distribution of carbon sequestration and transport by the circulation and the biological pump in this region. Moreover, the salinity and
<table><caption>**Table 3.** For each cruise all necessary information on data access (citations) are listed. The last column contains the link to the cruise summary reports (CSR), which describe the cruise details, but do not contain any data.</caption><thead><tr><th>Cruise name (dates)</th><th>physics (CTD, Conductivity - Temperature - Depth probe)</th><th>chemistry (BTL, bottle)</th><th>LADCP (Lowered Acoustic Doppler Current Profiler)</th><th>VM-ADCP (Vessel-mounted Acoustic Doppler Current Profiler)</th><th>thermo-salinograph</th><th>NCEI OCADS link (National Centers for Environmental Information Ocean Carbon and Acidification Data Portal)</th><th>CSR link (Cruise Summary Report)</th></tr></thead><tbody><tr><td>TAIPro2016 (19.08.2016-28.08.2016)</td><td>83</td><td>84</td><td>85</td><td>86</td><td>87</td><td>88</td><td>https://CSR.seadatanet.org/report/20166062</td></tr><tr><td>ESAW2015 (10.12.2015-16.12.2015)</td><td>90</td><td>91</td><td>not collected</td><td>92</td><td>not collected</td><td>not participating</td><td>https://nodc.ogs.it/catalogs/csrdetails/2774?3</td></tr><tr><td>ESAW2016 (05.04.2016-10.04.2016)</td><td>93</td><td>94</td><td>not collected</td><td>95</td><td>not collected</td><td>not participating</td><td>https://nodc.ogs.it/catalogs/csrdetails/2191?2</td></tr><tr><td>CRELEV2016 (02.06.2016-10.06.2016)</td><td>96</td><td>89</td><td>https://CSR.seadatanet.org/report/20165899</td><td></td><td></td><td></td><td></td></tr></tbody></table>

temperature of the deep Mediterranean have shown an increasing trend over time²⁵,⁴⁶, which implies significant changes in the water mass formation and thermohaline circulation. These changes need to be documented and understood in terms of their causes and consequences. The winter convection episodes in the WMED have become weaker over time, leading to a decrease of dissolved oxygen concentrations in the intermediate waters⁴⁷. This could affect the carbon export and the mesopelagic marine ecosystems, resulting in new balances within the trophic chain (e.g.⁴⁸⁻⁵⁴).

**Cruise summaries.** *TAIPro2016 Cruise.* The TAIPro2016 cruise (Palermo, Italy - Barcelona, Spain, RV ANGELES ALVARINO, 18.08.2016 - 29.08.2016) occupied three hydrographic transects, one across the relatively flat and deep Algero-Provençal basin, one across the deep and rugged Tyrrhenian Sea, and an additional one across the Balearic Sea (not foreseen in the Med-SHIP design). A total of 43 stations were realised (grey circles in Fig. 1), where physical (temperature, salinity, pressure, velocity) and biogeochemical parameters (oxygen, nutrients, carbonate system) and transient tracers (CFCs and SF₆) were measured.

The main focus of the cruise was the north-south transects across the WMED, between its northern and southern shores. Along the different sections, CTD stations including sampling of chemical parameters were conducted approximately every 10–30 nm. In addition, underway thermosalinograph and vessel-mounted ADCP measurements were performed between CTD stations.

The specific objectives of the TAIPro2016 cruise were (i) to observe the state of the Western Mediterranean Transition, (ii) to quantify the ventilation of different parts of the WMED, (iii) to assess the rate of warming and salinification, mainly of the intermediate and deep layer, (iv) to monitor the biogeochemical status of the basin, and (v) to update the quantification of the uptake of anthropogenic carbon in the deep basin.

*CRELEV2016 Cruise.* The CRELEV2016 cruise (Piraeus - Piraeus, Greece, RV AEGAEO, 02.06.2016 – 10.06.2016) was carried out in the Cretan Sea and in the Cretan Passage, south of Crete. The cruise was restricted within the limits of the Greek exclusive economic zone. In total 29 stations were sampled where CTD multiparameter profiles and samples for dissolved oxygen and nutrients were collected (black circles in Fig. 1), according to the GO-SHIP recommendations. On a subset of stations other samples were collected for biogeochemical (carbonate system, dissolved and particulate carbon, nitrogen and phosphorus) and biological (chlorophyll, plankton abundance and biomass, pigment analysis, others) parameters.

The main objective of the CRELEV2016 cruise was to examine a comprehensive set of parameters in the EMED in order to quantify variability and trends of physical and biogeochemical properties. The specific objectives of the cruise were (i) to observe the present post-convective state of the Cretan Sea, in relation to the state of the upper layer Ionian circulation, (ii) to observe the present state and to quantify changes of the hydrographic and biogeochemical properties of the Cretan Sea, the Cretan Strait and the Levantine Basin, (iii) to quantify the uptake of anthropogenic carbon in the EMED, (iv) to quantify changes in water formation rates and circulation in the EMED.

*ESAW Cruises.* The ESAW cruises (Split - Split, Croatia, RV BIOS DVA, 10.12.2015 – 15.12.2015 and 05.04.2016 – 10.04.2016) were conducted to obtain the multidisciplinary characterization of the Adriatic waters before and after the wintertime conditions, i.e., December 2015 and April 2016, hereinafter referred to as ESAW2015 and ESAW2016. The cruises were carried out in the middle and southern Adriatic Sea.

A detailed dataset was collected to compare physical and biogeochemical characteristics during pre- and post-wintertime conditions. The Adriatic Sea is characterised by a high variability due to the influence of local climate and meteorological conditions, as well as due to the remote effects of the lateral water exchange within the Adriatic and between the Adriatic and Ionian seas. Data were collected along two shore-to-shore key transects (Gargano-Split and Bari-Dubrovnik) and in the mid Adriatic depressions. A total of 55 hydrographic stations were completed (24 during ESAW2015 and 31 during ESAW2016, see red + orange + yellow circles
in Fig. 1). The scientific interest was focused on (i) the regeneration processes occurring in the deep waters in the middle and southern Adriatic pits, and (ii) on the distribution of the EIW in the Adriatic Sea, which has an important role in preconditioning the dense water formation, and in imprinting biogeochemical and biological processes.

**Data provenance.**

(a) CTD multiparametric probe and underway acquisitions.

(b) Discrete samples.

**Methods**

**Conductivity-temperature-depth (CTD) measurements.** *CTD setup and configuration.* During all cruises, the CTD-rosette systems were equipped with a sonar altimeter to reach the deepest point at each station maintaining a safe distance from the bottom of about 5–10 m, depending on the sea conditions. The CTDs were left soaking at the surface for a few minutes before the beginning of the downcast in order for the water pump to activate and until the conductivity cells were stabilized.

Details on Sea-Bird SBE911 + CTD sensor configurations used during the TAIPro2016 cruise are shown in the cruise report²⁰. The 24-bottles rosette system SBE32 was equipped with 12 litres Niskin bottles. The rosette was coupled with an SBE43 oxygen sensor, as well as with WETLabs Fluorometer and Turbidimeter.

During the CRELEV2016 cruise a Sea-Bird SBE911 + CTD system attached to a SBE32 carouse water sampler was used, equipped with 12 Niskin bottles (12 litres capacity). An SBE43 oxygen sensor was also used, as well as a Chelsea AquaTracka III fluorometer and a Chelsea AlphaTracka Mk II transmissometer.

During the ESAW cruises the CTD used was a Sea-Bird SBE911 +, coupled with an SBE43 oxygen sensor, a WETLabs FLNTURTD sensor for fluorescence and turbidity, and a WET Labs C-Star transmissometer. The CTD was attached to a rosette sampler holding 12 Niskin bottles (8 litres capacity) for water sampling. A few additional CTD casts were conducted using a separate self-recording SBE25 probe (pressure, temperature and conductivity sensors), sampling at 8 Hz.

**CTD operations and data processing.** For all cruises, the vertical profiles of all parameters were obtained by sampling the signals at 24 Hz, with the CTD/rosette going down at a speed of about 1 m/s. The data were transmitted on-line to the ships' computer and viewed in real time using Sea-Bird SeaSave software (https://www.seabird.com/asset-get.download.jsa?id=54627862734). Data were pre-processed on board using SBE Data Processing software (https://www.seabird.com/asset-get.download.jsa?id=54627862733) in the following steps according to the CTD manufacturer's recommendations: data conversion, low-pass filtering of pressure readings, temporal alignment of CTD's oxygen sensor, cell thermal mass correction, computing of derived variables, bin averaging at 1 dbar bins and splitting of downcast and upcast. Derived parameters, such as the potential temperature and potential density anomaly referred to 0 dbar, were calculated from original *in situ* data using the TEOS-10 Gibbs function approach⁵⁵. Details on CTD calibration and corrections can be found in Section "Technical Validation".

**Lowered acoustic doppler current profiler.** The model RDI Workhorse 300 kHz Lowered Acoustic Doppler Current Profiler (LADCP) was used along with CTD casts to measure current velocity and magnitude. During the TAIPro2016 cruise, two LADCP (a master facing down and a slave facing up, see Fig. 2) were used, while during the CRELEV2016 cruise a single downward-looking LADCP was used. For data post-processing the LDEO LADCP (version 10.16 and version IX, respectively for TAIPro2016 and CRELEV2016) software⁵⁶ was used, which converts the raw LADCP data, combines them with processed CTD data and navigational data, to give the u- and v- velocities profiles from surface to bottom. No LADCP data were collected during the ESAW cruises.

**Discrete water sample analyses.** *Salinity (measurement and calibration).* Water samples were taken at most stations to measure conductivity/salinity and to identify any bias or drift in the conductivity sensors. Standard sample procedure was used⁵⁷. Conductivity was measured with a Guildline Autosal 8400B (TAIPro2022 and ESAW cruises) and 8410 A (CRELEV2016). Once collected, the samples were brought into the room hosting the salinometer to give them time to adjust to the room temperature. The temperature of the water in which the conductivity cell is bathed must be as close as possible to the room temperature and needs to be around 21 °C for optimal results.

The instrument is standardized at the beginning and the end of a set of measurements with Standard IAPSO Seawater. For each sample, three readings of the conductivity ratio were taken (or as many samples as needed to obtain sufficient stability). Salinometer and CTD salinities were compared, new calibration coefficients were computed, and applied to the CTD sensor-derived raw values to obtain corrected salinity values.

*Dissolved oxygen (measurement and calibration).* Dissolved oxygen (DO) concentrations were measured at all CTD casts with a SBE43 sensor (Sea-Bird Clark electrode polarographic sensor), which has an accuracy of 2% DO saturation (around 5 µmol kg⁻¹ at 13 °C with salinity 38). To adjust the sensor response from the drift and offset, seawater samples were collected at almost all stations from surface to deep waters for discrete Winkler analysis performed onboard.
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Slave (upward looking, left picture) and master (downward looking, right picture) LADCP set up on the rosette during the TALPro2016 cruise.</figcaption></figure>

Winkler-titrated water samples were used as the DO reference concentrations⁵⁷. Sampling was done with the recommended precautions to prevent any biological activity and gas exchanges with the atmosphere⁵⁸. Temperature of the water in the sampling bottle was measured during sampling with a digital thermometer. Samples were then immediately fixed after collection with the reagents (#1: manganous chloride and #2: sodium hydroxide – sodium iodine) added with semi-automated dispensers according to⁵⁹,⁶⁰. The bottles were shaken for about 1 minute to bring each molecule of dissolved oxygen into contact with manganese (II) hydroxide. After fixation the precipitate was allowed to settle below the half of the bottle (min. 4 hours). The bottles were kept in the dark. Before titration, the precipitated hydroxides were dissolved with sulfuric acid. During TALPro2016 titration was performed directly in the sample bottle, while during CRELEV2016 the solution was transferred carefully into a titration beaker. This was done because of technical problems with the titration equipment, which prevented us from performing the titration in the sample bottle. The titration was carried out with a standardized thiosulphate solution using a Dosimat Metrohm (for TALPro2016 and CRELEV2016). The thiosulfate solution was calibrated by titrating it against a potassium iodate certified standard solution of 0.0100 N (CSK standard solution WAKO).

During ESAW cruises, water samples for DO determination were collected in 50 mL calibrated BOD (Biological Oxygen Demand) bottles and concentrations were measured by the Winkler iodometric titration using a Mettler-DL21 titrator for the determination of the equivalence point⁶¹. Temperature of the water in the sampling bottle was measured during sampling with a GMH 3475 digital thermometer. The precision of the measurements evaluated on five determinations on subsamples collected from the same Niskin bottle was 0.1%. After the cruise, these data have been used to check the sensor's response. After a thorough inspection and comparison, the SBE43 sensor data were retained good enough and were kept without change.

**Dissolved inorganic and organic nutrients.** During TALPro2016, samples for dissolved inorganic nutrients were collected at all stations from surface to deep waters (from 9 to 24 depth levels each cast). These unfiltered samples were collected through 60 mL polyethylene flasks, poisoned with 100 µl HgCl₂ (6 g/L for 100 µl of solution for 20 ml of samples, i.e. 20 µg L⁻¹) and stored in a dark room at 5 °C until analysis⁶². In previous works, comparisons have been made between freezing and poisoning, with no significant differences found (SOMLIT Intercomparisons, Service d'Observation en Milieu Littoral, https://www.somlit.fr/en/intercomparisons/). Samples can be stored for up to a year in the dark at room temperature⁶². All nutrient samples were analysed by a standard automated colorimetric system, using a Seal Analytical continuous flow AutoAnalyser III (AA3) at IMEV (Institut de la Mer de Villefranche, Sorbonne University & CNRS). Nitrate (NO₃⁻) ions were analysed according to the protocol⁶³ with a detection limit of 0.01 µmol L⁻¹. Phosphate (PO₄³⁻) ions were analysed according to the protocol⁶⁴ with a detection limit of 0.02 µmol L⁻¹. Silicate (Si(OH)₄) ions were analysed according to the protocol⁵⁸ with a detection limit of 0.02 µmol L⁻¹.

During CRELEV2016 nutrients were sampled at 29 stations at all the depths of closure of the Niskin bottles. Overall, 561 unfiltered samples were collected. Two subsamples were collected for each sampling depth. For the onboard analyses, seawater was collected in acid washed polycarbonate 10 mL vials and preserved at 4 °C until analysis, another subsample was collected in acid washed 20 mL HDPE vial and was immediately stored at -20 °C. The determination for nitrate, nitrite, phosphate and silicate were performed onboard by a colorimetric four-channel Continuous Flow Analyzer QuAAtro Seal Analytical following the methods reported in⁶⁵. In addition, during CRELEV2016 cruise, the MAGIC25 procedure⁶⁶ was followed for the determination of low phosphate concentrations. In the ultra-oligotrophic euphotic zone of the EMED the determination of phosphate with the standard phosphomolybdic blue method is difficult, because phosphate concentrations are typically at or below the detection limit of the standard method⁵⁸. MAGIC25 procedure was used instead of the classic phosphomolybdic blue method when the phosphate concentration was below 200 nM, in the layer 0–200 m, at all 29 stations. According to the MAGIC25 method, phosphates were pre-concentrated by 25 times, thus the procedure allowed detection of nanomolar concentration in seawater (detection limit ~1 nM). In total MAGIC25 measurements were performed at 263 samples in duplicates, for a total of 526 phosphate analyses
on board. For the determination of phosphates with MAGIC25 procedure, 250 mL of seawater was sampled in duplicates from the Niskin bottles at high density polycarbonate bottles of 250 mL. Turbidity blanks in duplicates were determined at several stations and depths. The bottles were then centrifuged with NaOH. The final volume of the MAGIC25 concentrate was 10 mL, which corresponded to a 25 times preconcentration. The colorimetric procedure was finally performed on board using a PERKIN-ELMER 25 Lambda spectrophotometer. Three calibration curves were performed on board, in duplicates, in seawater taken from the sampling area at 50 m depth and were used for the calculation of the phosphate concentration in the seawater samples. Moreover, during CRELEV2016 a total of 276 samples for total dissolved nitrogen and phosphorus were collected at 16 stations at standard depths. Samples for the analyses of Total Dissolved Nitrogen (TDN) and Phosphorus (TDP) were filtered on 0.2 µm filters in a Nalgene polysulfone 0.5 litres filtration system under mild vacuum and collected in 125 mL HDPE bottles from the Niskin bottle. The filtrate was collected in acid washed polyethylene vials rinsed with filtered seawater and immediately frozen (–20 °C) until laboratory analysis. The analysis was performed by wet persulfate oxidation, according to⁶⁵. Dissolved organic nitrogen (DON) and phosphorus (DOP) have been computed as follows: DON = TDN – (N-NO₂ + N-NO₃) and DOP = TDP – P-PO₄.

During ESAW cruises samples for nitrate, nitrite, ammonium, phosphate and silicate analysis were collected with prefiltration on GF/F (glass fiber) Whatman filters in polyethylene 100 mL capacity vials and were immediately stored at –20 °C. The analyses have been performed in the onshore National Institute of Oceanography and Applied Geophysics (OGS, Italy) laboratory by a colorimetric four-channel continuous Flow QuAAstro Seal Analyzer following the methods reported in⁶⁵. Dissolved organic nitrogen (DON) and phosphorus (DOP) have been computed as follows: DON = TDN – (NO₂⁻ + NO₃⁻ + NH₄⁺) and DOP = TDP – PO₄³⁻.

**CO₂ system variables.** During the TA!Pro2016 cruise four CO₂ variables were measured: pH, Total Alkalinity (TA), Dissolved Inorganic Carbon (DIC) and carbonate ion concentration (CO₃²⁻). During CRELEV2015 three CO₂ variables were measured (pH, TA, DIC) and during ESAW cruises two variables (pH and TA) were measured. A different subsection is devoted for each cruise to take into account some differences in the methodologies, although there are overlaps in the methods that were used.

**TA!Pro2016 cruise.** Due to on board time limitations, a complete sampling scheme for all stations and depths was unachievable: no CO₂ data was sampled in the Balearic Sea, while in the other two sections (Fig. 1) DIC, pH, TA, and CO₃²⁻ covered respectively 55%, 90%, 31% and 34% of the total discrete depths sampled. In the following sections a brief description of the methods and quality control schemes for the different CO₂ variables is given, while more information can be found in the TA!Pro2016 cruise report²⁰.

**DIC determination.** DIC samples were analysed on board with a VINDTA (Versatile INstrument for the Determination of Total inorganic carbon and titration Alkalinity) 3D system (www.marianda.com) coupled to a UIC 5011 coulometer, briefly, seawater CO₂ was extracted adding phosphoric acid within a stripping chamber, the generated CO₂ was carried by a pure N₂ gas stream into a coulometer cell where the coulometric titration was performed⁶⁷. Samples for DIC were collected after transient tracers and dissolved oxygen, preferably in 500 mL borosilicate bottles, alternatively in 200 mL topaz bottles. Bottles were rinsed and filled smoothly from the bottom till overflowing, leaving a headspace of 1% the bottle volume. Samples were stored at room temperature in the dark until analysis, maximum 3 days after sampling.

**pH determination.** Spectrophotometric pH in seawater was measured following⁶⁸ using a manual approach. Samples were collected in cylindrical optical glass 10-cm pathlength cells, which were filled to overflowing and immediately stopped. All the absorbance measurements were obtained in the thermostated (25 ± 0.2 °C) cell compartment of a BECKMAN DU800 single beam spectrophotometer. After blanking with the sampled seawater without dye, 50 µL of the unpurified m-cresol purple (2 mM, Sigma Aldrich) solution was added to each sample. The absorbance was measured at four fixed wavelengths (434, 487.6, 578 and 730 nm). pH in the total hydrogen ion concentration scale, is calculated using the formula by⁶⁸ and reported at 25 °C, hereinafter pH25T. The injection of the indicator slightly perturbs the sample pH, the double addition correction was applied. The magnitude of this correction over our pH range is small, ranging from 0.004 to 0.005 pH units, which are added to the uncorrected pH values.

**TA determination.** TA was analysed following the double end point potentiometric technique by⁶⁹, further improved by⁷⁰. This technique is faster than the whole curve titration, but comparable⁷¹. Seawater samples for TA were collected after pH samples in 600 mL borosilicate bottles, filled to overflowing and immediately stopped. Samples were stored in the dark at laboratory temperature until analysis, usually no later than 2 days. The TA potentiometric titrations were done with an automatic potentiometric titrator combined with a glass electrode and a Pt-1000 temperature probe. The system was coupled with a 5 mL exchangeable unit filled with 0.1 N hydrochloric acid to a final pH of 4.40⁶⁹. The electrode was standardised with a 4.41 pH phthalate buffer made in CO₂ free seawater⁷⁰. Salinity variations after the titration were considered in the final TA calculation. Each sample was analysed twice and the reported values are means of two values.

**CO₃²⁻ determination.** The carbonate ion concentration was determined manually following the spectrophotometric method as in⁷² with the modifications proposed by⁷³. CO₃²⁻ samples were collected after TA in cylindrical optical quartz 10-cm pathlength cuvettes, filled to overflowing and immediately stopped. After sampling the cells are immediately stabilised at 25 °C using the same temperature control as for pH. All the absorbance measurements were obtained in the thermostated (25 ± 0.2 °C) cell compartment of a BECKMAN DU800 single beam spectrophotometer. After blanking with the sampled seawater without dye, 20 µL of the titrant Pb(ClO₄)₂
solution (0.022 M, Fisher Scientific, 99.99% purity dissolved in distilled water) was added to each sample. The absorbance was measured at three wavelengths (234, 250 and 350 nm). The final $CO_3^{2-}$ concentration at 25 °C was calculated with the$^{73}$ formula which includes a titrant induced perturbation term.

**CRELEV2016 cruise.** DIC determination. Seawater samples for the determination DIC were drawn into 500–1000 mL borosilicate glass bottles sealed with glass stoppers according to the collection methods of$^{74}$. The samples were poisoned with saturated mercuric chloride solution to prevent further biological activity, stored at room temperature away from light and were analysed onshore within 6 months. Totally, 184 samples from 13 stations were collected for the determination of DIC. Samples were analysed by a VINDTA 3 C system based on a coulometric procedure using 10% phosphoric acid and a continuous nitrogen carrier gas supply with adjusted pressure to 1.5 bar, in a closed cell. The sea water sample (20 mL) is acidified with phosphoric acid (10% in 0.7 M NaCl), which converts all carbonate species into $CO_2$ gas. The generated $CO_2$ is carried into the coulometric cell using an inert gas (N$_2$), and titrated coulometrically using internally prepared anode and cathode solutions. The sampling, titration and DIC calculation of the sea water samples is controlled by LabVIEW TM software (https://www.ni.com/it-it/shop/labview.html). All samples were measured at 25 °C with a temperature-controlled water-bath. As a guidance for the laboratory work and preparation of internal standards and other material, the operating manual for the VINDTA 3 C by$^{75}$ was used.

**pH determination.** Samples were collected directly in the spectrophotometric cylindrical cells with 10 cm path length by rinsing twice and filling by overflowing twice the cell volume and stoppering with Teflon caps. All pH spectrophotometric measurements were performed onboard within 3 hours from subsampling. The samples of station 3 and station 5 were maintained and measured at 20 °C, while samples of stations 7, 9, 11, 13 and 15 were maintained and measured at 18 °C, closer to *in situ* temperature. But, because of condensation on the cell walls and of thermoregulation problems, from station 18 the samples were maintained and measured at 25 °C, to be closer to ambient temperature. Overall, 247 pH measurements were performed onboard. Samples were collected at standard depths at 16 stations. Samples were measured using a single beam spectrophotometer (Cary 50 Scan UV-Visible) with a thermostatic cell holder and a purified indicator dye m-cresol purple$^{76}$ following the Standard Operating Procedure (SOP6b$^{74}$). The temperature in the cell holder was maintained constant by a cryothermostat Lauda ECO Silver and the samples were pre-equilibrated at established temperature in a thermostat Pol-ECO thermostat.

**TA determination.** Samples were collected into 250 mL narrow-necked borosilicate bottles filling twice the volume from the bottom. Each bottle was poisoned with 100 µL of saturated mercuric chloride ($HgCl_2$), to halt biological activity, and immediately stored at 4 °C in the dark until analysis. Total alkalinity was determined within 6 months in an onshore OGS laboratory by potentiometric titration in an open cell (SOP 3b$^{74}$) on 35 samples collected at stations 23, 25, 27 and 28. TA analyses were also performed on the samples collected at the other stations by VINDTA3C instrument at HCMR laboratory. This instrument, allows the determination of TA using an open cell potentiometric acid titration method$^{74}$, however, the results from these analyses were discarded due to technical issues encountered at the electrode in the alkalinity cell, therefore only TA values measured on a 35 samples at OGS laboratories by open cell titration were kept.

**ESAW cruises.** pH determination. The pH was measured using a single-beam spectrophotometer (Cary 50 Scan UV-Visible) with a thermostated cell holder and a purified indicator dye m-cresol purple following the SOP3b$^{74}$. The temperature in the cell holder was maintained at 14.0 °C by a cryothermostat Lauda ECO Silver. To avoid $CO_2$ gain or loss, unfiltered seawater samples were collected directly into quartz cuvettes with 10 cm path-length, leaving no head space and pre-equilibrated at 14.0 °C in a Pol-ECO thermostat. pH was always measured within 1–3 hours from subsampling.

**TA determination.** For the total alkalinity, samples were pre-filtered on glass-fibre filters (Whatman GF/F) into 250 mL narrow-necked PVC bottles. Filtration was performed to remove phytoplankton cells and particles of $CaCO_3$, derived from calcifying organisms, which respectively can interact with HCl titrant solution and dissolve during measurements due to acid additions, inducing a relevant error in the exact estimation of the TA concentration$^{77}$. Each bottle was poisoned with 100 µL of saturated mercuric chloride ($HgCl_2$) to halt biological activity and stored at 4 °C in the dark until analysis. TA was determined within 6 months in an onshore OGS laboratory by potentiometric titration in an open cell$^{74}$.

**Dichlorodifluoromethane and sulfur hexafluoride.** Measurements of the transient tracers dichlorodifluoromethane (CFC-12 or freon-12) and sulfur hexafluoride (SF6) are used to characterize ventilation in the Mediterranean, and particularly temporal changes in ventilation (e.g.$^{44,78,79}$). During the TA1Pro2016 and CRELEV2016 cruises, the gas chromatograph/purge and trap (GC/PT) systems were used for the measurements of the transient tracers CFC-12 and SF6. The systems are modified versions of the set-up normally used for the analysis of CFCs$^{80,81}$. Samples were collected in 250 mL ground glass syringes. An aliquot of about 200 mL of the samples was injected into the analytical systems. The analytes were stripped out of the water phase by a flow (120 mL/min) of ultra-clean N2 during 5–6 minutes to the trap. The trap consists of 100 cm of 1/16” tubing packed with 70 cm Heysep D and is kept at –60 °C during the trapping phase. The trap was desorbed at 120 °C and the analytes passed on to the GC. The GC was setup with a 1/8” main column packed with 180 cm Carbograph 1AC (60–80 mesh) and a 20 cm Molsieve 5 A postcolumn, kept isothermal at 50 °C. The pre-column was packed with 10 cm Porasil C and 20 cm Molsieve 5 A in a 1/8” stainless steel column but had to be shortened
<figure><img src="image_6.png" /><figcaption>**Fig. 3** (left) Sea Surface Temperature from Copernicus Marine Environment Service (CMEMS) reanalysis (August 2016) and temperature measured by the thermosalinograph during TALPro2016; (right) salinity measured by the thermosalinograph during CRELEV2016.</figcaption></figure>

to 15 cm Porasil C and 20 cm Molsieve 5 A. Detection was performed on an Electron Capture Detector (ECD). This set-up allowed efficient analysis of SF6 and CFC-12.

During the ESAW cruises water samples for the determination of CFC-12 and SF₆ concentrations were collected from the Niskin bottles using 300 mL glass ampules. The ampules were directly attached to the Niskin bottles with a stainless-steel tubing system to prevent contact with the atmosphere during the sampling process (e.g.⁸²). The ampules were flushed with 3 times the volume during sampling. The samples were cooled in a water bath at ~10 °C to prevent outgassing of the trace gases before they were flame sealed and stored in aluminium boxes during the cruise; flame sealing was always performed immediately after sampling. The measurements of the samples were conducted onshore in the tracer lab at GEOMAR Helmholtz Centre for Ocean Research Kiel within one year since sampling. The two transient tracers were simultaneously analysed using a purge and trap GC-ECD system as described in⁸³.

**Additional variables.** A number of additional variables, not strictly included in the Med-SHIP/GO-SHIP sampling strategy, which are not further discussed here, have been collected: nitrogen and oxygen isotopes during TALPro2016 and particulate organic carbon (POC), particular nitrogen (PN), carbohydrates (total, mono- and poly-saccharides), and chlorophyll-a biomass during CRELEV2016, chlorophyll-a biomass, phytoplankton population structure (microplankton and nanoplankton), microbiology sampling, zooplankton sampling and ben-thos sampling during ESAW. Only the additional variables of the CRELEV2016 cruise are included in the dataset submitted to the repository listed in Table 3.

**Underway measurements.** *Thermosalinograph.* During TALPro2016 and CRELEV2016 a thermosalinograph collected underway temperature and salinity data (plus fluorescence for TALPro2016) along the cruise routes (Fig. 3). No thermosalinograph data were collected during the ESAW cruises. An SBE21 thermosalinograph was used during CRELEV2016, while during TALPro2016 the thermosalinograph system consisted of a SBE21 together with a SBE38 thermometer. Both systems worked independent from each other throughout the cruise. While temperature is taken at the water inlet in about 5 m depth (SBE38), salinity is estimated within the interior thermosalinograph from conductivity and interior temperature (SBE21). A Turner fluorometer (10-Au-005) was connected to the system.

*Vessel-mounted acoustic doppler current profiler.* The hydrographic data set comprises current measurements by hull-mounted ADCPs (Fig. 4). During TALPro2016 a RDI Ocean Surveyor (150 KHz) has worked along the ship's route. The depth range of the current profiler is about 300 m (number of cells was set to 60, bin size to 8 m). During CRELEV2016 a RDI Ocean Surveyor (75 KHz) was used to monitor ocean currents up to depths of 700 m (number of bins was set to 50, bin size to 15 m). Both systems ran in narrowband mode. During the ESAW cruises, a RDI broadband WorkHorse Mariner (300 KHz) registered horizontal velocities with a range of about 90 m (number of cells was set to 50, bin size to 4 m).

All systems were connected to the Global Positioning System (GPS) and controlled by computers using the conventional RDI VMDAS (Vessel-mounted Data Acquisition System, TRDI) software package, installed under a MS Windows system with a pinging set to be as fast as possible. No interferences with other used acoustical instruments were observed.

## Data Records

The datasets described in this paper have been submitted to repositories, ensuring accessibility and long-term preservation, as detailed in Table 3. The data are available at the following repositories: PANGAEA (83–87), NCEI/NOAA (88,89 and NODC-OGS) (90–96).

The details on file formats and types are described in the following sections. It is essential to underline that the primary objective of this paper is to offer a comprehensive depiction of cruises conducted under the Med-SHIP initiative. These cruises took place on diverse vessels, led by independent teams, each operating autonomously.
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>Fig. 4 Currents measured during TAIPRO2016 at 52 m depth (upper panel), during CREVLEV2016 at 53 m depth (middle panel), during ESAW2015 (bottom left panel) and ESAW2016 (bottom right panel) at 50.5 m depth.</figcaption></figure>

Consequently, the data records from these cruises exhibit a degree of heterogeneity and inconsistencies, evident in variations such as dataset naming conventions and repositories used. Through this paper, we underscore the imperative for improved coordination in forthcoming basin-wide efforts within the Med-SHIP framework. This coordination should aim for a high degree of consistency in the provision and description of collected data.

**Physics data of TAIPRO2016 cruise.** This is a tab-delimited dataset containing the data description<sup>83</sup>, the events (station listing and details on coordinates and time), the parameter list and corresponding columns (depth, pressure, temperature from 2 sensors, conductivity from 2 sensors, dissolved oxygen, turbidity, fluorescence,
salinity from 2 sensors, density, sigma-theta, potential temperature, absolute salinity, conservative temperature) and the licence.

**Chemistry data of TAIPro2016 cruise.** This is a tab-delimited dataset containing the data description⁸⁴, the events (station listing and details on coordinates and time), the parameter list and corresponding columns (pressure, salinity, temperature, sulfur hexafluoride, conductivity, freon-12, dissolved oxygen from sample, dissolved oxygen from CTD, nitrite, nitrate, phosphate, silicate, dissolved inorganic carbon, total alkalinity, pH, carbonate ion, turbidity, fluorescence, salinity, density, sigma-theta, potential temperature, absolute salinity, conservative temperature), parameter flags and the licence.

**Lowered acoustic doppler current profiler data of TAIPro2016 cruise.** This is a tab-delimited dataset containing the data description⁸⁵, the events (station listing and details on coordinates and time), the parameter list and corresponding columns (depth, east-west current velocity, north-south current velocity, vertical current velocity, elapsed time) and the licence.

**Vessel-mounted acoustic doppler current profiler data of TAIPro2016 cruise.** This is a netCDF file with an extensive header describing metadata⁸⁶, parameters names and units, as well as the processing steps, containing the following parameters: depth, zonal velocity component, meridional velocity component, received signal strength, percent good pings, editing flags, ship heading, ADCP transducer temperature, number of pings averaged per ensemble, ship zonal velocity component, ship meridional velocity component.

**Thermosalinograph data of TAIPro2016 cruise.** This is a tab-delimited dataset containing the data description⁸⁷, the events (with details on coordinates and time), the parameter list and corresponding columns (depth, elapsed time, temperature from 2 sensors, salinity from 2 sensors, fluorescence, conductivity from 2 sensors) and the licence.

**Physics data of ESAW cruises.** These are comma-separated values (csv) datasets in Ocean Data View (ODV) data format (https://odv.awi.de)⁹⁰,⁹³, containing the cruise details, the events (station details on coordinates and time), the parameter columns (bottom depth, pressure, conductivity from 2 sensors, temperature from 2 sensors, salinity from 2 sensors, dissolved oxygen, fluorescence, turbidity from 2 sensors, distance from the bottom, smoothed dissolved oxygen, smoothed fluorescence, smoothed turbidity).

**Chemistry data of ESAW cruises.** These are comma-separated values (csv) datasets containing the cruise details⁹¹,⁹⁴, the events (station details on coordinates, time and bottle position), the parameter columns (bottom depth, pressure, depth, temperature, salinity from CTD, salinity from sample, dissolved oxygen from CTD, dissolved oxygen from sample, nitrite, nitrate, phosphate, ammonium, silicate, total dissolved nitrogen, total dissolved phosphorus, total alkalinity, pH on total scale, in situ pH, dissolved organic carbon, sulfur hexafluoride, freon-12).

**Vessel-mounted acoustic doppler current profiler data of ESAW cruises.** These are comma-separated values (csv) datasets in Ocean Data View (ODV) data format (https://odv.awi.de)⁹²,⁹⁵, containing a header with the cruise and station details, and the following parameters: depth, eastward horizontal velocity, northward horizontal component.

**Physics data of CRELEV2016 cruise.** The single compressed file in⁹⁶ contains all cruise files⁹⁶. The physical data are contained in the file named CTD_dataset CRELEV-2016.txt. This is a text ASCII file dataset in Ocean Data View (ODV) data format (https://odv.awi.de) containing a header with the cruise details, information on each station, and the following parameters: bottom depth, pressure, temperature, dissolved oxygen, fluorescence, turbidity, potential temperature, salinity, potential density anomaly, corrected dissolved oxygen.

**Chemistry data of CRELEV2016 cruise.** The single compressed file in⁹⁶ also contains the chemical data, in the file named CRELEV-2016_bottle.txt⁹⁶. This is a text ASCII file dataset in Ocean Data View (ODV) data format (https://odv.awi.de) containing a header with the cruise details, information on each station, and the following parameters: bottom depth, pressure, temperature, salinity, turbidity, fluorescence, chlorophyll, dissolved oxygen from CTD, dissolved oxygen from samples, nitrite, nitrite + nitrate, nitrate, phosphate, phosphate with MAGIC procedure, silicate, total dissolved nitrogen, total dissolved phosphate, dissolved inorganic carbon, total alkalinity, in situ pH, freon-12, sulfur hexafluoride, dissolved organic carbon, particulate nitrogen, particulate organic carbon, monosaccharides, total saccharides, polysaccharides.

**Lowered acoustic doppler current profiler data of CRELEV2016 cruise.** The single compressed file in⁹⁶ also contains the LADCP data, in a folder named ADCP-2016⁹⁶, and this folder contains one file for each LADCP/CTD station (station_name.txt). These are text ASCII files containing a header with the station details and the following columns: depth, east-west current velocity, north-south current velocity, error velocity.

**Vessel-mounted Acoustic Doppler Current Profiler data of CRELEV2016 cruise.** The single compressed file in⁹⁶ also contains the VM-ADCP data⁹⁶, in a folder named ADCP-route, and this folder contains the whole route data acquisition divided in 31 files (route#.txt) of different lengths. These files are text ASCII files that are organized as a matrix where the first 3 rows contain details on time and coordinates and below there are the following rows (corresponding to each time step and lon/lat pair): depth, east-west current velocity, north-south current velocity.
<figure><img src="image_11.png" /><figcaption>**Fig. 5** Scatter plot of potential temperature and salinity in the WMED (a), in the Adriatic Sea (b) and in the EMED (c). Grey contour lines indicate potential density anomaly isolines (kg m⁻³).</figcaption></figure>

**Thermosalinograph data of CRELEV2016 cruise.** The single compressed file in⁹⁶ also contains the thermosalinograph data⁹⁶, in the file named CRELEV-2016_Thermosalinograph.txt. This is a text ASCII file dataset containing details on time and coordinates as well as the following parameters: temperature, salinity

## Technical Validation

To ensure the quality in terms of precision and accuracy of the collected data, we performed several checks and validations on the data before processing and analysing them. In particular, we checked and quality controlled the chemistry data from Niskin bottles using various methods. In general accuracy is ensured by following best practices recommendations⁵⁷. The different methods are described for each measurement in the following paragraphs.

**Temperature and salinity data.** For all cruises, we followed the best practices procedures⁵⁷ and recommendations from the CTD manufacturer. Conductivity sensors have been recently calibrated: for TalPro2016 cruise all sensors mounted on the Sea-Bird SBE911 + CTD system were calibrated in May 2016 at the manufacturer's; for ESAW cruises, pre- and post-cruise calibration have been performed at the Centre for Oceanographic Calibration and Metrology at OGS; for CRELEV2016 cruise the sensors were calibrated between December 2015 and January 2016 at the Sea-Bird facility in Germany (temperature and conductivity) and in the US (oxygen).

The CTD model used during all cruises (i.e. Sea-Bird Electronics, Inc. SBE911 + plus) has the following initial accuracies: conductivity ± 0.0003 S/m, temperature ± 0.001 °C, pressure ± 0.015% of full scale range.

Salinity data were corrected after the cruises using the conductivity measurements obtained from water samples and analysed with a salinometer or performing a laboratory sensor calibration (see subsection “Salinity (measurement and calibration)”). Accuracy for salinometer data is of the order of 0.002 in salinity (precision 0.0001).

We checked for any inconsistencies or anomalies in the thermohaline data using potential temperature vs. salinity diagrams (Fig. 5), which describe the wide range of temperature and salinity values found in the upper and deep layers within the Mediterranean basins. We identified and removed any outliers or spurious data points that deviated from the expected potential temperature vs. salinity relationships.

Salinity in the upper layers of the WMED and of the Adriatic ranges between 37 and 38.5. It is distinctively lower than in the EMED, where surface values are higher (around 39). The saltiest intermediate water is the EIW, its properties decrease in temperature and salinity, from 16 °C and >39 in the EMED, 15.5 °C and 38.9 in the Adriatic, to about 15 °C and 38.8 in the WMED. In the Adriatic Sea a patch of salty but cooler water was observed below the EIW. The deep layers in the three basins are distinguished according to their local potential temperature vs. salinity properties, subject both to advective and local convective processes, related to the dense water formation. Generally, the coolest and freshest are deep WMED waters, and warmest and saltiest, those in the EMED. In the latter, the deep waters of the Cretan Sea are saltier and warmer than those in the Cretan Passage.

**Lowered acoustic doppler current profiler data.** We checked and calibrated the LADCP data following the recommendations by⁵⁶. To maximize LADCP profile accuracy, depth and sound-speed data, which are usually derived from pressure, temperature and salinity measurements, were used during the post-processing. While some ADCPs are equipped with ancillary sensors for measuring pressure and temperature, LADCP systems more typically use data from high-quality CTDs for this purpose. We also used precision navigation from the shipboard GPS receiver, which is crucial for determining the ship's drift during the LADCP cast. The ship drift provides a valuable velocity-referencing constraint, which requires that the deployment/recovery positions are known as precisely as possible. Accuracy of LADCP data is typically between 3 and 4 cm s⁻¹, an estimate based
on comparisons with other independent velocity measurements, such as shipboard ADCP, moored ADCP, or
mechanical current meters. No independent assessment for the data of the cruises presented here has been done.

**Dissolved oxygen data.** Dissolved oxygen data from the CTD sensors have been either calibrated or checked after each cruise using the measurements done with the Winkler titration method (see subsection “Dissolved oxygen (measurement and calibration)”), which is a standard technique for measuring dissolved oxygen concentration in seawater samples and which sees the use of solutions that are calibrated against certified standard solutions. The reported precision that can be associated with dissolved oxygen data when following the described procedure is 2.2 µmol O₂ kg⁻¹. Furthermore, for TalPro2016, a similar method to the one described by⁴⁷ has been used to calibrate the SBE43 data. We matched the data from the upcast profiles with the Winkler samples based on pressure and minimized the residual of the squared differences between the calibrated and water sample DO data. To do this, we adjusted the calibration coefficients of the SBE43 sensor. These coefficients are the DO signal slope, SOC, the voltage at zero DO signal, Voffset, and the pressure factor correction E (see⁴⁷ for details). To ensure high robustness of the DO data adjustment, the number of samples used to determine each set of coefficients was maximised by using the largest possible number of consecutive profiles (i.e. samples with a constant DO_SBE43 vs. DO_Winkler ratio). The total number of samples used to calibrate the SBE43 sensors were 270 for TalPro2016 and 560 for CRELEV2016.

**Nutrient data.** For nutrients we used segmented flow analysis methods with Auto-Analyzers using standard methods⁵⁸,⁶³,⁶⁴. Furthermore, for CRELEV2016, the precision of the measurement was evaluated on five subsamples collected from the same Niskin bottle closed at 125 meters depth. The resulting standard deviations for nitrate + nitrite, phosphate and silicate were 0.011, 0.008, 0.023 µM, respectively. The limit of detection (LOD), based on the standard deviation of the response, σ, and slope, s (LOD = 3σ/s) were 0.04, 0.008 and 0.02 µM for nitrate + nitrite, phosphate and silicate, respectively. For ESAW, the detection limits were 0.006, 0.02, 0.03, 0.01 and 0.01 µM for nitrites, nitrates, ammonium, phosphates and silicates, respectively.

**Dissolved inorganic carbon data.** For DIC, since no calibration unit was available for the system during TalPro2016, samples of CO₂ certified reference material (CRM) from Prof. Dickson (Scripps Institution of Oceanography, San Diego, USA, batch #147) were analysed at the beginning and end of each batch of DIC analysis. The ratio between the CRM DIC nominal and measured value was used as an initial calibration factor, while the CRM DIC measured at the end of the batch was used to correct small temporal drifts. Due to time limitations, no typical precision exercise was done for DIC analysis. Some collected DIC bottles were analysed twice to check the reproducibility. Differences were usually less than 2 µmol kg⁻¹. Considering on board temperature problems (lab temperatures were high, and not stable, as described in subsection “Comments on quality control and internal consistency of CO₂ data”) affecting the conservation of the samples until analysis, and probably also the CO₂ extraction and coulometric reaction, our estimation of the DIC accuracy is 4 µmol kg⁻¹. Also during CRELEV2016, the accuracy of the DIC was assessed by measurements of CRM from Prof. Dickson (Scripps Institution of Oceanography, San Diego, USA). The precision achieved was 2 µmol kg⁻¹, and based on the CRMs used, the accuracy was ± 12 µmol kg⁻¹.

**pH data.** For pH, we used spectrophotometric methods and compared the results with CRM. During TalPro2016, the pH accuracy was checked measuring replicates from CO₂ CRM (Prof. Dickson batch #147). Two exercises were done: 10–12 samples were carefully drawn from each CRM bottle to avoid bubbles. The corresponding theoretical pH25T value for this batch using the dissociation constants from⁹⁷,⁹⁸ is 7.9197. The measured pH25T minus CRM pH25T values were −0.0032 ± 0.0008 and −0.0060 ± 0.0019 pH units. Please note that this information in the cruise report²⁰ is mistaken. No corrections were applied. To verify the precision of the pH measurements, replicate samples from the same Niskin bottle were analyzed twice during the cruise. The standard deviation of those analyses was ± 0.0015 pH units which could be considered as the reproducibility of the pH measurements during the TalPro2016 cruise. During CRELEV2016 and ESAW, the precision of the measurements, evaluated from five determinations on subsamples collected from different Niskin bottles closed at the same depth, were 0.0004 and 0.0008 pH units, respectively.

**Total alkalinity data.** During TalPro2016, CO₂ CRM (batch #147, the same as for DIC and pH analysis) samples were analysed to control the TA accuracy for every batch. As explained in⁴⁴ and in the TalPro2016 cruise report²⁰, the final pH of every batch of analyses was corrected to match the certified CRM TA value on the corresponding analysis. In addition, any temporal drift within each batch of analyses was controlled by titrating filtered and stabilized surface seawater at the beginning and at the end of each batch. TA accuracy is ±2 µmol kg⁻¹. Regarding the precision, due to time limitations, just once during the cruise replicate samples from the same Niskin bottle were titrated, the standard deviation of the 9 replicates was ±1 µmol kg⁻¹. For CRELEV2016, some technical issues occurred at the electrodes of the VINDTA instrument, and many TA data were discarded. For the data that were kept, it has been estimated that the accuracy was ± 4 µmol kg⁻¹, and the precision was <2 µmol kg⁻¹. During ESAW, as estimated by the analysis of certified reference seawater (CRM, provided by A.G. Dickson, Scripps Institution of Oceanography, USA), the accuracy was better than ±4 µmol kg⁻¹, and the precision was ± 2 µmol kg⁻¹.

**Ion carbonate (CO₃²⁻) data.** The accuracy was not checked as the presence of mercury in the CRM bottles interferes with the lead reagent. Regarding the reproducibility, we analysed replicate samples collected from the same Niskin bottle, coming up with a standard deviation of ±6 µmol kg⁻¹ (about ± 3% precision).
**Comments on quality control and internal consistency of CO₂ data.** During TAlPro2016 the laboratory temperature was high, varying from ≈30 °C in the first week to ≈27 °C in the second one, mainly due to calm weather, high air and water temperatures, and an air conditioning failure on board. Therefore, careful attention was paid to the temperature control in the thermostatic baths for DIC, pH and CO₃²⁻. However, we recognise that the final quality of the whole CO₂ system measurements is not ideal, as those in⁴⁴. A first inspection of the CO₂ measurements was done using typical oceanographical criteria, using property-to-property plots, with pressure and salinity being the independent variables for reference. Quality flags were assigned, using just three values, flag 2 for acceptable, 3 for non-acceptable and 9 for not measured or unavailable. An alternative way to check the precision, but not the accuracy, of the CO₂ measurements is inspecting the coherence of the results in waters with a small temporal variability and/or high ventilation time scales, where a higher homogeneity is expected. The former conditions are difficult to find in the Mediterranean Sea where changes in deep and bottom waters were evidenced (and indeed are to be monitored by the Med-SHIP program). During TAlPro2016 we selected samples deeper than 2500 dbar in the Tyrrhenian Sea, where physical and chemical vertical gradients are small and the transient tracer concentrations are low, finding standard deviations of concentrations that compared well with the typical precision analysis from replicate samples (if available) given in the above sections.

The lack of CRM for pH measurements and difficulties in pH metrology to relate measured pH to a traceable International System of Units leads to a confusion between pH accuracy and precision. Consequently, internal consistency studies are usually used for a general estimation of the CO₂ precision and accuracy. Recent publications deal with the lack of consistency between DIC, TA and pH, especially in regions, or waters with low pH⁹⁹. There is no clear consensus on the reasons behind this mismatch⁹⁹,¹⁰⁰. pH method details regarding the equipment, the dye characterization (purified or unpurified), equation for pH calculation and corresponding valid ranges of temperature and salinity. In addition, several CO₂ dissociation constants are available with no clear definition of the best recommended. A recent publication¹⁰¹ recommends⁹⁷ refitted by⁹⁸, or¹⁰² combined with the total boron concentration from¹⁰³, in combination with recommendations in⁴⁴ for the Mediterranean Sea. Using those constants, the pH, DIC and TA mean and standard deviation residuals (measured – calculated) values are respectively, -0.0029 ± 0.0085 pH units, -1.7 ± 5.1 µmol kg⁻¹, 1.9 ± 5.5 µmol kg⁻¹ for the whole water column and -0.0031 ± 0.0062 pH units, -1.8 ± 3.7 µmol kg⁻¹, 1.9 ± 3.9 µmol kg⁻¹ for deep Tyrrhenian Sea waters. Therefore, given the thresholds for consistency in GLODAP¹⁸ (i.e., 0.01–0.02 pH units for pH, and 4 µmol kg⁻¹ for DIC and TA), the TAlPro2016 CO₂ system measurements can be considered internally consistent.

Regarding CO₃²⁻, the original spectrophotometric method by⁷² using PbCl₂ evolved to¹⁰⁴ with a small readjustment in the equations relating absorbance measurements and the final concentration referred to at 25 °C. Later⁷³, changed the reagent to Pb(ClO₄)₂ with new equations. More recently¹⁰⁵, readjusted the equations correcting inaccuracies in the wavelength accuracy of the spectrophotometer equipment. Finally¹⁰⁶, widens the range of salinity and temperature equations validity. Those methods claim to be precise and accurate to about 2% the expected concentration. The scarcity of groups using these methods makes any assessment of a best method choice rather unclear. The CO₃²⁻ concentrations were calculated and reported with⁷³ equations, but the absorbance ratios and reference temperature can be used with more recent equations if required. With the same mentioned CO₂ constants, the mean and standard deviation of measured minus calculated CO₃²⁻ is 35 ± 14 µmol kg⁻¹, with a clear increase of the residuals with increasing calculated CO₃²⁻. Deep Tyrrhenian Sea water residuals are 38 ± 6 µmol kg⁻¹, so the measurements seem to be neither accurate nor precise nor internally consistent with the other CO₂ measurements. Further elaborating on this result is out of the scope of the current manuscript and is commented in a review work about the CO₃²⁻ spectrophotometric technique¹⁰⁷.

**Dichlorodifluoromethane and sulfur hexafluoride data.** Standardization was performed by injecting small volumes of gaseous standard containing SF₆ and CFC-12. This working standard was prepared by the company Dueste-Steiniger (Germany). The concentrations in the standard have been calibrated vs. a reference standard obtained from R.F. Weiss group at Scripps Institution of Oceanography, and the CFC-12 data are reported on the SIO98 scale. One calibration curve was measured to characterize the non-linearity of the system, and point calibrations were always performed between stations to determine the short-term drift of the detector. For TalPro2016, 19 replicate measurements were taken on 11 stations and the determined values for precision are 2.8% for SF₆ and 1.0% for CFC-12. During CRELEV2016 and ESAW2016 no replicate measurements were done, but it is likely that precision is similar, given that samples were measured by the same operators and instruments. During ESAW2015 the Niskin bottles had rubber springs that caused a large contamination for SF₆, so that all those samples had to be discarded.

**Thermosalinograph underway data.** The SBE21 thermosalinograph that was used during CRELEV2016 was calibrated in January 2016. The SBE38 and SBE21 that were used during TAlPro2016 were calibrated in July 2008 and in January 2016, respectively. Initial accuracies for these two models are conductivity ± 0.001 S/m, temperature ± 0.01 °C (SBE21) and ± 0.001 °C (SBE38). The Turner fluorometer that was connected to the system was also calibrated in January 2016. To further validate temperature data, we made a graphic check of the horizontal surface patterns with sea surface temperature data (see Fig. 3), as well as checking for any noise or spikes in the underway data.

**Vessel-mounted acoustic doppler current profiler data.** We checked and calibrated the underway data from ADCP using bottom-tracking methods to correct for ship motion and heading errors. Typical velocity accuracy is ±1.0% of measured velocity, or ±0.5 cm/s, for the Ocean Surveyor, and ±0.5% of measured velocity,
or ±0.5 cm/s, for the WorkHorse Mariner. For CRELEV2016 the ADCP data have been post-processed with the CASCADE software package¹⁰⁸, while for TALPro2016 and ESAW the post-processing has been done with the CODAS Software System (https://currents.soest.hawaii.edu/docs/adcp_doc/), which allows extracting data, assigning coordinates, editing and correcting velocity data. Moreover, the data were corrected for errors in the value of sound velocity in water, and misalignment of the instrument with respect to the axis of the ship (about 0.3° for TALPro2016, 1° for ESAW). After additional editing, to eliminate suspect values, and those where certain thresholds were surpassed, the output files were created by CODAS.

## Usage Notes

The aim of the paper is to provide an exhaustive description of the datasets that were collected during Med-SHIP 2016 cruises. It is thus not the aim to go into detail on oceanographic and biogeochemical aspects that came out from the data analysis. These have been (and will be further in future) the aims of specific scientific papers (e.g.²¹,²²,⁷⁹,¹⁰⁹,¹¹⁰). The reader is invited to refer to those sources to get insight into the quality of a key multidisciplinary database that sets the baseline for future repetitions in the framework of the Med-SHIP initiative.

## Code availability

No code was used to generate this dataset. Figure 5 was plotted by using Ocean Data View software (odv.awi.de), while the others have been generated using MATLAB ©.

Received: 25 January 2023; Accepted: 8 December 2023;

Published online: 09 January 2024

## References

1. Talley, L. D. et al. Changes in Ocean Heat, Carbon Content, and Ventilation: A Review of the First Decade of GO-SHIP Global Repeat Hydrography. *Ann. Rev. Mar. Sci.* 8, 185–215 (2016).

2. Bethoux, J. P. et al. The Mediterranean Sea: A miniature ocean for climatic and environmental studies and a key for the climatic functioning of the North Atlantic. *Prog. Oceanogr.* 44, 131–146 (1999).

3. Wüst, G. On the vertical circulation of the Mediterranean Sea. *J. Geophys. Res.* 66, 3261–3271 (1961).

4. Kinder, T. H. & Bryden, H. L. Gibraltar experiment: summary of the field program and initial results of the Gibraltar experiment. *Gibraltar Exp. Summ. F. Progr. Initial results Gibraltar Exp.* https://doi.org/10.1575/1912/7488 (1988).

5. Malanotte-Rizzoli, P. et al. *The Oceanography Report.* 1, 69 (1988).

6. Cheney, R., Douglas, B., Miller, L. & Milbert, D. *The Oceanography Report.* 67, 8–10 (1986).

7. Lehucher, P. M. et al. Progress from 1989 to 1992 in understanding the circulation of the Western Mediterranean Sea. *Oceanol. Acta* 18, 255–271 (1995).

8. Roether, W. et al. Recent changes in eastern Mediterranean deep waters. *Science* (80-.), 271, 333–335 (1996).

9. Roether, W. & Schlitzer, R. Eastern Mediterranean deep water renewal on the basis of chlorofluoromethane and tritium data. *Dyn. Atmos. Ocean.* 15, 333–354 (1991).

10. Roether, W., Beitzel, V., Sültenfuß, J. & Putzka, A. The Eastern Mediterranean tritium distribution in 1987. *J. Mar. Syst.* 20, 49–61 (1999).

11. Roether, W., Klein, B., Manca, B. B., Theocharis, A. & Kioroglou, S. Transient Eastern Mediterranean deep waters in response to the massive dense-water output of the Aegean Sea in the 1990s. *Prog. Oceanogr.* 74, 540–571 (2007).

12. Roether, W. & Lupton, J. E. Tracers confirm downward mixing of Tyrrhenian Sea upper waters associated with the Eastern Mediterranean Transient. *Ocean Sci.* 7, 91–99 (2011).

13. Key, R. M. et al. A global ocean carbon climatology: Results from Global Data Analysis Project (GLODAP). *Global Biogeochem. Cycles* 18, 1–23 (2004).

14. Sabine, C. L. et al. Global Ocean Data Analysis Project (GLODAP): Results and Data. (2005).

15. CIESM, 2012. Designing Med-SHIP: a Program for repeated oceanographic surveys. N° 43 in CIESM Workshop Monographs. *CIESM Work. Monogr.* 43, (2011).

16. Schroeder, K. et al. Mediterranean Sea Ship-based Hydrographic Investigations Program (Med-SHIP). *Oceanography* 28, 3, https://doi.org/10.5670/oceanog.2015.71 (2015).

17. Hainbucher, D. et al. Physical and biogeochemical parameters of the Mediterranean Sea during a cruise with RV Maria S. Merian in March 2018. *Earth Syst. Sci. Data* 12, 2747–2763 (2020).

18. Olsen, A. et al. An updated version of the global interior ocean biogeochemical data product, GLODAPv2.2020. *Earth Syst. Sci. Data* 12, 3653–3678 (2020).

19. Sloyan, B. M. et al. The global ocean ship-base hydrographic investigations program (GO-SHIP): A platform for integrated multidisciplinary ocean science. *Front. Mar. Sci.* 6, 1–21 (2019).

20. Jullion, L. TALPro2016: A Tyrrhenian Sea & Alger-Provençal component of the MedSHIP Programme, RV Angeles Alvariño, 18/08/16 – 29/08/16, Palermo (Italy) – Barcelona (Spain), Bremerhaven, EUROFLEETS2 Cruise Summary Report, https://epic.awi.de/id/eprint/49725/ (2016).

21. Velaoras, D., Papadopoulos, V. P., Kontoyiannis, H., Cardin, V. & Civitarese, G. Water masses and hydrography during April and June 2016 in the Cretan Sea and Cretan Passage (Eastern Mediterranean Sea). *Deep. Res. Part II Top. Stud. Oceanogr.* 164, 25–40 (2019).

22. Šantić, D. et al. Picoplankton distribution and activity in the deep waters of the Southern Adriatic Sea. *Water (Switzerland)* 11 (2019).

23. Tanhua, T. et al. The Mediterranean Sea system: A review and an introduction to the special issue. *Ocean Sci.* 9, 789–803 (2013).

24. Schroeder, K. TALPro2022 CRUISE REPORT R/V BELGICA Cruise n. 2022/12 (Version 1). *Zenodo*. https://doi.org/10.5281/zenodo.6918731 (2022).

25. Schroeder, K., Chiggiato, J., Bryden, H. L. & Borghini, M. & Ben Ismail, S. Abrupt climate shift in the Western Mediterranean Sea. *Sci. Rep.* 6, 1–7 (2016).

26. Astraldi, M. et al. Water mass properties and chemical signatures in the central Mediterranean region. *J. Mar. Syst.* 33–34, 155–177 (2002).

27. Millot, C. Circulation in the Western Mediterranean Sea. *J. Mar. Syst.* 20, 423–442 (1999).

28. Schroeder, K. et al. Abrupt warming and salting of the Western Mediterranean Deep Water after 2005: Atmospheric forcings and lateral advection. *J. Geophys. Res. Ocean.* 115, (2010).

29. Grignon, L., Smeed, D. A., Bryden, H. L. & Schroeder, K. Importance of the variability of hydrographic preconditioning for deep convection in the Gulf of Lion, NW Mediterranean. *Ocean Sci.* 6, 573–586 (2010).

30. Millot, C. Some features of the Algerian current. *J. Geophys. Res.* 90, 7169–7176 (1985).
31. Millot, C., Taupier-Letage, I. Circulation in the Mediterranean Sea, in: The Mediterranean Sea, the Handbook of Environmental Chemistry, Springer, Berlin/Heidelberg, 29–66 (2005).

32. Manca, B. B., Kovacevic, V., Gai, M. & Viezzoli, D. Dense water formation in the Southern Adriatic Sea and spreading into the Ionian Sea in the period 1997-1999. *J. Mar. Syst.* **33–34**, 133–154 (2002).

33. Gačić, M. *et al*. Salinity in the Sicily Channel corroborates the role of the Adriatic–Ionian Bimodal Oscillating System (BiOS) in shaping the decadal variability of the Mediterranean overturning circulation. *Ocean Sci.* **9**, 83–90 (2013).

34. García-Monteiro, S., Sobrino, J. A., Julien, Y., Sória, G. & Skokovic, D. Surface Temperature trends in the Mediterranean Sea from MODIS data during years 2003–2019. *Reg. Stud. Mar. Sci.* **49**, (2022).

35. Grbec, B. *et al*. Long-Term Trends, Variability and Extremes of *In Situ* Sea Surface Temperature Measured Along the Eastern Adriatic Coast and its Relationship to Hemispheric Processes. *Pure Appl. Geophys.* **175**, 4031–4046 (2018).

36. Mihanović, H. *et al*. Observation, Preconditioning and Recurrence of Exceptionally High Salinities in the Adriatic Sea. *Front. Mar. Sci.* **8**, 1–22 (2021).

37. Klein, B. *et al*. The large deep water transient in the Eastern Mediterranean. *Deep. Res. Part I Oceanogr. Res. Pap.* **46**, 371–414 (1999).

38. Klein, B. *et al*. Is the Adriatic returning to dominate the production of Eastern Mediterranean deep water? *Geophys. Res. Lett.* **33**, 2–5 (2006).

39. Hainbucher, D., Rubino, A. & Klein, B. Water mass characteristics in the deep layers of the western Ionian Basin observed during May 2003. *Geophys. Res. Lett.* **33**, 2–5 (2006).

40. Rubino, A. & Hainbucher, D. A large abrupt change in the abyssal water masses of the eastern Mediterranean. *Geophys. Res. Lett.* **34**, 1–5 (2007).

41. Bensi, M., Rubino, A., Cardin, V., Hainbucher, D. & Mancero-Mosquera, I. Structure and variability of the abyssal water masses in the Ionian Sea in the period 2003–2010. *J. Geophys. Res. Ocean.* **118**, 931–943 (2013).

42. Cardin, V., Civitarese, G., Hainbucher, D., Bensi, M. & Rubino, A. Thermohaline properties in the Eastern Mediterranean in the last three decades: Is the basin returning to the pre-EMT situation? *Ocean Sci.* **11**, 53–66 (2015).

43. Gačić, M. *et al*. Impact of dense-water flow over a sloping bottom on open-sea circulation: Laboratory experiments and an Ionian Sea (Mediterranean) example. *Ocean Sci.* **17**, 975–996 (2021).

44. Schneider, A., Tanhua, T., Roether, W. & Steinfeldt, R. Changes in ventilation of the mediterranean Sea during the past 25 year. *Ocean Sci.* **10**, 1–16 (2014).

45. Álvarez, M. *et al*. The CO2 system in the Mediterranean Sea: A basin wide perspective. *Ocean Sci.* **10**, 69–92 (2014).

46. Schroeder, K. *et al*. Rapid response to climate change in a marginal sea. *Sci. Rep.* **7**, 1–7 (2017).

47. Coppola, L. *et al*. Seasonal and inter-annual variations of dissolved oxygen in the northwestern Mediterranean Sea (DYFAMED site). *Prog. Oceanogr.* **162**, 187–201 (2018).

48. Skejić, S. *et al*. Coccolithophore diversity in open waters of the middle Adriatic Sea in pre- and post-winter periods. *Mar. Micropaleontol.* **143**, 30–45 (2018).

49. Batistić, M. *et al*. Occurrence of winter phytoplankton bloom in the open southern Adriatic: Relationship with hydroclimatic events in the Eastern Mediterranean. *Cont. Shelf Res.* **174**, 12–25 (2019).

50. Jasprica, N. *et al*. Phytoplankton distribution related to different winter conditions in 2016 and 2017 in the open southern Adriatic Sea (eastern Mediterranean). *J. Mar. Syst.* **226**, 103665, https://doi.org/10.1016/j.jimarsys.2021.103665 (2022).

51. Njire, J., Batistić, M., Kovačević, V., Garić, R. & Bensi, M. Tintinnid ciliate communities in pre- and post-winter conditions in the Southern Adriatic Sea (NE Mediterranean). *Water* **11**, 1–19 (2019).

52. Hure, M., Batistić, M., Kovačević, V., Bensi, M. & Garić, R. Copepod community structure in pre-and post-winter conditions in the southern Adriatic sea (NE Mediterranean). *J. Mar. Sci. Eng.* **8**, 1–23 (2020).

53. Tselepides, A., Zervakis, V., Polychronaki, T., Danovaro, R. & Chronis, G. Distribution of nutrients and particulate organic matter in relation to the prevailing hydrographic features of the Cretan Sea (NE Mediterrarean). *Prog. Oceanogr.* **46**, 113–142 (2000).

54. Kessouri, F. *et al*. Vertical mixing effects on phytoplankton dynamics and organic carbon export in the western Mediterranean Sea. *J. Geophys. Res. Ocean.* **123**, 1647–1669 (2018).

55. IOC, SCOR & IAPSO. The international thermodynamic equation of seawater – 2010: Calculation and use of thermodynamic properties. *Intergov. Oceanogr. Comm. Manuals Guid. No.* 56 196 (2010).

56. Turnherr, A. M. *How To Process LADCP Data With the LDEO Software (Versions IX.7 – IX.10).* The “GO-SHIP” manual for LADCP data acquisition ftp://ftp.ldeo.columbia.edu/pub/LADCP/UserManuals (2014).

57. Hood, E.M., Sabine, C.L. & Sloyan, B.M. The GO-SHIP Repeat Hydrography Manual: A Collection of Expert Reports and Guidelines. IOCCP Report Number 14, ICPO Publication Series Number 134. Available online at http://www.go-ship.org/ HydroMan.html (2010).

58. Strickland, J. D. H. & Parsons, T. R. A practical handbook of seawater analysis. *A Pract. Handb. seawater Anal.* (The Alger Press Ltd., 1972).

59. Carpenter, J. H. The Chesapeake Bay Institute technique for the Winkler dissolved oxygen method. *Limnol. Oceanogr.* **10**, 141–143 (1965).

60. Carpenter, J. H. The accuracy of the Winkler method for the dissolved oxygen. *Limnol. Oceanogr.* **10**, 135–140 (1965).

61. Oudot, C., Gerard, R., Morin, P. & Gningue, I. Precise shipboard determination of dissolved oxygen (Winkler procedure) for productivity studies with a commercial system. *Limnol. Oceanogr.* **33**, 146–150 (1988).

62. Kirkwood, D. S. Stability of solutions of nutrient salts during storage. *Mar. Chem.* **38**, 151–164 (1992).

63. Bendschneider, K. & Robinson, R. J. A new spectrophotometric method for the determination of nitrite in sea water. Technical Report No. 8. *University of Washington* (1952).

64. Murphy, R. & Riley, J. P. A Modified Single Solution Method for the Determination of Phosphate in Natural Waters. *Anal. Chim. Acta* **27**, 31–36 (1962).

65. Hansen, H. P. & Koroleff, F. Determination of nutrients. In *Methods of Seawater Analysis* (eds Grasshoff, K., Kremling, K. and Ehrhardt, M.). https://doi.org/10.1002/9783527613984.ch10 (1999).

66. Rimmelin, P. & Moutin, T. Re-examination of the MAGIC method to determine low orthophosphate concentration in seawater. *Anal. Chim. Acta* **548**, 174–182 (2005).

67. Johnson, K. M., Wills, K. D., Butler, D. B., Johnson, W. K. & Wong, C. S. Coulometric total carbon dioxide analysis for marine studies: maximizing the performance of an automated gas extraction system and coulometric detector. *Mar. Chem.* **44**, 167–187 (1993).

68. Clayton, T. D. & Byrne, H. B. Spectrophotometric seawater pH measurements: total hydrogen ion concentration scale calibration of m-cresol purple and at-sea results. *Deep. Res. Part I Oceanogr. Res. Pap.* **40**, 2115–2129 (1993).

69. Pérez, F. F. & Fraga, F. A precise and rapid analytical procedure for alkalinity determination. *Mar. Chem.* **21**, 169–182 (1987).

70. Pérez, F. F., Ríos, A. F., Rellán, T. & Álvarez, M. Improvements in a fast potentiometric seawater alkalinity determination. *Ciencias Mar.* **26**, 463–478 (2000).

71. Mintrop, L., Pérez, F. F., González-Dávila, M., Santana-Casiano, J. M. & Körtzinger, A. Alkalinity determination by potentiometry: Intercalibration using three different methods. *Ciencias Mar.* **26**, 23–37 (2000).

72. Byrne, R. H. & Yao, W. Procedures for measurement of carbonate ion concentrations in seawater by direct spectrophotometric observations of Pb(II) complexation. *Mar. Chem.* **112**, 128–135 (2008).
73. Patsavas, M. C. et al. Procedures for direct spectrophotometric determination of carbonate ion concentrations: Measurements in US Gulf of Mexico and East Coast waters. *Mar. Chem.* **168**, 80–85 (2015).

74. Dickson, A. G., Sabine, C. L. & Christian, J. R. Guide to best practices for ocean CO2 measurements. *PICES Spec. Publ.* **3**, 191 (2007).

75. Hartman, S.E., Dumousseaud, C., Roberts, A. Operating manual for the Marianda (Versatile Instrument for the Determination of Titration Alkalinity) VINDTA 3C for the laboratory based determination of Total and Total Dissolved Inorganic Carbon in seawater. *National Oceanography Center. Internal Document 2011*, pp66 (2011).

76. Liu, X., Patsavas, M. C. & Byrne, R. H. Purification and characterization of meta-cresol purple for spectrophotometric seawater pH measurements. *Environ. Sci. Technol.* **45**, 4862–4868 (2011).

77. Gattuso, J. P., Lee, K., Rost, B. & Schulz, K. G. Approaches and tools to manipulate the carbonate chemistry. Guide for Best Practices in Ocean Acidification Research and Data Reporting. In: Riebesell, U., Fabry, V. J., Hansson, L., Gattuso, J.-P. (Eds.), (2010).

78. Stöven, T. & Tanhua, T. Ventilation of the mediterranean sea constrained by multiple transient tracer measurements. *Ocean Sci.* **10**, 439–457 (2014).

79. Li, P. & Tanhua, T. Recent Changes in Deep Ventilation of the Mediterranean Sea; Evidence From Long-Term Transient Tracer Observations. *Front. Mar. Sci.* **7**, 1–23 (2020).

80. Bullister, J. L. & Weiss, R. F. Determination of CCl3F and CCl2F in seawater and air. Deep Sea Res. Part A, Oceanogr. Res. Pap. **35**, 839–853 (1988).

81. Bullister, J. L. & Wisegarver, D. P. The shipboard analysis of trace levels of sulfur hexafluoride, chlorofluorocarbon-11 and chlorofluorocarbon-12 in seawater. *Deep. Res. Part I Oceanogr. Res. Pap.* **55**, 1063–1074 (2008).

82. Stöven, T. Ventilation processes of the Mediterranean Sea based on CFC-12 and SF6 measurements. *Ifm-Geomar Diploma th*, **54** (2011).

83. Tanhua, T. Physical oceanography during MedSHIP cruise Talpro. *PANGAEA*. https://doi.org/10.1594/PANGAEA.902330 (2019a).

84. Tanhua, T. Hydrochemistry of water samples during MedSHIP cruise Talpro. *PANGAEA*. https://doi.org/10.1594/PANGAEA.902293 (2019b).

85. Schroeder, K. LADCP current profiles collected during MedSHIP cruise TALPro2016. *PANGAEA*. https://doi.org/10.1594/PANGAEA.932316 (2021a).

86. Schroeder, K. Vessel-mounted ADCP data (current profiles) collected during MedSHIP cruise TALPro2016. *PANGAEA*. https://doi.org/10.1594/PANGAEA.933628 (2021b).

87. Schroeder, K. Continuous thermosalinograph oceanography along MedSHIP cruise TALPro2016 on RV *Angeles Alvariño*. *PANGAEA*. https://doi.org/10.1594/PANGAEA.946435 (2022).

88. Jullion, L., Álvarez, M., Tanhua, T.: Dissolved inorganic carbon (DIC), total alkalinity, pH on total scale, chlorofluorocarbon 12 (CFC-12), temperature, salinity and other measurements collected from discrete samples and profile observations during the R/V *Ángeles Alvarino* TALPRO cruise (EXPOCODE 29AJ20160818) in the Mediterranean Sea from 2016-08-18 to 2016-08-29 (NCEI Accession 0213116). *NOAA National Centers for Environmental Information*. https://doi.org/10.25921/8f6a-zq94 (2020).

89. Civitarese, G., Álvarez, M., Tanhua, T. Dissolved inorganic carbon (DIC), total alkalinity, pH on total scale, chlorofluorocarbons (CFC-11, CFC-12), water temperature, salinity, dissolved oxygen, nutrients and other measurements collected from discrete samples and profile observations during the R/V *Aegaeo* CRELEV cruise (EXPOCODE 36AE20160602) in the Aegean Sea and Levantine Basin from 2016-06-02 to 2016-06-10 (NCEI Accession 0217938). *NOAA National Centers for Environmental Information*. https://doi.org/10.25921/c01n-sv64 (2020).

90. Kovačević, V., Bensi, M., Muslim, S., Mihanović H. Physical oceanography of the ESAW Leg 1 cruise in 2015. *NODC-OGS*. https://doi.org/10.13120/A0ACDE7F-AF65-4633-A021-7D15F52BC434 (2021a).

91. Kovačević, V., Giani, M. & Santinelli, C. Hydrochemistry of the water samples collected during ESAW leg 1 cruise in 2015. *NODC-OGS*. https://doi.org/10.13120/p0e3-3m55 (2021b).

92. Kovačević, V., Bensi, M., Muslim, S., Mihanović, H. & Matic, F. Vessel-mounted ADCP data collected during ESAW leg 1 cruise in 2015. *NODC-OGS*. https://doi.org/10.13120/61ee5a96-359f-42ab-a42c-5f92b40d41a3 (2022a).

93. Kovačević, V., Bensi, M., Muslim, S., Mihanović, H., Matic F. Physical oceanography of the ESAW Leg 2 cruise in 2016. *NODC-OGS*. https://doi.org/10.13120/2ffaa5c1-0f52-4168-80ab-0043d1ed38cf (2022b).

94. Kovačević, V., Giani, M. & Santinelli, C. Hydrochemistry of the water samples collected during ESAW leg 2 cruise in 2016. *NODC-OGS*. https://doi.org/10.13120/f938681f-6062-4a99-8125-cce9afd929db (2022c).

95. Kovačević, V., Ursella, L. & Mihanović, H. Vessel-mounted ADCP data collected during ESAW leg 2 cruise in 2016. *NODC-OGS*. https://doi.org/10.13120/a4881c9-7e9f-46fe-8276-981f83521d2f (2022d).

96. Civitarese, G. et al. CRELEV 2016. *NODC-OGS*. https://doi.org/10.13120/9a8c8fe1-5438-4854-8a57-a13339a654b0 (2022).

97. Mehrbach, C., Culberson, C. H., Hawley, J. E. & Pytkowicz, R. M. Measurement of the Apparent Dissociation Constants of Carbonic Acid in Seawater At Atmospheric Pressure. *Limnol. Oceanogr.* **18**, 897–907 (1973).

98. Dickson, A. G. & Millero, F. J. A comparison of the equilibrium constants for the dissociation of carbonic acid in seawater media. *Deep Sea Res. Part A, Oceanogr. Res. Pap.* **34**, 1733–1743 (1987).

99. Álvarez, M. et al. Global Ocean Spectrophotometric pH Assessment: Consistent Inconsistencies. *Environ. Sci. Technol.* **54**, 10977–10988 (2020).

100. Takeshita, Y. et al. Assessment of pH dependent errors in spectrophotometric pH measurements of seawater. *Mar. Chem.* **223**, 103801 (2020).

101. Woosley, R. J. Evaluation of the temperature dependence of dissociation constants for the marine carbon system using pH and certified reference materials. *Mar. Chem.* **229**, 103914 (2021).

102. Lueker, T. J., Dickson, A. G. & Keeling, C. D. Ocean pCO2 calculated from DIC, TA, and the Mehrbach equations for K1 and K2: Validation using laboratory measurements of CO2 in gas and seawater at equilibrium. *Mar. Chem.* **70**, 105–119 (2000).

103. Lee, K. et al. The universal ratio of boron to chlorinity for the North Pacific and North Atlantic oceans. *Geochim. Cosmochim. Acta* **74**, 1801–1811 (2010).

104. Easley, R. A. et al. Spectrophotometric measurement of calcium carbonate saturation states in seawater. *Environ. Sci. Technol.* **47**, 1468–1477 (2013).

105. Sharp, J. D. et al. Spectrophotometric Determination of Carbonate Ion Concentrations: Elimination of Instrument-Dependent Offsets and Calculation of *in Situ* Saturation States. *Environ. Sci. Technol.* **51**, 9127–9136 (2017).

106. Sharp, J. D. & Byrne, R. H. Carbonate ion concentrations in seawater: Spectrophotometric determination at ambient temperatures and evaluation of propagated calculation uncertainties. *Mar. Chem.* **209**, 70–80 (2019).

107. Guallart, F. et al. Spectrophotometric Measurement of Carbonate Ion in Seawater over a Decade: Dealing with Inconsistencies. *Environ. Sci. Technol.* **56**, 7381–7395 (2022).

108. Le Bot, P., Kermabon, C., Lherminier, P. & Gaillard, F. Chaine Automatisée de Suivi des Courantomètres Acoustiques Doppler Embarqués. CASCADE V6.1: Logiciel de validation et de visualisation des mesures ADCP de coque. *Rapp. OPS/LPO 11-01* **93** (2011).

109. Protopapa, M., Zervoudaki, S., Assimakopoulou, G., Velaoras, D. & Koppelmann, R. Mesozooplankton community structure in the Eastern Mediterranean Sea. *J. Mar. Syst.* **211**, 103401 (2020).

110. Belgacem, M. et al. Dissolved inorganic nutrients in the western Mediterranean Sea (2004–2017). *Earth Syst. Sci. Data* **12**, 1985–2011 (2020).

111. Schroeder K. *Current Systems in the Mediterranean Sea, In: Encyclopedia of Ocean Sciences - 3: Ocean Dynamics*, edited by: Cochran J.K., Bokuniewicz H.J., Yager P.L., Elsevier, https://doi.org/10.1016/B978-0-12-409548-9.11296-5 (2019).
## Acknowledgements

We are grateful to the captains, officers and crews of the RVs ANGELES ALVARIÑO, AEGAEO and BIOS DVA for their relentless effort and their dedication. These cruises would not have been successful without their support. We acknowledge the Grant agreement no. 312762 of the EUROFLEETS2 Project funded by the European Commission, which made the ship time available to conduct these cruises. For the DOC and biogeochemistry data collection during ESAW cruises thanks go also to Stefano Vestri and Yuri Galletti (from CNR-IBF, Pisa, Italy), Cinzia Fabbro and Walter Dellisanti (from OGS, Trieste, Italy). An internal IEO project supported the Inorganic Chemical Oceanography (INOCEN) group (M. Álvarez, M. Castaño-Carrera & I. Rodríguez) participation in TALPro2016. Finally, the authors dedicate this work to the memory of their friend and colleague, and pioneering oceanographer of the Mediterranean Sea, Wolfgang Roether (1935–2022). His work and dedication to understand the circulation, ventilation and oceanography of the Mediterranean Sea has been a source of great inspiration for the authors and for their dedication for systematic oceanographic investigations in the Mediterranean Sea.

## Author contributions

All authors participated in at least one of the cruises mentioned in the text or to the sample analyses in the onshore labs. K.S. and T.T. elaborated the concept of the paper, K.S. prepared the draft of the paper, with the help of T.T., M.A., V.K., M.B., D.V., G.C. and L.U. All other authors contributed to the parts that are devoted to the specific parameters they have collected, processed, analysed or interpreted.

## Competing interests

The authors declare no competing interests.

## Additional information

Correspondence and requests for materials should be addressed to K.S.

Reprints and permissions information is available at www.nature.com/reprints.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_12.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024