This paper is the outcome of a workshop held in Rome in November 2011 on the occasion of the 25th anniversary of the POEM (Physical Oceanography of the Eastern Mediterranean) program. In the workshop discussions, a number of unresolved issues were identified for the physical and biogeochemical properties of the Mediterranean Sea as a whole, i.e., comprising the Western and Eastern sub-basins. Over the successive two years, the related ideas were discussed among the group of scientists who participated in the workshop and who have contributed to the writing of this paper. Three major topics were identified, each of them being the object of a section divided into a number of different sub-sections, each addressing a specific physical, chemical or biological issue: 1. Assessment of basin-wide physical/biochemical properties, of their variability and interactions. 2. Relative importance of external forcing functions (wind stress, heat/moisture fluxes, forcing through straits) vs. internal variability. 3. Shelf/deep sea interactions and exchanges of physical/biogeochemical properties and how they affect the sub-basin circulation and property distribution. Furthermore, a number of unresolved scientific/methodological issues were also identified and are reported in each sub-section after a short discussion of the present knowledge. They represent the collegial consensus of the scientists contributing to the paper. Naturally, the unresolved issues presented here constitute the choice of the authors and therefore they may not be exhaustive and/or complete. The overall goal is to stimulate a broader interdisciplinary discussion among the scientists of the Mediterranean oceanographic community, leading to enhanced collaborative efforts and exciting future discoveries.

Microplastic research has mainly concentrated on open seas, while riverine plumes remain largely unexplored despite their hypothesized importance as a microplastic source to coastal waters. This work aimed to model coastal accumulation of microplastic particles (1-5 mm) emitted by the Po River over 1.5 years. We posit that river-induced microplastic accumulation on adjacent coasts can be predicted using (1) hydrodynamic-based and (2) remote sensing-based modelling. Model accumulation maps were validated against sampling at nine beaches, with sediment microplastic concentrations up to 78 particles/kg (dry weight). Hydrodynamic modelling revealed that discharged particle amount is only semi-coupled to beaching rates, which are strongly mouth dependent and occur within the first ten days. Remote sensing modelling was found to better capture river mouth relative strength, and accumulation patterns were found consistent with hydrodynamic modelling. This methodology lays groundwork for developing an operational monitoring system to assess microplastic pollution emitted by a major river.

In November 2011, an Atlantic depression affected the Mediterranean basin, eventually evolving into a Tropical-Like Cyclone (TLC or Mediterranean Hurricane, usually designated as Medicane). In the region affected by the Medicane, mean sea level pressures down to 990 hPa, wind speeds of hurricane intensity close to the eye (around 115 km/h) and intense rainfall in the prefrontal zone were reported. The intensity of this event, together with its long permanence over the sea, suggested its suitability as a paradigmatic case for investigating the sensitivity of a numerical modeling system to different configurations, air-sea interface parameterizations and coupling approaches. Toward this aim, a set of numerical experiments with different parameterization schemes and levels of coupling complexity was carried out within the Coupled Ocean Atmosphere Wave Sediment Transport System (COAWST), which allows the description of air-sea dynamics by coupling an atmospheric model (WRF), an ocean circulation model (ROMS), and a wave model (SWAN). The sensitivity to different initialization times and Planetary Boundary Layer (PBL) parameterizations was firstly investigated by running a set of WRF standalone (atmospheric-only) simulations. In order to better understand the effect of coupling on the TLC formation, intensification and trajectory, different configurations of atmosphere-ocean coupling were subsequently tested, eventually including the full coupling among atmosphere, ocean and waves, also changing the PBL parameterization and the formulation of the surface roughness. Results show a strong sensitivity of both the trajectory and the intensity of this TLC to the initial conditions, while the tracks and intensities provided by the coupled modeling approaches explored in this study do not introduce drastic modifications with respect to those resulting from a fine-tuned standalone atmospheric run, though they provide by definition a better physical and energetic consistency. Nevertheless; the use of different schemes for the calculation of the surface roughness from wave motion, which reflects the description of air-sea interface processes, can significantly affect the results in the fully coupled runs.

In this paper we document dense water formation throughout the Adriatic shelf and coastal area in January/February 2012, resulting in record-breaking densities observed during and after the event. The unprecedented dense water generation was preconditioned by a dry and warm year which resulted in a significant reduction of coastal freshwaters, superimposed on a long-term basin-wide salinity increase. The final event that triggered the dense water formation was an extended period of cold weather with strong and severe winds. Record-breaking potential density anomalies (above 30 kg m⁻³) were measured at several formation sites. Accumulated surface net heat and water losses in some coastal regions exceeded 1.5 GJ m⁻² and 250 kg m⁻² over 21 days, respectively. Excessiveness, importance of shelf-type dense water formation and effects on the thermohaline circulation and deep aquatic systems are discussed.

Sea surface temperature (SST) and turbidity (T) derived from Landsat 8 (L8) imagery were used to characterize river plumes in the northern Adriatic Sea (NAS) during a significant flood event in November 2014. Circulation patterns and sea surface salinity (SSS) from an operational coupled ocean-wave model supported the interpretation of the plumes' interaction with the receiving waters and among them. There was a good agreement of the SSS, T, and SST fields at the sub-mesoscale and mesoscale delineation of the major river plumes. L8 30 m resolution also enabled the description of smaller plume structures. The different plumes' reflectance spectra were related to the lithological fingerprint of the sediments in the river catchments. Sharp fronts in T and SST delimited each single river plume. The isotherms and turbidity isolines' coupling varied among the plumes due to differences in particle loads and surface temperatures in the discharged waters. The surface expressions of all the river plumes occurring in NAS were classified based on the occurrence of the plume dynamical regions in the L8 30 m resolution imagery.