'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' The indicator of the Kuroshio extension phase variations is based on the standardized high frequency altimeter Eddy Kinetic Energy (EKE) averaged in the area 142-149°E and 32-37°N and computed from the DUACS (https://duacs.cls.fr) delayed-time (reprocessed version DT-2021, CMEMS SEALEVEL_GLO_PHY_L4_MY_008_047) and near real-time (CMEMS SEALEVEL_GLO_PHY_L4_NRT_OBSERVATIONS_008_046) altimeter sea level gridded products. The change in the reprocessed version (previously DT-2018) and the extension of the mean value of the EKE (now 27 years, previously 20 years) induce some slight changes not impacting the general variability of the Kuroshio extension (correlation coefficient of 0.988 for the total period, 0.994 for the delayed time period only). '''CONTEXT''' The long-term mean and trends alone do not give a complete view of the likely changes in position of unstable western boundary current extensions (Kelly et al., 2010). The Kuroshio Extension is an eastward-flowing current in the subtropical western North Pacific after the Kuroshio separates from the coast of Japan at 35°N, 140°E. Being the extension of a wind-driven western boundary current, the Kuroshio Extension is characterized by a strong variability and is rich in large-amplitude meanders and energetic eddies (Niiler et al., 2003; Qiu, 2003, 2002). The Kuroshio Extension region has the largest sea surface height variability on sub-annual and decadal time scales in the extratropical North Pacific Ocean (Jayne et al., 2009; Qiu and Chen, 2010, 2005). Prediction and monitoring of the path of the Kuroshio are of huge importance for local economies as the position of the Kuroshio extension strongly determines the regions where phytoplankton and hence fish are located. '''CMEMS KEY FINDINGS''' The different states of the Kuroshio extension phase have been presented and validated by Bessières et al. (2013) and further reported by Drévillon et al. (2018) in the Copernicus Ocean State Report #2. Two rather different states of the Kuroshio extension are observed: an ‘elongated state’ (also called ‘strong state’) corresponding to a narrow strong steady jet, and a ‘contracted state’ (also called ‘weak state’) in which the jet is weaker and more unsteady, spreading on a wider latitudinal band. When the Kuroshio Extension jet is in a contracted (elongated) state, the upstream Kuroshio Extension path tends to become more (less) variable and regional eddy kinetic energy level tends to be higher (lower). In between these two opposite phases, the Kuroshio extension jet has many intermediate states of transition and presents either progressively weakening or strengthening trends. In 2018, the indicator reveals an elongated state followed by a weakening neutral phase since then. '''DOI (product):''' https://doi.org/10.48670/moi-00222

'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' Experimental altimeter satellite along-track sea surface heights anomalies (SLA) computed with respect to a twenty-year [1993, 2012] mean with a 5Hz (~1.3km) sampling. All the missions are homogenized with respect to a reference mission (see QUID document or http://duacs.cls.fr [http://duacs.cls.fr] pages for processing details). The product gives additional variables (e.g. Mean Dynamic Topography, Dynamic Atmosphic Correction, Ocean Tides, Long Wavelength Errors, Internal tide, …) that can be used to change the physical content for specific needs This product was generated as experimental products in a CNES R&D context. It was processed by the DUACS multimission altimeter data processing system. '''DOI (product) :''' https://doi.org/10.48670/moi-00137

'''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, distributes Level-4 product including the daily interpolated chlorophyll field with no data voids starting from the multi-sensor (MODIS-Aqua, NOAA-20-VIIRS, NPP-VIIRS, Sentinel3A-OLCI at 300m of resolution) (at 1 km resolution) and the monthly averaged chlorophyll concentration for the multi-sensor (at 1 km resolution) and Sentinel-OLCI Level-3 (at 300m resolution). Chlorophyll field are obtained by means of the Mediterranean regional algorithms: an updated version of the MedOC4 (Case 1 waters, Volpe et al., 2019, with new coefficients) and AD4 (Case 2 waters, Berthon and Zibordi, 2004). Discrimination between the two water types is performed by comparing the satellite spectrum with the average water type spectral signature from in situ measurements for both water types. Reference insitu dataset is MedBiOp (Volpe et al., 2019) where pure Case II spectra are selected using a k-mean cluster analysis (Melin et al., 2015). Merging of Case I and Case II information is performed estimating the Mahalanobis distance between the observed and reference spectra and using it as weight for the final merged value. The interpolated gap-free Level-4 Chl concentration is estimated by means of a modified version of the DINEOF algorithm by GOS (Volpe et al., 2018). DINEOF is an iterative procedure in which EOF are used to reconstruct the entire field domain. As first guess, it uses the SeaWiFS-derived daily climatological values at missing pixels and satellite observations at valid pixels. The other Level-4 dataset is the time averages of the L3 fields and includes the standard deviation and the number of observations in the monthly 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. Level-4 includes both monthly time averages and the daily-interpolated fields. Time averages are computed on the delayed-time data. The interpolated product starts from the L3 products at 1 km resolution. At the first iteration, DINEOF procedure uses, as first guess for each of the missing pixels the relative daily climatological pixel. A procedure to smooth out spurious spatial gradients is applied to the daily merged image (observation and climatology). From the second iteration, the procedure uses, as input for the next one, the field obtained by the EOF calculation, using only a number of modes: that is, at the second round, only the first two modes, at the third only the first three, and so on. At each iteration, the same smoothing procedure is applied between EOF output and initial observations. The procedure stops when the variance explained by the current EOF mode exceeds that of noise. '''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 pre+D2sence 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-chl-multi-l4-chl_1km_monthly-rt-v02 *dataset-oc-med-chl-multi-l4-interp_1km_daily-rt-v02 *dataset-oc-med-chl-olci-l4-chl_300m_monthly-rt-v02 '''Files format:''' *CF-1.4 *INSPIRE compliant '''DOI (product) :''' https://doi.org/10.48670/moi-00113

'''This product has been archived''' '''DEFINITION''' The linear change of zonal mean subsurface temperature over the period 1993-2019 at each grid point (in depth and latitude) is evaluated to obtain a global mean depth-latitude plot of subsurface temperature trend, expressed in °C. The linear change is computed using the slope of the linear regression at each grid point scaled by the number of time steps (27 years, 1993-2019). A multi-product approach is used, meaning that the linear change is first computed for 5 different zonal mean temperature estimates. The average linear change is then computed, as well as the standard deviation between the five linear change computations. The evaluation method relies in the study of the consistency in between the 5 different estimates, which provides a qualitative estimate of the robustness of the indicator. See Mulet et al. (2018) for more details. '''CONTEXT''' Large-scale temperature variations in the upper layers are mainly related to the heat exchange with the atmosphere and surrounding oceanic regions, while the deeper ocean temperature in the main thermocline and below varies due to many dynamical forcing mechanisms (Bindoff et al., 2019). Together with ocean acidification and deoxygenation (IPCC, 2019), ocean warming can lead to dramatic changes in ecosystem assemblages, biodiversity, population extinctions, coral bleaching and infectious disease, change in behavior (including reproduction), as well as redistribution of habitat (e.g. Gattuso et al., 2015, Molinos et al., 2016, Ramirez et al., 2017). Ocean warming also intensifies tropical cyclones (Hoegh-Guldberg et al., 2018; Trenberth et al., 2018; Sun et al., 2017). '''CMEMS KEY FINDINGS''' The results show an overall ocean warming of the upper global ocean over the period 1993-2019, particularly in the upper 300m depth. In some areas, this warming signal reaches down to about 800m depth such as for example in the Southern Ocean south of 40°S. In other areas, the signal-to-noise ratio in the deeper ocean layers is less than two, i.e. the different products used for the ensemble mean show weak agreement. However, interannual-to-decadal fluctuations are superposed on the warming signal, and can interfere with the warming trend. For example, in the subpolar North Atlantic decadal variations such as the so called ‘cold event’ prevail (Dubois et al., 2018; Gourrion et al., 2018), and the cumulative trend over a quarter of a decade does not exceed twice the noise level below about 100m depth. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00244

'''This product has been archived''' This dataset provide a times series of gap free map of Sea Surface Temperature (SST) foundation at high resolution on a 0.10 x 0.10 degree grid (approximately 10 x 10 km) for the Global Ocean, every 24 hours. Whereas along swath observation data essentially represent the skin or sub-skin SST, the Level 4 SST product is defined to represent the SST foundation (SSTfnd). SSTfnd is defined within GHRSST as the temperature at the base of the diurnal thermocline. It is so named because it represents the foundation temperature on which the diurnal thermocline develops during the day. SSTfnd changes only gradually along with the upper layer of the ocean, and by definition it is independent of skin SST fluctuations due to wind- and radiation-dependent diurnal stratification or skin layer response. It is therefore updated at intervals of 24 hrs. SSTfnd corresponds to the temperature of the upper mixed layer which is the part of the ocean represented by the top-most layer of grid cells in most numerical ocean models. It is never observed directly by satellites, but it comes closest to being detected by infrared and microwave radiometers during the night, when the previous day's diurnal stratification can be assumed to have decayed. The processing combines the observations of multiple polar orbiting and geostationary satellites, embedding infrared of microwave radiometers. All these sources are intercalibrated with each other before merging. A ranking procedure is used to select the best sensor observation for each grid point. An optimal interpolation is used to fill in where observations are missing. '''DOI (product) :''' https://doi.org/10.48670/mds-00321

'''This product has been archived''' '''DEFINITION''' The temporal evolution of thermosteric sea level in an ocean layer is obtained from an integration of temperature driven ocean density variations, which are subtracted from a reference climatology to obtain the fluctuations from an average field. The regional thermosteric sea level values are then averaged from 60°S-60°N aiming to monitor interannual to long term global sea level variations caused by temperature driven ocean volume changes through thermal expansion as expressed in meters (m). '''CONTEXT''' Most of the interannual variability and trends in regional sea level is caused by changes in steric sea level. At mid and low latitudes, the steric sea level signal is essentially due to temperature changes, i.e. the thermosteric effect (Stammer et al., 2013, Meyssignac et al., 2016). Salinity changes play only a local role. Regional trends of thermosteric sea level can be significantly larger compared to their globally averaged versions (Storto et al., 2018). Except for shallow shelf sea and high latitudes (> 60° latitude), regional thermosteric sea level variations are mostly related to ocean circulation changes, in particular in the tropics where the sea level variations and trends are the most intense over the last two decades. '''CMEMS KEY FINDINGS''' Significant (i.e. when the signal exceeds the noise) regional trends for the period 2005-2019 from the Copernicus Marine Service multi-ensemble approach show a thermosteric sea level rise at rates ranging from the global mean average up to more than 8 mm/year. There are specific regions where a negative trend is observed above noise at rates up to about -8 mm/year such as in the subpolar North Atlantic, or the western tropical Pacific. These areas are characterized by strong year-to-year variability (Dubois et al., 2018; Capotondi et al., 2020). Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00241

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''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-2020 period was 0.59% 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 hemisphehres. The significant increases in chlorophyll reported in 2016-2017 (Sathyendranath et al., 2018b) for the Atlantic and Pacific oceans at high latitudes continued to be observed after the 2020 extension, as well as the negative trends over the equatorial Pacific and the Indian Ocean Gyre. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''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:''' The KD490 product identifies the turbidity of the water column, i.e., how visible light in the blue-green region of the spectrum penetrates within the water column. It is directly related to the presence of scattering particles in the water column. This product is derived from OLCI and remapped at nominal 300m spatial resolution using cylindrical equirectangular projection. '''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. By comparing reflectances at different wavelengths and calibrating the result against in-situ measurements, an estimate of in water absorption parameters can be derived. '''Quality / Accuracy / Calibration information:''' Detailed description of cal/val is given in the relevant QUID, associated validation reports and quality documentation. '''DOI (product) :''' https://doi.org/10.48670/moi-00078

'''DEFINITION''' The Copernicus Marine IBI_OMI_seastate_extreme_var_swh_mean_and_anomaly OMI indicator is based on the computation of the annual 99th percentile of Significant Wave Height (SWH) from model data. Two different CMEMS products are used to compute the indicator: The Iberia-Biscay-Ireland Multi Year Product (IBI_MULTIYEAR_WAV_005_006) and the Analysis product (IBI_ANALYSISFORECAST_WAV_005_005). Two parameters have been considered for this OMI: * Map of the 99th mean percentile: It is obtained from the Multi-Year Product, the annual 99th percentile is computed for each year of the product. The percentiles are temporally averaged in the whole period (1980-2023). * Anomaly of the 99th percentile in 2024: The 99th percentile of the year 2024 is computed from the Analysis product. The anomaly is obtained by subtracting the mean percentile to the percentile in 2024. This indicator is aimed at monitoring the extremes of annual significant wave height and evaluate the spatio-temporal variability. The use of percentiles instead of annual maxima, makes this extremes study less affected by individual data. This approach was first successfully applied to sea level variable (Pérez Gómez et al., 2016) and then extended to other essential variables, such as sea surface temperature and significant wave height (Pérez Gómez et al 2018 and Álvarez-Fanjul et al., 2019). Further details and in-depth scientific evaluation can be found in the CMEMS Ocean State report (Álvarez- Fanjul et al., 2019). '''CONTEXT''' The sea state and its related spatio-temporal variability affect dramatically maritime activities and the physical connectivity between offshore waters and coastal ecosystems, impacting therefore on the biodiversity of marine protected areas (González-Marco et al., 2008; Savina et al., 2003; Hewitt, 2003). Over the last decades, significant attention has been devoted to extreme wave height events since their destructive effects in both the shoreline environment and human infrastructures have prompted a wide range of adaptation strategies to deal with natural hazards in coastal areas (Hansom et al., 2015). Complementarily, there is also an emerging question about the role of anthropogenic global climate change on present and future extreme wave conditions (Young and Ribal, 2019). The Iberia-Biscay-Ireland region, which covers the North-East Atlantic Ocean from Canary Islands to Ireland, is characterized by two different sea state wave climate regions: whereas the northern half, impacted by the North Atlantic subpolar front, is of one of the world’s greatest wave generating regions (Mørk et al., 2010; Folley, 2017), the southern half, located at subtropical latitudes, is by contrast influenced by persistent trade winds and thus by constant and moderate wave regimes. The North Atlantic Oscillation (NAO), which refers to changes in the atmospheric sea level pressure difference between the Azores and Iceland, is a significant driver of wave climate variability in the Northern Hemisphere. The influence of North Atlantic Oscillation on waves along the Atlantic coast of Europe is particularly strong in and has a major impact on northern latitudes wintertime (Gleeson et al., 2017; Martínez-Asensio et al. 2016; Wolf et al., 2002; Bauer, 2001; Kushnir et al., 1997; Bouws et al., 1996; Bacon and Carter, 1991). Swings in the North Atlantic Oscillation index produce changes in the storms track and subsequently in the wind speed and direction over the Atlantic that alter the wave regime. When North Atlantic Oscillation index is in its positive phase, storms usually track northeast of Europe and enhanced westerly winds induce higher than average waves in the northernmost Atlantic Ocean. Conversely, in the negative North Atlantic Oscillation phase, the track of the storms is more zonal and south than usual, with trade winds (mid latitude westerlies) being slower and producing higher than average waves in southern latitudes (Marshall et al., 2001; Wolf et al., 2002; Wolf and Woolf, 2006). Additionally, a variety of previous studies have uniquevocally determined the relationship between the sea state variability in the IBI region and other atmospheric climate modes such as the East Atlantic pattern, the Arctic Oscillation, the East Atlantic Western Russian pattern and the Scandinavian pattern (Izaguirre et al., 2011, Martínez-Asensio et al., 2016). In this context, long‐term statistical analysis of reanalyzed model data is mandatory not only to disentangle other driving agents of wave climate but also to attempt inferring any potential trend in the number and/or intensity of extreme wave events in coastal areas with subsequent socio-economic and environmental consequences. '''CMEMS KEY FINDINGS''' The climatic mean of 99th percentile (1980-2023) reveals a north-south gradient of Significant Wave Height with the highest values in northern latitudes (above 8m) and lowest values (2-3 m) detected southeastward of Canary Islands, in the seas between Canary Islands and the African Continental Shelf. This north-south pattern is the result of the two climatic conditions prevailing in the region and previously described. The 99th percentile anomalies in 2024 show that during this period, virtually the entire IBI region was affected by positive anomalies in maximum SWH, which exceeded the standard deviation of the historical record in the waters west of the Iberian Peninsula, the Spanish coast of the Bay of Biscay, and the African coast south of Cape Ghir. Anomalies reaching twice the standard deviation of the time series were also observed in coastal regions of the English Channel. '''DOI (product):''' https://doi.org/10.48670/moi-00249

'''Short description: ''' For the Global Ocean - In-situ observation yearly delivery in delayed mode of Ocean surface currents. '''Detailed description: ''' The In Situ delayed mode product designed for reanalysis purposes integrates the best available version of in situ data for Ocean surface currents. The data are collected from the Surface Drifter Data Assembly Centre (SD-DAC at NOAA AOML) completed by European data provided by EUROGOOS regional systems and national systems by the regional INS TAC components. All surface drifters data have been processed to check for drogue loss. Drogued and undrogued drifting buoy surface ocean currents are provided with a drogue presence flag as well as a wind slippage correction for undrogued buoy. '''Processing information: ''' From the near real time INS TAC product validated on a daily and weekly basis for forecasting purposes, and from the SD-DAC quality controlled dataset a scientifically validated product is created . It s a """"reference product"""" updated on a yearly basis. This product has been processed using a method that checks for drogue loss. Altimeter and wind data have been used to extract the direct wind slippage from the total drifting buoy velocities. The obtained wind slippage values have then been analyzed to identify probable undrogued data among the drifting buoy velocities dataset. A simple procedure has then been applied to produce an updated dataset including a drogue presence flag as well as a wind slippage correction. '''Suitability, Expected type of users / uses: ''' The product is designed to be assimilated into or for validation purposes of operational models operated by ocean forecasting centers for reanalysis purposes or for research community. These users need data aggregated and quality controlled in a reliable and documented manner.