ALBERT

Description: A coordinated effort involving trailblazing science — and icebreaking ships — from many nations is needed to fill gaps in our understanding of the Arctic Ocean and how it’s changing.

Description: We investigate the role of the tropics, the stratosphere, and atmosphere‐ocean coupling for seasonal forecasts of strong, potentially damaging, Northern Hemisphere extratropical winter wind storm frequencies. This is done by means of relaxation experiments with the European Centre for Medium‐Range Weather Forecasts model, which allow us to prescribe perfect forecasts for specific parts of the coupled atmosphere‐ocean system. We find that perfect predictions of the Northern Hemisphere stratosphere significantly enhance winter storm predictive skill between eastern Greenland and Northern Europe. Correct seasonal predictions of the occurrence of stratospheric sudden warmings play a decisive role. The importance of correctly predicting the tropics and of two‐way atmosphere‐ocean coupling, both for forecasting stratospheric sudden warming risk and, correspondingly, severe winter storm frequency, is noted.

Description: The northern part of the South China Sea is characterized by widespread occurrence of bottom simulating reflectors (BSR) indicating the presence of marine gas hydrate. Because the area covers both a tectonically inactive passive margin and the termination of a subduction zone, the influence of tectonism on the dynamics of gas hydrate systems can be studied in this region. Geophysical data show that there are multiple thrust faults on the active margin while much fewer and smaller faults exist in the passive margin. This tectonic difference matches with a difference in the geophysical characteristics of the gas hydrate systems. High hydrate saturation derived from ocean bottom seismometer data and controlled source electromagnetic data and conspicuous high‐amplitude reflections in P‐Cable 3D seismic data above the BSR are found in the anticlinal ridges of the active margin. In contrast all geophysical evidence for the passive margin points to normal to low hydrate saturations. Geochemical analyses of gas samples collected at seep sites on the active margin show methane with heavy δ13C isotope composition, while gas collected at the passive margin shows light carbon isotope composition. Thus, we interpret the passive margin as a typical gas hydrate province fuelled by biogenic production of methane and the active margin gas hydrate system as a system that is fuelled not only by biogenic gas production but also by additional advection of thermogenic methane from the subduction system.

Description: The mechanisms controlling the variability of oxygen levels in the ocean are poorly quantified. We focus here on the impact of wind synoptic variability associated with tropical convective regions and extra‐tropical storms. Removing the wind higher frequencies of variability (2 days – 1 month) in an atmosphere reanalysis used to force an ocean model decreases wind stress by up to 20% in the tropics and 50% in the mid‐latitudes, weakening wind‐driven ocean circulation by 20%. Oxygen levels decrease by up to 10 mmol.m‐3 in tropical oceans and 30 mmol.m‐3 in subtropical gyres mainly due to changes in advective processes. While a large part of the tropical oxygen anomaly has local origins, changes in oxygen levels in the subtropical gyres modulate tropical oxygen distribution. Our study suggests that the “storminess” of the ocean is an important parameter that could determine the future evolution of poorly oxygenated regions.

Description: The Benguela Upwelling system (BUS) is the most productive of all eastern boundary upwelling ecosystems and it hosts a well‐developed oxygen minimum zone. As such, the BUS is a potential hotspot for production of N2O, a potent greenhouse gas derived from microbially‐driven decay of sinking organic matter. Yet, the extent at which near‐surface waters emit N2O to the atmosphere in the BUS is highly uncertain. Here we present the first high‐resolution surface measurements of N2O across the northern part of the BUS (nBUS). We found strong gradients with a three‐fold increase in N2O concentrations near the coast as compared with open ocean waters. Our observations show enhanced sea‐to‐air fluxes of N2O (up to 1.67 nmol m−2 s−1) in association with local upwelling cells. Based on our data we suggest that the nBUS can account for 13% of the total coastal upwelling source of N2O to the atmosphere.

Description: Hydroxylamine (NH 2 OH), a short-lived intermediate in the nitrogen cycle, is a potential precursor of nitrous oxide (N 2 O) in the ocean. However, measurements of NH 2 OH in the ocean are sparse. Here we present a data set of depth profiles of NH 2 OH from the equatorial Atlantic Ocean and the eastern tropical South Pacific and compare it to N 2 O, nitrate, and nitrite profiles under varying oxygen conditions. The presence of NH 2 OH in surface waters points toward surface nitrification in the upper 100 m. Overall, we found a ratio of 1:3 between NH 2 OH and N 2 O in open ocean areas when oxygen concentrations were 〉50 μmol/L. In the equatorial Atlantic Ocean and the open ocean eastern tropical South Pacific, where nitrification is the dominant N 2 O production pathway, stepwise multiple regressions demonstrated that N 2 O, NH 2 OH, and nitrate concentrations were highly correlated, suggesting that NH 2 OH is a potential indicator for nitrification.

Description: Radiative forcing from volcanic aerosol impacts surface temperatures; however, the background climate state also affects the response. A key question thus concerns whether constraining forcing estimates is more important than constraining initial conditions for accurate simulation and attribution of posteruption climate anomalies. Here we test whether different realistic volcanic forcing magnitudes for the 1815 Tambora eruption yield distinguishable ensemble surface temperature responses. We perform a cluster analysis on a superensemble of climate simulations including three 30-member ensembles using the same set of initial conditions but different volcanic forcings based on uncertainty estimates. Results clarify how forcing uncertainties can overwhelm initial-condition spread in boreal summer due to strong direct radiative impact, while the effect of initial conditions predominate in winter, when dynamics contribute to large ensemble spread. In our setup, current uncertainties affecting reconstruction-simulation comparisons prevent conclusions about the magnitude of the Tambora eruption and its relation to the “year without summer.”

Description: Prediction and Research Moored Array in the Tropical Atlantic (PIRATA) is a multinational program initiated in 1997 in the tropical Atlantic to improve our understanding and ability to predict ocean-atmosphere variability. PIRATA consists of a network of moored buoys providing meteorological and oceanographic data transmitted in real time to address fundamental scientific questions as well as societal needs. The network is maintained through dedicated yearly cruises, which allow for extensive complementary shipboard measurements and provide platforms for deployment of other components of the Tropical Atlantic Observing System. This paper describes network enhancements, scientific accomplishments and successes obtained from the last 10 years of observations, and additional results enabled by cooperation with other national and international programs. Capacity building activities and the role of PIRATA in a future Tropical Atlantic Observing System that is presently being optimized are also described.

Description: Seafloor massive sulphide (SMS) deposits are of increasing economic interest in order to satisfy the relentless growth in worldwide metal demand. The Trans‐Atlantic Geotraverse (TAG) hydrothermal field at 26°N on the Mid‐Atlantic Ridge hosts several such deposits. This study presents new controlled source electromagnetic (CSEM), bathymetric and magnetic results from the TAG field. Potential SMS targets were selected based on their surface expressions in high‐resolution bathymetric data. High‐resolution reduced‐to‐the‐pole magnetic data show negative anomalies beneath and surrounding the SMS deposits, revealing large areas of hydrothermal alteration. CSEM data, sensitive to the electrical conductivity of SMS mineralization, further reveal a maximum thickness of up to 80 m and conductivities of up to 5 S/m. SMS samples have conductivities of up to a few thousand S/m, suggesting that remotely inferred conductivities represent an average of metal sulphide ores combined with silicified and altered host basalt that likely dominates at greater depths.

Description: The climate of the last two millennia was characterised by decadal to multi‐centennial variations which were recorded in terrestrial records and had important societal impacts. The cause of these climatic events is still under debate but changes in the North Atlantic circulation have often been proposed to play an important role. In this review we compile available high‐resolution paleoceanographic datasets from the northern North Atlantic and Nordic Seas. The records are grouped into regions related to modern ocean conditions and their variability is discussed. We additionally discuss our current knowledge from modelling studies, with a specific focus on the dynamical changes that are not well inferred from the proxy records. An illustration is provided through the analysis of two climate model ensembles and an individual simulation of the last millennium. This review thereby provides an up‐to‐date paleo‐perspective on the North Atlantic multidecadal to multi‐centennial ocean variability across the last two millennia.

Description: Key Points: • A regional ocean model is used to examine multidecadal shelf temperature changes on the Agulhas Bank • There are distinct shelf temperature regime changes in 1966 and 1996 • These regime shifts are caused by changes in coastal upwelling linked to large-scale wind variability The Agulhas Bank is an important area for the spawning of small pelagic fish and other species. Here, within a NEMO ocean model, we investigate changes in temperature over the Bank on multidecadal time scales. In agreement with previous observational studies, a shift to colder temperatures is found in 1997. The model also simulates an earlier shift from colder to warmer temperatures in 1966. These shifts are coastally confined and shown, using a climatologically forced model run as a control, to be driven by a north‐south migration in the large‐scale wind belts, rather than by changes in downward heat fluxes or changes in the Agulhas Current itself. The zonal wind changes on the Agulhas Bank show a significant relationship with the Southern Annular Mode, showing some promise for future predictability of cold and warm regimes on the Agulhas Bank. Thus, while the Agulhas Current has been shown in previous work to have a large impact on intra‐annual and interannual temperature variability, this work shows that multidecadal variability in temperature on the shelf is likely to be wind forced.

Description: Gravitational sinking of photosynthetically fixed particulate organic carbon (POC) constitutes a key component of the biological carbon pump. The fraction of POC leaving the surface ocean depends on POC sinking velocity (SV) and remineralization rate (Cremin), both of which depend on plankton community structure. However, the key drivers in plankton communities controlling SV and Cremin are poorly constrained. In fall 2014, we conducted a 6 weeks mesocosm experiment in the subtropical NE Atlantic Ocean to study the influence of plankton community structure on SV and Cremin. Oligotrophic conditions prevailed for the first 3 weeks, until nutrient‐rich deep water injected into all mesocosms stimulated diatom blooms. SV declined steadily over the course of the experiment due to decreasing CaCO3 ballast and – according to an optical proxy proposed herein – due to increasing aggregate porosity mostly during an aggregation event after the diatom bloom. Furthermore, SV was positively correlated with the contribution of picophytoplankton to the total phytoplankton biomass. Cremin was highest during a Synechococcus bloom under oligotrophic conditions and in some mesocosms during the diatom bloom after the deep‐water addition while it was particularly low during harmful algal blooms. The temporal changes were considerably larger in Cremin (max. 15‐fold) than in SV (max. 3‐fold). Accordingly, estimated POC transfer efficiency to 1000 m was mainly dependent on how the plankton community structure affected Cremin. Our approach revealed key players and interactions in the plankton food web influencing POC export efficiency thereby improving our mechanistic understanding of the biological carbon pump. Key points Sinking velocity was higher during oligotrophy than during blooms which is linked to ballast, porosity, and phytoplankton size structure Remineralization was highly variable but tended to be higher during Synechococcus or diatom blooms and lower during harmful algal blooms Plankton community structure had a considerably larger influence on particle remineralization rate than on sinking velocity

Description: Mid‐ocean ridge axes are marked by segmentation of the axes and underlying magmatic systems. Fine‐scale segmentation has mainly been studied along fast‐spreading ridges. Here we offer insight into the third‐ and fourth‐order segmentation of intermediate‐spreading ridges and their temporal evolution. The Alarcón Rise and the Endeavour Segment have similar spreading rates (49 and 52.5 mm/year, respectively) but contrasting morphologies that vary from an axial high with a relatively narrow axial summit trough to an axial valley. One‐meter resolution bathymetry acquired by autonomous underwater vehicles, lava geochemistry, and ages from sediment cores is combined with available seismic reflection profiles to analyze variations in (1) geometry and orientation of the axial summit trough or valley, (2) seafloor depth near the axis, and (3) distribution of hydrothermal vents, (4) lava chemistry, and (5) flow ages between contiguous axes. Along both intermediate‐spreading segments, third‐ and fourth‐order discontinuities and associated segments are similar in dimension to what has been observed along fast‐spreading ridges. The Alarcón Rise and the Endeavour Segment also allow the study of the evolution of fine‐scale segmentation over periods of 300 to 4,000 years. Comparison between old and young axes reveals that the evolution of fine‐scale segmentation depends on the intensity of the magmatic activity. High magmatic periods are associated with rapid evolution of third‐order segments, while low magmatic activity periods, dominated by tectonic deformation and/or hydrothermal activity, are associated with little to no change in segmentation.

Description: We study the structure and tectonics of the collision zone between the Nazca Ridge (NR) and the Peruvian margin constrained by seismic, gravimetric, bathymetric, and natural seismological data. The NR was formed in an on-ridge setting, and it is characterized by a smooth and broad shallow seafloor (swell) with an estimated buoyancy flux of ~7 Mg/s. The seismic results show that the NR hosts an oceanic lower crust 10–14 km thick with velocities of 7.2–7.5 km/s suggesting intrusion of magmatic material from the hot spot plume to the oceanic plate. Our results show evidence for subduction erosion in the frontal part of the margin likely enhanced by the collision of the NR. The ridge-trench collision zone correlates with the presence of a prominent normal scarp, a narrow continental slope, and (uplifted) shelf. In contrast, adjacent of the collision zone, the slope does not present a topographic scarp and the continental slope and shelf become wider and deeper. Geophysical and geodetic evidence indicate that the collision zone is characterized by low seismic coupling at the plate interface. This is consistent with vigorous subduction erosion enhanced by the subducting NR causing abrasion and increase of fluid pore pressure at the interplate contact. Furthermore, the NR has behaved as a barrier for rupture propagation of megathrust earthquakes (e.g., 1746 Mw 8.6 and 1942 Mw 8.1 events). In contrast, for moderate earthquakes (e.g., 1996 Mw 7.7 and 2011 Mw 6.9 events), the NR has behaved as a seismic asperity nucleating at depths 〉20 km.

Description: The variability of the Atlantic Meridional Overturning Circulation (AMOC) and its governing processes during the Last Glacial Maximum (LGM) is investigated in the Kiel Climate Model (KCM). Under LGM conditions, multidecadal AMOC variability is mainly forced by the surface heat flux variability linked to the East Atlantic pattern (EAP). In contrast, the multidecadal AMOC variability under preindustrial conditions is mainly driven by the surface heat flux variability associated with the North Atlantic Oscillation (NAO). Stand‐alone atmosphere model experiments show that relative to preindustrial conditions, the change in AMOC forcing under LGM conditions is tightly linked to the differences in topography. Key Points Multidecadal AMOC variability during the LGM and its associated physical processes have been investigated by means of a climate model Multidecadal AMOC variability during the LGM is mainly driven by surface heat flux variability linked to the East Atlantic pattern as opposed to the North Atlantic Oscillation under preindustrial conditions Change in topography during the LGM is responsible for the change in AMOC forcing

Description: Dissolved organic carbon (DOC) plays critical roles in marine carbon cycling, but its sources and sinks remain uncertain. In this study, we monitored DOC exudation rates of Sargassum natans under visible light (lambda 〉 390 nm) and solar radiation. DOC release rates ranged from 7 to 10 mu g C g(biomass)(-1) hr(-1) (wet weight) under visible light, but increased to 23 to 41 mu g C g(biomass)(-1) hr(-1) when exposed to natural sunlight. Results indicate that DOC released by Sargassum could amount to 0.3 to 1.2 Tg C/year, potentially contributing significantly to the marine DOC pool in the Gulf of Mexico and Western North Atlantic. We employed the Folin-Ciocalteu phenolic content method, nuclear magnetic resonance (NMR) spectroscopy, and ultrahigh resolution Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS) to characterize the diverse pool of organic compounds exuded from Sargassum. Results from these complementary methods showed that Sargassum release large quantities of phlorotannins, a class of polyphenols that have very similar properties to terrestrial DOC. These phlorotannins and their oxygenated phenolic derivatives exhibit a high hydrogen deficiency and functionalization (i.e., 4 to 6 oxygen atoms per aromatic ring), representing 5 to 18% of the released DOC isolated by solid phase extraction. Thus, Sargassum is the largest biological source of open ocean polyphenols recorded to date. The amount of polyphenolic DOC released by Sargassum challenges previous beliefs that all polyphenols found within the oceans are remnants of terrestrial organic matter, although the stability of phlorotannins and their derivatives needs to be further evaluated.

Description: We present a three-dimensional gas hydrate systems model of the southern Hikurangi subduction margin in eastern New Zealand. The model integrates thermal and microbial gas generation, migration, and hydrate formation. Modeling these processes has improved the understanding of factors controlling hydrate distribution. Three spatial trends of concentrated hydrate occurrence are predicted. The first trend (I) is aligned with the principal deformation front in the overriding Australian plate. Concentrated hydrate deposits are predicted at or near the apexes of anticlines and to be mainly sourced from focused migration and recycling of microbial gas generated beneath the hydrate stability zone. A second predicted trend (II) is related to deformation in the subducting Pacific plate associated with former Mesozoic subduction beneath Gondwana and the modern Pacific-Australian plate boundary. This trend is enhanced by increased advection of thermogenic gas through permeable layers in the subducting plate and focused migration into the Neogene basin fill above Cretaceous-Paleogene structures. The third trend (III) follows the northern margin of the Hikurangi Channel and is related to the presence of buried strata of the Hikurangi Channel system. The predicted trends are consistent with pronounced seismic reflection anomalies related to free gas in the pore space and strength of the bottom-simulating reflection. However, only trend I is also associated with clear and widespread seismic indications of concentrated gas hydrate. Total predicted hydrate masses at the southern Hikurangi Margin are between 52,800 and 69,800 Mt. This equates to 3.4–4.5 Mt hydrate/km2, containing 6.33 × 108–8.38 × 108 m3/km2 of methane.

Description: Global climatic changes during the last Glacial and Deglacial have been related to variations of the Atlantic Meridional Overturning Circulation (AMOC). Here, we present new and refined 231Pa/230Th down‐core profiles extending back to 30 ka BP from the northwestern Atlantic along the Atlantic Deep Western Boundary Current (DWBC), which is the main component of the southward deep backflow of the AMOC. Besides the well‐known Bermuda Rise records, available high‐resolution 231Pa/230Th data in the northwestern Atlantic are still sparse. Our new records along with reconstructions of deep water provenance from Nd isotopes constrain the timing and magnitude of past changes in AMOC from an additional northwestern Atlantic region forming a depth transect between 3000 and 4760 m water depth. Our extended and improved dataset confirms the weakening of the AMOC during deglacial cold spells such as Heinrich Event 1 and the Younger Dryas interrupted by a reinvigoration during the Bølling‐Allerød interstadial as seen in the prominent 231Pa/230Th records from the Bermuda Rise. However, in contrast to the Bermuda Rise records we find a clearly reduced circulation strength during the Last Glacial Maximum in the deep Atlantic.

Description: Gelatinous zooplankton (Cnidaria, Ctenophora, and Urochordata, namely, Thaliacea) are ubiquitous members of plankton communities linking primary production to higher trophic levels and the deep ocean by serving as food and transferring “jelly‐carbon” (jelly‐C) upon bloom collapse. Global biomass within the upper 200 m reaches 0.038 Pg C, which, with a 2–12 months life span, serves as the lower limit for annual jelly‐C production. Using over 90,000 data points from 1934 to 2011 from the Jellyfish Database Initiative as an indication of global biomass (JeDI: http://jedi.nceas.ucsb.edu, http://www.bco‐dmo.org/dataset/526852), upper ocean jelly‐C biomass and production estimates, organism vertical migration, jelly‐C sinking rates, and water column temperature profiles from GLODAPv2, we quantitatively estimate jelly‐C transfer efficiency based on Longhurst Provinces. From the upper 200 m production estimate of 0.038 Pg C year−1, 59–72% reaches 500 m, 46–54% reaches 1,000 m, 43–48% reaches 2,000 m, 32–40% reaches 3,000 m, and 25–33% reaches 4,500 m. This translates into ~0.03, 0.02, 0.01, and 0.01 Pg C year−1, transferred down to 500, 1,000, 2,000, and 4,500 m, respectively. Jelly‐C fluxes and transfer efficiencies can occasionally exceed phytodetrital‐based sediment trap estimates in localized open ocean and continental shelves areas under large gelatinous blooms or jelly‐C mass deposition events, but this remains ephemeral and transient in nature. This transfer of fast and permanently exported carbon reaching the ocean interior via jelly‐C constitutes an important component of the global biological soft‐tissue pump, and should be addressed in ocean biogeochemical models, in particular, at the local and regional scale.

Description: Sea surface salinity (SSS) is an important variable in the global ocean circulation. However, decadal to interdecadal changes in SSS are not well understood due to the lack of instrumental data. Here we reconstruct SSS from a paired, bimonthly resolved coral δ18O and Sr/Ca record from La Reunion Island that extends from 1913 to 1995. Coral Sr/Ca correlates with regional sea surface temperature (SST) back to 1966, when instrumental coverage is good, while coral δ18O does not. The slope of the monthly (annual mean) coral Sr/Ca-SST regression is −0.040 mmol/mol per 1 °C (−0.068 mmol/mol per 1 °C) consistent with published estimates of the Sr/Ca-SST relationship. Coral Sr/Ca suggest a warming of 0.39 °C since 1913. δ18O seawater is calculated by subtracting the temperature component from measured coral δ18O, using coral Sr/Ca as well as historical SST products. The derived δ18O seawater reconstructions are correlated (r 〉 0.6), and all show a significant shift in the midtwentieth century (−0.17‰ to −0.19‰), indicating a freshening of SSS by 0.7 psu. However, the timing of this shift depends on the temperature component and varies from 1947 (δ18O seawater calculated with historical SST) to the late 1950s (δ18O seawater calculated with coral Sr/Ca). Coral Sr/Ca shows warm temperature anomalies in the mid-1950s, while historical SST products show warm anomalies from 1940 to 1945 followed by cooling in the 1950s, a pattern typical for the World War II bias. This suggests that historical SST may bias reconstructions of δ18O seawater and SSS from corals.

Description: Evidence from proxy records indicates that millennial‐scale abrupt climate shifts, called Dansgaard‐Oeschger events, happened during past glacial cycles. Various studies have been conducted to uncover the physical mechanism behind them, based on the assumption that climate mean state determines the variability. However, our study shows that the Dansgaard‐Oeschger events can regulate the mean state of the Northern Hemisphere ice sheets. Sensitivity experiments show that the simulated mean state is influenced by the amplitude of the climatic noise. The most likely cause of this phenomenon is the nonlinear response of the surface mass balance to temperature. It could also cause the retreat processes to be faster than the buildup processes within a glacial cycle. We propose that the climate variability hindered ice sheet development and prevented the Earth system from entering a full glacial state from Marine Isotope Stage 4 to Marine Isotope Stage 3 about 60,000 years ago.

Description: The interaction between the atmosphere, specifically the North Atlantic Oscillation (NAO), and the North Atlantic ocean circulation on sub‐decadal timescale is analyzed in a subset of models participating in the Coupled Model Intercomparison Project phase 5 (CMIP5). From preindustrial control runs of at least 500 years length, we derive anomaly patterns in the atmospheric and ocean circulation and of air‐sea heat exchange. All models simulate a distinct dipolar oceanic overturning anomaly at the sub‐decadal timescale, with centers at 30° N and 55° N. The dipolar overturning anomaly goes along with marked anomalies in the North Atlantic sea surface temperature and gyre circulation. Lag‐regression analyses demonstrate, with relatively small ensemble spread, how the atmosphere and the ocean circulation interact. The dipolar anomalies in the overturning are forced by NAO‐related wind stress curl anomalies. Anomalous surface heat fluxes in concert with anomalous vertical motions drive a meridional dipolar heat content anomaly in the upper ocean, and it is this dipolar heat content anomaly which carries the coupled system from one phase of the sub‐decadal cycle to the other by reversing the tendencies in the overturning circulation. The coupled sub‐decadal variability derived from the CMIP5 models is characterized by three elements: a wind‐driven part steering the dipolar overturning anomaly, surface heat flux anomalies that support a heat build‐up in the subpolar gyre region, and the heat storage memory which is instrumental in the phase reversal of the NAO.

Description: Key Points: - Rapid subsurface oceanographic change in the tropical W Atlantic reflect shifting Subtropical Gyre - Subsurface warming responds to deglacial AMOC perturbations (Heinrich Stadials 2, 1, and the Younger Dryas) - Southward propagation of Salinity Maximum Water during Northern Hemisphere cold spells shift the mixing zone of tropical and subtropical waters During times of deglacial Atlantic Meridional Overturning Circulation (AMOC) perturbations, the tropical Atlantic experienced considerable warming at subsurface levels. Coupled ocean‐atmosphere simulations corroborate the tight teleconnection between the tropical Atlantic and climate change at high northern latitudes, but still underestimate the relevance of the subsurface N Atlantic Subtropical Gyre (STG) for heat and salt storage and its sensitivity to rapid climatic change. We here reconstruct vertical and lateral temperature and salinity gradients in the tropical W Atlantic and the Caribbean over the last 30 kyrs, based on planktic deep and shallow dwelling foraminiferal Mg/Ca and δ18O‐records. The rapid and large amplitude subsurface changes illustrate a dynamic STG associated with abrupt shifts of North Atlantic hydrographic and atmospheric regimes. During full glacial conditions, the STG has been shifted southward while intensified Ekman‐downwelling associated to strengthened trade winds fostered the formation of warm and saline Salinity Maximum Water (SMW). The southward propagation of SMW was facilitated by the glacially eastward deflected North Brazil Current. During periods of significant AMOC perturbations (Heinrich Stadials 1, and the Younger Dryas), extreme subsurface warming by ~6°C led to diminished lateral subsurface temperature gradients. Coevally, a deep thermocline suggests that SMW fully occupied the subsurface tropical W Atlantic and that the STG reached its southernmost position. During the Holocene, modern‐like conditions gradually developed with the northward retreat of SMW and the development of a strong thermocline ridge between the Subtropical Gyre and the tropical W Atlantic.

Description: Quaternary East Asian winter monsoon (EAWM) evolution has long been attributed to high-latitude Northern Hemisphere climate change. However, it cannot explain the distinct relationships of the EAWM in the northern and southern East Asian marginal sea in paleoclimatic records. Here we present an EAWM record of the northern East China Sea over the past 300 ka and a transient climate simulation with the Kiel Climate Model through the Holocene. Both proxy record and simulation suggest anticorrelated long-term EAWM evolution between the northern East China Sea and the South China Sea. We suggest that this spatial discrepancy of EAWM can be interpreted as El Niño–Southern Oscillation (ENSO)-like controlling, which generates cyclonic/anticyclonic wind anomalies in the northern/southern East Asian marginal sea. This research explains much of the controversy in nonorbital scale variability of Quaternary EAWM records in the East Asian marginal sea and supports a potent role of tropical forcing in East Asian winter climate change.

Description: New marine geophysical data acquired across the partly ice‐covered northern East Greenland continental margin highlight a complex interaction between tectonic and magmatic events. Breakup‐related lava flows are imaged in reflection seismic data as seaward dipping reflectors (SDRs), which are found to decrease in size both northwards and southwards from a central point at 75° N. We provide evidence that the magnetic anomaly pattern in the shelf area is related to volcanic phases and not to the presence of oceanic crust. The remnant magnetization of the individual lava flows is used to deduce a relative timing of the emplacement of the volcanic wedges. We find that the SDRs have been emplaced over a period of 2‐4 Ma progressively from north to south and from landward to seaward. The new data indicate a major post‐middle Eocene magmatic phase around the landward termination of the West Jan Mayen Fracture Zone. This post‐40 Ma volcanism likely was associated with the progressive separation of the Jan Mayen microcontinent from East Greenland. The break‐up of the Greenland Sea started at several isolated seafloor spreading cells whose location was controlled by rift structures and led to the present‐day segmentation of the margin. The original rift basins were subsequently connected by steady‐state seafloor spreading that propagated southwards, from the Greenland Fracture Zone to the Jan Mayen Fracture Zone. Key Points Polyphase Cenozoic volcanic rifting and consecutive emplacement of breakup‐related lava flows units along the northern East Greenland margin Breakup along restricted margin segments is followed by north to south directed progressive opening of the Greenland Sea Widespread post‐middle Eocene (〈 40 Ma) offshore magmatism, associated with the breakup of the Jan Mayen microcontinent from East Greenland

Description: The Eastern Tropical North Pacific Ocean hosts one of the world's largest oceanic oxygen deficient zones (ODZs). Hot spots for reactive nitrogen (N-r) removal processes, ODZs generate conditions proposed to promote N-r inputs via dinitrogen (N-2) fixation. In this study, we quantified N-2 fixation rates by N-15 tracer bioassay across oxygen, nutrient, and light gradients within and adjacent to the ODZ. Within subeuphotic oxygen-deplete waters, N-2 fixation was largely undetectable; however, addition of dissolved organic carbon stimulated N-2 fixation in suboxic (〈20 mu mol/kg O-2) waters, suggesting that diazotroph communities are likely energy limited or carbon limited and able to fix N-2 despite high ambient concentrations of dissolved inorganic nitrogen. Elevated rates (〉9 nmol N center dot L-1 center dot day(-1)) were also observed in suboxic waters near volcanic islands where N-2 fixation was quantifiable to 3,000 m. Within the overlying euphotic waters, N-2 fixation rates were highest near the continent, exceeding 500 mu mol N center dot m(-2)center dot day(-1) at one third of inshore stations. These findings support the expansion of the known range of diazotrophs to deep, cold, and dissolved inorganic nitrogen-replete waters. Additionally, this work bolsters calls for the reconsideration of ocean margins as important sources of N-r. Despite high rates at some inshore stations, regional N-2 fixation appears insufficient to compensate for N-r loss locally as observed previously in the Eastern Tropical South Pacific ODZ.

Description: During R/V Meteor cruise 141/1, pore fluids of near surface sediments were investigated to find indications for hydrothermal activity in the Terceira Rift (TR), a hyper‐slow spreading center in the Central North Atlantic Ocean. To date, submarine hydrothermal fluid venting in the TR has only been reported for the D. João de Castro seamount, which presently seems to be inactive. Pore fluids sampled close to a volcanic cone at 2800 m water depth show an anomalous composition with Mg, SO4, and total alkalinity (TA) concentrations significantly higher than seawater and a nearby reference core. The most straightforward way of interpreting these deviations is the dissolution of the hydrothermally formed mineral caminite (MgSO4 0.25Mg(OH)2 0.2H2O). This interpretation is corroborated by a thorough investigation of fluid isotope systems (δ26Mg, δ30Si, δ34S, δ44/42Ca, and 87Sr/86Sr). Caminite is known from mineral assemblages with anhydrite, and forms in hydrothermal recharge zones only under specific conditions such as high fluid temperatures and in altered oceanic crust, which are conditions generally met at the TR. We hypothesize that caminite was formed during hydrothermal activity and is now dissolving during the waning state of the hydrothermal system, so that caminite mineralization is shifted out of its stability zone. Ongoing fluid circulation through the basement is transporting the geochemical signal via slow advection towards the seafloor.

Description: Since 2009, unprecedented comprehensive long-term gas hydrate observations have become available from Ocean Networks Canada's NEPTUNE cabled ocean observatory at the northern Cascadia margin. Several experiments demonstrate the scientific importance of permanent power and Internet connectivity to the ocean floor as they have advanced the field of gas hydrate related research. One example is the cabled crawler Wally at Barkley Canyon, enabling live in situ exploration of the hydrate mounds and its associated benthic communities through the crawler's mobility and permanent accessibility throughout the year. Another example is a bubble-imaging sonar at Clayoquot Slope, revealing the strong relationship between ebullition of natural gas and tidal pressure, without apparent correlation to earthquakes, storms, or temperature fluctuations, in year-long continuous recordings. Finally, regular observatory maintenance cruises allow additional science sampling including echo-sounder surveys to extend the observatory footprint. Long-term trends in the data are not yet apparent but can also become evident from continuous measurements, as ocean observatories such as NEPTUNE are built for a 25-year lifetime, and expansion of the observatory networks makes these findings comparable and testable. Plain Language Summary Natural gas near the ocean floor creates a rapidly changing environment where it is important to collect data continuously in order to determine the magnitude, speed, and potential mechanism of change. This long-standing challenge of year-round access to the deep ocean has been tackled by Ocean Networks Canada through cabling the northern Cascadia seafloor, providing power and Internet communication-ideal for power-hungry instruments, large data volumes, and real-time access. The presence of gas influences the shape of the seafloor, animal activity, and potential escape of methane, a potent greenhouse gas. A seafloor crawler Wally was operated around deep canyon mounds of gas hydrate (a solid gas-water composite) since 2009 and helped discover environmental changes influencing sea life. Further along the continental slope, an acoustic sonar monitored rising methane bubbles where the bubbling appears to be controlled neither by earthquakes, winter storms, nor subtle temperature changes but actually strongly by tidal pressure. Regular maintenance of the observatory by ship allows more data to be collected near the cabled seafloor sites, extending the observations to a larger area. Ocean observatories are built to last decades and therefore more data for more research can be collected, potentially detecting relatively slow processes as well.

Description: There is a converging body of evidence supporting a measurable slowdown of the Atlantic Meridional Overturning Circulation (AMOC) as climate warms and Northern Hemisphere ice sheets inexorably shrink. Within this context, we assess the variability of the AMOC during the Holocene based on a marine sediment core retrieved from the deep northwest Atlantic, which sensitively recorded large‐scale deglacial transitions in deep water circulation. While there is a diffuse notion of Holocene variability in Labrador and Nordic Seas overturning, we report a largely invariable deep water circulation for the last ~11,000 years, even during the meltwater pulse associated with the 8.2‐ka event. Sensitivity tests along with high‐resolution 231Pa/230Th data constrain the duration and the magnitude of possible Holocene AMOC variations. The generally constant baseline during the Holocene suggests attenuated natural variability of the large‐scale AMOC on submillennial timescales and calls for compensating effects involving the upstream components of North Atlantic Deep Water.

Description: Numerical models have been highly successful in simulating global carbon and nutrient cycles in today's ocean, together with observed spatial and temporal patterns of chlorophyll and plankton biomass at the surface. With this success has come some confidence in projecting the century-scale response to continuing anthropogenic warming. There is also increasing interest in using such models to understand the role of plankton ecosystems in past oceans. However, today's marine environment is the product of billions of years of continual evolution—a process that continues today. In this paper, we address the questions of whether an assumption of species invariance is sufficient, and if not, under what circumstances current model projections might break down. To do this, we first identify the key timescales and questions asked of models. We then review how current marine ecosystem models work and what alternative approaches are available to account for evolution. We argue that for timescales of climate change overlapping with evolutionary timescales, accounting for evolution may to lead to very different projected outcomes regarding the timescales of ecosystem response and associated global biogeochemical cycling. This is particularly the case for past extinction events but may also be true in the future, depending on the eventual degree of anthropogenic disruption. The discipline of building new numerical models that incorporate evolution is also hugely beneficial in itself, as it forces us to question what we know about adaptive evolution, irrespective of its quantitative role in any specific event or environmental changes.

Description: Five initialization and ensemble generation methods are investigated with respect to their impact on the prediction skill of the German decadal prediction system "Mittelfristige Klimaprognose" (MiKlip). Among the tested methods, three tackle aspects of model‐consistent initialization using the ensemble Kalman filter (EnKF), the filtered anomaly initialization (FAI) and the initialization method by partially coupled spin‐up (MODINI). The remaining two methods alter the ensemble generation: the ensemble dispersion filter (EDF) corrects each ensemble member with the ensemble mean during model integration. And the bred vectors (BV) perturb the climate state using the fastest growing modes. The new methods are compared against the latest MiKlip system in the low‐resolution configuration (Preop‐LR), which uses lagging the climate state by a few days for ensemble generation and nudging toward ocean and atmosphere reanalyses for initialization. Results show that the tested methods provide an added value for the prediction skill as compared to Preop‐LR in that they improve prediction skill over the eastern and central Pacific and different regions in the North Atlantic Ocean. In this respect, the EnKF and FAI show the most distinct improvements over Preop‐LR for surface temperatures and upper ocean heat content, followed by the BV, the EDF and MODINI. However, no single method exists that is superior to the others with respect to all metrics considered. In particular, all methods affect the Atlantic Meridional Overturning Circulation in different ways, both with respect to the basin‐wide long‐term mean and variability, and with respect to the temporal evolution at the 26° N latitude.

Description: Ground-breaking measurements from the ocean observatories initiative Irminger Sea surface mooring (60°N, 39°30′W) are presented that provide the first in situ characterization of multiwinter surface heat exchange at a high latitude North Atlantic site. They reveal strong variability (December 2014 net heat loss nearly 50% greater than December 2015) due primarily to variations in frequency of intense short timescale (1–3 days) forcing. Combining the observations with the new high resolution European Centre for Medium Range Weather Forecasts Reanalysis 5 (ERA5) atmospheric reanalysis, the main source of multiwinter variability is shown to be changes in the frequency of Greenland tip jets (present on 15 days in December 2014 and 3 days in December 2015) that can result in hourly mean heat loss exceeding 800 W/m2. Furthermore, a new picture for atmospheric mode influence on Irminger Sea heat loss is developed whereby strongly positive North Atlantic Oscillation conditions favor increased losses only when not outweighed by the East Atlantic Pattern.

Description: The Bjerknes feedback is the dominant positive feedback in the equatorial ocean basins. To examine the seasonality, symmetry, and stationarity of the Pacific and Atlantic Bjerknes feedbacks we decompose them into three feedback elements that relate thermocline depth, sea surface temperature (SST), and western basin wind stress variability to each other. We partition feedback elements into composites associated with positive or negative anomalies. Using robust regression, we diagnose the strength of each composite. For the recent period 1993‐2012, composites of the Pacific Bjerknes feedback elements agree well with previous work. Positive composites are generally stronger than negative composites, and all feedback elements are weakest in late boreal spring. In the Atlantic, differences between positive and negative composites are less consistent across feedback elements. Specifically, wind variability seems to play a less important role in shaping atmosphere‐ocean coupling in the Atlantic when compared to the Pacific. However, a clear seasonality emerges: Feedback elements are generally strong in boreal summer and, for the negative composites, again in boreal winter. The Atlantic Bjerknes feedback is dominated by subsurface‐surface coupling. Applying our analysis to overlapping 25‐year periods for 1958‐2009 shows that the strengths of feedback elements in both ocean basins vary on decadal time scales. While the overall asymmetry of the Pacific Bjerknes feedback is robust, the strength and symmetry of Atlantic feedback elements vary considerably between decades. Our results indicate that the Atlantic Bjerknes feedback is non‐stationary on decadal time scales.

Description: Sea surface waves significantly affect the wind energy input to the Ekman layer in the upper ocean. In the study, we first incorporated the wave-induced Coriolis-Stokes forcing, the reduction of wind stress caused by wave generation, and wave dissipation into the classical Ekman model to investigate the kinetic energy balance in the wave-affected Ekman layer. Then, both the theoretical steady state solution for the idealized condition and the nonsteady state solution for the realistic ocean were derived. Total energy input to the wave-affected Ekman layer includes the wind stress energy input and the wave-induced energy input. Based on the WAVEWATCH III model, the wave spectrum was simulated to represent realistic random directional wave conditions. The wind stress energy input and the wave-induced energy input to the wave-affected Ekman layer in the Antarctic Circumpolar Current in the period from 1988 to 2010 were then calculated. The annual mean total energy input in the Antarctic Circumpolar Current region was 402.5 GW and the proportions of the wind stress energy input and the wave-induced energy input were, respectively, 85% and 15%. Particularly, total energy input in the Antarctic Circumpolar Current in the wave-affected Ekman layer model was 59.8% lower than that in the classical Ekman model. We conclude that surface waves play a significant role in the wind energy input to the Ekman layer.

Description: The export of organic carbon from the surface ocean forms the basis of the biological carbon pump, an important planetary carbon flux. Typically, only a small fraction of primary productivity (PP) is exported (quantified as the export efficiency: export/PP). Here we assemble a global data synthesis to reveal that very high export efficiency occasionally occurs. These events drive an apparent inverse relationship between PP and export efficiency, which is opposite to that typically used in empirical or mechanistic models. At the global scale, we find that low PP, high export efficiency regimes tend to occur when macrozooplankton and bacterial abundance are low. This implies that a decoupling between PP and upper ocean remineralization processes can result in a large fraction of PP being exported, likely as intact cells or phytoplankton-based aggregates. As the proportion of PP being exported declines, macrozooplankton and bacterial abundances rise. High export efficiency, high PP regimes also occur infrequently, possibly associated with nonbiologically mediated export of particles. A similar analysis at a biome scale reveals that the factors affecting export efficiency may be different at regional and global scales. Our results imply that the whole ecosystem structure, rather than just the phytoplankton community, is important in setting export efficiency. Further, the existence of low PP, high export efficiency regimes imply that biogeochemical models that parameterize export efficiency as increasing with PP may underestimate export flux during decoupled periods, such as at the start of the spring bloom.

Description: A pronounced warm anomaly occurred at the Peruvian coast in early 2017. This “Coastal Niño” caused heavy rainfalls, leading to flooding in Peru and Ecuador. At the same time, neutral conditions prevailed in the equatorial Pacific. Using observational sea surface temperature data sets and an ocean reanalysis product for the time period 1900 to 2010, previous similar events are investigated. Eighteen coastal warming events without corresponding equatorial Pacific warming are identified. Further analysis shows, however, that only four of these events are not connected to the central equatorial Pacific. All other periods of strong coastal warm anomalies are directly followed or preceded by El Niño‐like conditions. The “stand‐alone” coastal warming events are characterized by comparatively low equatorial heat content. We thus hypothesize that the depleted heat content in the equatorial Pacific in the wake of the strong 2015/2016 El Niño prevented the warming to spread westward in 2017.

Description: Southeast Greenland has been a major participant in the ice sheet mass loss over the last several decades. Interpreting the evolution of glacier fronts requires information about their depth below sea level and ocean thermal forcing, which are incompletely known in the region. Here, we combine airborne gravity and multibeam echo sounding data from the National Aeronautics and Space Administration's Oceans Melting Greenland (OMG) mission with ocean probe and fishing boat depth data to reconstruct the bathymetry extending from the glacier margins to the edges of the continental shelf. We perform a three‐dimensional inversion of the gravity data over water and merge the solution with a mass conservation reconstruction of bed topography over land. In contrast with other parts of Greenland, we find few deep troughs connecting the glaciers to the sources of warm Atlantic Water, amidst a relatively uniform, shallow (350 m) continental shelf. The deep channels include the Kangerlugssuaq, Sermilik, Gyldenløve, and Tingmiarmiut Troughs.

Description: International Ocean Discovery Program Expedition 352 recovered sedimentary‐volcaniclastic successions and extensional structures (faults and extensional veins) that allow the reconstruction of the Izu‐Bonin forearc tectonic evolution using a combination of shipboard core data, seismic reflection images, and calcite vein microstructure analysis. The oldest recorded biostratigraphic ages within fault‐bounded sedimentary basins (Late Eocene to Early Oligocene) imply a ~15 Ma hiatus between the formation of the igneous basement (52 to 50 Ma) and the onset of sedimentation. At the upslope sites (U1439 and U1442) extension led to the formation of asymmetric basins reflecting regional stretch of ~16–19% at strain rates of ~1.58 × 10−16 to 4.62 × 10−16 s−1. Downslope Site U1440 (closer to the trench) is characterized by a symmetric graben bounded by conjugate normal faults reflecting regional stretch of ~55% at strain rates of 4.40 × 10−16 to 1.43 × 10−15 s−1. Mean differential stresses are in the range of ~70–90 MPa. We infer that upper plate extension was triggered by incipient Pacific Plate rollback ~15 Ma after subduction initiation. Extension was accommodated by normal faulting with syntectonic sedimentation during Late Eocene to Early Oligocene times. Backarc extension was assisted by magmatism with related Shikoku and Parece‐Vela Basin spreading at ~25 Ma, so that parts of the arc and rear arc, and the West Philippine backarc Basin were dismembered from the forearc. This was followed by slow‐rift to postrift sedimentation during the transition from forearc to arc rifting to spreading within the Shikoku‐Parece‐Vela Basin system.

Description: Understanding oceanic cadmium (Cd) cycling is paramount due to its micronutrient‐like behavior in seawater, which has been inferred from its similarity to phosphate (PO4) cycling. Cadmium concentrations follow a nutrient‐like consumption‐regeneration cycle in the top of the water column and are mainly controlled by water mass mixing and circulation in the deep ocean. However, an additional scavenging mechanism through cadmium sulfide (CdS) precipitates, occurring within sinking biogenic particles in oxygen deficient zones (ODZ), has been proposed. In this study, we report Cd stable isotope and concentration data for seven vertical seawater profiles sampled during GEOTRACES cruise GA08 in the northern Cape and Angola Basins, which feature a significant ODZ along their eastern margins. Outside the ODZ, Cd cycling is similar to that previously reported for the South Atlantic. While water mass mixing largely controls deep ocean Cd isotope signatures, Cd isotope fractionation in surface waters can be modelled as an open system at steady‐state buffered by organic ligand complexation. In the ODZ, stronger Cd depletion relative to PO4 is associated with a shift in δ114Cd towards heavier values, which is indicative of CdS precipitation. Our interpretation is supported by experimental CdS precipitation data and a size‐resolved particle model involving bacterial sulfate reduction as a precursor of CdS. Our estimates of the CdS flux to the seafloor (107 to 109 mol yr‐1) indicate that CdS precipitation is a significant process of Cd removal and constitutes a non‐negligible Cd sink that needs to be better quantified by Cd isotope analyses of marine sediments.

Description: Ultraslow spreading ridges are poorly understood plate boundaries consisting of magmatic and amagmatic segments that expose mostly mantle peridotite and only traces of basalt and gabbro. The slowest part of the global spreading system is represented by the eastern Gakkel Ridge in the Central Arctic Ocean, where crustal accretion is characterized by extreme focusing of melt to discrete magmatic centers. Close to its eastern tip lies the unusual 5,310 m deep Gakkel Rift Deep (GRD) with limited sediment infill, which is in strong contrast to the broader sediment‐filled rift valleys to the east and west. Here, we report an 40Ar/39Ar age of 3.65±0.01 Ma for a pillow basalt from a seamount located on the rim the GRD confirming ultraslow spreading rates of ~7 mm/yr close to the Laptev Sea as suggested from aeromagnetic data. Its geochemistry points to an alkaline lava, attributed to partial melting of a source that underwent prior geochemical enrichment. We note that the GRD extracts compositionally similar melts as the sparsely magmatic zone further west but at much slower spreading velocities of only ~6‐7 mm/yr, indicating the widespread occurrence of similarly fertile mantle in the High Arctic. This enriched source differs from sub‐continental lithospheric mantle that influences magmatism along the Western Volcanic Zone (Goldstein et al. 2008) and is similar to metasomatized mantle ‐ shown to influence melt genesis along the Eastern Volcanic Zone.

Description: We present a 2‐D P wave velocity model and a coincident multichannel seismic reflection profile characterizing the structure of the southern Costa Rica margin and incoming Cocos Ridge. The seismic profiles image the ocean and overriding plates from the trench across the entire offshore margin, including the structures involved in the 2002 Osa earthquake. The overriding plate consists of three domains: Domain I displays thin‐skinned deformation of an imbricate thrust system composed of fractured rocks. Domain II shows ~15‐km‐long landward dipping reflection packages and active deformation of the shelf sediment. Domain III is little fractured and appears to be dominated by elastic deformation, overlain by ~2‐km‐thick landward dipping strata. The velocity structure supports the argument that the bulk of the margin is highly consolidated rock. Thick‐skinned tectonics probably causes the uplift of Domains II and III. The oceanic plate shows crustal thickness variations from ~14 km at the trench (Cocos Ridge) to 6–7 km beneath the shelf. We combine (1) interplate geometry and fracturing degree, (2) tectonic stresses and brittle strain, and (3) earthquake locations, to investigate relationships between structure and earthquake generation. The 2002 Osa sequence nucleated at the leading flank of subducting seamounts in the area of highest tectonic overpressure. Both estimated rock fracturing and modeled brittle strain steadily increase from the leading flank of the subducting seamounts to their top, reflecting the progressive damage caused by the seamount. Therefore, the seismicity and structural‐mechanical evolution of the upper plate reflect the downward propagation of the leading edge of seamounts.

Description: We examine the intra‐arc crustal seismicity of the Southern Andes Volcanic Zone (SVZ). Our aim is to resolve inter‐seismic deformation in an active magmatic arc dominated by both margin‐parallel (Liquiñe‐Ofqui fault system, (LOFS)) and Andean transverse faults (ATF). Crustal seismicity provides information about the schizosphere tectonic state, delineating the geometry and kinematics of high strain domains driven by oblique‐subduction. Here, we present local seismicity based on 16‐months data collected from 34 seismometers monitoring a ~200 km long section of the Southern Volcanic Zone, including the Lonquimay and Villarrica volcanoes. We located 356 crustal events with magnitudes between Mw 0.6 and Mw 3.6. Local seismicity occurs at depths down to 40 km in the forearc and consistently shallower than 12 km beneath the volcanic chain, suggesting a convex shape of the crustal seismogenic layer bottom. Focal mechanisms indicate strike‐slip faulting consistent with ENE‐WSW shortening in line with the long‐term deformation history revealed by structural geology studies. However, we find regional to local‐scale variations in the shortening axes orientation as revealed by the nature and spatial distribution of microseismicity, within three distinctive latitudinal domains. In the northernmost domain, seismicity is consistent with splay faulting at the northern termination of the LOFS; in the central domain, seismicity distributes along ENE‐ and WNW‐striking discrete faults, spatially associated with, hitherto seismic ATF. The southernmost domain, in turn, is characterized by activity focused along a N15°E striking master branch of the LOFS. These observations indicate a complex strain compartmentalization pattern within the intra‐arc crust, where variable strike‐slip faulting dominates over dip‐slip movements.

Description: The North Equatorial Undercurrent (NEUC) has been suggested to act as an important oxygen supply route towards the oxygen minimum zone in the Eastern Tropical North Atlantic. Observational estimates of the mean NEUC strength are uncertain due to the presence of elevated mesoscale activities, and models have difficulties in simulating a realistic NEUC. Here we investigate the interannual variability of the NEUC and its impact onto oxygen based on the output of a high‐resolution ocean general circulation model (OGCM) and contrast the results with an unique data set of 21 ship sections along 23° W and a conceptual model. We find that the interannual variability of the NEUC in the OGCM is related to the Atlantic Meridional Mode (AMM) with a stronger and more northward NEUC during negative AMM phases. Discrepancies between OGCM and observations suggest a different role of the NEUC in setting the regional oxygen distribution. In the model a stronger NEUC is associated with a weaker oxygen supply towards the east. We attribute this to a too strong recirculation between the NEUC and the northern branch of the South Equatorial Current (nSEC) in the OGCM. Idealized experiments with the conceptual model support the idea that the impact of NEUC variability on oxygen depends on the source water pathway. A strengthening of the NEUC supplied out of the western boundary acts to increase oxygen levels within the NEUC. A strengthening of the recirculations between NEUC and the nSEC results in a reduction of oxygen levels within the NEUC.

Description: The South-East Madagascar Bloom occurs in an oligotrophic region of the southwest Indian Ocean. Phase locked to austral summer, this sporadic feature exhibits substantial temporal and spatial variability. Several studies, with different hypotheses, have focused on the initiation mechanism triggering the bloom, but none has been as yet clearly substantiated. With 19 years of ocean color data set available as well as in situ measurements (Argo profiles), the time is ripe to review this feature. The bloom is characterized in a novel manner, and a new bloom index is suggested, yielding 11 bloom years, including 3 major bloom years (1999, 2006, and 2008). Spatially, the bloom varies from a mean structure (22–32°S; 50–70°E) both zonally and meridionally. A colocation analysis of Argo profiles and chlorophyll-a data revealed a bloom occurrence in a shallow-stratified layer, with low-salinity water in the surface layers. Additionally, a quantitative assessment of the previous hypotheses is performed and bloom occurrence is found to coincide with La Niña events and reduced upwelling intensity south of Madagascar. A stronger South-East Madagascar Current during La Niña may support a detachment of the current from the coasts, dampening the upwelling south of Madagascar, and feeding low-salinity waters into the Madagascar Basin, hence increasing stratification. Along with abundance of light, these provide the right conditions for a nitrogen-fixing cyanobacterial phytoplankton bloom onset

Description: Determination of Mg/Ca in foraminifer shells as a proxy of seawater temperature is of particular interest in paleoclimate reconstruction. Here we show that femtosecond–200 nm–laser ablation–inductively coupled plasma–mass spectrometry is a suitable technique to precisely and accurately determine Mg/Ca in the micrometer-sized calcareous chambers of foraminifers. At low fluence (0.3–0.6 J/cm 2 ) the double-charged 44 Ca ++ and the single-charged 25 Mg + ions are measured nearly simultaneously. Integrated single-shot measurements using a pulse repetition rate of 1 Hz enable precise analyses at a depth resolution of about 50–100 nm/pulse corresponding to an ablated material of 0.3–0.6 ng calcite/pulse for a spot size of 55 μm. High-resolution analyses can be performed until a depth of 10–20 μm and thus particularly suitable for thin-shelled foraminifers. Reproducibility (relative standard deviation) is about 5% as approved by homogeneous reference materials. Calibration is performed with the microanalytical synthetic reference material MACS-3. Magnesium and Ca data of different carbonate and silicate reference materials agree within uncertainties with reference values. The procedure has been successfully applied for detailed analyses of single chambers and shell-depth profiles of live individuals and empty planktic and benthic foraminifer tests from different ocean basins.

Description: Interannual variability of Antarctic Intermediate Water (AAIW) in the tropical North Atlantic is investigated using the GECCO2 ocean state estimate and Argo data. AAIW salinity variability near the western boundary is highly correlated with the transport along the western boundary on interannual timescales. Northward propagating anomalies are associated with the western boundary transport variability that, to some extent, is related to the large‐scale wind stress curl forcing by means of the Sverdrup balance. AAIW anomalies also propagate westward with the speed of baroclinic Rossby waves, indicating that the displacement of the meridional salinity gradient by westward propagation of baroclinic Rossby waves plays a role in the variability of AAIW characteristics. Slower eastward spreading of AAIW anomalies is identified on decadal timescales likely associated with the advection of salinity anomalies by weak eastward current bands. Understanding the observed interannual and decadal variability of AAIW salinity is important to properly interpret salinity changes reported in response to changes in the hydrological cycle.

Description: Springtime stratospheric final warming (SFW) variability has been suggested to be linked to the tropospheric circulation, particularly over the North Atlantic sector. These findings, however, are based on reanalysis data that cover a rather short period of time (1979 to present). The present work aims to improve the understanding of drivers, trends and surface impact of dynamical variability of boreal SFWs using chemistry‐climate models. We use multidecadal integrations of the fully coupled chemistry‐climate models Community Earth System Model version 1 (Whole Atmosphere Community Climate Model) and ECHAM/Modular Earth Submodel System Atmospheric Chemistry‐O. Four sensitivity experiments are analyzed to assess the impact of external factors; namely, the quasi‐biennial oscillation, sea surface temperature (SST) variability, and anthropogenic emissions. SFWs are classified into two types with respect to their vertical development; that is, events which occur first in the midstratosphere (10‐hPa first SFWs) or first in the upper stratosphere (1‐hPa first SFWs). Our results confirm previous reanalysis results regarding the differences in the time evolution of stratospheric conditions and near‐surface circulation between 10 and 1‐hPa first SFWs. Additionally, a tripolar SST pattern is, for the first time, identified over the North Atlantic in spring months related to the SFW variability. Our analysis of the influence of remote modulators on SFWs revealed that the occurrence of major warmings in the previous winter favors the occurrence of 10‐hPa first SFWs later on. We further found that quasi‐biennial oscillation and SST variability significantly affect the ratio between 1‐hPa first and 10‐hPa first SFWs. Finally, our results suggest that ozone recovery may impact the timing of the occurrence of 1‐hPa first SFWs.

Description: Mesoscale dynamics of the Agulhas Current system determine the exchange between the Indian and Atlantic oceans, thereby influencing the global overturning circulation. Using a series of ocean model experiments compared to observations, we show that the representation of mesoscale eddies in the Agulhas ring path improves with increasing resolution of submesoscale flows. Simulated submesoscale dynamics are validated with time‐mean horizontal‐wavenumber spectra from satellite sea surface temperature measurements and mesoscale dynamics with spectra from sea surface height. While the Agulhas ring path in a nonsubmesoscale‐resolving (1/20)° configuration is associated with too less power spectral densities on all scales and too steep spectral slopes, the representation of the mesoscale dynamics improves when the diffusion and the dissipation of the model are reduced and some small‐scale features are resolved. Realistic power spectral densities over all scales are achieved when additionally the horizontal resolution is increased to (1/60)° and a larger portion of the submesoscale spectrum is resolved. Results of an eddy detection algorithm applied to the model outputs as well as to a gridded sea surface height satellite product show that in particular strong cyclones are much better represented when submesoscale flows are resolved by the model. The validation of the submesoscale dynamics with sea surface temperature spectra provides guidance for the choice of advection schemes and explicit diffusion and dissipation as well as for further subgrid‐scale parameterizations. For the Agulhas ring path, the use of upstream biased advection schemes without explicit diffusion and dissipation is found to be associated with realistically simulated submesoscales.

Description: Large reservoirs of methane present in Arctic marine sediments are susceptible to rapid warming, promoting increasing methane emissions. Gas bubbles in the water column can be detected, and flow rates can be quantified using hydroacoustic survey methods, making it possible to monitor spatiotemporal variability. We present methane (CH4) bubble flow rates derived from hydroacoustic data sets acquired during 11 research expeditions to the western Svalbard continental margin (2008-2014). Three seepage areas emit in total 725-1,125 t CH4/year, and bubble fluxes are up to 2 kg.m(-2).year (-1). Bubble fluxes vary between different surveys, but no clear trend can be identified. Flux variability analyses suggest that two areas are geologically interconnected, displaying alternating flow changes. Spatial migration of bubble seepage was observed to follow seasonal changes in the theoretical landward limit of the hydrate stability zone, suggesting that formation/dissociation of shallow hydrates, modulated by bottom water temperatures, influences seafloor bubble release. Plain Language Summary It has been speculated that the release of methane (a potent greenhouse gas) from the seafloor in some Arctic Ocean regions is triggered by warming seawater. Emissions of gas bubbles from the seafloor can be detected by ship-mounted sonars. In 2008, a methane seepage area west of Svalbard was hydroacoustically detected for the first time. This seepage was hypothesized to be caused by dissociation of hydrates (ice-like crystals consisting of methane and water) due to ocean warming. We present an analysis of sonar data from 11 surveys conducted between 2008 and 2014. This study is the first comparison of methane seepage-related hydroacoustic data over such a long period. The hydroacoustic mapping and quantification method allowed us to assess the locations and intensity of gas bubble release, and how these parameters change over time, providing necessary data for numerical flux and climate models. No trend of increasing gas flow was identified. However, we observed seasonal variations potentially controlled by seasonal formation and dissociation of shallow hydrates. The hydrate formation/dissociation process is likely controlled by changes of bottom water temperatures. Alternating gas emissions between two neighboring areas indicate the existence of fluid pathway networks within the sediments.

Description: Free gas migration through the gas hydrate stability zone (GHSZ) and subsequent gas seepage at the seabed are characteristic features in marine gas hydrate provinces worldwide. The biogenic or thermogenic gas is typically transported along faults from deeper sediment strata to the GHSZ. Several mechanisms have been proposed to explain free gas transport through the GHSZ. While inhibition of hydrate formation by elevated salinities and temperatures have been addressed previously in studies simulating unfocused, area-wide upward advection of gas, which is not adequately supported by field observations, the role of focused gas flow through chimney-like structures has been underappreciated in this context. Our simulations suggest that gas migration through the GHSZ is, fundamentally, a result of methane gas supply in excess of its consumption by hydrate formation. The required high gas flux is driven by local overpressure, built up from gas accumulating below the base of the GHSZ that fractures the overburden when exceeding a critical pressure, thereby creating the chimney-like migration pathway. Initially rapid hydrate formation raises the temperature in the chimney structure, thereby facilitating further gas transport through the GHSZ. As a consequence, high hydrate saturations form preferentially close to the seafloor, where temperatures drop to bottom water values, producing a prominent subsurface salinity peak. Over time, hydrates form at a lower rate throughout the chimney structure, while initial temperature elevation and salinity peak dissipate. Thus, our simulations suggest that the near-surface salinity peak and elevated temperatures are a result of transient high-flux gas migration through the GHSZ.

Description: The presence of gas hydrates (GHs) increases the stiffness and strength of marine sediments. In elasto‐plastic constitutive models, it is common to consider GH saturation (Sh) as key internal variable for defining the contribution of GHs to composite soil mechanical behavior. However, the stress‐strain behavior of GH‐bearing sediments (GHBS) also depends on the microscale distribution of GH and on GH‐sediment fabrics. A thorough analysis of GHBS is difficult, because there is no unique relation between Sh and GH morphology. To improve the understanding of stress‐strain behavior of GHBS in terms of established soil models, this study summarizes results from triaxial compression tests with different Sh, pore fluids, effective confining stresses, and strain histories. Our data indicate that the mechanical behavior of GHBS strongly depends on Sh and GH morphology, and also on the strain‐induced alteration of GH‐sediment fabrics. Hardening‐softening characteristics of GHBS are strain rate‐dependent, which suggests that GH‐sediment fabrics dynamically rearrange during plastic yielding events. We hypothesize that rearrangement of GH‐sediment fabrics, through viscous deformation or transient dissociation and reformation of GHs, results in kinematic hardening, suppressed softening, and secondary strength recovery, which could potentially mitigate or counteract large‐strain failure events. For constitutive modeling approaches, we suggest that strain rate‐dependent micromechanical effects from alterations of the GH‐sediment fabrics can be lumped into a nonconstant residual friction parameter. We propose simple empirical evolution functions for the mechanical properties and calibrate the model parameters against the experimental data. Plain Language Summary Gas hydrates (GHs) are crystalline‐like solids, which are formed from natural gas molecules and water at high pressure and low temperature. GHs, and particularly methane hydrates, are naturally abundant in marine sediments. It is known that the presence of GH increases the mechanical stiffness and strength of sediments, and there is strong effort in analyzing and quantifying these effects in order to understand potential risks of sediment destabilization or slope failure. Based on our experimental results from high‐pressure geotechnical studies, we show that not only the initial amount and distribution of GH are important for the increased strength of GH‐bearing sediments but also the dynamic rearrangement of GH‐sediment fabrics during deformation characterizes the stress‐strain response and enables strength recovery after failure. We propose that different microstructural mechanisms contribute to this rearrangement and strength recovery of GH sediment. However, we consider these complicated processes in a simplified manner in an improved numerical model, which can be applied for geotechnical risk assessment on larger scales.

Description: The isotopic composition of neodymium dissolved in seawater consists of a distal, advected component that reflects water mass mixing and circulation, but near land can also contain a large local component originating from terrestrial sources such as aeolian or fluvial material. In order to use Nd isotopes to reconstruct paleocirculation, it is important to detect any local influences on the seawater signal recorded in deep sea sediments. Here we present rare earth element (REE) and Nd isotope (εNd) records from the deep Caribbean for two well‐studied time intervals in the Late Pliocene and Early Pleistocene. We measured trace element and REE compositions of weakly cleaned foraminifera to investigate if the Nd isotope signal from the same samples contained a local component. We find distinct changes in REE compositions across glaciations that are consistent with increases in the supply of local terrestrial material to the basin likely the results of glacially driven changes in sea level. Despite these larger terrestrial inputs, the Ce anomaly (Ce/Ce*) became more pronounced during glaciations indicating a better deep Caribbean ventilation. Short negative Nd isotope excursions occurred during three of the four studied glaciations, independently of any other proxy indicators for changes in ocean circulation suggesting that inputs from local terrigenous sources of Nd controlled the signal. We recommend that studies that aim to use εNd as a paleocirculation tracer routinely measure REE compositions of the authigenic phase to identify any possible terrestrial influence on the signal.

Description: A systematic study of Benguela Nino and Benguela Nina events during 1958 to 2015 including those that developed before the satellite era (1982) is carried out using an ocean general circulation model in combination with a linear equatorial model. Altogether, 21 strong warm and cold anomalous coastal events are identified among which 6 undocumented extreme coastal events are reported. Results suggest that most of these extreme coastal events including the newly identified ones are linked to remote equatorial forcing via mode 2 equatorial Kelvin waves. The latter propagates after approaching the African coast poleward as coastally trapped waves leading surface temperature anomalies along the Angola-Benguela current system by one month. One to two months before the peak of Benguela Ninos or Ninas usually occurring in March-April, a large-scale wind stress forcing is observed with both local (variations of alongshore coastal wind stress) and remote forcing developing simultaneously. Results further suggest that surface temperature anomalies off Southern Angola and in the Angola-Benguela Front are associated with equatorial dynamics and meridional wind stress fluctuations off the southwestern African coast north of 15 degrees S. Similar mechanisms are observed for Northern Namibia in combination with forcing by local meridional wind stress variations. Plain Language Summary The Benguela upwelling system located in the southeastern Atlantic Ocean supports a large marine ecosystem due to upwelling conditions. Every few years, anomalous warm and cold coastal events occur in the southeastern Atlantic and are detrimental for Angola, Namibia, and South Africa, as they affect fisheries and rainfall like El Nino phenomenon in the Pacific. To study these coastal events from 1958 to 2015, we use the output from a tropical Atlantic simulation in combination with the solution of a simple linear equatorial model. We study the anomalous coastal events including the ones that occurred before the satellite era (before 1982) and examine the role of the local wind forcing and the remote forcing associated with equatorial variability. We describe so far undocumented extreme events occurring from 1958 to 2015. Results suggest that most of the extreme coastal warm and cold events are associated with the propagation of equatorial Kelvin waves along the equatorial waveguide which trigger poleward-propagating coastal trapped waves along the southwestern African coast. One to two months before the peak season (usually March-April) of the anomalous coastal events, a large-scale wind pattern is observed, encompassing both variations of alongshore coastal wind in the southeastern Atlantic and zonal wind along the equatorial Atlantic.

Description: Astronomical tuning in the Mediterranean region is primarily based on organically‐mediated proxies, such as cyclicity of organic rich layers or changes in foraminiferal assemblages. Both during and post deposition, organic proxies can be affected by complex processes not immediately related to the changes in precession (insolation) they are assumed to reflect. Here we present an isotopic proxy which exhibits precessional cyclicity yet is inorganic. Seawater lead (Pb) isotope records over four precessional cycles between 6.6 and 6.5 Ma, from bulk sediment leachates of three Messinian, circum‐Mediterranean marginal locations, show variations consistent with precessional cyclicity. During insolation minima, the Pb isotope signatures from all three sites converge to similar values, suggesting a regional process is affecting all three locations at that time. Data from the marginal sites are compared with new data from ODP Site 978 and published data from a variety of geological archives from the Mediterranean region to determine the mechanism(s) causing the observed variability. While the comparisons are not fully conclusive, the timing of events suggest that increased dust production from North Africa during insolation minima is the most likely control. This hypothesis implies that authigenic marine Pb isotope records have the potential to provide a reliable inorganic tie point for Mediterranean cyclostratigraphy where sub‐precessional resolution is required. An inorganic tie point could also provide the means to resolve long‐standing problems in Mediterranean stratigraphy on precessional and sub‐precessional timescales which have been obscured due to post‐depositional changes (e.g., sapropel burn‐down) or suboptimal ecological conditions (e.g., the Messinian Salinity Crisis).

Description: We investigate the origin of fresh water on the shelves near Cape Farewell (south Greenland) using sections of three hydrographic cruises in May (HUD2014007) and June 2014 (JR302 and Geovide). We partition the fresh water between meteoric water sources and sea ice melt or brine formation using the δ18O of sea water. The sections illustrate the presence of the East Greenland Coastal Current (EGCC) close to shore east of Cape Farewell. West of Cape Farewell, it partially joins the shelf break, with a weaker near‐surface remnant of the EGCC observed on the shelf southwest and west of Cape Farewell. The EGCC traps the freshest waters close to Greenland and carries a brine signature below 50‐m depth. The cruises illustrate a strong increase in meteoric water of the shelf upper layer (by more than a factor 2) between early May and late June, likely to result from East and South Greenland spring melt. There was also a contribution of sea ice melt near the surface but with large variability both spatially and also between the two June cruises. Furthermore, gradients in the freshwater distribution and its contributions are larger east of Cape Farewell than west of Cape Farewell, which is related to the EGCC being more intense and closer to the coast east of Cape Farewell than west of it. Large temporal variability in the currents is found between different sections to the east and southeast of Cape Farewell, likely related to changes in wind conditions.

Description: A new global climate model setup using FESOM2.0 for the sea ice‐ocean component and ECHAM6.3 for the atmosphere and land surface has been developed. Replacing FESOM1.4 by FESOM2.0 promises a higher efficiency of the new climate setup compared to its predecessor. The new setup allows for long‐term climate integrations using a locally eddy‐resolving ocean. Here it is evaluated in terms of (1) the mean state and long‐term drift under preindustrial climate conditions, (2) the fidelity in simulating the historical warming, and (3) differences between coarse and eddy‐resolving ocean configurations. The results show that the realism of the new climate setup is overall within the range of existing models. In terms of oceanic temperatures, the historical warming signal is of smaller amplitude than the model drift in case of a relatively short spin‐up. However, it is argued that the strategy of “de‐drifting” climate runs after the short spin‐up, proposed by the HighResMIP protocol, allows one to isolate the warming signal. Moreover, the eddy‐permitting/resolving ocean setup shows notable improvements regarding the simulation of oceanic surface temperatures, in particular in the Southern Ocean.

Description: Sections PDFPDF Tools Share Abstract Long‐term data characterizing the oceans' biological carbon pump are essential for understanding impacts of climate variability on marine ecosystems. The “Bakun upwelling intensification hypothesis” suggests intensified coastal upwelling due to a greater land‐sea temperature gradient influenced by global warming. We present long time series of bathypelagic (approximately 1,200–3,600 m) particle fluxes from a coastal (CBeu: 2003–2016) and an offshore (CBmeso: 1988–2016) sediment trap setting located in the Canary Current upwelling. Organic carbon (Corg) and biogenic opal (BSi, diatoms) fluxes were twofold to threefold higher at the coastal upwelling site compared to the offshore site, respectively, and showed higher seasonality with flux maxima in spring. A relationship between winter and spring BSi fluxes to the North Atlantic Oscillation index was best expressed at the offshore site CBmeso. Lithogenic (dust) fluxes regularly peaked in winter when frequent low‐altitude dust storms and deposition occurred, decreasing offshore by about threefold. We obtained a high temporal match of short‐term peaks of BSi and dust fluxes in winter to spring at the inner site CBeu. We found synchronous flux variations at both sites and an anomalous year 2005, characterized by high BSi and Corg fluxes under a low North Atlantic Oscillation. Corg and BSi fluxes revealed a decreasing trend from 2006 to 2016 at the coastal site CBeu, pointing to coastal upwelling relaxation during the last two decades. The permanent offshore upwelling zone of the deflected Canary Current represented by the flux record of CBmeso showed no signs of increasing upwelling as well which contradicts the Bakun hypothesis.

Description: Magma emplacement in organic‐rich sedimentary basins is a main driver of past environmental crises. Using a 2D numerical model, we investigate the process of thermal cracking in contact aureoles of cooling sills and subsequent transport and emission of thermogenic methane by hydrothermal fluids. Our model includes a Mohr‐Coulomb failure criterion to initiate hydrofracturing and a dynamic porosity/permeability. We investigate the Karoo Basin, taking into account host‐rock material properties from borehole data, realistic total organic carbon content, and different sill geometries. Consistent with geological observations, we find that thermal plumes quickly rise at the edges of saucer‐shaped sills, guided along vertically fractured high permeability pathways. Contrastingly, less focused and slower plumes rise from the edges and the central part of flat‐lying sills. Using a novel upscaling method based on sill‐to‐sediment ratio we find that degassing of the Karoo Basin occurred in two distinct phases during magma invasion. Rapid degassing triggered by sills emplaced within the top 1.5 km emitted ~1.6·103 Gt of thermogenic methane, while thermal plumes originating from deeper sills, carrying a 12‐times greater mass of methane, may not reach the surface. We suggest that these large quantities of methane could be re‐mobilized by the heat provided by neighboring sills. We conclude that the Karoo LIP may have emitted as much as ~22.3·103 Gt of thermogenic methane in the half million years of magmatic activity, with emissions up to 3 Gt/year. This quantity of methane and the emission rates can explain the negative δ13C excursion of the Toarcian environmental crisis. Key Points Sill geometry and emplacement depth as well as intruded host rock type are the main factors controlling methane mobilization and degassing Dehydration‐related porosity increase and pore‐pressure‐induced hydrofracturing are important mechanisms for a quick transport of methane from sill to the surface The Karoo Basin may have degassed ~22.3·103 Gt of thermogenic methane in the half million years of magmatic activity

Description: The Atlantic Subtropical Cells (STCs) are shallow wind‐driven overturning circulations connecting the tropical upwelling areas to the subtropical subduction regions. In both hemispheres they are characterized by equatorward transport at thermocline level, upwelling at the equator and poleward Ekman transport in the surface layer. This study uses recent data from Argo oats complemented by ship sections at the western boundary as well as reanalysis products to estimate the meridional water mass transports and to investigate the vertical and horizontal structure of the STCs from an observational perspective. The seasonally varying depth of meridional velocity reversal is used as the interface between the surface poleward ow and the thermocline equatorward ow. The latter is bounded by the 26.0 kg m‐3 isopycnal at depth. We find that the thermocline layer convergence is dominated by the southern hemisphere water mass transport (9.0 ±1.1 Sv from the southern hemisphere compared to 2.9 ±1.3 Sv from the northern hemisphere) and that this transport is mostly confined to the western boundary. Compared to the asymmetric convergence at thermocline level, the wind‐driven Ekman divergence in the surface layer is more symmetric, being 20.4 ±3.1 Sv between 10°N and 10°S. The net poleward transports (Ekman minus geostrophy) in the surface layer concur with values derived from reanalysis data (5.5 ±0.8 Sv at 10°S and 6.4 ±1.4 Sv at 10°N). A diapycnal transport of about 4 Sv across the 26.0 kg m‐3 isopycnal is required in order to maintain the mass balance in the STC circulation.

Description: Precise knowledge of wintertime sea ice production in Arctic polynyas is not only required to enhance our understanding of atmosphere‐sea ice‐ocean interactions but also to verify frequently utilized climate and ocean models. Here, a high‐resolution (2‐km) Moderate Resolution Imaging Spectroradiometer (MODIS) thermal infrared satellite data set featuring spatial and temporal characteristics of 17 Arctic polynya regions for the winter seasons 2002/2003 to 2017/2018 is directly compared to an akin low‐resolution Advanced Microwave Scanning Radiometer‐EOS (AMSR‐E) passive microwave data set for 2002/2003 to 2010/2011. The MODIS data set is purely based on a 1‐D energy‐balance model, where thin‐ice thicknesses (≤ 20 cm) are directly derived from ice‐surface temperature swath data and European Centre for Medium‐Range Weather Forecasts Re‐Analysis‐Interim atmospheric reanalysis data on a quasi‐daily basis. Thin‐ice thicknesses in the AMSR‐E data set are derived empirically. Important polynya properties such as areal extent and potential thermodynamic ice production can be estimated from both pan‐Arctic data sets. Although independently derived, our results show that both data sets feature quite similar spatial and temporal variations of polynya area (POLA) and ice production (IP), which suggests a high reliability. The average POLA (average accumulated IP) for all Arctic polynyas combined derived from both MODIS and AMSR‐E are 1.99×105 km2 (1.34×103 km3) and 2.29×105 km2 (1.31×103 km3), respectively. Narrow polynyas in areas such as the Canadian Arctic Archipelago are notably better resolved by MODIS. Analysis of 16 winter seasons provides an evaluation of long‐term trends in POLA and IP, revealing the significant increase of ice formation in polynyas along the Siberian coast.

Description: A new release of the Max Planck Institute for Meteorology Earth System Model version 1.2 (MPI-ESM1.2) is presented. The development focused on correcting errors in and improving the physical processes representation, as well as improving the computational performance, versatility, and overall user friendliness. In addition to new radiation and aerosol parameterizations of the atmosphere, several relatively large, but partly compensating, coding errors in the model's cloud, convection, and turbulence parameterizations were corrected. The representation of land processes was refined by introducing a multilayer soil hydrology scheme, extending the land biogeochemistry to include the nitrogen cycle, replacing the soil and litter decomposition model and improving the representation of wildfires. The ocean biogeochemistry now represents cyanobacteria prognostically in order to capture the response of nitrogen fixation to changing climate conditions and further includes improved detritus settling and numerous other refinements. As something new, in addition to limiting drift and minimizing certain biases, the instrumental record warming was explicitly taken into account during the tuning process. To this end, a very high climate sensitivity of around 7 K caused by low-level clouds in the tropics as found in an intermediate model version was addressed, as it was not deemed possible to match observed warming otherwise. As a result, the model has a climate sensitivity to a doubling of CO2 over preindustrial conditions of 2.77 K, maintaining the previously identified highly nonlinear global mean response to increasing CO2 forcing, which nonetheless can be represented by a simple two-layer model.

Description: The Weddell Gyre (WG) is one of the main oceanographic features of the Southern Ocean south of the Antarctic Circumpolar Current which plays an influential role in global ocean circulation as well as gas exchange with the atmosphere. We review the state‐of‐the art knowledge concerning the WG from an interdisciplinary perspective, uncovering critical aspects needed to understand this system's role in shaping the future evolution of oceanic heat and carbon uptake over the next decades. The main limitations in our knowledge are related to the conditions in this extreme and remote environment, where the polar night, very low air temperatures, and presence of sea ice year‐round hamper field and remotely sensed measurements. We highlight the importance of winter and under‐ice conditions in the southern WG, the role that new technology will play to overcome present‐day sampling limitations, the importance of the WG connectivity to the low‐latitude oceans and atmosphere, and the expected intensification of the WG circulation as the westerly winds intensify. Greater international cooperation is needed to define key sampling locations that can be visited by any research vessel in the region. Existing transects sampled since the 1980s along the Prime Meridian and along an East‐West section at ~62°S should be maintained with regularity to provide answers to the relevant questions. This approach will provide long‐term data to determine trends and will improve representation of processes for regional, Antarctic‐wide, and global modeling efforts—thereby enhancing predictions of the WG in global ocean circulation and climate.

Description: Global biogeochemical ocean models rely on many parameters, which govern the interaction between individual components, and their response to the physical environment. They are often assessed/calibrated against quasi-synoptic data sets of dissolved inorganic tracers. However, a good fit to one observation might not necessarily imply a good match to another. We investigate whether two different metrics—the root-mean-square error to nutrients and oxygen and a metric measuring the overlap between simulated and observed oxygen minimum zones (OMZs)—help to constrain a global biogeochemical model in different aspects of performance. Three global model optimizations are carried out. Two single-objective optimizations target the root-mean-square metric and a sum of both metrics, respectively. We then present and explore multiobjective optimization, which results in a set of compromise solutions. Our results suggest that optimal parameters for denitrification and nitrogen fixation differ when applying different metrics. Optimization against observed OMZs leads to parameters that enhance fixed nitrogen cycling; this causes too low nitrate concentrations and a too high global pelagic denitrification rate. Optimization against nutrient and oxygen concentrations leads to different parameters and a lower global fixed nitrogen turnover; this results in a worse fit to OMZs. Multiobjective optimization resolves this antagonistic effect and provides an ensemble of parameter sets, which help to address different research questions. We finally discuss how systematic model calibration can help to improve models used for projecting climate change and its effect on fisheries and climate gas emissions.

Description: Marine sediments host large amounts of methane (CH4), which is a potent greenhouse gas. Quantitative estimates for methane release from marine sediments are scarce, and a poorly constrained temporal variability leads to large uncertainties in methane emission scenarios. Here, we use 2D and 3D seismic reflection, multibeam bathymetric, geochemical and sedimentological data to (I) map and describe pockmarks in the Witch Ground Basin (central North Sea), (II) characterize associated sedimentological and fluid migration structures, and (III) analyze the related methane release. More than 1500 pockmarks of two distinct morphological classes spread over an area of 225 km2. The two classes form independently from another and are corresponding to at least two different sources of fluids. Class 1 pockmarks are large in size (〉 6 m deep, 〉 250 m long, and 〉 75 m wide), show active venting, and are located above vertical fluid conduits that hydraulically connect the seafloor with deep methane sources. Class 2 pockmarks, which comprise 99.5 % of all pockmarks, are smaller (0.9‐3.1 m deep, 26‐140 m long, and 14‐57 m wide) and are limited to the soft, fine‐grained sediments of the Witch Ground Formation and possibly sourced by compaction‐related dewatering. Buried pockmarks within the Witch Ground Formation document distinct phases of pockmark formation, likely triggered by external forces related to environmental changes after deglaciation. Thus, greenhouse gas emissions from pockmark fields cannot be based on pockmark numbers and present‐day fluxes but require an analysis of the pockmark forming processes through geological time.

Description: The Lagrangian method-where current location and intensity are determined by tracking the movement of flow along its path-is the oldest technique for measuring the ocean circulation. For centuries, mariners used compilations of ship drift data to map out the location and intensity of surface currents along major shipping routes of the global ocean. In the mid-20th century, technological advances in electronic navigation allowed oceanographers to continuously track freely drifting surface buoys throughout the ice-free oceans and begin to construct basin-scale, and eventually global-scale, maps of the surface circulation. At about the same time, development of acoustic methods to track neutrally buoyant floats below the surface led to important new discoveries regarding the deep circulation. Since then, Lagrangian observing and modeling techniques have been used to explore the structure of the general circulation and its variability throughout the global ocean, but especially in the Atlantic Ocean. In this review, Lagrangian studies that focus on pathways of the upper and lower limbs of the Atlantic Meridional Overturning Circulation (AMOC), both observational and numerical, have been gathered together to illustrate aspects of the AMOC that are uniquely captured by this technique. These include the importance of horizontal recirculation gyres and interior (as opposed to boundary) pathways, the connectivity (or lack thereof) of the AMOC across latitudes, and the role of mesoscale eddies in some regions as the primary AMOC transport mechanism. There remain vast areas of the deep ocean where there are no direct observations of the pathways of the AMOC.

Description: To maintain the chance of keeping the average global temperature increase below 2 degrees C and to limit long-term climate change, removing carbon dioxide from the atmosphere (carbon dioxide removal, CDR) is becoming increasingly necessary. We analyze optimal and cost-effective climate policies in the dynamic integrated assessment model (IAM) of climate and the economy (DICE2016R) and investigate (1) the utilization of (ocean) CDR under different climate objectives, (2) the sensitivity of policies with respect to carbon cycle feedbacks, and (3) how well carbon cycle feedbacks are captured in the carbon cycle models used in state-of-the-art IAMs. Overall, the carbon cycle model in DICE2016R shows clear improvements compared to its predecessor, DICE2013R, capturing much better long-term dynamics and also oceanic carbon outgassing due to excess oceanic storage of carbon from CDR. However, this comes at the cost of a (too) tight short-term remaining emission budget, limiting the model suitability to analyze low-emission scenarios accurately. With DICE2016R, the compliance with the 2 degrees C goal is no longer feasible without negative emissions via CDR. Overall, the optimal amount of CDR has to take into account (1) the emission substitution effect and (2) compensation for carbon cycle feedbacks.

Description: During the summer monsoon, the western tropical Indian Ocean is predicted to be a hot spot for dimethylsulfide emissions, the major marine sulfur source to the atmosphere, and an important aerosol precursor. Other aerosol relevant fluxes, such as isoprene and sea spray, should also be enhanced, due to the steady strong winds during the monsoon. Marine air masses dominate the area during the summer monsoon, excluding the influence of continentally derived pollutants. During the SO234-2/235 cruise in the western tropical Indian Ocean from July to August 2014, directly measured eddy covariance DMS fluxes confirm that the area is a large source of sulfur to the atmosphere (cruise average 9.1 μmol m−2 d−1). The directly measured fluxes, as well as computed isoprene and sea spray fluxes, were combined with FLEXPART backward and forward trajectories to track the emissions in space and time. The fluxes show a significant positive correlation with aerosol data from the Terra and Suomi-NPP satellites, indicating a local influence of marine emissions on atmospheric aerosol numbers.

Description: Anthropogenic activities have resulted in enhanced lead (Pb) emissions to the environment over the past century, mainly through the combustion of leaded gasoline. Here, we present the first combined dissolved (DPb), labile (LpPb) and particulate (PPb) Pb dataset from the Northeast Atlantic (Celtic Sea) since the phasing out of leaded gasoline in Europe. Concentrations of DPb in surface waters have decreased by 4-fold over the last four decades. We demonstrate that anthropogenic Pb is transported from the Mediterranean Sea over long distances (〉2500 km). Benthic DPb fluxes exceeded the atmospheric Pb flux in the region, indicating the importance of sediments as a contemporary Pb source. A strong positive correlation between DPb, PPb and LpPb indicates a dynamic equilibrium between the phases and the potential for particles to ‘buffer’ the DPb pool. This study provides insights into Pb biogeochemical cycling and demonstrates the potential of Pb in constraining ocean circulation patterns.

Description: The variability of the Atlantic Meridional Overturning Circulation (AMOC) may play a role in sea surface temperature predictions on seasonal to decadal time scales. Therefore, AMOC seasonal cycles are a potential baseline for interpreting predictions. Here we present estimates for the seasonal cycle of transports of volume, temperature, and freshwater associated with the upper limb of the AMOC in the eastern subpolar North Atlantic on the Extended Ellett Line hydrographic section between Scotland and Iceland. Due to weather, ship‐based observations are primarily in summer. Recent glider observations during other seasons present an opportunity to investigate the seasonal variability in the upper layer of the AMOC. First, we document a new method to quality control and merge ship, float, and glider hydrographic observations. This method accounts for the different spatial sampling rates of the three platforms. The merged observations are used to compute seasonal cycles of volume, temperature, and freshwater transports in the Rockall Trough. These estimates are similar to the seasonal cycles in two eddy‐resolving ocean models. Volume transport appears to be the primary factor modulating other Rockall Trough transports. Finally, we show that the weakest transports occur in summer, consistent with seasonal changes in the regional‐scale wind stress curl. Although the seasonal cycle is weak compared to other variability in this region, the amplitude of the seasonal cycle in the Rockall Trough, roughly 0.5–1 Sv about a mean of 3.4 Sv, may account for up to 7–14% of the heat flux between Scotland and Greenland.

Description: We reanalyze existing paleodata of global mean surface temperature ΔTg and radiative forcing ΔR of CO2 and land ice albedo for the last 800,000 years to show that a state‐dependency in paleoclimate sensitivity S, as previously suggested, is only found if ΔTg is based on reconstructions, and not when ΔTg is based on model simulations. Furthermore, during times of decreasing obliquity (periods of land‐ice sheet growth and sea level fall) the multi‐millennial component of reconstructed ΔTg diverges from CO2, while in simulations both variables vary more synchronously, suggesting that the differences during these times are due to relatively low rates of simulated land ice growth and associated cooling. To produce a reconstruction‐based extrapolation of S for the future we exclude intervals with strong ΔTg‐CO2 divergence and find that S is less state‐dependent, or even constant (state‐independent), yielding a mean equilibrium warming of 2–4 K for a doubling of CO2.

Description: Provenance studies of widely distributed tephras, integrated within a well-defined temporal framework, are important to deduce systematic changes in the source, scale, distribution and changes in regional explosive volcanism. Here, we establish a robust tephro-chronostratigraphy for a total of 157 marine tephra layers collected during IODP Expedition 352. We infer at least three major phases of highly explosive volcanism during Oligocene to Pleistocene time. Provenance analysis based on glass composition assigns 56 of the tephras to a Japan source, including correlations with 12 major and widespread tephra layers resulting from individual eruptions in Kyushu, Central Japan and North Japan between 115 ka and 3.5 Ma. The remaining 101 tephras are assigned to four source regions along the Izu-Bonin arc. One, of exclusively Oligocene age, is proximal to the Bonin Ridge islands; two reflect eruptions within the volcanic front and back-arc of the central Izu-Bonin arc, and a fourth region corresponds to the Northern Izu-Bonin arc source. First-order volume estimates imply eruptive magnitudes ranging from 6.3 to 7.6 for Japan-related eruptions and between 5.5 and 6.5 for IBM eruptions. Our results suggest tephras between 30 and 22 Ma that show a subtly different Izu-Bonin chemical signature compared to the recent arc. After a ∼11 m.y. gap in eruption, tephra supply from the Izu-Bonin arc predominates from 15 to 5 Ma, and finally a subequal mixture of tephra sources from the (palaeo)Honshu and Izu-Bonin arcs occurs within the last ∼5 Ma.

Description: Subsurface coherent eddies are well-known features of ocean circulation, but the sparsity of observations prevents an assessment of their importance for biogeochemistry. Here, we use a global eddying (0.1° ) ocean-biogeochemical model to carry out a census of subsurface coherent eddies originating from eastern boundary upwelling systems (EBUS), and quantify their biogeochemical effects as they propagate westward into the subtropical gyres. While most eddies exist for a few months, moving over distances of 100s of km, a small fraction (〈 5%) of long-lived eddies propagates over distances greater than 1000km, carrying the oxygen-poor and nutrient-rich signature of EBUS into the gyre interiors. In the Pacific, transport by subsurface coherent eddies accounts for roughly 10% of the offshore transport of oxygen and nutrients in pycnocline waters. This "leakage" of subsurface waters can be a significant fraction of the transport by nutrient-rich poleward undercurrents, and may contribute to the well-known reduction of productivity by eddies in EBUS. Furthermore, at the density layer of their cores, eddies decrease climatological oxygen locally by close to 10%, thereby expanding oxygen minimum zones. Finally, eddies represent low-oxygen extreme events in otherwise oxygenated waters, increasing the area of hypoxic waters by several percent and producing dramatic short-term changes that may play an important ecological role. Capturing these non-local effects in global climate models, which typically include non-eddying oceans, would require dedicated parameterizations.

Description: Climate models depict large diversity in the strength of the El Niño/Southern Oscillation (ENSO) (ENSO amplitude). Here we investigate ENSO-amplitude diversity in the Coupled Model Intercomparison Project Phase 5 (CMIP5) by means of the linear recharge oscillator model, which reduces ENSO dynamics to a two-dimensional problem in terms of eastern equatorial Pacific sea surface temperature anomalies (T) and equatorial Pacific upper ocean heat content anomalies (h). We find that a large contribution to ENSO-amplitude diversity originates from stochastic forcing. Further, significant interactions exist between the stochastic forcing and the growth rates of T and h with competing effects on ENSO amplitude. The joint consideration of stochastic forcing and growth rates explains more than 80% of the ENSO-amplitude variance within CMIP5. Our results can readily explain the lack of correlation between the Bjerknes Stability index, a measure of the growth rate of T, and ENSO amplitude in a multimodel ensemble.

Description: Monthly mean sea level anomalies in the tropical Pacific for the period 1961-2002 are reconstructed using a linear, multi-mode model driven by monthly mean wind stress anomalies from the NCEP/NCAR and ERA-40 reanalysis products. Overall, the sea level anomalies reconstructed by both wind stress products agree well with the available tide gauge data, although with poor performance at Kanton Island in the western-central equatorial Pacific and reduced amplitude at Christmas Island. The reduced performance is related to model error in locating the pivot point in sea level variability associated with the so-called “tilt” mode. We present evidence that the pivot point was further west during the period 1993-2014 than during the period 1961-2002 and attribute this to a persistent upward trend in the zonal wind stress variance along the equator west of 160° W throughout the period 1961-2014. Experiments driven by the zonal component of the wind stress alone reproduce much of the trend in sea level found in the experiments driven by both components of the wind stress. The experiments show an upward trend in sea level in the eastern tropical Pacific over the period 1961-2002, but with a much stronger upward trend when using the NCEP/NCAR product. We argue that the latter is related to an overly strong eastward trend in zonal wind stress in the eastern-central Pacific that is believed to be a spurious feature of the NCEP/NCAR product.

Description: Continental margin ecosystems in the western North Pacific Ocean are subject to strong climate forcing and anthropogenic impacts. To evaluate mechanisms controlling phytoplankton biomass and community structure variations in marginal sea‐open ocean boundary regions, brassicasterol, dinosterol and C37 alkenones were measured in suspended particles in summer and autumn from 2012 to 2013 in the northeastern East China Sea and the western Tsushima Strait (NEECS‐WTS). In summer, the concentrations of brassicasterol (40 ‐ 1535 ng L‐1) and dinosterol (4.2 ‐ 94 ng L‐1) were higher in the southwest of Cheju Island, while C37 alkenones (0 ‐ 30 ng L‐1) were higher in the south of Cheju Island. In autumn, brassicasterol (12 ‐ 106 ng L‐1), dinosterol (2.4 ‐ 21 ng L‐1) and C37 alkenones (0.7 – 7.0 ng L‐1) were higher in the southwest of Cheju Island and the WTS, and higher C37 alkenones also occurred in the Okinawa Trough. Correlation analysis of biomarkers and environmental conditions (temperature, salinity and inorganic nutrient concentrations) clearly demonstrated that phytoplankton biomass and community structure variations can be well elucidated by water masses as indexed by temperature and salinity. High nutrients from the Changjiang River were the main cause of high biomass in summer, while nutrients from subsurface water were likely the key factor regulating phytoplankton biomass in open ocean water stations in autumn. This study indicates that mechanisms controlling phytoplankton biomass in marginal sea‐open ocean boundary regions should be classified by various water masses with different nutrient concentrations, instead of by geography.

Description: A submesoscale coherent vortex (SCV) with a low oxygen core is characterized from underwater glider and mooring observations from the eastern tropical North Atlantic, north of the Cape Verde Islands. The eddy crossed the mooring with its center and a 1 month time series of the SCV's hydrographic and upper 100 m currents structure was obtained. About 45 days after, and ∼100 km west, the SCV frontal zone was surveyed in high temporal and spatial resolution using an underwater glider. Satellite altimetry showed the SCV was formed about 7 months before at the Mauritanian coast. The SCV was located at 80-100 m depth, its diameter was ∼100 km and its maximum swirl velocity ∼0.4 m s-1. A Burger number of 0.2 and a vortex Rossby number 0.15 indicate a flat lens in geostrophic balance. Mooring and glider data show in general comparable dynamical and thermohaline structures, the glider in high spatial resolution, the mooring in high temporal resolution. Surface maps of chlorophyll concentration suggest high productivity inside and around the SCV. The low potential vorticity (PV) core of the SCV is surrounded by filamentary structures, sloping down at different angles from the mixed layer base and with typical width of 10-20 km and a vertical extent of 50-100 m.

Description: Earth is 180 Myr into the current supercontinent cycle, and the next supercontinent is predicted to form in 250 Myr. The continuous changes in continental configuration can move the ocean between resonant states, and the semidiurnal tides are currently large compared to the past 252 Myr due to tidal resonance in the Atlantic. This leads to the hypothesis that there is a “supertidal” cycle linked to the supercontinent cycle. Here this is tested using new tectonic predictions for the next 250 Myr as bathymetry in a numerical tidal model. The simulations support the following hypothesis: a new tidal resonance will appear 150 Myr from now, followed by a decreasing tide as the supercontinent forms 100 Myr later. This affects the dissipation of tidal energy in the oceans, with consequences for the evolution of the Earth‐Moon system, ocean circulation and climate, and implications for the ocean's capacity of hosting and evolving life.

Description: The Antarctic Circumpolar Current (ACC) has highly energetic mesoscale phenomena, but their impacts on phytoplankton biomass, productivity, and biogeochemical cycling are not understood well. We analyze satellite observations and an eddy‐rich ocean model to show that they drive chlorophyll anomalies of opposite sign in winter versus summer. In winter, deeper mixed layers in positive sea surface height (SSH) anomalies reduce light availability, leading to anomalously low chlorophyll concentrations. In summer with abundant light, however, positive SSH anomalies show elevated chlorophyll concentration due to higher iron level, and an iron budget analysis reveals that anomalously strong vertical mixing enhances iron supply to the mixed layer. Features with negative SSH anomalies exhibit the opposite tendencies: higher chlorophyll concentration in winter and lower in summer. Our results suggest that mesoscale modulation of iron supply, light availability and vertical mixing plays an important role in causing systematic variations in primary productivity over the seasonal cycle.

Description: Halogen- and sulfur-containing compounds are supersaturated in the surface ocean, which results in their emission to the atmosphere. These compounds can be transported to the stratosphere, where they impact ozone, the background aerosol layer, and climate. In this study we calculate the seasonal and interannual variability of transport from the West Indian Ocean (WIO) surface to the stratosphere for 2000-2016 with the Lagrangian transport model FLEXPART using ERA-Interim meteorological fields. We investigate the transport relevant for very short lived substances (VSLS) with tropospheric lifetimes corresponding to dimethylsulfide (1 day), methyl iodide (CH3I, 3.5 days), bromoform (CHBr3, 17 days), and dibromomethane (CH2Br2, 150 days). The stratospheric source gas injection of VSLS tracers from the WIO shows a distinct annual cycle associated with the Asian monsoon. Over the 16-year time series, a slight increase in source gas injection from the WIO to the stratosphere is found for all VSLS tracers and during all seasons. The interannual variability shows a relationship with sea surface temperatures in the WIO as well as the El Niño-Southern Oscillation. During boreal spring of El Niño, enhanced stratospheric injection of VSLS from the tropical WIO is caused by positive sea surface temperature anomalies and enhanced vertical uplift above the WIO. During boreal fall of La Niña, strong injection is related to enhanced atmospheric upward motion over the East Indian Ocean and a prolonged Indian summer monsoon season. Related physical mechanisms and uncertainties are discussed in this study

Description: The notion of a shallow northern sourced intermediate water mass is a well evidenced feature of the Atlantic circulation scheme of the Last Glacial Maximum (LGM). However, recent observations from stable carbon isotopes (δ13C) at the Corner Rise in the deep northwest Atlantic suggested a significant contribution of a Northern Component Water mass to the abyssal northwest Atlantic basin that has not been described before. Here we test the hypothesis of this northern sourced water mass underlying the southern sourced glacial Antarctic Bottom Water by measuring the authigenic neodymium (Nd) isotopic composition from the same sediments from 5,010-m water depth. Neodymium isotopes act as a semiconservative water mass tracer capable of distinguishing between Northern and Southern Component Waters at the northwest Atlantic. Our new Nd isotopic record resolves various water mass changes from the LGM to the early Holocene in agreement with existing Nd-based reconstructions from across the west Atlantic Ocean. Especially pronounced are the Younger Dryas and Bølling-Allerød with unprecedented changes in the Nd isotopic composition. For the LGM we found Nd isotopic evidence for a northern sourced water mass contributing to abyssal depths, thus being in agreement with previous δ13C data from Corner Rise. Overall, however, the deep northwest Atlantic was still dominated by southern sourced water, since we found signatures that are intermediate between northern and southern end member compositions. Furthermore, this new record indicates that C and Nd isotopes were partly decoupled, pointing to nonconservative behavior of one or more likely of both water mass proxies during the LGM.

Description: The MARCAN project, launched last January, is working to fill a gap in our knowledge of how freshwater flowing underground shapes and alters the continental margins.

Description: Recent evidence from mooring data in the equatorial Atlantic reveals that semi-annual and longer time scale ocean current variability is close to being resonant with equatorial basin modes. Here we show that intraseasonal variability, with time scales of 10's of days, provides the energy to maintain these resonant basin modes against dissipation. The mechanism is analogous to that by which storm systems in the atmosphere act to maintain the atmospheric jet stream. We demonstrate the mechanism using an idealised model set-up that exhibits equatorial deep jets. The results are supported by direct analysis of available mooring data from the equatorial Atlantic Ocean covering a depth range of several thousand meters. The analysis of the mooring data suggests that the same mechanism also helps maintain the seasonal variability.

Description: Reading the sediment record in terms of past climates is challenging since linking climate change to the associated responses of sedimentary systems is not always straightforward. Here we analyze the erosional response of landscapes on the Tibetan Plateau to interglacial climate forcing. Using the theory of dynamical systems on Holocene time series of geochemical proxies, we derive a sedimentary response model that accurately simulates observed proxy variation in three lake records. The model suggests that millennial variations in sediment composition reflect a self‐organization of landscapes in response to abrupt climate change between 11.6 and 11.9 ka BP. The self‐organization is characterized by oscillations in sediment supply emerging from a feedback between physical and chemical erosion processes, with estimated response times between 3,000 to 18,000 years depending on catchment topography. The implications of our findings emphasize the need for landscape response models to decipher the paleoclimatic code in continental sediment records.

Description: Extensive field campaigns in the Mauritanian upwelling region between 2005 and 2016 provide the database for analyzing the seasonal variability of the eastern boundary circulation (EBC) and associated water mass distribution at 18°N. The data set includes shipboard upper ocean current, hydrographic, and oxygen measurements from nine research cruises conducted during upwelling (December to April) and relaxation (May to July) seasons. During the upwelling season, the EBC closely resembles a classical eastern boundary current regime, with a poleward undercurrent flowing beneath an equatorward coastal jet. In contrast, elevated poleward flow exceeding 30 cm/s and extending from the surface down to 250-m depth is observed during the relaxation season. The pronounced seasonal variability of the across-shore structure of the EBC can be related to local wind forcing and is in general agreement with Sverdrup balance. The EBC transport is correlated to the wind stress curl leading the transport by 7 days. The short lead time suggests a fast response of locally forced waves adjusting the EBC to wind forcing. The seasonal and vertical water mass distribution is presented based on hydrographic observations. The meridional oxygen distribution and corresponding water mass partitioning into South and North Atlantic Central Water masses reveal a possible northerly ventilation pathway in the deeper layers of the central water stratum. Our results suggest that the poleward surface flow and the poleward undercurrent both are a consequence of the cyclonic wind stress curl forcing and thus propose to name it the Mauritania Current.

Description: The South China Sea (SCS) is the largest semienclosed marginal sea in the western Pacific (WP) and connects to the west Pacific through the Luzon Strait (LU). In this study, we use the observation of transient tracer chlorofluorocarbon‐12 (CFC‐12) to calculate the ventilation time scales of the SCS, LU, and WP. The CFC‐12 and oxygen data are used together to identify the sandwiched structure vertically of the flows across the LU. The CFC‐12 and oxygen distributions reveal a pronounced decrease westward across the LU and a slight decrease southward in the transport of the SCS. The mean age gradient of the salinity minimum (Smin) water between the WP and the northern SCS could be a consequence of intensive mixing and entrainment of the inflow water from the WP. An expected difference in age between the LU and SSCS is verified to reflect the transit time for the given water layers in the SCS. Thus, a mean transit time of 77 ± 20 years is estimated for the intermediate water in the SCS interior.

Description: The ocean is estimated to contribute up to ~20% of global fluxes of atmospheric nitrous oxide (N2O), an important greenhouse gas and ozone depletion agent. Marine oxygen minimum zones contribute disproportionately to this flux. To further understand the partition of nitrification and denitrification and their environmental controls on marine N2O fluxes, we report new relationships between oxygen concentration and rates of N2O production from nitrification and denitrification directly measured with 15N tracers in the Eastern Tropical Pacific. Highest N2O production rates occurred near the oxic‐anoxic interface, where there is strong potential for N2O efflux to the atmosphere. The dominant N2O source in oxygen minimum zones was nitrate reduction, the rates of which were 1 to 2 orders of magnitude higher than those of ammonium oxidation. The presence of oxygen significantly inhibited the production of N2O from both nitrification and denitrification. These experimental data provide new constraints to a multicomponent global ocean biogeochemical model, which yielded annual oceanic N2O efflux of 1.7–4.4 Tg‐N (median 2.8 Tg‐N, 1 Tg = 1012 g), with denitrification contributing 20% to the oceanic flux. Thus, denitrification should be viewed as a net N2O production pathway in the marine environment.

Description: Meter‐scale AUV bathymetric mapping and ROV sampling of the entire 47 km‐long Alarcon Rise between the Pescadero and Tamayo transforms show that the shallowest inflated portion of the segment hosts all four active hydrothermal vent fields and the youngest, hottest, and highest effusion rate lava flows. This shallowest inflated part is located ∼1/3 of the way between the Tamayo and Pescadero transforms and is paved by a 16 km2 channelized flow that erupted from 9 km of en echelon fissures and is larger than historic flows on the East Pacific Rise or on the Gorda and Juan de Fuca Ridges. Starting ∼5 km south of the Pescadero transform, 6.5 km of the Alarcon Rise is characterized by faulted ridges and domes of fractionated lavas ranging from basaltic andesite to rhyolite with up to 77.3 wt % SiO2. These are the first known rhyolites from the submarine global mid‐ocean ridge system. Silicic lavas range from 〉11.7 ka, to as young as 1.1 ka. A basalt‐to‐basaltic andesite sequence and an andesite‐to‐dacite‐to‐rhyolite sequence are consistent with crystal fractionation but some intermediate basaltic andesite and andesite formed by mixing basalt with dacite or rhyolite. Magmatism occurred along the bounding Tamayo and Pescadero transforms as extensive channelized flows. The flows erupted from ring faults surrounding uplifted sediment hills inferred to overlie sills. The transforms are transtensional to accommodate magma migration from the adjacent Alarcon Rise.

Description: We analyze the contribution of the Agulhas Current on the central water masses of the Benguela upwelling system (BUS) over the last decades in a high-resolution ocean simulation driven by atmospheric reanalysis. The BUS is an Eastern Boundary Upwelling System where upwelling of cold nutrient-rich water favors biomass growth. The two distinct subregions, North and South Benguela, differ in nutrient and oxygen properties of the upwelling water mass. Our analysis indicates that the contribution of Agulhas water to the upwelling is very strong in both subregions. Although the water masses feeding the upwelling have a common origin, their pathways are distinct in both regions. Whereas for the central waters of South Benguela the path is rather direct from where it is formed, the central waters of North Benguela takes a longer route through the equatorial current system. Not only the travel time from the Agulhas Current to the BUS is longer but also the central water mass is twice as old for the northern part when compared to the southern. Our analysis traces the pathways, history, and origin of the central water masses feeding upwelling in the BUS and emphasizes the direct impact of the Agulhas Current on the upwelling region. The variability of that link between the Indian Ocean and the South Atlantic is likely to change the nutrient and oxygen content, as well as temperature and salinity of the water masses in the upwelling region.

Description: 234Th measurements are widely used to estimate the downward carbon flux of particles via the oceanic Biological Pump. Carbon export is evaluated from 234Th-238U disequilibrium assuming either steady state (SS) conditions, or including a non-SS (NSS) correction. We use a novel stochastic simulation to quantify the temporal variation of vertical carbon and 234Th (dissolved and particulate) concentration profiles with high temporal resolution. We calculate seasonal export as if in situ measured with sediment trap and SS- and NSS-234Th approaches and quantify the periods of validity for SS/NSS conditions defined in previous works. The SS approach is valid throughout the entire season in oligotrophic regions. In temperate regions, the SS introduces a bias in the export if sampling takes place outside specific temporal windows. Windows of validity range from days in short blooms of ~15-day duration to weeks in blooms longer than ~30 days.

Description: The classic scenario for the generation of Dansgaard‐Oeschger (DO) events assumes a link to changes in the Atlantic Meridional Overturning Circulation (AMOC) induced by North Atlantic freshwater perturbations. Recent proxy data emphasize the existence of leads and lags between DO fingerprints in Greenland and Antarctic records, highlighting the potential of a Southern Hemisphere control on these events. Investigating this possibility, we provide a conceptual model resulting from phase space reconstructions based on the northern and southern ice core records. The resulting patterns closely resemble AMOC hysteresis, consistent with a northern abrupt warming linked to gradual global temperature changes. This suggests that rapid DO warmings associated with abrupt AMOC transitions from a relatively weak (cold stadial) state to a stronger (warm interstadial) state can be controlled by global forcing that can be linked to the Southern Hemisphere, rather than by the end of a local temporary forcing in the North Atlantic.

Description: We present a dynamically consistent gridded data set of the global, monthly mean oxygen isotope ratio of seawater ( urn:x-wiley:jgrc:media:jgrc23118:jgrc23118-math-0001). The data set was created from an optimized simulation of an ocean general circulation model constrained by global monthly urn:x-wiley:jgrc:media:jgrc23118:jgrc23118-math-0002 data collected from 1950 to 2011 and climatological salinity and temperature data collected from 1951 to 1980. The optimization was obtained using the adjoint method for variational data assimilation, which yields a simulation that is consistent with the observational data and the physical laws embedded in the model. Our data set performs equally well as a previous data set in terms of model‐data misfit but brings an improvement in terms of the seasonal cycle and physical consistency. As a result the data set does not show any sharp transitions between water masses or in areas where the data coverage is low. The data assimilation method shows high potential for interpolating sparse data sets in a physically meaningful way. Comparatively big errors, however, are found in our data set in the surface levels in the Arctic Ocean mainly because the influence of isotopically highly depleted precipitation is not preserved in the sea ice model, and the low model resolution of about 285 km horizontally. The data set is publicly available, and it is anticipated to be useful for a large range of applications in (paleo‐) oceanographic studies.

Description: The spatial distribution of subsurface parameters such as permeability are increasingly relevant for regional to global climate, land surface, and hydrologic models that are integrating groundwater dynamics and interactions. Despite the large fraction of unconsolidated sediments on Earth's surface with a wide range of permeability values, current global, high‐resolution permeability maps distinguish solely fine‐grained and coarse‐grained unconsolidated sediments. Representative permeability values are derived for a wide variety of unconsolidated sediments and applied to a new global map of unconsolidated sediments to produce the first geologically constrained, two‐layer global map of shallower and deeper permeability. The new mean logarithmic permeability of the Earth's surface is −12.7 ± 1.7 m2 being 1 order of magnitude higher than that derived from previous maps, which is consistent with the dominance of the coarser sediments. The new data set will benefit a variety of scientific applications including the next generation of climate, land surface, and hydrology models at regional to global scales.

Description: The ability of state‐of‐the‐art regional climate models to simulate cyclone activity in the Arctic is assessed based on an ensemble of 13 simulations from 11 models from the Arctic‐CORDEX initiative. Some models employ large‐scale spectral nudging techniques. Cyclone characteristics simulated by the ensemble are compared with the results forced by four reanalyses (ERA‐Interim, National Centers for Environmental Prediction‐Climate Forecast System Reanalysis, National Aeronautics and Space Administration‐Modern‐Era Retrospective analysis for Research and Applications Version 2, and Japan Meteorological Agency‐Japanese 55‐year reanalysis) in winter and summer for 1981–2010 period. In addition, we compare cyclone statistics between ERA‐Interim and the Arctic System Reanalysis reanalyses for 2000–2010. Biases in cyclone frequency, intensity, and size over the Arctic are also quantified. Variations in cyclone frequency across the models are partly attributed to the differences in cyclone frequency over land. The variations across the models are largest for small and shallow cyclones for both seasons. A connection between biases in the zonal wind at 200 hPa and cyclone characteristics is found for both seasons. Most models underestimate zonal wind speed in both seasons, which likely leads to underestimation of cyclone mean depth and deep cyclone frequency in the Arctic. In general, the regional climate models are able to represent the spatial distribution of cyclone characteristics in the Arctic but models that employ large‐scale spectral nudging show a better agreement with ERA‐Interim reanalysis than the rest of the models. Trends also exhibit the benefits of nudging. Models with spectral nudging are able to reproduce the cyclone trends, whereas most of the nonnudged models fail to do so. However, the cyclone characteristics and trends are sensitive to the choice of nudged variables.