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''' Variations of the Mediterranean Outflow Water at 1000 m depth are monitored through area-averaged salinity anomalies in specifically defined boxes. The salinity data are extracted from several CMEMS products and averaged in the corresponding monitoring domain: * IBI-MYP: IBI_MULTIYEAR_PHY_005_002 * IBI-NRT: IBI_ANALYSISFORECAST_PHYS_005_001 * GLO-MYP: GLOBAL_REANALYSIS_PHY_001_030 * CORA: INSITU_GLO_TS_REP_OBSERVATIONS_013_002_b * ARMOR: MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012 The anomalies of salinity have been computed relative to the monthly climatology obtained from IBI-MYP. Outcomes from diverse products are combined to deliver a unique multi-product result. Multi-year products (IBI-MYP, GLO,MYP, CORA, and ARMOR) are used to show an ensemble mean and the standard deviation of members in the covered period. The IBI-NRT short-range product is not included in the ensemble, but used to provide the deterministic analysis of salinity anomalies in the most recent year. '''CONTEXT''' The Mediterranean Outflow Water is a saline and warm water mass generated from the mixing processes of the North Atlantic Central Water and the Mediterranean waters overflowing the Gibraltar sill (Daniault et al., 1994). The resulting water mass is accumulated in an area west of the Iberian Peninsula (Daniault et al., 1994) and spreads into the North Atlantic following advective pathways (Holliday et al. 2003; Lozier and Stewart 2008, de Pascual-Collar et al., 2019). The importance of the heat and salt transport promoted by the Mediterranean Outflow Water flow has implications beyond the boundaries of the Iberia-Biscay-Ireland domain (Reid 1979, Paillet et al. 1998, van Aken 2000). For example, (i) it contributes substantially to the salinity of the Norwegian Current (Reid 1979), (ii) the mixing processes with the Labrador Sea Water promotes a salt transport into the inner North Atlantic (Talley and MacCartney, 1982; van Aken, 2000), and (iii) the deep anti-cyclonic Meddies developed in the African slope is a cause of the large-scale westward penetration of Mediterranean salt (Iorga and Lozier, 1999). Several studies have demonstrated that the core of Mediterranean Outflow Water is affected by inter-annual variability. This variability is mainly caused by a shift of the MOW dominant northward-westward pathways (Bozec et al. 2011), it is correlated with the North Atlantic Oscillation (Bozec et al. 2011) and leads to the displacement of the boundaries of the water core (de Pascual-Collar et al., 2019). The variability of the advective pathways of MOW is an oceanographic process that conditions the destination of the Mediterranean salt transport in the North Atlantic. Therefore, monitoring the Mediterranean Outflow Water variability becomes decisive to have a proper understanding of the climate system and its evolution (e.g. Bozec et al. 2011, Pascual-Collar et al. 2019). The CMEMS IBI-OMI_WMHE_mow product is aimed to monitor the inter-annual variability of the Mediterranean Outflow Water in the North Atlantic. The objective is the establishment of a long-term monitoring program to observe the variability and trends of the Mediterranean water mass in the IBI regional seas. To do that, the salinity anomaly is monitored in key areas selected to represent the main reservoir and the three main advective spreading pathways. More details and a full scientific evaluation can be found in the CMEMS Ocean State report Pascual et al., 2018 and de Pascual-Collar et al. 2019. '''CMEMS KEY FINDINGS''' The absence of long-term trends in the monitoring domain Reservoir (b) suggests the steadiness of water mass properties involved on the formation of Mediterranean Outflow Water. Results obtained in monitoring box North (c) present an alternance of periods with positive and negative anomalies. The last negative period started in 2016 reaching up to the present. Such negative events are linked to the decrease of the northward pathway of Mediterranean Outflow Water (Bozec et al., 2011), which appears to return to steady conditions in 2020 and 2021. Results for box West (d) reveal a cycle of negative (2015-2017) and positive (2017 up to the present) anomalies. The positive anomalies of salinity in this region are correlated with an increase of the westward transport of salinity into the inner North Atlantic (de Pascual-Collar et al., 2019), which appear to be maintained for years 2020-2021. Results in monitoring boxes North and West are consistent with independent studies (Bozec et al., 2011; and de Pascual-Collar et al., 2019), suggesting a westward displacement of Mediterranean Outflow Water and the consequent contraction of the northern boundary. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00258

'''This product has been archived''' '''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 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''' Regional trends for the period 2005-2019 from the Copernicus Marine Service multi-ensemble approach show warming at rates ranging from the global mean average up to more than 8 W/m2 in some specific regions (e.g. northern hemisphere western boundary current regimes). There are specific regions where a negative trend is observed above noise at rates up to about -5 W/m2 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-00236

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

'''DEFINITION''' The indicator of Volume Transport Anomaly in Selected Vertical Sections in the Iberia–Biscay–Ireland (IBI) region (OMI_CIRCULATION_VOLTRANS_IBI_section_integrated_anomalies) is defined as the time series of annual mean volume transport calculated across a set of vertical ocean sections. These sections have been selected to represent the temporal variability of key ocean currents within the IBI domain. The monitored ocean currents include the transport towards the North Sea through the Rockall Trough (RTE) (Holliday et al., 2008; Lozier and Stewart, 2008), the Canary Current (CC) (Knoll et al., 2002; Mason et al., 2011), the Azores Current (AC) (Mason et al., 2011), the Algerian Current (ALG) (Tintoré et al., 1988; Benzohra and Millot, 1995; Font et al., 1998), and the net transport along the 48° N latitude parallel (N48) (see OMI figure). To produce ensemble-based results, six datasets provided by the Copernicus Marine Service have been used: * '''IBI-REA''' & '''IBI-INT''': IBI_MULTIYEAR_PHY_005_002 (reanalysis and interim datasets) * '''GLO-REA''': GLOBAL_MULTIYEAR_PHY_001_030 (reanalysis) * '''ARMOR''': MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012 (reprocessed observations) * '''MED-REA''': MEDSEA_MULTIYEAR_PHY_006_004 (reanalysis) * '''NWS-REA''': NWSHELF_MULTIYEAR_PHY_004_009 (reanalysis) The time series displays the ensemble mean (blue line), the ensemble spread (grey shaded area), and the mean transport with reversed sign (red dashed line), which indicates the threshold of anomaly values corresponding to a reversal in the direction of the current transport. In addition, the trend analysis at the 95% confidence level is shown in the bottom-right corner of each diagram. Further details on the product are provided in the corresponding Product User Manual (de Pascual-Collar et al., 2026a) and Quality Information Document (de Pascual-Collar et al., 2026b), as well as in de Pascual-Collar et al., 2024. '''CONTEXT''' The IBI area is a highly complex region characterized by a remarkable variety of ocean currents. Among them, we can highlight those that originate as a result of the closure of the North Atlantic Drift (Mason et al., 2011; Holliday et al., 2008; Peliz et al., 2007; Bower et al., 2002; Knoll et al., 2002; Pérez et al., 2001; Jia, 2000); the subsurface currents flowing northward along the continental slope (de Pascual-Collar et al., 2019; Pascual et al., 2018; Machín et al., 2010; Fricourt et al., 2007; Knoll et al., 2002; Mazé et al., 1997; White & Bowyer, 1997); and the exchange currents occurring in the Strait of Gibraltar and the Alboran Sea (Sotillo et al., 2016; Font et al., 1998; Benzohra & Millot, 1995; Tintoré et al., 1988). The variability of ocean currents in the IBI domain is relevant to the global thermohaline circulation and other climatic and environmental processes. For example, as discussed by Fasullo and Trenberth (2008), subtropical gyres play a crucial role in the meridional energy balance. The poleward salt transport of Mediterranean water, driven by subsurface slope currents, has significant implications for salinity anomalies in the Rockall Trough and the Nordic Seas, as studied by Holliday (2003), Holliday et al. (2008), and Bozec et al. (2011). The Algerian Current serves as the only pathway for Atlantic Water to reach the Western Mediterranean. '''CMEMS KEY FINDINGS''' The volume transport time series reveal periods during which the monitored currents exhibited notably high or low variability. Specifically, the RTE current shows pronounced variability in 2010 and during 2014–2015; the N48 section between 2012 and 2014; the ALG current in 2006 and 2017; the AC current between 2005–2007 and in 2021; and the CC current between 2005–2007. Furthermore, certain periods display anomalies of sufficient magnitude (in absolute value) to indicate a reversal in the net transport direction of the current. This is the case for the ALG current in 2017 and 2024 (with net transport towards the west), and for the CC current in 2010 (with net transport towards the north). Trend analysis over the period 1993–2023 does not reveal any statistically significant trends for the monitored currents. However, the confidence interval for the trend in the ALG section is close to rejecting the null hypothesis of no trend. '''DOI (product):''' https://doi.org/10.48670/mds-00351

'''DEFINITION''' The temporal evolution of thermosteric sea level in an ocean layer (here: 0-700m) is obtained from an integration of temperature driven ocean density variations, which are subtracted from a reference climatology (here 1993-2014) to obtain the fluctuations from an average field. The annual mean thermosteric sea level of the year 2017 is substracted from a reference climatology (1993-2014) at each grid point to obtain a global map of thermosteric sea level anomalies in the year 2017, expressed in millimeters per year (mm/yr). '''CONTEXT''' Most of the interannual variability and trends in regional sea level is caused by changes in steric sea level (Oppenheimer et al., 2019). 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''' Higher-than-average thermosteric sea level is reported over most areas of the global ocean and the European regional seas in 2018. In some areas – e.g. the western boundary current regions of the Pacific and Atlantic Ocean in both hemispheres reach values of more than 0.2 m. There are two areas of lower-than-average thermosteric sea level, which stand out from the generally higher-than-average conditions: the western tropical Pacific, and the subpolar North Atlantic. The latter is linked to the so called “North Atlantic cold event” which persists since a couple of years (Dubois et al., 2018). However, its signature has significantly reduced compared to preceding years.

'''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''' The global mean sea level is reflecting changes in the Earth’s climate system in response to natural and anthropogenic forcing factors such as ocean warming, land ice mass loss and changes in water storage in continental river basins. Thermosteric sea-level variations result from temperature related density changes in sea water associated with volume expansion and contraction. Global thermosteric sea level rise caused by ocean warming is known as one of the major drivers of contemporary global mean sea level rise (Cazenave et al., 2018; Oppenheimer et al., 2019). '''CMEMS KEY FINDINGS''' Since the year 2005 the upper (0-2000m) near-global (60°S-60°N) thermosteric sea level rises at a rate of 1.3±0.2 mm/year. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00240

'''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''' The global mean sea level is reflecting changes in the Earth’s climate system in response to natural and anthropogenic forcing factors such as ocean warming, land ice mass loss and changes in water storage in continental river basins. Thermosteric sea-level variations result from temperature related density changes in sea water associated with volume expansion and contraction. Global thermosteric sea level rise caused by ocean warming is known as one of the major drivers of contemporary global mean sea level rise (Cazenave et al., 2018; Oppenheimer et al., 2019). '''CMEMS KEY FINDINGS''' Since the year 2005 the upper (0-700m) near-global (60°S-60°N) thermosteric sea level rises at a rate of 0.9±0.1 mm/year. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00239

'''DEFINITION''' Time mean meridional Eulerian streamfunctions are computed using the velocity field estimate provided by the Copernicus Marine Mediterranean Sea reanalysis over the period from 1987 to the year preceding the current one [-1Y], operationally extended yearly. The Eulerian meridional streamfunction is evaluated by integrating meridional velocity daily data first in a vertical direction, then in a meridional direction, and finally averaging over the reanalysis period. The Mediterranean overturning indices are derived for the eastern and western Mediterranean Sea by computing the annual streamfunction in the two areas separated by the Strait of Sicily around 36.5°N, and then considering the associated maxima. In each case a geographical constraint focused the computation on the main region of interest. For the western index, we focused on deep-water formation regions, thus excluding both the effect of shallow physical processes and the Gibraltar net inflow. For the eastern index, we investigate the Levantine and Cretan areas corresponding to the strongest meridional overturning cell locations, thus only a zonal constraint is defined. Time series of annual mean values is provided for the Mediterranean Sea using the Mediterranean 1/24o eddy resolving reanalysis (Escudier et al., 2020, 2021). More details can be found in the Copernicus Marine Ocean State Report issue 4 (OSR4, von Schuckmann et al., 2020) Section 2.4 (Lyubartsev et al., 2020) and in the QUID. '''CONTEXT''' The western and eastern Mediterranean clockwise meridional overturning circulation is connected to deep-water formation processes. The Mediterranean Sea 1/24o eddy resolving reanalysis (MEDSEA_MULTIYEAR_PHY_006_004, Escudier et al., 2020, 2021) is used to show the interannual variability of the Meridional Overturning Index. Details on the product are delivered in the PUM and QUID of this OMI. The Mediterranean Meridional Overturning Index is defined here as the maxima of the clockwise cells in the eastern and western Mediterranean Sea and is associated with deep and intermediate water mass formation processes that occur in specific areas of the basin: Gulf of Lion, Southern Adriatic Sea, Cretan Sea and Rhodes Gyre (Pinardi et al., 2015). As in the global ocean, the overturning circulation of the western and eastern Mediterranean are paramount to determine the stratification of the basins (Cessi, 2019). In turn, the stratification and deep water formation mediate the exchange of oxygen and other tracers between the surface and the deep ocean (e.g., Johnson et al., 2009; Yoon et al., 2018). In this sense, the overturning indices are potential gauges of the ecosystem health of the Mediterranean Sea, and in particular they could instruct early warning indices for the Mediterranean Sea to support the Sustainable Development Goal (SDG) 13 Target 13.3. '''CMEMS KEY FINDINGS''' The western and eastern Mediterranean overturning indices (WMOI and EMOI) are synthetic indices of changes in the thermohaline properties of the Mediterranean basin related to changes in the main drivers of the basin scale circulation. The western sub-basin clockwise overturning circulation is associated with the deep-water formation area of the Gulf of Lion, while the eastern clockwise meridional overturning circulation is composed of multiple cells associated with different intermediate and deep-water sources in the Levantine, Aegean, and Adriatic Seas. On average, the EMOI shows higher values than the WMOI indicating a more vigorous overturning circulation in eastern Mediterranean. The difference is mostly related to the occurrence of the eastern Mediterranean transient (EMT) climatic event, and linked to a peak of the EMOI in 1992 (Roether et al. 1996, 2014, Gertman et al. 2006). In 1999, the difference between the two indices started to decrease because EMT water masses reached the Sicily Strait flowing into the western Mediterranean Sea (Schroeder et al., 2016). The western peak in 2006 is discussed to be linked to anomalous deep-water formation during the Western Mediterranean Transition (Smith, 2008; Schroeder et al., 2016). Thus, the WMOI and EMOI indices are a useful tool for long-term climate monitoring of overturning changes in the Mediterranean Sea. '''DOI (product):''' https://doi.org/10.48670/mds-00317

'''DEFINITION''' The Mediterranean water mass formation rates are evaluated in 4 areas as defined in the Ocean State Report issue 2 (OSR2, von Schuckmann et al., 2018) section 3.4 (Simoncelli and Pinardi, 2018): (1) the Gulf of Lions for the Western Mediterranean Deep Waters (WMDW); (2) the Southern Adriatic Sea Pit for the Eastern Mediterranean Deep Waters (EMDW); (3) the Cretan Sea for Cretan Intermediate Waters (CIW) and Cretan Deep Waters (CDW); (4) the Rhodes Gyre, the area of formation of the so-called Levantine Intermediate Waters (LIW) and Levantine Deep Waters (LDW). Annual water mass formation rates have been computed using daily mixed layer depth estimates (density criteria Δσ = 0.01 kg/m3, 10 m reference level) considering the annual maximum volume of water above mixed layer depth with potential density within or higher the specific thresholds specified in Table 1 then divided by seconds per year. Annual mean values are provided using the Mediterranean 1/24o eddy resolving reanalysis (Escudier et al. 2020, 2021). Time spans from 1987 to the year preceding the current one [-1Y], operationally extended yearly. '''CONTEXT''' The formation of intermediate and deep water masses is one of the most important processes occurring in the Mediterranean Sea, being a component of its general overturning circulation. This circulation varies at interannual and multidecadal time scales and it is composed of an upper zonal cell (Zonal Overturning Circulation) and two main meridional cells in the western and eastern Mediterranean (Pinardi and Masetti 2000). The objective is to monitor the main water mass formation events using the eddy resolving Mediterranean Sea Reanalysis (MEDSEA_MULTIYEAR_PHY_006_004, Escudier et al. 2020, 2021) and considering Pinardi et al. (2015) and Simoncelli and Pinardi (2018) as references for the methodology. The Mediterranean Sea Reanalysis can reproduce both Eastern Mediterranean Transient and Western Mediterranean Transition phenomena and catches the principal water mass formation events reported in the literature. This will permit constant monitoring of the open ocean deep convection process in the Mediterranean Sea and a better understanding of the multiple drivers of the general overturning circulation at interannual and multidecadal time scales. Deep and intermediate water formation events reveal themselves by a deep mixed layer depth distribution in four Mediterranean areas: Gulf of Lions, Southern Adriatic Sea Pit, Cretan Sea and Rhodes Gyre. '''KEY FINDINGS''' The Western Mediterranean Deep Water (WMDW) formation events in the Gulf of Lion appear to be larger after 1999 consistently with Schroeder et al. (2006, 2008) related to the Eastern Mediterranean Transient event. This modification of WMDW after 2005 has been called Western Mediterranean Transition. WMDW formation events are consistent with Somot et al. (2016) and the event in 2009 is also reported in Houpert et al. (2016). The Eastern Mediterranean Deep Water (EMDW) formation in the Southern Adriatic Pit region displays a period of water mass formation between 1988 and 1993, in agreement with Pinardi et al. (2015), in 1996, 1999 and 2000 as documented by Manca et al. (2002). Weak deep water formation in winter 2006 is confirmed by observations in Vilibić and Šantić (2008). An intense deep water formation event is detected in 2012-2013 (Gačić et al., 2014). Last years are characterized by large events starting from 2017 (Mihanovic et al., 2021). Cretan Intermediate Water formation rates present larger peaks between 1989 and 1993 with the ones in 1992 and 1993 composing the Eastern Mediterranean Transient phenomena. The Cretan Deep Water formed in 1992 and 1993 is characterized by the highest densities of the entire period in accordance with Velaoras et al. (2014). The Levantine Deep Water formation rate in the Rhode Gyre region presents the largest values between 1992 and 1993 in agreement with Kontoyiannis et al. (1999). '''DOI (product):''' https://doi.org/10.48670/mds-00318