Search

Submesoscale features with profound impact on ocean dynamics and climate-relevant fluxes are frequently observed in the upper ocean including Arctic region. Yet, modelling these features remains a challenge due to the difficulties in the parameterization of submesoscale processes and high resolution required, in particular, in the polar regions. The most effective way to study such phenomena is joint modelling and observational work. Several autonomous observation platforms have been deployed as part of Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) experiment within an approximately 50 km radius around the central observatory. Data from these buoys in combination with data from the central observatory provide a unique opportunity to reconstruct 3D water properties and velocity by constraining a numerical model that resolves the dynamics of the (sub-)mesoscale. It turns out that a minimum root mean square error between results of an optimal interpolation and observations indicates a characteristic length scale of about 7.5 km, corresponding approximately the first-mode barolinic Rossby radius in the area of investigation. However, results of the interpolation are questionable at the sub-mesoscale due to the distribution of the buoy observations in time and horizontal space. In order to describe the in-situ data to achieve a better characterization and understanding of (sub-)mesoscale dynamics we developed and applied a modification of the 3D regional model FESOM-C. The observed temperature and salinity were used to nudge the model to obtain an optimized solution at the resolution of the models. A series of simulations with different horizontal resolutions and model parameters make it possible to analyze the ability of models of this type to reproduce the observed dynamics, to estimate eddy kinetic energy and power spectra, and to compare findings with the observations used to nudge the model. We will show the eddy-induced fluxes and characteristics of eddies along the track of the beginning winter MOSAiC drift.

Using the depth (z) and density (ϱ) frameworks, we analyze local contributions to AMOC variability in a 900‐year simulation with the AWI climate model. Both frameworks reveal a consistent interdecadal variability; however, the correlation between their maxima deteriorates on year‐to‐year scales. We demonstrate the utility of analyzing the spatial patterns of sinking and diapycnal transformations through depth levels and isopycnals. The success of this analysis relies on the spatial binning of these maps which is especially crucial for the maps of vertical velocities which appear to be too noisy in the main regions of upwelling and downwelling because of stepwise bottom topography. Furthermore, we show that the AMOC responds to fast (annual or faster) fluctuations in atmospheric forcing associated with the NAO. This response is more obvious in the ϱ than in the z framework. In contrast, the link between AMOC and deep water production south of Greenland is found for slower fluctuations and is consistent between the frameworks.

Recent studies on seafloor mapping have presented different modelling methods for the automatic classification of seafloor sediments. However, most of these studies have applied these models to seafloor data with appropriate number of ground-truth samples, which raises the question whether these methods are applicable to studies with smaller numbers of ground-truth data. In this study, we aim to address this issue by conducting sediment class-specific predictions using ensemble modelling to map areas with limited or without ground-truth data and combined with hydro-acoustic datasets. The resulting class-specific maps were then assembled into one map, where the most probable class was assigned to the appropriate location. Our approach was able to predict sediment classes without bias to the class with more ground-truth data and produced reliable seafloor sediment distributions maps that can be used for seafloor monitoring. Sediment shifts of a heterogenous seafloor in the Sylt Outer Reef, German North Sea were also assessed to understand the sediment dynamics in the area. The analyses of sediment shifts showed that the western area of the Sylt Outer Reef is highly active, and the results of the analyses assisted in providing recommendations on future seafloor monitoring activities.

The dynamics of the boundary layer of the ocean significantly affect the interaction between ocean and atmosphere and, as a result, global climate. The sub-ice boundary layer of the ocean and its dynamics have not been thoroughly studied because of the extremely difficult conditions for observation, in particular during winter. Current understanding of spatial-temporal variability of (sub)mesoscales of the upper Arctic Ocean is extremely limited.At the same time, one of the most important features of the upper ocean layers are the small-scale processes that influence and possibly determine the vertical and horizontal transport of heat, salt, and biologically relevant substances. As a consequence, mathematical models, in particular climate models, experience serious difficulties in parameterization of processes not resolved by the models because of the lack sufficient knowledge to detail the spatial variability at the (sub-)mesoscale.To a better characterization and understanding of (sub)mesoscale dynamics and its role in vertical transport of energy and mass we apply a 3D regional ocean model FESOM-C. The observed vertical hydrological structure and a corresponding reconstructed horizontal temperature and salinity fields were imposed as a part of the forcing for the numerical model. These fields and information about the vertical hydrological structure were utilized by the model as initial conditions and for constraining (nudging) during the spin-up period. After the initial spin-up period, once the model had adjusted to our initial conditions, we performed several free runs.We expect that our 3D numerical studies of eddy properties will contribute to a better characterisation and understanding of (sub)mesoscale dynamics in the Arctic Ocean and its role in the vertical transport of energy and mass.

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.

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 a secular linear trend 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 understanding, estimating and predicting the individual contributions to local relative sea level, with focus on the most recent publications. The current best estimate of historical sea-level rise in Venice, based on tide-gauge data after removal of subsidence effects, is 1.23 ± 0.13 mm/year (period from 1872 to 2019). Subsidence thus contributed to about half of the observed relative sea-level rise over the same period. A higher – yet more uncertain – rate of sea-level rise is observed during recent decades, estimated from tide-gauge data to be about 2.76 ± 1.75 mm/year in the period 1993–2019 for the climatic component alone. An unresolved issue is the contrast between the observational capacity of tide gauges and satellite altimetry, with the latter tool not covering the Venice Lagoon. Water mass exchanges through the Gibraltar Strait currently constitute a source of substantial uncertainty for estimating future deviations of the Mediterranean mean sea-level trend from the global-mean value. Subsidence and regional atmospheric and oceanic circulation mechanisms can deviate Venetian relative sea-level trends from the global mean values for several decades. Regional processes will likely continue to determine significant interannual and interdecadal variability of Venetian sea level with magnitude comparable to that observed in the past, as well as non-negligible differential trends. Our estimate of the likely range of mean sea-level rise in Venice by 2100 due to climate change is presently estimated between 11 and 110 centimetres. An improbable yet possible high-end scenario linked to strong ice-sheet melting yields about 170 centimetres of mean sea-level rise in Venice by 2100. Projections of natural and human induced vertical land motions are currently not available, but historical evidence demonstrates that they can produce a significant contribution to the relative sea-level rise in Venice, further increasing the hazard posed by climatically-induced sea-level changes.

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 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, and the Højer and Rømø Deep areas. 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. 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 fulfillment of the energy balance. 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 are necessary and useful benchmarks for further studies in the area, including baroclinic and sediment dynamics investigations.

Detecting changes of sediment boundaries on the seafloor is important for a better understanding of sediment dynamics and related impacts to benthic habitats. Side-scan sonars (SSS) perform more cost-effectively in shallow waters than other acoustic systems because of their larger swath widths, and the resolution of its images does not change with varying water depth. However, as they are generally towed behind the survey vessel, they tend to have lower positioning accuracy, which makes them unreliable for change detection analyses. In this study, we present a workflow that processes SSS data in a way that makes them fit for change detection analyses. To test the capacity of SSS mosaics for change detection, we used a free software called “Digital Shoreline Analysis System”, which was developed by the United States Geological Survey for ArcGIS version 10.4 onwards. The methods were applied in three areas in the Sylt Outer Reef, German Bight, North Sea. Our results showed that with appropriate processing, SSS mosaics could be used for change detection of sharp sediment boundaries. We found a common trend in the sediment distribution patterns of coarse sediments by monitoring the movement of their boundaries. The boundaries moved in northeast-southwest direction and boundary movements of less than 20 m were typically observed. The methods presented here are semi-automated, repeatable, and replicable, which has potential for wide-scale monitoring of sediment distribution patterns.

A newly developed coastal model FESOM-C based on three-dimensional unstructured meshes and finite volume is applied to simulate dynamics of the southeastern part of the North Sea. Variable horizontal resolution enables using meshes that are coarse in the open sea but refined in the shallow areas (which include the Wadden Sea and the estuaries) to resolve important small-scale process (such as wetting and drying, sub-mesoscales eddies and dynamics of steep coastal fronts). Model results for the simulation for the period between January 2010 and December 2014 agree reasonably well with data from numerous autonomous observation stations with high temporal and spatial resolution, located in the region, data from ferry boxes and glider expeditions. The analysis of numerical solution convergence on meshes with different horizontal resolutions allows identifying areas where high mesh resolution (wetting and drying zones, shallow areas) and low mesh resolution (open boundary, open sea, and deep regions) are optimal for numerical simulations.

The Lena River (Lena R.) heat flux affects the Laptev Sea hydrology. Published long-term estimates range from 14.0 to 15.7 EJ·a−1, based on data from Kyusyur, at the river outlet. A novel daily stream temperature (Tw) dataset was used to evaluate contemporary Lena R. heat flux, which is 16.4 ± 2.7 EJ·a−1 (2002–2011), confirming upward trends in both Tw and water runoff. Our field data from Kyusyur, however, reveal a significant negative bias, −0.8 °C in our observations, in observed Tw values from Kyusyur compared to the cross-section average Tw. Minor Lena R. tributaries discharge colder water during July–September, forming a cold jet affecting Kyusyur Tw data. Major Tw negative peaks mostly coincide with flood peaks on the Yeremeyka River, one of these tributaries. This negative bias was accounted for in our reassessment. Revised contemporary Lena R. heat flux is 17.6 ± 2.8 EJ·a−1 (2002–2011) and is constrained from above at 26.9 EJ·a−1 using data from Zhigansk, approximately 500 km upstream Kyusyur. Heat flux is controlled by stream temperature in June, during the freshet period, while from late July to mid-September, water runoff is a dominant factor.

The long-wave dynamics of the Lombok Strait, which is the most important link of the West Indonesian throughflow connecting the Pacific and Indian Ocean waters, was simulated and analyzed. A feature of the strait is its extremely complex relief, on which water transport creates a field of pronounced vertical velocities, which requires consideration of the nonhydrostatic component of pressure. The work presents a 3-D nonhydrostatic model in curvilinear coordinates, which is verified on a test problem. Particular attention is paid to the method of solving the 3-D elliptical solver for a nonhydrostatic problem in boundary-matched coordinates and a vertical &sigma, level. The difference in transport through the Lombok Strait is determined by the difference in atmospheric pressure over the Pacific and Indian Oceans. Based on the results of the global simulation, the role of these factors in terms of their variability is analyzed, and the value of nonhydrostatic pressure in the dynamics of the Lombok Strait is revealed and evaluated. The vertical dynamics of the Lombok Strait are considered in detail based on hydrostatic and nonhydrostatic approaches.

Monitoring changes in eukaryotic microbial communities is critical for understanding ecosystem dynamics, trophic interactions and the impacts of climate change. Long-term time series are an important tool for monitoring changes in ecological communities, but time series from a single location may not be representative of regional dynamics. In the German Bight, the Helgoland Roads time series is such a long-term series. Here, we consider the spatial dynamics of the eukaryotic microbes as an indicator of the representativeness of the Helgoland Roads site for the coastal German Bight, which is located in the North Sea. The eukaryotic microbial community in the German Bight was analysed at Helgoland Roads and two coastal stations (Cuxhaven and Wilhelmshaven) between March and October 2016 using metabarcoding. In addition, an oceanographical model was used to check for potential hydrological connectivity between the stations during the sampling period. Our results showed that the communities were different at the three stations. Helgoland was dominated by dinoflagellates, whereas the coastal stations had more diverse communities. Furthermore, differences were observed in the dinoflagellate and diatom communities between the three stations. Lagrangian particle tracking applied to the model results, showed limited connectivity between Helgoland and the coastal stations in 2016. The differences between Helgoland and the coastal stations were correlated with the different hydrological regimes and associated nutrient contents. Our observations suggest the presence of different eukaryotic microbial communities separated by complex hydrological conditions in the coastal German Bight.

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 by the implementation of time stepping, the use of terrain-following vertical coordinate and 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 differt resolution or instead of strongly deformed quadrilateral cells. The description of the model numerical part is complemented by test cases illustrating the model performance.

The Lena plume dynamics in the Lena Delta region of the Laptev Sea are explored in simulations performed with the Finite Volume Coastal Ocean Model (FVCOM) on a mesh with the horizontal resolution 0.4–5 km. The impact of wind and tides on the Lena plume propagation is analysed based on simulations for the summer season of 2008 and also on idealised experiments. All main Lena River freshwater channels (Trofimovskaya, Bykovskaya, Tumatskaya and Olenekskaya) produce buoyant outflows in the summer season. The surface plume buoyancy signature proves to be highly variable in time, especially in case of upwelling favourable wind events. Winds stronger than 6 m s−1 can already turn the dynamics of flows from all main freshwater channels to the wind-driven state. During the summer season, the bulk of freshwater from the Lena River stays in the eastern Laptev Sea because of location of the main Lena River freshwater channels, their large Kelvin numbers and light summer winds. Westward and northward plume excursions are wind-driven, and the model skill in simulating them depends on the available wind forcing. The main mechanism of tidal influence in the freshwater plume zone is through tidally induced mixing, except for the northern vicinity of the delta, where residual circulation may contribute to the plume eastward transport significantly.

Tidal processes play an important role in the dynamics of shelf circulation in the Laptev Sea. The Unstructured Grid Finite Volume Coastal Ocean Model (FVCOM) is used to simulate the tidal dynamics in the Lena Delta region of the Laptev Sea in ice-free barotropic case. The grid element size ranges from 400 m to 5 km. The major semidiurnal tidal waves M2M2 and S2S2 are investigated with the M2M2 being the most important in generating large sea level amplitudes and currents over shallow areas. A correction to the tidal elevation at the open boundary is proposed, which minimizes the discrepancy between the model prediction and observations. The observations include both recent mooring data and the standard set of tide gauge measurements used in previous studies. The comparison of results to known tidal solutions is carried out. The paper also discusses the residual circulation and energy fluxes and assesses the impact of additional bathymetric information. Highlights • Our simulations reproduce the semidiurnal tidal waves M2M2 and S2S2 in the Laptev Sea region. • We develop special procedure for the construction of optimal open boundary conditions for tidal elevation for M2M2 and S2S2 constituents. • The simulated tidal maps show an improved agreement with observations. • We analyze barotropic currents, residual circulation and evolution of energy fluxes in the region. • We consider the energy balance for the M2M2 and S2S2 waves and the sensitivity to the bathymetry.

The water temperature characteristics of the Lena River at basin outlet during the summer season (June–September) are considered. The analysis is based on a long-term data series covering the period from the beginning of observation (1936) to the present time (2012) at Kusur (Kyusyur) gauging station and complementary data at several stations downstream and one station upstream. These additional data are rarely used, but their analysis is important for understanding processes in the basin outlet area. The differences between the stream surface temperatures at Kusur station and 200 km downstream to the north at Habarova (Khabarova) station have almost always been an anomalously large and negative for the considered period since the beginning of observation during open water season from July to September. The description of this difference and its analysis are presented. To sort the problem out, we consider the observational data in terms of the hydrology and morphology of the Lena River delta and main channel area and apply statistical and deterministic modelling approaches. The inability of water temperature observational data at Kusur station to represent the mean cross-sectional temperature is addressed. The analysis of the water temperature trends at the Kusur and Habarova stations is also presented.

The data series from several gauging stations on the Lena River are considered. The difference of the water temperature between Kusur station and 200 kilometers downstream to the north Habarova station often becomes anomaly large during open water season (June-September). The analysis of the observational data is presented. The non-representativeness of measurements on Kusur station is supposed, reconstruction of the real temperature regime on the base of the statistical and deterministic modeling as well as the optimization tools is carried out.