Description: Connectivity is a fundamental process driving the persistence of marine populations and their adaptation potential in response to environmental change. In this study, we analysed the population genetics of two morphologically highly similar deep-sea sponge clades (Phakellia hirondellei and the 'Topsentia-and-Petromica (TaP)' clade) at three locations in the Cantabrian Sea. Sponge taxonomy was assessed by spicule analyses, as well as by 18S sequencing and COI sequencing. The corresponding host microbiome was analysed by 16S rRNA gene sequencing. In addition we set up an oceanographic modelling framework, for which we used seawater flow cytometry data (derived from bottom depths of CTD casts) as ground-truthing data.

Description: A projection of larval dispersal patterns of Atlantic cold seep mussels Gigantidas childressi, G. mauritanicus, Bathymodiolus heckerae and B. boomerang was carried out for the next 50 years under the constraint of global warming predicted by the IPCC for the most pessismistic scenario. Simulations were run at +00 years, +25 years and +50 years from initial years of 2014 to 2019 (+00Y) at 21 locations on the US, European and African coasts using the VIKING20X model, in which the Atlantic water temperatures predicted by the FOCI model were forced to the future dates. The dataset consists of a number of 5775 simulations carried out over 5 years X 5 spawning dates per prediction period (+00Y, +25Y, +50Y) with, for predictions at +25Y and +50Y, a repetition of simulations per quantile (0.025, 0.16, 0.5, 0.67 and 0.975) to take into account for the most extreme variations in water mass temperatures predicted by the FOCI model for a given date.

Description: These data aim at evaluating the hypothesis of long-distance dispersal across the North Atlantic and the Equatorial Atlantic belt for the cold seep mussels Gigantidas childressi, G. mauritanicus, Bathymodiolus heckerae and B. boomerang. We combined mitochondrial Cox1 barcoding of some mussel specimens from both sides of the Atlantic (American vs European/African margins) with larval dispersal trajectories simulated from the VIKING20X model of the Atlantic circulation at a spatial scale not yet investigated. Larval dispersal modelling data correspond to transports of larvae over one year in surface waters from 21 geographic localities over 5 consecutive years (2015, 2016, 2017, 2018 and 2019) and 5 spawning dates (November, December, January, February and March) per year. Genetic data correspond to the geo-referenced sequences obtained for the 4 mussel species from some of the localities where larvae have been released during the modelling approach.

Description: This data set presents the results of an automated cluster analysis using Gaussian mixture models of the entire Atlantic seafloor environment. The analysis was based on eight global datasets and their derivatives: Bathymetry, slope, terrain ruggedness index, topographic position index, sediment thickness, POC flux, salinity, dissolved oxygen, temperature, current velocity, and phytoplankton abundance in surface waters along with seasonal variabilities (see Source data set). We obtained nine seabed areas (SBAs) that portray the Atlantic seafloor that are shown as polygons in the data set. The attribute table holds short descriptions of each SBA as well as about the colours used in the accompanying paper publication. Data sets like this can be used for further analysis like e.g. for landscape ecology metrics to identify regions of interest. The compressed file further contains a style file that can be used to directly load the correct style in the QGIS software package.

Description: Highlights: • OGCM simulations of the AMOC are highly sensitive to the subarctic freshwater forcing. • Trends in the simulated AMOC are linked to the salinity of the DSOW. • DSOW salinity trends can be traced back to the freshwater transport by the NAC. • The NAC freshwater budget is highly affected by the salinity restoring used in OGCMs. • Modifications in the subarctic precipitation can help to minimize the restoring flux. Global ocean sea-ice models with an atmospheric forcing based on bulk formulations of the air-sea fluxes exhibit spurious trends in key flow indices like the Atlantic Meridional Overturning Circulation (AMOC), constraining their use in investigations of multi-decadal ocean variability. To identify the critical model factors affecting the temporal evolution of the AMOC on time scales of up to 60 years, a series of experiments with both eddy-permitting (0.25°) and non-eddying (0.5°) ocean-ice models has been performed, focusing on the influence of artificial choices for the freshwater forcing, in particular the restoring of sea surface salinity towards climatological values. The atmospheric forcing builds on the proposal for Coordinated Ocean-ice Reference Experiments (CORE), utilizing the refined atmospheric reanalysis products for 1948–2006 compiled by Large and Yeager. Sensitivity experiments with small variations in precipitation (within the observational uncertainty) and sea surface salinity restoring in the subarctic Atlantic produce a wide range of AMOC transports, between upward drifts to more than 22 Sv and nearly-collapsed states with less than 7 Sv, reflecting the excessive role of the salinity feedback in such simulations. In all cases the AMOC is tightly related to the density of the Denmark Strait overflow; changes in that density are governed by the salinity in the Nordic Seas; and in turn, that salinity is strongly affected by the properties of the inflowing North Atlantic water.

Description: Highlights: • We combine high-resolution ocean models with population genetics • Variation in wind-driven ocean currents mediates the collapse of A. anguilla • Female eels are philopatric within the Sargasso Sea, while males maintain gene flow • We present first evidence of the role of ocean currents in shaping species’ evolution Summary: Worldwide, exploited marine fish stocks are under threat of collapse [1]. Although the drivers behind such collapses are diverse, it is becoming evident that failure to consider evolutionary processes in fisheries management can have drastic consequences on a species’ long-term viability [2]. The European eel (Anguilla anguilla; Linnaeus, 1758) is no exception: not only does the steep decline in recruitment observed in the 1980s [ 3 and 4] remain largely unexplained, the punctual detection of genetic structure also raises questions regarding the existence of a single panmictic population [ 5, 6 and 7]. With its extended Transatlantic dispersal, pinpointing the role of ocean dynamics is crucial to understand both the population structure and the widespread decline of this species. Hence, we combined dispersal simulations using a half century of high-resolution ocean model data with population genetics tools. We show that regional atmospherically driven ocean current variations in the Sargasso Sea were the major driver of the onset of the sharp decline in eel recruitment in the beginning of the 1980s. The simulations combined with genotyping of natural coastal eel populations furthermore suggest that unexpected evidence of coastal genetic differentiation is consistent with cryptic female philopatric behavior within the Sargasso Sea. Such results demonstrate the key constraint of the variable oceanic environment on the European eel population.

Description: Highlights: • Phase II of the Coordinated Ocean-ice Reference Experiments (CORE-II) is introduced. • Solutions from CORE-II simulations from eighteen participating models are presented. • Mean states in the North Atlantic with a focus on AMOC are examined. • The North Atlantic solutions differ substantially among the models. • Many factors, including parameterization choices, contribute to these differences. Simulation characteristics from eighteen global ocean–sea-ice coupled models are presented with a focus on the mean Atlantic meridional overturning circulation (AMOC) and other related fields in the North Atlantic. These experiments use inter-annually varying atmospheric forcing data sets for the 60-year period from 1948 to 2007 and are performed as contributions to the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II). The protocol for conducting such CORE-II experiments is summarized. Despite using the same atmospheric forcing, the solutions show significant differences. As most models also differ from available observations, biases in the Labrador Sea region in upper-ocean potential temperature and salinity distributions, mixed layer depths, and sea-ice cover are identified as contributors to differences in AMOC. These differences in the solutions do not suggest an obvious grouping of the models based on their ocean model lineage, their vertical coordinate representations, or surface salinity restoring strengths. Thus, the solution differences among the models are attributed primarily to use of different subgrid scale parameterizations and parameter choices as well as to differences in vertical and horizontal grid resolutions in the ocean models. Use of a wide variety of sea-ice models with diverse snow and sea-ice albedo treatments also contributes to these differences. Based on the diagnostics considered, the majority of the models appear suitable for use in studies involving the North Atlantic, but some models require dedicated development effort.

Description: Highlights: • Global mean sea level simulated in interannual CORE simulations. • Regional sea level patterns simulated in interannual CORE simulations. • Theoretical foundation for analysis of global mean sea level and regional patterns. Abstract: We provide an assessment of sea level simulated in a suite of global ocean-sea ice models using the interannual CORE atmospheric state to determine surface ocean boundary buoyancy and momentum fluxes. These CORE-II simulations are compared amongst themselves as well as to observation-based estimates. We focus on the final 15 years of the simulations (1993–2007), as this is a period where the CORE-II atmospheric state is well sampled, and it allows us to compare sea level related fields to both satellite and in situ analyses. The ensemble mean of the CORE-II simulations broadly agree with various global and regional observation-based analyses during this period, though with the global mean thermosteric sea level rise biased low relative to observation-based analyses. The simulations reveal a positive trend in dynamic sea level in the west Pacific and negative trend in the east, with this trend arising from wind shifts and regional changes in upper 700 m ocean heat content. The models also exhibit a thermosteric sea level rise in the subpolar North Atlantic associated with a transition around 1995/1996 of the North Atlantic Oscillation to its negative phase, and the advection of warm subtropical waters into the subpolar gyre. Sea level trends are predominantly associated with steric trends, with thermosteric effects generally far larger than halosteric effects, except in the Arctic and North Atlantic. There is a general anti-correlation between thermosteric and halosteric effects for much of the World Ocean, associated with density compensated changes.

Description: Estimating the magnitude of Agulhas leakage, the volume flux of water from the Indian to the Atlantic Ocean, is difficult because of the presence of other circulation systems in the Agulhas region. Indian Ocean water in the Atlantic Ocean is vigorously mixed and diluted in the Cape Basin. Eulerian integration methods, where the velocity field perpendicular to a section is integrated to yield a flux, have to be calibrated so that only the flux by Agulhas leakage is sampled. Two Eulerian methods for estimating the magnitude of Agulhas leakage are tested within a high-resolution two-way nested model with the goal to devise a mooring-based measurement strategy. At the GoodHope line, a section halfway through the Cape Basin, the integrated velocity perpendicular to that line is compared to the magnitude of Agulhas leakage as determined from the transport carried by numerical Lagrangian floats. In the first method, integration is limited to the flux of water warmer and more saline than specific threshold values. These threshold values are determined by maximizing the correlation with the float-determined time series. By using the threshold values, approximately half of the leakage can directly be measured. The total amount of Agulhas leakage can be estimated using a linear regression, within a 90% confidence band of 12 Sv. In the second method, a subregion of the GoodHope line is sought so that integration over that subregion yields an Eulerian flux as close to the float-determined leakage as possible. It appears that when integration is limited within the model to the upper 300 m of the water column within 900 km of the African coast the time series have the smallest root-mean-square difference. This method yields a root-mean-square error of only 5.2 Sv but the 90% confidence band of the estimate is 20 Sv. It is concluded that the optimum thermohaline threshold method leads to more accurate estimates even though the directly measured transport is a factor of two lower than the actual magnitude of Agulhas leakage in this model.

Description: Model drift in the Labrador Sea in eddy permitting model simulations is examined using a series of configurations based on the NEMO numerical framework. There are two phases of the drift that we can identify, beginning with an initial rapid 3-year period, associated with the adjustment of the model from its initial conditions followed by an extended model drift/adjustment that continued for at least another decade. The drift controlled the model salinity in the Labrador Sea, over-riding the variability. Thus, during this initial period, similar behavior was observed between the inter-annually forced experiments as with perpetual year forcing. The results also did not depend on whether the configuration was global, or regional North Atlantic Ocean. The inclusion of an explicit sea-ice component did not seem to have a significant impact on the interior drift. Clear cut evidence for the drift having an advective nature was shown, based on two separate currents/flow regimes. We find, as expected, the representation of freshwater in the sub-polar gyre’s boundary currents important. But this study also points out another, equally important process and pathway: the input of high salinity mode water from the subtropical North Atlantic. The advective regime is dependent on the details of the model, such as the representation of the freshwater transport in the model’s East Greenland Current being very sensitive to the strength of the local sea surface salinity restoring (and the underlying field that the model is being restored to).