he Global ARMOR3D L4 Reprocessed dataset is obtained by combining satellite (Sea Level Anomalies, Geostrophic Surface Currents, Sea Surface Temperature) and in-situ (Temperature and Salinity profiles) observations through statistical methods. References : - ARMOR3D: Guinehut S., A.-L. Dhomps, G. Larnicol and P.-Y. Le Traon, 2012: High resolution 3D temperature and salinity fields derived from in situ and satellite observations. Ocean Sci., 8(5):845–857. - ARMOR3D: Guinehut S., P.-Y. Le Traon, G. Larnicol and S. Philipps, 2004: Combining Argo and remote-sensing data to estimate the ocean three-dimensional temperature fields - A first approach based on simulated observations. J. Mar. Sys., 46 (1-4), 85-98. - ARMOR3D: Mulet, S., M.-H. Rio, A. Mignot, S. Guinehut and R. Morrow, 2012: A new estimate of the global 3D geostrophic ocean circulation based on satellite data and in-situ measurements. Deep Sea Research Part II : Topical Studies in Oceanography, 77–80(0):70–81.

'''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

'''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), the monthly averaged chlorophyll concentration for the multi-sensor and climatological fields, all at 1 km 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. Monthly 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. SeaWiFS daily climatology provides reference for the calculation of Quality Indices (QI) for Chl observations. '''Processing information:''' Multi-sensor product is 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. Single sensor NASA Level-2 data are destriped and then all Level-2 data are remapped at 1 km spatial resolution 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) 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. The entire data set is consistent and processed in one-shot mode (with an unique software version and identical configurations). For the climatology Rrs data are derived from the latest SeaWiFS NASA reprocessing (R2018.0) and converted to chlorophyll concentrations with the same algorithm as the one used for other L4 and/or L3 products. '''Description of observation methods/instruments:''' Ocean color technique exploits the emerging electromagnetic radiation from the sea surface in different wavelengths. The spectral variability of this signal defines the so called ocean color which is affected by the presence of phytoplankton. '''Quality / Accuracy / Calibration information:''' A detailed description of both the cal/val and a more in depth view of the method is given in QUID-009-038to045-071-073-078-079-095-096.pdf over Copernicus 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-rep-v02 *dataset-oc-med-chl-multi-l4-interp_1km_daily-rep *dataset-oc-med-chl-seawifs-l4-chl_1km_daily-climatology-v02 '''Files format:''' *CF-1.4 *INSPIRE compliant. '''DOI (product) :''' https://doi.org/10.48670/moi-00114

'''Short description''': The data are provided weekly over a regular grid at 1/4° horizontal resolution, from the surface to 1500 m depth (representative of each Wednesday). The velocities are obtained by solving a diabatic formulation of the Omega equation, starting from ARMOR3D data (MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012 ) and ERA5 surface fluxes. '''DOI (product) :''' https://doi.org/10.48670/moi-00053

'''Short description:''' Altimeter satellite along-track sea surface heights anomalies (SLA) computed with respect to a twenty-year [1993, 2012] mean with a 1Hz (~7km) sampling. It serves in delayed-time applications. This product is processed by the DUACS multimission altimeter data processing system. It processes data from all altimeter missions available (e.g. Sentinel-6A, Jason-3, Sentinel-3A, Sentinel-3B, Saral/AltiKa, Cryosat-2, Jason-1, Jason-2, Topex/Poseidon, ERS-1, ERS-2, Envisat, Geosat Follow-On, HY-2A, HY-2B, etc.). The system exploits the most recent datasets available based on the enhanced GDR/NTC production. All the missions are homogenized with respect to a reference mission. Part of the processing is fitted to the Global ocean. (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 Atmospheric Correction, Ocean Tides, Long Wavelength Errors) that can be used to change the physical content for specific needs (see PUM document for details) '''Associated products''' A time invariant product https://resources.marine.copernicus.eu/product-detail/SEALEVEL_GLO_PHY_NOISE_L4_STATIC_008_033/INFORMATION describing the noise level of along-track measurements is available. It is associated to the sla_filtered variable. It is a gridded product. One file is provided for the global ocean and those values must be applied for Arctic and Europe products. For Mediterranean and Black seas, one value is given in the QUID document. '''DOI (product)''': https://doi.org/10.48670/moi-00146

'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' The ocean monitoring indicator on regional mean sea level is derived from the DUACS delayed-time (DT-2021 version) altimeter gridded maps of sea level anomalies based on a stable number of altimeters (two) in the satellite constellation. These products are distributed by the Copernicus Climate Change Service and the Copernicus Marine Service (SEALEVEL_GLO_PHY_CLIMATE_L4_MY_008_057). The mean sea level evolution estimated in the Irish-Biscay-Iberian (IBI) region is derived from the average of the gridded sea level maps weighted by the cosine of the latitude. The annual and semi-annual periodic signals are removed (least square fit of sinusoidal function) and the time series is low-pass filtered (175 days cut-off). The curve is corrected for the regional mean effect of the Glacial Isostatic Adjustment (GIA) using the ICE5G-VM2 GIA model (Peltier, 2004). During 1993-1998, the Global men sea level (hereafter GMSL) has been known to be affected by a TOPEX-A instrumental drift (WCRP Global Sea Level Budget Group, 2018; Legeais et al., 2020). This drift led to overestimate the trend of the GMSL during the first 6 years of the altimetry record (about 0.04 mm/y at global scale over the whole altimeter period). A correction of the drift is proposed for the Global mean sea level (Legeais et al., 2020). Whereas this TOPEX-A instrumental drift should also affect the regional mean sea level (hereafter RMSL) trend estimation, currently this empirical correction is currently not applied to the altimeter sea level dataset and resulting estimated for RMSL. Indeed, the pertinence of the global correction applied at regional scale has not been demonstrated yet and there is no clear consensus achieved on the way to proceed at regional scale. Additionally, the estimation of such a correction at regional scale is not obvious, especially in areas where few accurate independent measurements (e.g. in situ)- necessary for this estimation - are available. The trend uncertainty is provided in a 90% confidence interval (Prandi et al., 2021). This estimate only considers errors related to the altimeter observation system (i.e., orbit determination errors, geophysical correction errors and inter-mission bias correction errors). The presence of the interannual signal can strongly influence the trend estimation considering to the altimeter period considered (Wang et al., 2021; Cazenave et al., 2014). The uncertainty linked to this effect is not taken into account. '''CONTEXT''' The indicator on area averaged sea level is a crucial index of climate change, and individual components contribute to sea level rise, including expansion due to ocean warming and melting of glaciers and ice sheets (WCRP Global Sea Level Budget Group, 2018). According to the recent IPCC 6th assessment report, global mean sea level (GMSL) increased by 0.20 (0.15 to 0.25) m over the period 1901 to 2018 with a rate 25 of rise that has accelerated since the 1960s to 3.7 (3.2 to 4.2) mm yr-1 for the period 2006–2018. Human activity was very likely the main driver of observed GMSL rise since 1970 (IPCC WGII, 2021). The weight of the different contributions evolves with time and in the recent decades the mass change has increased, contributing to the on-going acceleration of the GMSL trend (IPCC, 2022a; Legeais et al., 2020; Horwath et al., 2022). At regional scale, sea level does not change homogenously, and RMSL rise can also be influenced by various other processes, with different spatial and temporal scales, such as local ocean dynamic, atmospheric forcing, Earth gravity and vertical land motion changes (IPCC WGI, 2021). Rising sea level can strongly affect population and infrastructures in coastal areas, increase their vulnerability and risks for food security, particularly in low lying areas and island states. Adverse impacts from floods, storms and tropical cyclones with related losses and damages have increased due to sea level rise, and increase their vulnerability and increase risks for food security, particularly in low lying areas and island states (IPCC, 2022b). Adaptation and mitigation measures such as the restoration of mangroves and coastal wetlands, reduce the risks from sea level rise (IPCC, 2022c). In IBI region, the RMSL trend is modulated by decadal variations. As observed over the global ocean, the main actors of the long-term RMSL trend are associated with anthropogenic global/regional warming. Decadal variability is mainly linked to the strengthening or weakening of the Atlantic Meridional Overturning Circulation (AMOC) (e.g. Chafik et al., 2019). The latest is driven by the North Atlantic Oscillation (NAO) (e.g. Delworth and Zeng, 2016). Along the European coast, the NAO also influences the along-slope winds dynamic which in return significantly contributes to the local sea level variability observed (Chafik et al., 2019). '''CMEMS KEY FINDINGS''' Over the [1993/01/01, 2021/08/02] period, the basin-wide RMSL in the IBI area rises at a rate of 3.8  0.82 mm/year. '''DOI (product):''' https://doi.org/10.48670/moi-00252

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' For the Global ocean, the ESA Ocean Colour CCI surface Chlorophyll (mg m-3, 4 km resolution) using the OC-CCI recommended chlorophyll algorithm is made available in CMEMS format. L3 products are daily files, while the L4 are monthly composites. ESA-CCI data are provided by Plymouth Marine Laboratory at 4km resolution. These are processed using the same in-house software as in the operational processing. Standard masking criteria for detecting clouds or other contamination factors have been applied during the generation of the Rrs, i.e., land, cloud, sun glint, atmospheric correction failure, high total radiance, large solar zenith angle (actually a high air mass cutoff, but approximating to 70deg zenith), coccolithophores, negative water leaving radiance, and normalized water leaving radiance at 555 nm 0.15 Wm-2 sr-1 (McClain et al., 1995). 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 chlorophyll content can be derived. A detailed description of calibration & validation is given in the relevant QUID, associated validation reports and quality documentation. '''Processing information:''' ESA-CCI data are provided by Plymouth Marine Laboratory at 4km resolution. These are processed using the same in-house software as in the operational processing. The entire CCI data set is consistent and processing is done in one go. Both OC CCI and the REP product are versioned. Standard masking criteria for detecting clouds or other contamination factors have been applied during the generation of the Rrs, i.e., land, cloud, sun glint, atmospheric correction failure, high total radiance, large solar zenith angle (actually a high air mass cutoff, but approximating to 70deg zenith), coccolithophores, negative water leaving radiance, and normalized water leaving radiance at 555 nm 0.15 Wm-2 sr-1 (McClain et al., 1995). '''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 chlorophyll content can be derived. '''Quality / Accuracy / Calibration information:''' Detailed description of cal/val is given in the relevant QUID, associated validation reports and quality documentation. '''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. '''DOI (product) :''' https://doi.org/10.48670/moi-00101

'''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 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). 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-2023 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 -5 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

'''DEFINITION''' OMI_CLIMATE_SST_BAL_area_averaged_anomalies product includes time series of monthly mean SST anomalies over the period 1982-2024, relative to the 1991-2020 climatology, averaged for the Baltic Sea. The SST Level 4 analysis products that provide the input to the monthly averages are taken from the reprocessed product SST_BAL_SST_L4_REP_OBSERVATIONS_010_016 with a recent update to include 2023. The product has a spatial resolution of 0.02 in latitude and longitude. The OMI time series runs from Jan 1, 1982 to December 31, 2024 and is constructed by calculating monthly averages from the daily level 4 SST analysis fields of the SST_BAL_SST_L4_REP_OBSERVATIONS_010_016 product . The climatology period from 1991 to 2020 (30 years) is selected according to WMO recommendations (WMO, 2017) and the most recent practice from the U.S. National Oceanic and Atmospheric Administration practice (https://wmo.int/media/news/updated-30-year-reference-period-reflects-changing-climate). See the Copernicus Marine Service Ocean State Reports (section 1.1 in Von Schuckmann et al., 2016; section 3 in Von Schuckmann et al., 2018) for more information on the OMI product. '''CONTEXT''' Sea Surface Temperature (SST) is an Essential Climate Variable (GCOS) that is an important input for initialising numerical weather prediction models and fundamental for understanding air-sea interactions and monitoring climate change (GCOS 2010). The Baltic Sea is a region that requires special attention regarding the use of satellite SST records and the assessment of climatic variability (Høyer and She 2007; Høyer and Karagali 2016). The Baltic Sea is a semi-enclosed basin with natural variability and it is influenced by large-scale atmospheric processes and by the vicinity of land. In addition, the Baltic Sea is one of the largest brackish seas in the world. When analysing regional-scale climate variability, all these effects have to be considered, which requires dedicated regional and validated SST products. Satellite observations have previously been used to analyse the climatic SST signals in the North Sea and Baltic Sea (BACC II Author Team 2015; Lehmann et al. 2011). Recently, Høyer and Karagali (2016) demonstrated that the Baltic Sea had warmed 1-2 oC from 1982 to 2012 considering all months of the year and 3-5 °C when only July-September months were considered. This was corroborated in the Ocean State Reports (section 1.1 in Von Schuckmann et al., 2016; section 3 in Von Schuckmann et al., 2018). '''CMEMS KEY FINDINGS''' The basin-average trend of SST anomalies for Baltic Sea region amounts to 0.039±0.003°C/year over the period 1982-2024 which corresponds to an average warming of 1.68°C. Adding the North Sea area, the average trend amounts to 0.026±0.002°C/year over the same period, which corresponds to an average warming of 1.19°C for the entire region since 1982. '''DOI (product):''' https://doi.org/10.48670/moi-00205

'''Short description:''' The Reprocessed (REP) Mediterranean (MED) dataset provides a stable and consistent long-term Sea Surface Temperature (SST) time series over the Mediterranean Sea (and the adjacent North Atlantic box) developed for climate applications. This product consists of daily (nighttime), optimally interpolated (L4), satellite-based estimates of the foundation SST (namely, the temperature free, or nearly-free, of any diurnal cycle) at 0.05° resolution grid covering the period from 1st January 1981 to present (approximately one month before real time). The MED-REP-L4 product is built from a consistent reprocessing of the collated level-3 (merged single-sensor, L3C) climate data record (CDR) v.3.0, provided by the ESA Climate Change Initiative (CCI) and covering the period up to 2021, and its interim extension (ICDR) that allows the regular temporal extension for 2022 onwards. '''DOI (product) :''' https://doi.org/10.48670/moi-00173