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: 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: 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: A new Earth system model, the Flexible Ocean and Climate Infrastructure (FOCI), is introduced. A first version of FOCI consists of a global high-top atmosphere (ECHAM6.3) and an ocean model (NEMO3.6) as well as sea ice (LIM2) and land surface model components (JSBACH), which are coupled through the OASIS3-MCT software package. FOCI includes a number of optional modules which can be activated depending on the scientific question of interest. In the atmosphere, interactive stratospheric chemistry can be used (ECHAM6-HAMMOZ) to study, for example, the effects of the ozone hole on the climate system. In the ocean, a biogeochemistry model (MOPS) is available to study the global carbon cycle. A unique feature of FOCI is the ability to explicitly resolve mesoscale ocean eddies in specific regions. This is realized in the ocean through nesting; first examples for the Agulhas Current and the Gulf Stream systems are described here. FOCI therefore bridges the gap between coarse-resolution climate models and global high-resolution weather prediction and ocean-only models. It allows to study the evolution of the climate system on regional and seasonal to (multi-) decadal scales. The development of FOCI resulted from a combination of the long-standing expertise in ocean and climate modeling in several research units and divisions at GEOMAR. FOCI will thus be used to complement and interpret long-term observations in the Atlantic, enhance the process understanding of the role of mesoscale oceanic eddies for large-scale oceanic and atmospheric circulation patterns, study feedback mechanisms with stratospheric processes, estimate future ocean acidification, improve the simulation of the Atlantic Meridional Overturning Circulation changes and their influence on climate, ocean chemistry and biology. In this paper we present both the scientific vision for the development of FOCI as well as some technical details. This includes a first validation of the different model components using several configurations of FOCI. Results show that the model in its basic configuration runs stably under pre-industrial control as well as under historical forcing, and produces a mean climate and variability which compares well with observations, reanalysis products and other climate models. The nested configurations reduce some long-standing biases in climate models and are an important step forward to include the atmospheric response in multi-decadal eddy-rich configurations.

Description: The North Atlantic Deep Water (NADW) is a crucial component of the Atlantic meridional overturning circulation and is therefore an important factor of the climate system. In order to estimate the mean relative contributions, sources, and pathways of the NADW at the southern exit of the Labrador Sea, a Lagrangian particle experiment is performed. The particles were seeded according to the strength of the velocity field along the 53∘ N section and traced 40 years backward in time in the three-dimensional velocity and hydrography field. The resulting transport pathways, their sources and corresponding transit timescales were inferred. Our experiment shows that, of the 30.1 Sv of NADW passing 53∘ N on average, the majority of this water is associated with a diapycnal mass flux without contact to the atmosphere, accounting for 14.3 Sv (48 %), where 6.2 Sv originate from the Labrador Sea, compared to 4.7 Sv from the Irminger Sea. The second-largest contribution originates from the mixed layer with 7.2 Sv (24 %), where the Labrador Sea contribution (5.9 Sv) dominates over the Irminger Sea contribution (1.0 Sv). Another 5.7 Sv (19 %) of NADW crosses the Greenland–Scotland Ridge within the NADW density class, where about two-thirds pass the Denmark Strait, while one-third crosses the Iceland–Scotland Ridge. The NADW exported at 53∘ N is hence dominated by entrainment through the diapycnal mass flux and mixed-layer origin in the Labrador Sea.

Description: With the increase in computational power, ocean models with kilometer-scale resolution have emerged over the last decade. These models have been used for quantifying the energetic exchanges between spatial scales, informing the design of eddy parametrizations, and preparing observing networks. The increase in resolution, however, has drastically increased the size of model outputs, making it difficult to transfer and analyze the data. It remains, nonetheless, of primary importance to assess more systematically the realism of these models. Here, we showcase a cloud-based analysis framework proposed by the Pangeo project that aims to tackle such distribution and analysis challenges. We analyze the output of eight submesoscale-permitting simulations, all on the cloud, for a crossover region of the upcoming Surface Water and Ocean Topography (SWOT) altimeter mission near the Gulf Stream separation. The cloud-based analysis framework (i) minimizes the cost of duplicating and storing ghost copies of data and (ii) allows for seamless sharing of analysis results amongst collaborators. We describe the framework and provide example analyses (e.g., sea-surface height variability, submesoscale vertical buoyancy fluxes, and comparison to predictions from the mixed-layer instability parametrization). Basin- to global-scale, submesoscale-permitting models are still at their early stage of development; their cost and carbon footprints are also rather large. It would, therefore, benefit the community to document the different model configurations for future best practices. We also argue that an emphasis on data analysis strategies would be crucial for improving the models themselves.

Description: A hierarchy of global 1/4° (ORCA025) and Atlantic Ocean 1/20° nested (VIKING20X) ocean/sea-ice models is described. It is shown that the eddy-rich configurations performed in hindcasts of the past 50–60 years under CORE and JRA55-do atmospheric forcings realistically simulate the large-scale horizontal circulation, the distribution of the mesoscale, overflow and convective processes, and the representation of regional current systems in the North and South Atlantic. The representation, and in particular the long-term temporal evolution, of the Atlantic Meridional Overturning Circulation (AMOC) strongly depends on numerical choices for the application of freshwater fluxes. The interannual variability of the AMOC instead is highly correlated among the model experiments and also with observations, including the 2010 minimum observed by RAPID at 26.5° N pointing at a dominant role of the forcing. Regional observations in western boundary current systems at 53° N, 26.5° N and 11° S are explored in respect to their ability to represent the AMOC and to monitor the temporal evolution of the AMOC. Apart from the basin-scale measurements at 26.5° N, it is shown that in particular the outflow of North Atlantic Deepwater at 53° N is a good indicator of the subpolar AMOC trend during the recent decades, if the latter is provided in density coordinates. The good reproduction of observed AMOC and WBC trends in the most reasonable simulations indicate that the eddy-rich VIKING20X is capable in representing realistic forcing-related and ocean-intrinsic trends.

Description: Observations of the eastern subpolar North Atlantic in the 2010s show exceptional freshening and cooling of the upper ocean, peaking in 2016 with the lowest salinities recorded for 120 years. Published theories for the mechanisms driving the freshening include: reduced transport of saltier, warmer surface waters northwards from the subtropics associated with reduced meridional overturning; shifts in the pathways of fresher, cooler surface water from the Labrador Sea driven by changing patterns of wind stress; and the eastward expansion of the subpolar gyre. Using output from a high-resolution hindcast model simulation, we propose that the primary cause of the exceptional freshening and cooling is reduced surface heat loss in the Labrador Sea. Tracking virtual fluid particles in the model backwards from the eastern subpolar North Atlantic between 1990 and 2020 shows the major cause of the freshening and cooling to be an increased outflow of relatively fresh and cold surface waters from the Labrador Sea; with a minor contribution from reduced transport of warmer, saltier surface water northward from the subtropics. The cooling, but not the freshening, produced by these changing proportions of waters of subpolar and subtropical origin is mitigated by reduced along-track heat loss to the atmosphere in the North Atlantic Current. We analyse modelled boundary exchanges and water mass transformation in the Labrador Sea to show that since 2000, while inflows of lighter surface waters remain steady, the increasing output of these waters is due to reduced surface heat loss in the Labrador Sea beginning in the early 2000s. Tracking particles further upstream reveals that the primary source of the increased volume of lighter water transported out of the Labrador Sea is increased recirculation of water, and therefore longer residence times, in the upper 500–1000 m of the subpolar gyre.

Description: The inflow of relatively warm and salty water from the Indian Ocean into the South Atlantic via Agulhas leakage is important for the global overturning circulation and the global climate. In this study, we analyse the robustness of Agulhas leakage estimates as well as the thermohaline property modifications of Agulhas leakage south of Africa. Lagrangian experiments with both the newly developed tool Parcels and the well established tool Ariane were performed to simulate Agulhas leakage in the eddy-rich ocean–sea-ice model INALT20 (1/20∘ horizontal resolution) forced by the JRA55-do atmospheric boundary conditions. The average transport, its variability, trend and the transit time from the Agulhas Current to the Cape Basin of Agulhas leakage is simulated comparably with both Lagrangian tools, emphasizing the robustness of our method. Different designs of the Lagrangian experiment alter in particular the total transport of Agulhas leakage by up to 2 Sv, but the variability and trend of the transport are similar across these estimates. During the transit from the Agulhas Current at 32∘ S to the Cape Basin, a cooling and freshening of Agulhas leakage waters occurs especially at the location of the Agulhas Retroflection, resulting in a density increase as the thermal effect dominates. Beyond the strong air–sea exchange around South Africa, Agulhas leakage warms and salinifies the water masses below the thermocline in the South Atlantic.