ALBERT

Description: Key Points: - We reconstruct the temporal evolution of seawater isotope ratios of boron, strontium, lithium, and osmium over the last 65 million years - The evolution of seawater boron isotope ratio shows similarity to the evolution of strontium, lithium and osmium isotope ratios - Randomly drawn, smooth time series are provided for use in uncertainty propagation in calculation of palaeo pH The boron isotope ratio of seawater (δ11Bsw) is a parameter which must be known to reconstruct palaeo pH and CO2 from boron isotope measurements of marine carbonates. Beyond a few million years ago, δ11Bsw is likely to have been different to modern. Palaeo δ11Bsw can be estimated by simultaneously constraining the vertical gradients in foraminiferal δ11B (Δδ11B) and pH (ΔpH). A number of subtly different techniques have been used to estimate ΔpH in the past, all broadly based on assumptions about vertical gradients in oxygen, and/or carbon, or other carbonate system constraints. In this work we pull together existing data from previous studies, alongside a constraint on the rate of change of δ11Bsw from modeling. We combine this information in an overarching statistical framework called a Gaussian Process. The Gaussian Process technique allows us to bring together data and constraints on the rate of change in δ11Bsw to generate random plausible evolutions of δ11Bsw. We reconstruct δ11Bsw, and by extension palaeo pH, across the last 65Myr using this novel methodology. Reconstructed δ11Bsw is compared to other seawater isotope ratios, namely ,87/86 Sr, 187/188 Os , and δ7Li, which we also reconstruct with Gaussian Processes. Our method provides a template for incorporation of future δ11Bsw constraints, and a mechanism for propagation of uncertainty in δ11Bsw into future studies.

Description: Air-sea interaction in late boreal winter is studied over the extratropical North Atlantic (NA) during 1960–2020 by examining the relationship between sea-surface temperature (SST) and total turbulent heat flux (THF). The two quantities are positively correlated on interannual timescales over the central-midlatitude and subpolar NA, suggesting the atmosphere on average drives SST and THF variability is independent of SST. On decadal timescales and over the central-midlatitude NA the correlation is negative, suggesting ocean processes on average drive SST and THF variability is sensitive to SST. The correlation is positive over the subpolar NA. There, interannual and decadal THF variability is governed by the North Atlantic Oscillation (NAO). During the major late 20th and early 21st century SST increase in the subpolar NA diminishing oceanic heat loss associated with a weakening NAO was observed. This study suggests that the atmosphere is more sensitive to SST over the central-midlatitude than subpolar NA. Key Points: - Regional variation in the nature of air-sea interaction over the extratropical North Atlantic (NA) north of 35°N - Timescale dependence in relationship between sea-surface temperature (SST) and turbulent heat flux over the central-midlatitude NA - The atmosphere is more sensitive to SST variability over the central-midlatitude than subpolar NA

Description: Scientists are using a variety of geochemical, geophysical, and numerical methods to study offshore freshened groundwater and better understand its role in the global water cycle.

Description: Low-level jets (LLJs), vertical profiles with a wind speed maxima in the lowest hundred meters of the troposphere, have multiple impacts in the Earth system, but a global present-day climatology based on contemporary data does not exist. We use the spatially and temporally complete data set from ERA5 reanalysis to compile a global climatology of LLJs for studying the formation mechanisms, characteristics, and trends during the period of 1992–2021. In the global mean, LLJs are detected 21% of the time with more cases over land (32%) than over the ocean (15%). We classified the LLJs into three categories: non-polar land (LLLJ), polar land (PLLJ), and coastal (CLLJ) LLJs. For LLLJ, the averaged frequency of occurrence is 20% and 75% of them are associated with a near-surface temperature inversion as a prerequisite for an inertial oscillation. PLLJs are also associated with a temperature inversion and occur even more frequently with 59% of the time. These are also the lowest and the strongest LLJs among the three categories. CLLJs are particularly frequent in some marine hotspots, situated along the west coast of continents, with neutral to unstable stratification close to the surfaces and a stably stratified layer aloft. We found distinct regional trends in both the frequency and intensity of LLJs over the past decades, which can have implications for the emission and transport of aerosols, and the transport of atmospheric moisture. Future studies could address changes in LLJs and the associated implications in more detail, based on the here released ERA5-based LLJ data. Key Points: - First global comprehensive low-level jet (LLJ) climatology using ERA5 - Polar LLJs are the strongest and most frequent among the detected types - Distinct past trends in regional LLJ frequency and intensity

Description: The Banda Sea is of crucial importance for the circulation of the world's oceans, as it is part of the connection between the Pacific to the Indian Ocean. One peculiarity of the upper ocean hydrography in the Banda Sea is the occurrence of barrier layers. The regionality and temporal variability of barrier layer thickness (BLT) in the Banda Sea are examined in this study utilizing in-situ observations and ocean reanalysis output. It is found that a barrier layer occurs in over 90 % of the observational data profiles, and in over 72 % of those profiles, the BLT is shallower than 10 m. Furthermore, we find a seasonal cycle in BLT with a maximum thickness of about 60 m occurring during austral autumn and winter and coinciding with the presence of low saline waters fed by the regional river discharge and rainfall from the Java Sea and Makassar Strait. In addition, we identify the existence of a quasi-permanent anticyclonic circulation cell in the Banda Sea that may support the trapping of surface freshwater by retention. The anticyclonic circulation is most likely wind-driven because it coincides with the regional Ekman pumping pattern. Modulation of the anticyclone is via seasonal variability in the wind stress curl which in turn may explain the efficiency of freshwater retention and thus the BLT. The annual mean BLT distribution in the Banda Sea shows a preferential region of thickened barrier layers around 6o-8oS and 124o-126oE and resampling the pattern of the monthly mean climatology. Key Points: - First study estimating barrier layer thickness (BLT) in the Banda Sea using comprehensive observations - A quasi-permanent barrier layer exists in the Banda Sea with seasonal variation in occurrence and thickness - The intrusion of low saline waters and anticyclonic circulation are identified as the main mechanisms for creating and modulating the local BLT

Description: The recent severe European summer heat waves of 2015 and 2018 co-occurred with cold subpolar North Atlantic (NA) sea surface temperatures (SSTs). However, a significant connection between this oceanic state and European heat waves was not yet established. We performed two AMIP-like model experiments: (a) employing daily 2018 SSTs as observed and (b) applying a novel approach to remove the negative NA SST anomaly, while keeping SST daily and small-scale variability. Comparing these experiments, we find that cold subpolar NA SSTs significantly increase heat wave duration and magnitude downstream over the European continent. Surface temperature and circulation anomalies are connected by the upper-tropospheric summer wave pattern of meridional winds over the North Atlantic European sector, which is enhanced with cold NA SSTs. Our results highlight the relevance of the subpolar NA region for European summer conditions, a region that is marked by large biases in current coupled climate model simulations. Key Points: - Model study designed to investigate the ocean impact on European heat waves by prescribing observed and realistic ocean surface conditions - Cold subpolar North Atlantic sea surface temperatures significantly enhance heat wave intensity and duration over the European continent - North Atlantic ocean and European surface temperature and circulation anomalies are bridged by the upper-tropospheric summer mean wave

Description: In geoscience and other fields, researchers use models as a simplified representation of reality. The models include processes that often rely on uncertain parameters that reduce model performance in reflecting real-world processes. The problem is commonly addressed by adapting parameter values to reach a good match between model simulations and corresponding observations. Different optimization tools have been successfully applied to address this task of model calibration. However, seeking one best value for every single model parameter might not always be optimal. For example, if model equations integrate over multiple real-world processes which cannot be fully resolved, it might be preferable to consider associated model parameters as random parameters. In this paper, a random parameter is drawn from a wide probability distribution for every singe model simulation. We developed an optimization approach that allows us to declare certain parameters random while optimizing those that are assumed to take fixed values. We designed a corresponding variant of the well known Covariance Matrix Adaption Evolution Strategy (CMA-ES). The new algorithm was applied to a global biogeochemical circulation model to quantify the impact of zooplankton mortality on the underlying biogeochemistry. Compared to the deterministic CMA-ES, our new method converges to a solution that better suits the credible range of the corresponding random parameter with less computational effort.

Description: Flow of dense shelf water provide an efficient mechanism for pumping CO 2 to the deep ocean along the continental shelf slope, particularly around the Antarctic bottom water (AABW) formation areas where much of the global bottom water is formed. However, the contribution of the formation of AABW to sequestering anthropogenic carbon ( C ant ) and its consequences remain unclear. Here, we show prominent transport of C ant (25.0 ± 4.7 Tg C yr −1 ) into the deep ocean (〉2,000 m) in four AABW formation regions around Antarctica based on an integrated observational data set (1974–2018). This maintains a lower C ant in the upper waters than that of other open oceans to sustain a stronger CO 2 uptake capacity (16.9 ± 3.8 Tg C yr −1 ). Nevertheless, the accumulation of C ant can further trigger acidification of AABW at a rate of −0.0006 ± 0.0001 pH unit yr −1 . Our findings elucidate the prominent role of AABW in controlling the Southern Ocean carbon uptake and storage to mitigate climate change, whereas its side effects (e.g., acidification) could also spread to other ocean basins via the global ocean conveyor belt. Plain Language Summary The Southern Ocean is thought to uptake and store a large amount of anthropogenic CO 2 ( C ant ), but little attention has been paid to the Antarctic coastal regions in the south of 60°S, mainly due to the lack of observations. Based on an integrated data set, we discovered the deep penetration of C ant and a visible pattern of relatively high concentration of C ant along the AABW formation pathway, and the concentration of C ant along the shelf‐slope is higher than that of other marginal seas at low‐mid latitudes, implying a highly effective C ant transport in AABW formation areas. We also found strong upper‐layer CO 2 uptake and a significant acidification rate in the deep waters of the Southern Ocean due to the AABW‐driven CO 2 transport, which is 3 times faster than those in other deep oceans. It is therefore crucial to understand how the Antarctic shelf regions affect the global carbon cycle through the uptake and transport of anthropogenic CO 2 , which also drives acidification in the other ocean basins. Key Points We show evidence for the accumulation of C ant along the Antarctic shelf‐slope into the deep ocean The process of AABW formation drives C ant downward transport at 25.0 ± 4.7 Tg C yr −1 , sustaining the CO 2 uptake in the surface ocean This further triggers acidification of AABW at a rate of −0.0006 ± 0.0001 pH unit yr −1 , which is faster than in other deep oceans

Description: To reach their net-zero targets, countries will have to compensate hard-to-abate CO2 emissions through carbon dioxide removal (CDR). Yet, current assessments rarely include socio-cultural or institutional aspects or fail to contextualize CDR options for implementation. Here we present a context-specific feasibility assessment of CDR options for the example of Germany. We assess 14 CDR options, including three chemical carbon capture options, six options for bioenergy combined with carbon capture and storage (BECCS), and five options that aim to increase ecosystem carbon uptake. The assessment addresses technological, economic, environmental, institutional, social-cultural and systemic considerations using a traffic-light system to evaluate implementation opportunities and hurdles. We find that in Germany CDR options like cover crops or seagrass restoration currently face comparably low implementation hurdles in terms of technological, economic, or environmental feasibility and low institutional or social opposition but show comparably small CO2 removal potentials. In contrast, some BECCS options that show high CDR potentials face significant techno-economic, societal and institutional hurdles when it comes to the geological storage of CO2. While a combination of CDR options is likely required to meet the net-zero target in Germany, the current climate protection law includes a limited set of options. Our analysis aims to provide comprehensive information on CDR hurdles and possibilities for Germany for use in further research on CDR options, climate, and energy scenario development, as well as an effective decision support basis for various actors. Key Points: - More context-specific assessments of carbon dioxide removal (CDR) options are needed to guide national net-zero decision making - Ecosystem-based CDR options with comparably low implementation hurdles in Germany show relatively small CO2 removal potentials - High CDR potential options in Germany face high institutional, technological and societal hurdles linked in many ways to geological storage

Description: The Northwest Tropical Atlantic (NWTA) is a region of complex surface ocean circulation. The most prominent feature is the North Brazil Current (NBC) and its retroflection at 8°N, which leads to the formation of numerous mesoscale eddies known as NBC rings. The NWTA also receives the outflow of the Amazon River, generating freshwater plumes that can extend up to 100,000 km2. We show that these two processes influence the spatial variability of the region's surface latent heat flux (LHF). On the one hand, the presence of surface freshwater modifies the vertical stratification of the ocean, the mixed layer heat budget, and thus the air-sea heat exchanges. On the other hand, NBC rings create a highly heterogeneous mesoscale sea surface temperature (SST) field that directly influences the near-surface atmospheric circulation. These effects are illustrated by observations from the ElUcidating the RolE of Cloud-Circulation Coupling in ClimAte - Ocean Atmosphere (EUREC4A-OA) and Atlantic Tradewind Ocean-Atmosphere Interaction Campaign (ATOMIC) experiments, satellite and reanalysis data. We decompose the LHF budget into several terms controlled by different atmospheric and oceanic processes to identify the mechanisms leading to LHF changes. We find LHF variations of up to 160 W m2, of which 100 W m2 are associated with wind speed changes and 40 W m2 with SST variations. Surface currents or heat release associated with stratification changes remain as second-order contributions with LHF variations of less than 10 W m2 each. This study highlights the importance of considering these three components to properly characterize LHF variability at different spatial scales, although it is limited by the scarcity of collocated observations. Key Points: - Latent heat flux (LHF) presents strong spatial variations in the northwest tropical Atlantic (NWTA), which has a complex ocean circulation - Surface winds and sea surface temperature are the major drivers of LHF changes. The Amazon plume remains as a second-order contributor - It is necessary to distinguish between spatial scales (mesoscale and below vs. large-scale) when assessing the ocean's influence on LHF

Description: Tropical sea surface temperature (SST) biases can cause atmospheric biases on global scales, hence SST needs to be represented well in climate models. A major source of uncertainties is the representation of turbulent mixing in the oceanic boundary layer, or mixed layer (ML). In the present study we focus on near-inertial wave (NIW) induced mixing. The performance of two mixing schemes, Turbulent Kinetic Energy and K-profile parameterization (KPP), is assessed at two sites (11.5°N, 23°W and 15°N, 38°W) in the tropical Atlantic. At 11.5°N, turbulence observations (eddy diffusivities, shear and stratification) are available for comparison. We find that the schemes differ in their representation of NIWs, but both under-represent the observed enhanced diffusivities below the observed ML. However, we find that the models do mix below the ML at 15°N when a storm passes nearby. The near-inertial oscillations remain below the ML for the following 10 days. Near-inertial kinetic energy (NIKE) biases in the models are not directly correlated with the wind speed, the MLD biases, or the stratification at the ML base. Instead, NIKE biases are sensitive to the vertical mixing scheme parameterization. NIKE biases are lowest when the KPP scheme is used. Key Points: - Observations of inertial oscillations are used to evaluate the performance of two vertical mixing schemes in two high-resolution models - Both the K-profile parameterization and the Turbulent Kinetic Energy closure underestimate the NIW-induced mixing - Near-inertial kinetic energy biases are sensitive to the vertical mixing parameterization

Description: Global coupled climate models are in continuous need for evaluation against independent observations to reveal systematic model deficits and uncertainties. Changes in terrestrial water storage (TWS) as measured by satellite gravimetry missions GRACE and GRACE-FO provide valuable information on wetting and drying trends over the continents. Challenges arising from a comparison of observed and modelled water storage trends are related to gravity observations including non-water related variations such as, for example, glacial isostatic adjustment (GIA). Therefore, correcting secular changes in the Earth's gravity field caused by ongoing GIA is important for the monitoring of long-term changes in terrestrial water from GRACE in particular in former ice-covered regions. By utilizing a new ensemble of 56 individual realizations of GIA signals based on perturbations of mantle viscosities and ice history, we find that many of those alternative GIA corrections change the direction of GRACE-derived water storage trends, for example, from gaining mass into drying conditions, in particular in Eastern Canada. The change in the sign of the TWS trends subsequently impacts the conclusions drawn from using GRACE as observational basis for the evaluation of climate models as it influences the dis-/agreement between observed and modelled wetting/drying trends. A modified GIA correction, a combined GRACE/GRACE-FO data record extending over two decades, and a new generation of climate model experiments leads to substantially larger continental areas where wetting/drying trends currently observed by satellite missions coincide with long-term predictions obtained from climate model experiments.

Description: The Cabo Verde Archipelago is related to a mantle plume located close to the rotational pole of the African Plate. It consists of islands and seamounts arranged in a horseshoe‐shaped pattern open to the west, thus forming two volcanic chains, each with a weak east‐west age progression. High‐resolution swath bathymetry of 12 Cabo Verde seamounts is used here to assign each seamount to its pre‐shield, shield or post‐shield evolutionary stage, respectively. The eastern seamounts exhibit degraded and partially eroded morphologies, and are mainly in their post‐shield stage. A new 40 Ar‐ 39 Ar date for Senghor Seamount at 14.872 ± 0.027 Ma supports old ages for the eastern seamounts. The western seamounts generally exhibit younger volcanic‐edifice‐construction morphologies, showing fresh effusive and explosive volcanics, including rarely observed deep‐water explosive volcanism in the Charles Darwin Volcanic Field. Furthermore, the two previously unknown seamounts Sodade and Tavares in the westernmost termini of both volcanic chains exhibit pristine volcanic morphologies, in agreement with present‐day volcanism and seismic activity recorded from the western seamounts. The islands and seamounts rest on three submarine platforms to the east, northwest and southwest, respectively. Taken together, the seamount and island data suggest a shift in igneous activity from the eastern to the other platforms at about 8–6 Ma. However, the complex evolution pattern for both volcanic chains includes the simultaneous occurrence of pre‐shield or shield edifices at any time, followed by erosional and rejuvenation stages. The new seamount data still demonstrate ongoing westward submarine‐growth in both volcanic chains. Plain Language Summary The Cabo Verde volcanic islands and seamounts are located in the central Atlantic Ocean, ∼570 km off the west coast of Africa. They form a horseshoe‐shaped archipelago with two volcanic chains, which were formed by the African plate moving very slowly over a mantle hotspot (the Cabo Verde Plume). Both the northern and southern volcanic chains show weak east‐to‐west age progressions from ∼26 million years to the present day. This study uses underwater topographic data and observations/rock sampling via remotely operated vehicles from 12 submarine volcanic seamounts, including two previously unknown seamounts, collected during four research cruises in the Cabo Verde Archipelago. Geomorphology is used to classify each seamount as being in its pre‐shield, shield or post‐shield evolutionary stage, respectively. Cabo Verde islands and seamounts rest on three submarine morphological platforms, reflecting westward jumps of the main igneous activity, and also confirming the westward migration of the Cabo Verde hotspot beneath both volcanic chains. Key Points We present bathymetrical maps of 12, in part previously uncharted Cabo Verde seamounts Geomorphology reflects various evolutionary seamount stages and relative ages. Four older seamounts indicate late Quaternary sea level lowstands Islands and seamounts rest on three morphological platforms, indicating westward jumps of the main igneous activity

Description: Mineral dust is one of the most abundant atmospheric aerosol species and has various far-reaching effects on the climate system and adverse impacts on air quality. Satellite observations can provide spatio-temporal information on dust emission and transport pathways. However, satellite observations of dust plumes are frequently obscured by clouds. We use a method based on established, machine-learning-based image in-painting techniques to restore the spatial extent of dust plumes for the first time. We train an artificial neural net (ANN) on modern reanalysis data paired with satellite-derived cloud masks. The trained ANN is applied to cloud-masked, gray-scaled images, which were derived from false color images indicating elevated dust plumes in bright magenta. The images were obtained from the Spinning Enhanced Visible and Infrared Imager instrument onboard the Meteosat Second Generation satellite. We find up to 15% of summertime observations in West Africa and 10% of summertime observations in Nubia by satellite images miss dust plumes due to cloud cover. We use the new dust-plume data to demonstrate a novel approach for validating spatial patterns of the operational forecasts provided by the World Meteorological Organization Dust Regional Center in Barcelona. The comparison elucidates often similar dust plume patterns in the forecasts and the satellite-based reconstruction, but once trained, the reconstruction is computationally inexpensive. Our proposed reconstruction provides a new opportunity for validating dust aerosol transport in numerical weather models and Earth system models. It can be adapted to other aerosol species and trace gases. Key Points: - We present the first fast reconstruction of cloud-obscured Saharan dust plumes through novel machine learning applied to satellite images - The reconstruction algorithm utilizes partial convolutions to restore cloud-induced gaps in gray-scaled Meteosat Second Generation-Spinning Enhanced Visible and Infrared Imager Dust RGB images - World Meteorological Organization dust forecasts for North Africa mostly agree with the satellite-based reconstruction of the dust plume extent

Description: Observation‐based quantification of ocean carbon dioxide (CO 2 ) uptake relies on synthesis data sets such as the Surface Ocean CO 2 ATlas (SOCAT). However, the data collection effort has dramatically declined and the number of annual data sets in SOCATv2023 decreased by ∼35% from 2017 to 2021. This decline has led to a 65% increase (from 0.15 to 0.25 Pg C yr −1 ) in the standard deviation of seven SOCAT‐based air‐sea CO 2 flux estimates. Reducing the availability of the annual data to that in the year 2000 creates substantial bias (50%) in the long‐term flux trend. The annual mean CO 2 flux is insensitive to the seasonal skew of the SOCAT data and to the addition of the lower accuracy data set available in SOCAT. Our study highlights the need for sustained data collection and synthesis, to inform the Global Carbon Budget assessment, the UN‐led climate negotiations, and measurement, reporting, and verification of ocean‐based CO 2 removal projects. Plain Language Summary The Surface Ocean CO 2 ATlas (SOCAT) data set plays a crucial role in estimating the ocean carbon sink component of the Global Carbon Budget. However, the number of data sets available in SOCAT each year has drastically decreased since 2017. This study shows that the uncertainty in the data‐based ocean CO 2 flux estimate has increased by 65% due to this decline in data availability. The estimated fluxes, especially the long‐term flux trend, are remarkably affected by the data availability in SOCAT, reducing the reliability of ocean CO 2 uptake estimates in years and regions with sparse observations. Key Points Lower surface ocean f CO 2 data availability leads to higher uncertainty in data‐based estimates of ocean CO 2 uptake The long‐term trend in the ocean CO 2 flux increases by 1.5 times for subsequent years if the data availability is reduced to that in 2000 The annual mean CO 2 flux is not sensitive to the seasonal skew in the data and to the addition of low accuracy data

Description: Freshwater input from Greenland ice sheet melt has been increasing in the past decades from warming temperatures. To identify the impacts from enhanced meltwater input into the subpolar North Atlantic from 1997 to 2021, we use output from two nearly identical simulations in the eddy-rich model VIKING20X (1/20°) only differing in the freshwater input from Greenland: one with realistic interannually varying runoff increasing in the early 2000s and the other with climatologically (1961–2000) continued runoff. The majority of the additional freshwater remains within the boundary current enhancing the density gradient toward the warm and salty interior waters yielding increased current velocities. The accelerated boundary current shows a tendency to enhanced, upstream shifted eddy shedding into the Labrador Sea interior. Further, the experiments allow to attribute higher stratification and shallower mixed layers southwest of Greenland and deeper mixed layers in the Irminger Sea, particularly in 2015–2018, to the runoff increase in the early 2000s. Key Points The West Greenland Current (WGC) freshens and cools with the observed recent increase in meltwater runoff from Greenland The density gradient across the boundary current intensifies, strengthening the WGC and increasing local eddy formation Enhanced meltwater runoff contributed to an eastward shift in deep convection towards the Irminger Sea (2015–2018) Plain Language Summary Global warming has accelerated the melting of the Greenland ice sheet over the past few decades resulting in enhanced freshwater input into the North Atlantic. The additional freshwater can potentially inhibit deep water formation and have future implications on ocean circulation. To determine the influence from Greenland melt, we compare two high-resolution model experiments all with the same forcing but differing input of Greenland freshwater fluxes from 1997 to 2021. We find that in the experiment with realistically increasing Greenland meltwater, the water becomes fresher and cooler along the continental shelf and boundary of the subpolar gyre. The density difference between the shelf and interior increases with more freshwater, resulting in faster West Greenland Current speeds and enhanced eddy formation. Deeper mixed layers are found in the eastern Irminger Sea, particularly in 2015–2018. From 2009 to 2013, there were shallower mixed layers in the Labrador Sea where less Greenland meltwater was mixed downwards and spread eastward, causing mixed layers to deepen in the Irminger Sea.

Description: We report the first catalog of low‐frequency earthquakes in the Hikurangi subduction zone, located beneath the Kaimanawa Range of the North Island at 50 km depth, downdip of regularly recurring (every 4–5 years) deep M7 slow slip events. To systematically detect low‐frequency earthquakes within the regional continuous seismic data, we utilized a matched‐filter approach with template waveforms derived from previous observations of tectonic tremor. We built our catalog of 36 low‐frequency earthquake sources, that produced almost 21,000 events over more than a decade, with two matched‐filter search iterations. In each iteration, the detections were gathered into families and their coherent waveforms processed and stacked to extract high‐quality waveforms, allowing us to pick seismic phase arrivals to locate the low‐frequency earthquakes. We highlight three characteristic features to validate that our detected events are indeed low‐frequency earthquakes: the eponymous deficit of high frequencies in their seismic waveforms, the episodic swarms of activity that define their activity through time, and their location at the plate boundary with a double‐couple source mechanism and geometry consistent with the subduction interface. Considering the observed low‐frequency earthquakes' relationship to neighboring slow slip, we observe the event swarms to occur much more frequently than the M7 slow slip events located just updip. Similar to other deep low‐frequency earthquakes in other subduction zones, we suggest that this characteristic clustering in time is driven by more frequent, smaller slow slip events that are not clearly observable at the surface. Plain Language Summary Slow slip is episodic fault slip that lasts days, weeks or months, rather than the rapid ruptures of regular earthquakes. Geodetic observations of the surface displacement produced by slow slip suggest that their timing and location influence the seismic cycle of nearby faults and may even trigger large earthquakes. Although slow slip does not produce seismic radiation itself, slow slip is often accompanied by tiny repetitive seismic signals. These tiny seismic events, called low‐frequency earthquakes, can act as a powerful indicator of when and where slow slip is happening. In this study, we develop a new approach to detect low‐frequency earthquakes within continuous seismic waveforms, revealing the first observations of low‐frequency earthquakes in the Hikurangi subduction zone beneath the North Island of New Zealand. Our catalog of low‐frequency earthquakes suggests a complex pattern of slow fault slip at depth, with more frequent activity than geodetic data alone would suggest. The observed low‐frequency earthquake activity in the Hikurangi subduction zone thus represents a unique opportunity to study the slip history at depth beneath the North Island of New Zealand. Key Points 36 low‐frequency earthquake sources are extracted from continuous waveforms through template matching, deblurring, and unsupervised learning Low‐frequency earthquake sources locate close to the plate boundary with source mechanisms consistent with the subduction interface Detected low‐frequency earthquakes are likely triggered by small, frequent, and deep slow slip not geodetically observable at the surface

Description: The potential for future earthquakes on faults is often inferred from inversions of geodetically derived surface velocities for locking on faults using kinematic models such as block models. This can be challenging in complex deforming zones with many closely spaced faults or where deformation is not readily described with block motions. Furthermore, surface strain rates are more directly related to coupling on faults than surface velocities. We present a methodology for estimating slip deficit rate directly from strain rate and apply it to New Zealand for the purpose of incorporating geodetic data in the 2022 revision of the New Zealand National Seismic Hazard Model. The strain rate inversions imply slightly higher slip deficit rates than the preferred geologic slip rates on sections of the major strike‐slip systems including the Alpine Fault, the Marlborough Fault System and the northern part of the North Island Fault System. Slip deficit rates are significantly lower than even the lowest geologic estimates on some strike‐slip faults in the southern North Island Fault System near Wellington. Over the entire plate boundary, geodetic slip deficit rates are systematically higher than geologic slip rates for faults slipping less than one mm/yr but lower on average for faults with slip rates between about 5 and 25 mm/yr. We show that 70%–80% of the total strain rate field can be attributed to elastic strain due to fault coupling. The remaining 20%–30% shows systematic spatial patterns of strain rate style that is often consistent with local geologic style of faulting. Plain Language Summary The potential for future earthquakes on faults is often inferred from velocities of the ground surface derived from satellite geodesy, but this approach can be challenging in complex deforming zones with many closely spaced faults. We present a new methodology for estimating the rate at which energy is accumulating on faults using measurements of surface strain rates. The method is applied to New Zealand for the purpose of incorporating geodetic data in the 2022 revision of the New Zealand National Seismic Hazard Model. We show that 70%–80% of the total deformation field can be attributed to energy accumulation on known active faults while the source of the remaining 20%–30% remains unknown. Along some of the major faults in New Zealand we find some important differences in rates of energy accumulation from what is expected from geologic data. Estimated rates are significantly lower than even the lowest geologic estimates on some faults in the fault system near highly‐populated Wellington. Key Points We develop a method to invert geodetically derived strain rates for slip deficit rates on faults We find small but systematic differences between slip deficit rates and geologic slip rates About 70%–80% of the surface strain can be attributed to elastic strain due to coupling on faults

Description: Oceanic detachment faulting, a major mode of seafloor accretion at slow and ultraslow spreading ridges, is thought to occur during magma‐poor phases and be abandoned when magmatism increases. In this framework, detachment faulting is the result of temporal variations in magma flux, which is inconsistent with recent geophysical observations at the Longqi segment on the Southwest Indian Ridge (49°42′E). In this paper, we focus on this sequentially active detachment faulting system that includes an old, inactive detachment fault and a younger, active detachment fault. We investigate the mechanisms controlling the temporal evolution of this tectonomagmatic system by using 2D mid‐ocean ridge spreading models that simulate faulting and magma intrusion into a visco‐elasto‐plastic continuum. Our models show that temporal variations in magma flux alone are insufficient to match the inferred temporal evolution of the sequentially active detachment system. Rather we find that sequentially active detachment faulting spontaneously occurs at the Longqi segment as a function of lithospheric thickness. This finding is in agreement with an analytical model, which shows that a thicker axial lithosphere results in a smaller fault heave and that a flatter angle in lithosphere thickening away from the accretion axis stabilizes the active fault. A thicker axial lithosphere and its flatter off‐axis angle combined have the potential to modulate sequentially active detachment faulting at the Longqi segment. Our results thus suggest that temporal changes of magmatism are not necessary for the development and abandonment of detachment faults at ultraslow spreading ridges.

Description: Mesoscale eddies are common in the subtropical Northwest Pacific, however, relatively little is known about their spatial variability and temporal evolution, and how these impact upper ocean biogeochemistry. Here we investigate these using observations of a cyclonic eddy carried out along four sequential transects. Consistent with previous observations of cyclonic eddies, the eddy core had doming isopycnals, bringing elevated nutrient waters nearer to the surface. However, we also found that the upper layer of the eddy above the nutricline had significantly lower phosphate concentrations within its core relative to its edge. We attributed this to elevated N 2 fixation within the eddy core, which was likely driven by enhanced subsurface iron supply, ultimately resulting in increased phosphate consumption. Eddy‐enhanced N 2 fixation was additionally supported by the elevation of nitrate + nitrite to phosphate ratios below the euphotic zone. Moreover, we observed that while the upward displacement of isopycnals within the eddy core led to an increase in phytoplankton biomass in the lower euphotic zone, there was no significant increase in total phytoplankton biomass across the entire euphotic zone. Cyclonic eddies in the subtropical North Pacific are projected to be becoming more frequent, implying that such dynamics could become increasingly important for regulating nutrient biogeochemistry and ultimately productivity of the region. Key Points Lower phosphate concentrations were observed above the nutricline within the eddy core in comparison to the edge Enhanced N2 fixation within the eddy core is proposed to have driven increased phosphate consumption No substantial total phytoplankton biomass increase was found within the eddy core

Description: The marine cyanobacterium Trichodesmium has the remarkable ability to interact with and utilize air‐borne dust as a nutrient source. However, dust may adversely affect Trichodesmium through buoyancy loss and exposure to toxic metals. Our study explored the effect of desert dust on buoyancy and mortality of natural Red Sea puff‐shaped Trichodesmium thiebautii . Sinking velocities and ability of individual colonies to stay afloat with increasing dust loads were studied in sedimentation chambers. Low dust loads of up to ∼400 ng per colony did not impact initial sinking velocity and colonies remained afloat in the chamber. Above this threshold, sinking velocity increased linearly with the colony dust load at a slope matching prediction based on Stoke's law. The potential toxicity of dust was assessed with regards to metal dissolution kinetics, differentiating between rapidly released metals, which may impact surface blooms, and gradually released metals that may impact dust‐centering colonies. Incubations with increasing dust concentrations revealed colony death, but the observed lethal dose far exceeded dust concentrations measured in coastal and open ocean systems. Removal of toxic particles as a mechanism to reduce toxicity was explored using SEM‐EDX imaging of colonies incubated with Cu‐minerals, yet observations did not support this pathway. Combining our current and former experiments, we suggest that in natural settings the nutritional benefits gained by Trichodesmium via dust collection outweigh the risks of buoyancy loss and toxicity. Our data and concepts feed into the growing recognition of the significance of dust for Trichodesmium 's ecology and subsequently to ocean productivity. Plain Language Summary Trichodesmium spp. are abundant cyanobacteria, forming extensive blooms in low latitude warm oceans, and contribute significantly to carbon (C) and nitrogen (N) fixation, recycling and export. Desert dust deposited on the ocean surface was shown to supply Trichodesmium with the scarce micronutrient iron. Spherical, millimeter‐sized colonies of Trichodesmium from different ocean basins were reported to actively accumulate dust in their cores. While dust accumulation likely helps Trichodesmium obtain nutrients, it may come at a cost. Metals released from dust may induce toxicity and the dust weight could send Trichodesmium to the ocean depth. Our experimental study with natural Red Sea colonies examined some trade‐offs of dust accumulation. Links between dust load and colony buoyancy were examined in sedimentation experiments. Toxicity thresholds for surface blooms and dust‐accumulating colonies were determined from mortality assays and dust dissolution measurements. We found that metal‐induced toxicity to Trichodesmium is unlikely at typical oceanic dust fluxes, and that dust‐containing colonies can remain buoyant. At high loads, dust weight determined the colony's sinking velocity. Our findings and concepts can be extended to additional aerosols and Trichodesmium ‐rich habitats, and may assist in assessing Trichodesmium 's distribution, ecophysiology, and contribution to C or N transport to the deep ocean. Key Points Dust collected by Trichodesmium colonies from seawater as a nutrient source may result in metal toxification and buoyancy loss At moderate dust loads, colonies kept their buoyancy, but above 400 ng, sinking velocities increased linearly with dust loads Desert dust induced Trichodesmium mortality through toxic metal release, yet the lethal dose far exceeded oceanic dust concentrations

Description: In the northeastern tropical Atlantic, a region of high potential vorticity (PV) determines the size of the exchange window for the interior thermocline flow of the subtropical cell via its variations in strength and extent. Variability of this PV barrier has the potential to impact the ventilation of the tropical Atlantic on decadal timescales. Here, the impact of the North Atlantic Oscillation (NAO) on the PV barrier related to isopycnals within the thermocline of the subtropical-tropical Atlantic Ocean is assessed from Argo observations for the time period of 2006-2022. Relative to the negative NAO phase (2009-2010), during the positive NAO phase (2014-2019), the North Atlantic subtropical high and the northeast trades are intensified. Satellite-derived wind stress curl shows increased upwelling/downwelling on the equatorward/poleward side of the trade wind zone, respectively. In the subtropical-tropical Atlantic, a symmetric pattern of isopycnal heave is observed: rising isopycnals within 20 degrees N and 20 degrees S and sinking poleward of that. With rising isopycnals, the PV barrier in the northeastern tropical Atlantic becomes stronger. Analyses of geostrophic velocities and the Sverdrup streamfunction show that during the positive NAO phase there are increased equatorward velocities at thermocline level along the western boundary and reduced velocities through the interior as a result of intensified northeast trades and therefore a strengthened PV barrier. Intensified trades lead to enhanced subduction of thermocline waters and, independent of that, to a strengthened Equatorial Undercurrent transport as observed at the mooring site at 0 degrees, 23 degrees W, likely via the pulling effect of the subtropical cells. In the North Atlantic Ocean, subducted water from the subtropics has two possible pathways within the thermocline toward the equatorial region: the interior pathway and the pathway along the western boundary. The size of the exchange window between subtropics and tropics depends on the extent of a barrier zone in the eastern part of the basin that is associated with wind-driven upwelling of density surfaces. The North Atlantic Oscillation (NAO) is the dominant atmospheric climate mode in the North Atlantic and in this study, we show how the NAO impacts the barrier for the equatorward thermocline flow in the tropical Atlantic Ocean. During positive NAO phases (e.g., 2014-2019), density surfaces become shallower and strengthen the barrier, while during negative NAO phases (e.g., 2009-2010) the barrier weakens. Geostrophic velocity analysis reveals that during positive NAO phases more thermocline water is transported equatorward via the western boundary and less via the interior pathway. Additionally, observations from a mooring site at 0 degrees, 23 degrees W show stronger Equatorial Undercurrent transport as a result of intensified trade winds during positive NAO phases. Trade winds in the northeastern tropical Atlantic strengthen during positive phases of the North Atlantic Oscillation (NAO+) Potential vorticity barrier for the interior equatorward thermocline flow of the North Atlantic Subtropical Cell strengthens during NAO+ Annual subduction of thermocline water and Equatorial Undercurrent transport increase simultaneously from 2008 to 2018

Description: Viscosity in the momentum equation is needed for numerical stability, as well as to arrest the direct cascade of enstrophy at grid scales. However, a viscous momentum closure tends to over-dissipate eddy kinetic energy. To return excessively dissipated energy to the system, the viscous closure is equipped with what is called dynamic kinetic energy backscatter. The amplitude of backscatter is based on the amount of unresolved kinetic energy (UKE). This energy is tracked through space and time via a prognostic equation. Our study proposes to add advection of UKE by the resolved flow to that equation to explicitly consider the effects of nonlocality on the subgrid energy budget. UKE can consequently be advected by the resolved flow before it is reinjected via backscatter. Furthermore, we suggest incorporating a stochastic element into the UKE equation to account for missing small-scale variability, which is not present in the purely deterministic approach. The implementations are tested on two intermediate complexity setups of the global ocean model FESOM2: an idealized channel setup and a double-gyre setup. The impacts of these additional terms are analyzed, highlighting increased eddy activity and improved flow characteristics when advection and carefully tuned, stochastic sources are incorporated into the UKE budget. Additionally, we provide diagnostics to gain further insights into the effects of scale separation between the viscous dissipation operator and the backscatter operator responsible for the energy injection. Oceanic swirls or "eddies" have a typical size of 10-100 km, which is close to the smallest scales that global ocean models commonly resolve. For physical and numerical reasons, these models require the addition of artificial terms that influence the flow near its smallest scales. Common approaches have the drawback of introducing systematic loss of kinetic energy contained in the eddies, which leads to errors that also affect the oceanic circulation on global scales. In our research, we compensate for this error by returning some of the missing energy back into the simulation, using a so-called kinetic energy backscatter scheme. In this work, we continue the development of an already existing and successful backscatter scheme, adding certain improvements to the way energy is budgeted and returned to the flow: we ensure that the local energy budget is attached to each fluid parcel as it is transported by the large-scale flow, and we also add a random forcing term that mimics unknown sources of such energy to bring its statistical properties closer to reality. We demonstrate that these modifications effectively improve the characteristics of the simulated flow. Extension of the subgrid energy equation of the kinetic energy backscatter parameterization by adding advection and a stochastic term Both additional terms improve several flow characteristics in two idealized test cases, a channel and a double-gyre Scale analysis reveals the necessity of sufficient scale separation between viscous energy dissipation and energy injection via backscatter. Key Points: - Extension of the subgrid energy equation of the kinetic energy backscatter parameterization by adding advection and a stochastic term - Both additional terms improve several flow characteristics in two idealized test cases, a channel and a double-gyre - Scale analysis reveals the necessity of sufficient scale separation between viscous energy dissipation and energy injection via backscatter

Description: We present a continuous ∼6.2 Ma long record of explosive activity from the Northwest Pacific volcanic arcs based on a composite tephra sequence derived from Ocean Drilling Program Sites 882A and 884B, and core MD01‐2416 on the Detroit Seamount. Geochemical fingerprinting of tephra glass using major and trace element analyses and correlations of tephra layers between the three cores allowed the identification of 119 unique tephras, suggesting eruptions of magnitude (M) of 5.8–7.8. Age estimates for all the identified eruptions were obtained with the help of published and further refined age models for the studied cores, direct 40 Ar/ 39 Ar dating of four ash layers, and Bayesian age modeling. The glass compositions vary from low‐ to high‐K 2 O basaltic andesite to rhyolite and exhibit typical subduction‐related affinity. The majority of the tephras originated from Kamchatka, only a few tephras—from the neighboring Kuril and Aleutian arcs. The glass compositions revealed no temporal trends but made it possible to identify their source volcanic zones in Kamchatka and, in some cases, to determine their source eruptive centers. Our data indicates episodes of explosive activity recorded in the Detroit tephra sequence at ∼6,200, 5,600–5,000, 4,300–3,700 ka, and almost continuous activity since ∼3,000 ka. Within the latter episode, the most active intervals can be identified at 1,700–1,600, 1,150–1,050, and 600–50 ka. Geochemically fingerprinted and dated Detroit tephra sequence form a framework for dating and correlating diverse paleoenvironmental archives across the Northwest Pacific and for studies of geochemical evolution of the adjacent volcanic arcs. Plain Language Summary Explosive volcanic eruptions produce defragmented material named tephra, which can be spread over large distances and form layers in sediments on ocean floor and continents. Long continuous tephra sequences preserved in marine sediments provide one of the best chronicles of the explosive eruptions, and allow detailed evaluation of their timing relative to climatic changes. We studied one of such natural records of explosive volcanism preserved in the sediments covering the Detroit Seamount in the Northwest Pacific. We identified 119 tephra layers, which have been buried in the sediments during the last 6.2 Ma and represent volcanic eruptions with ≥7 km 3 tephra volume. We analyzed geochemical composition and determined age of each tephra. Most tephras were found to originate from volcanoes in Kamchatka, a few from the Kuril and Aleutian volcanoes. We found that the explosive activity recorded in the Detroit tephra sequence was not uniform over time. It peaked at ∼6,200, 5,600–5,000, 4,300–3,700, has continued since ∼3,000 thousand years ago until present. All tephra layers from our study can be used as unique isochrons for dating and correlating paleoenvironmental archives across the Northwest Pacific and for the reconstruction of the detailed volcanic record in the Earth history. Key Points We report age and composition for 119 tephras from sediment cores representing ∼6.2 Ma record of explosive volcanism in the NW Pacific The tephras have subduction‐related origin and mostly originate from volcanic eruptions with magnitude (M) of 5.8–7.8 in Kamchatka The data indicates episodes of explosive activity at ∼6,200, 5,600–5,000, 4,300–3,700 ka, and almost continuous activity since ∼3,000 ka

Description: Due to the complexity of 2D magnetic anomaly maps north of 18°S and the sparsity of seismic data, the tectonic evolution of the northern Lau Basin has not yet been unraveled. We use a multi-method approach to reconstruct the formation of the basin at ∼16°S by compiling seismic, magnetic, gravimetric and geochemical data along a 185 km-long crustal transect. We identified a crustal zonation which preserves the level of subduction input at the time of the crust's formation. Paired with the seafloor magnetization, the crustal zonation enabled us to qualitatively approximate the dynamic spreading history of the region. Further assessment of the recent tectonic activity and the degree of tectonic overprinting visible in the crust both suggest a complex tectonic history including a dynamically moving spreading center and the reorganizing of the local magma supply. Comparing the compiled data sets has revealed substantial differences in the opening mechanisms of the two arms of the Overlapping Spreading Center (OSC) that is made up by the northernmost tip of the Fonualei Rift and Spreading Center in the east and the southernmost segment of the Mangatolu Triple Junction in the west. The observed transition from a predominantly tectonic opening mechanism at the eastern OSC arm to a magmatic opening mechanism at the western OSC arm coincides with an equally sharp transition from and strongly subduction influenced crust to a crust with virtually no subduction input. The degree of subduction input alters the geochemical composition, as well as the lithospheric stress response. Key Points Oceanic crust in the north-eastern Lau Basin formed at the now reorganized FRSC-MTJ system The position and the opening mechanisms of back-arc basin spreading center's change more dynamically at mid-ocean ridges Different opening mechanisms at the southern Mangatolu Triple Junction and northern Fonualei Rift Spreading Center despite their proximity

Description: Submarine groundwater discharge (SGD) is a globally important process supplying nutrients and trace elements to the coastal environment, thus playing a pivotal role in sustaining marine primary productivity. Along with nutrients, groundwater also contains allochthonous microbes that are discharged from the terrestrial subsurface into the sea. Currently, little is known about the interactions between groundwater‐borne and coastal seawater microbial populations, and groundwater microbes' role upon introduction to coastal seawater populations. Here, we investigated seawater microbial abundance, activity and diversity in a site strongly influenced by SGD. In addition, through laboratory‐controlled bottle incubations, we mimicked different mixing scenarios between groundwater and seawater. Our results demonstrate that the addition of 0.1 μm filtered groundwater stimulated heterotrophic activity and increased microbial abundance compared to control coastal seawater, whereas 0.22 μm filtration treatments induced primary productivity and Synechococcus growth. 16S rRNA gene sequencing showed a strong shift from a SAR11‐rich community in the control samples to Rhodobacteraceae dominance in the 〈0.1 μm treatment, in agreement with Rhodobacteraceae enrichment in the SGD field site. These results suggest that microbes delivered by SGD may affect the abundance, activity and diversity of intrinsic microbes in coastal seawater, highlighting the cryptic interplay between groundwater and seawater microbes in coastal environments, which has important implications for carbon cycling. Plain Language Summary Submarine groundwater discharge (SGD) is an important process where groundwater flows into the ocean along the coast. When the groundwater mixes with seawater, the microbes from both sources interact with each other, which can impact the diversity, activity, and amount of microbes in the coastal environment. Currently, little is known about how groundwater‐borne microbes affect marine microbial populations. Our research shows that when groundwater microbes are removed before mixing groundwater with seawater, the abundance and activity of certain microbes that consume organic matter significantly increase. Additionally, we noticed a significant difference in the types of microbes present between the sites where SGD occurs versus background (uninfluenced) coastal water, especially in terms of the microbes that consume organic matter. Overall, this study suggests that there is a connection between groundwater and seawater microbes, which can influence the delicate balance between organisms that produce carbon and those that consume it. This has important implications for how carbon cycles globally. Key Points Groundwater discharge into the coastal zone delivers both nutrients and allochthonous microbes Groundwater microbes interact with seawater populations, by which affecting the delicate autotroph‐heterotroph balance Subterranean microbial processes are key drivers of food webs, potentially affecting biogenic carbon fluxes in the ocean

Description: With almost 700 Pg of carbon, marine dissolved organic carbon (DOC) stores more carbon than all living biomass on Earth combined. However, the controls behind the persistence and the spatial patterns of DOC concentrations on the basin scale remain largely unknown, precluding quantitative assessments of the fate of this large carbon pool in a changing climate. Net removal rates of DOC along the overturning circulation suggest lifetimes of millennia. These net removal rates are in stark contrast to the turnover times of days to weeks of heterotrophic microorganisms, which are the main consumers of organic carbon in the ocean. Here, we present a dynamic “MICrobial DOC” model (MICDOC) with an explicit representation of picoheterotrophs to test whether ecological mechanisms may lead to observed decadal to millennial net removal rates. MICDOC is in line with 〉40,000 DOC observations. Contrary to other global models, the reactivity of DOC fractions is not prescribed, but emerges from a dynamic feedback between microbes and DOC governed by carbon and macronutrient availability. A colimitation of macronutrients and organic carbon on microbial DOC uptake explains 〉70% of the global variation of DOC concentrations, and governs characteristic features of its distribution. Here, decadal to millennial net removal rates emerge from microbial processes acting on time scales of days to weeks, suggesting that the temporal variability of the marine DOC inventory may be larger than previously thought. With MICDOC, we provide a foundation for assessing global effects on DOC related to changes in heterotrophic microbial communities in a future ocean Plain Language Summary The ocean stores more carbon as dissolved organic compounds (DOC) than all animals and plants on land and the oceans combined. However, numerical models used for future climate scenarios lack an implementation of processes transforming DOC back to CO 2 by marine microorganisms. Here, we present a global dynamical ocean model that explicitly considers the processes of DOC degradation by marine microorganisms. In the present ocean, the availability of organic carbon but also nitrogen and phosphorus control the amount of carbon stored as DOC, as the lack of these nutrients inhibits its degradation by bacteria. The identification of these ecological controls allows a quantitative assessment of the fate of this large carbon reservoir in the future. The findings indicate that the marine DOC reservoir is potentially more dynamic than previously thought, since decadal to millennial scale net removal rates might be a result of microbial processes acting on shorter time scales Key Points A model to reconcile millennial‐scale bulk dissolved organic carbon degradation rates and short‐term microbial turnover times is presented Macronutrient colimitation can explain observed concentration patterns of dissolved organic carbon in the surface ocean Continuous microbial reworking suggests a higher temporal variability of the marine dissolved organic matter inventory than previously thought

Description: The study of offshore freshened groundwater (OFG) is gaining importance due to population growth and environmental pressure on coastal water resources. Marine controlled source electromagnetic (CSEM) methods can effectively map the spatial extent of OFG systems using electrical resistivity as a proxy. Integrating these resistivity models with sub-surface properties, such as host-rock porosity, allows for estimates of pore-water salinity. However, evaluating the uncertainty in pore-water salinity using resistivity models obtained from deterministic inversion approaches presents challenges, as they provide only one best-fit model, with no associated estimate of uncertainty. To address this limitation, we employ trans-dimensional Markov-Chain Monte-Carlo inversion on marine time-domain CSEM data, acquired in the Canterbury Bight, New Zealand. We integrate resistivity posterior probability distributions with borehole and seismic reflection data to estimate pore-water salinity with corresponding uncertainty estimates. The results highlight a low-salinity groundwater body in the center of the survey area, hosted by consecutive silty- and fine-sand layers approximately 20–60 km from the coast. The posterior probability distribution of resistivity models indicates freshening of the OFG body toward the shoreline within a permeable, coarse-sand layer 40–150 m beneath the seafloor, suggesting an active connection between the OFG body and the terrestrial groundwater system. The approach demonstrates how Bayesian inversion constrains the uncertainties in resistivity models and subsequently in pore-water salinity estimates. Our findings highlight the potential of Bayesian inversion to enhance our understanding of OFG systems and provide uncertainty constraints for hydrogeological modeling, thereby contributing to sustainable water resource development. Key Points A Bayesian workflow is employed to evaluate uncertainty in pore-water salinity estimates Offshore groundwater in Canterbury Bight stores freshened pore-water in fine-grained sediments, likely extending from the onshore aquifer Correlation between pore-water salinities and seismic-derived stratigraphy provides boundary conditions for hydrogeological modeling Plain Language Summary Geophysical methods that measure the electromagnetic properties of the Earth are effective in investigating freshwater sources beneath the seafloor. By combining the geophysical and geological information, we can better assess the quality of this groundwater. In this study, we develop a workflow that uses statistical methods to integrate electromagnetic observations with borehole and acoustic measurements along the eastern coast of the South Island of New Zealand. We aim to improve our understanding of the groundwater quality beneath the seafloor. Our research confirms the presence of freshened groundwater within the sandy seafloor up to 60 km from the coastline. Importantly, our observations indicate that the groundwater quality increases toward the coast. These findings are significant as they enhance the hydrogeological modeling of the groundwater system and suggest its potential as a source of freshwater.

Description: Detecting phase arrivals and pinpointing the arrival times of seismic phases in seismograms is crucial for many seismological analysis workflows. For land station data, machine learning methods have already found widespread adoption. However, deep learning approaches are not yet commonly applied to ocean bottom data due to a lack of appropriate training data and models. Here, we compiled an extensive and labeled ocean bottom seismometer (OBS) data set from 15 deployments in different tectonic settings, comprising ∼90,000 P and ∼63,000 S manual picks from 13,190 events and 355 stations. We propose PickBlue, an adaptation of the two popular deep learning networks EQTransformer and PhaseNet. PickBlue joint processes three seismometer recordings in conjunction with a hydrophone component and is trained with the waveforms in the new database. The performance is enhanced by employing transfer learning, where initial weights are derived from models trained with land earthquake data. PickBlue significantly outperforms neural networks trained with land stations and models trained without hydrophone data. The model achieves a mean absolute deviation of 0.05 s for P-waves and 0.12 s for S-waves, and we apply the picker on the Hikurangi Ocean Bottom Tremor and Slow Slip OBS deployment offshore New Zealand. We integrate our data set and trained models into SeisBench to enable an easy and direct application in future deployments. Key Points We assembled a database of ocean Bottom Seismometer (OBS) waveforms and manual P and S picks, on which we train PickBlue, a deep learning picker Our picker significantly outperforms pickers trained with land-based data with confidence values reflecting the likelihood of outlier picks The picker and database are available in the SeisBench platform, allowing easy and direct application to OBS traces and hydrophone records

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: The Hikurangi Margin east of New Zealand's North Island hosts an extensive gas hydrate province with numerous gas hydrate accumulations related to the faulted structure of the accretionary wedge. One such hydrate feature occurs in a small perched upper‐slope basin known as Urutī Basin. We investigated this hydrate accumulation by combining a long‐offset seismic line (10‐km‐long receiver array) with a grid of high‐resolution seismic lines acquired with a 600‐m‐long hydrophone streamer. The long‐offset data enable quantitative velocity analysis, while the high‐resolution data constrain the three‐dimensional geometry of the hydrate accumulation. The sediments in Urutī Basin dip landward due to ongoing deformation of the accretionary wedge. These strata are clearly imaged in seismic data where they cross a distinct bottom simulating reflection (BSR) that dips counterintuitively in the opposite direction to the regional dip of the seafloor. BSR‐derived heat flow estimates reveal a distinct heat flow anomaly that coincides spatially with the upper extent of a landward‐verging thrust fault. We present a conceptual model of this gas hydrate system that highlights the roles of fault‐controlled fluid flow at depth merging into strata‐controlled fluid flow into the hydrate stability zone. The result is a layer‐constrained accumulation of concentrated gas hydrate in the dipping strata. Our study provides new insight into the interplay between deep faulting, fluid flow and gas hydrate formation within an active accretionary margin. Plain Language Summary Gas hydrates are ice‐like substances in which natural gas molecules are trapped in a cage of water molecules. They exist where the pressure is high, temperature is cold, and enough methane is present. These conditions exist in the marine environment at water depths greater than 300–500 m near sediment‐rich continental margins and in polar regions. It is important to study gas hydrates because they represent a significant part of the Earth's carbon budget and influence the flow of methane into the oceans and atmosphere. In this study, we use the seismic reflection method to generate images of gas‐hydrate‐bearing marine sediments east of New Zealand. Our data reveal an intriguing relationship between deep‐sourced fluid flow upward along a tectonic fault, and shallower flow through dipping sediments. This complex fluid flow pattern has led to disruption of the gas hydrate system and the formation of concentrated gas hydrate deposits within the dipping sediments. Our study highlights the relationships between relatively deep tectonic processes (faulting and fluid flow) and the shallow process of gas hydrate formation in an active subduction zone. Key Points A distinct gas‐hydrate to free‐gas transition is mapped using high‐ and low‐frequency seismic data Gas and hydrate accumulations in the Urutī Basin are controlled by the structural setting, ongoing deep‐sourced fluid flow, and near surface stratigraphy Regions of high modeled heat flow can be directly related to accumulations of gas and gas hydrates

Description: We conducted two‐dimensional numerical simulations to investigate the mechanisms underlying the strong spatiotemporal correlation observed between submarine landslides and gas hydrate dissociation due to glacial sea‐level drops. Our results suggest that potential plastic deformation or slip could occur at localized and small scales in the shallow‐water portion of the gas hydrate stability zone (GHSZ). This shallow‐water portion of the GHSZ typically lies within the area enclosed by three points: the BGHSZ–seafloor intersection, the seafloor at ∼600 m below sea level (mbsl), and the base of the GHSZ (BGHSZ) at ∼1,050 mbsl in low‐latitude regions. The deep BGHSZ (〉1,050 mbsl) could not slip; therefore, the entire BGHSZ was not a complete slip surface. Glacial hydrate dissociation alone is unlikely to cause large‐scale submarine landslides. Observed deep‐water (much greater than 600 mbsl) turbidites containing geochemical evidence of glacial hydrate dissociation potentially formed from erosion or detachment in the GHSZ pinch‐out zone. Plain Language Summary Many submarine landslides spatiotemporally correlate with gas hydrate dissociation. However, direct mechanical evidence supporting whether the overpressure and deformation due to glacial sea‐level drop‐induced hydrate dissociation are adequate for triggering submarine landslides is lacking. Here, we present two‐dimensional thermal‐hydraulic‐chemical and geomechanical models of a gas‐hydrate system in response to glacial sea‐level drops and conduct sensitivity analyses of the model behavior under a wide range of key conditions from a global perspective. Our simulations suggest that glacial hydrate dissociation might induce plastic deformation or slip at localized and small scales only possibly within the shallow‐water portion of the hydrate stability zone. The deep part (〉1,050 m below sea level) of the bottom boundary of the hydrate stability zone could not slip; therefore, the entire bottom boundary of the hydrate stability zone was not a complete slip surface. We demonstrate that glacial hydrate dissociation alone is unlikely to trigger large‐scale submarine landslides. Our work highlights the vicinity of the upper limit of the hydrate stability zone (where the base of the hydrate stability zone intersects the seafloor) as an important area for investigating overpressure and focused fluid flow, localized plastic deformation or slip, and downslope sediment transport related to glacial hydrate dissociation. Key Points Glacial hydrate dissociation might cause potential plastic deformation or slip at localized and small scales in shallow parts of the GHSZ The large deformation surface at the BGHSZ boundary of the potential plastic deformation zone was not a complete slip surface Glacial sea‐level drop‐induced gas hydrate dissociation alone is unlikely to have caused large‐scale submarine landslides

Description: Carbon disulfide (CS2) has recently gained attention as an important precursor for the atmospheric trace gas carbonyl sulfide (OCS), which delivers sulfur to the stratospheric sulfur layer and impacts the radiative budget of the Earth. CS2 is naturally produced in the ocean and emitted to the atmosphere. However, the magnitude of its marine emissions is only poorly constrained due to lacking understanding of its production and consumption processes. Here, we present incubation experiments with and without UV light treatment and provide evidence for a previously not considered UV-light-driven degradation process of CS2 in seawater, following first-order kinetics. In addition to its already known photochemical production process, CS2 production is found in the dark, depending on the amount of dissolved organic sulfur present in seawater. We provide novel production and consumption rates of CS2 in seawater that pave the way toward mechanistically quantifying marine emissions of this important trace gas. Key Points: - Carbon disulfide in seawater is degraded by UV light at time scales of days - Carbon disulfide is produced in seawater without UV light at rates comparable to photochemical production - Carbon disulfide dark production is limited by dissolved organic sulfur

Description: Current earthquake forecasting approaches are mainly based on probabilistic assumptions, as earthquakes seem to occur randomly. Such apparent randomness can however be caused by deterministic chaos, rendering deterministic short‐term forecasts possible. Due to the short historical and instrumental record of earthquakes, chaos detection has proven challenging, but more frequently occurring slow slip events (SSE) are promising candidates to probe for determinism. Here, we characterize the SSE signatures obtained from GNSS position time series in the Hikurangi Subduction Zone (New Zealand) to investigate whether the seemingly random SSE occurrence is governed by chaotic determinism. We find evidence for deterministic chaos for stations recording shallow SSEs, suggesting that short‐term deterministic forecasting of SSEs, similar to weather forecasts, might indeed be possible over timescales of a few weeks. We anticipate that our findings could open the door for next‐generation SSE forecasting, adding new tools to existing probabilistic approaches. Plain Language Summary Since earthquakes appear to occur randomly, the currently available probabilistic predictions are based on past earthquake records. These predictions estimate the likelihood of an earthquake of a given magnitude occurring within a defined time period. In contrast to such probabilistic approaches, deterministic systems are fully predictable, albeit often confined to short time scales due to their potential chaotic behavior. Probing for deterministic predictability in the earthquake cycle is intractable due to the limited historical instrumental record. However, frequently occurring slow slip events ‐ captured by transient GNSS displacements that can last several weeks ‐ provide a unique opportunity to explore deterministic predictability in these types of slow earthquakes. By studying GNSS time series from various stations on New Zealand’s North Island, we have discovered evidence suggesting that these irregularly occurring slow slip events might be governed by chaotic determinism. This implies the potential to forecast both timing and magnitude of slow slip events a few weeks in advance using deterministic methods, much like we predict weather patterns. Consequently, our theoretical findings could therefore pave the way for innovative approaches to short‐term slow slip forecasting. Key Points Nonlinear analysis of GNSS displacement time series unveils evidence for deterministic chaos in slow slip events in New Zealand Our theoretical findings imply that irregularly occurring slow slip events could potentially be forecasted a few weeks in advance

Description: The Sr/Ca ratio of coral aragonite is used to reconstruct past sea surface temperature (SST). Twentyone laboratories took part in an interlaboratory study of coral Sr/Ca measurements. Results show interlaboratory bias can be significant, and in the extreme case could result in a range in SST estimates of 7°C. However, most of the data fall within a narrower range and the Porites coral reference material JCp- 1 is now characterized well enough to have a certified Sr/Ca value of 8.838 mmol/mol with an expanded uncertainty of 0.089 mmol/mol following International Association of Geoanalysts (IAG) guidelines. This uncertainty, at the 95% confidence level, equates to 1.5°C for SST estimates using Porites, so is approaching fitness for purpose. The comparable median within laboratory error is 〈0.5°C. This difference in uncertainties illustrates the interlaboratory bias component that should be reduced through the use of reference materials like the JCp-1. There are many potential sources contributing to biases in comparative methods but traces of Sr in Ca standards and uncertainties in reference solution composition can account for half of the combined uncertainty. Consensus values that fulfil the requirements to be certified values were also obtained for Mg/Ca in JCp-1 and for Sr/Ca and Mg/Ca ratios in the JCt-1 giant clam reference material. Reference values with variable fitness for purpose have also been obtained for Li/Ca, B/Ca, Ba/Ca, and U/Ca in both reference materials. In future, studies reporting coral element/Ca data should also report the average value obtained for a reference material such as the JCp-1.

Description: Carbonate buildups and mounds are impressive biogenic structures throughout Earth history. In the recent NE Atlantic, cold-water coral (CWC) reefs form giant carbonate mounds of up to 300 m of elevation. The expansion of these coral carbonate mounds is paced by climatic changes during the past 2.7 Myr. Environmental control on their development is directly linked to controls on its main constructors, the reef-building CWCs. Seawater density has been identified as one of the main controlling parameter of CWC growth in the NE Atlantic. One possibility is the formation of a pycnocline above the carbonate mounds, which is increasing the hydrodynamic regime, supporting elevated food supply, and possibly facilitating the distribution of coral larvae. The potential to reconstruct past seawater densities from stable oxygen isotopes of benthic foraminifera has been further developed: a regional equation gives reliable results for three different settings, peak interglacials (e.g., Holocene), peak glacials (e.g., Last Glacial Maximum), and intermediate setting (between the two extremes). Seawater densities are reconstructed for two different NE Atlantic CWC carbonate mounds in the Porcupine Seabight indicating that the development of carbonate mounds is predominantly found at a seawater density range between 27.3 and 27.7 kg m−3 (σΘ notation). Comparable to recent conditions, we interpret the reconstructed density range as a pycnocline serving as boundary layer, on which currents develop, carrying nutrition and possibly coral larvae. The close correlation of CWC reef growth with reconstructed seawater densities through the Pleistocene highlights the importance of pycnoclines and intermediate water mass dynamics.

Description: The complex deglacial to Holocene oceanographic development in the Gulf of Guayaquil (Eastern Equatorial Pacific) is reconstructed for sea surface and subsurface ocean levels from (isotope) geochemical proxies based on marine sediment cores. At sea surface, southern sourced Cold Coastal Water and tropical Equatorial Surface Water/Tropical Surface Water are intimately related. In particular since ~10 ka, independent sea surface temperature proxies capturing different seasons emphasize the growing seasonal contrast in the Gulf of Guayaquil, which is in contrast to ocean areas further offshore. Cold Coastal Water became rapidly present in the Gulf of Guayaquil during the austral winter season in line with the strengthening of the Southeast Trades, while coastal upwelling off Peru gradually intensified and expanded northward in response to a seasonally changing atmospheric circulation pattern affecting the core locations intensively since 4 ka BP. Equatorial Surface Water, instead, was displaced and Tropical Surface Water moved northward together with the Equatorial Front. At subsurface, the presence of Equatorial Under Current-sourced Equatorial Subsurface Water was continuously growing, prominently since ~10–8 ka B.P. During Heinrich Stadial 1 and large parts of the Bølling/Allerød, and similarly during short Holocene time intervals at ~5.1–4 ka B.P. and ~1.5–0.5 ka B.P., the admixture of Equatorial Subsurface Water was reduced in response to both short-term weakening of Equatorial Under Current strength from the northwest and emplacement by tropical Equatorial Surface Water, considerably warming the uppermost ocean layers.

Description: High-latitude cold-water coral (CWC) reefs are particularly susceptible due to enhanced CO2 uptake in these regions. Using precisely dated (U/Th) CWCs (Lophelia pertusa) retrieved during research cruise POS 391 (Lopphavet 70.6°N, Oslofjord 59°N) we applied boron isotopes (δ11B), Ba/Ca, Li/Mg and U/Ca ratios to reconstruct the environmental boundary conditions of CWC reef growth. The sedimentary record from these CWC reefs reveals a lack of corals between ∼ 6.4 and 4.8 ka. The question remains if this phenomenon is related to changes in the carbonate system or other causes. The initial postglacial setting had elevated Ba/Ca ratios, indicative of meltwater fluxes showing a decreasing trend towards cessation at 6.4 ka with a oscillation pattern similar to continental glacier fluctuations. Downcore U/Ca ratios reveal an increasing trend, which is outside the range of modern U/Ca variability in L. pertusa, suggesting changes of seawater pH near 6.4 ka. The reconstructed BWT at Lopphavet reveals a striking similarity to Barent Sea-Surface and sub-Sea-Surface-Temperature records. We infer that meltwater pulses weakened the North Atlantic Current system resulting in southward advances of cold and CO2 rich Arctic waters. A corresponding shift in the δ11B record from ∼ 25.0‰ to ∼ 27.0 ‰ probably implies enhanced pH-up regulation of the CWCs due to the higher pCO2 concentrations of ambient seawater, which hastened Mid-Holocene CWC reef decline on the Norwegian Margin.

Description: “Flip‐flop” detachment mode represents an endmember type of lithosphere‐scale faulting observed at almost amagmatic sections of ultraslow‐spreading mid‐ocean ridges. Recent numerical experiments using an imposed steady temperature structure show that an axial temperature maximum is essential to trigger flip‐flop faults by focusing flexural strain in the footwall of the active fault. However, ridge segments without significant melt budget are more likely to be in a transient thermal state controlled, at least partly, by the faulting dynamics themselves. Therefore, we investigate which processes control the thermal structure of the lithosphere and how feedbacks with the deformation mechanisms can explain observed faulting patterns. We present results of 2‐D thermo‐mechanical numerical modeling including serpentinization reactions and dynamic grain size evolution. The model features a novel form of parametrized hydrothermal cooling along fault zones as well as the thermal and rheological effects of periodic sill intrusions. We find that the interplay of hydrothermal fault zone cooling and periodic sill intrusions in the footwall facilitates the flip‐flop detachment mode. Hydrothermal cooling of the fault zone pushes the temperature maximum into the footwall, while intrusions near the temperature maximum further weaken the rock and promote the formation of new faults with opposite polarity. Our model allows us to put constraints on the magnitude of two processes, and we obtain most reasonable melt budgets and hydrothermal heat fluxes if both are considered. Furthermore, we frequently observe two other faulting modes in our experiments complementing flip‐flop faulting to yield a potentially more robust alternative interpretation for existing observations. Plain Language Summary At mid‐ocean ridges, two plates diverge and new seafloor is created. The nature and appearance of this new seafloor strongly depend on spreading velocity and the availability of magmatic melts. At one of the melt‐poorest and slowest‐spreading ridges, a special form of large‐scale tectonic faults, so‐called flip‐flop detachments, can be observed. Tectonic faults can act as pathways for fluids circulating through the seafloor, which provides a significant cooling effect for the young plate. The interplay of magmatic activity, faulting and fluid circulation is evident at many different ridges with different magmatic activity and spreading rates. Flip‐flop faulting is restricted to only a few ridge sections worldwide, and we here investigate the prerequisites for this special spreading mode. To do so, we set up a computer model of an ultraslow‐spreading mid‐ocean ridge including the effects of sparse magmatism as well as the cooling effect associated with fluid circulation. We find that feedbacks between faulting dynamics, hydrothermal cooling and magmatic activity control the magnitude and spatial location of each individual process. Seafloor and subsurface observations are best explained by calculations with moderate melt input and hydrothermal circulation acting together. Key Points We implemented hydrothermal cooling and magmatic intrusion in a thermo‐mechanical model to explain detachment faulting at ultraslow ridges Stable flip‐flop detachment faulting is observed for setups considering both melt input and hydrothermal heat fluxes at realistic magnitudes Two other faulting modes frequently observed in our model offer potential alternative interpretations for existing seafloor observations

Description: Submarine landslides pose a hazard to coastal communities and critical seafloor infrastructure, occurring on all of the world's continental margins, from coastal zones to hadal trenches. Offshore monitoring has been limited by the largely unpredictable occurrence of submarine landslides and the need to cover large regions. Recent subsea monitoring has provided new insights into the preconditioning and run-out of submarine landslides using active geophysical techniques. However, these tools measure a small spatial footprint and are power- and memory-intensive, thus limiting long-duration monitoring. Most landslide events remain unrecorded. In this chapter, we first show how passive acoustic and seismologic techniques can record acoustic emissions and ground motions created by terrestrial landslides. This terrestrial-focused research has catalyzed advances in characterizing submarine landslides using onshore and offshore networks of broadband seismometers, hydrophones, and geophones. We discuss new insights into submarine landslide preconditioning, timing, location, velocity, and down-slope evolution arising from these advances. Finally, we outline challenges, emphasizing the need to calibrate seismic and acoustic signals generated by submarine landslides. Passive seismic and acoustic sensing has a strong potential to enable more complete hazard catalogs to be built and open the door to emerging techniques (such as fiber-optic sensing) to fill key knowledge gaps.

Description: Key Points: - North Atlantic biases are alleviated by an eddying nested ocean configuration embedded in a global climate model, FOCI-VIKING10 - It is indicated that reduction of the North Atlantic biases could improve the representation of NAO sub-decadal (8 years) variability - For detecting weak external imprints with limited computational resources, an ensemble with a coarse-resolution model is favorable Increasing the horizontal resolution of an ocean model is frequently seen as a way to reduce the model biases in the North Atlantic, but we are often limited by computational resources. Here, a two-way nested ocean model configuration (VIKING10) that consists of a high-resolution (1/10°) component and covers the northern North Atlantic, is embedded in a 1/2° ocean grid as part of the global chemistry-climate model, FOCI (called FOCI-VIKING10). This configuration yields a significantly improved path of the North Atlantic current (NAC), which here reduces the North Atlantic cold bias by ∼50%. Compared with the coarse-resolution, non-eddying model, the improved thermal state of upper ocean layers and surface heat fluxes in a historical simulation based on FOCI-VIKING10 are beneficial for simulating the subdecadal North Atlantic Oscillation (NAO) variability (i.e., a period of 8 years). A northward drift of the NAO-forced ocean thermal anomalies as seen in observations and the eddying FOCI-VIKING10, provide a lagged ocean feedback to the NAO via changes in the net surface heat flux, leading to the NAO periodicity of 8 years. This lagged feedback and the 8 years variability of the NAO cannot be captured by the non-eddying standard FOCI historical simulation. Furthermore, the argumentative responses of the North Atlantic to the 11-year solar cycle are re-examined in this study. The reported solar-induced NAO-like responses are confirmed in the 9-member ensemble mean based on FOCI but with low robustness among individual members. A lagged NAO-like response is only found in the nested eddying simulation but absent from the non-eddying reference simulation, suggesting North Atlantic biases importantly limit climate model capability to realistically solar imprints in North Atlantic climate.

Description: The past ∼200 million years of Earth's geomagnetic field behavior have been recorded within oceanic basalts, many of which are only accessible via scientific ocean drilling. Obtaining the best possible paleomagnetic measurements from such valuable samples requires an a priori understanding of their magnetic mineralogies when choosing the most appropriate protocol for stepwise demagnetization experiments (either alternating field or thermal). Here, we present a quick, and non‐destructive method that utilizes the amplitude‐dependence of magnetic susceptibility to screen submarine basalts prior to choosing a demagnetization protocol, whenever conducting a pilot study or other detailed rock‐magnetic characterization is not possible. We demonstrate this method using samples acquired during International Ocean Discovery Program Expedition 391. Our approach is rooted in the observation that amplitude‐dependent magnetic susceptibility is observed in basalt samples whose dominant magnetic carrier is multidomain titanomagnetite (∼TM 60–65 , (Ti 0.60–0.65 Fe 0.35–0.40 )Fe 2 O 4 ). Samples with low Ti contents within titanomagnetite or samples that have experienced a high degree of oxidative weathering do not display appreciable amplitude dependence. Due to their low Curie temperatures, basalts that possess amplitude‐dependence should ideally be demagnetized either using alternating fields or via finely‐spaced thermal demagnetization heating steps below 300°C. Our screening method can enhance the success rate of paleomagnetic studies of oceanic basalt samples. Plain Language Summary Oceanic basalts are ideal recorders of the Earth's magnetic field. To decipher magnetic histories recorded in rocks, paleomagnetists need to isolate the magnetization directions and intensities within rocks by one of two possible methods. One method typically involves progressively heating the samples to high temperatures. The other method involves exposing samples to alternating magnetic fields with increasing peak field intensities. Both of these methods are ultimately destructive to the original magnetization preserved within rocks. However, without knowledge of a given rock's magnetic mineralogy, randomly choosing thermal or alternating field demagnetization methods may result in high failure rates. We developed a pre‐screening method to help decide which cleaning method will likely be more successful for a given sample based on low‐field magnetic susceptibility measurements. These measurements do not affect the original magnetic information recorded in a rock, thereby permitting subsequent paleomagnetic studies on the same sample. Our technique can be performed as rapidly as 2 min per sample, is non‐destructive, and does not require complicated sample preparation. Key Points Paleomagnetic studies utilize either alternating field or thermal demagnetization, but it is difficult to choose the best protocol a priori Amplitude‐dependence of magnetic susceptibility measurements permits preliminary magnetic mineralogy characterization in submarine basalts Rapid amplitude‐dependence measurements may aid in deciding upon the best demagnetization protocol for submarine basalt samples

Description: Marginal seas influenced by large rivers are characterized by complex hydrodynamic and organic matter cycling processes. However, the impacts of hydrodynamics on the composition and reactivity of particulate organic matter (POM) remain unclear. Here we conducted a comprehensive study on the bulk, molecular and biological properties of suspended POM in the Changjiang Estuary and adjacent area subjected to strong currents, eddies as well as typhoons during spring and autumn. D/L‐enantiomers of particulate amino acids (PAA) were analyzed to evaluate the bioreactivity of POM and quantify bacterial‐derived organic carbon. We found that POM bioavailability as indicated by carbon‐normalized yields of PAA (PAA‐C%) reflected the ecosystem productivity. Relatively high PAA‐C% values (20−35%) were observed in productive areas influenced by Changjiang River plume, cyclonic eddies and typhoons, likely related to the enhanced nutrient availability arising from hydrodynamic processes. In contrast, the oligotrophic Taiwan Warm Current‐influenced regions featured relatively low POM bioavailability (PAA‐C% 〈 10%) despite typhoons facilitating water mixing. The PAA‐C% values showed a significant positive correlation with extracellular enzyme activity, indicating that bioavailable POM can rapidly stimulate heterotrophic transformation. Hot spots of elevated bioavailable POM showed high contributions of bacterial organic carbon. A large portion (∼2/3) of bacterial organic carbon was present in the form of bacterial detritus, suggesting that patches of these biological hot spots represent important sites of carbon sequestration. Together, our findings indicate that fresh POM production is largely controlled by nutrient supply driven by hydrodynamic processes, with important implications for carbon sequestration in the dynamic ocean margins. Plain Language Summary Marginal seas are subject to complex hydrodynamic processes and play an important role in carbon sequestration. Disentangling the linkages between hydrodynamics and organic carbon reactivity and composition is crucial to understanding the regional carbon cycle. Here we collected suspended particulate organic matter (POM) in the Changjiang Estuary and adjacent coastal areas. Based on the biomarker D/L‐amino acids, we assessed the bioavailability of POM and quantified the organic carbon originating from bacteria. We found that high bioactivity of POM occurred in productive Changjiang River plume, cyclonic eddy, and typhoon influenced areas. These hydrodynamic processes appear to increase nutrient availability, therefore promoting phytoplankton growth. Bioavailable POM can rapidly stimulate heterotrophic activity and facilitate the transformation of algal‐derived organic carbon to bacterial detritus, thus contributing to carbon sequestration. Our findings suggest that the production of bioavailable POM is largely controlled by hydrodynamically driven nutrient supply. Key Points We use D/L‐amino acids to assess the bioreactivity and bacterial origins of particulate organic matter (POM) in the dynamic Changjiang Estuary and adjacent area High bioavailability of POM occurs in productive regions affected by Changjiang River plume, cyclonic eddies and typhoons Hot spots of bioavailable POM represent important sites for carbon sequestration

Description: The Congo River supplies vast quantities of trace metals (TMs) to the South Atlantic Ocean, but TM budgets for the Congo plume derived using radium isotopes for GEOTRACES cruise GA08 suggest additional input other than the river outflow. Considering the tight correlations between most dissolved TMs and salinity in the plume and the high rainfall during the wet season over the Congo shelf, we hypothesized that wet atmospheric deposition is a TM source to the Congo plume. Observed TM concentrations in rainwaters across the Congo shelf were mostly comparable to values from previous work in the North Atlantic and Mediterranean Sea. Wet deposition contributed the equivalent of 43% dCd, 21% dCu, 20% dPb and 68% dZn of the Congo River fluxes. Our findings show an important role of wet deposition in supplying TMs to the South Atlantic overlapping with the region that receives substantial TM fluxes from the Congo River. Key Points The Congo River is an important source of trace metals (TMs) to the South Atlantic Ocean revealed by data from GEOTRACES cruise GA08 Wet deposition (rainfall) is identified as an additional TM source to the Congo plume by concurrently considering river and rain data Rainfall supplies anthropogenic dTMs (Cd, Cu, Pb and Zn) with fluxes equivalent to 20%–68% of those from the Congo River on the Congo shelf Plain Language Summary The Congo River has the second largest freshwater discharge volume globally and creates an extensive near-equatorial plume into the Atlantic Ocean. The Congo plume constitutes an important source of trace metals (TMs) to the ocean, which impacts biogeochemical cycles in the tropical and subtropical ocean. However, existing work suggests a discrepancy within the TM budgets in the Congo plume and points to unknown source other than the Congo River or shelf sediments. Most TM concentrations across the Congo plume remain tightly correlated with salinity, suggesting that any additional sources are likely also freshwater-derived or enter the ocean at the river mouth coincidently with direct riverine TM inputs. Here, TM concentrations in ocean, river and rainwater collected during the GEOTRACES GA08 cruise are combined to suggest that wet deposition augmented some Congo TM fluxes to the ocean. Fluxes of anthropogenic Cd, Cu, Pb and Zn to the Congo shelf from wet deposition are of the same order of magnitude as the Congo River. Concentrations of these elements in rainwater are similar to prior observations reported for the North Atlantic and Mediterranean Sea, suggesting that a large fraction of the global range of rainwater concentrations over the ocean has been captured in our observations.

Description: Iodine cycling in the ocean is closely linked to productivity, organic carbon export, and oxygenation. However, iodine sources and sinks at the seafloor are poorly constrained, which limits the applicability of iodine as a biogeochemical tracer. We present pore water and solid phase iodine data for sediment cores from the Peruvian continental margin, which cover a range of bottom water oxygen concentrations, organic carbon rain rates and sedimentation rates. By applying a numerical reaction‐transport model, we evaluate how these parameters determine benthic iodine fluxes and sedimentary iodine‐to‐organic carbon ratios (I:C org ) in the paleo‐record. Iodine is delivered to the sediment with organic material and released into the pore water as iodide (I − ) during early diagenesis. Under anoxic conditions in the bottom water, most of the iodine delivered is recycled, which can explain the presence of excess dissolved iodine in near‐shore anoxic seawater. According to our model, the benthic I − efflux in anoxic areas is mainly determined by the organic carbon rain rate. Under oxic conditions, pore water dissolved I − is oxidized and precipitated at the sediment surface. Much of the precipitated iodine re‐dissolves during early diagenesis and only a fraction is buried. Particulate iodine burial efficiency and I:C org burial ratios do increase with bottom water oxygen. However, multiple combinations of bottom water oxygen, organic carbon rain rate and sedimentation rate can lead to identical I:C org , which limits the utility of I:C org as a quantitative oxygenation proxy. Our findings may help to better constrain the ocean's iodine mass balance, both today and in the geological past. Key Points The impact of early diagenesis on benthic iodine fluxes and iodine burial was quantitatively evaluated using a reaction‐transport model Dissolved iodine anomalies in the water column are indicative of benthic efflux from anoxic sediments with high organic carbon turnover Not only bottom water oxygen but also organic carbon delivery and sedimentation rate determine sedimentary iodine‐to‐organic carbon ratios

Description: The understanding of silicate weathering and its role as a sink for atmospheric CO 2 is important to get a better insight into how the Earth shifts from warm to cool climates. The lithium isotope composition (δ 7 Li) of marine carbonates can be used as a proxy to track the past chemical weathering of silicates. A high‐resolution δ 7 Li record would be helpful to evaluate the role of silicate weathering during the late Cretaceous climate cooling. Here, we assess chalk as a potential archive for reconstructing Late Cretaceous seawater Li isotope composition by comparing Maastrichtian chalk from Northern Germany (Hemmoor, Kronsmoor) to a Quaternary coccolith ooze from the Manihiki Plateau (Pacific Ocean) as a lithological analog to modern conditions. We observe a negative offset of 3.9 ± 0.6‰ for the coccolith ooze relative to the modern seawater Li isotope composition (+31.1 ± 0.3‰; 2SE; n = 54), a value that falls in the range of published offsets for modern core‐top samples and for brachiopod calcite. Further, the negative offset between the Li isotope compositions of Manihiki coccolith ooze and modern planktonic foraminifera is 2.3 ± 0.6‰. Although chalk represents a diagenetically altered modification of pelagic nannofossil ooze, manifested by changes in the composition of trace elements, we observe a consistent offset of Li isotope data between Maastrichtian chalk and Maastrichtian planktonic foraminiferal data (−1.4 ± 0. 5‰) that lies within the uncertainty of modern values. We therefore suggest that chalk can be used as a reliable archive for δ 7 Li reconstructions. Key Points Chalk is a reliable archive for the Li isotope composition of seawater Coccolith ooze has a negative offset of 3.9 ± 0.6‰ from modern seawater for Li isotope ratios The estimated mean value for the late Maastrichtian seawater Li isotope composition is +27.5 ± 1.0‰

Description: The coastal ocean contributes to regulating atmospheric greenhouse gas concentrations by taking up carbon dioxide (CO2) and releasing nitrous oxide (N2O) and methane (CH4). In this second phase of the Regional Carbon Cycle Assessment and Processes (RECCAP2), we quantify global coastal ocean fluxes of CO2, N2O and CH4 using an ensemble of global gap-filled observation-based products and ocean biogeochemical models. The global coastal ocean is a net sink of CO2 in both observational products and models, but the magnitude of the median net global coastal uptake is similar to 60% larger in models (-0.72 vs. -0.44 PgC year-1, 1998-2018, coastal ocean extending to 300 km offshore or 1,000 m isobath with area of 77 million km2). We attribute most of this model-product difference to the seasonality in sea surface CO2 partial pressure at mid- and high-latitudes, where models simulate stronger winter CO2 uptake. The coastal ocean CO2 sink has increased in the past decades but the available time-resolving observation-based products and models show large discrepancies in the magnitude of this increase. The global coastal ocean is a major source of N2O (+0.70 PgCO2-e year-1 in observational product and +0.54 PgCO2-e year-1 in model median) and CH4 (+0.21 PgCO2-e year-1 in observational product), which offsets a substantial proportion of the coastal CO2 uptake in the net radiative balance (30%-60% in CO2-equivalents), highlighting the importance of considering the three greenhouse gases when examining the influence of the coastal ocean on climate. The coastal ocean regulates greenhouse gases. It acts as a sink of carbon dioxide (CO2) but also releases nitrous oxide (N2O) and methane (CH4) into the atmosphere. This synthesis contributes to the second phase of the Regional Carbon Cycle Assessment and Processes (RECCAP2) and provides a comprehensive view of the coastal air-sea fluxes of these three greenhouse gases at the global scale. We use a multi-faceted approach combining gap-filled observation-based products and ocean biogeochemical models. We show that the global coastal ocean is a net sink of CO2 in both observational products and models, but the coastal uptake of CO2 is similar to 60% larger in models than in observation-based products due to model-product differences in seasonality. The coastal CO2 sink is strengthening but the magnitude of this strengthening is poorly constrained. We also find that the coastal emissions of N2O and CH4 counteract a substantial part of the effect of coastal CO2 uptake in the atmospheric radiative balance (by 30%-60% in CO2-equivalents), highlighting the need to consider these three gases together to understand the influence of the coastal ocean on climate. We synthesize air-sea fluxes of CO2, nitrous oxide and methane in the global coastal ocean using observation-based products and ocean models The coastal ocean CO2 sink is 60% larger in ocean models than in observation-based products due to systematic differences in seasonality Coastal nitrous oxide and methane emissions offset 30%-60% of the CO2 coastal uptake in the net radiative balance

Description: Accessible seafloor minerals located near mid‐ocean ridges are noticed to mitigate the projected metal demands of the net‐zero energy transition, promoting growing interest in quantifying the global distributions of seafloor massive sulfides (SMS). Mineral potentials are commonly estimated using geophysical and geological data that lastly rely on additional confirmation studies using sparsely available, locally limited, seafloor imagery, grab samples, and coring data. This raises the challenge of linking in situ confirmation data to geophysical data acquired at disparate spatial scales to obtain quantitative mineral predictions. Although multivariate data sets for marine mineral research are incessantly acquired, robust, integrative data analysis requires cumbersome workflows and experienced interpreters. We introduce an automated two‐step machine learning approach that integrates the mound detection through image segmentation with geophysical data. SMS predictors are subsequently clustered into distinct classes to infer marine mineral potentials that help guide future exploration. The automated workflow employs a U‐Net convolutional neural network to identify mound structures in bathymetry data and distinguishes different mound classes through the classification of mound architectures and magnetic signatures. Finally, controlled source electromagnetic data are utilized together with in situ sampling data to reassess predictions of potential SMS volumes. Our study focuses on the Trans‐Atlantic Geotraverse area, which is among the most explored SMS areas worldwide and includes 15 known SMS sites. The automated workflow classifies 14 of the 15 known mounds as exploration targets of either high or medium priority. This reduces the exploration area to less than 7% of the original survey area from 49 to 3.1 km 2 . Key Points A two‐step machine learning workflow identifies mound structures in bathymetry data and classifies their origins based on auxiliary data Significant increase in potential seafloor massive sulfides (SMS) edifices detected within the trans‐Atlantic geo‐traverse hydrothermal field distributed within latitudinal bands SMS mineral potential is likely lower than previously assumed due to heterogeneously distributed mineralization within mounds

Description: Burial driven recycling is an important process in the natural gas hydrate (GH) systems worldwide, characterized by complex multiphysics interactions like gas migration through an evolving gas hydrate stability zone (GHSZ), competing gas-water-hydrate (i.e., fluid-fluid-solid) phase transitions, locally appearing and disappearing phases, and evolving sediment properties (e.g., permeability, reaction surface area, and capillary entry pressure). Such a recycling process is typically studied in homogeneous or layered sediments. However, there is mounting evidence that structural heterogeneity and anisotropy linked to normal and inclined fault systems or anomalous sediment layers have a strong impact on the GH dynamics. Here, we consider the impacts of such a structurally complex media on the recycling process. To capture the properties of the anomalous layers accurately, we introduce a fully mass conservative, high-order, discontinuous Galerkin (DG) finite element based numerical scheme. Moreover, to handle the rapidly switching thermodynamic phase states robustly, we cast the problem of phase transitions as a set of variational inequalities, and combine our DG discretization scheme with a semi-smooth Newton solver. Here, we present our new simulator, and demonstrate using synthetic geological scenarios, (a) how the presence of an anomalous high-permeability layer, like a fracture or brecciated sediment, can alter the recycling process through flow-localization, and more importantly, (b) how an incorrect or incomplete approximation of the properties of such a layer can lead to large errors in the overall prediction of the recycling process. Key Points Structural heterogeneity linked to inclined fault systems or anomalous sediment layers have a strong impact on the gas hydrate dynamics The presence of anomalous high-permeability layers within gas hydrate stability zone alters the recycling process through flow-localization The presented discontinuous Galerkin scheme is able to accurately capture the gas hydrate recycling processes through strongly anisotropic materials

Description: The Nazca Ridge (NR) was formed near the interaction of a hotspot mantle plume and an active spreading center. We use active-source wide-angle seismic data to obtain 2-D Vp and Vs tomographic models, and hence the Poisson's ratio (ν) structure beneath the NR. Results show a ∼2 km thick seismic layer 2A with ν values of 0.25–0.32 in the uppermost crust interpreted as pillow basalts with a low degree of fracturing and/or hydrothermal alteration. The 2A/B boundary layer presents ν values of 0.27–0.29 consistent with pillow basalts/sheeted dykes units. A ∼3 km layer 2B overlies a ∼10 km layer 3 with ν values of 0.24–0.3 at the 2/3 boundary layer. The lowermost layer 3 presents ν values of 0.28 ± 0.02 suggesting an increase in Mg content (≥10% wt). The NR crust (∼15 km thick) requires an increment of the asthenospheric mantle potential temperature in ∼100°C formed by passive adiabatic decompression melting.

Description: Upwelling of subsurface waters injects macronutrients (fixed N, P and Si) and micronutrient trace metals (TMs) into surface waters supporting elevated primary production in Eastern Boundary Upwelling Regions (EBUR). The eastern South Atlantic features a highly productive shelf sea transitioning to a low productivity N-Fe (co)limited open ocean. Whilst a gradient in most TM concentrations is expected in any off-shelf transect, the factors controlling the magnitude of cross-shelf TM fluxes are poorly constrained. Here, we present dissolved TM concentrations of Fe, Co, Mn, Cd, Ni and Cu within the Benguela Upwelling System (BUS) from the coastal section of the GEOTRACES GA08 cruise. Elevated dissolved Fe, Co, Mn, Cd, Ni, Cu and macronutrient concentrations were observed near shelf sediments. Benthic sources supplied 2.22 ± 0.99 μmol Fe m-2 d-1, 0.05 ± 0.03 μmol Co m-2 d-1, 0.28 ± 0.11 μmol Mn m-2 d-1 and were found to be the dominant source to shallow shelf waters compared to atmospheric depositions. Similarly, off-shelf transfer was a more important source of TMs to the eastern South Atlantic Ocean compared to atmospheric deposition. Assessment of surface (shelf, upper 200 m) and subsurface (shelf edge, 200 - 500 m) fluxes of Fe and Co indicated TM fluxes from subsurface were 2 - 5 times larger than those from surface into the eastern South Atlantic Ocean. Under future conditions of increasing ocean deoxygenation, these fluxes may increase further, potentially contributing to a shift towards more extensive regional limitation of primary production by fixed N availability. Key Points Shelf sediments release redox-sensitive trace metals (TMs) to overlying oxygen-depleted waters in the Benguela Upwelling System (BUS) Sediment-derived TMs are upwelled and laterally transported constituting a major source to shelf waters and to the eastern South Atlantic Subsurface fluxes of dissolved Fe and Co from the shelf edge play an important role in supplying Fe and Co to the eastern South Atlantic

Description: The Baltic Basin is known for its numerous Paleozoic hydrocarbon reservoirs. There is published evidence that hydrocarbons are leaking from the seafloor, however, little is known about the hydrocarbon migration pathways from Paleozoic source and reservoir rocks toward the seafloor and their escape structures. To investigate these processes, we utilize a new set of multibeam, parametric sediment sub-bottom profiler and 2D seismic reflection data. The integrated analysis of seismic profiles, diffraction imaging and bathymetric maps allow to identify a hydrocarbon migration system within Silurian and Devonian strata that consists of layer parallel and updip migration beneath sealing layers, migration across seals along faults, and seafloor escape structures in form of elongated depressions. The general migration trend is directed updip, from the Paleozoic reservoirs below the southeastern Baltic Sea toward the Gotland Depression in the northwest. The locations of the hydrocarbon escape structures at the seafloor and their elongated shape are mainly controlled by the regional geological setting of outcropping Paleozoic layers. In addition, iceberg scouring may have facilitated hydrocarbon migration through the Quaternary deposits. The description of this hydrocarbon migration system fills the gap between the known reservoirs and the observed hydrocarbon accumulations and seepages. With regard to potential Carbon Capture and Storage projects, the identification of this hydrocarbon migration system is of great importance, as potential storage sites may be leaking. Key Points Numerous elongated fluid escape depressions are observed at the eastern margin of the Gotland Deep, central Baltic Sea First evidence for fluid migration pathways from Paleozoic toward Quaternary strata in the region Locations of fluid escape is controlled by the regional tectonic setting

Description: Sub-seabed fluid flow, gas hydrate accumulation and seafloor methane seepage are tightly interwoven processes with implications for marine biodiversity, ocean chemistry and seafloor stability. We combine long-offset seismic reflection data with high-resolution seismic data to investigate shallow structural deformation and its relationship to focused gas migration and hydrate accumulation in the southern Hikurangi subduction wedge. Anticlines, effective traps for focusing free gas, are characterized by both normal faults and vertical zones of hydraulic fracturing within the hydrate stability zone. The normal faults form as a result of sediment layer folding and gravitational collapse of ridges during uplift. We document both longitudinal (ridge-parallel) and transverse (ridge-perpendicular) extensional structures (normal faults and elongated hydraulic fracture zones) in the sub-seafloor of anticlinal ridges. Intriguingly, gas flow through ridges close to the deformation front of the wedge exploits longitudinal structures, while ridges further inboard are characterized by gas flow along transverse structures. This highlights pronounced changes in the shallow deformation of ridges in different parts of the wedge, associated with a switching of the least and intermediate principal stress directions. It is critical to understand these shallow stress fields because they control fluid flow patterns and methane seepage out of the seafloor. Key Points Gas migration through ridges occurs along both longitudinal (ridge-parallel) and transverse (ridge-perpendicular) zones of fracturing Shallow stress fields differ significantly between ridges, reflecting differences in ridge evolution and deformation Seismic reflection images of the base of gas hydrate stability and gas-water contacts are strongly affected by seismic frequency content

Description: Polar marine ecosystems are particularly vulnerable to the effects of climate change. Warming temperatures, freshening seawater, and disruption to sea-ice formation potentially all have cascading effects on food webs. New approaches are needed to better understand spatiotemporal interactions among biogeochemical processes at the base of Southern Ocean food webs. In marine systems, isoscapes (models of the spatial variation in the stable isotopic composition) of carbon and nitrogen have proven useful in identifying spatial variation in a range of biogeochemical processes, such as nutrient utilization by phytoplankton. Isoscapes provide a baseline for interpreting stable isotope compositions of higher trophic level animals in movement, migration, and diet research. Here, we produce carbon and nitrogen isoscapes across the entire Southern Ocean (〉40°S) using surface particulate organic matter isotope data, collected over the past 50 years. We use Integrated Nested Laplace Approximation-based approaches to predict mean annual isoscapes and four seasonal isoscapes using a suite of environmental data as predictor variables. Clear spatial gradients in δ13C and δ15N values were predicted across the Southern Ocean, consistent with previous statistical and mechanistic views of isotopic variability in this region. We identify strong seasonal variability in both carbon and nitrogen isoscapes, with key implications for the use of static or annual average isoscape baselines in animal studies attempting to document seasonal migratory or foraging behaviors.

Description: Magmatic activity that accompanied the collision between Arabia and Eurasia at ∼27 Ma, provides unique opportunities for understanding the triggers and magma reservoirs for collisional magmatism and its different styles in magmatic fronts and back-arcs. We present new ages and geochemical-isotopic results for magmatic rocks that formed during the collision between Arabia and Eurasia in NE Iran, which was a back-arc region to the main magmatic arcs of Iran. Our new zircon U-Pb ages indicate that collisional magmatism began at ∼24 Ma in the NE Iran back-arc, although magmatism in this area started in the Late Cretaceous time and continued until the Pleistocene. The collisional igneous rocks are characteristically bimodal, and basaltic-andesitic and dacitic-rhyolitic components show significant isotopic differences; εNd(t) = +4.4 to +7.4 and εHf(t) = +5.4 to +9.5 for mafic rocks and εNd(t) = +0.2 to +8.4 and εHf(t) = +3.4 to +12.3 for silicic rocks. The isotopic values and modeling suggest that fractional crystallization and assimilation-fractional crystallization played important roles in the genesis of felsic rocks in the NE Iran collisional zone. Trace element and isotopic modeling further emphasize that the main triggers of the magmatism in NE Iran comprise a depleted to the enriched mantle and the Cadomian continental crust of Iran. Our results also emphasize the temporal magmatic variations in the NE Iran back-arc from Late Cretaceous to Pleistocene.

Description: Sources of uncertainty (i.e., internal variability, model and scenario) in Atlantic Niño variability projections were quantified in 49 models participating in the Coupled Model Intercomparison Phases 5 (CMIP5) and 6 (CMIP6). By the end of the twenty‐first century, the ensemble mean change in Atlantic Niño variability is −0.07 ± 0.10˚C, with 80% of CMIP models projecting a decrease, and representing a 16% reduction relative to the 1981–2005 ensemble mean. Models' projections depict a large spread, with variability changes ranging from 0.23˚C to −0.50˚C. Internal variability is the main source of uncertainty until 2045 but model uncertainty dominates thereafter, eventually explaining up to 80% of the total uncertainty. The scenario uncertainty remains low (〈1%) throughout the twenty‐first century. The total uncertainty on Atlantic Niño variability projections is not improved when considering only CMIP models with a realistic zonal equatorial Atlantic sea surface temperature gradient. Plain Language Summary Sources of uncertainty (i.e., internal variability, model and scenario) in future projections of the Atlantic Niño variability were evaluated in global coupled models participating in the Coupled Model Intercomparison Phases 5 (CMIP5) and 6 (CMIP6). Relative to 1981–2005, models' projections depict a large spread, ranging from increasing Atlantic Niño variability by up to 0.23˚C to decreasing by up to −0.50˚C. By the end of the twenty‐first century, the ensemble mean Atlantic Niño variability change is −0.07 ± 0.10˚C with 80% of the global coupled models simulating a decrease. This change in the ensemble mean Atlantic Niño variability, relative to the period 1981–2005, represents a 16% reduction. During the first four decades of projection, the internal variability is the main contributor to the total uncertainty; thereafter model uncertainty dominates and explains up to 80% of the total uncertainty at the end of the twenty‐first century. The scenario uncertainty remains low (〈1%) throughout the twenty‐first century. The total uncertainty on Atlantic Niño variability projections is not improved when considering only CMIP models with a realistic zonal equatorial Atlantic sea surface temperature gradient. Key Points 80% of the CMIP models simulate a decrease of the Atlantic Niño variability at the end of the 21st century The model uncertainty explains about 80% of the total uncertainty on Atlantic Niño variability projections at the end of the 21st century Global warming signal is not detectable throughout scenarios due to large internal variability and model uncertainties

Description: Glacial-isostatic adjustment (GIA) is the key process controlling relative sea-level (RSL) and paleo-topography. The viscoelastic response of the solid Earth is controlled by its viscosity structure. Therefore, the appropriate choice of Earth structure for GIA models is still an important area of research in geodynamics. We construct 18 3D Earth structures that are derived from seismic tomography models and are geodynamically constrained. We consider uncertainties in 3D viscosity structures that arise from variations in the conversion from seismic velocity to temperature variations (factor r) and radial viscosity profiles (RVP). We apply these Earth models to a 3D GIA model, VILMA, to investigate the influence of such structure on RSL predictions. The variabilities in 3D Earth structures and RSL predictions are investigated for globally distributed sites and applied for comparisons with regional 1D models for ice center (North America, Antarctica) and peripheral regions (Central Oregon Coast, San Jorge Gulf). The results from 1D and 3D models reveal substantial influence of lateral viscosity variations on RSL. Depending on time and location, the influence of factor r and/or RVP can be reverse, for example, the same RVP causes lowest RSL in Churchill and largest RSL in Oregon. Regional 1D models representing the structure beneath the ice and 3D models show similar influence of factor r and RVP on RSL prediction. This is not the case for regional 1D models representing the structure beneath peripheral regions indicating the dependence on the 3D Earth structure. The 3D Earth structures of this study are made available.

Description: The Los Chocoyos (14.6°N, 91.2°W) supereruption happened ∼75,000 years ago in Guatemala and was one of the largest eruptions of the past 100,000 years. It emitted enormous amounts of sulfur, chlorine, and bromine, with multi‐decadal consequences for the global climate and environment. Here, we simulate the impact of a Los Chocoyos‐like eruption on the quasi‐biennial oscillation (QBO), an oscillation of zonal winds in the tropical stratosphere, with a comprehensive aerosol chemistry Earth System Model. We find a ∼10‐year disruption of the QBO starting 4 months post eruption, with anomalous easterly winds lasting ∼5 years, followed by westerlies, before returning to QBO conditions with a slightly prolonged periodicity. Volcanic aerosol heating and ozone depletion cooling leads to the QBO disruption and anomalous wind regimes through radiative changes and wave‐mean flow interactions. Different model ensembles, volcanic forcing scenarios and results of a second model back up the robustness of our results.

Description: The Havre Trough backarc basin in the southwest Pacific is in the rifting stage of development. We distinguish five types of basalt there based on their amount and kind of slab component: backarc basalts (BAB) with little or no slab component; modified BAB (mBAB) with slight amounts; reararc (RA) with more; remnants of the pre‐existing arc (Colville Ridge horsts, CRH); and arc front volcanoes within the Havre Trough. Previous sub‐arc mantle is quickly removed and replaced by more fertile mantle with less slab component. The ambient mantle is “Pacific” isotopically, and more enriched in Nb/Yb and Nd and Hf isotope ratios north of the Central Kermadec Discontinuity at 32°S than to the south. The contrast may reflect inheritance in the south of mantle that was depleted during spreading that formed the southern South Fiji Basin, and a higher degree of melting because of a wetter slab‐derived flux. The slab component also differs along strike, more like a dry melt in the north and a super‐critical fluid in the south. The mass fraction of slab component increases southward in the backarc as well as the arc front. Reararc volcanoes have the most slab component (1‐2%) and form indistinct ridges at high angles to, and 〈50 km behind, frontal volcanoes. Backarc basalts have less and occur throughout the basin. Slab components are distributed further into the backarc, and more irregularly, during the rifting than spreading stage of backarc basin development. The rifting stage is disorganized geochemically as well as spatially.

Description: We have successfully constructed and tested a new, portable, Hybrid Lister‐Outrigger (HyLO) probe designed to measure geothermal gradients in submarine environments. The lightweight, low‐cost probe is 1‐3 m long, contains 4‐12 semiconductor temperature sensors that have a temperature resolution of 0.002 oC, a sample rate of 〈2 seconds, and a maximum working depth of ~2100 meters below sea level (mbsl). Probe endurance is continuous via ship‐power to water depths of ~700 mbsl, or up to ~1 week on batteries in depths 〉500 mbsl. Data are saved on solid‐state disks, transferred directly to the ship during deployment via a data cable, or transmitted via Bluetooth when the probe is at the sea surface. The probe contains an accelerometer to measure tilt, and internal pressure, temperature, and humidity gauges. Key advantages of this probe include (1) near‐real time temperature measurements and data transfer; (2) a low‐cost, transportable, and lightweight design; (3) easy and rapid two‐point attachment to a gravity corer, (4) short (3‐5 minute) thermal response times; (5) high temporal/spatial resolution and (6) longer deployment endurance compared to traditional methods. We successfully tested the probe both in lakes and during sea trials in May 2019 offshore Montserrat during the R/V Meteor Cruise 154/2. Probe‐measured thermal gradients were consistent with seafloor ocean‐drilling temperature measurements. Ongoing probe improvements include the addition of real‐time bottom‐camera feeds and long‐term (6‐12 month) deployment for monitoring. Key Points - We have designed, developed, and tested a low‐cost, portable hybrid Lister‐type probe to measure shallow thermal gradients - The probe consists of lightweight, quickly interchangeable/expendable components deployable to 2100 meters depth - The probe provides high vertical and temporal temperature resolution and rapid data transmission, reducing down‐time

Description: The giant tsunami that swept the Pacific from Alaska to Antarctica in 1946, was generated along one of three Alaska Trench instrumentally recorded aftershock areas following great and giant earthquakes. Aftershock areas were investigated during the past decade with multibeam bathymetry, OBS wide‐angle seismic, reprocessed legacy and new seismic reflection images. Summarized and updated here are previous papers and additional data. Tectonic structures collocated with aftershock area boundaries indicate possible lengths of rupture in future great earthquakes. NE aftershock area boundaries relate to subducted lower plate structures whereas the SW zone upper plate retains Beringian structural relicts. The lower to middle slope transition separating a stronger continental framework rock from a weaker accreted prism occurs along splay fault zones previously interpreted as backstops in seismic images. Damage zones along splay faults are generally 1 km wide dipping typically 21°. Splays form slip paths from the plate interface to the seafloor much shorter than the 3° to 4° dipping plate interface beneath the frontal prism. Associated seafloor vent structures indicate overpressured fluids at depth. Splay fault dip and its rigid hanging wall impart greater seafloor uplift than the accreted prism per unit of slip making them effective tsunami generators. Backstop splay fault zones run along the entire Alaska Trench. Beneath the frontal prism, active bend faults add rugosity to the plate interface and km high relief is commonly imaged in reprocessed legacy and new seismic data. The 1946 Unimak great (M8.6) earthquake epicenter is located near the backstop splay fault zone.

Description: Submarine landslides can be several orders of magnitude larger than their terrestrial counterparts and can pose significant hazards across entire ocean basins. The landslide failure mechanism strongly controls the associated tsunami hazard. The Tampen Slide offshore Norway is one of the largest landslides on Earth but remains poorly understood due to its subsequent burial beneath up to 450 m of sediments. Here, we use laterally extensive (16,000 km2), high‐resolution processed 3D seismic reflection data to characterize the upper Tampen Slide. We identify longitudinal (downslope, movement‐parallel) chutes and ridges that are up‐to‐40 m high, as well as extensional and compressional (cross‐slope) ridges. This is the first time that longitudinal ridges of such size have been imaged in a deep marine setting. The first phase of the Tampen Slide involved the simultaneous translation of over 720 km3 of sediments along a single failure plane. This was followed by spreading along the head‐ and sidewall, and the formation of a retrogressive debris flow and slump, the volumes of which are insignificant compared to the first failure. The process responsible for movement of such a large sediment volume along a single glide plane differs significantly from that of other passive margin megaslides, which typically comprise numerous smaller landslides that fail retrogressively along multiple glide planes. The trigger mechanism (e.g. an earthquake), the presence of mechanically strong obstructions (e.g. igneous topographical high), and the number and location of weak layers may be key factors that determine whether megaslides develop along a single plane or retrogressively.

Description: The development of the South Asian monsoon (SAM) and Himalaya‐Tibetan Plateau uplift were closely intertwined with some studies suggesting that uplift initiated the monsoon whereas others link tectonics with monsoon‐controlled exhumation. Silicate weathering controls atmospheric CO2 on geological timescales resulting in a large potential for monsoon strength and the Himalayan orogeny to influence global climate but detailed records of SAM‐induced weathering on million year (Myr) timescales are lacking. Here, we present radiogenic Sr, Nd, and Pb isotope compositions of clay minerals produced by silicate weathering and transported to the central Bay of Bengal. The radiogenic isotope data exhibit a relatively small range and demonstrate a remarkably consistent mixture of sources dominated by Himalayan rocks and the Indo‐Burman ranges, which consist of sediments derived from the Himalayas. This suggests that the spatial pattern of regional weathering, which today is highest in the regions of strongest monsoon rains, has persisted in a similar form for the last 27 Myrs. A pronounced increase in primary clay mineral abundance (from 9% to 22%) coincident with global cooling 13.9 Myrs ago points to a shift in the weathering regime given that the clay provenance did not change dramatically. Relatively weaker chemical weathering intensity during the mid and late Miocene cooling suggests increased aridity and changes in the large scale atmospheric circulation in the SAM domain. The establishment of the dry winter monsoon season during the mid and late Miocene may have caused this shift in the weathering regime and can reconcile much of the contrasting evidence for SAM initiation.

Description: The aftershock distribution of the 2014 Mw 8.1 Iquique earthquake offshore northern Chile, identified from a long‐term deployment of ocean bottom seismometers installed eight months after the mainshock, in conjunction with seismic reflection imaging, provides insights into the processes regulating the up‐dip limit of coseismic rupture propagation. Aftershocks up‐dip of the mainshock hypocenter frequently occur in the upper plate and are associated with normal faults identified from seismic reflection data. We propose that aftershock seismicity near the plate boundary documents subduction erosion that removes mass from the base of the wedge and results in normal faulting in the upper plate. The combination of very little or no sediment accretion and subduction erosion over millions of years has resulted in a very weak and aseismic frontal wedge. Our observations thus link the shallow subduction zone seismicity to subduction erosion processes that control the evolution of the overriding plate. Key Points: - We investigate structure and seismicity at the up-dip end of the 2014 Iquique earthquake rupture using amphibious seismic data. - Seismicity up-dip of the 2014 Iquique earthquake occurs over a broad range likely interpreted to be related to the basal erosion processes. - Coseismic stress changes and aftershocks activate extensional faulting of the upper plate and subduction erosion.

Description: Key Points: • diatom sedaDNA composition is concordant with reconstructed sea-ice dynamics, SSTs and subsurface salinities over the past millenia • loss of diatom richness at ~11.1 cal kyr BP is possibly a consequence of increased freshwater input from Kamchatka • shifts of potential ecotypes of species from the genus Chaetoceros with changing environmental conditions We traced diatom composition and diversity through time using diatom derived sedimentary ancient DNA (sedaDNA) from eastern continental slope sediments off Kamchatka (North Pacific) by applying a short, diatom‐specific marker on 63 samples in a DNA metabarcoding approach. The sequences were assigned to diatoms that are common in the area and characteristic of cold water. SedaDNA allowed us to observe shifts of potential lineages from species of the genus Chaetoceros that can be related to different climatic phases, suggesting that pre‐adapted ecotypes might have played a role in the long‐term success of species in areas of changing environmental conditions. These sedaDNA results complement our understanding of the long‐term history of diatom assemblages and their general relationship to environmental conditions of the past. Sea‐ice diatoms (Pauliella taeniata (Grunow) Round & Basson, Attheya septentrionalis (Østrup) R.M.Crawford and Nitzschia frigida (Grunow)) detected during the late glacial and Younger Dryas are in agreement with previous sea‐ice reconstructions. A positive correlation between pennate diatom richness and the sea‐ice proxy IP25 suggests that sea ice fosters pennate diatom richness, whereas a negative correlation with June insolation and temperature points to unfavorable conditions during the Holocene. A sharp increase in proportions of freshwater diatoms at ∼11.1 cal kyr BP implies the influence of terrestrial runoff and coincides with the loss of 42% of diatom sequence variants. We assume that reduced salinity at this time stabilized vertical stratification which limited the replenishment of nutrients in the euphotic zone.

Description: An extensive data set of biogenic silica (BSi) fluxes is presented for the Peruvian oxygen minimum zone (OMZ) at 11ºS and 12ºS. Each transect extends from the shelf to the upper slope (∼1000 m) and dissects the permanently anoxic waters between ∼200 – 500m water depth. BSi burial (2100 mmol m‐2 yr‐1) and recycling fluxes (3300 mmol m‐2 yr‐1) were highest on the shelf with mean preservation efficiencies (34±15%) that exceed the global mean of 10 – 20%. BSi preservation was highest on the inner shelf (up to 56%), decreasing to 7% and 12% under anoxic waters and below the OMZ, respectively. The data suggest that the main control on BSi preservation is the rate at which reactive BSi is transported away from undersaturated surface sediments by burial and bioturbation to the underlying saturated sediment layers where BSi dissolution is thermodynamically and/or kinetically inhibited. BSi burial across the entire Peruvian margin between 3ºS to 15ºS and down to 1000m water depth is estimated to be 0.1 – 0.2 Tmol yr‐1; equivalent to 2 – 7% of total burial on continental margins. Existing global data permit a simple relationship between BSi rain rate to the seafloor and the accumulation of unaltered BSi, giving the possibility to reconstruct rain rates and primary production from the sediment archive in addition to benthic Si turnover in global models.

Description: Crystallographic preferred orientation (CPO) and the associated seismic anisotropy of serpentinites are important factors for the understanding of tectonic settings involving hydrated Earth´s mantle, for example, at slow-spreading mid-ocean ridges. CPO of lizardite and magnetite in low-grade metamorphic serpentinites from the Atlantis Massif oceanic core complex (Mid-Atlantic Ridge, 30°N) were determined using synchrotron high energy X-ray diffraction in combination with Rietveld texture analysis. Serpentinite mesh structures show weak CPO while deformed samples show a single (0001) maximum perpendicular to the foliation. Seismic anisotropies calculated from CPO show up to 〉11% anisotropy for compressional waves (Vp) and shear wave splitting up to 0.38 km/s in the deformed samples. This indicates that deformation in shear zones controls elastic anisotropy and highlights its importance in defining the seismic signature of hydrated upper mantle.

Description: Taupō volcano, New Zealand, is a large caldera volcano that has been highly active through the Holocene. It most recently erupted ∼1800 years ago but there have been multiple periods of historic volcanic unrest. We use seismological and geodetic analysis to show that in 2019 Taupō underwent a period of unrest characterised by increased seismic activity through multiple swarms and was accompanied by ground deformation within the caldera. The earthquakes, which include non-double-couple events, serve to outline an aseismic zone beneath the most recent eruptive vents. This aseismic zone is coincident with an inflating source, based on forward modelling of ground deformation data. We infer that this aseismic and deforming region delineates the location of the present day magma reservoir that is ≥250 km3 in volume and has a melt fraction of 〉 20–30%, inhibiting seismic activity. Our analysis shows that the 2019 unrest at Taupō was volcanic in nature and origin, demonstrating that this is an active and potentially hazardous volcano, and that improving our monitoring and understanding of its behaviour is important.

Description: Microbial degradation of dissolved organic matter (DOM) contributes to the formation and preservation of oxygen minimum zones (OMZs) in the ocean, but information on the spatial distribution and molecular composition of DOM in OMZ regions is scarce. We quantified molecular components of DOM that is, dissolved amino acids (DAA) and dissolved combined carbohydrates (DCCHO), in the upwelling region off Peru. We found the highest concentrations of DCCHO in fully oxygenated surface waters steeply declining at shallow depth. The highest DAA concentrations were observed close to the surface also, but attenuation of DAA concentration over depth was less pronounced. Compositional changes of DCCHO were strongest within more oxygenated waters. Compositional changes of DAA were also evident under suboxic conditions (〈5 µmol O2 kg−1) and indicated bacterial peptide degradation. Moreover, specific free amino acids (alanine and threonine) were enhanced within suboxic waters, pointing to a potential production of dissolved organic nitrogen under suboxic conditions. Our results therewith suggest that deoxygenation supports a spatial decoupling of DCCHO and DAA production and degradation dynamics and give new insights to carbon and nitrogen cycling in the OMZ off Peru.

Description: We use a 30-year time series (1986–2016) of dichlorodifluoromethane (CFC-12) concentrations with a refined transit time distribution (TTD) method, to estimate the temporal variation of anthropogenic carbon (Cant) in the Central Labrador Sea. We determined that the saturation of CFC-12 and sulfur hexafluroide (SF6) in newly-formed Labrador Sea Water had departed significantly from 100% and varied systematically with time. Multiple linear regression of the time-varying saturation, with the tracer's atmospheric growth rate and the wintertime mixed layer depth as independent variables, allowed reconstruction of the saturation history of CFC-12 and SF6 in wintertime surface waters, which was implemented in the TTD method. Use of the time-varying saturation for CFC-12 gave Cant concentrations ∼7 μmol kg−1 larger than estimates obtained assuming a constant saturation of 100%. The resulting Cant column inventories were ∼20% larger and displayed lower interannual variability compared to conventional TTD-based estimates. The column inventory of Cant increased at an average rate of 1.8 mol m−2 y−1 over the 30-year period. However, the accumulation rate of Cant was higher than this average in the early 1990s and since 2013, whereas inventories remained almost unchanged between 2003 and 2012. The variation in the Cant accumulation rate is shown to be linked to temporal variability in the relative layer thickness of the annually ventilated Labrador Sea Water and the underlying Deep Intermediate Water. The non-steady Cant accumulation highlights the importance of sampling frequency, especially in regions of variable deep mixing and high carbon inventories, and potential misinterpretation of Cant dynamics

Description: The presence of clouds in the Arctic regulates the surface energy budget (SEB) over the sea-ice surface and the ice-free ocean. Following several previous field campaigns, the cloud-radiation relationship, including cloud vertical structure and phase, has been elucidated; however, modeling of this relationship has matured slowly. In recognition of the recent decline in the Arctic sea-ice extent, representation of the cloud system in numerical models should consider the effects of areas covered by sea ice and ice-free areas. Using an in situ stationary meteorological observation data set obtained over the ice-free Arctic Ocean by the Japanese Research Vessel Mirai (September 2014), coordinated evaluation of six regional climate models (RCMs) with nine model runs was performed by focusing on clouds and the SEB. The most remarkable findings were as follows: (1) reduced occurrence of unstable stratification with low-level cloud water in all models in comparison to the observations, (2) significant differences in cloud water representations between single- and double-moment cloud schemes, (3) extensive differences in partitioning of hydrometeors including solid/liquid precipitation, and (4) pronounced lower-tropospheric air temperature biases. These issues are considered as the main sources of SEB uncertainty over ice-free areas of the Arctic Ocean. The results from a coupled RCM imply that the SEB is constrained by both the atmosphere and the ocean (and sea ice) with considerable feedback. Coordinated improvement of both stand-alone atmospheric and coupled RCMs would promote a more comprehensive and improved understanding of the Arctic air-ice-sea coupled system.

Description: Carbonate lithologies host considerable quantities of the Earth’s freshwater resources and partially supply a quarter of the global population with drinkable water. Carbonates constitute substantial amounts of the global coastlines, yet it is not known if and how they can sustain freshened groundwater offshore. Here, we use controlled source electromagnetic, seismic reflection, and core sample data to derive a lithological model for the eastern margin of the Maltese Islands and identify four distinct resistivity anomalies within the Upper Coralline Limestone, Globigerina Limestone, and Blue Clay formations. The anomalies hosted in the former are likely associated to low porosities, whereas the anomaly within the latter is indicative of pore fluid freshening. Hydrogeological modeling suggests that freshened pore fluids, emplaced during sea-level lowstands and preserved in low permeability units, are potentially still found within carbonate shelves. However, resource potential is low due to its relict nature and low permeability host environment. Key points ● Geophysical data and hydrogeological modeling are applied to detect offshore freshened groundwater in a semi-arid carbonate setting ● Globigerina Limestone and Blue Clay located offshore SE Malta likely host a disconnected offshore freshened groundwater body ● The resistive anomalies within the Upper Coralline Limestone are interpreted as localized porosity variations.Accepted Article This article is protected by copyright. All rights reserved. ● This OFG was emplaced during sea-level lowstands and preserved in low permeability units

Description: The Arctic Ocean receives a large supply of dissolved organic matter (DOM) from its catchment and shelf sediments, which can be traced across much of the basin’s upper waters. This signature can potentially be used as a tracer. On the shelf, the combination of river discharge and sea-ice formation, modifies water densities and mixing considerably. These waters are a source of the halocline layer that covers much of the Arctic Ocean, but also contain elevated levels of DOM. Here we demonstrate how this can be used as a supplementary tracer and contribute to evaluating ocean circulation in the Arctic. A fraction of the organic compounds that DOM consists of fluoresce and can be measured using in-situ fluorometers. When deployed on autonomous platforms these provide high temporal and spatial resolution measurements over long periods. The results of an analysis of data derived from several Ice Tethered Profilers (ITPs) offer a unique spatial coverage of the distribution of DOM in the surface 800m below Arctic sea-ice. Water mass analysis using temperature, salinity and DOM fluorescence, can clearly distinguish between the contribution of Siberian terrestrial DOM and marine DOM from the Chukchi shelf to the waters of the halocline. The findings offer a new approach to trace the distribution of Pacific waters and its export from the Arctic Ocean. Our results indicate the potential to extend the approach to separate freshwater contributions from, sea-ice melt, riverine discharge and the Pacific Ocean. Key Points: Arctic surface waters with comparable temperature and salinity have contrasting in situ dissolved organic matter fluorescence. Organic matter fluorescence can tracklow salinity waters feeding into the Transpolar Drift and haloclinelayers. Siberian and Chukchishelf waters can be separated based on their fluorescence to salinity relationship

Description: Biological productivity in the ocean directly influences the partitioning of carbon between the atmosphere and ocean interior. Through this carbon cycle feedback, changing ocean productivity has long been hypothesized as a key pathway for modulating past atmospheric carbon dioxide levels and hence global climate. Because phytoplankton preferentially assimilate the light isotopes of carbon and the major nutrients nitrate and silicic acid, stable isotopes of carbon (C), nitrogen (N), and silicon (Si) in seawater and marine sediments can inform on ocean carbon and nutrient cycling, and by extension the relationship with biological productivity and global climate. Here, we compile water column C, N, and Si stable isotopes from GEOTRACES-era data in four key ocean regions to review geochemical proxies of oceanic carbon and nutrient cycling based on the C, N, and Si isotopic composition of marine sediments. External sources and sinks as well as internal cycling (including assimilation, particulate matter export, and regeneration) are discussed as likely drivers of observed C, N, and Si isotope distributions in the ocean. The potential for C, N, and Si isotope measurements in sedimentary archives to record aspects of past ocean C and nutrient cycling is evaluated, along with key uncertainties and limitations associated with each proxy. Constraints on ocean C and nutrient cycling during late Quaternary glacial-interglacial cycles and over the Cenozoic are examined. This review highlights opportunities for future research using multielement stable isotope proxy applications and emphasizes the importance of such applications to reconstructing past changes in the oceans and climate system.

Description: Focused gas migration through the gas hydrate stability zone in vertical gas conduits is a global phenomenon. The process can lead to concentrated gas hydrate formation and seafloor gas seepage, which influences seafloor biodiversity and ocean biogeochemistry. However, much is unknown about how gas and gas hydrate co-exist within and around gas conduits. We present seismic imaging of the gas hydrate system beneath a four-way closure anticlinal ridge at New Zealand's southern Hikurangi subduction margin. Gas has accumulated beneath the base of gas hydrate stability to a thickness of up to ∼240 m, which has ultimately led to hydraulic fracturing and propagation of a vertical gas conduit to the seafloor. Despite the existence of an array of normal faults beneath the ridge, these structures are not exploited as long-range gas flow conduits. Directly beneath the conduit, and extending upward from the regional base of gas hydrate stability, is a broad zone characterized by both negative- and positive-polarity reflections. We interpret this zone as a volume of sediment hosting both gas hydrate and free gas, that developed due to partial gas trapping beneath a mass transport deposit. Similar highly reflective zones have been identified at the bases of other gas conduits, but they are not intrinsic to all gas conduits through gas hydrate systems. We suggest that pronounced intervening sealing units within the gas hydrate stability zone determine whether or not they form.

Description: Cyclonic ocean eddies drive upwelling of deep waters enhanced in nutrients, which can elevate phytoplankton productivity. At mid‐latitudes in the North Atlantic, satellite images show enhanced chlorophyll‐a associated with eddies. However, surface macronutrient concentrations are often not fully depleted in this region, implying enhanced macronutrient supply is not the primary control. We conducted high resolution sampling through two mid‐latitude Atlantic eddies in late spring, located 800 and 350 km east of the Newfoundland Grand Banks. Waters outside of both eddies had unused residual macronutrients, low dissolved iron, and iron‐stressed phytoplankton. Inside both eddies, plankton biomass was higher and macronutrient concentrations lower. However, full macronutrient drawdown and an absence of iron stress were only present in the eddy nearer the continental shelf. From these two examples, iron supply and proximity to shelf iron sources appear to be important factors regulating productivity and macronutrient utilization in mid‐latitude North Atlantic cyclonic eddies.

Description: A GEOMAR (Kiel, Germany) research team has developed a passive electric field acquisition system for Autonomous Underwater Vehicles (AUVs) to optimize seafloor massive sulfides exploration. This sensor was made of two perpendicular and horizontal pairs of electrodes, and was successfully tested over active basalt-hosted hydrothermal site TAG (26°N, Mid-Atlantic Ridge) and several inactive sites in its vicinity. The resulting data underline the efficiency of combining deep-sea electric and magnetic measurements for searching for active and inactive hydrothermal vent fields. With these datasets, it becomes possible to determine the geological nature of the targets and to constrain the characteristics of fluid circulation at depth without involving costly and invasive underwater tools such as Remotely Operated Vehicles or even manned submersibles to collect samples. Data analysis also revealed that AUV attitude variations induce distortions of the electric signal. These distortions start prevailing for dives at altitudes higher than 90 m above the seafloor, as the distance between the AUV becomes too important to guarantee that the signal produced by the geological target still dominates. To improve the acquisition system and reduce the overall noise, we discuss solutions that limit the impact of such attitude variations. These solutions consist of minor adjustments, such as masts at AUVs stern to tow damping electrodes arrays. In such configurations, we believe that deep-sea passive electric measurements combined with high-resolution magnetic measurements can become a highly efficient seafloor exploration tool, including for sulfide deposits associated with inactive hydrothermal systems.

Description: The Milos, Christiana-Santorini-Kolumbo (CSK) and Kos-Yali-Nisyros (KYN) volcanic complexes of the Aegean Volcanic Arc have repeatedly produced highly explosive eruptions from at least ∼360 ka into historic times and still show recent unrest. We present the marine tephra record from an array of 50, up to 7.4 m long, sediment cores along the arc collected in 2017 during RV Poseidon cruise POS513, which complements earlier work on distal to ultra-distal eastern Mediterranean sediment cores. A unique set of glass-shard trace element (LA-ICPMS) compositions complements our major element (EMP) data on 220 primary ash layers and 40 terrestrial samples to support geochemical fingerprinting for correlations with 19 known tephras from all three volcanic complexes and with the 39 ka Campanian Ignimbrite from the Campi Flegrei, Italy. The correlations include eleven eruptions from CSK (Kameni, Kolumbo 1650, Minoan, Cape Riva, Cape Tripiti, Upper Scoriae 1 and 2, Middle Pumice, Cape Thera, Lower Pumice, Cape Therma 3). We identify a previously unknown widespread tephra from a plinian eruption on Milos (Firiplaka Tephra). Near the KYN we correlate marine tephras with the Kos Plateau Tuff, the Yali 1 and Yali 2 tephras, and the Upper and Lower Pumice on Nisyros. Between these two major tephras, we found two tephras from Nisyros not yet observed on land. The four Nisyros tephras form a systematic trend toward more evolved magma compositions. In the companion paper we use the tephrostratigraphic framework established here to constrain new eruption ages and magnitudes as a contribution to volcanic hazard assessment.

Description: Arctic sea ice cover has been steeply declining since the onset of satellite observations in the late 1970s. However, the available annually resolved sea ice data prior to this time are limited. Here, we evaluated the suitability of annual trace element (Mg/Ca) ratios and growth increments from the long-lived annual increment-forming benthic coralline red alga, Clathromorphum compactum, as high-resolution sea ice cover archive. It has previously been shown that growth of C. compactum is strongly light controlled and therefore greatly limited during polar night and underneath sea ice cover. We compare algal data from 11 sites collected throughout the Canadian Arctic, Greenland and Svalbard, with satellite sea ice data. Our results suggested that algal growth anomalies most often produced better correlations to sea ice concentration than Mg/Ca ratios or when averaging growth and Mg/Ca anomalies. High Arctic regions with persistently higher sea ice concentrations and shorter ice-free seasons showed strongest correlations between algal growth anomalies and satellite sea ice concentration over the study period (1979-2015). At sites where ice breakup took place prior to the return of sufficient solar irradiance, algal growth was most strongly tied to a combination of solar irradiance and other factors such as temperature, suspended sediments, phytoplankton blooms and cloud cover. These data are the only annually resolved in situ marine proxy data known to date and are of utmost importance to gain a better understanding of the sea ice system and to project future sea ice conditions.

Description: The West Iberia margin is the focus of intense research since the 1980s, with some of the most exemplary geophysical cross-sections and drilling expeditions. Those data sets have been used to create conceptual models of rifting used as a template to interpret margins worldwide. We present two collocated ∼350 km long lines of multi-channel seismic (MCS) streamer data and wide-angle seismic (WAS) data collected across the Tagus Abyssal Plain (TAP). We use travel-times of first arrivals identified at WAS and reflected seismic phases identified at both WAS and MCS records to jointly invert for the P wave velocity (Vp) distribution and the geometry of a sediment unconformity, the top of the basement, and the Moho boundary. The Vp model shows that the TAP basement is more complex than previously inferred, presenting abrupt boundaries between five domains. Domain I under the foot of the slope and Domain III under the abyssal plain display Vp values and gradients of thin continental crust. In between, Domain II displays a steep Vp gradient and high Vp values at shallow depth that support that basement is made of exhumed partly serpentinized mantle. Domain IV and Domain V, further oceanward, have oceanic crust Vp structure. The new results support an unanticipated complex rift history during the initial separation of Iberia and America. We propose a geodynamic scenario characterized by two phases of extension separated by a jump of the locus of extension, caused by the northward propagation of the oceanic spreading center during the J-anomaly formation, which terminated continental rifting.

Description: When organic matter from thawed permafrost is released, the sources and sinks of greenhouse gases (GHGs), like carbon dioxide (CO2) and methane (CH4) in Arctic rivers will be influenced in the future. However, the temporal variation, environmental controls, and magnitude of the Arctic riverine GHGs are largely unknown. We measured in situ high temporal resolution concentrations of CO2, CH4, and oxygen (O2) in the Ambolikha River in northeast Siberia between late June and early August 2019. During this period, the largely supersaturated riverine CO2 and CH4 concentrations decreased steadily by 90% and 78%, respectively, while the O2 concentrations increased by 22% and were driven by the decreasing water temperature. Estimated gas fluxes indicate that during late June 2019, significant emissions of CO2 and CH4 were sustained, possibly by external terrestrial sources during flooding, or due to lateral exchange with gas-rich downstream-flowing water. In July and early August, the river reversed its flow constantly and limited the water exchange at the site. The composition of dissolved organic matter and microbial communities analyzed in discrete samples also revealed a temporal shift. Furthermore, the cumulative total riverine CO2 emissions (36.8 gC-CO2 m−2) were nearly five times lower than the CO2 uptake at the adjacent floodplain. Emissions of riverine CH4 (0.21 gC-CH4 m−2) were 16 times lower than the floodplain CH4 emissions. Our study revealed that the hydraulic connectivity with the land in the late freshet, and reversing flow directions in Arctic streams in summer, regulate riverine carbon replenishment and emissions.

Description: Although submarine landslides have been studied for decades, a persistent challenge is the integration of diverse geoscientific datasets to characterise failure processes. We present a core-log-seismic integration study of the Tuaheni Landslide Complex to investigate intact sediments beneath the undeformed seafloor as well as post-failure landslide deposits. Beneath the undeformed seafloor are coherent reflections underlain by a weakly-reflective and chaotic seismic unit. This chaotic unit is characterised by variable shear strength that correlates with density fluctuations. The basal shear zone of the Tuaheni landslide likely exploited one (or more) of the low shear strength intervals. Within landslide deposits is a widespread “Intra-debris Reflector”, previously interpreted as the landslide’s basal shear zone. This reflector is a subtle impedance drop around the boundary between upper and lower landslide units. However, there is no pronounced shear strength change across this horizon. Rather, there is a pronounced reduction in shear strength ∼10-15 m above the Intra-debris Reflector that presumably represents an induced weak layer that developed during failure. Free gas accumulates beneath some regions of the landslide and is widespread deeper in the sedimentary sequence, suggesting that free gas may have played a role in pre-conditioning the slope to failure. Additional pre-conditioning or failure triggers could have been seismic shaking and associated transient fluid pressure. Our study underscores the importance of detailed core-log-seismic integration approaches for investigating basal shear zone development in submarine landslides.

Description: In this study, we used stable isotopes of oxygen (δ18O), deuterium (δD), and dissolved inorganic carbon (δ13CDIC) in combination with temperature, salinity, oxygen and nutrient concentrations to characterize the coastal (71-78 °W) and an oceanic (82-98 °W) water masses (SAAW-Subantarctic Surface Water; STW-Subtropical Water; ESSW-Equatorial Subsurface water; AAIW-Antarctic Intermediate Water; PDW-Pacific Deep Water) of the Southeast Pacific (SEP). The results show that δ18O and δD can be used to differentiate between SAAW-STW, SAAW-ESSW and ESSW-AAIW. δ13CDIC signatures can be used to differentiate between STW-ESSW (oceanic section), SAAW-ESSW, ESSW-AAIW and AAIW-PDW. Compared with the oceanic section, our new coastal section highlights differences in both the chemistry and geometry of water masses above 1000 m. Previous paleoceanographic studies using marine sediments from the SEP continental margin used the present-day hydrological oceanic transect to compare against, as the coastal section was not sufficiently characterized. We suggest that our new results of the coastal section should be used for past characterizations of the SEP water masses that are usually based on continental margin sediment samples.

Description: Abrupt fluid emissions from shallow marine sediments pose a threat to seafloor installations like wind farms and offshore cables. Quantifying such fluid emissions and linking pockmarks, the seafloor manifestations of fluid escape, to flow in the sub-seafloor remains notoriously difficult due to an incomplete understanding of the underlying physical processes. Here, using a compositional multi-phase flow model, we test plausible gas sources for pockmarks in the south-eastern North Sea, which recent observations suggest have formed in response to major storms. We find that the mobilization of pre-existing gas pockets is unlikely because free gas, due to its high compressibility, damps the propagation of storm-induced pressure changes deeper into the subsurface. Rather, our results point to spontaneous appearance of a free gas phase via storm-induced gas exsolution from pore fluids. This mechanism is primarily driven by the pressure-sensitivity of gas solubility, and the appearance of free gas is largely confined to sediments in the vicinity of the seafloor. We show that in highly permeable sediments containing gas-rich pore fluids, wave-induced pressure changes result in the appearance of a persistent gas phase. This suggests that seafloor fluid escape structures are not always proxies for overpressured shallow gas and that periodic seafloor pressure changes can induce persistent free gas phase to spontaneously appear. Key Points - Storm-induced pressure changes can lead to spontaneous appearance of free gas phase near the seafloor - This process is driven by pressure-sensitive phase instabilities - This mechanism could help explain elusive gas sources in recently observed pockmarks in the North Sea

Description: Controls on the deformation pattern (shortening mode and tectonic style) of orogenic forelands during lithospheric shortening remain poorly understood. Here, we use high-resolution 2D thermomechanical models to demonstrate that orogenic crustal thickness and foreland lithospheric thickness significantly control the shortening mode in the foreland. Pure-shear shortening occurs when the orogenic crust is not thicker than the foreland crust or thick, but the foreland lithosphere is thin (〈70–80 km, as in the Puna foreland case). Conversely, simple-shear shortening, characterized by foreland underthrusting beneath the orogen, arises when the orogenic crust is much thicker. This thickened crust results in high gravitational potential energy in the orogen, which triggers the migration of deformation to the foreland under further shortening. Our models present fully thick-skinned, fully thin-skinned, and intermediate tectonic styles in the foreland. The first tectonics forms in a pure-shear shortening mode whereas the others require a simple-shear mode and the presence of thick (〉∼4 km) sediments that are mechanically weak (friction coefficient 〈∼0.05) or weakened rapidly during deformation. The formation of fully thin-skinned tectonics in thick and weak foreland sediments, as in the Subandean Ranges, requires the strength of the orogenic upper lithosphere to be less than one-third as strong as that of the foreland upper lithosphere. Our models successfully reproduce foreland deformation patterns in the Central and Southern Andes and the Laramide province.

Description: Germany 2050: For the first time Germany reached a balance between its sources of anthropogenic CO2 to the atmosphere and newly created anthropogenic sinks. This backcasting study presents a fictional future in which this goal was achieved by avoiding (∼645 Mt CO2), reducing (∼50 Mt CO2) and removing (∼60 Mt CO2) carbon emissions. This meant substantial transformation of the energy system, increasing energy efficiency, sector coupling, and electrification, energy storage solutions including synthetic energy carriers, sector-specific solutions for industry, transport, and agriculture, as well as natural-sink enhancement and technological carbon dioxide options. All of the above was necessary to achieve a net-zero CO2 system for Germany by 2050.

Description: High-resolution velocity models developed using full-waveform inversion (FWI) can image fine details of the nature and structure of the subsurface. Using a 3D FWI velocity model of hyper-thinned crust at the Deep Galicia Margin (DGM) west of Iberia, we constrain the nature of the crust at this margin by comparing its velocity structure with those in other similar tectonic settings. Velocities representative of both the upper and lower continental crust are present, but there is no clear evidence for distinct upper and lower crustal layers within the hyper-thinned crust. Our velocity model supports exhumation of the lower crust under the footwalls of fault blocks to accommodate the extension. We used our model to generate a serpentinization map for the uppermost mantle at the DGM, at a depth of 100 ms (∼340 m) below the S-reflector, a low-angle detachment that marks the base of the crust at this margin. We find a good alignment between serpentinized areas and the overlying major block bounding faults on our map, suggesting that those faults played an important role in transporting water to the upper mantle. Further, we observe a weak correlation between fault heaves and serpentinization beneath the hanging-wall blocks, indicating that serpentinization was controlled by complex faulting during rifting. A good match between topographic highs of the S and local highly serpentinized areas of the mantle suggests that the morphology of the S was affected by the volume-increasing process of serpentinization and deformation of the overlying crust. Key Points Exhumation of the lower crust under the footwall of the normal faults to accommodate extension Overlying faults in the crust control water transport to the mantle Topography of the S-reflector is affected by the serpentinization process and deformations of the overlying crust

Description: Stable barium isotopes are a potential proxy for riverine inputs into the ocean that reflect monsoon variability and climate change. However, dissolved Ba isotope (δ138BaDBa) geochemistry in river estuaries, a dynamic land to ocean transition zone, has rarely been systematically examined to date. Here, we show that significant Ba isotope fractionation occurs at near-zero salinities in the Yangtze and Pearl River Estuary, whereas conservative mixing dominates δ138BaDBa distributions beyond low salinities, which are well predicted by an ion exchange model. Elevated δ138BaDBa in the river endmember results from preferential removal of light Ba isotopes by adsorption to fluvial particles. Subsequently, δ138BaDBa rapidly drops to minimum signatures at increased salinities indicating particle desorption of isotopically light Ba. Nevertheless, the apparently conservative δ138BaDBa-salinity relationship beyond the low-salinity minimum in both estuaries provides a modern calibration for using Ba isotopes as a proxy for paleosalinity and river water inputs into the ocean.

Description: The SW Iberian margin is one of the most seismogenic and tsunamigenic areas in W-Europe, where large historical and instrumental destructive events occurred. To evaluate the sensitivity of the tsunami impact on the coast of SW Iberia and NW Morocco to the fault geometry and slip distribution for local earthquakes, we carried out a set of tsunami simulations considering some of the main known active crustal faults in the region: the Gorringe Bank (GBF), Marquês de Pombal (MPF), Horseshoe (HF), North Coral Patch (NCPF) and South Coral Patch (SCPF) thrust faults, and the Lineament South (LSF) strike-slip fault. We started by considering for all of them relatively simple planar faults featuring with uniform slip distribution; we then used a more complex 3D fault geometry for the faults constrained with a large 2D multichannel seismic dataset (MPF, HF, NCPF, and SCPF); and finally, we used various heterogeneous slip distributions for the HF. Our results show that using more complex 3D fault geometries and slip distributions, the peak wave height at the coastline can double compared to simpler tsunami source scenarios from planar fault geometries. Existing tsunami hazard models in the region use homogeneous slip distributions on planar faults as initial conditions for tsunami simulations and therefore underestimate tsunami hazard. Complex 3D fault geometries and non-uniform slip distribution should be considered in future tsunami hazard updates. The tsunami simulations also support the finding that submarine canyons attenuate the wave height reaching the coastline, while submarine ridges and shallow shelves have the opposite effect.

Description: We present two ∼150-km-long orthogonal 2D P-wave tomographic velocity models across and along the ridge axis of the ultraslow-spreading Southwest Indian Ridge at 64°30′E. Here, detachment faults largely accommodate seafloor accretion by mantle exhumation. The velocity models are constructed by inverting first arrival traveltimes recorded by 32 ocean bottom seismometers placed on the two profiles. The velocities increase rapidly with depth, from 3 to 3.5 km/s at the seafloor to 7 km/s at depths ranging from 1.5 to 6 km below the seafloor. The vertical gradient decreases for velocities 〉7 km/s. We suggest that changes in velocity with depth are related to changes in the degree of serpentinization and interpret the lithosphere to be composed of highly fractured and fully serpentinized peridotites at the top with a gradual downward decrease in serpentinization and pore space to unaltered peridotites. One active and five abandoned detachment faults are identified on the ridge-perpendicular profile. The active axial detachment fault (D1) shows the sharpest lateral change (horizontal gradient of ∼1 s–1) and highest vertical gradient (∼2 s–1) in the velocities. In the western section of the ridge-parallel profile, the lithosphere transitions from non-volcanic to volcanic over a distance of ∼10 km. The depth extent of serpentinization on the ridge-perpendicular profile ranges from ∼2 to 5 km, with the deepest penetration at the D1 hanging wall. On the ridge-parallel profile, this depth (∼2.5–4 km) varies less as the profile crosses the D1 hanging wall at ∼5–9 km south of the ridge axis.

Description: Porewater dissolved silicic acid (DSi) concentrations and stable Si isotope compositions (δ30Si) together with biogenic silica (bSiO2) contents of sediments in five sediment cores collected from the southern Mariana Trench are presented. These data suggest the occurrence of bSiO2 dissolution and concomitant authigenic clay formation in three bSiO2-bearing cores. A reaction-transport model constrained by the measured geochemical data was applied to quantify the rates of Si turnover. Model results predicted the greatest rates of both bSiO2 dissolution and authigenic clay formation at the trench axis core that displayed low bSiO2 contents and abundant detrital materials, suggesting that detrital materials may be a limiting factor for bSiO2 diagenesis. Model results further predicted that ∼40%–70% of DSi generated by bSiO2 dissolution is consumed by authigenic clay formation. This is the first study that demonstrates active silica diagenesis in the hadal realm and has implications for understanding benthic Si cycling in deep-sea settings. Key Points - Biogenic silica diagenesis was examined for the first time in hadal trench sediments - Availability of detrital materials may be a limiting factor for biogenic silica (bSiO2) diagenesis in bSiO2-rich sediments of the Mariana Trench - ∼40%–80% of dissolved silicic acid generated by bSiO2 dissolution is fixed by authigenic silicate formation

Description: Early Pleistocene Marine Isotope Stage (MIS)-31 (1.081–1.062 Ma) is a unique interval of extreme global warming, including evidence of a West Antarctic Ice Sheet (WAIS) collapse. Here we present a new 1,000-year resolution, spanning 1.110–1.030 Ma, diatom-based reconstruction of primary productivity, relative sea surface temperature changes, sea-ice proximity/open ocean conditions and diatom species absolute abundances during MIS-31, from the Scotia Sea (59°S) using deep-sea sediments collected during International Ocean Discovery Program (IODP) Expedition 382. The lower Jaramillo magnetic reversal (base of C1r.1n, 1.071 Ma) provides a robust and independent time-stratigraphic marker to correlate records from other drill cores in the Antarctic Zone of the Southern Ocean (AZSO). An increase in open ocean species Fragilariopsis kerguelensis in early MIS-31 at 53°S (Ocean Drilling Program Site 1,094) correlates with increased obliquity forcing, whereas at 59°S (IODP Site U1537; this study) three progressively increasing, successive peaks in the relative abundance of F. kerguelensis correlate with Southern Hemisphere-phased precession pacing. These observations reveal a complex pattern of ocean temperature change and sustained sea surface temperature increase lasting longer than a precession cycle within the Atlantic sector of the AZSO. Timing of an inferred WAIS collapse is consistent with delayed warmth (possibly driven by sea-ice dynamics) in the southern AZSO, supporting models that indicate WAIS sensitivity to local sub-ice shelf melting. Anthropogenically enhanced impingement of relatively warm water beneath the ice shelves today highlights the importance of understanding dynamic responses of the WAIS during MIS-31, a warmer than Holocene interglacia

Description: Along ultraslow spreading ridges melt is distributed unequally, but melt focusing guides melt away from amagmatic segments toward volcanic centers. An interplay of tectonism and magmatism is thought to control melt ascent, but the detailed process of melt extraction is not yet understood. We present a detailed image of the seismic velocity structure of the Logachev volcanic center and adjacent region along the Knipovich Ridge. With travel times of P- and S-waves of 3,959 earthquakes we performed a local earthquake tomography. We simultaneously inverted for source locations, velocity structure and the Vp/Vs-ratio. An extensive low velocity anomaly coincident with high Vp/Vs-ratios 〉1.9 lies underneath the volcanic center at depths of 10 km below sea level in an aseismic area. More shallow, tightly clustered earthquake swarms connect the anomaly to a shallow anomaly with high Vp/Vs-ratio beneath the basaltic seafloor. We consider the deep low-velocity anomaly to represent an area of partial melt from which melts ascent vertically to the surface and northwards into the adjacent segment. By comparing tomographic studies of the Logachev and Southwest Indian Ridge Segment-8 volcano we conclude that volcanic centers of ultraslow spreading ridges host spatially confined, circular partial melt areas below 10 km depth, in contrast to the shallow extended melt lenses along fast spreading ridges. Lateral feeding over distances of 35 km is possible at orthogonal spreading segments, but limited at the obliquely spreading Knipovich Ridge. Key Points - Active volcanic centers at ultraslow spreading ridges host deeper and more confined partial melt areas than faster spreading ridges - Earthquake swarms delineate melt ascent paths from the partial melt area to the surface - Lateral feeding at shallow depths into subordinate segments is prevented by ridge obliquity