• Energy Research
• Open Access close
• Restricted close
• ES close
• PT close
• IT close
• European Marine Science close

Species distribution models (SDMs) are becoming an important tool for marine conservation and management. Yet while there is an increasing diversity and volume of marine biodiversity data for training SDMs, little practical guidance is available on how to leverage distinct data types to build robust models. We explored the effect of different data types on the fit, performance and predictive ability of SDMs by comparing models trained with four data types for a heavily exploited pelagic fish, the blue shark (Prionace glauca), in the Northwest Atlantic: two fishery-dependent (conventional mark-recapture tags, fisheries observer records) and two fishery-independent (satellite-linked electronic tags, pop-up archival tags). We found that all four data types can result in robust models, but differences among spatial predictions highlighted the need to consider ecological realism in model selection and interpretation regardless of data type. Differences among models were primarily attributed to biases in how each data type, and the associated representation of absences, sampled the environment and summarized the resulting species distributions. Outputs from model ensembles and a model trained on all pooled data both proved effective for combining inferences across data types and provided more ecologically realistic predictions than individual models. Our results provide valuable guidance for practitioners developing SDMs. With increasing access to diverse data sources, future work should further develop truly integrative modeling approaches that can explicitly leverage strengths of individual data types while statistically accounting for limitations, such as sampling biases.  Please see the README document ("README.md") and the accompanying published article: Braun, C. D., M. C. Arostegui, N. Farchadi, M. Alexander, P. Afonso, A. Allyn, S. J. Bograd, S. Brodie, D. P. Crear, E. F. Culhane, T. H. Curtis, E. L. Hazen, A. Kerney, N. Lezama-Ochoa, K. E. Mills, D. Pugh, N. Queiroz, J. D. Scott, G. B. Skomal, D. W. Sims, S. R. Thorrold, H. Welch, R. Young-Morse, R. Lewison. In press. Building use-inspired species distribution models: using multiple data types to examine and improve model performance. Ecological Applications. Accepted. DOI: < article DOI will be added when it is assigned >

[CASTELLÀ] En este trabajo de final de grado hemos realizado un prototipo mecánico capaz de captar la energía que se pierde de las olas del mar (energía undimotriz) en el espigón de Vilanova i la Geltrú. Estudiando y analizando los mecanismos y prototipos de las empresas actuales en el mercado, obtenemos unas ideas para nuestro dispositivo. A parte hemos analizado el comportamiento del oleaje en nuestro punto de estudio, que es en el espigón de la playa del faro (Vilanova i la Geltrú) y teniendo en cuenta el oleaje y las mareas, hemos diseñado un prototipo más idóneo para esta ubicación. Nuestra idea final trataría de poner varios dispositivos a lo largo del espigón para así de esta manera poder captar la máxima energía en esta zona. [ANGLÈS] In this final degree project we have made a mechanical prototype capable of capturing the energy that is lost from the waves of the sea (wave energy) in the breakwater of Vilanova i la Geltrú. By studying and analysing the mechanisms and prototypes of current companies in the market, we obtain some ideas for our device. We have also analysed the behaviour of the waves at our point of study, which is at the pier of the beach of the lighthouse (Vilanova i la Geltrú) and taking into account the waves and tides, we have designed a more suitable prototype for this location. Our final idea would try to put several devices along the jetty in order to capture the maximum energy in this area.

The dramatic decline of European eel (Anguilla anguilla) populations over recent decades has attracted considerable attention and concern. Furthermore, little is known concerning the sensitivity of eel’s early stages to projected future environmental change. Here we investigated, for the first time, the potential combined effects of ocean warming (OW; ∆ + 4oC; 18oC) and acidification (OA; ∆ - 0.4 pH units) on the survival and migratory behaviour of A. anguilla glass eels, namely their preference towards riverine cues (freshwater and geosmin). Recently arrived individuals were exposed to isolated and combined OW- and OA–conditions for 100 days, adjusting for the salinity gradients associated with upstream migration. A two-choice test was used to investigate migratory activity and shifts in preference towards freshwater environments. While OW decreased survival and increased migratory activity, OA appears to hinder migratory response, reducing the preference for riverine cues. Our results suggest that future conditions could potentially favor an early settlement of glass eels, reducing the proportion of fully-migratory individuals. Further research into the effects of climate change on eel migration and habitat selection is needed to implement efficient conservation plans for this critically endangered species. Santos&Borges_etal_DataIn this file, the data from our study on the influence of ocean warming and ocean acidification on the riverward migration of glass eels is present. These include survival and behavioural outputs as well as a code sheet for easier table read. The tables are R software ready and can be exported as '.csv' in order to run in the Rscript.BorgesSantos_etal_Data.xlsxSantos&Borges_etal_RscriptThis is the code for the Rscript of the statistical analysis performed in our study.BorgesSantos_etal_Rscript.txt

Dataset abstract This dataset contains quality-checked meteorological observations of air temperature, relative humidity, dew point, barometric pressure and observations of downwelling solar radiation and ultraviolet radiation. Further it contains the wind speed and direction relative to the ship but not corrected for air-flow distortion, and translated into the earth reference frame. For each of these variables observations are available from a portside and starboard side sensor. The dataset also contains, cloud base height and sky cover at three levels measured with a Ceilometer. As additional information the solar azimuth and altitude angle have been calculated for the ship’s position every five minutes and have been added as a one-minute time series using the nearest value. The ship’s position, heading, course and speed over ground are also provided. The wind speed measurements were made at a height of approximately 30.5 meters above sea level. The measurement height of the temperature and humidity probes is 23.7 meters above sea level. The barometric pressure was measured at 20 meters above sea level. The observations have been screened for implausible values and on some occasions despiking based on visual inspection and a rolling interquartile range filter have been applied. Solar radiation measurements are affected by shadowing of the ship, and the air temperature and humidity by the heating of air that passes over the ship. Masks are provided to flag affected observations. The wind speed readings are affected by airflow distortion and should be used with consideration until a dataset of corrected wind speeds is published. More details on airflow distortion can be requested from the contact person. Dataset contents ACE_filtered_meteorological_data_1min.csv, data file, comma-separated values diff_TA1_TA3_WDR2_5min_1.png, metadata, portable network graphics ratio_SR1_SR3_solangle_5min_1.png, metadata, portable network graphics data_file_header.txt, metadata, text README.txt, metadata, text ace_filtered_meteorological_data_change_log.txt, metadata, text Change log v1.1 - The range check for skycover (SC) and cloudlevel (CL) was added to the quality-checking routines. 53 data points violated the range check for these variables: these have now been marked as NaN. v1.0 - Initial release of verified meteorological data. Dataset license This meteorological dataset is made available under the Creative Commons Attribution 4.0 International License (CC BY 4.0) whose full description can be found at https://creativecommons.org/licenses/by/4.0/ {"references": ["Smith, Shawn R., Mark A. Bourassa, and Ryan J. Sharp. Establishing More Truth in True Winds. JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY 16 (1999): 14.", "Stull, Roland. Wet-Bulb Temperature from Relative Humidity and Air Temperature. 2011. Journal of Applied Meteorology and Climatology 50, no. 11 (9 September 2011): 2267\u201369. https://doi.org/10.1175/JAMC-D-11-0143.1."]} The Antarctic Circumnavigation Expedition was made possible by funding from the Swiss Polar Institute and Ferring Pharmaceuticals. Antarctic Circumnavigation Expedition – Delivering Added value To Antarctica (ACE-DATA) is funded by the Swiss Data Science Center as Project number 17-02.

# Long-term changes in taxonomic and functional composition of European marine fish communities The GitHub linked repository is here: [European_demersal_fish_assemblages (](https://github.com/auroreRECE/European_demersal_fish_assemblages)DOI [10.5281/zenodo.11190119](https://zenodo.org/doi/10.5281/zenodo.11190119)) ## Overview This project is dedicated to studying the influence of environmental conditions and fishing on the functional and taxonomic structure of a demersal fish community in Europe. This GitHub repository provides the code of the Receveur et al. (2024) publication in Ecography. ## Data files description ### df\_MFA.csv This file contains the coordinates resulting from the Multiple Factor Analysis (MFA): * X : the row numbers ; * ID_unique : a unique ID number corresponding to the trawls ; * Dim.1 : the coordinate of each trawl on the first MFA dimension ; * Dim.2 : the coordinate of each trawl on the second MFA dimension ; * Dim.3 : the coordinate of each trawl on the third MFA dimension ; ### df\_PCA.csv This file contains the coordinates * X : the row numbers ; * ID_unique : a unique ID number corresponding to the trawls ; * Dim.1 : the coordinate of each trawl on the first PCA dimension ; * Dim.2 : the coordinate of each trawl on the second PCA dimension ; * Dim.3 : the coordinate of each trawl on the third PCA dimension ; ### df\_env.csv This file contains the following environmental parameters: * X : the row numbers ; * ID_unique : a unique ID number corresponding to the trawls ; * Year : the Year of each trawl ; * Quarter : the Quarter of each trawl ; * Ecoregion : the Ecoregion where each trawl has been done; * Survey : the name of the Survey ; * x_my_spatial_id : the longitude of the ICES rectangle where the trawl has been done ; * y_my_spatial_id : the latitude of the ICES rectangle where the trawl has been done ; * my_spatial_id : an ID for the ICES rectangle where the trawl has been done ; * depth : the bottom depth (meters) ; * depth_span : the bottom depth variability (maximum depth of the ICES cell - minimum depth) (meters) ; * chloro_mea: the mean chlorophyll-a concentration (mg/m³) ; * mlotst_mea : the mean mixed layer depth (meters) ; * oxy_bottom_mea : the mean bottom dissolved oxygen (umol/l) ; * oxy_surf_mea : the mean surface dissolved oxygen (umol/l) ; * temp_bottom_mea : the mean bottom temperature (°C) ; * temp_surf_mea : the mean surface temperature (°C) ; * curr_surf_mea : the mean surface current strength (m/s) ; * curr_bottom_mea : the mean bottom current strength (m/s) ; * sal_surf_mea : the mean surface salinity (PSU) ; * chloro_std : the standard deviation of chlorophyll-a concentration (mg/m³) ; * mlotst_std : the standard deviation of mixed layer depth (meters) ; * oxy_bottom_std : the standard deviation of bottom dissolved oxygen (umol/l) ; * oxy_surf_std : the standard deviation of surface dissolved oxygen (umol/l) ; * temp_bottom_std : the standard deviation of bottom temperature (°C) ; * temp_surf_std : the standard deviation of surface temperature (°C) ; * curr_surf_std : the standard deviation of surface current strength (m/s) ; * curr_bottom_std : the standard deviation of bottom current strength (m/s) ; * sal_surf_std : the standard deviation of surface salinity (PSU). ## Raw Data sources ### Biological data Trawls content is publicly available for the North East Atlantic (DATRAS database). Mediterranean data (MEDITS database) are available upon request to Maritime Affairs and Fisheries (MARE DATACOLLECTIONFRAMEWORK). The project uses the following surveys: | Survey Code | Survey name | Area | Period | References | | :---------- | :----------------------------------------------------- | :------------------------------------- | :-------: | :--------: | | BITS | Baltic International Trawl Survey | Baltic Sea | 1994-2019 | 4 | | BTS | Beam Trawl Survey | Celtic Sea; English Channel; North Sea | 1997-2019 | 7 | | BTS-VIII | Beam Trawl Survey – Bay of Biscay | Bay of Biscay | 2011-2019 | 7 | | DWS | Deepwater Survey | Irish Sea | 2006-2007 | 8 | | DYFS | Inshore Beam Trawl Survey | Southern North Sea | 2002-2019 | 7 | | EVHOE | French Southern Atlantic Bottom trawl Survey | Bay of Biscay and Celtic Sea | 2003-2019 | 1 | | FR-CGFS | French Channel ground Survey | English Channel | 1997-2019 | 2 | | IE-IAMS | Irish Anglerfish and megrim Survey | Scottish rockall and Irish Sea | 2016-2019 | 2 | | IE-IGFS | Irish Groundfish | Ireland Shelf Sea | 2003-2019 | 2 | | MEDITS | International bottom trawl survey in the Mediterranean | Mediterranean Sea | 1994-2018 | 9 | | NIGFS | Northern Ireland Groundfish Survey | Irish Sea | 2009-2019 | 2 | | NS-IBTS | North Sea International Bottom Trawl Survey | North Sea | 1997-2019 | 2 | | PT-IBTS | Portuguese International Bottom Trawl Survey | Portugal Shelf Sea | 2003-2017 | 2 | | ROCKALL | Scottish Rockall Survey (until 2010) | Rockall plateau | 2003-2009 | 2 | | SCOROC | Scottish Rockall Survey (from 2011) | Scottish plateau | 2011-2019 | 2 | | SCOWCGFS | Scottish West Coast Groundfish Survey | Scottish west coast | 2011-2019 | 2 | | SNS | Sole Net Survey | Southern North Sea | 2002-2019 | 7 | | SP-ARSA | Spanish Gulf of Cadiz Bottom Trawl Survey | Spain | 2003-2019 | 6 | | SP-NORTH | Spanish North Bottom Trawl Survey | North of Spain | 2003-2019 | 2 | | SP-PORC | Spanish Porcupine Bottom Trawl Survey | Irish Sea | 2003-2019 | 5 | | SWC-IBTS | Scottish West Coast International Bottom Trawl Survey | Scotland Shelf Sea | 1999-2010 | 2 | ### Trait data The complete traits data table is available upon request. It is a combination of the publicly available PANGAEA database, Fishbase information, and inference based on the FISHLIFE project. ### Environmental variables The data used are all publicly available on the Copernicus website. ### Fishing data The data used are all publicly available on the Global Fishing Watch website. ## Recommended Citation Please use the following citation: Receveur, A., Leprieur F., Ellingsen K., Keith D., Kleisner K., McLean M., Mérigot B., Mills K., Mouillot D., Rufino M., Trindade-Santos I., Van Hoey G., Albouy C., Auber A. Data for “Long-term changes in taxonomic and functional composition of European marine fish communities.” Dryad Digital Repository. (2024). doi.org/10.5061/dryad.x69p8czsj ## Acknowledgments This research is a product of the MAESTRO group funded by the synthesis center CESAB of the French Foundation for Research on Biodiversity (FRB). We thank France Filière Pêche (FFP) who founded the MAESTRO project. We also warmly thank all those who have contributed in any way to the scientific surveys and data collection/provision (European Institutions and scientists implicated in DATRAS-BTS, MEDITS, and DCF). ## References 1. ICES. The EVHOE survey (France). ICES Documents. (1997). Available at: https://archimer.ifremer.fr/doc/00036/14707/12013.pdf 2. ICES. Manual of the IBTS North Eastern Atlantic Surveys. Series of ICES Survey Protocols SISP 15 (2017). doi:10.17895/ices.pub.3519 3. ICES. Manual for the International Bottom Trawl Surveys Revision VIII. Series of ICES Survey Protocols SISP 10 - IBTS IX. (2015). 4. https://ices-library.figshare.com/articles/report/SISP_7_-*Manual_for_the_Baltic_International_Trawl_Surveys_BITS*/19050986 5. https://gis.ices.dk/geonetwork/srv/api/records/ce94a257-c8b3-44f7-9fd0-6bd7449ce073 6. http://ices.dk/sites/pub/CM%20Doccuments/2002/D/D0302A.pdf 7. https://ices-library.figshare.com/articles/report/SISP_14_-*Manual_for_the_Offshore_Beam_Trawl_Surveys_WGBEAM*/19051328 8. https://gis.ices.dk/geonetwork/srv/api/records/936b4fb7-9baa-4dbc-abd0-b1b7bda16406 9. https://archimer.ifremer.fr/doc/00117/22783/20585.pdf Evidence of large-scale biodiversity degradation in marine ecosystems has been reported worldwide, yet most research has focused on few species of interest or on limited spatiotemporal scales. Here we assessed the spatial and temporal changes in the taxonomic and functional composition of fish communities in European seas over the last 25 years (1994-2019). We then explored how these community changes were linked to environmental gradients and fishing pressure. We show that the spatial variation in fish species composition is more than two times higher than the temporal variation, with a marked spatial continuum in taxonomic composition and a more homogenous pattern in functional composition. The regions warming the fastest are experiencing an increasing dominance and total abundance of r-strategy fish species (lower age of maturity). Conversely, regions warming more slowly show an increasing dominance and total abundance of K-strategy species (high trophic level and late reproduction). Among the considered environmental variables, sea surface temperature, surface salinity, and chlorophyll-a most consistently influenced communities’ spatial patterns, while bottom temperature and oxygen had the most consistent influence on temporal patterns. Changes in communities’ functional composition were more closely related to environmental conditions than taxonomic changes. Our study demonstrates the importance of integrating community-level species traits across multi-decadal scales and across a large region to better capture and understand ecosystem-wide responses and provides a different lens on community dynamics that could be used to support sustainable fisheries management.

[Experimental design: thermal performance experiments] All experiments were run in climate-controlled incubation facilities of the Institut Mediterrani d’Estudis Avançats (Mallorca, Spain). Following 48 hrs under ambient (collection site) conditions, samples were transferred to individual experimental aquaria, which consisted of a double layered transparent plastic bag filled with 2 L of filtered seawater (60 μm) (following Savva et al. 2018). 16 experimental bags were suspended within 80L temperature-controlled baths. In total, ten baths were used, one for each experimental temperature treatment. Bath temperatures were initially set to the acclimatization temperature (i.e. in situ temperatures) and were subsequently increased or decreased by 1 °C every 24 hours until the desired experimental temperature was achieved. Experimental temperatures were: 15, 18, 21, 24, 26, 28, 30, 32, 34 and 36°C (Table S2). For each species, four replicate aquarium bags were used for each temperature treatment with three individually marked seagrass shoots or three algal fragments placed into each bag. For P. oceanica, each marked plant was a single shoot including leaves, vertical rhizome and roots. For C. nodosa, each marked individual consisted of a 10 cm fragment of horizontal rhizome containing three vertical shoots. Individually marked seaweeds contained the holdfast, and 4-5 fronds of P. pavonica (0.98 ± 0.06 g FW; mean ± SE) or a standardised 5-8 cm fragment with meristematic tip for C. compressa (3.67 ± 0.1 g FW; mean ± SE). Experimental plants were cleaned of conspicuous epiphytes. Once the targeted temperatures were reached in all of the baths, experiments ran for 14 days for the algal species and 21 days for seagrasses to allow for measurable growth in all species at the end of the experiment. Experiments were conducted inside a temperature-controlled chamber at constant humidity and air temperature (15 °C). Bags were arranged in a 4x4 grid within each bath, enabling four species/population treatments to be run simultaneously. Bags were mixed within each bath so that one replicate bag was in each row and column of the grid, to minimise any potential within bath effects of bag position. Replicate bags were suspended with their surface kept open to allow gas exchange and were illuminated with a 14h light:10h dark photoperiod through fluorescent aquarium growth lamps. The water within the bags were mixed with aquaria pumps. The light intensity within each bag was measured via a photometric bulb sensor (LI-COR) and ranged between 180-258 μmol m-2 s-1. Light intensity was constant between experiments and did not significantly differ between experimental treatments (p > 0.05). The temperature in the baths was controlled and recorded with an IKS-AQUASTAR system, which was connected to heaters and thermometers. The seawater within the bags was renewed every 72 hrs and salinity was monitored daily with an YSI multi-parameter meter. Distilled water was added when necessary to ensure salinity levels remained within the range of 36-39 PSU, typical of the study region. Carbon and Nitrogen concentrations in the leaf tissue were measured at the end of the experiment for triplicates of the 24ºC treatment for each species and location (Fig. S2) at Unidade de Técnicas Instrumentais de Análise (University of Coruña, Spain) with an elemental analyser FlashEA112 (ThermoFinnigan). [Species description and distribution] The species used in this study are all common species throughout the Mediterranean Sea, although differ in their biological traits, evolutionary histories and thermo-geographic affinities (Fig. S1). P. oceanica is endemic to the Mediterranean Sea with the all other Posidonia species found in temperate Australia (Aires et al. 2011). The distribution of P. oceanica is restricted to the Mediterranean, spanning from Gibraltar in the west to Cyprus in the east and north into the Aegean and Adriatic seas (Telesca et al. 2015) (Fig. S1A). C. nodosa distribution extends across the Mediterranean Sea and eastern Atlantic Ocean, where it is found from south west Portugal, down the African coast to Mauritania and west to Macaronesia (Alberto et al. 2008) (Fig. S1B). Congeneric species of C. nodosa are found in tropical waters of the Red Sea and Indo-Pacific, suggesting origins in the region at least prior to the closure of the Suez Isthmus, approximately 10Mya. Like C. nodosa, Cystoseira compressa has a distribution that extends across the Mediterranean and into the eastern Atlantic, where it is found west to Macaronesia and south to northwest Africa (Fig. S1C). The genus Cystoseira has recently been reclassified to include just four species with all congeneric Cystoseira spp. having warm-temperate distributions from the Mediterranean to the eastern Atlantic (Orellana et al. 2019). The distribution of Padina pavonica is conservatively considered to resemble C. nodosa and C. compressa, spanning throughout the Mediterranean and into the eastern Atlantic. We considered the poleward distribution limit of P. pavonica to be the British Isles 50ºN (Herbert et al. 2016). P. pavonica was previously thought to have a global distribution, but molecular analysis of the genus has found no evidence to support this (Silberfeld et al. 2013). Instead it has been suggested that P. pavonica was potentially misclassified outside of the Mediterranean, due to morphological similarity with congeneric species (Silberfeld et al. 2013). Padina is a monophyletic genus with a worldwide distribution from tropical to cold temperate waters (Silberfeld et al. 2013). Most species have a regional distribution, with few confirmed examples of species spanning beyond a single marine realm (sensu Spalding et al. 2007). [Metabolic rates] Net production (NP), gross primary production (GPP) and respiration (R) were measured for all species from the four sites for five different experimental temperatures containing the in-situ temperature during sampling up to a 6ºC warming (see SM Table S3 for details). Individuals of the different species were moved to methacrylate cylinders containing seawater treated with UV radiation to remove bacteria and phytoplankton, in incubation tanks at the 5 selected temperatures. Cylinders were closed using gas-tight lids that prevent gas exchange with the atmosphere, containing an optical dissolved oxygen sensor (ODOS® IKS), with a measuring range from 0-200 % saturation and accuracy at 25ºC of 1% saturation, and magnetic stirrers inserted to ensure mixing along the height of the core. Triplicates were measured for each species and location, along with controls consisting in cylinders filled with the UV-treated seawater, in order to account for any residual production or respiration derived from microorganisms (changes in oxygen in controls was subtracted from treatments). Oxygen was measured continuously and recorded every 15 minutes for 24 hours. Changes in the dissolved oxygen (DO) were assumed to result from the biological metabolic processes and represent NP. During the night, changes in DO are assumed to be driven by R, as in the absence of light, no photosynthetic production can occur. R was calculated from the rate of change in oxygen at night, from half an hour after lights went off to half an hour before light went on (NP in darkness equalled R). NP was calculated from the rate of change in DO, at 15 min intervals, accumulated over each 24 h period. Assuming that daytime R equals that during the night, GPP was estimated as the sum of NP and R. To derive daily metabolic rates, we accumulated individual estimates of GPP, NP, and R resolved at 15 min intervals over each 24 h period during experiments and reported them in mmol O2 m−3 day−1. A detailed description of calculation of metabolic rates can be found at Vaquer-Sunyer et al. (Vaquer-Sunyer et al. 2015). [Thermal distribution and thermal safety margins] We estimated the realised thermal distribution for the four experimental species by downloading occurrence records from the Global Biodiversity Information Facility (GBIF.org (11/03/2020) GBIF Occurrence Download). Occurrence records from GBIF were screened for outliers and distributions were verified from the primary literature (Alberto et al. 2008, Draisma et al. 2010, Ni-Ni-Win et al. 2010, Silberfeld et al. 2013, Telesca et al. 2015, Orellana et al. 2019) and Enrique Ballesteros (pers. comms) (Fig. S1). Mean, 1st and 99th percentiles of daily SST’s were downloaded for each occurrence site for the period between 1981-2019 using the SST products described above (Table S4). Thermal range position of species at each experimental site were standardised by their global distribution using a Range Index (RI; Sagarin & Gaines 2002). Median SST at the experimental collection sites were standardized relative to the thermal range observed across a species realized distribution, using the equation: RI = 2(SM- DM)/DB where SM = the median temperature at the experimental collection site, Dm = the thermal midpoint of the species global thermal distribution and DB = range of median temperatures (ºC) that a species experiences across its distribution. The RI scales from -1 to 1, whereby ‘-1’ represents the cool, leading edge of a species distribution, ‘0’ represents the thermal midpoint of a species distribution and ‘1’ represents the warm, trailing edge of a species distribution (Sagarin & Gaines 2002). Thermal safety margins for each population were calculated as the difference between empirically derived upper thermal limits for each population and the maximum long term habitat temperatures recorded at collection sites. Each population’s thermal safety margin was plotted against its range position to examine patterns in thermal sensitivity across a species distribution. [Growth measurements and statistical analyses] Net growth rate of seagrass shoots was measured using leaf piercing-technique (Short & Duarte 2001). At the beginning of the experiment seagrass shoots were pierced just below the ligule with a syringe needle and shoot growth rate was estimated as the elongation of leaf tissue in between the ligule and the mark position of all leaves in a shoot at the end of the experiment, divided by the experimental duration. Net growth rate of macroalgae individuals was measured as the difference in wet weight at the end of the experiment from the beginning of the experiment divided by the duration of the experiment. Moisture on macroalgae specimens was carefully removed before weighing them. Patterns of growth in response to temperature were examined for each experimental population using a gaussian function: g = ke[-0.5(TMA-μ)2/σ2], where k = amplitude, μ = mean and σ = standard deviation of the curve. Best fit values for each parameter were determined using a non-linear least squares regression using the ‘nlstools’ package (Baty et al. 2015) in R (Team 2020). 95% CI for each of the parameters were calculated using non-parametric bootstrapping of the mean centred residuals. The relationship between growth metrics and the best-fit model was determined by comparing the sum of squared deviations (SS) of the observed data from the model, to the SS of 104 randomly resampled datasets. Growth metrics were considered to display a significant relationship to the best-fit model if the observed SS was smaller than the 5th percentile of randomised SS. Upper thermal limits were defined as the optimal temperature + 2 standard deviations (95th percentile of curve) or where net growth = 0. Samples that had lost all pigment or structural integrity by the end of the experiment were considered dead and any positive growth was treated as zero. Comparative patterns in thermal performance between populations have fundamental implications for a species thermal sensitivity to warming and extreme events. Despite this, within-species variation in thermal performance is seldom measured. Here we compare thermal performance between-species variation within communities, for two species of seagrass (Posidonia oceanica and Cymodocea nodosa) and two species of seaweed (Padina pavonica and Cystoseira compressa). Experimental populations from four locations spanning approximately 75% of each species global distribution and a 6ºC gradient in summer temperatures were exposed to 10 temperature treatments (15ºC to 36ºC), reflecting median, maximum and future temperatures. Experimental thermal performance displayed the greatest variability between species, with optimal temperatures differing by over 10ºC within the same location. Within-species differences in thermal performance were also important for P. oceanica which displayed large thermal safety margins within cool and warm-edge populations and small safety margins within central populations. Our findings suggest patterns of thermal performance in Mediterranean seagrasses and seaweeds retain deep ‘pre-Mediterranean’ evolutionary legacies, suggesting marked differences in sensitivity to warming within and between benthic marine communities. [Sample collection] Sample collections were conducted at two sites, separated by approximately 1 km, within each location. Collections were conducted at the same depth (1-3 m) at each location and were spaced across the reef or meadow to try and minimise relatedness between shoots or fragments. Upon collection, fragments were placed into a mesh bag and transported back to holding tanks in cool, damp, dark conditions (following Bennett et al. 2021). Fragments were kept in aerated holding tanks in the collection sites at ambient seawater temperature and maintained under a 14:10 light-dark cycle until transport back to Mallorca, where experiments were performed. Prior to transport, P. oceanica shoots were clipped to 25 cm length (from meristem to tip), to standardise initial conditions and remove old tissue for transport. For transport back to Mallorca, fragments were packed in layers within cool-boxes. Cool-packs were wrapped in damp tea towels (rinsed in seawater) and placed between layers of samples. Samples from Catalonia, Crete and Cyprus experienced approximately 12hrs of transit time. On arrival at the destination, samples were returned to holding tanks with aerated seawater and a 14:10 light-dark cycle. [Sea temperature measurements and reconstruction] Sea surface temperature data for each collection site were based on daily SST maps with a spatial resolution of 1/4°, obtained from the National Center for Environmental Information (NCEI, https://www.ncdc.noaa.gov/oisst (Reynolds et al. 2007). These maps have been generated through the optimal interpolation of Advanced Very High Resolution Radiometer (AVHRR) data for the period 1981-2019. Underwater temperature loggers (ONSET Hobo pro v2 Data logger) were deployed at each site and recorded hourly temperatures throughout one year. In order to obtain an extended time series of temperature at each collection site, a calibration procedure was performed comparing logger data with sea surface temperature from the nearest point on SST maps. In particular, SST data were linearly fitted to logger data for the common period. Then, the calibration coefficients were applied to the whole SST time series to obtain corrected-SST data and reconstruct daily habitat temperatures from 1981-2019. [Field collections] Thermal tolerance experiments were conducted on two seagrass species (P. oceanica and Cymodocea nodosa) and two brown seaweed species (Cystoseira compressa and P. pavonica) from four locations spanning 8 degrees in latitude and 30 degrees in longitude across the Mediterranean (Fig. 1, Table S1). These four species were chosen as they are dominant foundation species and cosmopolitan across the Mediterranean Sea. Thermal performance experiments from Catalonia and Mallorca were conducted simultaneously in June 2016 for seaweeds (P. pavonica and C. compressa) and in August 2016 for seagrasses (P. oceanica and C. nodosa). Experiments for all four species were conducted in July 2017 for Crete and in September 2017 for Cyprus. Horizon 2020 Framework Programme, Award: 659246; Juan de la Cierva Formacion, Award: FJCI-2016-30728; Spanish Ministry of Economy, Industry and Competitiveness, Award: MedShift, CGL2015-71809-P; Spanish Ministry of Science, Innovation and Universities, Award: SUMAECO, RTI2018-095441-B-C21

Project: GCOS Earth Heat Inventory - A study under the Global Climate Observing System (GCOS) concerted international effort to update the Earth heat inventory (EHI), and presents an updated international assessment of ocean warming estimates, and new and updated estimates of heat gain in the atmosphere, cryosphere and land over the period from 1960 to present. Summary: The file “GCOS_EHI_1960-2020_Earth_Heat_Inventory_Ocean_Heat_Content_data.nc” contains a consistent long-term Earth system heat inventory over the period 1960-2020. Human-induced atmospheric composition changes cause a radiative imbalance at the top-of-atmosphere which is driving global warming. Understanding the heat gain of the Earth system from this accumulated heat – and particularly how much and where the heat is distributed in the Earth system - is fundamental to understanding how this affects warming oceans, atmosphere and land, rising temperatures and sea level, and loss of grounded and floating ice, which are fundamental concerns for society. This dataset is based on a study under the Global Climate Observing System (GCOS) concerted international effort to update the Earth heat inventory published in von Schuckmann et al. (2020), and presents an updated international assessment of ocean warming estimates, and new and updated estimates of heat gain in the atmosphere, cryosphere and land over the period 1960-2020. The dataset also contains estimates for global ocean heat content over 1960-2020 for different depth layers, i.e., 0-300m, 0-700m, 700-2000m, 0-2000m, 2000-bottom, which are described in von Schuckmann et al. (2022). This version includes an update of heat storage of global ocean heat content, where one additional product (Li et al., 2022) had been included to the initial estimate. The Earth heat inventory had been updated accordingly, considering also the update for continental heat content (Cuesta-Valero et al., 2023).

Ocean acidification and warming (OA-W) result mainly from the absorption of carbon dioxide and heat by the oceans, altering its physical and chemical properties and affecting carbonate secretion by marine calcifiers such as gastropods. These processes are ongoing, and the projections of their aggravation are not encouraging. This work assesses the concomitant effect of the predicted pH decrease and temperature rise on early life stages of the neogastropod Tritia reticulata (L.), a common scavenger of high ecological importance on coastal ecosystems of the NE Atlantic. Veligers were exposed for 14 days to 12 OA-W experimental scenarios generated by a factorial design of three pH levels (targeting 8.1, 7.8 and 7.5) at four temperatures (16, 18, 20 and 22 °C). Results reveal effects of both pH and temperature (T °C) on larval development, growth, shell integrity and survival, individually or interactively at different exposure times. All endpoints were initially driven by pH, with impaired development and high mortalities being recorded in the first week, constrained by the most acidic scenarios (pHtarget 7.5). Development was also significantly driven by T °C, and its acceleration with warming was observed for the remaining exposure time. Still, by the end of this 2-weeks trial, larval performance and survival were highly affected by the interaction between pH and T °C: growth under warming was evident but only for T °C ≤ 20 °C and carbonate saturation (pHtarget ≥ 7.8). In fact, carbonate undersaturation rendered critical larval mortality (100%) at 22 °C, and the occurrence of extremely vulnerable, unshelled specimens in all other tested temperatures. As recruitment cohorts are the foundation for future populations, our results point towards the extreme vulnerability of this species in case tested scenarios become effective that, according to the IPCC, are projected for the northern hemisphere, where this species is ubiquitous, by the end of the century. Increased veliger mortality associated with reduced growth rates, shell dissolution and loss under OA-W projected scenarios will reduce larval performance, jeopardizing T. reticulata subsistence. In order to allow full comparability with other ocean acidification data sets, the R package seacarb (Gattuso et al, 2019) was used to compute a complete and consistent set of carbonate system variables, as described by Nisumaa et al. (2010). In this dataset the original values were archived in addition with the recalculated parameters (see related PI). The date of carbonate chemistry calculation by seacarb is 2020-06-12.