Search

Abstract Marine heatwaves (MHWs) represent anomalously warm temperature conditions of seawater that may affect marine life and ocean biogeochemistry. Under such conditions, phytoplankton communities may modify their structure and functions, and their resilience is not assured. This study characterizes the impact of MHWs on the phytoplankton spring bloom in the North‐Western Mediterranean Sea. Here, we synergistically combine autonomous observations from BioGeoChemical‐Argo floats, satellite‐based and marine ecosystem model data, and show that MHW events occurring during winter drastically inhibit phytoplankton carbon biomass in spring by up to 70%. Such reduction is related to the enhanced stratification of the water column under MHWs which hinders the renewal of nutrients from deep‐ocean reservoirs, thus preventing surface phytoplankton from blooming. This process negatively impacts particulate organic carbon stocks within the mixed layer, while severe events cause an earlier shift of phytoplankton phenology that provokes changes in zooplankton biomass distribution.

The response of aquatic and terrestrial organisms to climate change can depend on biological sex. A key challenge is to unravel the interactive effects of sex and climate change at the individual and population levels and the cascading effects on communities. This new understanding is essential to improve climate adaptation and mitigation strategies.

Atmospheric temperature is a key variable to detect and attribute climate change. Due to the relative sparseness of ground-based observations and heterogeneity of satellite data, global atmospheric reanalysis products are considered valuable datasets for studying and monitoring the climate, since these usually ensure spatially complete and continuous temporal coverage. Consequently, evaluating differences among the existing reanalyses is key to identifying inconsistencies. To this aim, the current study intercompares the climatological mean, variability, and linear trends for upper air temperature provided from four recent atmospheric reanalysis products (ERA5, ERA-Interim, MERRA-2, and JRA-55) The Reanalysis Multi-Model Ensemble-mean (RMME) is used as a comparator. Radiosonde observations are included for comparison on the regional scale (tropics). The results reveal that all evaluated reanalyses provide a consistent reproduction of the upper-air temperature profile. Temperature differences from RMME within ±0.25 K are found in both the troposphere and lower stratosphere, except for a few specific regions. Larger differences (>± 1.5 K) and discrepancies among the datasets are found in the upper stratosphere. Agreement between reanalyses increased after 1998. Differences in the temperature time series and seasonal cycle at the regional scales are smaller in the Northern Hemisphere mid-latitudes than in the tropics and Southern Hemisphere. A robust cooling trend in the lower stratosphere during the period of largest ozone depletion (1980-1997) and a warming trend in the troposphere for the period 1998-2018 are clearly found, at the global level, in all the datasets. Temperature trends and variability in the tropics are consistent in all reanalysis datasets with the homogenized radiosonde records from the lower to middle troposphere and in the lower stratosphere. However, large differences are found in the upper troposphere, tropical tropopause layer (TTL), and middle stratosphere. The well-known temperature variability in the lower stratosphere associated with Quasi-Biennial Oscillation (QBO) is captured in both reanalyses and observational datasets. Among the reanalyses, ERA5 shows the smallest temperature difference with respect to homogenized upper-air radiosoundings records.

This paper presents major gaps and challenges for implementing the UN Decade of Ocean Science for Sustainable Development (2021-2030) in the Mediterranean region. The authors make recommendations on the scientific knowledge needs and co-design actions identified during two consultations, part of the Decade preparatory-phase, framing them in the Mediterranean Sea’s unique environmental and socio-economic perspectives. According to the ‘Mediterranean State of the Environment and Development Report 2020’ by the United Nations Environment Programme Mediterranean Action Plan and despite notable progress, the Mediterranean region is not on track to achieve and fully implement the Sustainable Development Goals of Agenda 2030. Key factors are the cumulative effect of multiple human-induced pressures that threaten the ecosystem resources and services in the global change scenario. The basin, identified as a climate change vulnerability hotspot, is exposed to pollution and rising impacts of climate change. This affects mainly the coastal zones, at increasing risk of extreme events and their negative effects of unsustainable management of key economic assets. Transitioning to a sustainable blue economy is the key for the marine environment’s health and the nourishment of future generations. This challenging context, offering the opportunity of enhancing the knowledge to define science-based measures as well as narrowing the gaps between the Northen and Southern shores, calls for a joint (re)action. The paper reviews the state of the art of Mediterranean Sea science knowledge, sets of trends, capacity development needs, specific challenges, and recommendations for each Decade’s societal outcome. In the conclusions, the proposal for a Mediterranean regional programme in the framework of the Ocean Decade is addressed. The core objective relies on integrating and improving the existing ocean-knowledge, Ocean Literacy, and ocean observing capacities building on international cooperation to reach the “Mediterranean Sea that we want”.

The role of eddies in the North Atlantic decadal variability is investigated in this study by using two ocean reanalyses, including an eddy permitting (or eddy poor) reanalysis with horizontal resolution of 0.25 degree and 75 vertical levels and an eddy resolving (or eddy rich) reanalysis with horizontal resolution of 1/12 degree and 50 vertical levels. The prominent mid-1990s warming and post-2005 cooling trend as part of the North Atlantic decadal variability is well displayed in both reanalyses with no significant difference between them. The main driver of the mid-1990s warming and post-2005 cooling trend is the increase and reduction of the meridional ocean heat transport showing similar patterns in both reanalyses. The relative contribution of the heat transport anomalies from eddies to the total heat transport anomalies is slightly larger in eddy resolving than in eddy permitting ocean reanalysis. However, the total mean ocean meridional heat transport increases by 10% in eddy resolving reanalysis with respect to eddy permitting reanalysis and is mainly due to the associated increase of the mean states (temperature and velocity). Therefore, the increase of eddy population due to the increase of horizontal resolution, found by comparing the two datasets, does not affect the MHT anomalies significantly and, consequently, the North Atlantic decadal variability. It is found that the importance of the model horizontal resolution for the North Atlantic decadal variability depends on the interaction between the eddies (small scale) and the mean state (large scale) at decadal time scales. Although the fast increase of computational power will allow soon for eddy-resolving predictions, the need to use high resolution modeling tools for decadal predictions depends on the importance of initialization methods and the interaction between small scale and large-scale variabilities. This study has pivotal implications for the development of North Atlantic decadal prediction systems.

Since May 2022, the Mediterranean Sea has been experiencing an exceptionally long marine heatwave event. Warm anomalies, mainly occurring in the Western basin, have persisted until boreal spring 2023, making this event the longest Mediterranean marine heat wave of the last four decades. In this work, the 2022/2023 anomaly is characterized, using in-situ and satellite measurements, together with state of the art reanalysis products. The role of atmospheric forcing is also investigated; the onset and growth of sea surface temperature anomalies is found to be related to the prevalence of anticyclonic conditions in the atmosphere, which have also caused severe droughts in the Mediterranean region over the same period. Analysis of in-situ observations from the Lampedusa station and of ocean reanalyzes reveals that wind-driven vertical mixing led to the penetration of the warm anomalies below the sea surface, where they have persisted for several months, particularly in the central part of the basin. The evolution of the 2022/23 event is compared with the severe 2003 event, to put recent conditions in the context of climate change.

Sea surface salinity (SSS) is one of the Essential Climate Variables (ECVs), defined by the Global Climate Observing System (GCOS). Salinity is modified by river discharge, land run-off, precipitation, and evaporation, and it is advected by oceanic currents. In turn, ocean circulation, the water cycle, and biogeochemistry are deeply impacted by salinity variations. The Mediterranean Sea represents a hot spot for the variability of salinity. Despite the ever-increasing number of moorings and floating buoys, in situ SSS estimates have low coverage, hindering the monitoring of SSS patterns. Conversely, satellite sensors provide SSS surface data at high spatial and temporal resolution, complementing the sparseness of in situ datasets. Here, we describe a multidimensional optimal interpolation algorithm, specifically configured to provide a new daily SSS dataset at 1/16° grid resolution, covering the entire Mediterranean Sea (Med L4 SSS). The main improvements in this regional algorithm are: the ingestion of satellite SSS estimates from multiple satellite missions (NASA’s Soil Moisture Active Passive (SMAP), ESA’s Soil Moisture and Ocean Salinity (SMOS) satellites), and a new background (first guess), specifically built to improve coastal reconstructions. The multi-sensor Med L4 SSS fields have been validated against independent in situ SSS samples, collected between 2010–2020. They have also been compared with global weekly Copernicus Marine Environment Monitoring Service (CMEMS) and Barcelona Expert Centre (BEC) regional products, showing an improved performance. Power spectral density analyses demonstrated that the Med L4 SSS field achieves the highest effective spatial resolution, among all the datasets analysed. Even if the time series is relatively short, a clear interannual trend is found, leading to a marked salinification, mostly occurring in the Eastern Mediterranean Sea.

The Mediterranean community represented in this paper is the result of more than 30 years of EU and nationally funded coordination, which has led to key contributions in science concepts and operational initiatives. Together with the establishment of operational services, the community has coordinated with universities, research centers, research infrastructures and private companies to implement advanced multi-platform and integrated observing and forecasting systems that facilitate the advancement of operational services, scientific achievements and mission-oriented innovation. Thus, the community can respond to societal challenges and stakeholders needs, developing a variety of fit-for-purpose services such as the Copernicus Marine Service. The combination of state-of-the-art observations and forecasting provides new opportunities for downstream services in response to the needs of the heavily populated Mediterranean coastal areas and to climate change. The challenge over the next decade is to sustain ocean observations within the research community, to monitor the variability at small scales, e.g., the mesoscale/submesoscale, to resolve the sub-basin/seasonal and inter-annual variability in the circulation, and thus establish the decadal variability, understand and correct the model-associated biases and to enhance model-data integration and ensemble forecasting for uncertainty estimation. Better knowledge and understanding of the level of Mediterranean variability will enable a subsequent evaluation of the impacts and mitigation of the effect of human activities and climate change on the biodiversity and the ecosystem, which will support environmental assessments and decisions. Further challenges include extending the science-based added-value products into societal relevant downstream services and engaging with communities to build initiatives that will contribute to the 2030 Agenda and more specifically to SDG14 and the UN's Decade of Ocean Science for sustainable development, by this contributing to bridge the science-policy gap. The Mediterranean observing and forecasting capacity was built on the basis of community best practices in monitoring and modeling, and can serve as a basis for the development of an integrated global ocean observing system.

Spectrally resolved water-leaving radiances (ocean colour) and inferred chlorophyll concentration are key to studying phytoplankton dynamics at seasonal and inter-annual scales, for a better understanding of the role of phytoplankton in marine biogeochemistry; the global carbon cycle; and the response of marine ecosystems to climate variability, change and feedback processes. Ocean colour data also have a critical role in operational observation systems monitoring coastal eutrophication, harmful algal blooms, and sediment plumes. The contiguous ocean-colour record reached 21 years in 2018; however, it is comprised of a number of one-off missions such that creating a consistent time-series of ocean-colour data requires merging of the individual sensors (including MERIS, Aqua-MODIS, SeaWiFS, VIIRS, and OLCI) with differing sensor characteristics, without introducing artefacts. By contrast, the next decade will see consistent observations from operational ocean colour series with sensors of similar design and with a replacement strategy. Also, by 2029 the record will start to be of sufficient duration to discriminate climate change impacts from natural variability, at least in some regions. This paper describes the current status and future prospects in the field of ocean colour focusing on large to medium resolution observations of oceans and coastal seas. It reviews the user requirements in terms of products and uncertainty characteristics and then describes features of current and future satellite ocean-colour sensors, both operational and innovative. The key role of in situ validation and calibration is highlighted as are ground segments that process the data received from the ocean-colour sensors and deliver analysis-ready products to end-users. Example applications of the ocean-colour data are presented, focusing on the climate data record and operational applications including water quality and assimilation into numerical models. Current capacity building and training activities pertinent to ocean colour are described and finally a summary of future perspectives is provided.

Sea surface temperature (SST) is a fundamental physical variable for understanding, quantifying and predicting complex interactions between the ocean and the atmosphere. Such processes determine how heat from the sun is redistributed across the global oceans, directly impacting large- and small-scale weather and climate patterns. The provision of daily maps of global SST for operational systems, climate modeling and the broader scientific community is now a mature and sustained service coordinated by the Group for High Resolution Sea Surface Temperature (GHRSST) and the CEOS SST Virtual Constellation (CEOS SST-VC). Data streams are shared, indexed, processed, quality controlled, analyzed, and documented within a Regional/Global Task Sharing (R/GTS) framework, which is implemented internationally in a distributed manner. Products rely on a combination of low-Earth orbit infrared and microwave satellite imagery, geostationary orbit infrared satellite imagery, and in situ data from moored and drifting buoys, Argo floats, and a suite of independent, fully characterized and traceable in situ measurements for product validation (Fiducial Reference Measurements, FRM). Research and development continues to tackle problems such as instrument calibration, algorithm development, diurnal variability, derivation of high-quality skin and depth temperatures, and areas of specific interest such as the high latitudes and coastal areas. In this white paper, we review progress versus the challenges we set out 10 years ago in a previous paper, highlight remaining and new research and development challenges for the next 10 years (such as the need for sustained continuity of passive microwave SST using a 6.9 GHz channel), and conclude with needs to achieve an integrated global high-resolution SST observing system, with focus on satellite observations exploited in conjunction with in situ SSTs. The paper directly relates to the theme of Data Information Systems and also contributes to Ocean Observing Governance and Ocean Technology and Networks within the OceanObs2019 objectives. Applications of SST contribute to all the seven societal benefits, covering Discovery; Ecosystem Health & Biodiversity; Climate Variability & Change; Water, Food, & Energy Security; Pollution & Human Health; Hazards and Maritime Safety; and the Blue Economy.

The Copernicus Marine Environment Monitoring Service (CMEMS) provides regular and systematic reference information on the physical and biogeochemical ocean and sea-ice state for the global ocean and the European regional seas. CMEMS serves a wide range of users (more than 15,000 users are now registered to the service) and applications. Observations are a fundamental pillar of the CMEMS value-added chain that goes from observation to information and users. Observations are used by CMEMS Thematic Assembly Centres (TACs) to derive high-level data products and by CMEMS Monitoring and Forecasting Centres (MFCs) to validate and constrain their global and regional ocean analysis and forecasting systems. This paper presents an overview of CMEMS, its evolution, and how the value of in situ and satellite observations is increased through the generation of high-level products ready to be used by downstream applications and services. The complementary nature of satellite and in situ observations is highlighted. Long-term perspectives for the development of CMEMS are described and implications for the evolution of the in situ and satellite observing systems are outlined. Results from Observing System Evaluations (OSEs) and Observing System Simulation Experiments (OSSEs) illustrate the high dependencies of CMEMS systems on observations. Finally future CMEMS requirements for both satellite and in situ observations are detailed.

Measuring the ocean surface currents at high spatio-temporal resolutions is crucial for scientific and socio-economic applications. Since the early 1990s, the synoptic and global-scale monitoring of the ocean surface currents has been provided by constellations of radar altimeters. By construction, altimeter constellations provide only the geostrophic component of the marine surface currents. In addition, given the effective spatial-temporal resolution of the altimeter-derived products (O (100 km) and O (10 days), respectively), only the largest ocean mesoscale features can be resolved. In order to enhance the altimeter system capabilities, we propose a synergistic use of high resolution sea surface Chlorophyll observations (Chl) and altimeter-derived currents’ estimates. The study is focused on the Mediterranean Sea, where the most energetic signals are found at spatio-temporal scales up to 10 km and a few days. The proposed method allows for inferring the marine surface currents from the evolution of the Chl field, relying on altimeter-derived currents as a first-guess estimate. The feasibility of this approach is tested through an Observing System Simulation Experiment, starting from biogeochemical model outputs distributed by the European Copernicus Marine Service. Statistical analyses based on the 2017 daily data showed that our approach can improve the altimeter-derived currents accuracy up to 50%, also enhancing their effective spatial resolution up to 30 km. Moreover, the retrieved currents exhibit larger temporal variability than the altimeter estimates over annual to weekly timescales. Our method is mainly limited to areas/time periods where/when Chl gradients are larger and are modulated by the marine currents’ advection. Its application is thus more efficient when the surface Chl evolution is not dominated by the biological activity, mostly occurring in the mid-February to mid-March time window in the Mediterranean Sea. Preliminary tests on the method applicability to satellite-derived data are also presented and discussed.

Air–sea heat fluxes are essential climate variables, required for understanding air–sea interactions, local, regional and global climate, the hydrological cycle and atmospheric and oceanic circulation. In situ measurements of fluxes over the ocean are sparse and model reanalysis and satellite data can provide estimates at different scales. The accuracy of such estimates is therefore essential to obtain a reliable description of the occurring phenomena and changes. In this work, air–sea radiative fluxes derived from the SEVIRI sensor onboard the MSG satellite and from ERA5 reanalysis have been compared to direct high quality measurements performed over a complete annual cycle at the ENEA oceanographic observatory, near the island of Lampedusa in the Central Mediterranean Sea. Our analysis reveals that satellite derived products overestimate in situ direct observations of the downwelling short-wave (bias of 6.1 W/m2) and longwave (bias of 6.6 W/m2) irradiances. ERA5 reanalysis data show a negligible positive bias (+1.0 W/m2) for the shortwave irradiance and a large negative bias (−17 W/m2) for the longwave irradiance with respect to in situ observations. ERA5 meteorological variables, which are needed to calculate the air–sea heat flux using bulk formulae, have been compared with in situ measurements made at the oceanographic observatory. The two meteorological datasets show a very good agreement, with some underestimate of the wind speed by ERA5 for high wind conditions. We investigated the impact of different determinations of heat fluxes on the near surface sea temperature (1 m depth), as determined by calculations with a one-dimensional numerical model, the General Ocean Turbulence Model (GOTM). The sensitivity of the model to the different forcing was measured in terms of differences with respect to in situ temperature measurements made during the period under investigation. All simulations reproduced the true seasonal cycle and all high frequency variabilities. The best results on the overall seasonal cycle were obtained when using meteorological variables in the bulk formulae formulations used by the model itself. The derived overall annual net heat flux values were between +1.6 and 40.4 W/m2, depending on the used dataset. The large variability obtained with different datasets suggests that current determinations of the heat flux components and, in particular, of the longwave irradiance, need to be improved. The ENEA oceanographic observatory provides a complete, long-term, high resolution time series of high quality in situ observations. In the future, more similar sites worldwide will be needed for model and satellite validations and to improve the determination of the air–sea exchange and the understanding of related processes.

Remote sensing data provide a huge number of sea surface observations, but cannot give direct information on deeper ocean layers, which can only be provided by sparse in situ data. The combination of measurements collected by satellite and in situ sensors represents one of the most effective strategies to improve our knowledge of the interior structure of the ocean ecosystems. In this work, we describe a Multi-Layer-Perceptron (MLP) network designed to reconstruct the 3D fields of ocean temperature and chlorophyll-a concentration, two variables of primary importance for many upper-ocean bio-physical processes. Artificial neural networks can efficiently model eventual non-linear relationships among input variables, and the choice of the predictors is thus crucial to build an accurate model. Here, concurrent temperature and chlorophyll-a in situ profiles and several different combinations of satellite-derived surface predictors are used to identify the optimal model configuration, focusing on the Mediterranean Sea. The lowest errors are obtained when taking in input surface chlorophyll-a, temperature, and altimeter-derived absolute dynamic topography and surface geostrophic velocity components. Network training and test validations give comparable results, significantly improving with respect to Mediterranean climatological data (MEDATLAS). 3D fields are then also reconstructed from full basin 2D satellite monthly climatologies (1998–2015) and resulting 3D seasonal patterns are analyzed. The method accurately infers the vertical shape of temperature and chlorophyll-a profiles and their spatial and temporal variability. It thus represents an effective tool to overcome the in-situ data sparseness and the limits of satellite observations, also potentially suitable for the initialization and validation of bio-geophysical models.

Measurements of ocean surface topography collected by satellite altimeters provide geostrophic estimates of the sea surface currents at relatively low resolution. The effective spatial and temporal resolution of these velocity estimates can be improved by optimally combining altimeter data with sequences of high resolution interpolated (Level 4) Sea Surface Temperature (SST) data, improving upon present-day values of approximately 100 km and 15 days at mid-latitudes. However, the combined altimeter/SST currents accuracy depends on the area and input SST data considered. Here, we present a comparative study based on three satellite-derived daily SST products: the Remote Sensing Systems (REMSS, 1/10 ∘ resolution), the UK Met Office OSTIA (1/20 ∘ resolution), and the Multiscale Ultra-High resolution SST (1/100 ∘ resolution). The accuracy of the marine currents computed with our synergistic approach is assessed by comparisons with in-situ estimated currents derived from a global network of drifting buoys. Using REMSS SST, the meridional currents improve up to more than 20% compared to simple altimeter estimates. The maximum global improvements for the zonal currents are obtained using OSTIA SST, and reach 6%. Using the OSTIA SST also results in slight improvements (≃1.3%) in the zonal flow estimated in the Southern Ocean (45 ∘ S to 70 ∘ S). The homogeneity of the input SST effective spatial resolution is identified as a crucial requirement for an accurate surface current reconstruction. In our analyses, this condition was best satisfied by the lower resolution SST products considered.

In the coastal environment the optical properties can vary on temporal scales that are shorter than the near-polar orbiting satellite temporal resolution (~1 image per day), which does not allow capturing most of the coastal optical variability. The objective of this work is to fill the gap between the near-polar orbiting and geostationary sensor temporal resolutions, as the latter sensors provide multiple images of the same basin during the same day. To do that, a Level 3 hyper-temporal analysis-ready Ocean Color (OC) dataset, named Virtual Geostationary Ocean Color Sensor (VGOCS), has been created. This dataset contains the observations acquired over the North Adriatic Sea by the currently functioning near-polar orbiting sensors, allowing approaching the geostationary sensor temporal resolution. The problem in using data from different sensors is that they are characterized by different uncertainty sources that can introduce artifacts between different satellite images. Hence, the sensors have different spatial and spectral resolutions, their calibration procedures can have different accuracies, and their Level 2 data can be retrieved using different processing chains. Such differences were reduced here by adjusting the satellite data with a multi-linear regression algorithm that exploits the Fiducial Reference Measurements data stream of the AERONET-OC water-leaving radiance acquired at the Acqua Alta Oceanographic Tower, located in the Gulf of Venice. This work aims to prove the suitability of VGOCS in analyzing the coastal optical variability, presenting the improvement brought by the adjustment on the quality of the satellite data, the VGOCS spatial and temporal coverage, and the inter-sensor differences. Hence, the adjustment will strongly increase the agreement between the satellite and in situ data and between data from different near-polar orbiting OC imagers; moreover, the adjustment will make available data traditionally masked in the standard processing chains, increasing the VGOCS spatial and temporal coverage, fundamental to analyze the coastal optical variability. Finally, the fulfillment by VGOCS of the three conditions for a hyper-temporal dataset will be demonstrated in this work.

In the context of the Copernicus Program, EUMETSAT prioritizes the creation of an ocean color infrastructure for system vicarious calibration (OC-SVC). This work aims to reply to this need by proposing the European Radiometry Buoy and Infrastructure (EURYBIA). EURYBIA is designed as an autonomous European infrastructure operating within the Marine Optical Network (MarONet) established by University of Miami (Miami, FL, USA) based on the Marine Optical Buoy (MOBY) experience and NASA support. MarONet addresses SVC requirements in different sites, consistently and in a traceable way. The selected EURYBIA installation is close to the Lampedusa Island in the central Mediterranean Sea. This area is widely studied and hosts an Atmospheric and Oceanographic Observatory for long-term climate monitoring. The EURYBIA field segment comprises off-shore and on-shore infrastructures to manage the observation system and perform routine sensors calibrations. The ground segment includes the telemetry center for data communication and the processing center to compute data products and uncertainty budgets. The study shows that the overall uncertainty of EURYBIA SVC gains computed for the Sentinel-3 OLCI mission under EUMETSAT protocols is of about 0.05% in the blue-green wavelengths after a decade of measurements, similar to that of the reference site in Hawaii and in compliance with requirements for climate studies.

Estimating long-term modifications of the sea surface temperature (SST) is crucial for evaluating the current state of the oceans and to correctly assess the impact of climate change at regional scales. In this work, we analyze SST variations within the Mediterranean Sea and the adjacent Northeastern Atlantic box (west of the Strait of Gibraltar) over the last 37 years, by using a satellite-based dataset from the Copernicus Marine Environment Monitoring Service (CMEMS). We found a mean warming trend of 0.041 ± 0.006 ∘ C/year over the whole Mediterranean Sea from 1982 to 2018. The trend has an uneven spatial pattern, with values increasing from 0.036 ± 0.006 ∘ C/year in the western basin to 0.048 ± 0.006 ∘ C/year in the Levantine−Aegean basin. The Northeastern Atlantic box and the Mediterranean show a similar trend until the late 1990s. Afterwards, the Mediterranean SST continues to increase, whereas the Northeastern Atlantic box shows no significant trend, until ~2015. The observed change in the Mediterranean Sea affects not only the mean trend but also the amplitude of the Mediterranean seasonal signal, with consistent relative increase and decrease of summer and winter mean values, respectively, over the period considered. The analysis of SST changes occurred during the “satellite era” is further complemented by reconstructions also based on direct in situ SST measurements, i.e., the Extended Reconstructed SST (ERSST) and the Hadley Centre Sea Ice and Sea Surface Temperature dataset (HadISST), which go back to the 19th century. The analysis of these longer time series, covering the last 165 years, indicates that the increasing Mediterranean trend, observed during the CMEMS operational period, is consistent with the Atlantic Multidecadal Oscillation (AMO), as it closely follows the last increasing period of AMO. This coincidence occurs at least until 2007, when the apparent onset of the decreasing phase of AMO is not seen in the Mediterranean SST evolution.

Monitoring sea surface salinity (SSS) and density variations is crucial to investigate the global water cycle and the ocean dynamics, and to analyse how they are impacted by climate change. Historically, ocean salinity and density have suffered a poor observational coverage, which hindered an accurate assessment of their surface patterns, as well as of associated space and time variability and trends. Different approaches have thus been proposed to extend the information obtained from sparse in situ measurements and provide gap-free fields at regular spatial and temporal resolution, based on the combination of in situ and satellite data. In the framework of the Copernicus Marine Environment Monitoring Service, a daily (weekly sampled) global reprocessed dataset at ¼° × ¼° resolution has been produced by modifying a multivariate optimal interpolation (OI) technique originally developed within MyOcean project. The algorithm has been applied to in situ salinity/density measurements covering the period from 1993 to 2016, using satellite sea surface temperature differences to constrain the surface patterns. This improved algorithm and the new dataset are described and validated here with holdout approach and independent data.

We present a methodology to detect oil spills using MODIS near-infrared sun glittered radiance imagery. The methodology was developed by using a set of seven MODIS images (training dataset) and validated using four other images (validation dataset). The method is based on the ratio image R = L'GN/LGN, where L'GN is the MODIS-retrieved normalized sun glint radiance image and LGN the same quantity, but obtained from the Cox and Munk isotropic (independent of wind direction) sun glint model. We show that in the R image, while clean water pixel values tend to one, oil spills stand out as anomalies. Moreover, we provide a criterion to distinguish between positive and negative oil-water contrast. A pixel in an R image is classified as a potential oil spill or water via a variable threshold Rs as a function of L'GN, where the threshold values are obtained from the slicks of our training dataset. Two different fitting curves are provided for Rs, according to the contrast sign. The selection of the correct fitting curve is based on the contrast type, resulting from the criterion above. Results indicate that the thresholding is able to isolate the spills and that the spills of the validation dataset are successfully detected. Spurious look-alike features, such as clouds, and other non-spill features, e.g., large areas at the glint region border, are also detected as oil spills and must be eliminated. We believe that our methodology represents a novel and promising, though preliminary, approach towards automatic oil spill detection in optical satellite images.

Regional relationships to estimate the main Phytoplankton Functional Types (PFTs) and Size Classes (PSCs) from satellite data are presented. Following the abundance-based approach and selecting the Total Chlorophyll a (TChla) as descriptor of the trophic status of the environment, empirical relations between the TChla concentration and seven accessory pigments, marker for the main algal groups, have been developed for the Mediterranean Sea. Using only in-situ data acquired in this basin, firstly a previous regional diagnostic pigment analysis has been conducted to evaluate the specific pigment ratios featuring the phytoplankton assemblage that occurs in the Mediterranean Sea. Secondly, the new regional PFT and PSC algorithms have been calibrated and validated on the in-situ dataset. The statistical analysis showed a very good predictive power for all the new regional models. A quantitative comparison with global abundance-based models applied to our validation dataset showed that the regionalization improves the uncertainty and the spread of about one order of magnitude for all the classes (e.g., in the nano class, where the mean bias error improves from −0.056 to 0.001 mg m−3). These results highlighted that a regionalization for the PSC and PFT estimates are required, to take into account the peculiar bio-optical properties of the Mediterranean Sea. Finally, the new regional equations have been applied to the Mediterranean TChla satellite (1998–2015) time series to estimate annual and monthly PFT and PSC climatology. The analysis of the climatological maps, relative to the phytoplankton assemblage distribution patterns, reveals that all the three size classes reach their maxima in the higher nutrient areas, with absolute values >3 mg m−3 of TChla for micro-, and about 1.6 and 0.4 mg m−3 for nano- and pico-phytoplankton, respectively. Moreover, the nano component shows intermediate percentage values in the whole basin, ranging from 30 to 40% of the TChla in the western basin, up to 45% in the more productive areas. In terms of chlorophyll concentration, in the coastal areas we find the predominance of the Diatoms and Haptophytes, while in the ultra-oligotrophic waters Prokaryotes predominates on the other groups, constituting the principal component of the pico-phytoplankton.

Sediment dynamics is the primary driver of the evolution of the coastal geomorphology and of the underwater shelf clinoforms. In this paper, we focus on mesoscale and sub-mesoscale processes, such as coastal currents and river plumes, and how they shape the sediment dynamics at regional or basin spatial scales. A new methodology is developed that combines observational data with numerical modelling: the aim is to pair satellite measurements of suspended sediment with velocity fields from numerical oceanographic models, to obtain an estimation of the sediment flux. A numerical divergence of this flux is then computed. The divergence field thus obtained shows how the aforementioned mesoscale processes distribute the sediments. The approach was applied and discussed on the Adriatic Sea, for the winter of 2012, using data provided by the ESA Coastcolour project and the output of a run of the MIT General Circulation Model.

We present a study on the potential of the Copernicus Imaging Microwave Radiometer (CIMR) mission for the global monitoring of Sea-Surface Salinity (SSS) using Level-4 (gap-free) analysis processing. Space-based SSS are currently provided by the Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) satellites. However, there are no planned missions to guarantee continuity in the remote SSS measurements for the near future. The CIMR mission is in a preparatory phase with an expected launch in 2026. CIMR is focused on the provision of global coverage, high resolution sea-surface temperature (SST), SSS and sea-ice concentration observations. In this paper, we evaluate the mission impact within the Copernicus Marine Environment Monitoring Service (CMEMS) SSS processing chain. The CMEMS SSS operational products are based on a combination of in situ and satellite (SMOS) SSS and high-resolution SST information through a multivariate optimal interpolation. We demonstrate the potential of CIMR within the CMEMS SSS operational production after the SMOS era. For this purpose, we implemented an Observing System Simulation Experiment (OSSE) based on the CMEMS MERCATOR global operational model. The MERCATOR SSSs were used to generate synthetic in situ and CIMR SSS and, at the same time, they provided a reference gap-free SSS field. Using the optimal interpolation algorithm, we demonstrated that the combined use of in situ and CIMR observations improves the global SSS retrieval compared to a processing where only in situ observations are ingested. The improvements are observed in the 60% and 70% of the global ocean surface for the reconstruction of the SSS and of the SSS spatial gradients, respectively. Moreover, the study highlights the CIMR-based salinity patterns are more accurate both in the open ocean and in coastal areas. We conclude that CIMR can guarantee continuity for accurate monitoring of the ocean surface salinity from space.

We present a method for the remote retrieval of the sea surface currents in the Mediterranean Sea. Combining the altimeter-derived currents with sea-surface temperature information, we created daily, gap-free high resolution maps of sea surface currents for the period 2012–2016. The quality of the new multi-sensor currents has been assessed through comparisons to other surface-currents estimates, as the ones obtained from drifting buoys trajectories (at the basin scale), or HF-Radar platforms and ocean numerical model outputs in the Malta–Sicily Channel. The study yielded that our synergetic approach can improve the present-day derivation of the surface currents in the Mediterranean area up to 30% locally, with better performances for the the meridional component of the motion and in the western section of the basin. The proposed reconstruction method also showed satisfying performances in the retrieval of the ageostrophic circulation in the Sicily Channel. In this area, assuming the High Frequency Radar-derived currents as reference, the merged multi-sensor currents exhibited improvements with respect to the altimeter estimates and numerical model outputs, mainly due to their enhanced spatial and temporal resolution.

In being at the base of the marine food web, phytoplankton is particularly important for marine ecosystem functioning (e.g., biodiversity). Strong anthropization, over-exploitation of natural resources, and climate change affect the natural amount of phytoplankton and, therefore, represent a continuous threat to the biodiversity in marine waters. In particular, a concerning risks for coastal waters is the increase in nutrient inputs of terrestrial/anthropogenic origin that can lead to undesirable modifications of phytoplankton concentration (i.e., eutrophication). Monitoring chlorophyll (Chl) concentration, which is a proxy of phytoplankton biomass, is an efficient tool for recording and understanding the response of the marine ecosystem to human pressures and thus for detecting eutrophication. Here, we compute Chl trends over the Mediterranean Sea by using satellite data, also highlighting the fact that remote sensing may represent an efficient and reliable solution to synoptically control the "good environmental status" (i.e., the Marine Directive to achieve Good Environmental Status of EU marine waters by 2020) and to assess the application of international regulations and environmental directives. Our methodology includes the use of an ad hoc regional (i.e., Mediterranean) algorithm for Chl concentration retrieval, also accounting for the difference between offshore (i.e., Case I) and coastal (i.e., Case II) waters. We apply the Mann-Kendall test and the Sens's method for trend estimation to the Chl concentration de-seasonalized monthly time series, as obtained from the X-11 technique. We also provide a preliminary analysis of some particular trends by evaluating their associated inter-annual variability. The high spatial resolution of our approach allows a clear identification of intense trends in those coastal waters that are affected by river outflows. We do not attempt to attribute the observed trends to specific anthropogenic events. However, the trends that we document are consistent with the findings of several previous studies.

Knowledge of the vertical structure of the bio-chemical properties of the ocean is crucial for the estimation of primary production, phytoplankton distribution, and biological modelling. The vertical profiles of chlorophyll-a (Chla) are available via in situ measurements that are usually quite rare and not uniformly distributed in space and time. Therefore, obtaining estimates of the vertical profile of the Chla field from surface observations is a new challenge. In this study, we employed an Artificial Neural Network (ANN) to reconstruct the 3-Dimensional (3D) Chla field in the Mediterranean Sea from surface satellite estimates. This technique is able to reproduce the highly nonlinear nature of the relationship between different input variables. A large in situ dataset of temperature and Chla calibrated fluorescence profiles, covering almost all Mediterranean Sea seasonal conditions, was used for the training and test of the network. To separate sources of errors due to surface Chla and temperature satellite estimates, from errors due to the ANN itself, the method was first applied using in situ surface data and then using satellite data. In both cases, the validation against in situ observations shows comparable statistical results with respect to the training, highlighting the feasibility of applying an ANN to infer the vertical Chla field from surface in situ and satellite estimates. We also analyzed the usefulness of our approach to resolve the Chla prediction at small temporal scales (e.g., day) by comparing it with the most widely used Mediterranean climatology (MEDATLAS). The results demonstrated that, generally, our method is able to reproduce the most reliable profile of Chla from synoptical satellite observations, thus resolving finer spatial and temporal scales with respect to climatology, which can be crucial for several marine applications. We demonstrated that our 3D reconstructed Chla field could represent a valid alternative to overcome the absence or discontinuity of in situ sampling.

Knowledge of the link between ocean hydrodynamics and distribution of small pelagic fish species is fundamental for the sustainable management of fishery resources. Both commercial and scientific communities are indeed seeking to provide services that could "connect the dots" among in situ and remote observations, numerical ocean modelling, and fisheries. In the Mediterranean Sea and, in particular, in the Sicily Channel the reproductive strategy of the European Anchovy (Engraulis encrasicolus) is strongly influenced by the oceanographic patterns, which are often visible in sea surface temperature satellite data. Based on these experimental evidences, we propose here a more general approach where the role of ocean currents, wind effects, and mesoscale activity are tied together. To investigate how these features affect anchovy larvae distribution, we pair ichthyoplankton observations to a wide remote sensing data set, and to Lagrangian numerical simulations for larval transport. Our analysis shows that while the wind-induced coastal current is able to transport anchovy larvae from spawning areas to the recruiting area off the Sicilian south-eastern tip, significant cross-shore transport due to the combination of strong northwesterly mistral winds and topographic effects delivers larvae away from the coastal conveyor belt. We then use a potential vorticity approach to describe the occurrence of larvae cross-shore transport. We conclude that monitoring and quantifying the upwelling on the southern Sicilian coast during the spawning season allows to estimate the cross-shore transport of larvae and the consequent decrease of individuals within the recruiting area.