The DTM is a compilation of multibeam echosounder surveys at Central Aegean Sea - around Santorini Island (Greece), acquired with R/V Aegaeo using Seabeam 2120 (20 kHz) echosounder during 2001, 2006, 2016 and 2017 cruises. Grid resolution is 1/32 arc-minutes processed for the purpose of the HR DTMs layer of EMODnet Bathymetry HRSM, October 2018.

The DTM is a compilation of multibeam echosounder surveys at Korinthiakos Gulf (Greece) acquired with R/V Aegaeo using Seabeam 2120 (20 kHz) echosounder from 2001 to 2004. Grid resolution is 1/32 arc-minutes processed for the purpose of the HR DTMs layer of EMODnet Bathymetry HRSM, October 2018.

Moving 10-years analysis of nitrite at Mediterranean Sea for each season : - winter (January-March), - spring (April-June), - summer (July-September), - autumn (October-December). Every year of the time dimension corresponds to the 10-year centered average of each season. Decades span from 1965-1974 until 2004-2013. Observational data span from 1961 to 2013. Depth range (IODE standard depths): -1500.0, -1400.0, -1300.0, -1200.0, -1100.0, -1000.0, -900.0, -800.0, -700.0, -600.0, -500.0, -400.0, -300.0, -250.0, -200.0, -150.0, -125.0, -100.0, -75.0, -50.0, -30.0, -20.0, -10.0, -5.0, -0.0. Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. Profiles were interpolated at standard depths using weighted parabolic interpolation algorithm (Reiniger and Ross, 1968). GEBCO 1min topography is used for the contouring preparation. Analyzed filed masked using relative error threshold 0.3 and 0.5. DIVA settings: A constant value for signal-to-noise ratio was used equal to 3. Correlation length was optimized and filtered vertically and a seasonally-averaged profile was used. Logarithmic transformation applied to the data prior to the analysis. Background field: the data mean value is subtracted from the data. Detrending of data: no. Advection constraint applied: no. Originators of Italian data sets-List of contributors: o Brunetti Fabio (OGS) o Cardin Vanessa, Bensi Manuel doi:10.6092/36728450-4296-4e6a-967d-d5b6da55f306 o Cardin Vanessa, Bensi Manuel, Ursella Laura, Siena Giuseppe doi:10.6092/f8e6d18e-f877-4aa5-a983-a03b06ccb987 o Cataletto Bruno (OGS) o Cinzia Comici Cinzia (OGS) o Civitarese Giuseppe (OGS) o DeVittor Cinzia (OGS) o Giani Michele (OGS) o Kovacevic Vedrana (OGS) o Mosetti Renzo (OGS) o Solidoro C.,Beran A.,Cataletto B.,Celussi M.,Cibic T.,Comici C.,Del Negro P.,De Vittor C.,Minocci M.,Monti M.,Fabbro C.,Falconi C.,Franzo A.,Libralato S.,Lipizer M.,Negussanti J.S.,Russel H.,Valli G., doi:10.6092/e5518899-b914-43b0-8139-023718aa63f5 o Celio Massimo (ARPA FVG) o Malaguti Antonella (ENEA) o Fonda Umani Serena (UNITS) o Bignami Francesco (ISAC/CNR) o Boldrini Alfredo (ISMAR/CNR) o Marini Mauro (ISMAR/CNR) o Miserocchi Stefano (ISMAR/CNR) o Zaccone Renata (IAMC/CNR) o Lavezza, R., Dubroca, L. F. C., Ludicone, D., Kress, N., Herut, B., Civitarese, G., Cruzado, A., Lefèvre, D., Souvermezoglou, E., Yilmaz, A., Tugrul, S., and Ribera d'Alcala, M.: Compilation of quality controlled nutrient profiles from the Mediterranean Sea, doi:10.1594/PANGAEA.771907, 2011.

Moving 10-years analysis of ammonium at Mediterranean Sea for each season : - winter (January-March), - spring (April-June), - summer (July-September), - autumn (October-December). Every year of the time dimension corresponds to the 10-year centered average of the winter season. Decades span from 1977-1986 until 2004-2013. Observational data span from 1970 to 2013. Depth range (IODE standard depths): -1000.0, -900.0, -800.0, -700.0, -600.0, -500.0, -400.0, -300.0, -250.0, -200.0, -150.0, -125.0, -100.0, -75.0, -50.0, -30.0, -20.0, -10.0, -5.0, -0.0. Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. Profiles were interpolated at standard depths using weighted parabolic interpolation algorithm (Reiniger and Ross, 1968). GEBCO 1min topography is used for the contouring preparation. Analyzed filed masked using relative error threshold 0.3 and 0.5. DIVA settings: A constant value for signal-to-noise ratio was used equal to 3. Correlation length was optimized and filtered vertically and a seasonally-averaged profile was used. Logarithmic transformation applied to the data prior to the analysis. Background field: the data mean value is subtracted from the data. Detrending of data: no. Advection constraint applied: no. Originators of Italian data sets-List of contributors: o Brunetti Fabio (OGS) o Cardin Vanessa, Bensi Manuel doi:10.6092/36728450-4296-4e6a-967d-d5b6da55f306 o Cardin Vanessa, Bensi Manuel, Ursella Laura, Siena Giuseppe doi:10.6092/f8e6d18e-f877-4aa5-a983-a03b06ccb987 o Cataletto Bruno (OGS) o Cinzia Comici Cinzia (OGS) o Civitarese Giuseppe (OGS) o DeVittor Cinzia (OGS) o Giani Michele (OGS) o Kovacevic Vedrana (OGS) o Mosetti Renzo (OGS) o Solidoro C.,Beran A.,Cataletto B.,Celussi M.,Cibic T.,Comici C.,Del Negro P.,De Vittor C.,Minocci M.,Monti M.,Fabbro C.,Falconi C.,Franzo A.,Libralato S.,Lipizer M.,Negussanti J.S.,Russel H.,Valli G., doi:10.6092/e5518899-b914-43b0-8139-023718aa63f5 o Celio Massimo (ARPA FVG) o Malaguti Antonella (ENEA) o Fonda Umani Serena (UNITS) o Bignami Francesco (ISAC/CNR) o Boldrini Alfredo (ISMAR/CNR) o Marini Mauro (ISMAR/CNR) o Miserocchi Stefano (ISMAR/CNR) o Zaccone Renata (IAMC/CNR) o Lavezza, R., Dubroca, L. F. C., Ludicone, D., Kress, N., Herut, B., Civitarese, G., Cruzado, A., Lefèvre, D., Souvermezoglou, E., Yilmaz, A., Tugrul, S., and Ribera d'Alcala, M.: Compilation of quality controlled nutrient profiles from the Mediterranean Sea, doi:10.1594/PANGAEA.771907, 2011.

Moving 10-years analysis of phosphates at Mediterranean Sea for each season : - winter (January-March), - spring (April-June), - summer (July-September), - autumn (October-December). Every year of the time dimension corresponds to the 10-year centered average of each season. Decades span from 1960-1969 until 2004-2013. Observational data span from 1960 to 2013. Depth range (IODE standard depths): -1500.0, -1400.0, -1300.0, -1200.0, -1100.0, -1000.0, -900.0, -800.0, -700.0, -600.0, -500.0, -400.0, -300.0, -250.0, -200.0, -150.0, -125.0, -100.0, -75.0, -50.0, -30.0, -20.0, -10.0, -5.0, -0.0. Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. Profiles were interpolated at standard depths using weighted parabolic interpolation algorithm (Reiniger and Ross, 1968). GEBCO 1min topography is used for the contouring preparation. Analyzed filed masked using relative error threshold 0.3 and 0.5. DIVA settings: A constant value for signal-to-noise ratio was used equal to 3. Correlation length was optimized and filtered vertically and a seasonally-averaged profile was used. Logarithmic transformation applied to the data prior to the analysis. Background field: the data mean value is subtracted from the data. Detrending of data: no. Advection constraint applied: no. Originators of Italian data sets-List of contributors • Brunetti Fabio (OGS) • Cardin Vanessa, Bensi Manuel doi:10.6092/36728450-4296-4e6a-967d-d5b6da55f306 • Cardin Vanessa, Bensi Manuel, Ursella Laura, Siena Giuseppe doi:10.6092/f8e6d18e-f877-4aa5-a983-a03b06ccb987 • Cataletto Bruno (OGS) • Cinzia Comici Cinzia (OGS) • Civitarese Giuseppe (OGS) • DeVittor Cinzia (OGS) • Giani Michele (OGS) • Kovacevic Vedrana (OGS) • Mosetti Renzo (OGS) • Solidoro C.,Beran A.,Cataletto B.,Celussi M.,Cibic T.,Comici C.,Del Negro P.,De Vittor C.,Minocci M.,Monti M.,Fabbro C.,Falconi C.,Franzo A.,Libralato S.,Lipizer M.,Negussanti J.S.,Russel H.,Valli G., doi:10.6092/e5518899-b914-43b0-8139-023718aa63f5 • Celio Massimo (ARPA FVG) • Malaguti Antonella (ENEA) • Fonda Umani Serena (UNITS) • Bignami Francesco (ISAC/CNR) • Boldrini Alfredo (ISMAR/CNR) • Marini Mauro (ISMAR/CNR) • Miserocchi Stefano (ISMAR/CNR) • Zaccone Renata (IAMC/CNR) • Lavezza, R., Dubroca, L. F. C., Ludicone, D., Kress, N., Herut, B., Civitarese, G., Cruzado, A., Lefèvre, D., Souvermezoglou, E., Yilmaz, A., Tugrul, S., and Ribera d’Alcala, M.: Compilation of quality controlled nutrient profiles from the Mediterranean Sea, doi:10.1594/PANGAEA.771907, 2011. Units: umol/l

Moving 10-years winter (January-March) analysis of total nitrogen at Mediterranean Sea for each season : - winter (January-March), - spring (April-June), - summer (July-September), - autumn (October-December). Every year of the time dimension corresponds to the 10-year centered average of each season. Decades span from 1990-1999 until 1995-2004. Observational data span from 1979 to 2013. Depth range (IODE standard depths): -1000.0, -900.0, -800.0, -700.0, -600.0, -500.0, -400.0, -300.0, -250.0, -200.0, -150.0, -125.0, -100.0, -75.0, -50.0, -30.0, -20.0, -10.0, -5.0, -0.0. Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. Profiles were interpolated at standard depths using weighted parabolic interpolation algorithm (Reiniger and Ross, 1968). GEBCO 1min topography is used for the contouring preparation. Analyzed filed masked using relative error threshold 0.3 and 0.5. DIVA settings: A constant value for signal-to-noise ratio was used equal to 3. Correlation length was optimized and filtered vertically and a seasonally-averaged profile was used. Logarithmic transformation applied to the data prior to the analysis. Background field: the data mean value is subtracted from the data. Detrending of data: no. Advection constraint applied: no. Originators of Italian data sets-List of contributors: o Brunetti Fabio (OGS) o Cardin Vanessa, Bensi Manuel doi:10.6092/36728450-4296-4e6a-967d-d5b6da55f306 o Cardin Vanessa, Bensi Manuel, Ursella Laura, Siena Giuseppe doi:10.6092/f8e6d18e-f877-4aa5-a983-a03b06ccb987 o Cataletto Bruno (OGS) o Cinzia Comici Cinzia (OGS) o Civitarese Giuseppe (OGS) o DeVittor Cinzia (OGS) o Giani Michele (OGS) o Kovacevic Vedrana (OGS) o Mosetti Renzo (OGS) o Solidoro C.,Beran A.,Cataletto B.,Celussi M.,Cibic T.,Comici C.,Del Negro P.,De Vittor C.,Minocci M.,Monti M.,Fabbro C.,Falconi C.,Franzo A.,Libralato S.,Lipizer M.,Negussanti J.S.,Russel H.,Valli G., doi:10.6092/e5518899-b914-43b0-8139-023718aa63f5 o Celio Massimo (ARPA FVG) o Malaguti Antonella (ENEA) o Fonda Umani Serena (UNITS) o Bignami Francesco (ISAC/CNR) o Boldrini Alfredo (ISMAR/CNR) o Marini Mauro (ISMAR/CNR) o Miserocchi Stefano (ISMAR/CNR) o Zaccone Renata (IAMC/CNR) o Lavezza, R., Dubroca, L. F. C., Ludicone, D., Kress, N., Herut, B., Civitarese, G., Cruzado, A., Lefèvre, D., Souvermezoglou, E., Yilmaz, A., Tugrul, S., and Ribera d'Alcala, M.: Compilation of quality controlled nutrient profiles from the Mediterranean Sea, doi:10.1594/PANGAEA.771907, 2011.

Moving 10-years analysis of silicates at Mediterranean Sea for each season : - winter (January-March), - spring (April-June), - summer (July-September), - autumn (October-December). Every year of the time dimension corresponds to the 10-year centered average of each season. Decades span from 1960-1969 until 2004-2013. Observational data span from 1960 to 2013. Depth range (IODE standard depths): -1500.0, -1400.0, -1300.0, -1200.0, -1100.0, -1000.0, -900.0, -800.0, -700.0, -600.0, -500.0, -400.0, -300.0, -250.0, -200.0, -150.0, -125.0, -100.0, -75.0, -50.0, -30.0, -20.0, -10.0, -5.0, -0.0. Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. Profiles were interpolated at standard depths using weighted parabolic interpolation algorithm (Reiniger and Ross, 1968). GEBCO 1min topography is used for the contouring preparation. Analyzed filed masked using relative error threshold 0.3 and 0.5. DIVA settings: A constant value for signal-to-noise ratio was used equal to 3. Correlation length was optimized and filtered vertically and a seasonally-averaged profile was used. Logarithmic transformation applied to the data prior to the analysis. Background field: the data mean value is subtracted from the data. Detrending of data: no. Advection constraint applied: no. Originators of Italian data sets-List of contributors • Brunetti Fabio (OGS) • Cardin Vanessa, Bensi Manuel doi:10.6092/36728450-4296-4e6a-967d-d5b6da55f306 • Cardin Vanessa, Bensi Manuel, Ursella Laura, Siena Giuseppe doi:10.6092/f8e6d18e-f877-4aa5-a983-a03b06ccb987 • Cataletto Bruno (OGS) • Cinzia Comici Cinzia (OGS) • Civitarese Giuseppe (OGS) • DeVittor Cinzia (OGS) • Giani Michele (OGS) • Kovacevic Vedrana (OGS) • Mosetti Renzo (OGS) • Solidoro C.,Beran A.,Cataletto B.,Celussi M.,Cibic T.,Comici C.,Del Negro P.,De Vittor C.,Minocci M.,Monti M.,Fabbro C.,Falconi C.,Franzo A.,Libralato S.,Lipizer M.,Negussanti J.S.,Russel H.,Valli G., doi:10.6092/e5518899-b914-43b0-8139-023718aa63f5 • Celio Massimo (ARPA FVG) • Malaguti Antonella (ENEA) • Fonda Umani Serena (UNITS) • Bignami Francesco (ISAC/CNR) • Boldrini Alfredo (ISMAR/CNR) • Marini Mauro (ISMAR/CNR) • Miserocchi Stefano (ISMAR/CNR) • Zaccone Renata (IAMC/CNR) • Lavezza, R., Dubroca, L. F. C., Ludicone, D., Kress, N., Herut, B., Civitarese, G., Cruzado, A., Lefèvre, D., Souvermezoglou, E., Yilmaz, A., Tugrul, S., and Ribera d’Alcala, M.: Compilation of quality controlled nutrient profiles from the Mediterranean Sea, doi:10.1594/PANGAEA.771907, 2011. Units: umol/l

EMODnet Chemistry aims to provide access to marine chemistry data sets and derived data products concerning eutrophication, ocean acidification and contaminants. The chemicals chosen reflect importance to the Marine Strategy Framework Directive (MSFD). This regional aggregated dataset contains all unrestricted EMODnet Chemistry data on contaminants (24 parameters), and covers the Mediterranean Sea with 4517 CDI records divided per matrices: 520 biota profiles, 560 water profiles, 3437 sediment (26 Vertical profiles and 3411 Time series). In the water dataset, the vertical profiles temporal range is from 1974-09-12 to 2015-09-14. In sediment dataset, vertical profiles temporal range is from 2010-08-02 to 2014-09-28 and time series temporal range is from 1981-06-27 to 2018-08-02. In biota time series temporal range is from 1979-03-29 to 2017-03-15. Data were aggregated and quality controlled by ‘Hellenic Centre for Marine Research, Hellenic National Oceanographic Data Centre (HCMR/HNODC)’ from Greece. Regional datasets concerning contaminants are automatically harvested. Parameter names in these datasets are based on P01, BODC Parameter Usage Vocabulary, which is available at: http://seadatanet.maris2.nl/bandit/browse_step.php . Each measurement value has a quality flag indicator. The resulting data collections for each Sea Basin are harmonised, and the collections are quality controlled by EMODnet Chemistry Regional Leaders using ODV Software and following a common methodology for all Sea Regions. Harmonisation means that: (1) unit conversion is carried out to express contaminant concentrations with a limited set of measurement units (according to EU directives 2013/39/UE; Comm. Dec. EU 2017/848) and (2) merging of variables described by different “local names” ,but corresponding exactly to the same concepts in BODC P01 vocabulary. The harmonised dataset can be downloaded as ODV spreadsheet (TXT file), which is composed of metadata header followed by tab separated values. This worksheet can be imported to ODV Software for visualisation (More information can be found at: https://www.seadatanet.org/Software/ODV ). The same dataset is offered also as XLSX file in a long/vertical format, in which each P01 measurement is a record line. Additionally, there are a series of columns that split P01 terms in subcomponents (measure, substance, CAS number, matrix...).This transposed format is more adapted to worksheet applications users (e.g. LibreOffice Calc). The 24 parameter names in this metadata record are based on P02, SeaDataNet Parameter Discovery Vocabulary, which is available at: http://seadatanet.maris2.nl/v_bodc_vocab_v2/vocab_relations.asp?lib=P02 . Detailed documentation will be published soon. The original datasets can be searched and downloaded from EMODnet Chemistry Download Service: https://emodnet-chemistry.maris.nl/search

Moving 10-years analysis of total phosphorus at Mediterranean Sea for each season : - winter (January-March), - spring (April-June), - summer (July-September), - autumn (October-December). Every year of the time dimension corresponds to the 10-year centered average of each season. Decades span from 1984-1993 until 1998-2007. Observational data span from 1978 to 2013. Depth range (IODE standard depths): -1000.0, -900.0, -800.0, -700.0, -600.0, -500.0, -400.0, -300.0, -250.0, -200.0, -150.0, -125.0, -100.0, -75.0, -50.0, -30.0, -20.0, -10.0, -5.0, -0.0. Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. Profiles were interpolated at standard depths using weighted parabolic interpolation algorithm (Reiniger and Ross, 1968). GEBCO 1min topography is used for the contouring preparation. Analyzed filed masked using relative error threshold 0.3 and 0.5. DIVA settings: A constant value for signal-to-noise ratio was used equal to 3. Correlation length was optimized and filtered vertically and a seasonally-averaged profile was used. Logarithmic transformation applied to the data prior to the analysis. Background field: the data mean value is subtracted from the data. Detrending of data: no. Advection constraint applied: no. Originators of Italian data sets-List of contributors: o Brunetti Fabio (OGS) o Cardin Vanessa, Bensi Manuel doi:10.6092/36728450-4296-4e6a-967d-d5b6da55f306 o Cardin Vanessa, Bensi Manuel, Ursella Laura, Siena Giuseppe doi:10.6092/f8e6d18e-f877-4aa5-a983-a03b06ccb987 o Cataletto Bruno (OGS) o Cinzia Comici Cinzia (OGS) o Civitarese Giuseppe (OGS) o DeVittor Cinzia (OGS) o Giani Michele (OGS) o Kovacevic Vedrana (OGS) o Mosetti Renzo (OGS) o Solidoro C.,Beran A.,Cataletto B.,Celussi M.,Cibic T.,Comici C.,Del Negro P.,De Vittor C.,Minocci M.,Monti M.,Fabbro C.,Falconi C.,Franzo A.,Libralato S.,Lipizer M.,Negussanti J.S.,Russel H.,Valli G., doi:10.6092/e5518899-b914-43b0-8139-023718aa63f5 o Celio Massimo (ARPA FVG) o Malaguti Antonella (ENEA) o Fonda Umani Serena (UNITS) o Bignami Francesco (ISAC/CNR) o Boldrini Alfredo (ISMAR/CNR) o Marini Mauro (ISMAR/CNR) o Miserocchi Stefano (ISMAR/CNR) o Zaccone Renata (IAMC/CNR) o Lavezza, R., Dubroca, L. F. C., Ludicone, D., Kress, N., Herut, B., Civitarese, G., Cruzado, A., Lefèvre, D., Souvermezoglou, E., Yilmaz, A., Tugrul, S., and Ribera d'Alcala, M.: Compilation of quality controlled nutrient profiles from the Mediterranean Sea, doi:10.1594/PANGAEA.771907, 2011.

Moving 10-years analysis of nitrates at Mediterranean Sea for each season : - winter (January-March), - spring (April-June), - summer (July-September), - autumn (October-December). Every year of the time dimension corresponds to the 10-year centered average of each season. Decades span from 1960-1969 until 2004-2013. Observational data span from 1960 to 2013. Depth range (IODE standard depths): -1500.0, -1400.0, -1300.0, -1200.0, -1100.0, -1000.0, -900.0, -800.0, -700.0, -600.0, -500.0, -400.0, -300.0, -250.0, -200.0, -150.0, -125.0, -100.0, -75.0, -50.0, -30.0, -20.0, -10.0, -5.0, -0.0. Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. Profiles were interpolated at standard depths using weighted parabolic interpolation algorithm (Reiniger and Ross, 1968). GEBCO 1min topography is used for the contouring preparation. Analyzed filed masked using relative error threshold 0.3 and 0.5. DIVA settings: A constant value for signal-to-noise ratio was used equal to 3. Correlation length was optimized and filtered vertically and a seasonally-averaged profile was used. Logarithmic transformation applied to the data prior to the analysis. Background field: the data mean value is subtracted from the data. Detrending of data: no. Advection constraint applied: no. Originators of Italian data sets-List of contributors • Brunetti Fabio (OGS) • Cardin Vanessa, Bensi Manuel doi:10.6092/36728450-4296-4e6a-967d-d5b6da55f306 • Cardin Vanessa, Bensi Manuel, Ursella Laura, Siena Giuseppe doi:10.6092/f8e6d18e-f877-4aa5-a983-a03b06ccb987 • Cataletto Bruno (OGS) • Cinzia Comici Cinzia (OGS) • Civitarese Giuseppe (OGS) • DeVittor Cinzia (OGS) • Giani Michele (OGS) • Kovacevic Vedrana (OGS) • Mosetti Renzo (OGS) • Solidoro C.,Beran A.,Cataletto B.,Celussi M.,Cibic T.,Comici C.,Del Negro P.,De Vittor C.,Minocci M.,Monti M.,Fabbro C.,Falconi C.,Franzo A.,Libralato S.,Lipizer M.,Negussanti J.S.,Russel H.,Valli G., doi:10.6092/e5518899-b914-43b0-8139-023718aa63f5 • Celio Massimo (ARPA FVG) • Malaguti Antonella (ENEA) • Fonda Umani Serena (UNITS) • Bignami Francesco (ISAC/CNR) • Boldrini Alfredo (ISMAR/CNR) • Marini Mauro (ISMAR/CNR) • Miserocchi Stefano (ISMAR/CNR) • Zaccone Renata (IAMC/CNR) • Lavezza, R., Dubroca, L. F. C., Ludicone, D., Kress, N., Herut, B., Civitarese, G., Cruzado, A., Lefèvre, D., Souvermezoglou, E., Yilmaz, A., Tugrul, S., and Ribera d’Alcala, M.: Compilation of quality controlled nutrient profiles from the Mediterranean Sea, doi:10.1594/PANGAEA.771907, 2011. Units: umol/l