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 global yearly ocean CO2 sink represents the ocean uptake of CO2 from the atmosphere computed over the whole ocean. It is expressed in PgC per year. The ocean monitoring index is presented for the period 1985 to year-1. The yearly estimate of the ocean CO2 sink corresponds to the mean of a 100-member ensemble of CO2 flux estimates (Chau et al. 2022). The range of an estimate with the associated uncertainty is then defined by the empirical 68% interval computed from the ensemble. '''CONTEXT''' Since the onset of the industrial era in 1750, the atmospheric CO2 concentration has increased from about 277±3 ppm (Joos and Spahni, 2008) to 412.44±0.1 ppm in 2020 (Dlugokencky and Tans, 2020). By 2011, the ocean had absorbed approximately 28 ± 5% of all anthropogenic CO2 emissions, thus providing negative feedback to global warming and climate change (Ciais et al., 2013). The ocean CO2 sink is evaluated every year as part of the Global Carbon Budget (Friedlingstein et al. 2022). The uptake of CO2 occurs primarily in response to increasing atmospheric levels. The global flux is characterized by a significant variability on interannual to decadal time scales largely in response to natural climate variability (e.g., ENSO) (Friedlingstein et al. 2022, Chau et al. 2022). '''CMEMS KEY FINDINGS''' The rate of change of the integrated yearly surface downward flux has increased by 0.04±0.01e-1 PgC/yr2 over the period 1985 to year-1. The yearly flux time series shows a plateau in the 90s followed by an increase since 2000 with a growth rate of 0.06±0.04e-1 PgC/yr2. In 2021 (resp. 2020), the global ocean CO2 sink was 2.41±0.13 (resp. 2.50±0.12) PgC/yr. The average over the full period is 1.61±0.10 PgC/yr with an interannual variability (temporal standard deviation) of 0.46 PgC/yr. In order to compare these fluxes to Friedlingstein et al. (2022), the estimate of preindustrial outgassing of riverine carbon of 0.61 PgC/yr, which is in between the estimate by Jacobson et al. (2007) (0.45±0.18 PgC/yr) and the one by Resplandy et al. (2018) (0.78±0.41 PgC/yr) needs to be added. A full discussion regarding this OMI can be found in section 2.10 of the Ocean State Report 4 (Gehlen et al., 2020) and in Chau et al. (2022). '''DOI (product):''' https://doi.org/10.48670/moi-00223

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' This RRS product is defined as the ratio of upwelling radiance and downwelling irradiance at 412, 443, 490, 510, 560 and 665 nm wavebands (corresponding to MERIS), and can also be expressed as the ratio of normalized water leaving Radiance (nLw) and the extra-terrestrial solar irradiance (F0). The ESA Climate Change Initiative is a 2-part programme aiming to produce “climate quality” merged data records from multiple sensors. The Ocean Colour project within this programme has a primary focus on chlorophyll in open oceans, using the highest quality Rrs merging process to date. This uses a combination of bandshifting to a reference sensor and temporally-weighted bias correction to align independent sensors into a coherent and minimally-biased set of reflectances. These are derived from level 2 data produced by SeaDAS l2gen (SeaWiFS) and Polymer (MODIS, VIIRS, MERIS and OLCI-3A) , and the resulting Rrs bias corrected. '''Processing information:''' ESA-CCI Rrs raw data are provided by Plymouth Marine Laboratory, currently at 4km resolution. These are processed to produce CMEMS representations 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 (70deg), large spacecraft zenith angle (56deg), coccolithophores, negative water leaving radiance, and normalized water leaving radiance at 560 nm 0.15 Wm-2 sr-1 (McClain et al., 1995). For the regional products, a variant of the OC-CCI chain is run to produce high resolution data at the 1km resolution necessary. '''DOI (product) :''' https://doi.org/10.48670/moi-00077

'''Short description:''' For the Global Ocean- In-situ observation yearly delivery in delayed mode. The In Situ delayed mode product designed for reanalysis purposes integrates the best available version of in situ data for temperature and salinity measurements. These data are collected from main global networks (Argo, GOSUD, OceanSITES, World Ocean Database) completed by European data provided by EUROGOOS regional systems and national system by the regional INS TAC components. It is updated on a yearly basis. The time coverage has been extended in the past by integration of EN4 data for the period 1950-1990. Acces through CMEMS Catalogue after registration: http://marine.copernicus.eu/ '''Detailed description: ''' Ocean circulation models need information on the interior of the ocean to be able to generate accurate forecast. This information is only available from in-situ measurement. However this information is acquired all around the world and not easily available to the operational users. Therefore, INS TAC , by connecting to a lot of international networks, collects, controls and disseminates the relevant in-situ data to operational users . For reanalysis purposes, operational centres needs to access to the best available datasets with the best possible coverage and where additional quality control procedures have been performed. This dataset suits research community needs Each year, a new release of this product is issued containing all the observations gathered by the INS TAC global component operated by Coriolis. '''Processing information:''' From the near real time INS TAC product validated on a daily and weekly basis for forecasting purposes, a scientifically validated product is created . It s a ""reference product"" updated on a yearly basis. This product has been controlled using an objective analysis (statistical tests) method and a visual quality control (QC). This QC procedure has been developed with the main objective to improve the quality of the dataset to the level required by the climate application and the physical ocean re-analysis activities. It provides T and S weekly gridded fields and individual profiles both on their original level and interpolated level. The measured parameters, depending on the data source, are : temperature, salinity. The reference level of measurements is immersion (in meters) or pressure (in decibars). The EN4 data were converted to the CORA NetCDF format without any additional validation. '''Quality/accuracy/calibration information:''' The process is done in two steps using two different time windows, corresponding to two runs of objective analysis, with an additional visual QC inserted between. The first run was done on a window of three weeks, to capture the most doubtful profiles which were then checked visually by an operator to decide whether or not it was bad data or real oceanic phenomena. The second run was done on a weekly basis to fit the modelling needs. '''Suitability, Expected type of users / uses: ''' The product is designed for assimilation into operational models operated by ocean forecasting centres for reanalysis purposes or for research community. These users need data aggregated and quality controlled in a reliable and documented manner.

'''DEFINITION''' The Iberia Biscay Ireland (IBI) Sea Surface Temperature extreme from Reanalysis ocean monitoring indicator (OMI) (OMI_CLIMATE_TEMPSAL_IBI_extreme_var_temp_mean_and_anomaly) is based on the computation of the annual 99th percentile of Sea Surface Temperature (SST) from model data. Two different Copernicus Marine products are used to compute the indicator: The IBI Reanalysis (IBI_MULTIYEAR_PHY_005_002) and the IBI Analysis product (IBI_ANALYSISFORECAST_PHY_005_001). Two parameters have been considered for this OMI: * '''Map of the 99th mean percentile''': It is obtained from the reanalysis product, the annual 99th percentile is computed for each year of the product. The percentiles are temporally averaged over the whole period (1993-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 from the 2024 percentile. This indicator is aimed at monitoring the extremes of sea surface temperature every year and at checking their variations in space. The use of percentiles instead of annual maxima, makes this extremes study less affected by individual data. This study of extreme variability was first applied to the 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 Alvarez Fanjul et al., 2019). More details and a full scientific evaluation can be found in the CMEMS Ocean State report (Alvarez Fanjul et al., 2019). '''CONTEXT''' The Sea Surface Temperature (SST) is one of the essential ocean variables, hence the monitoring of this variable is of key importance, since its variations can affect the ocean circulation, marine ecosystems, and ocean-atmosphere exchange processes. As the oceans continuously interact with the atmosphere, trends of sea surface temperature can also have an effect on the global climate. While the global-averaged sea surface temperatures have increased since the beginning of the 20th century (Hartmann et al., 2013) in the North Atlantic, anomalous cold conditions have also been reported since 2014 (Mulet et al., 2018; Dubois et al., 2018). The IBI area is a complex dynamic region with a remarkable variety of ocean physical processes and scales involved. The SST field in the region is strongly dependent on latitude, with higher values towards the South (Locarnini et al. 2013). This latitudinal gradient is supported by the presence of the eastern part of the North Atlantic subtropical gyre that transports cool water from the northern latitudes towards the equator. Additionally, the IBI region is under the influence of the Sea Level Pressure dipole established between the Icelandic low and the Bermuda high. Therefore, the interannual and interdecadal variability of the surface temperature field may be influenced by the North Atlantic Oscillation pattern (Czaja and Frankignoul, 2002; Flatau et al., 2003). Upwelling processes, taking place in the coastal margins, are also relevant in the IBI region. The most referenced one is the eastern boundary coastal upwelling system off the African and western Iberian coast (Sotillo et al., 2016), although other smaller upwelling systems have also been described in the northern coast of the Iberian Peninsula (Alvarez et al., 2011), the south-western Irish coast (Edwars et al., 1996) and the European Continental Slope (Dickson, 1980). '''CMEMS KEY FINDINGS''' In the IBI region, the 99th mean percentile for 1993-2023 shows a north-south pattern driven by the climatological distribution of temperatures in the North Atlantic. In the coastal regions of Africa and the Iberian Peninsula, the mean values are influenced by the upwelling processes (Sotillo et al., 2016). These results are consistent with the ones presented in Álvarez Fanjul (2019) for the period 1993-2016. The analysis of the 99th percentile SST anomaly for the year 2024 reveals that the northeastern Atlantic region, between latitudes 36° N and 48° N, experienced thermal anomalies exceeding twice the standard deviation. Similar anomalies are also observed near the northeastern Iberian Peninsula, suggesting that inshore and coastal areas may have been affected as well. In contrast, the upwelling region west of the Iberian Peninsula shows negative anomalies in maximum SST, indicating an intensification of upwelling processes in this area. '''DOI (product):''' https://doi.org/10.48670/moi-00254

'''DEFINITION''' This product includes the Mediterranean Sea satellite chlorophyll trend map based on regional chlorophyll reprocessed (MY) product as distributed by CMEMS OC-TAC (OCEANCOLOUR_MED_BGC_L3_NRT_009_141). This dataset, derived from multi-sensor (SeaStar-SeaWiFS, AQUA-MODIS, NOAA20-VIIRS, NPP-VIIRS, Envisat-MERIS and Sentinel3-OLCI) (at 1 km resolution) 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. (2023). The trend map is obtained by applying Colella et al. (2016) methodology, where the Mann-Kendall test (Mann, 1945; Kendall, 1975) and Sens’s method (Sen, 1968) are applied on deseasonalized monthly time series, as obtained from the X-11 technique (see e. g. Pezzulli et al. 2005), to estimate, trend magnitude and its significance. The trend is expressed in % per year that represents the relative changes (i.e., percentage) corresponding to the dimensional trend [mg m-3 y-1] with respect to the reference climatology (1997-2014). Only significant trends (p < 0.05) are included. This OMI has been introduced since the 2nd issue of Ocean State Report in 2017. '''CONTEXT''' Phytoplankton are key actors in the carbon cycle and, as such, recognised as an Essential Climate Variable (ECV). 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). The Mediterranean Sea is an oligotrophic basin, where chlorophyll concentration decreases following a specific gradient from West to East (Colella et al. 2016). The highest concentrations are observed in coastal areas and at the river mouths, where the anthropogenic pressure and nutrient loads impact on the eutrophication regimes (Colella et al. 2016). The the use of long-term time series of consistent, well-calibrated, climate-quality data record is crucial for detecting eutrophication. 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. '''KEY FINDINGS''' The chlorophyll trend in the Mediterranean Sea for the 1997–2024 period confirms the findings of the previous release, with predominantly negative values observed across most of the basin. On average, the region shows a trend of approximately -0.77% per year, slightly more negative than the overall trend reported previously. As in earlier assessments, weak positive trends persist in specific areas such as the northern Aegean Sea and the Sicily Channel. Compared to the 1997–2023 period, new positive trends are now evident in the Gulf of Lion. In contrast to the findings of Salgado-Hernanz et al. (2019), which were based on satellite observations from 1998 to 2014, this analysis does not reveal a distinct difference between the western and eastern Mediterranean basins. Notably, the Ligurian Sea now exhibits a negative trend, diverging from the positive trends identified by Colella et al. (2016) for the 1998–2009 period and by Salgado-Hernanz et al. (2019) for 1998–2014. Similarly, the waters of the Northern Adriatic Sea show weak positive trends, differing from the strong negative trend previously reported by Colella et al. (2016), and also representing a shift from the positive values observed by Salgado-Hernanz et al. (2019). '''DOI (product):''' https://doi.org/10.48670/moi-00260

'''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''' The CMEMS NORTHWESTSHELF_OMI_tempsal_extreme_var_temp_mean_and_anomaly OMI indicator is based on the computation of the annual 99th percentile of Sea Surface Temperature (SST) from model data. Two different CMEMS products are used to compute the indicator: The North-West Shelf Multi Year Product (NWSHELF_MULTIYEAR_PHY_004_009) and the Analysis product (NORTHWESTSHELF_ANALYSIS_FORECAST_PHY_004_013). Two parameters are included on 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 over the whole period (1993-2019). * Anomaly of the 99th percentile in 2020: The 99th percentile of the year 2020 is computed from the Analysis product. The anomaly is obtained by subtracting the mean percentile from the 2020 percentile. This indicator is aimed at monitoring the extremes of sea surface temperature every year and at checking their variations in space. The use of percentiles instead of annual maxima, makes this extremes study less affected by individual data. This study of extreme variability was first applied to the 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 Alvarez Fanjul et al., 2019). More details and a full scientific evaluation can be found in the CMEMS Ocean State report (Alvarez Fanjul et al., 2019). '''CONTEXT''' This domain comprises the North West European continental shelf where depths do not exceed 200m and deeper Atlantic waters to the North and West. For these deeper waters, the North-South temperature gradient dominates (Liu and Tanhua, 2021). Temperature over the continental shelf is affected also by the various local currents in this region and by the shallow depth of the water (Elliott et al., 1990). Atmospheric heat waves can warm the whole water column, especially in the southern North Sea, much of which is no more than 30m deep (Holt et al., 2012). Warm summertime water observed in the Norwegian trench is outflow heading North from the Baltic Sea and from the North Sea itself. '''CMEMS KEY FINDINGS''' The 99th percentile SST product can be considered to represent approximately the warmest 4 days for the sea surface in Summer. Maximum anomalies for 2020 are up to 4oC warmer than the 1993-2019 average in the western approaches, Celtic and Irish Seas, English Channel and the southern North Sea. For the atmosphere, Summer 2020 was exceptionally warm and sunny in southern UK (Kendon et al., 2021), with heatwaves in June and August. Further north in the UK, the atmosphere was closer to long-term average temperatures. Overall, the 99th percentile SST anomalies show a similar pattern, with the exceptional warm anomalies in the south of the domain. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product)''' https://doi.org/10.48670/moi-00273

'''This product has been archived''' "''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-2022), following the Vantrepotte and Mélin (2009) method. The trend estimate is not shown because the length of the time series does not allow to completely differentiate the climate trend to the natural variability of the primary production. 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 presents 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 demonstrates 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''' Global primary production does not show specific trend and remain relatively constant over the archive 1998-2022. The temporal variability of the primary production appears to be mainly driven by the seasonal variation. However, some specific inter-annual event may induce noticeable increase or decrease in primary production, as for example in the second part of 2011. '''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:''' 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. Processing of these data was mainly carried out in the OC-CCI framework producing a climate-quality consistent dataset using the latest and most complete knowledge of satellite sensor calibration, characterization and attitude, complete (as far as possible) ancillary data sets, latest versions of models and algorithms etc. The data were then repackaged, using custom software, to suit the requirements of CMEMS. The remote sensing of Ocean Colour represents a measure of the spectral variations in the light leaving the water surface, subsequently interpreted in terms of concentrations of optically-significant constituents in the water. The electromagnetic signal collected by the sensor on-board the satellite is largely determined by photons that have never reached the water surface, but have been backscattered within the atmosphere through multiple interactions between gas molecules and aerosols. After removing the atmospheric contribution, the water leaving radiance recorded at a given time by the satellite reflects the optical properties of the water which, in turn, mirrors a specific structure and biogeochemical composition of the marine waters. A detailed description of calibration & validation is given in the relevant QUID, associated validation reports and quality documentation. '''How to reference product:''' The User will ensure that the original product OCEANCOLOUR_GLO_OPTICS_L3_REP_OBSERVATIONS_009_064 -or value added products or derivative works developed from it including pictures- shall credit CMEMS and ESA/CCI by explicitly making mention of the originator in the following manner: ""Generated using Copernicus Marine and ESA/CCI Product"". For publication purposes, the User shall ensure that the credits mention CMEMS and ESA/CCI in the following manner: ""This study has been conducted using Copernicus Marine and ESA/CCI Product"". For all detailed information concerning the use of this product, see the Service Commitments and Licence on the Copernicus Marine website. '''Processing information:''' Processing of these data was mainly carried out in the OC-CCI framework producing a climate-quality consistent dataset using the latest and most complete knowledge of satellite sensor calibration, characterization and attitude, complete (as far as possible) ancillary data sets, latest versions of models and algorithms etc. The data were then repackaged, using custom software, to suit the requirements of CMEMS. '''Description of observation methods/instruments:''' The remote sensing of Ocean Colour represents a measure of the spectral variations in the light leaving the water surface, subsequently interpreted in terms of concentrations of optically-significant constituents in the water. The electromagnetic signal collected by the sensor on-board the satellite is largely determined by photons that have never reached the water surface, but have been backscattered within the atmosphere through multiple interactions between gas molecules and aerosols. After removing the atmospheric contribution, the water leaving radiance recorded at a given time by the satellite reflects the optical properties of the water which, in turn, mirrors a specific structure and biogeochemical composition of the marine waters. '''Quality / Accuracy / Calibration information:''' The user is referred to the QUID 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-00103