These gridded products are produced from the following upstream data: - for satellites SARAL/AltiKa, Cryosat-2, HaiYang-2B, Jason-3, Copernicus Sentinel-3A/B, Sentinel-6 MF, SWOT Nadir => NRT (Near-Real-Time) Nadir along-track (or Level-3) SEA LEVEL products (DOI: https://doi.org/10.48670/moi-00147) delivered by the Copernicus Marine Service (http://marine.copernicus.eu/ ). The gridded product is based on near-real-time (NRT) Level-3 Nadir datasets for the period from July 7, 2025, to December 31, 2025. => MY (Multi-Year) Nadir along-track (or Level-3) SEA LEVEL products (DOI: https://doi.org/10.48670/moi-00146 ) delivered by the Copernicus Marine Service (CMEMS, http://marine.copernicus.eu/ ). The gridded product is based on MY Level-3 Nadir datasets for the period from March 28, 2023, to July 6, 2025. - for SWOT KaRIn : the L3_LR_SSH Expert v3.0 product distributed by AVISO (DOI: https://doi.org/10.24400/527896/A01-2023.018) from March 28, 2023 to December 31, 2025. One mapping algorithm is proposed: the MIOST approach which provides which provides global Sea Surface Height (SSH) solutions. The MIOST method is capable of accounting for various modes of ocean surface topography variability (e.g., geostrophic, barotropic, equatorial wave dynamics) by constructing multiple independent components within a predefined covariance model.

These gridded products are produced from the along-track (or Level-3) SEA LEVEL products (DOI: doi.org/10.48670/moi-00147) delivered by the Copernicus Marine Service (CMEMS, marine.copernicus.eu) for satellites SARAL/AltiKa, Cryosat-2, HaiYang-2B, Jason-3, Copernicus Sentinel-3A/B, Sentinel-6 MF, SWOT nadir, and SWOT Level-3 KaRIn sea level products (DOI: https://doi.org/10.24400/527896/A01-2023.018). Three mapping algorithms are proposed: MIOST, 4DvarNET, 4DvarQG: - the MIOST approach which give the global SSH solutions: the MIOST method is able of accounting for various modes of variability of the ocean surface topography (e.g., geostrophic, barotrope, equatorial waves dynamic …) by constructing several independent components within an assumed covariance model. - the 4DvarNET approach for the regional SSH solutions: the 4DvarNET mapping algorithm is a data-driven approach combining a data assimilation scheme associated with a deep learning framework. - the 4DvarQG approach for the regional SSH solutions: the 4DvarQG mapping technique integrates a 4-Dimensional variational (4DVAR) scheme with a Quasi-Geostrophic (QG) model.

The Sentinel-6 Level-2P skewness products was developed to estimate the skewness from Sentinel-6 LR (Low Resolution Mode) and HR (High Resolution Mode) acquisitions. That demonstration product is generated by different retracking processes, provides an initial estimation of such a phenomenon and allows a finer description of the sea state.

'''DEFINITION''' The trend map is derived from version 5 of the global climate-quality chlorophyll time series produced by the ESA Ocean Colour Climate Change Initiative (ESA OC-CCI, Sathyendranath et al. 2019; Jackson 2020) and distributed by CMEMS. The trend detection method is based on the Census-I algorithm as described by Vantrepotte et al. (2009), where the time series is decomposed as a fixed seasonal cycle plus a linear trend component plus a residual component. The linear trend is expressed in % year -1, and its level of significance (p) calculated using a t-test. Only significant trends (p < 0.05) are included. '''CONTEXT''' Phytoplankton are key actors in the carbon cycle and, as such, recognised as an Essential Climate Variable (ECV). Chlorophyll concentration is the most widely used measure of the concentration of phytoplankton present in the ocean. Drivers for chlorophyll variability range from small-scale seasonal cycles to long-term climate oscillations and, most importantly, anthropogenic climate change. Due to such diverse factors, the detection of climate signals requires a long-term time series of consistent, well-calibrated, climate-quality data record. Furthermore, chlorophyll analysis also demands the use of robust statistical temporal decomposition techniques, in order to separate the long-term signal from the seasonal component of the time series. '''CMEMS KEY FINDINGS''' The average global trend for the 1997-2021 period was 0.51% per year, with a maximum value of 25% per year and a minimum value of -6.1% per year. Positive trends are pronounced in the high latitudes of both northern and southern hemispheres. The significant increases in chlorophyll reported in 2016-2017 (Sathyendranath et al., 2018b) for the Atlantic and Pacific oceans at high latitudes appear to be plateauing after the 2021 extension. The negative trends shown in equatorial waters in 2020 appear to be remain consistent in 2021. '''DOI (product):''' https://doi.org/10.48670/moi-00230

'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' This product is a REP L4 global total velocity field at 0m and 15m. It consists of the zonal and meridional velocity at a 3h frequency and at 1/4 degree regular grid. These total velocity fields are obtained by combining CMEMS REP satellite Geostrophic surface currents and modelled Ekman currents at the surface and 15m depth (using ECMWF ERA5 wind stress). 3 hourly product, daily and monthly means are available. This product has been initiated in the frame of CNES/CLS projects. Then it has been consolidated during the Globcurrent project (funded by the ESA User Element Program). '''DOI (product) :''' https://doi.org/10.48670/moi-00050 '''Product Citation:''' Please refer to our Technical FAQ for citing products: http://marine.copernicus.eu/faq/cite-cmems-products-cmems-credit/?idpage=169.

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' In wavenumber spectra, the 1hz measurement error is the noise level estimated as the mean value of energy at high wavenumbers (below 20km in term of wave length). The 1hz noise level spatial distribution follows the instrumental white-noise linked to the Surface Wave Height but also connections with the backscatter coefficient. The full understanding of this hump of spectral energy (Dibarboure et al., 2013, Investigating short wavelength correlated errors on low-resolution mode altimetry, OSTST 2013 presentation) still remain to be achieved and overcome with new retracking, new editing strategy or new technology. '''DOI (product) :''' https://doi.org/10.48670/moi-00144

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' Marine primary production corresponds to the amount of inorganic carbon which is converted into organic matter during the photosynthesis, and which feeds upper trophic layers. The daily primary production is estimated from satellite observations with the Antoine and Morel algorithm (1996). This algorithm modelized the potential growth in function of the light and temperature conditions, and with the chlorophyll concentration as a biomass index. The monthly area average is computed from monthly primary production weighted by the pixels size. The trend is computed from the deseasonalised time series (1998-2019), following the Vantrepotte and Mélin method. More details are provided in the Ocean State Reports 4 (Cossarini et al. ,2020). '''CONTEXT''' Marine primary production is at the basis of the marine food web and produce about 50% of the oxygen we breath every year (Behrenfeld et al., 2001). Study primary production is of paramount importance as ocean health and fisheries are directly linked to the primary production (Pauly and Christensen, 1995, Fee et al., 2019). Changes in primary production can have consequences on biogeochemical cycles, and specially on the carbon cycle, and impact the biological carbon pump intensity, and therefore climate (Chavez et al., 2011). Despite its importance for climate and socio-economics resources, primary production measurements are scarce and do not allow a deep investigation of the primary production evolution over decades. Satellites observations and modelling can fill this gap. However, depending of their parametrisation, models can predict an increase or a decrease in primary production by the end of the century (Laufkötter et al., 2015). Primary production from satellite observations present therefore the advantage to dispose an archive of more than two decades of global data. This archive can be assimilated in models, in addition to direct environmental analysis, to minimise models uncertainties (Gregg and Rousseaux, 2019). In the Ocean State Reports 4, primary production estimate from satellite and from modelling are compared at the scale of the Mediterranean Sea. This demonstrate the ability of such a comparison to deeply investigate physical and biogeochemical processes associated to the primary production evolution (Cossarini et al., 2020) '''CMEMS KEY FINDINGS''' The trend for the global ocean is negative over the period 1998-2019 with a decline in primary production of about 0.67 mgC.m-2.yr-1 or equivalently 0.2 %.yr-1. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00225

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' The Global Ocean Satellite monitoring and marine ecosystem study group (GOS) of the Italian National Research Council (CNR), in Rome operationally produces Level-4 product includes monthly averaged datasets of the diffuse attenuation coefficient of light at 490 nm (Kd490) for multi-sensor (MODIS-AQUA, NOAA20-VIIRS, NPP-VIIRS, Sentinel3A-OLCI at 300m of resolution) (at 1 km resolution) and Sentinel3A-OLCI observations (at 300m resolution). Kd490 is the diffuse attenuation coefficient of light at 490 nm, and is a measure of the turbidity of the water column, i.e., how visible light in the blue-green region of the spectrum penetrates the water column. It is directly related to the presence of absorbing and scattering matter in the water column and is estimated through the ratio between Rrs at 490 and 555 nm. For the multi-sensor dataset, single sensor Rrs fields are band-shifted, over the SeaWiFS native bands (using the QAAv6 model, Lee et al., 2002) and merged with a technique aimed at smoothing the differences among different sensors. This technique is developed by the GOS. The QAA allows the inversion of the radiative transfer equations to compute the Inherent Optical Properties. Level-4 product includes monthly averages along with the standard deviation and the number of observations in the period of integration. '''Processing information:''' Multi-sensor products are constituted by MODIS-AQUA, NOAA20-VIIRS, NPP-VIIRS and Sentinel3A-OLCI. For consistency with NASA L2 dataset, BRDF correction was applied to Sentinel3A-OLCI prior to band shifting and multi sensor merging. Hence, the single sensor OLCI data set is also distributed after BRDF correction. Single sensor NASA Level-2 data are destriped and then all Level-2 data are remapped at 1 km spatial resolution (300m for Sentinel3A-OLCI) using cylindrical equirectangular projection. Afterwards, single sensor Rrs fields are band-shifted, over the SeaWiFS native bands (using the QAAv6 model, Lee et al., 2002) and merged with a technique aimed at smoothing the differences among different sensors. This technique is developed by The Global Ocean Satellite monitoring and marine ecosystem study group (GOS) of the Italian National Research Council (CNR, Rome). Then geophysical fields (i.e. chlorophyll, kd490, bbp, aph and adg) are estimated via state-of-the-art algorithms for better product quality. Time averages are computed on the delayed-time data. '''Description of observation methods/instruments:''' Ocean colour technique exploits the emerging electromagnetic radiation from the sea surface in different wavelengths. The spectral variability of this signal defines the so-called ocean colour which is affected by the presence of phytoplankton. '''Quality / Accuracy / Calibration information:''' A detailed description of the calibration and validation activities performed over this product can be found on the CMEMS web portal. '''Suitability, Expected type of users / uses:''' This product is meant for use for educational purposes and for the managing of the marine safety, marine resources, marine and coastal environment and for climate and seasonal studies. '''Dataset names :''' *dataset-oc-med-opt-multi-l4-kd490_1km_monthly-rt-v02 *dataset-oc-med-opt-olci-l4-kd490_300m_monthly-rt '''Files format:''' *CF-1.4 *INSPIRE compliant '''DOI (product) :''' https://doi.org/10.48670/moi-00117

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' The time series are derived from the regional chlorophyll reprocessed (REP) product as distributed by CMEMS. This dataset, derived from multi-sensor (SeaStar-SeaWiFS, AQUA-MODIS, NOAA20-VIIRS, NPP-VIIRS, Envisat-MERIS and Sentinel3A-OLCI) Rrs spectra produced by CNR using an in-house processing chain, is obtained by means of the Mediterranean Ocean Colour regional algorithms: an updated version of the MedOC4 (Case 1 (off-shore) waters, Volpe et al., 2019, with new coefficients) and AD4 (Case 2 (coastal) waters, Berthon and Zibordi, 2004). The processing chain and the techniques used for algorithms merging are detailed in Colella et al. (2021). Monthly regional mean values are calculated by performing the average of 2D monthly mean (weighted by pixel area) over the region of interest. The deseasonalized time series is obtained by applying the X-11 seasonal adjustment methodology on the original time series as described in Colella et al. (2016), and then the Mann-Kendall test (Mann, 1945; Kendall, 1975) and Sens’s method (Sen, 1968) are subsequently applied to obtain the magnitude of trend. '''CONTEXT''' Phytoplankton and chlorophyll concentration as a proxy for phytoplankton respond rapidly to changes in environmental conditions, such as light, temperature, nutrients and mixing (Colella et al. 2016). The character of the response depends on the nature of the change drivers, and ranges from seasonal cycles to decadal oscillations (Basterretxea et al. 2018). Therefore, it is of critical importance to monitor chlorophyll concentration at multiple temporal and spatial scales, in order to be able to separate potential long-term climate signals from natural variability in the short term. In particular, phytoplankton in the Mediterranean Sea is known to respond to climate variability associated with the North Atlantic Oscillation (NAO) and El Niño Southern Oscillation (ENSO) (Basterretxea et al. 2018, Colella et al. 2016). '''CMEMS KEY FINDINGS''' In the Mediterranean Sea, the trend average for the 1997-2020 period is slightly negative (-0.580.62% per year). Due to the change in processing techniques and chlorophyll retrieval, this trend estimate cannot be compared directly to those previously reported. The observations time series (in grey) shows minima values have been quite constant until 2015 and then there is a little decrease up to 2020, when an absolute minimum occurs with values lower than 0.04 mg m-3. Throughout the time series, maxima are variable year by year (with absolute maximum in 2015, >0.14 mg m-3), showing an evident reduction since 2016. In the last years of the series, the decrease of chlorophyll concentrations is also observed in the deseasonalized timeseries (in green) with a marked step in 2020. This attenuation of chlorophyll values in the last years results in an overall negative trend for the Mediterranean Sea. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00259

'''This product has been archived''' For operational and online products, please visit https://marine.copernicus.eu '''Short description:''' For the Global Ocean - the OSTIA diurnal skin Sea Surface Temperature product provides daily gap-free maps of: *Hourly mean skin Sea Surface Temperature at 0.25° x 0.25° horizontal resolution, using in-situ and satellite data from infra-red radiometers. The Operational Sea Surface Temperature and Ice Analysis (OSTIA) system is run by the Met Office. A 1/4° (approx. 28 km) hourly analysis of skin Sea Surface temperature (SST) is produced daily for the global ocean. The skin temperature of the ocean is the temperature measured by satellite infra-red radiometers and can experience a large diurnal cycle. The skin SST L4 product is created by combining: 1. the OSTIA foundation SST analysis which uses in-situ and satellite observations; 2. the OSTIA diurnal warm layer analysis which uses satellite observations; and 3. a cool skin model. OSTIA uses satellite data provided by the GHRSST project. '''DOI (product) :''' https://doi.org/10.48670/moi-00167