Moving 6-year analysis and visualization of Water body silicate in the North Sea. Four seasons (December-February, March-May, June-August, September-November). Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVAnd (Data-Interpolating Variational Analysis) tool, version 2.7.12. results were subjected to the minfield option in DIVAnd to avoid negative/underestimated values in the interpolated results; error threshold masks L1 (0.3) and L2 (0.5) are included as well as the unmasked field. The depth dimension allows visualizing the gridded field at various depths.

'''DEFINITION''' Ocean heat content (OHC) is defined here as the deviation from a reference period (1993-2014) and is closely proportional to the average temperature change from z1 = 0 m to z2 = 700 m depth: OHC=∫_(z_1)^(z_2)ρ_0 c_p (T_yr-T_clim )dz [1] with a reference density of = 1030 kgm-3 and a specific heat capacity of cp = 3980 J kg-1 °C-1 (e.g. von Schuckmann et al., 2009). Time series of annual mean values area averaged ocean heat content is provided for the Mediterranean Sea (30°N, 46°N; 6°W, 36°E) and is evaluated for topography deeper than 300m. '''CONTEXT''' Knowing how much and where heat energy is stored and released in the ocean is essential for understanding the contemporary Earth system state, variability and change, as the oceans shape our perspectives for the future. The quality evaluation of MEDSEA_OMI_OHC_area_averaged_anomalies is based on the “multi-product” approach as introduced in the second issue of the Ocean State Report (von Schuckmann et al., 2018), and following the MyOcean’s experience (Masina et al., 2017). Six global products and a regional (Mediterranean Sea) product have been used to build an ensemble mean, and its associated ensemble spread. The reference products are: • The Mediterranean Sea Reanalysis at 1/24 degree horizontal resolution (MEDSEA_MULTIYEAR_PHY_006_004, DOI: https://doi.org/10.25423/CMCC/MEDSEA_MULTIYEAR_PHY_006_004_E3R1, Escudier et al., 2020) • Four global reanalyses at 1/4 degree horizontal resolution (GLOBAL_MULTIYEAR_PHY_ENS_001_031): GLORYS, C-GLORS, ORAS5, FOAM • Two observation based products: CORA (INSITU_GLO_PHY_TS_OA_MY_013_052) and ARMOR3D (MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012). Details on the products are delivered in the PUM and QUID of this OMI. '''CMEMS KEY FINDINGS''' The ensemble mean ocean heat content anomaly time series over the Mediterranean Sea shows a continuous increase in the period 1993-2022 at rate of 1.38±0.08 W/m2 in the upper 700m. After 2005 the rate has clearly increased with respect the previous decade, in agreement with Iona et al. (2018). '''DOI (product):''' https://doi.org/10.48670/moi-00261

Seasonal climatology of Water body silicate for Loire river for the period 1950-2024 and for the following seasons: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Observation data span from 1950 to 2024. Depth levels (m): [0.0, 2.0, 4.0, 6.0, 8.0, 10.0, 15.0, 20.0, 25.0, 30.0, 35.0, 40.0, 45.0, 50.0, 60.0, 70.0, 80.0, 90.0, 100.0, 110.0, 120.0, 130.0]. Data sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVAnd analysis: the computation was done with DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.12, using GEBCO 15 sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps is 0.01 degrees. Horizontal correlation length is defined seasonally (in meters): 133000 (winter), 180000 (spring), 150000 (summer), 210000 (autumn). Vertical correlation length was optimized and vertically filtered and a seasonally-averaged profile was used (DIVAnd.fitvertlen). Signal-to-noise ratio was fixed to 1 for vertical profiles and 0.1 for time series to account for the redundancy in the time series observations. A logarithmic transformation (DIVAnd.loglin) was applied to the data prior to the analysis. Background field: the vertically-filtered data mean profile is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: umol/l.

As part of the European Horizon Europe FOCCUS project (https://foccus-project.eu/), the metadata inventory of European coastal platforms has been extracted. The inventory was based on the following History and Latest products, downloaded from the CMEMS website (https://marine.copernicus.eu/fr/acces-donnees) at: 1) Global Ocean-In-Situ Near-Real-Time Observation, 2) Atlantic Iberian Biscay Irish Ocean-In-Situ Near Real Time Observations, 3) Mediterranean Sea-In-Situ Near Real Time Observations, 4) Atlantic-European North West Shelf-Ocean In-Situ Near Real Time Observations. To carry out this inventory, it was decided to target only coastal platforms, located less than 200km from the coast and at a depth of less than 400m. For mobile platforms, it was also decided to focus only on the first position in the file. This data must be located within 200 km of the coast and at a depth of less than 400 m. In this inventory, FerryBox platforms have all been considered as coastal platforms. The following platforms were extracted from the products: BO (Bottles), CT (CTD), DB (Drifting Buoys), FB (Ferry Box), GL (Gliders), HF (High Frequency Radar), MO (Mooring), PF (Profiling Float), TG (Tide Gauge) and XB (XBT). Once the metadata had been extracted from the files, duplicates were removed (files with the same names). Duplicate platforms of type _TS_ and _WS_ were merged (date and parameters). Latest‘ files have been merged with ’History" files. Missing metadata have been replaced in the Excel file by ‘Missing Data’. Some old dates were also revised by hand because they had been badly extracted, as well as some institution names that included special characters. Platforms located on estuaries/rivers/lakes/ponds have also been removed by hand. This inventory identified a total of 10,479 coastal platforms.

'''DEFINITION''' Estimates of Ocean Heat Content (OHC) are obtained from integrated differences of the measured temperature and a climatology along a vertical profile in the ocean (von Schuckmann et al., 2018). The products used include three global reanalyses: GLORYS, C-GLORS, ORAS5 (GLOBAL_MULTIYEAR_PHY_ENS_001_031) and two in situ based reprocessed products: CORA5.2 (INSITU_GLO_PHY_TS_OA_MY_013_052) , ARMOR-3D (MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012). Additionally, the time series based on the method of von Schuckmann and Le Traon (2011) has been added. The regional OHC values are then averaged from 60°S-60°N aiming i) to obtain the mean OHC as expressed in Joules per meter square (J/m2) to monitor the large-scale variability and change. ii) to monitor the amount of energy in the form of heat stored in the ocean (i.e. the change of OHC in time), expressed in Watt per square meter (W/m2). Ocean heat content is one of the six Global Climate Indicators recommended by the World Meterological Organisation for Sustainable Development Goal 13 implementation (WMO, 2017). '''CONTEXT''' Knowing how much and where heat energy is stored and released in the ocean is essential for understanding the contemporary Earth system state, variability and change, as the ocean shapes our perspectives for the future (von Schuckmann et al., 2020). Variations in OHC can induce changes in ocean stratification, currents, sea ice and ice shelfs (IPCC, 2019; 2021); they set time scales and dominate Earth system adjustments to climate variability and change (Hansen et al., 2011); they are a key player in ocean-atmosphere interactions and sea level change (WCRP, 2018) and they can impact marine ecosystems and human livelihoods (IPCC, 2019). '''CMEMS KEY FINDINGS''' Since the year 2005, the near-surface (0-300m) near-global (60°S-60°N) ocean warms at a rate of 0.4 ± 0.1 W/m2. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00233

This dataset comprises energy density and proximal composition (water, ash, lipid and protein contents) for anchovy (Engraulis encrasicolus) and sardine (Sardina pilchardus) from the Bay of Biscay, the English Channel and the southern North Sea between 2014 and 2017. Fish were sampled throughout various seasons thanks to the PELGAS (May, Bay of Biscay), EVHOE (October-November, Bay of Biscay), CGFS and CAMANOC (September-October, English Channel) surveys conducted by IFREMER on the RV “Thalassa”, to the JUVENA (September-October, Bay of Biscay) survey conducted by AZTI, and from commercial landings within the European Data Collection Framework (DCF) during the CAPTAIN project (France Filière Pêche). During the surveys, pelagic (PELGAS and JUVENA) or demersal (EVHOE and CGFS) trawl hauls are undertaken to identify species and measure individual fish traits. Professional sampling was performed from pelagic trawl or purse-seine catches. From the various surveys, a sub-sampling of the trawls was performed to cover as much as possible the spatial extent of the surveys along the french coast. From the various selected trawls, a sub-sampling of 5 fish per size class (when possible) was performed to cover the size range of each species, based on the following size classes : sardine (1 : <15 cm ; 2 : 15-20 cm ; 3 : >20cm), anchovy (1 : <10cm ; 2 : 10-14 cm ; 3 : >14cm). Each fish was individually measured to the nearest tenth of a centimeter and weighted to the nearest tenth of a gram. These measurements were taken either at sea or later in the laboratory. The collected fish were frozen individually at -20°C before laboratory processing. In the laboratory, maturity stages were determined following ICES guidelines (ICES, 2008) based on macroscopic gonads observation and using a six-stage key as follows: stages 1 & 2 indicate immature and developing individuals, stages 3–5 indicate three steps of increasing gonad development and the spawning period (stage 3: partial pre-spawning; stage 4: spawning (hydrated); stage 5: partial post-spawning), and stage 6 features the final post-spawning period. Fish characterised by maturity stages 3, 4 or 5 were considered as being in an active reproductive period as opposed to fish in stages 1, 2 or 6. Fish were then ground and freeze-dried during at least 48 hours. Water content of the entire fish was determined from dry mass and wet mass ratio. Then, fish were ground again to obtain fine homogeneous dry powder for subsequent analysis. Energy density measurements were performed following the protocols of Gatti et al. (2018). Two subsamples of fish powder were placed in an adiabatic bomb calorimeter (IKA C-4000 adiabatic bomb calorimeter, IKA-WerkeGmbh & co. KG) for energy measurements. The energy density (ED, kJ.g-1 dry mass) was determined by measuring the heat released through the combustion of a small subsample, approximately 200 mg. If the coefficient of variation between the two measurements exceeded 3%, a third measurement was made. Finally, ED subsamples measurements were averaged and assigned to each individual fish. Energy density analyses were conducted on 503 individuals for anchovy and 976 individuals for sardine. Ash content was determined gravimetrically by combusting dried tissue in a muffle furnace at 550°C for six hours. Lipids and proteins were analysed by a certified laboratory (Labocea, Plouzané, France). Protein content was estimated using the Kjeldahl method. It consists in first determining the quantity of nitrogen contained in the sample, and to convert it into protein content using a conversion factor (6.25 here). Lipid content was determined by hydrolysis, using petroleum ether as an organic solvent. Carbohydrates represent less than 1% of fish mass and were thus neglected. Protein, lipid and ash content did not exactly sum to 1 in DW (anchovy: mean = 0.91, sd = 0.04; sardine: mean = 0.90, sd = 0.04). This discrepancy may arise from residual water, measurement uncertainties, or to a lesser extent the exclusion of carbohydrates. Body component contents have been normalised by dividing each component by the sum of lipids, proteins and ash content, to sum to one, enabling comparisons between fishes, assuming proportional errors across the components. A total of 116 and 104 proximate composition analyses were performed for anchovy and sardine, respectively.

Moving 6-year analysis of Water body dissolved oxygen concentration in the 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 6-year centered average of the season. 6-years periods span from 1970-1975 until 2019-2024. Description of DIVA analysis: The computation was done with the DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.12, using GEBCO 30sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps grids is dx=dy=0.125 degrees (around 13.5km and 10.9km accordingly). The vertical resolution is 27 depth levels: [0.,5.,10.,20.,30.,50.,75.,100.,125.,150.,200.,250.,300.,400.,500.,600.,700.,800.,900.,1000.,1100.,1200.,1300.,1400.,1500.,1750.,2000.]. The horizontal correlation length is 200km. The vertical correlation length (in meters) was set twices the vertical resolution: [10.,10.,20.,20.,40.,50.,50.,50.,50.,100.,100.,100.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.,500.,500.,500.]. Duplicates check was performed using the following criteria for space and time: dlon=0.001deg., dlat=0.001deg., ddepth=1m, dtime=1hour, dvalue=0.1. The error variance (epsilon2) was set equal to 1 for profiles and 10 for time series to reduce the influence of close data near the coasts. An anamorphosis transformation was applied to the data (function DIVAnd.Anam.loglin) to avoid unrealistic negative values: threshold value=200. A background analysis field was used for all years (1970-2024) with correlation length equal to 600km and error variance (epsilon2) equal to 20. Quality control of the observations was applied using the interpolated field (QCMETHOD=3). Residuals (differences between the observations and the analysis (interpolated linearly to the location of the observations) were calculated. Observations with residuals outside the minimum and maximum values of the 99% quantile were discarded from the analysis. 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.

This visualization product displays the spatial distribution of seafloor litter density per trawl. EMODnet Chemistry included the collection of marine litter in its 3rd phase. Since the beginning of 2018, data of seafloor litter collected by international fish-trawl surveys have been gathered and processed in the EMODnet Chemistry Marine Litter Database (MLDB). The harmonization of all the data has been the most challenging task considering the heterogeneity of the data sources, sampling protocols (OSPAR and MEDITS protocols) and reference lists used on a European scale. Moreover, within the same protocol, different gear types are deployed during bottom trawl surveys. In cases where the wingspread and/or number of items were/was unknown, it was not possible to use the data because these fields are needed to calculate the density. Data collected before 2011 are concerned by this filter. When the distance reported in the data was null, it was calculated from: - the ground speed and the haul duration using the following formula: Distance (km) = Haul duration (h) * Ground speed (km/h); - the trawl coordinates if the ground speed and the haul duration were not filled in. The swept area was calculated from the wingspread (which depends on the fishing gear type) and the distance trawled: Swept area (km²) = Distance (km) * Wingspread (km) Densities were calculated on each trawl and year using the following computation: Density (number of items per km²) = ∑Number of items / Swept area (km²) Then a grid with 30km x 30km cells was used to calculate the weighted mean of densities in each cell from the formula : Weighted mean (number of items per km²) = ∑ (Distance (km) * Density (number of items per km²)) / ∑ Distance (km) Percentiles 50, 75, 95 & 99 were calculated taking into account data for all years. More information on data processing and calculation are detailed in the attached methodology document. Warning: the absence of data on the map does not necessarily mean that they do not exist, but that no information has been entered in the Marine Litter Database for this area. This work is based on the work presented in the following scientific article: O. Gerigny, M. Brun, M.C. Fabri, C. Tomasino, M. Le Moigne, A. Jadaud, F. Galgani, Seafloor litter from the continental shelf and canyons in French Mediterranean Water: Distribution, typologies and trends, Marine Pollution Bulletin, Volume 146, 2019, Pages 653-666, ISSN 0025-326X, https://doi.org/10.1016/j.marpolbul.2019.07.030.

Sardine physiological measurments from september to november 2020

Ensemble simulations of the ecosystem model Apecosm (https://apecosm.org) forced by the IPSL-CM6-LR climate model with the climate change scenario SSP1-2.6. The output files contain yearly mean biomass density for 3 communities (epipelagic, mesopelagic migratory and mesopelagic redidents) and 100 size classes (ranging from 0.12cm to 1.96m) The model grid file is also provided. Units are in J/m2 and can be converted in kg/m2 by dividing by 4e6. These outputs are associated with the "Assessing the time of emergence of marine ecosystems from global to local scales using IPSL-CM6A-LR/APECOSM climate-to-fish ensemble simulations" paper from the Earth's Future "Past and Future of Marine Ecosystems" Special Collection.