The city of Venice and the surrounding lagoonal ecosystem are highly vulnerable to variations in relative sea level. In the past ∼150 years, this was characterized by an average rate of relative sea-level rise of about 2.5 mm/year resulting from the combined contributions of vertical land movement and sea-level rise. This literature review reassesses and synthesizes the progress achieved in quantification, understanding and prediction of the individual contributions to local relative sea level, with a focus on the most recent studies. Subsidence contributed to about half of the historical relative sea-level rise in Venice. The current best estimate of the average rate of sea-level rise during the observational period from 1872 to 2019 based on tide-gauge data after removal of subsidence effects is 1.23 ± 0.13 mm/year. A higher – but more uncertain – rate of sea-level rise is observed for more recent years. Between 1993 and 2019, an average change of about +2.76 ± 1.75 mm/year is estimated from tide-gauge data after removal of subsidence. Unfortunately, satellite altimetry does not provide reliable sea-level data within the Venice Lagoon. Local sea-level changes in Venice closely depend on sea-level variations in the Adriatic Sea, which in turn are linked to sea-level variations in the Mediterranean Sea. Water mass exchange through the Strait of Gibraltar and its drivers currently constitute a source of substantial uncertainty for estimating future deviations of the Mediterranean mean sea-level trend from the global-mean value. Regional atmospheric and oceanic processes will likely contribute significant interannual and interdecadal future variability in Venetian sea level with a magnitude comparable to that observed in the past. On the basis of regional projections of sea-level rise and an understanding of the local and regional processes affecting relative sea-level trends in Venice, the likely range of atmospherically corrected relative sea-level rise in Venice by 2100 ranges between 32 and 62 cm for the RCP2.6 scenario and between 58 and 110 cm for the RCP8.5 scenario, respectively. A plausible but unlikely high-end scenario linked to strong ice-sheet melting yields about 180 cm of relative sea-level rise in Venice by 2100. Projections of human-induced vertical land motions are currently not available, but historical evidence demonstrates that they have the potential to produce a significant contribution to the relative sea-level rise in Venice, exacerbating the hazard posed by climatically induced sea-level changes.

We describe FESOM-C, the coastal branch of the Finite-volumE Sea ice – Ocean Model (FESOM2), which shares with FESOM2 many numerical aspects, in particular its finite-volume cell-vertex discretization. Its dynamical core differs in the implementation of time stepping, the use of a terrain-following vertical coordinate, and the formulation for hybrid meshes composed of triangles and quads. The first two distinctions were critical for coding FESOM-C as an independent branch. The hybrid mesh capability improves numerical efficiency, since quadrilateral cells have fewer edges than triangular cells. They do not suffer from spurious inertial modes of the triangular cell-vertex discretization and need less dissipation. The hybrid mesh capability allows one to use quasi-quadrilateral unstructured meshes, with triangular cells included only to join quadrilateral patches of different resolution or instead of strongly deformed quadrilateral cells. The description of the model numerical part is complemented by test cases illustrating the model performance.

Background: Blooms of marine microalgae play a pivotal role in global carbon cycling. Such blooms entail successive blooms of specialized clades of planktonic bacteria that collectively remineralize gigatons of algal biomass on a global scale. This biomass is largely composed of distinct polysaccharides, and the microbial decomposition of these polysaccharides is therefore a process of prime importance.Results: In 2020, we sampled a complete biphasic spring bloom in the German Bight over a 90-day period. Bacterioplankton metagenomes from 30 time points allowed reconstruction of 251 metagenome-assembled genomes (MAGs). Corresponding metatranscriptomes highlighted 50 particularly active MAGs of the most abundant clades, including many polysaccharide degraders. Saccharide measurements together with bacterial polysaccharide utilization loci (PUL) expression data identifiedβ-glucans (diatom laminarin) andα-glucans as the most prominent and actively metabolized dissolved polysaccharide substrates. Both substrates were consumed throughout the bloom, withα-glucan PUL expression peaking at the beginning of the second bloom phase shortly after a peak in flagellate and the nadir in bacterial total cell counts.Conclusions: We show that the amounts and composition of dissolved polysaccharides, in particular abundant storage polysaccharides, have a pronounced influence on the composition of abundant bacterioplankton members during phytoplankton blooms, some of which compete for similar polysaccharide niches. We hypothesize that besides the release of algal glycans, also recycling of bacterial glycans as a result of increased bacterial cell mortality can have a significant influence on bacterioplankton composition during phytoplankton blooms.

This study is dedicated to the tidal dynamics in the Sylt-Rømø Bight with a focus on the non-linear processes. The FESOM-C model was used as the numerical tool, which works with triangular, rectangular or mixed grids and is equipped with a wetting/drying option. As the model's success at resolving currents largely depends on the quality of the bathymetric data, we have created a new bathymetric map for an area based on recent studies of Lister Deep, Lister Ley, Højer Deep and Rømø Deep. This new bathymetric product made it feasible to work with high-resolution grids (up to 2 m in the wetting/drying zone). As a result, we were able to study the tidal energy transformation and the role of higher harmonics in the domain in detail. For the first time, the tidal ellipses, maximum tidally induced velocities, energy fluxes and residual circulation maps were constructed and analysed for the entire bight. Additionally, tidal asymmetry maps were introduced and constructed. The full analysis was performed on two grids with different structures and showed a convergence of the results as well as fulfilment of the energy balance. A great deal of attention has been paid to the selection of open-boundary conditions, model validation against tide gauges and recent in situ current data. The tidal residual circulation and asymmetric tidal cycles largely define the circulation pattern, transport and accumulation of sediment, and the distribution of bedforms in the bight; therefore, the results presented in the article are necessary and useful benchmarks for further studies in the area, including baroclinic and sediment dynamics investigations.

A newly developed coastal model, FESOM-C, based on three-dimensional unstructured meshes and finite volume, is applied to simulate the dynamics of the southeastern North Sea. Variable horizontal resolution enables coarse meshes in the open sea with refined meshes in shallow areas including the Wadden Sea and estuaries to resolve important small-scale processes such as wetting and drying, sub-mesoscale eddies, and the dynamics of steep coastal fronts. Model results for a simulation of the period from January 2010 to December 2014 agree reasonably well with data from numerous regional autonomous observation stations with high temporal and spatial resolutions, as well as with data from FerryBoxes and glider expeditions. Analyzing numerical solution convergence on meshes of different horizontal resolutions allows us to identify areas where high mesh resolution (wetting and drying zones and shallow areas) and low mesh resolution (open boundary, open sea, and deep regions) are optimal for numerical simulations.

The Atlantic meridional overturning circulation (AMOC) is one of the most important characteristics of an ocean model run. Using the depth (z) and density frameworks, we analyze how the sinking and diapycnal transformations defining the AMOC as well as AMOC strength and variability react to mesh refinement from low to higher resolution in two model runs driven by the CORE‐II forcing. Both runs can represent the key locations of sinking and diapycnal transformations behind AMOC, that is, northeastern North Atlantic. Although their spatial patterns do not change significantly with resolution in both frameworks as the consequence of the same atmospheric forcing, the quantitative differences, reaching several sverdrups, are seen in different locations between two model runs for both frameworks. In particular, the refinement leads to the strongest differences in the vertical transport and diapycnal transformations in the latitude range between 30°N and 55°N. The z framework emphasizes the role of localized upwelling around the Gulf Stream separation site, whereas the density framework emphasizes the contribution of (spurious) diapycnal mixing around the Grand Banks. Both effects are reduced in the higher‐resolution run, leading to higher AMOC south of 26°N as compared to the low‐resolution run, despite the AMOC maxima, located at high latitudes, are higher in the low‐resolution run. We suggest that both AMOC frameworks should be used routinely in standard analyses, including forthcoming intercomparison projects.