(Supplementary Material 2) Landmark coordinates of the analyzed specimens with associated metadata

The use of planktonic foraminifera in paleoceanographic studies relies on the assumption that morphospecies represent biological species with ecological preferences that are stable through time and space. However, genetic surveys unveiled a considerable level of diversity in most morphospecies of planktonic foraminifera. This diversity is significant for paleoceanographic applications because cryptic species were shown to display distinct ecological preferences that could potentially help refine paleoceanographic proxies. Subtle morphological differences between cryptic species of planktonic foraminifera have been reported, but so far their applicability within paleoceanographic studies remains largely unexplored. Here we show how information on genetic diversity can be transferred to paleoceanography using Globorotalia inflata as a case study. The two cryptic species of G. inflata are separated by the Brazil-Malvinas Confluence (BMC), a major oceanographic feature in the South Atlantic. Based on this observation, we developed a morphological model of cryptic species detection in core top material. The application of the cryptic species detection model to Holocene samples implies latitudinal oscillations in the position of the confluence that are largely consistent with reconstructions obtained from stable isotope data. We show that the occurrence of cryptic species in G. inflata, can be detected in the fossil record and used to trace the migration of the BMC. Since a similar degree of morphological separation as in G. inflata has been reported from other species of planktonic foraminifera, the approach presented in this study can potentially yield a wealth of new paleoceanographical proxies.

Terrain model of the eastern slope of the Porcupine Seabight (50 m grid size)

The topography of the eastern margin of the Porcupine Seabight was surveyed in June 2000 utilizing swath bathymetry. The survey was carried out during RV Polarstern cruise ANT XVII/4 as part of the GEOMOUND project. The main objective was to map and investigate the seafloor topography of this region. The investigated area contains a variability of morphological features such as deep sea channels and giant mounds. The survey was planned and realized on the basis of existing data so as to guarantee the complete coverage of the margin. In order to achieve a resolution of the final digital terrain model (DTM) that meets the project demands, data processing was adjusted accordingly. The grid spacing of the DTM was set to 50 m and an accuracy better than 1% of the water depth was achieved for 96% of the soundings.

Terrain model of the Håkon Mosby Mud Volcano (10 m grid size)

The Håkon Mosby Mud Volcano is a natural laboratory to study geological, geochemical, and ecological processes related to deep-water mud volcanism. High resolution bathymetry of the Håkon Mosby Mud Volcano was recorded during RV Polarstern expedition ARK-XIX/3 utilizing the multibeam system Hydrosweep DS-2. Dense spacing of the survey lines and slow ship speed (5 knots) provided necessary point density to generate a regular 10 m grid. Generalization was applied to preserve and represent morphological structures appropriately. Contour lines were derived showing detailed topography at the centre of the Håkon Mosby Mud Volcano and generalized contours in the vicinity. We provide a brief introduction to the Håkon Mosby Mud Volcano area and describe in detail data recording and processing methods, as well as the morphology of the area. Accuracy assessment was made to evaluate the reliability of a 10 m resolution terrain model. Multibeam sidescan data were recorded along with depth measurements and show reflectivity variations from light grey values at the centre of the Håkon Mosby Mud Volcano to dark grey values (less reflective) at the surrounding moat.

Medem Channel profiles based on stand-alone model simulation and corrections with nudging

State-of-the-art process-based models have shown to be applicable to the simulation and prediction of coastal morphodynamics. On annual to decadal temporal scales, these models may show limitations in reproducing complex natural morphological evolution patterns, such as the movement of bars and tidal channels, e.g. the observed decadal migration of the Medem Channel in the Elbe Estuary, German Bight. Here a morphodynamic model is shown to simulate the hydrodynamics and sediment budgets of the domain to some extent, but fails to adequately reproduce the pronounced channel migration, due to the insufficient implementation of bank erosion processes. In order to allow for long-term simulations of the domain, a nudging method has been introduced to update the model-predicted bathymetries with observations. The model-predicted bathymetry is nudged towards true states in annual time steps. Sensitivity analysis of a user-defined correlation length scale, for the definition of the background error covariance matrix during the nudging procedure, suggests that the optimal error correlation length is similar to the grid cell size, here 80-90 m. Additionally, spatially heterogeneous correlation lengths produce more realistic channel depths than do spatially homogeneous correlation lengths. Consecutive application of the nudging method compensates for the (stand-alone) model prediction errors and corrects the channel migration pattern, with a Brier skill score of 0.78. The proposed nudging method in this study serves as an analytical approach to update model predictions towards a predefined 'true' state for the spatiotemporal interpolation of incomplete morphological data in long-term simulations.

Water level time series and sediment grain size model

The environment of ebb-tidal deltas between barrier island systems is characterized by a complex morphology with ebb- and flood-dominated channels, shoals and swash bars connecting the ebb-tidal delta platform to the adjacent island. These morphological features reveal characteristic surface sediment grain-size distributions and are subject to a continuous adaptation to the prevailing hydrodynamic forces. The mixed-energy tidal inlet Otzumer Balje between the East Frisian barrier islands of Langeoog and Spiekeroog in the southern North Sea has been chosen here as a model study area for the identification of relevant hydrodynamic drivers of morphology and sedimentology. We compare the effect of high-energy, wave-dominated storm conditions to mid-term, tide-dominated fair-weather conditions on tidal inlet morphology and sedimentology with a process-based numerical model. A multi-fractional approach with five grain-size fractions between 150 and 450 µm allows for the simulation of corresponding surface sediment grain-size distributions. Net sediment fluxes for distinct conditions are identified: during storm conditions, bed load sediment transport is generally onshore directed on the shallower ebb-tidal delta shoals, whereas fine-grained suspended sediment bypasses the tidal inlet by wave-driven currents. During fair weather the sediment transport mainly focuses on the inlet throat and the marginal flood channels. We show how the observed sediment grain-size distribution and the morphological response at mixed-energy tidal inlets are the result of both wave-dominated less frequent storm conditions and mid-term, tide-dominant fair-weather conditions.

Morphological variability of Globorotalia menardii in DSDP Sites 68-502 and 68-503

Variability in the test of Globorotalia menardii during the past 8 million years has been investigated at DSDP Site 502A (Caribbean Sea) and DSDP Site 503A (Eastern Equatorial Pacific). Measurements were made of spire height (delta x), maximum diameter (delta y), the tangent angles of the upper and lower peripheral keels (phi 1, phi 2, respectively), the number of chambers in the final whorl, and the area of the silhouette in keel view. Four morphotypes alpha, beta, gamma, and delta were distinguished. Morphotype alpha was found in strata ranging in age from the Late Miocene through the Holocene. It shows a continuous increase in delta x and delta y until the Late Pleistocene. During and after the final closure of the ancient Central American Seaway (between 2.4 Ma and 1.8 Ma) there was a rapid increase in the area of the test in keel view. At the Caribbean Sea site, morphotype beta evolved during the past 0.22 Ma. It is less inflated than alpha and has a more delicate test. In the morphospace of delta x vs. delta y, morphotypes alpha and beta can be distinguished by a separation line delta y = 3.2 * delta x - 160 (delta x and delta y in µm). Plots of morphotype alpha are below that line, those of beta are above it. Morphotype alpha is taken to be Globorotalia menardii menardii Parker, Jones & Brady (1865) and includes G. menardii 'A' Bolli (1970). Morphotype beta is identified as G. menardii cultrata (d'Orbigny). Morphotypes gamma and delta are extinct Upper Miocene to Pliocene forms which evolved from morphotype alpha. They have a narrower phi 1 angle and more chambers (>=7) than morphotype alpha commonly with 5 to 6 chambers (7 in transitional forms). In contemporaneous samples morphotype delta can be distinguished from gamma by a smaller value of phi 1 and 8 or more chambers in the final whorl. Morphotype gamma is taken to be G. limbata (Fornasini, 1902) and includes the junior synonym G. menardii 'B' Bolli (1970). Morphotype delta is G. multicamerata Cushman & Jarvis (1930). With the exception of the Late Pleistocene development of G. menardii cultrataonly in the Caribbean the morphological changes of G. menardii at DSDP Sites 502A and 503A are similar. The development from the ancestral G. menardii menardii of the G. limbata - G. multicamerata lineage during the Pliocene and of G. menardii cultrata during the Late Pleistocene suggests responses at the two sites to a changing palaeoceanography during and after the formation of the Isthmus of Panama.

Results of analysis and model simulation of the bed profiles

Bedforms such as dunes and ripples are ubiquitous in rivers and coastal seas, and commonly described as triangular shapes from which height and length are calculated to estimate hydrodynamic and sediment dynamic parameters. Natural bedforms, however, present a far more complicated morphology; the difference between natural bedform shape and the often assumed triangular shape is usually neglected, and how this may affect the flow is unknown. This study investigates the shapes of natural bedforms and how they influence flow and shear stress, based on four datasets extracted from earlier studies on two rivers (the Rio Paraná in Argentina, and the Lower Rhine in The Netherlands). The most commonly occurring morphological elements are a sinusoidal stoss side made of one segment and a lee side made of two segments, a gently sloping upper lee side and a relatively steep (6 to 21°) slip face. A non-hydrostatic numerical model, set up using Delft3D, served to simulate the flow over fixed bedforms with various morphologies derived from the identified morphological elements. Both shear stress and turbulence increase with increasing slip face angle and are only marginally affected by the dimensions and positions of the upper and lower lee side. The average slip face angle determined from the bed profiles is 14°, over which there is no permanent flow separation. Shear stress and turbulence above natural bedforms are higher than above a flat bed but much lower than over the often assumed 30° lee side angle.