Description: Rapid increases in upper 700‐m Indian Ocean heat content (IOHC) since the 2000s have focused attention on its role during the recent global surface warming hiatus. Here, we use ocean model simulations to assess distinct multidecadal IOHC variations since the 1960s and explore the relative contributions from wind stress and buoyancy forcing regionally and with depth. Multidecadal wind forcing counteracted IOHC increases due to buoyancy forcing from the 1960s to the 1990s. Wind and buoyancy forcing contribute positively since the mid‐2000s, accounting for the drastic IOHC change. Distinct timing and structure of upper ocean temperature changes in the eastern and western Indian Ocean are linked to the pathway how multidecadal wind forcing associated with the Interdecadal Pacific Oscillation is transmitted and affects IOHC through local and remote winds. Progressive shoaling of the equatorial thermocline—of importance for low‐frequency variations in Indian Ocean Dipole occurrence—appears to be dominated by multidecadal variations in wind forcing.

Description: Mesoscale eddies can be strengthened by the absorption of submesoscale eddies resulting from mixed-layer baroclinic instabilities. This is shown for mesoscale eddies in the Agulhas Current system by investigating the kinetic energy cascade with a spectral and a coarse-graining approach in two model simulations of the Agulhas region. One simulation resolves mixed-layer baroclinic instabilities and one does not. When mixed-layer baroclinic instabilities are included, the largest submesoscale near-surface fluxes occur in winter-time in regions of strong mesoscale activity for upscale as well as downscale directions. The forward cascade at the smallest resolved scales occurs mainly in frontogenetic regions in the upper 30 m of the water column. In the Agulhas ring path, the forward cascade changes to an inverse cascade at a typical scale of mixed-layer eddies (15 km). At the same scale, the largest sources of the upscale flux occur. After the winter, the maximum of the upscale flux shifts to larger scales. Depending on the region, the kinetic energy reaches the mesoscales in spring or early summer aligned with the maximum of mesoscale kinetic energy. This indicates the importance of submesoscale flows for the mesoscale seasonal cycle. A case study shows that the underlying process is the mesoscale absorption of mixed-layer eddies. When mixed-layer baroclinic instabilities are not included in the simulation, the open-ocean upscale cascade in the Agulhas ring path is almost absent. This contributes to a 20 %-reduction of surface kinetic energy at mesoscales larger than 100 km when submesoscale dynamics are not resolved by the model.

Description: Warming of the North Atlantic region in climate history often was associated with massive melting of the Greenland Ice Sheet. To identify the meltwater’s impacts and isolate these from internal variability and other global warming factors, we run single-forcing simulations including small ensembles using three complex climate models differing only in their ocean components. In 200-year long pre-industrial climate simulations, we identify robust consequences of abruptly increasing Greenland runoff by 0.05 Sv: sea-level rise of 44±10 cm, subpolar North Atlantic surface cooling of 0.7˚C and a moderate AMOC decline of 1.1–2.0 Sv. The latter two emerge in under three decades—and reverse on the same timescale after the perturbation ends in year 100. The ocean translates the step-change perturbation into a multi-decadal to centennial signature in the deep overturning circulation. In all simulations, internal variability creates notable uncertainty in estimating trends, time of emergence and duration of the response.

Description: In highly fragmented and relatively stable cold-seep ecosystems, species are expected to exhibit high migration rates and long-distance dispersal of long-lived pelagic larvae to maintain genetic integrity over their range. Accordingly, several species inhabiting cold seeps are widely distributed across the whole Atlantic Ocean, with low genetic divergence between metapopulations on both sides of the Atlantic Equatorial Belt (AEB, i.e. Barbados and African/European margins). Two hypotheses may explain such patterns: (i) the occurrence of present-day gene flow or (ii) incomplete lineage sorting due to large population sizes and low mutation rates. Here, we evaluated the first hypothesis using the cold seep mussels Gigantidas childressi, G. mauritanicus, Bathymodiolus heckerae and B. boomerang. We combined COI barcoding of 763 individuals with VIKING20X larval dispersal modelling at a large spatial scale not previously investigated. Population genetics supported the parallel evolution of Gigantidas and Bathymodiolus genera in the Atlantic Ocean and the occurrence of a 1-3 Million-year-old vicariance effect that isolated populations across the Caribbean Sea. Both population genetics and larval dispersal modelling suggested that contemporary gene flow and larval exchanges are possible across the AEB and the Caribbean Sea, although probably rare. When occurring, larval flow was eastward (AEB - only for B. boomerang) or northward (Caribbean Sea - only for G. mauritanicus). Caution is nevertheless required since we focused on only one mitochondrial gene, which may underestimate gene flow if a genetic barrier exists. Non-negligible genetic differentiation occurred between Barbados and African populations, so we could not discount the incomplete lineage sorting hypothesis. Larval dispersal modelling simulations supported the genetic findings along the American coast with high amounts of larval flow between the Gulf of Mexico (GoM) and the US Atlantic Margin, although the Blake Ridge population of B. heckerae appeared genetically differentiated. Overall, our results suggest that additional studies using nuclear genetic markers and population genomics approaches are needed to clarify the evolutionary history of the Atlantic bathymodioline mussels and to distinguish between ongoing and past processes.

Description: The Makassar Strait, the main passageway of the Indonesian Throughflow (ITF), is an important component of Indo-Pacific climate through its inter-basin redistribution of heat and freshwater. Observational studies suggest that wind-driven freshwater advection from the marginal seas into the Makassar Strait modulates the strait's surface transport. However, direct observations are too short (〈15 years) to resolve variability on decadal timescales. Here we use a series of global ocean simulations to assess the advected freshwater contributions to ITF transport across a range of timescales. The simulated seasonal and interannual freshwater dynamics are consistent with previous studies. On decadal timescales, we find that wind-driven advection of South China Sea (SCS) waters into the Makassar Strait modulates upper-ocean ITF transport. Atmospheric circulation changes associated with Pacific decadal variability appear to drive this mechanism via Pacific lower-latitude western boundary current interactions that affect the SCS circulation. Key Points: - A global ocean model is used to show how freshwater impacts the decadal variability of transport through the main Indonesian Throughflow pathway - Wind-driven advection of South China Sea freshwater induces an upstream pressure gradient that reduces transport - Freshwater input is modulated by atmospheric circulation changes associated with Pacific decadal variability

Description: Deep convection and associated deep water formation are key processes for climate variability, since they impact the oceanic uptake of heat and trace gases and alter the structure and strength of the global overturning circulation. For long, deep convection in the subpolar North Atlantic was thought to be confined to the central Labrador Sea in the western subpolar gyre (SPG). However, there is increasing observational evidence that deep convection also has occurred in the eastern SPG south of Cape Farewell and in the Irminger Sea, in particular, in 2015–2018. Here we assess this recent event in the context of the temporal evolution of spatial deep convection patterns in the SPG since the mid-twentieth century, using realistic eddy-rich ocean model simulations. These reveal a large interannual variability with changing contributions of the eastern SPG to the total deep convection volume. Notably, in the late 1980s to early 1990s, the period with highest deep convection intensity in the Labrador Sea related to a persistent positive phase of the North Atlantic Oscillation, the relative contribution of the eastern SPG was small. In contrast, in 2015–2018, deep convection occurred with an unprecedented large relative contribution of the eastern SPG. This is partly linked to a smaller north-westward extent of deep convection in the Labrador Sea compared to previous periods of intensified deep convection, and may be a first fingerprint of freshening trends in the Labrador Sea potentially associated with enhanced Greenland melting and the oceanic advection of the 2012–2016 eastern North Atlantic fresh anomaly.

Description: Landscape maps based on multivariate cluster analyses provide an objective and comprehensive view on the (marine) environment. They can hence support decision making regarding sustainable ocean resource handling and protection schemes. Across a large number of scales, input parameters and classification methods, numerous studies categorize the ocean into seascapes, hydro-morphological provinces or clusters. Many of them are regional, however, while only a few are on a basin scale. This study presents an automated cluster analysis of the entire Atlantic seafloor environment, 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. As a result, we obtained nine seabed areas (SBAs) that portray the Atlantic seafloor. Some SBAs have a clear geological and geomorphological nature, while others are defined by a mixture of terrain and water body characteristics. The majority of the SBAs, especially those covering the deep ocean areas, are coherent and show little seasonal and hydrographic variation, whereas other, nearshore SBAs, are smaller sized and dominated by high seasonal changes. To demonstrate the potential use of the marine landscape map for marine spatial planning purposes, we mapped out local SBA diversity using the patch richness index developed in landscape ecology. It identifies areas of high landscape diversity, and is a practical way of defining potential areas of interest, e.g. for designation as protected areas, or for further research. Clustering probabilities are highest (100%) in the center of SBA patches and decrease towards the edges (〈 98%). On the SBA point cloud which was reduced for probabilities 〈98%, we ran a diversity analysis to identify and highlight regions that have a high number of different SBAs per area, indicating the use of such analyses to automatically find potentially delicate areas. We found that some of the highlights are already within existing EBSAs, but the majority is yet unexplored.

Description: Earth System Sciences have been generating increasingly larger amounts of heterogeneous data in recent years. We identify the need to combine Earth System Sciences with Data Sciences, and give our perspective on how this could be accomplished within the sub-field of Marine Sciences. Marine data hold abundant information and insights that Data Science techniques can reveal. There is high demand and potential to combine skills and knowledge from Marine and Data Sciences to best take advantage of the vast amount of marine data. This can be accomplished by establishing Marine Data Science as a new research discipline. Marine Data Science is an interface science that applies Data Science tools to extract information, knowledge, and insights from the exponentially increasing body of marine data. Marine Data Scientists need to be trained Data Scientists with a broad basic understanding of Marine Sciences and expertise in knowledge transfer. Marine Data Science doctoral researchers need targeted training for these specific skills, a crucial component of which is co-supervision from both parental sciences. They also might face challenges of scientific recognition and lack of an established academic career path. In this paper, we, Marine and Data Scientists at different stages of their academic career, present perspectives to define Marine Data Science as a distinct discipline. We draw on experiences of a Doctoral Research School, MarDATA, dedicated to training a cohort of early career Marine Data Scientists. We characterize the methods of Marine Data Science as a toolbox including skills from their two parental sciences. All of these aim to analyze and interpret marine data, which build the foundation of Marine Data Science.

Description: While forced ocean hindcast simulations are useful for a wide range of applications, a key limitation is their inability to simulate ocean-atmosphere feedbacks. As a consequence, they need to rely on artificial choices such as sea surface salinity restoring and other corrections affecting the surface freshwater fluxes. Fully coupled models overcome these limitations, but lack the correct timing of variability due to weaker observational constraints. This leads to a mismatch between forced and coupled models on interannual to decadal timescales. A possibility to combine the advantages of both modelling strategies is to apply a partial coupling (PCPL), i.e. replacing the surface winds stress in the ocean component by wind stress derived from reanalysis. To identify the capabilities, limitations and possible use cases of partial coupling, we perform a fully coupled, two partially coupled and an ocean-only experiment using an all-Atlantic nested ocean configuration at eddying resolution in a global climate model. We show that the correct timing of Atlantic Meridional Overturning Circulation (AMOC) variability in PCPL experiments is robust on timescales below 5 years. Mid-latitude wind stress curl changes contribute to decadal AMOC variability, but North Atlantic buoyancy fluxes are not significantly altered by incorporating reanalysed wind stress anomalies, limiting the success of PCPL on this timescale. Long term trends of the AMOC in PCPL mode are consistent with fully coupled model experiments under historic atmospheric boundary conditions, suggesting that a partially coupled model is still able to simulate the important ocean-atmosphere feedbacks necessary to maintain a stable AMOC.

Description: Marine Heatwaves (MHWs) are ocean extreme events, characterized by anomalously high temperatures, which can have significant ecological impacts. The Northeast U.S. continental shelf is of great economical importance as it is home to a highly productive ecosystem. Local warming rates exceed the global average and the region experienced multiple MHWs in the last decade with severe consequences for regional fisheries. Due to the lack of subsurface observations, the depth-extent of MHWs is not well-known, which hampers the assessment of impacts on pelagic and benthic ecosystems. This study utilizes a global ocean circulation model with a high-resolution (1/20°) nest in the Atlantic to investigate the depth structure of MHWs and associated drivers on the Northeast U.S. continental shelf. It is shown that MHWs exhibit varying spatial extents, with some only occurring at depth. The highest intensities are found around 100 m depth with temperatures exceeding the climatological mean by up to 7°C, while surface intensities are typically smaller (around 3°C). Distinct vertical structures are associated with different spatial MHW patterns and drivers. Investigation of the co-variability of temperature and salinity reveals that over 80% of MHWs at depth (〉50 m) coincide with extreme salinity anomalies. Two case studies provide insight into opposing MHW patterns at the surface and at depth, being forced by anomalous air-sea heat fluxes and Gulf Stream warm core ring interaction, respectively. The results highlight the importance of local ocean dynamics and the need to realistically represent them in climate models.