Coastal zones connect terrestrial and marine ecosystems forming a unique environment that is under increasing anthropogenic pressure. Rising sea levels, sinking coasts, and changing precipitation patterns modify hydrodynamic gradients and may enhance sea-land exchange processes in both tidal and non-tidal systems. Furthermore, the removal of flood protection structures as restoration measure contributes locally to the changing coastlines. A detailed understanding of the ecosystem functioning of coastal zones and the interactions between connected terrestrial and marine ecosystems is still lacking. Here, we propose an interdisciplinary approach to the investigation of interactions between land and sea at shallow coasts, and discuss the advantages and the first results provided by this approach as applied by the research training group Baltic TRANSCOAST. A low-lying fen peat site including the offshore shallow sea area on the southern Baltic Sea coast has been chosen as a model system to quantify hydrophysical, biogeochemical, sedimentological, and biological processes across the land-sea interface. Recently introduced rewetting measures might have enhanced submarine groundwater discharge (SGD) as indicated by distinct patterns of salinity gradients in the near shore sediments, making the coastal waters in front of the study site a mixing zone of fresh-and brackish water. High nutrient loadings,

Microplastics (MP) data collection from the aquatic environment is a challenging endeavour that sets apparent limitations to regional and global MP quantification. Expensive data collection causes small sample sizes and oftentimes existing data sets are compared without accounting for natural variability due to hydrodynamic processes governing the distribution of particles. In Warnow estuarine sediments (Germany) we found significant correlations between high-density polymer size fractions (≥500 µm) and sediment grain size. Among potential predictor variables (source and environmental terms) sediment grain size was the critical proxy for MP abundance. The MP sediment relationship can be explained by the force necessary to start particle transport: at the same level of fluid motion, transported sediment grains and MP particles are offset in size by one to two orders of magnitude. Determining grain-size corrected MP abundances by fractionated granulometric normalisation is recommended as a basis for future MP projections and identification of sinks and sources.

The well-known Knudsen relations and the total exchange flow (TEF) analysis framework provide quantifications of exchange flow across an open boundary to the adjacent ocean in terms of bulk values (Knudsen theory: inflow and outflow volume or salinity) or with resolution in salinity space (TEF: profiles of volume and salt flux in salinity coordinates). In the present study, these theories are extended toward mixing of salinity, defined as the decay of salinity variance due to turbulent mixing. In addition to the advective fluxes, diffusive fluxes across the boundary are also considered now. These new Knudsen and TEF relations for mixing are derived by applying Gauss’s theorem to the salinity square and salinity variance equations. As a result of the analysis, four different Knudsen relations for the mixing in estuaries are derived. The first one is exact and considers nonperiodicity as well as nonconstancy of the inflow and outflow salinities. The other three formulations are approximate only, in the sense that either nonperiodicity or nonconstancy or both are relaxed. The simplest of those formulations has recently been derived by MacCready et al. and estimates the estuarine mixing as the product of inflow salinity, outflow salinity, and time-averaged river runoff. These four mixing estimates are systematically assessed by means of a number of idealized estuarine test cases. For periodic tidal flow, the simplest estimate still predicts the effective (physical plus numerical) mixing within an error of about 10%.

In straight tidal estuaries, residual overturning circulation results mainly from a competition between gravitational forcing, wind forcing, and friction. To systematically investigate this for tidally energetic estuaries, the dynamics of estuarine cross sections is analyzed in terms of the relation between gravitational forcing, wind stress, and the strength of estuarine circulation. A system-dependent basic Wedderburn number is defined as the ratio between wind forcing and opposing gravitational forcing at which the estuarine circulation changes sign. An analytical steady-state solution for gravitationally and wind-driven exchange flow is constructed, where tidal mixing is parameterized by parabolic eddy viscosity. For this simple but fundamental situation, is calculated, meaning that the up-estuary wind forcing needs to be 15% of the gravitational forcing to invert estuarine circulation. In three steps, relevant physical processes are added to this basic state: (i) tidal dynamics are resolved by a prescribed semidiurnal tide, leading to caused by tidal straining; (ii) lateral circulation is added by introducing cross-channel bathymetry, smoothly increasing from 0.47 (flat bed) to 1.3 (parabolic bed) due to an increasing effect of lateral circulation on estuarine circulation; and (iii) full dynamics of a real tidally energetic inlet with highly variable forcing, where results from a two-dimensional linear regression.

The present study aims to estimate effective diahaline turbulent salinity fluxes and diffusivities in numerical model simulations of estuarine scenarios. The underlying method is based on a quantification of salinity mixing per salinity class, which is shown to be twice the turbulent salinity transport across the respective isohaline. Using this relation, the recently derived universal law of estuarine mixing, predicting that average mixing per salinity class is twice the respective salinity times the river run‐off, can be directly derived. The turbulent salinity transport is accurately decomposed into physical (due to the turbulence closure) and numerical (due to truncation errors of the salinity advection scheme) contributions. The effective diahaline diffusivity representative for a salinity class and an estuarine region results as the ratio of the diahaline turbulent salinity transport and the respective (negative) salinity gradient, both integrated over the isohaline area in that region and averaged over a specified period. With this approach, the physical (or numerical) diffusivities are calculated as half of the product of physical (or numerical) mixing and the isohaline volume, divided by the square of the isohaline area. The method for accurately calculating physical and numerical diahaline diffusivities is tested and demonstrated for a three‐dimensional idealized exponential estuary. As a major product of this study, maps of the spatial distribution of the effective diahaline diffusivities are shown for the model estuary.