Abstract Ocean water clarity, influenced by marine chlorophyll concentration, significantly alters the distribution of shortwave radiation in the water column. This work aims to assess the effects of varying chlorophyll on the upper-ocean physical properties and their subsequent impact on the atmosphere, using a coupled ocean-atmosphere regional model for the Mediterranean and Black Seas. We performed 11-year (2011–2021) twin-simulation experiments based on different chlorophyll concentrations to estimate the penetration of solar radiation in the ocean. The first simulation used a monthly climatology field of chlorophyll concentrations derived from satellite observations, while in the second experiment, the chlorophyll concentration was kept constant at 0.05 \$\mathrm{m}\mathrm{g}\ {\mathrm{m}}^{-3}\$, representing clear water conditions. Results show that radiative heating driven by chlorophyll amplifies the seasonal cycle of temperature in the upper layers, leading to increased surface warming in summer and surface cooling in winter. Also, higher surface chlorophyll contributes to cooling in subsurface layers throughout the year due to its shading effect. The temperature response to chlorophyll variations is controlled by the mixed layer depth and a balance between (a) direct near-surface radiative heating due to the chlorophyll absorption and (b) indirect cooling resulting from vertical turbulent mixing processes with subsurface waters. The atmosphere moderates the seasonal sea surface temperature (SST) response caused by chlorophyll differential heating primarily through changes in latent heat flux. Ultimately, our simulations suggest that increased surface chlorophyll concentrations enhance the Mediterranean overturning circulation, highlighting the necessity of incorporating realistic optical forcing into regional climate modeling studies.

The contribution of surface and lateral forcing to the observed Arabian Gulf warming trends is studied based on the results of a high-resolution (1/100°, 60 vertical layers) MIT general circulation model (MITgcm) covering the period 1993-2021. The model validation against available observations reveals that the simulation satisfactorily reproduces the main features of the Arabian Gulf's dynamics and their variability. We show that the heat content of the Arabian Gulf generally follows the reported variability of sea surface temperature, with significant increasing trends of 0.1 × 10 7 J m −3 and 0.2°C per decade. The interannual variability of the heat content is dominated by the surface heat fluxes, while the long-term warming of the basin is primarily driven by lateral fluxes. The analyses of the heat exchanges through the Strait of Hormuz indicate a pronounced upward trend in the transported heat toward the Arabian Gulf, which is associated with an increase in both the volume and temperature of the exchanged waters. Considering the inflow and outflow in the Strait separately, the temperature increase is more prominent in the inflowing waters; however, the dominant factor driving the rising trend in heat content exchanges is the increase in the volume of waters being exchanged. This implies that the observed warming of the Arabian Gulf during the investigated period is directly related to the acceleration of its overturning circulation.

Accurate modelling of air-sea processes is essential for reliable forecasts of Mediterranean tropical-like cyclones (also known as “Medicanes”). Medicanes occasionally develop in the Mediterranean causing extreme weather conditions with catastrophic potential due to excessive precipitation, windstorms, and coastal flooding. In this work, we investigate how the complexity of ocean-wave-atmosphere coupling and model initialization affect the simulated track and intensity of the Medicane Ianos (2020). Results indicate that the model's initial conditions and the cyclone's development stage are the main drivers of track position errors, while ocean and wave feedback have a significant impact on the intensity and evolution of the cyclone. Compared with an atmosphere-only simulation, an atmosphere-ocean coupled system reproduces the cyclone's SST cooling effect (up to 3.7 °C), in agreement also with the satellite observations thus, reducing the cyclone intensity, as estimated by the minimum MSLP, the 10-m wind speed and the surface enthalpy flux. Adding a wave model to the coupled system, further increases the magnitude of ocean cooling (by about 1.2 °C), due to increased sea surface roughness leading to increased wind stress and enhanced upper ocean mixing. Overall, surface waves are shown to have competing effects on cyclone intensity i.e., negative feedback via increasing the surface momentum flux and positive feedback via increasing the enthalpy flux, the latter being more sensitive to surface roughness rather than to SST modifications brought by the wave coupled system. The turbulent air-sea fluxes under high winds, appear to be very sensitive to sea-state patterns resolved by the coupled models, highlighting the need to improve forecasting systems for extreme weather events in the Mediterranean.

A CTD/ADCP/surface-drifter survey in fall 2004 reveals the behaviour of a mesoscale unidirectional flow coming from the Cretan sea in the south with depths ∼1000 m and entering a channel-like area of the Cycladic shelf in the north, that forms a contraction which leads to a bottom elevation (sill depth ∼100 m), and finally returning into the Cretan Sea in the lee-side of the sill. The flow decelerates/accelerates upstream/downstream of the sill. The along-stream density contours near the sill bottom are raised prior to reaching the sill, while they deepen in the lee side of it indicating supercriticality. The long-wavelength internal wave speeds with realistic stratification and no-rotation are higher than the section averaged flow speeds and indicate subcriticality. A key element in this apparent paradox is the large height of the sill that potentially increases the body (drag) force exerted on the flow by the sill while flow blocking is also observed upstream of the sill.

This study examines the impact of different turbidity products on the Aegean Sea surface physical characteristics, by performing twin-experiment simulations using a high-resolution regional ocean model. The turbidity products used include an in-situ based diffuse attenuation coefficient dataset at 490 nm (kd490, in m- 1) and a satellite derived kd490 product. Satellite turbidity products are broadly used in ocean simulations due to their spatiotemporal coverage and algorithm universality. Their validation and empirical components are trained mainly in phytoplankton driven regions and this may cause systematic differences in oligotrophic areas of variable optical properties’ composition. In the Aegean Sea, the in-situ based turbidity product accounts for the contribution of suspended particles in the solar heating profile, having further implications in the surface characteristics. The Aegean Sea upper-ocean thermohaline characteristics and general circulation patterns, reveal distinct differences between the twin-experiment simulations, showcasing mesoscale to locally induced impact of the turbidity variations. The turbidity impact on the air-sea interaction fluxes affects both thermodynamic processes i.e., solar radiation penetration and absorption in the water column, as well as dynamic processes i.e., momentum fluxes due to changes of the sea surface temperature and subsequently to the momentum drag coefficient. The Aegean Sea surface characteristics in the in-situ based turbidity product simulation, show a stronger decoupling between the North and the South Aegean Sea, when compared with the satellite derived turbidity product simulation. These results highlight the importance of incorporating more realistic turbidity products in ocean models, especially for optically complex regions such as the Aegean Sea.

The dynamics of fluid flows give rise to robust, persistent circulation features that underpin the flow and exert strong control over the advection of water masses, either enhancing it or suppressing it, collectively known as lagrangian coherent structures. Lagrangian approaches and metrics have been shown to be better suited than eulerian ones at locating and delineating such structures and capturing the effect they have on the formation and dispersion of water masses, particularly at the smaller scales. In this paper, we use the framework of lagrangian coherent structures to analyse the ocean velocity fields over a climatological year obtained from a high-resolution eddy-resolving model in order to investigate the lagrangian regimes that affect the motion, separation and mixing of water masses in the Mediterranean Sea. The lagrangian regimes that develop in each sub-basin over the course of the year are characterised and regions of persistent lagrangian activity and coherent structure formation and presence are identified. A quantitative picture of the seasonal variability of the lagrangian coherent structure-induced horizontal mixing and vortex formation is obtained.

The diffuse attenuation coefficient is an indicator of light availability in the surface layer, and is used in a broad range of applications, including numerical simulations, for the parameterization of the light transmission in the water column. In this study, a new dataset of the diffuse attenuation coefficient for the Eastern Mediterranean Sea test case is developed using an existing optical dataset of 2614 beam attenuation coefficient profiles. This method introduces a way of overcoming the difficulty of measuring the diffuse attenuation coefficient in-situ by utilizing the most routinely measured variable, the beam attenuation coefficient. The proposed approach uses existing semi-analytical relationships and a neural network. The neural network, a multi-layer perceptron regression model, is trained and validated with a dataset of 29398 concurrent bio-optical in-situ measurements from the PROSOPE cruise and remotely sensed surface variables. The model is applied to the Eastern Mediterranean dataset and the results are interpolated into a gridded gap-free field, with a grid resolution of 0.0416° x 0.0416°, which is assessed and compared with a satellite-derived product, investigating their significant differences. The resulting field's mean value is slightly reduced with respect to the satellite product, showing regions of higher turbidity, with the most prominent located in the northern Aegean Sea in regions of excess colored dissolved organic matter and around mesoscale features and in the Cretan and Levantine Sea in regions of higher mesoscale activity.

Summary Global warming affects the aquatic ecosystems, accelerating pathogenic microorganisms' and toxic microalgae's growth and spread in marine habitats, and in bivalve molluscs. New parasite invasions are directly linked to oceanic warming. Consumption of pathogen-infected molluscs impacts human health at different rates, depending, inter alia, on the bacteria taxa. It is therefore necessary to monitor microbiological and chemical contamination of food. Many global cases of poisoning from bivalve consumption can be traced back to Mediterranean regions. This article aims to examine the marine bivalve's infestation rate within the scope of climate change, as well as to evaluate the risk posed by climate change to bivalve welfare and public health. Biological and climatic data literature review was performed from international scientific sources, Greek authorities and State organizations. Focusing on Greek aquaculture and bivalve fisheries, high-risk index pathogenic parasites and microalgae were observed during summer months, particularly in Thermaikos Gulf. Considering the climate models that predict further temperature increases, it seems that marine organisms will be subjected in the long term to higher temperatures. Due to the positive linkage between temperature and microbial load, the marine areas most affected by this phenomenon are characterized as ‘high risk’ for consumer health.

Research infrastructures have been established throughout Europe in order to create robust organizations that will facilitate and enhance research and innovation processes and will advance society with innovative products and services. The Hellenic Integrated Marine Observing, Forecasting and Technology System (component of HIMIOFoTS RI) has been implemented in the framework of the National Roadmap for Research Infrastructures to form a large-scale infrastructure for the marine environment in Greece. It links together ocean observing and forecasting systems, coastal zone monitoring and management practices, as well as ocean engineering testing facilities. The overarching framework of the system supports the coordination of five organizations with expertise in the field of marine science and technology, the central management of research activities, and the common development of services and products. It comprises facilities and resources while it provides open access to research communities (academia, industry) to support the scientific advancements and innovation in their fields. The Hellenic Marine Observing, Forecasting and Technology System was further enhanced during its implementation through significant upgrades and developments in order to extend its observing capacity and the forecasting and technological abilities, while advancing the provided services and products.

The interannual variability of the Mediterranean overturning circulation is investigated using a high-resolution (1/36°) ocean model. As the overturning circulation regulates the replenishment and ventilation of the deep layers, we study the spatiotemporal scales of the maximum value of the overturning streamfunction over three main sub-basins of dense water formation (Aegean Sea, Adriatic, and the northwestern Mediterranean). The variability of the zonal overturning is also discussed. The spectrum analysis shows that the overturning variability has its largest signal on annual timescales in all sub-basins, explained by perpetual winter formation. On shorter frequencies (decadal) there are marked differences observed, due to regional processes of the overturning cells, led by buoyancy flux long-term variability in each sub-basin. The decomposition of the total overturning circulation into barotropic, geostrophic shear, and Ekman components revealed weakening and strengthening for the Aegean and Adriatic Sea total overturning, respectively, with opposite trends for the barotropic and geostrophic shear components. The simultaneous contribution of the Ekman and geostrophic component to the total overturning differentiates the variability of zonal overturning circulation from the local meridional overturning circulation of the three sub-basins. The cross spectra between the maximum overturning value and the buoyancy fluxes also revealed that the system keeps the “memory” of this forcing and shows annual variability.

In this article, Part 1 of a two-part series, we run and evaluate the skill of a regional physical–biogeochemical stochastic ocean model based on NEMO. The domain covers the Bay of Biscay at 1/36° resolution, as a case study for open-ocean and coastal shelf dynamics. We generate model ensembles based on assumptions about errors in the atmospheric forcing, the ocean model parameterizations and in the sources and sinks of the biogeochemical variables. The resulting errors are found to be mainly driven by the wind forcing uncertainties, with the rest of the perturbed forcing and parameters locally influencing the ensemble spread. Biogeochemical uncertainties arise from intrinsic ecosystem model errors and from errors in the physical state. Uncertainties in physical forcing and parameterization are found to have a larger impact on chlorophyll spread than uncertainties in ecosystem sources and sinks. The ensembles undergo quantitative verification with respect to observations, focusing on upper-ocean properties. Despite a tendency for ensembles to be generally under-dispersive, they appear to be reasonably consistent with respect to sea surface temperature data. The largest statistical sea-level biases are observed in coastal regions. These biases hint at the presence of high-frequency error sources currently unaccounted for, and suggest that the ensemble-based uncertainties are unfit to model error covariances for assimilation. Model ensembles for chlorophyll appear to be consistent with ocean colour data only at times. The stochastic model is qualitatively evaluated by analysing its ability at generating consistent multivariate incremental model corrections. Corrections to physical properties are associated with large-scale biases between model and data, with diverse characteristics in the open-ocean and the shelves. Mesoscale features imprint their signature on temperature and sea-level corrections, as well as on chlorophyll corrections due to the vertical velocities associated with vortices. Small scale local corrections are visible over the shelves. Chlorophyll information has measurable impact on physical variables.

In this Part 2 article of a two-part series, observations based on satellite missions were used to evaluate the empirical consistency of model ensembles generated via stochastic modelling of ocean physics and biogeochemistry. A high-resolution Bay of Biscay configuration was used as a case study to explore the model error subspace in both the open and coastal ocean. In Part 1 of this work, three experiments were carried out to generate model ensembles by perturbing only physics, only biogeochemistry, and both of them simultaneously. In Part 2 of this work, empirical consistency was checked, first by means of rank histograms projecting the data onto the model ensemble classes, and second, by pattern-selective consistency criteria in the space of “array modes” defined as eigenvectors of the representer matrix. Rank histograms showed large dependency on geographical region and on season for sea surface temperature (SST), sea-level anomaly (SLA), and phytoplankton functional types (PFT), shifting from consistent model-data configurations to large biases because of model ensemble underspread. Consistency for SST array modes was found to be verified at large, small and coastal scales soon after the ensemble spin-up. Array modes for the along-track sea-level showed useful consistent information at large scales and at the mesoscale; for the gridded SLA was verified only at large scale. Array modes showed that biogeochemical model uncertainties generated by stochastic physics, were effectively detected by PFT measurements at large scales, as well as at mesoscale and small-scale. By contrast, perturbing only biogeochemistry, with an identical physical forcing across the ensemble, limits the potential of PFT measurements at detecting and possibly correcting small-scale biogeochemical model errors. When an ensemble was found to be inconsistent with observations along a particular direction (here, an array mode), a plausible reason is that other error processes must have been active in the model, in addition to the ones at work across the ensemble.

A rich data set on particulate matter optical properties and parameters (beam attenuation coefficient, volume concentration, particle size and PSD slope), accompanied by measurements of biochemical indices (particulate organic carbon, particulate nitrogen and their stable isotopic composition) was obtained from the surface to deep waters across the Mediterranean Sea, in March-April 2018. A decrease of beam attenuation coefficients, total particle volume concentrations, particulate organic carbon and nitrogen concentrations was noted towards the eastern Mediterranean Sea (EMed) in comparison to the western Mediterranean Sea (WMed). LISST-derived optical properties were significantly correlated with water mass characteristics. Overall, the most turbid water mass identified in the Mediterranean Sea was the Surface Atlantic water (AW), and the most transparent was the Transitional Mediterranean Water (TMW) in the Cretan Sea, whereas a general decrease in particulate matter concentration is observed from the WMed towards the EMed. Relatively depleted δ13C-POC values in the particle pool of the open Mediterranean Sea can be attributed to contribution from terrestrial inputs, mainly via atmospheric deposition. Throughout the entire water column, a significant positive correlation between particle beam attenuation coefficient and particulate organic carbon concentration is observed in the open Mediterranean Sea. Such relationship suggests the predominance of organic particles with biogenic origin. POC concentration and particle median diameter D50 are significantly and negatively correlated both in the WMed and the EMed Sea, confirming that small particles are POC-rich. At depth, a prominent decrease of most measured parameters was observed, with the exception of particle median diameter that increased substantially in the EMed towards the deep sea, suggesting potentially enhanced aggregation processes. The low particle size distribution slope ξ observed in the EMed, corresponding to larger particle populations, supports the above notion. Basin-wide Rayleigh-type isotopic fractionation in vertical profiles of δ15N-PN across the Mediterranean Sea, underlines the differences in the trophic characters of the two sub-basins and highlights the role of circulation changes on biogeochemical parameters and the redistribution of particulate matter as a source of nutrients in the water column.