ALBERT

Description: In the boreal summer of 2021, the equatorial Atlantic experienced the strongest warm event, that is, Atlantic Niño, since the beginning of satellite observations in the 1970s. Such events have far‐reaching impacts on large‐scale wind patterns and rainfall over the surrounding continents. Yet, developing a paradigm of how Atlantic Niño interacts with the upper‐ocean currents and intraseasonal waves remains elusive. Here we show that the equatorial Kelvin wave associated with the onset of the 2021 Atlantic Niño modulated both the background flow and the eddy flux of the equatorial upper‐ocean circulation, causing an extremely weak and delayed tropical instability wave (TIW) season. TIW‐induced variations of sea surface temperature (SST), sea surface salinity, sea surface height, and eddy temperature advection were exceptionally weak during May to July, the climatological peak of TIW activity, but rebounded in August when higher than normal variability was observed. Moored velocity data at 23°W show that during the peak of the 2021 Atlantic Niño from June to August, the Equatorial Undercurrent was deeper and stronger than usual. An anomalously weak eddy momentum flux strongly suppressed barotropic energy conversion north of the equator from May to July, likely contributing to low TIW activity. Reduced baroclinic energy conversion also might have played a role, as the meridional gradient of SST was sharply reduced during the Atlantic Niño. Despite extremely weak TIW velocities, modest intraseasonal variability of chlorophyll‐a (Chl‐ a ) was observed during the Atlantic Niño, due to pronounced meridional Chl‐ a gradients that partly compensated for the weak TIWs. Plain Language Summary Every few years the eastern equatorial Atlantic Ocean is significantly warmer than usual during boreal summer. Such warm events are referred to as Atlantic Niño events, and share similarities with El Niño events in the Pacific. In 2021, the strongest Atlantic Niño in at least four decades was observed in the equatorial Atlantic. This study is the first that investigates the complex interaction between Atlantic Niño, tropical Atlantic upper‐ocean currents, and equatorial waves based on various observational data sets. We show that the developing 2021 Atlantic Niño weakened both the background flow and the variability of near‐surface currents in May, which in turn largely reduced the strength of intraseasonal (20–50 days) waves that are usually generated by instability of the upper‐ocean zonal currents. As a consequence, the cooling effect that these waves usually have north of the equator and the warming effect along the equator vanished from May to July 2021. Interestingly, variability of chlorophyll concentration was enhanced, suggesting that enhanced meridional chlorophyll gradients compensated for reduced wave activity. Key Points The developing 2021 Atlantic Niño led to weaker equatorial surface currents and reduced vertical shear of upper‐ocean horizontal velocity Strong reduction of the surface flow, eddy flux, and meridional temperature gradient in May caused extremely weak and delayed tropical instability wave (TIW) season Reduced meridional TIW advection contributed to sharpen the north equatorial Chl‐ a front resulting in modest intraseasonal Chl‐ a variability

Description: Phytoplankton primary productivity (PP) varies significantly over environmental gradients, particularly in physically‐dynamic systems such as estuaries and coastal seas. During summer, runoff peaks in the Changjiang River driving large environmental gradients in both the Changjiang estuary and adjacent East China Sea (ECS), likely driving significant variability in PP. As satellite models of PP often underperform in coastal waters, we aimed to develop a novel approach for assessing net PP variability in such a dynamic environment. Parallel in situ measurements of Fast Repetition Rate (FRR) fluorometry and carbon (C) uptake rates were conducted for the first time in this region during two summer cruises in 2019 and 2021. A series of 13 C‐incubations ( n = 31) were performed, with measured PP ranging from ∼6 to 1,700 mgC m −3 d −1 . Net PP values were significantly correlated with salinity ( r = 0.45), phytoplankton chlorophyll a (Chl‐ a , r = 0.88), Photosystem II (PSII) functional absorption cross‐section ( σ PSII , r = −0.76) and maximum PSII quantum yield ( F v / F m , r = 0.59). Stepwise regression analysis showed that Chl‐ a and σ PSII were the strongest predictors of net PP. A generalized additive model (GAM) was also used to estimate net PP considering nonlinear effects of Chl‐ a and σ PSII . We demonstrate that GAM outperforms linear modeling approaches in estimating net PP in this study, as evidenced by a lower root mean square error (∼140 vs. 250 mgC m −3 d −1 ). Our novel approach provides a valuable tool to examine carbon cycling dynamics in this important region. Plain Language Summary The East China Sea has a complex current system that creates a highly dynamic physical environment for phytoplankton, particularly during the summer months. Net primary productivity (PP) is highly variable in this region, yet characterizing these spatial patterns in PP is difficult due to the lack of a high‐resolution data collecting method. Therefore, a strong need exists for a quick and easily implemented method for monitoring PP in this dynamic system. Based on parallel measurements of phytoplankton biomass and photophysiology, we present a novel approach that allows us to rapidly and easily assess regional PP at a high resolution. The high data volume potentially afforded by our net PP estimation method could not only contribute to a better understanding of PP variations in such a dynamic environment, but also help fill the large gaps in field data needed for validating satellite‐based PP models. Key Points Parallel in situ measurements of net primary productivity (PP) and Fast Repetition Rate fluorometry were conducted in the Changjiang estuary Productivity was highest at stations with high Chl and low σ PSII , typically located along the Chiangjiang river plume front A generalized additive model was developed to estimate net PP, providing an approach for assessing regional C‐cycling dynamics

Description: Hydrothermal fluids in ultramafic‐hosted hydrothermal systems associated with oceanic detachment faults can be more oxidizing compared to mafic‐hosted vent sites. These fluids form a mineral assemblage of pyrite, magnetite and hematite. At 13°30′N on the Mid‐Atlantic Ridge, chlorite‐quartz breccias recovered from an exposed fault scarp contain pyrite, with abundant magnetite and hematite, indicating that the redox of the fluids was variable. In primary micron‐scale zonations in pyrite, Ni, Co, and Se have a decoupled relationship, recording fluctuations in the chemical composition and temperature of hydrothermal fluid as the grains grew. Secondary zonations that erase and overprint primary zonations are limited to the grain margin and permeable regions within the grain core. Secondary zonations formed via two processes: (a) grain dissolution followed by overgrowth, and (b) remobilization of metals during oxidizing fluid flow events. In both instances, Ni and Co have been mobilized and concentrated, and are not lost to the hydrothermal fluid. Superimposed on these features is evidence of grain scale deformation related to periods of fault movement along the detachment surface. Sulfur isotope ratios (δ 34 S) in pyrite systematically decrease from the grain margin to the grain core, indicating that increased amounts of sulfur were derived from thermochemical sulfate reduction of seawater. Thus, pyrite records the evolution of fluid flow and deformation events during exhumation along the detachment surface from ∼1 to 2 km below the seafloor at the base of the lava pile, with temporal fluctuations in fluid redox identified as an important process in controlling Ni and Co enrichment in pyrite. Plain Language Summary Detachment faults are long lived faults that can expose ultramafic rocks at the seafloor. We aim to investigate the links between hydrothermal activity and detachment fault formation. To do this we use pyrite as a tape recorder for past fluid flow events. Across individual mineral grains, distinct zonations in metal content and sulfur isotope ratios show that the incursion of seawater occurred periodically during pyrite growth, increasing during fault movement events that lead to changes in the temperature and pH of the fluids in the fault zone. These changes concentrated metals toward the center of individual mineral grains. Zonations were then overprinted by later deformation‐related events, providing evidence that the samples formed at deeper crustal levels below the seafloor and were progressively exhumed at the seafloor over time. Key Points Microtextural, geochemical, and isotopic variations in subseafloor pyrite record the history of sample exhumation along a detachment fault Nickel and Co are remobilized and concentrated in pyrite across individual mineral grains in response to fluctuating fluid redox conditions Evidence of pyrite deformation and alteration mineralogy of samples indicates sample exhumation from a depth of 1–2 km

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: 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: 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: 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: “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: 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: 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 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: 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: 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: 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: 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: 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: 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: 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: Various nutrient sources in the upper waters of oceanic subtropical gyres, which are the Earth's largest oligotrophic ecosystems, play a crucial role in governing the sequestration of atmospheric CO2

Description: Key Points: - A new CHBr3 emission inventory based on natural and anthropogenic sources suggests that the latter account for 12%–28% of the global emissions - In the NH, new anthropogenic estimates increase known natural CHBr3 emissions by up to 70.5%, leading to higher atmospheric CHBr3 levels - At the NH extratropical tropopause, CHBr3 is enhanced by 0.9 ppt Br due to anthropogenic sources thus doubling natural CHBr3 abundances Bromoform (CHBr3) contributes to stratospheric ozone depletion but is not regulated under the Montreal Protocol due to its short lifetime and large natural sources. Here, we show that anthropogenic sources contribute significantly to the amount of CHBr3 transported into the Northern Hemisphere (NH) extratropical stratosphere. We present a new CHBr3 emission inventory comprised of natural and anthropogenic sources, with the latter estimated from ship ballast, power plant cooling and desalination plant brine water. Including anthropogenic sources in the new inventory increases CHBr3 emissions by up to 31.5% globally and 70.5% in the NH. In consequence, atmospheric CHBr3 is also significantly higher, especially over the NH extratropics during boreal winter. Here anthropogenic sources enhance bromine at the tropopause by 0.9 ppt Br, thus doubling natural CHBr3 abundances. For some latitudes, tropopause bromine increases by 2.4 ppt Br suggesting significant contributions of anthropogenic CHBr3 to the NH lowermost stratosphere.

Description: Valdivia Bank (VB) is a Late Cretaceous oceanic plateau formed by volcanism from the Tristan-Gough hotspot at the Mid-Atlantic Ridge (MAR). To better understand its origin and evolution, magnetic data were used to generate a magnetic anomaly grid, which was inverted to determine crustal magnetization. The magnetization model reveals quasi-linear polarity zones crossing the plateau and following expected MAR paleo-locations, implying formation by seafloor spreading over ∼4 Myr during the formation of anomalies C34n-C33r. Paleomagnetism and biostratigraphy data from International Ocean Discovery Program Expedition 391 confirm the magnetic interpretation. Anomaly C33r is split into two negative bands, likely by a westward ridge jump. One of these negative anomalies coincides with deep rift valleys, indicating their age and mechanism of formation. These findings imply that VB originated by seafloor spreading-type volcanism during a plate reorganization, not from a vertical stack of lava flows as expected for a large volcano. Key Points - Valdivia Bank is characterized by quasi-linear magnetic anomalies that are parallel to the inferred paleo-Mid-Atlantic Ridge - Magnetic anomalies imply that the plateau becomes younger E-W consistent with formation via seafloor spreading during anomalies C34n-C33r - Rift valleys, division of C33r, and anomaly curvature imply complex ridge tectonics and a ridge jump

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

Description: The Radiative Forcing Model Intercomparison Project (RFMIP) allows estimates of effective radiative forcing (ERF) in the Coupled Model Intercomparison Project phase six (CMIP6). We analyze the RFMIP output, including the new experiments from models that use the same parameterization for anthropogenic aerosols (RFMIP-SpAer), to characterize and better understand model differences in aerosol ERF. We find little changes in the aerosol ERF for 1970–2014 in the CMIP6 multi-model mean, which implies greenhouse gases primarily explain the positive trend in the total anthropogenic ERF. Cloud-mediated effects dominate the present-day aerosol ERF in most models. The results highlight a regional increase in marine cloudiness due to aerosols, despite suppressed cloud lifetime effects in that RFMIP-SpAer experiment. Negative cloud-mediated effects mask positive direct effects in many models, which arise from strong anthropogenic aerosol absorption. The findings suggest opportunities to better constrain simulated ERF by revisiting the optical properties and long-range transport of aerosols. Key Points: - Coupled Model Intercomparison Project phase six (CMIP6) averaged trend in aerosol effective radiative forcing (ERF) is small for 1970–2014 and weakly positive for 2000–2014 - Positive direct aerosol radiative effects in CMIP6 models are associated with strong aerosol absorption - Diverse and often strong cloud-mediated effects primarily determine the magnitude of aerosol ERF in CMIP6

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: The oceanic uptake and resulting storage of the anthropogenic CO2 (Cant) that humans have emitted into the atmosphere moderates climate change. Yet our knowledge about how this uptake and storage has progressed in time remained limited. Here, we determine decadal trends in the storage of Cant by applying the eMLR(C*) regression method to ocean interior observations collected repeatedly since the 1990s. We find that the global ocean storage of Cant grew from 1994 to 2004 by 29 ± 3 Pg C dec−1 and from 2004 to 2014 by 27 ± 3 Pg C dec−1 (±1σ). The storage change in the second decade is about 15 ± 11% lower than one would expect from the first decade and assuming proportional increase with atmospheric CO2. We attribute this reduction in sensitivity to a decrease of the ocean buffer capacity and changes in ocean circulation. In the Atlantic Ocean, the maximum storage rate shifted from the Northern to the Southern Hemisphere, plausibly caused by a weaker formation rate of North Atlantic Deep Waters and an intensified ventilation of mode and intermediate waters in the Southern Hemisphere. Our estimates of the Cant accumulation differ from cumulative net air-sea flux estimates by several Pg C dec−1, suggesting a substantial and variable, but uncertain net loss of natural carbon from the ocean. Our findings indicate a considerable vulnerability of the ocean carbon sink to climate variability and change. Key Points: - The global ocean storage of anthropogenic carbon grew by 29 ± 3 and 27 ± 3 Pg C dec−1 from 1994 to 2004 and 2004 to 2014, respectively - The change in oceanic storage of anthropogenic carbon relative to the atmospheric CO2 growth decreased by 15 ± 11% from the first to the second decade - This reduction is attributed to a decrease of the ocean buffer capacity and changes in ocean circulation

Description: Marine nitrogen (N2) fixation supports significant primary productivity in the global ocean. However, in one of the most productive regions of the world ocean, the northern Humboldt Upwelling System (HUS), the magnitude and spatial distribution of this process remains poorly characterized. This study presents a spatially resolved dataset of N2 fixation rates across six coastal transects of the northern HUS off Peru (8°S – 16°S) during austral summer. N2 fixation rates were detected throughout the waters column including within the OMZ between 12°S and 16°S. N2 fixation rates were highest where the subsurface Oxygen Minimum Zone (OMZ, O2 〈20 µmol L-1) was most intense and estimated nitrogen (N) loss was highest. There, rates were measured throughout the water column. Hence the vertical and spatial distribution of rates indicates colocation of N2 fixation with N loss in the coastal productive waters of the northern HUS. Despite high phosphate and total dissolvable iron (TdFe) concentrations throughout the study area, N2 fixation was still generally low (1.19 ± 3.81 nmol L-1 d-1) and its distribution could not be directly explained by these two factors. Our results suggest that the distribution was likely influenced by a complex interplay of environmental factors including phytoplankton biomass and organic matter availability, and potentially iron, or other trace metal (co)-limitation of both N2 fixation and primary production. In general, our results support previous conclusions that N2 fixation in the northern HUS plays a minor role as a source of new N and to replenish the regional N loss. Key Points: A north-to-south pattern in N2 fixation rates was observed implying increased N turnover between 12°S and 16°S where N loss was pronounced Highest N2 fixation rates were measured in coastal productive waters above and within the OMZ, showing no clear relationship with Fe or P The magnitude of N2 fixation was low compared to predictions, estimated to account for ∼0.3% of primary production and 〈2% of local N loss

Description: We introduce a new hypothesis concerning the role of internal climate dynamics in the non-linear transitions from interglacial to glacial (IG-G) state since the Mid Pleistocene Transition (MPT). These transitions encompass large and abrupt changes in atmospheric CO2, ice volume, and temperature that we suggest involve critical interactions between insolation and high amplitude oscillations in ocean/atmosphere circulation patterns. Specifically, we highlight the large amplitude of millennial-scale climate oscillations across the transition from Marine Isotope Stage (MIS) 5 to 4, which we argue led to amplified cooling of the deep ocean and we demonstrate that analogous episodes of extreme cooling systematically preceded glacial periods of the last 800 kyr. We suggest that such cooling necessitates a reconfiguration of the deep ocean to avoid a density paradox between northern and southern-sourced deep waters (SSW), which could be accomplished by increasing the relative volume and or salinity of SSW, thus providing the necessary storage capacity for the subsequent (delayed) and relatively abrupt drawdown of CO2. We therefore explain the transient decoupling of Antarctic temperature from CO2 across MIS 5/4 as a direct consequence of millennial activity at that time. We further show that similar climatic decoupling typically occurred during times of low obliquity and was a ubiquitous feature of IG-G transitions over the past 800 kyr, producing the appearance of bimodality in records of CO2, benthic δ18O and others. Finally we argue that the apparent lack of bimodality in the pre-MPT record of benthic δ18O implies that the dynamics associated with IG-G transitions changed across the MPT

Description: The Atlantic Meridional Overturning Circulation (AMOC) is a key feature of the North Atlantic with global ocean impacts. The AMOC's response to past changes in forcings during the Holocene provides important context for the coming centuries. Here, we investigate AMOC trends using an emerging set of transient simulations using multiple global climate models for the past 6,000 years. Although some models show changes, no consistent trend in overall AMOC strength during the mid-to-late Holocene emerges from the ensemble. We interpret this result to suggest no overall change in AMOC, which fits with our assessment of available proxy reconstructions. The decadal variability of the AMOC does not change in ensemble during the mid- and late-Holocene. There are interesting AMOC changes seen in the early Holocene, but their nature depends a lot on which inputs are used to drive the experiment.

Description: We calculate the depth to magnetic basement and the average crustal magnetic susceptibility, which is sensitive to the presence of iron-rich minerals, to interpret the present structure and the tecto-magmatic evolution in the Central Tethyan belt. Our results demonstrate exceptional variability of crustal magnetization with smooth, small-amplitude anomalies in the Gondwana realm and short-wavelength high-amplitude variations in the Laurentia realm. Poor correlation between known ophiolites and magnetization anomalies indicates that Tethyan ophiolites are relatively poorly magnetized, which we explain by demagnetization during recent magmatism. We analyze regional magnetic characteristics for mapping previously unknown oceanic fragments and mafic intrusions, hidden beneath sedimentary sequences or overprinted by tectono-magmatic events. By the style of crustal magnetization, we distinguish three types of basins and demonstrate that many small-size basins host large volumes of magmatic rocks within or below the sedimentary cover. We map the width of magmatic arcs to estimate paleo-subduction dip angle and find no systematic variation between the Neo-Tethys and Paleo-Tethys subduction systems, while the Pontides magmatic arc has shallow (∼15o) dip in the east and steep (∼50-55o) dip in the west. We recognize an unknown, buried 450 km-long magmatic arc along the western margin of the Kırşehir massif formed above steep (55o) subduction. We propose that lithosphere fragmentation associated with Neo-Tethys subduction systems may explain high-amplitude, high-gradient crustal magnetization in the Caucasus Large Igneous Province. Our results challenge conventional regional geological models, such as Neo-Tethyan subduction below the Greater Caucasus, and call for reevaluation of the regional paleotectonics. Key Points: Magnetic regionalization does not fully match regional geological models in the Central Tethyan Belt We identify previously unknown magmatic arcs and ocean relics Magnetization is weak in Gondwana and strong in Laurentia terranes: Kirsehir massif has Laurentia affinity

Description: Offshore meteoric groundwater (OMG) has long been hypothesised to be a driver of seafloor geomorphic processes in continental margins worldwide. Testing this hypothesis has been challenging because of our limited understanding of the distribution and rate of OMG flow and seepage, and their efficacy as erosive/destabilising agents. Here we carry out numerical simulations of groundwater flow and slope stability using conceptual models and evolving stratigraphy - for passive siliciclastic and carbonate margin cases – to assess whether OMG and its evolution during a late Quaternary glacial cycle can generate the pore pressures required to trigger mechanical instabilities on the seafloor. Conceptual model results show that mechanical instabilities by OMG flow are most likely to occur in the outer shelf to upper slope, at or shortly before the Last Glacial Maximum sea level lowstand. Models with evolving stratigraphy show that OMG flow is a key driver of pore pressure development and instability in the carbonate margin case. In the siliciclastic margin case, OMG flow plays a secondary role in preconditioning the slope to failure. The higher degree of spatial/stratigraphic heterogeneity of carbonate margins, lower shear strengths of their sediments, and limited generation of overpressures by sediment loading may explain the higher susceptibility of carbonate margins, in comparison to siliciclastic margins, to mechanical instability by OMG flow. OMG likely played a more significant role in carbonate margin geomorphology (e.g. Bahamas, Maldives) than currently thought. Key Points Offshore meteoric groundwater (OMG) flow can drive mechanical instabilities in the outer shelf to upper slope Such instabilities occur at, or shortly after, the Last Glacial Maximum sea level lowstand Carbonate margins are more susceptible to mechanical instability by OMG than siliciclastic margins

Description: The current narrative of artificial upwelling (AU) is to translocate nutrient rich deep water to the ocean surface, thereby stimulating the biological carbon pump (BCP). Our refined narrative takes the response of the solubility pump and the CO2 emission scenario into account. Using global ocean-atmosphere model experiments we show that the effectiveness of a hypothetical maximum AU deployment in all ocean areas where AU is predicted to lower surface pCO2, the draw down of CO2 from the atmosphere during years 2020–2100 depends strongly on the CO2 emission scenario and ranges from 1.01 Pg C/year (3.70 Pg CO2/year) under RCP 8.5 to 0.32 Pg C/year (1.17 Pg CO2/year) under RCP 2.6. The solubility pump becomes equally effective compared to the BCP under the highest emission scenario (RCP 8.5), but responds with CO2 outgassing under low CO2 emission scenarios. Key Points: - Artificial upwelling (AU) effectiveness to draw down CO2 from the atmosphere is strongly dependent on the future CO2 emission scenario - The solubility pump becomes as effective as the biological carbon pump under high emission scenarios - Organic matter transfer efficiency decreases under AU, likely due to higher water temperatures below the ocean's surface

Description: An insufficient supply of the micronutrient iron (Fe) limits phytoplankton growth across large parts of the ocean. Ambient Fe speciation and solubility are largely dependent on seawater physico-chemical properties. We calculated the apparent Fe solubility (SFe(III)app) at equilibrium for ambient conditions, where SFe(III)app is defined as the sum of aqueous inorganic Fe(III) species and Fe(III) bound to organic matter formed at a free Fe3+ concentration equal to the solubility of Fe hydroxide. We compared the SFe(III)app to measured dissolved Fe (dFe) in the Atlantic and Pacific Oceans. The SFe(III)app was overall ∼2 to 4-fold higher than observed dFe at depths less than 1000 m, ∼2-fold higher than the dFe between 1000-4000 m and ∼3-fold higher than dFe below 4000 m. Within the range of used parameters, our results showed that there was a similar trend in the vertical distributions of horizontally averaged SFe(III)app and dFe. Our results suggest that vertical dFe distributions are underpinned by changes in SFe(III)app which are driven by relative changes in ambient pH and temperature. Since both pH and temperature are essential parameters controlling ambient Fe speciation, these should be accounted for in investigations of changing Fe dynamics, particularly in the context of ocean acidification and warming. Key Points Apparent iron solubility is driven by ambient pH, temperature (T) and dissolved organic carbon (DOC), and showed a 6-fold variation between surface (pH= 8.05 on the total scale, DOC= 71.8 µmol L-1, T= 20.4 °C) and deep oceanic waters (pH= 7.82, DOC= 38.6 µmol L-1, T= 1.1°C). Higher values of apparent iron solubility were determined for deep Atlantic and Pacific waters, with lower values in subtropical gyres. Calculated apparent iron solubility showed a similar trend in vertical distribution to dissolved iron, highlighting the importance of considering the impact of changes in ambient physico-chemical conditions on seawater iron chemistry.

Description: The correlation between concentrations of dissolved barium (dBa) and silicon (dSi) in the modern ocean supports the use of Ba as a paleoceanographic proxy. However, the mechanisms behind their linkage and the exact processes controlling oceanic Ba cycling remain enigmatic. To discern the extent to which this association arises from biogeochemical processes versus physical mixing, we examine the behavior of Ba and Si at the Congo River-dominated Southeast Atlantic margin where active biological processes and large boundary inputs override the large-scale ocean circulation. Here we present the first combined measurements of dissolved stable Ba (δ138Ba) and Si (δ30Si) isotopes as well as Ba and Si fluxes estimated based on 228Ra from the Congo River mouth to the northern Angola Basin. In the surface waters, river-borne particle desorption or dissolution and shelf inputs lead to non-conservative additions of both dBa and dSi to the Congo-shelf-zone, with the Ba flux increasing more strongly than that of Si across the shelf. In the epipelagic and mesopelagic layers, Ba and Si are decoupled likely due to different depths of in situ barite precipitation and biogenic silica production. In the deep waters of the northern Angola Basin, we observe large enrichment of dBa, likely originating from high benthic inputs from the Congo deep-sea fan sediments. Our results reveal different mechanisms controlling the biogeochemical cycling of Ba and Si and highlight a strong margin influence on marine Ba cycling. Their close association across the global ocean must therefore mainly be a consequence of the large-scale ocean circulation. Key Points Stronger enrichment of dissolved barium (dBa) than silicon (dSi) observed in the shelf-zone of the Congo plume Diatom silica production has negligible effect on dissolved Ba isotopic compositions in large river plumes Strong dBa enrichment (up to 24 nM) in the deep water of the northern Angola Basin likely originates from high benthic inputs

Description: This contribution to the RECCAP2 (REgional Carbon Cycle Assessment and Processes) assessment analyzes the processes that determine the global ocean carbon sink, and its trends and variability over the period 1985-2018, using a combination of models and observation-based products. The mean sea-air CO2 flux from 1985 to 2018 is -1.6 +/- 0.2 PgC yr(-1) based on an ensemble of reconstructions of the history of sea surface pCO(2) (pCO(2) products). Models indicate that the dominant component of this flux is the net oceanic uptake of anthropogenic CO2, which is estimated at -2.1 +/- 0.3 PgC yr(-1) by an ensemble of ocean biogeochemical models, and -2.4 +/- 0.1 PgC yr(-1) by two ocean circulation inverse models. The ocean also degasses about 0.65 +/- 0.3 PgC yr(-1) of terrestrially derived CO2, but this process is not fully resolved by any of the models used here. From 2001 to 2018, the pCO2 products reconstruct a trend in the ocean carbon sink of -0.61 +/- 0.12 PgC yr(-1) decade(-1), while biogeochemical models and inverse models diagnose an anthropogenic CO2-driven trend of -0.34 +/- 0.06 and -0.41 +/- 0.03 PgC yr(-1) decade(-1), respectively. This implies a climate-forced acceleration of the ocean carbon sink in recent decades, but there are still large uncertainties on the magnitude and cause of this trend. The interannual to decadal variability of the global carbon sink is mainly driven by climate variability, with the climate-driven variability exceeding the CO2-forced variability by 2-3 times. These results suggest that anthropogenic CO2 dominates the ocean CO2 sink, while climate-driven variability is potentially large but highly uncertain and not consistently captured across different methods.

Description: Geochemical analyses of trace elements in the ocean water column have suggested that pelagic clay‐rich sediments are a major source of various elements to bottom‐waters. However, corresponding high‐quality measurements of trace element concentrations in porewaters of pelagic clay‐rich sediments are scarce, making it difficult to evaluate the contributions from benthic processes to global oceanic cycles of trace elements. To bridge this gap, we analyzed porewater and bulk sediment concentrations of vanadium, chromium, cobalt, nickel, copper, arsenic, molybdenum, barium and uranium, as well as concentrations of the major oxidants nitrate, manganese, iron, and sulfate in the top 30 cm of cores collected along a transect from Hawaii to Alaska. The data show large increases in porewater concentrations of vanadium, manganese, cobalt, nickel, copper, and arsenic within the top cm of the sediment, consistent with the release of these elements from remineralized organic matter. The sediments are a sink for sulfate, uranium, and molybdenum, even though conditions within the sampled top 30 cm remain aerobic. Porewater chromium concentrations generally increase with depth due to release from sediment particles. Extrapolated to the global aerial extent of pelagic clay sediment, the benthic fluxes in mol yr −1 are Ba 3.9 ± 3.6 × 10 9 , Mn 3.4 ± 3.5 × 10 8 , Co 2.6 ± 1.3 × 10 7 , Ni 9.6 ± 8.6 × 10 8 , Cu 4.6 ± 2.4 × 10 9 , Cr 1.7 ± 1.1 × 10 8 , As 6.1 ± 7.0 × 10 8 , V 6.0 ± 2.5 × 10 9 . With the exception of vanadium, calculated fluxes across the sediment–water interface are consistent with the variability in bottom‐water concentrations and ocean residence time of the studied elements.

Description: The additional water from the Antarctic ice sheet and ice shelves due to climate‐induced melt can impact ocean circulation and global climate. However, the major processes driving melt are not adequately represented in Coupled Model Intercomparison Project phase 6 (CMIP6) models. Here, we analyze a novel multi‐model ensemble of CMIP6 models with consistent meltwater addition to examine the robustness of the modeled response to meltwater, which has not been possible in previous single‐model studies. Antarctic meltwater addition induces a substantial weakening of open‐ocean deep convection. Additionally, Antarctic Bottom Water warms, its volume contracts, and the sea surface cools. However, the magnitude of the reduction varies greatly across models, with differing anomalies correlated with their respective mean‐state climatology, indicating the state‐dependency of the climate response to meltwater. A better representation of the Southern Ocean mean state is necessary for narrowing the inter‐model spread of response to Antarctic meltwater. Plain Language Summary The melting of the Antarctic ice sheet and ice shelves can have significant impacts on ocean circulation and thermal structure, but current climate models do not fully capture these effects. In this study, we analyze seven climate models to understand how they respond to the addition of meltwater from Antarctica. We find that the presence of Antarctic meltwater leads to a significant weakening of deep convection in the open ocean. The meltwater also causes Antarctic Bottom Water to warm and its volume to decrease, while the sea surface cools and sea ice expands. However, the magnitude of the response to meltwater varies across models, suggesting that the mean‐state conditions of the Southern Ocean play a role. A better representation of the mean state and the inclusion of Antarctic meltwater in climate models will help reduce uncertainties and improve our understanding of the impact of Antarctic meltwater on climate. Key Points Antarctic meltwater substantially reduces the strength of simulated Southern Ocean deep convection in climate models The additional meltwater induces Antarctic Bottom Water warming and contraction, with dense water classes converting to lighter ones Differences in the magnitude of these responses between models can be partly attributed to their different base states

Description: The Java ‐ Lesser Sunda margin, which features multi‐scale subducting oceanic basement relief, is classified as neutral (Lombok and Sumbawa) to erosional (Central Java to Bali) in comparison to its accretionary counterpart offshore Sumatra. However, a comprehensive analysis of how plate boundary and upper plate structure across the neutral to erosional transition are modulated by the subduction of oceanic basement relief is lacking to date. To shed light on the tectonic parameters that push the margin into the neutral or erosional domain, we combine multi‐channel reflection seismic images derived through a grid‐based P‐wave velocity inversion, and multibeam bathymetric maps. The data document how different scales of subducting topography modify seafloor morphology, upper plate structure, and décollement position. Large‐scale subducting features cause a landward shift of the deformation front, shortening of the accretionary wedge, and seafloor steepening at the relief's trailing edge. Small‐scale subducting ridges primarily impact the frontal prism resulting in over‐steepening at the trench and localized slope failure. Ahead of subducting relief, deformation of the accretionary wedge encompasses enhanced compression and a reduction in seafloor slope but appears independent of the size of the relief. Ridge and seamount subduction induce frontal erosion and basal erosion offshore Lombok and Bali, respectively. Our P‐wave velocity models indicate that the rigidity of the upper plate's base along the eastern Sunda margin is significantly lower than the worldwide trend. We conclude that this favors the genesis of tsunami earthquakes that have occurred on the Java margin. Plain Language Summary The convergence of the tectonic plates drives a wide variety of geological processes along the plate margins, including the formation of the forearc accretionary wedge, volcanic activities, and megathrust earthquakes. Over the past 40 years, marine research shows that different sizes of oceanic reliefs (seamounts and ridges) are widely distributed over the seafloor, approaching the trench, and eventually subducted underneath the overriding plate. An accurate observation of the subducted reliefs and their tectonic impact on the overriding plate depends on different observation approaches, data processing methods, and the evolutionary history of the forearc. In the Java margin, the oceanic seafloor features massive seamounts with different scales and shapes, and the bathymetry of the overriding plate is highly disturbed. This provides us with the best opportunity of studying the rugged seafloor's seismogenic and geological impacts. By using state‐of‐the‐art seismic imaging techniques, we image the subsurface structures, obtain the forearc velocity, identify the seamounts, and discuss the seamounts' effect on structural deformation and megathrust earthquake occurrence. Distinctively, the marine forearc gets shortened and thickened significantly by seamount subduction. Structural images sharply illustrate different deformation patterns and stress regimes at the seamount's different flanks and reveal the possible process of subduction erosion. Key Points Upper plate deformation scales with variable subducting relief, as observed along the Java Trench in seismic sections and bathymetry Subduction of seafloor topography induces progression from an accretion‐dominated domain toward a phase of subduction erosion The overall low rigidity of the upper plate's base may contribute to the Java margin earthquake's tsunami‐genesis

Description: Accurately predicting future ocean acidification (OA) conditions is crucial for advancing OA research at regional and global scales, and guiding society's mitigation and adaptation efforts. This study presents a new model-data fusion product covering 10 global surface OA indicators based on 14 Earth System Models (ESMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6), along with three recent observational ocean carbon data products. The indicators include fugacity of carbon dioxide, pH on total scale, total hydrogen ion content, free hydrogen ion content, carbonate ion content, aragonite saturation state, calcite saturation state, Revelle Factor, total dissolved inorganic carbon content, and total alkalinity content. The evolution of these OA indicators is presented on a global surface ocean 1° × 1° grid as decadal averages every 10 years from preindustrial conditions (1750), through historical conditions (1850–2010), and to five future Shared Socioeconomic Pathways (2020–2100): SSP1-1.9, SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5. These OA trajectories represent an improvement over previous OA data products with respect to data quantity, spatial and temporal coverage, diversity of the underlying data and model simulations, and the provided SSPs. The generated data product offers a state-of-the-art research and management tool for the 21st century under the combined stressors of global climate change and ocean acidification. The gridded data product is available in NetCDF at the National Oceanic and Atmospheric Administration (NOAA) National Centers for Environmental Information: https://www.ncei.noaa.gov/data/oceans/ncei/ocads/metadata/0259391.html, and global maps of these indicators are available in jpeg at: https://www.ncei.noaa.gov/access/ocean-carbon-acidification-data-system/synthesis/surface-oa-indicators.html. Key Points: - This study presents the evolution of 10 ocean acidification (OA) indicators in the global surface ocean from 1750 to 2100 - By leveraging 14 Earth System Models (ESMs) and the latest observational data, it represents a significant advancement in OA projections - This inter-model comparison effort showcases the overall agreements among different ESMs in projecting surface ocean carbon variables

Description: We evaluate the decadal evolution of ventilation and anthropogenic carbon (C-ant) in the Nordic Seas between 1982 and the 2010s. Ventilation changes on decadal timescale are identified by evaluating decadal changes in mean ages and apparent oxygen utilization in each of the four main basins of the Nordic Seas (the Greenland and Iceland Seas, and the Norwegian and Lofoten Basins). The ages are derived from the transient time distribution approach, based on the transient tracers chlorofluorocarbon-12 (CFC-12) and sulfur hexafluoride (SF6). The different decades show different phases in ventilation, with the 2000s being overall better ventilated than the 1990s in all basins. For the Greenland Sea, we also show that the 2010s are better ventilated than the 2000s, with a clear shift in hydrographic properties. The evolution of concentrations and inventory of C-ant is linked to the ventilation state. The deep waters get progressively older over the analyzed period, which is connected to the increased fraction of deep water from the Arctic Ocean.Plain Language Summary The ocean region between Greenland, Iceland, and Norway, called the Nordic Seas, is a main site of deep-water formation. This process produces dense waters and brings surface waters to larger depths, thereby ventilating the water below. This transports, among other things, man-made CO2 (anthropogenic carbon; C-ant) and oxygen from the atmosphere into the interior ocean, thereby reducing the amount of CO2 stored in the atmosphere. This study investigates how the ventilation has changed in the Nordic Seas from 1982 to the 2010s. We find that the ventilation has changed with time, from a rather well-ventilated state in 1982, to a reduced ventilation in the 1990s, and then a restrengthened ventilation from the 2000s.

Description: We present the first observational evidence for convectively generated cold pools (CP) as driving mechanism for low-level jets (LLJ). Our findings are based on a unique campaign data set that allowed us to perform a systematic assessment of the process. During the three-month campaign in Germany, 6.8% of all identified LLJ profiles were connected to a CP (CPLLJ). Most measured CPLLJs appeared with the CP front and lasted for up to two hours. Moreover, we have observed a CP favoring the formation of a several-hours long LLJ. In that case, a strong LLJ and cooling of the atmosphere between the surface and at least 400 m a.g.l. were seen when the density current reached the measurement site. The development led to the formation of a near-surface temperature inversion during daytime as a prerequisite for the LLJ, not unlike the mechanism of nocturnal LLJs. Key Points: - Observed low-level jets connected to convective cold pools were about 7% of all jet profiles during summer campaign in Germany - Convective cold pools favored reduced frictional coupling of the wind field as a prerequisite for generating low-level jets during daytime - Low-level jets connected to convective cold pools were on average weaker but gustier than nocturnal jets

Description: The development of stable barium (Ba) isotope measurements provides a novel tool to investigate the geochemical cycling of Ba in the ocean and its sediments. In sediment pore waters, gradients of dissolved Ba concentrations result from various diagenetic processes. The distribution and fractionation of Ba isotopes in the pore waters are expected to further improve our understanding of these early diagenetic control mechanisms. Here, we present pore water profiles of dissolved stable Ba isotopic signatures (δ138Bapw) from shallow water sediments covering the entire Pearl River Estuary (PRE) in Southern China. We find pronounced depth-dependent Ba isotope variations generally showing a shift from heavy to light δ138Bapw signatures from the sediment surface down to 15 cm depth. These gradients are well reproduced by a diffusion-reaction model, which generates an apparent fractionation factor (138ε) of −0.60 ± 0.10‰ pointing to preferential removal of low-mass Ba isotopes from the pore water during solution-solid phase interactions. Consequently, the combined diagenetic processes induce the highest δ138Bapw values of +0.5 to +0.7‰ in the pore waters of the topmost sediment layer. Although the detrital fraction dominates the Ba content in the PRE surface sediments, the determined gradients of pore water Ba isotopes, together with concentration variations of Ba and other redox-sensitive elements such as manganese (Mn), show that non-detrital excess Ba carriers including Mn oxides and authigenic barite clearly affect the post-depositional Ba dynamics. Stable Ba isotopes are thus a potentially powerful tracer of Ba geochemistry during early sediment diagenesis in estuarine depositional environments. Key Points We present a data set of dissolved stable Ba isotopic compositions in surface sediment pore waters of a large river estuary Pore water Ba isotope values generally decrease with increasing sediment depth, reflecting post-depositional Ba isotope fractionation A diffusion-reaction model predicts the distribution and fractionation of stable Ba isotopes in the sediment pore waters well

Description: Subtropical gyres cover 26%-29% of the world's surface ocean and are conventionally regarded as ocean deserts due to their permanent stratification, depleted surface nutrients, and low biological productivity. Despite tremendous advances over the past three decades, particularly through the Hawaii Ocean Time-series and the Bermuda Atlantic Time-series Study, which have revolutionized our understanding of the biogeochemistry in oligotrophic marine ecosystems, the gyres remain understudied. We review current understanding of upper ocean biogeochemistry in the North Pacific Subtropical Gyre, considering other subtropical gyres for comparison. We focus our synthesis on spatial variability, which shows larger than expected dynamic ranges of properties such as nutrient concentrations, rates of N-2 fixation, and biological production. This review provides new insights into how nutrient sources drive community structure and export in upper subtropical gyres. We examine the euphotic zone (EZ) in subtropical gyres as a two-layered vertically structured system: a nutrient-depleted layer above the top of the nutricline in the well-lit upper ocean and a nutrient-replete layer below in the dimly lit waters. These layers vary in nutrient supply and stoichiometries and physical forcing, promoting differences in community structure and food webs, with direct impacts on the magnitude and composition of export production. We evaluate long-term variations in key biogeochemical parameters in both of these EZ layers. Finally, we identify major knowledge gaps and research challenges in these vast and unique systems that offer opportunities for future studies. Key Points Subtropical gyres display larger spatiotemporal dynamics in biogeochemical properties than previously considered An improved two-layer framework is proposed for the study of nutrient-driven and biologically mediated carbon export in the euphotic zone Future research will benefit from high-resolution samplings, improved sensitivity of nutrient analyses, and advanced modeling capabilities

Description: Because new observations have revealed that the Labrador Sea is not the primary source for waters in the lower limb of the Atlantic Meridional Overturning Circulation (AMOC) during the Overturning in the Subpolar North Atlantic Programme (OSNAP) period, it seems timely to re-examine the traditional interpretation of pathways and property variability for the AMOC lower limb from the subpolar gyre to 26.5 degrees N. In order to better understand these connections, Lagrangian experiments were conducted within an eddy-rich ocean model to track upper North Atlantic Deep Water (uNADW), defined by density, between the OSNAP line and 26.5 degrees N as well as within the Labrador Sea. The experiments reveal that 77% of uNADW at 26.5 degrees N is directly advected from the OSNAP West section along the boundary current and interior pathways west of the Mid-Atlantic Ridge. More precisely, the Labrador Sea is a main gateway for uNADW sourced from the Irminger Sea, while particles connecting OSNAP East to 26.5 degrees N are exclusively advected from the Iceland Basin and Rockall Trough along the eastern flank of the Mid-Atlantic Ridge. Although the pathways between OSNAP West and 26.5 degrees N are only associated with a net formation of 1.1 Sv into the uNADW layer, they show large density changes within the layer. Similarly, as the particles transit through the Labrador Sea, they undergo substantial freshening and cooling that contributes to further densification within the uNADW layer. Key Points: - The large majority of upper North Atlantic Deep Water (uNADW) sourced from the Irminger Sea transits through the Labrador Sea before reaching 26.5°N - Interior pathways along the eastern flank of the Mid-Atlantic Ridge connect the Iceland Basin and Rockall Trough to 26.5°N - Though uNADW is mainly sourced in the eastern subpolar gyre, its transit in the Labrador Sea is associated with further property changes

Description: Mesoscale eddies are frequently observed in the Eastern Tropical North Atlantic (ETNA), yet their effects on the transport and distribution of biogeochemical solutes, and specifically on the production and remineralization of dissolved organic matter (DOM) remain difficult to elucidate. Here, we investigated the submesoscale variability of chromophoric DOM (CDOM) and fluorescent DOM (FDOM) together with microbial production and remineralization processes in two cyclonic eddies (CEs) in the ETNA during summer and winter 2019. One CE, formed near the coast off Mauritania during the post-upwelling season, was sampled along a ∼900 km zonal corridor between Mauritania and the Cape Verde Islands. The other CE, formed nearby Brava Island, was out of coastal influence. Four fluorescent components were identified with parallel factor analysis, two humic-like, and two protein-like components. Humic-like FDOM components correlated to optode-based community respiration and were also good indicators of upwelling associated with the Brava Island CE as they correlated to physical parameters (e.g., temperature) and to dissolved inorganic nitrogen. The tryptophan-like FDOM components correlated with the carbon and nitrogen content of semi-labile DOM, phytoplankton biomass, community respiration, and bacterial production. Overall, our study revealed that DOM optical properties are suitable for tracing freshly produced organic matter and the transport of remineralized DOM within offshore eddies. Key Points: - Four fluorescent dissolved organic matter (FDOM) components were studied in two cyclonic eddies (CEs) in the Eastern Tropical North Atlantic - Tryptophan-like FDOM was an indicator of the CEs' productivity as it correlated with semi-labile dissolved organic matter and microbial metabolic activities - Humic-like FDOM was a by-product of microbial respiration; its distribution within an offshore CE covaried with nutrient upwelling

Description: Key Points: - Glacier-derived particles release 2–46% of labile particulate lead (Pb) upon mixing with seawater - Pb dynamics in glacier fjords are characterized by a rapid release of dissolved Pb followed by readsorption on a timescale of hours-to-days - Dissolved Pb release from the Greenland Ice Sheet is likely within the range 0.2–1 Mmol yr−1 Higher than expected concentrations of dissolved lead (dPb) have been consistently observed along glaciated coastlines and it is widely hypothesized that there is a net release of dPb from glacier-derived sediments. Here we further corroborate that dPb concentrations in diverse locations around west Greenland (3.2–252 pM) and the Western Antarctic Peninsula (7.7–107 pM) appear to be generally higher than can be explained by addition of dPb from glacier-derived freshwater. The distribution of dPb across the salinity gradient is unlike any other commonly studied trace element (e.g., Fe, Co, Ni, Cu, Mn, and Al) implying a dynamic, reversible exchange between dissolved and labile particulate Pb. Incubating a selection of glacier-derived particles from SW Greenland (Ameralik and Nuup Kangerlua) and Svalbard (Kongsfjorden), with a range of labile particulate Pb (LpPb) content (11–113 nmol g−1), the equivalent of 2–46% LpPb was released as dPb within 24 hr of addition to Atlantic seawater. Over longer time periods, the majority of this dPb was typically readsorbed. Sediment loading was the dominant factor influencing the net release of dPb into seawater, with a pronounced decline in net dPb release efficiency when sediment load increased from 20 to 500 mg L−1. Yet temperature also had some effect with 68 ± 22% higher dPb release at 11°C compared to 4°C. Future regional changes in dPb dynamics may therefore be more sensitive to short-term suspended sediment dynamics, and potentially temperature changes, than to changing interannual runoff volume.

Description: Mechanisms related to sub-seabed fluid flow processes are complex and inadequately understood. Petrophysical properties, availability of gases, topography, stress directions, and various geological parameters determine the location and intensity of leakage which change over time. From tens of seafloor pockmarks mapped along Vestnesa Ridge on the west-Svalbard margin, only six show persistent present-day seepage activity in sonar data. To investigate the causes of such restricted gas seepage, we conducted a study of anisotropy within the conduit feeding one of these active pockmarks (i.e., Lunde Pockmark). Lunde is ∼400–500 m in diameter, and atop a ∼300–400 m wide seismic chimney structure. We study seismic anisotropy using converted S-wave data from 22 ocean-bottom seismometers (OBSs) located in and around the pockmark. We investigate differences in symmetry plane directions in anisotropic media using null energy symmetries in transverse components. Subsurface stress distribution affects fault/fracture orientations and seismic anisotropy, and we use S-wave and high-resolution 3D seismic data to infer stress regimes in and around the active seep site and study the effect of stresses on seepage. We observe the occurrence of changes in dominant fault/fracture and horizontal stress orientations in and around Lunde Pockmark and conclude minimum (NE-SW) and maximum (SE-NW) horizontal stress directions. Our analysis indicates a potential correlation between hydrofractures and horizontal stresses, with up to a ∼32% higher probability of alignment of hydrofractures and faults perpendicular to the inferred minimum horizontal stress direction beneath the Lunde Pockmark area. Key Points The S-wave analysis using ocean-bottom seismic (OBS) data indicates seismic anisotropy around a seeping pockmark on the W-Svalbard Margin The occurrence and orientation of symmetry planes in shallow anisotropic sediments vary across the pockmark Combined analyses using S-wave and 3-D seismic data suggest that preferred fault and fracture orientations follow local stress conditions

Description: Key Points: - During 1993–2019, the East Greenland Coastal Current is freshest in 2010 and 2012 notably matching years of exceptional Greenland runoff - Freshwater anomalies from sea-ice melt and Arctic export advected along east Greenland are of similar magnitudes as those linked to runoff - Simulation of fresh coastal waters requires improved surface boundary conditions and/or models capable of representing mesoscale dynamics Accelerated melting of the Greenland Ice Sheet is considered a tipping element in the freshwater balance of the subpolar North Atlantic (SPNA). The East Greenland Current (EGC) and Coastal Current (EGCC) are the major conduits for transporting Arctic-sourced and Greenland glacial freshwater. Understanding freshwater changes in the EGC system and drivers thereof is crucial for connecting tipping elements in the SPNA. Using the eddy-rich model VIKING20X (1/20°) and Copernicus GLORYS12 (1/12°), we find that from 1993 to 2019 freshwater remains close to the shelf with interannual extremes in freshwater content (FWC) attributable to the imprint of Greenland melt only in years 2010 and 2012. Runoff increased significantly from 1995 to 2005 and Arctic freshwater export after 2005. Overall, regional wind patterns, sea ice melt and increasingly glacial ice and snow meltwater runoff along with the Arctic-sourced Polar Water set interannual FWC variations in the EGC system. We emphasize that these freshwater sources have different seasonal timing. South of 65°N sea ice melts year round and retreats to north of 65°N, where melt in summer prevails. Greenland runoff peaks in June–August with only some locations of year round discharge. Alongshore winds intensify in fall and winter where reduced onshore Ekman transport allows for freshwater to spread laterally in the EGC. We show that sea ice melt, runoff and wind can cause interannual variations of comparable magnitude. All of which makes attributing ocean freshening events to Greenland meltwater inflow at current magnitudes a major challenge.

Description: The Christiana‐Santorini‐Kolumbo (CSK) volcanic field has hosted more than 100 explosive eruptions in the past 250,000 years, including the 1650 CE eruption of Kolumbo Volcano. Previous studies have established a link between regional tectonics and volcanism in the CSK volcanic field. While 2D seismic reflection data have given valuable insight into regional faulting, detailed fault zone characterization has been precluded by the sparsely spaced profiles. Using 3D seismic reflection data around Kolumbo Volcano, we provide the first 3D characterization of fault zones in the CSK volcanic field. Beneath the volcano's northwestern flank, and farther to the northwest, normal faults are predominantly NE‐SW trending, with mean fault trends between 044° and 049°. Normal faults beneath the southeastern flank are slightly more north‐oriented, with mean fault trends between 028° and 038°. Our detailed fault zone analysis reveals clear NW‐SE directed extension around the volcano, consistent with published focal mechanisms from microseismicity. The Kolumbo Fault Zone, ∼6 km northwest of Kolumbo Volcano, is characterized by distinct relay ramps between major overstepping normal faults. Regional 2D seismic profiles reveal a previously undocumented volcanic cone directly above the fault zone. Magma ascent to this cone has likely exploited enhanced vertical permeability associated with distributed deformation within a relay ramp. We suggest that fault relay structures may play an important role, over a range of spatial scales, in focusing magma ascent within the CSK volcanic field. Plain Language Summary In the last 250,000 years, more than 100 explosive eruptions have occurred in the “Christiana‐Santorini‐Kolumbo” volcanic field in the Aegean Sea. Eruptions like these represent a serious natural hazard for the region. In this study, we explored how tectonic processes are related to volcanic activity. We did this by studying tectonic deformation around the submarine Kolumbo Volcano, which last erupted violently in 1650 CE. We used three‐dimensional (3D) seismic reflection data, which provide high‐resolution imagery of the seafloor and underlying sediments. The data set shows how the sediments beneath the seafloor have been disrupted by tectonic faults, which have formed as the crust is being slowly pulled apart (extended). The orientations of the faults show that extension in and around the volcano is happening along a northwest to southeast orientation. Based on our new data, we suggest that the movement of magma through the crust might occur preferentially through structural features called “relay ramps.” Relay ramps are regions of complex tectonic deformation that exist between overlapping extensional faults. Our 3D imagery of fault zones in this volcanic field gives a better understanding of how tectonic and volcanic processes interact with each other. Key Points 3D seismic data reveal unprecedented detail of normal faulting around the submarine Kolumbo Volcano, Aegean Sea Long-term extension (NW-SE oriented) around Kolumbo Volcano is consistent with previous studies of seismicity and field mapping on Santorini Relay ramps accommodate strain in step-overs between normal faults and may be exploited as permeable zones for vertical magma ascent

Description: Marine isoprene plays a crucial role in the formation of secondary organic aerosol within the remote marine boundary layer. Due to scarce field measurements of oceanic isoprene and limited laboratory-based studies of isoprene production, assessing the importance of marine isoprene on atmospheric chemistry and climate is challenging. Calculating in-field isoprene production rates is a crucial step to predict marine isoprene concentrations and the subsequent emissions to the atmosphere. The distribution, sources, and dominant environmental factors of isoprene were determined in the Northwest Pacific Ocean in 2019. The nutrient enrichment in the Kuroshio Oyashio Extension (KOE) surface seawater, driven by the upwelling and atmospheric deposition, promoted the growth of phytoplankton and elevated the isoprene concentration. This was confirmed by observed responses of isoprene to nutrients and aerosol dust additions in a ship-based incubation experiment, where the isoprene concentrations increased by 70% (t = 4.417, p 〈 0.001) and 35% (t = 2.387, p 〈 0.05), respectively. Biogenic isoprene production rates in the deck incubation experiments were positively related to chlorophyll a, temperature, and solar radiation, with an average production of 7.33 +/- 4.27 pmol L (-1) day (-1). Photochemical degradation of dissolved organic matter was likely an abiotic source of isoprene, contributing to approximately 14% of the total production. Driven by high isoprene production and extreme physical disturbance, the KOE showed very high emissions of isoprene of 46.0 +/- 13.0 nmol m(-2) day (-1), which led to a significant influence on the oxidative capacity of the local atmosphere.

Description: During the Cretaceous, there were two factors that had important influences on the East Asian climate, the East Asian coastal mountains and Earth's orbital cycling. An important question is how the coastal mountains modulated the variability of East Asian climate over orbital timescales. Here, we perform simulations with the coastal mountains of 0, 2, and 4 km high and three orbital configurations to answer the question. Our results show that a mountain range at the East Asian coast can amplify the impacts of orbital forcing on East Asian climate. Specifically, precipitation over the Songliao Basin in Northeastern China has significant changes as the coastal mountain range is about 4 km high. Combining our simulation results with orbitally‐controlled sedimentary deposits from the Songliao Basin, we conclude that the altitude of the coastal mountain range was very likely higher than 2 km in the Late Cretaceous. Plain Language Summary Tectonic events and solar insolation are the two important factors impacting variations of the climate system in the geological past. Regional climate responses to variations in the radiation from the sun over 10 4 –10 5 years were often magnified or dampened by tectonic events. Cretaceous sedimentary records in East Asia suggest that East Asian climate was influenced by the solar insolation. Geological evidence showed that a mountain range existed along the East Asian coast then. Would this mountain range modulate impacts of solar insolation on East Asian climate? Our modeling results show that the influence of solar insolation on East Asian climate can be amplified by the coastal mountain range, depending on the mountain elevation. When the coastal mountain range is ∼2 km high, the amplification effects become significant. When its altitude reaches ∼4 km, the response of East Asian climate to solar insolation is considerably strengthened, and such a condition is supported by the rhythm induced by the climate variation due to solar insolation archived in the Cretaceous strata in the Songliao Basin. Thus, we speculate that the East Asian coastal mountains might have reached an altitude more than 2 km in the Late Cretaceous. Key Points East Asian climate was sensitive to orbital forcing in the Late Cretaceous East Asian coastal mountains amplified orbital forcing on East Asian climate variability East Asian coastal mountains were likely higher than 2 km in the Late Cretaceous

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: Arc‐backarc systems are inherently shaped by subduction, representing an essential window into processes acting in the Earth's interior such as the recycling of subducted slabs. Furthermore, they are setting where new crust is formed and are believed to be sites where juvenile continental crust emerges. We present a seismic refraction and wide‐angle velocity model across the Izu arc‐backarc system, and use its characteristic features to constrain geochemically and petrologically different compartments, revealing processes governing crustal formation overlying subduction zones. Our result delineates the Izu arc with a maximum thickness of ∼20 km and the Shikoku Basin with thicknesses of ∼7 to 11 km. In the volcanic arc, the middle crust of the felsic to intermediate tonalitic layer (6.0–6.5 km/s) is remarkably thicker beneath the basalt‐dominated area than in the rhyolite‐dominated area, indicating that basaltic volcanism is indispensable in the transformation process from arc to continental crust. However, rhyolitic volcanism may relate to the juvenile stage of arc evolution or the remelting of middle crust due to the insufficient supply of basaltic magma from the mantle. The mafic restite and cumulates, which used to be part of the arc crustal material, are delaminated and foundered into the mantle, forming extremely low mantle velocities (〈7.5 km/s). In the Shikoku Basin, our result supports a fertile mantle source with passive upwelling and normal temperature during the opening process, but the lack of high velocity in the lower crust rules out hydrous melts entrained from the subducting slab or anomalous mantle trapped during subduction zone reconfiguration. Plain Language Summary As a vital factor in supporting the conditions for the evolution of life and ecosystems, the origin and evolution of the continents are still enigmatic. Volcanic arcs are generally seen as a place for creating continental crust while recycling the incoming subducting slab. In this study, we present a seismic velocity structure model across the Izu arc and Shikoku Basin, offshore south of Japan, to demonstrate the rules contained behind the transformation from arc to continental crust. Our results support that basaltic volcanism in the volcanic arc nurtures the generation of felsic to intermediate rocks, which provides the bulk of the continental crust. During this process, other anti‐continent materials, like mafic rocks, tend to be foundered into the mantle. Therefore, we propose that constant basaltic volcanism is critical in transferring arc crust to continental crust. Key Points A long seismic refraction and wide‐angle profile presents the seismic structure across the Izu arc and Shikoku Basin The transformation from arc to continental crust is closely associated with basaltic volcanism from the rear arc to volcanic front Passive melting of a fertile mantle source under normal temperature governs the opening of the Shikoku Basin

Description: At the Blanco transform fault system (BTFS) off Oregon, 138 local earthquakes and 84 double-couple focal mechanisms from ocean-bottom-seismometer recordings jointly discussed with bathymetric features reveal a highly segmented transform system without prominent fracture zone traces longer than 100 km. In the west, seismicity is focused at deep troughs (i.e., the West and East Blanco, and Surveyor Depressions). In the east, the BTFS lacks a characteristic transform valley and instead developed the Blanco Ridge, which is the most seismically active feature, showing strike-slip and dip-slip faulting. Sandwiched between the two main segments of the BTFS is the Cascadia depression, representing a short intra-transform spreading segment. Seismic slip vectors reveal that stresses at the eastern BTFS are roughly in line with plate motion. In contrast, stresses to the west are clockwise skewed, indicating ongoing reorganization of the OTF system. As we observed no prominent fracture zones at the BTFS, plate tectonic reconstructions suggest that the BTFS developed from non-transform offsets rather than pre-existing transform faults during a series of ridge propagation events. Our observations suggest that the BTFS can be divided into two oceanic transform systems. The eastern BTFS is suggested to be a mature transform plate boundary since ∼0.6 Ma. In contrast, the western BTFS is an immature transform system, which is still evolving to accommodate far-field stress change. The BTFS acts as a natural laboratory to yield processes governing the development of oceanic transform faults. Key Points Local seismicity of the Blanco transform fault system (BTFS) reveals along-strike variations dominated by strike-slip and oblique dip-slip The BTFS developed from non-transform offsets rather than discrete transform faults in response to plate rotation and ridge propagation The BTFS consists of a mature plate boundary in the east and an immature system in the west, separated by a central spreading center

Description: Quantifying changes in oceanic aerobic respiration is essential for understanding marine deoxygenation. Here we use an Earth system model to investigate if and to what extent oxygen utilization rate (OUR) can be used to track the temporal change of true respiration (Rtrue). Rtrue results from the degradation of particulate and dissolved organic matter in the model ocean, acting as ground truth to evaluate the accuracy of OUR. Results show that in thermocline and intermediate waters of the North Atlantic Subtropical Gyre (200–1,000 m), vertically integrated OUR and Rtrue both decrease by 0.2 molO2/m2/yr from 1850 to 2100 under global warming. However, in the mesopelagic Tropical South Atlantic, integrated OUR increases by 0.2 molO2/m2/yr, while the Rtrue integral decreases by 0.3 molO2/m2/yr. A possible reason for the diverging OUR and Rtrue is ocean mixing, which affects water mass composition and maps remote respiration changes to the study region. Key Points: - Our model study confirms earlier findings that oxygen utilization rate (OUR) underestimates true respiration (Rtrue) in mesopelagic ocean - Despite OUR underestimate Rtrue, OUR can adequately estimate long-term changes in Rtrue in the mesopelagic North Atlantic subtropical gyre - OUR cannot adequately estimate climate-driven changes in Rtrue in the mesopelagic tropical South Atlantic where different water masses mix

Description: We present a model of thermal lithospheric thickness (the depth where the geotherm reaches a temperature of 1300°C) and surface heat flow in Tibet and adjacent regions based on a new thermal-isostasy method. The method accounts for crustal density heterogeneity, is free from any assumption of a steady-state lithosphere thermal regime, and assumes that deviations from crustal Airy-type isostasy are caused by lithosphere thermal heterogeneity. We observe a highly variable lithospheric thermal structure which we interpret as representing longitudinal variations in the northern extent of the subducting Indian plate, southward subduction of the Asian plate beneath central Tibet, and possible preservation of fragmented Tethyan paleo-slabs. Cratonic-type cold and thick lithosphere (200-240 km) with a predicted surface heat flow of 40-50 mW/m2 typifies the Tarim Craton, the northwest Yangtze Craton, and most of the Lhasa Block that is likely refrigerated by underthrusting Indian lithosphere. We identify a ‘North Tibet anomaly’ with thin (〈80 km) lithosphere and high surface heat flow (〉80-100 mW/m2). We interpret this anomaly as the result of removal of lithospheric mantle and asthenospheric upwelling at the junction of the Indian and Asian slabs with opposite subduction polarities. Other parts of Tibet typically have intermediate lithosphere thickness of 120-160 km and a surface heat flow of 45-60 mW/m2, with patchy anomalies in eastern Tibet. While different uplift mechanisms for Tibet predict different lithospheric thermal regimes, our results in terms of a highly variable thermal structure beneath Tibet suggest that topographic uplift is caused by an interplay of several mechanisms. Key Points Thick Tibetan lithosphere defines the longitudinally variable northern extent of the Indian plate The ‘North Tibet Anomaly’ with a hot mantle marks the junction of the Indian and Asian slabs Tethyan paleo-slabs explain variable lithosphere thickness in eastern Tibet

Description: Export of sinking particles from the surface ocean is critical for carbon sequestration and to provide energy to the deep biosphere. The magnitude and spatial patterns of this export have been estimated in the past by in situ particle flux observations, satellite-based algorithms, and ocean biogeochemical models; however, these estimates remain uncertain. Here, we use a recent machine learning reconstruction of global ocean particle size distributions (PSDs) from Underwater Vision Profiler 5 measurements to estimate carbon fluxes by sinking particles (35 μm–5 mm equivalent spherical diameter) from the surface ocean. We combine global maps of PSD properties with empirical relationships constrained against in situ flux observations to calculate particulate carbon export from the euphotic zone (5.8 ± 0.1 Pg C y−1) and annual maximum mixed layer depths (6.1 ± 0.1 Pg C y−1). The new flux reconstructions suggest a less variable seasonal cycle in the tropical ocean and a more persistent export in the Southern Ocean than previously recognized. Smaller particles (less than 418 μm) contribute most of the flux globally, while larger particles become more important at high latitudes and in tropical upwelling regions. Export from the annual maximum mixed layer exceeds that from the euphotic zone over most of the low-latitude ocean, suggesting shallow particle recycling and net heterotrophy in the deep euphotic zone. These estimates open the way to fully three-dimensional global reconstructions of particle fluxes in the ocean, supported by the growing database of in situ optical observations.

Description: Until now, few offshore seismic studies have acquired simultaneously P- and S- wave data to derive in detail the seismic structure of the oceanic crust. We present 2-D Vp and Vs models using wide-angle seismic data at the Indian basin adjacent to the NinetyEast Ridge. Here, an outcrop basement located at the middle of the seismic line presents uppermost crustal Poisson's ratios (ν) of 0.28–0.29 (Vp ∼ 4.2 km/s and Vs ∼ 2.3 km/s). At the flanks of the outcrop basement, the sediment cover is 200–300 m thick and ν values are similar (0.28–0.3), but Vp and Vs values are higher (4.5–4.8 and 2.4–2.6 km/s, respectively). We interpret the relatively lower Vp and Vs around the basement outcrop in terms of hydrothermal alteration, while at the flanks of the basement outcrop, hydrothermal alteration has most likely ceased by sedimentation and compaction processes. Across the seismic layer 2, the Vp–Vs trend is linear and follows a ν value of 0.28–0.29, however, at the seismic layer 2/3 transition, the Vp–Vs trend abruptly changes following a ν value of 0.25–0.26. These reduced observed ν values at the layer 2/3 transition are lower than those reported by laboratory measurements for gabbro (ν ∼ 0.293) and are interpreted in terms of epidotization at the dike-gabbro contact and/or crack-change properties around the lower part of the intrusive sheeted dike section. Key Points We obtain 2-D Vp and Vs models from active seismic data for the Indian oceanic crust The seismic models suggest hydrothermal alteration near a basement outcrop Poisson's ratios change at the layer 2/3 transition from 0.28–0.29 to 0.25–0.26

Description: Bacterial metabolism largely drives the sequestration of refractory organic matter in the ocean. However, a lack of understanding exists regarding the abundance and reactivity of bacterial particulate organic matter (POM). Here we report the bacterial contributions to suspended POM collected in the oligotrophic Western Pacific Warm Pool (WPWP). Around 27% of particulate organic carbon (POC) and ∼39% of particulate nitrogen (PN) in the surface ocean were derived from bacteria. Most of the bacterial POM (∼87%) was labile or semi-labile, and ∼85% of bacterial POM was removed between depths of ∼100–300 m. Bacterial POM constituted only ∼8% and ∼13% of refractory POC and PN, respectively. The rapid cycling of bacterial POM in upper waters was likely related to oligotrophic conditions and facilitated by higher temperatures in the WPWP. Taken together, these observations indicate that bacterial POM plays a crucial role in supplying energy for bacterial respiration. Key Points We assess bacterial contributions to suspended particulate organic matter (POM) in the Western Pacific Warm Pool on the basis of D-amino acid biomarkers Bacterial organics constitute 27% of surface ocean particulate organic carbon (POC) and 39% of particulate nitrogen (PN), but majority (∼87%) is labile or semi-labile Rapid cycling of bacterial POM in the upper ocean results in a contribution of only ∼8% to refractory POC and ∼13% to PN

Description: Understanding the future fate of the Greenland Ice Sheet (GIS) in the context of anthropogenic CO2 emissions is crucial to predict sea level rise. With the fully coupled Earth system model of intermediate complexity CLIMBER-X, we study the stability of the GIS and its transient response to CO2 emissions over the next 10 Kyr. Bifurcation points exist at global temperature anomalies of 0.6 and 1.6 K relative to pre-industrial. For system states in the vicinity of the equilibrium ice volumes corresponding to these temperature anomalies, mass loss rate and sensitivity of mass loss to cumulative CO2 emission peak. These critical ice volumes are crossed for cumulative emissions of 1,000 and 2,500 GtC, which would cause long-term sea level rise by 1.8 and 6.9 m respectively. In summary, we find tipping of the GIS within the range of the temperature limits of the Paris agreement. Key Points Bifurcation points exist at global mean temperature anomalies of 0.6 and 1.6 K relative to pre-industrial Mass loss rate and sensitivity to cumulative CO2 emission peak near the equilibrium ice volumes belonging to these temperature anomalies Substantial long-term mass loss of the Greenland ice sheet for cumulative emissions larger than 1,000 Gt carbon

Description: This study addresses the evolution of global tidal dynamics since the Last Glacial Maximum focusing on the extraction of tidal levels that are vital for the interpretation of geologic sea-level markers. For this purpose, we employ a truly-global barotropic ocean tide model which considers the non-local effect of Self-Attraction and Loading. A comparison to a global tide gauge data set for modern conditions yields agreement levels of 65%–70%. As the chosen model is data-unconstrained, and the considered dissipation mechanisms are well understood, it does not have to be re-tuned for altered paleoceanographic conditions. In agreement with prior studies, we find that changes in bathymetry during glaciation and deglaciation do exert critical control over the modeling results with minor impact by ocean stratification and sea ice friction. Simulations of 4 major partial tides are repeated in time steps of 0.5–1 ka and augmented by 4 additional partial tides estimated via linear admittance. These are then used to derive time series from which the tidal levels are determined and provided as a global data set conforming to the HOLSEA format. The modeling results indicate a strengthened tidal resonance by M2, but also by O1, under glacial conditions, in accordance with prior studies. Especially, a number of prominent changes in local resonance conditions are identified, that impact the tidal levels up to several meters difference. Among other regions, resonant features are predicted for the North Atlantic, the South China Sea, and the Arctic Ocean. Key Points Evolution of four major partial tides from Last Glacial Maximum until present times Validation of the employed ocean tide model with present-day tide gauge data and dissipation rates Diligent derivation of global tidal levels for the interpretation of sea level indexpoints

Description: Ecological interactions among phytoplankton occur in a moving fluid environment. Oceanic flows can modulate the competition and coexistence between phytoplankton populations, which in turn can affect ecosystem function and biogeochemical cycling. We explore the impact of submesoscale velocity gradients on phytoplankton ecology using observations, simulations, and theory. Observations reveal that the relative abundance of Synechoccocus oligotypes varies on 1–10 km scales at an ocean front with submesoscale velocity gradients at the same scale. Simulations in realistic flow fields demonstrate that regions of divergence in the horizontal flow field can substantially modify ecological competition and dispersal on timescales of hours to days. Regions of positive (negative) divergence provide an advantage (disadvantage) to local populations, resulting in up to ∼20% variation in community composition in our model. We propose that submesoscale divergence is a plausible contributor to observed taxonomic variability at oceanic fronts, and can lead to regional variability in community composition.

Description: Key Points: The high-resolution eddy's 3D structures are estimated from the satellite observations based on the interior + surface quasigeostrophic (isQG) method Introducing the higher-resolution sea surface temperature observation to isQG method enhances the effective resolution of reconstruction A correction scheme in density reconstruction within the mixed layer is suggested to offset the absence of mixed layer dynamics By incorporating the high-resolution satellite remote sensing sea surface temperature (SST) with low-resolution sea level anomaly and sea surface salinity (SSS), this study explores the reconstructability of the three-dimensional (3D) eddy structures via the “interior + surface quasigeostrophic” (isQG) method in the South Indian Ocean. We apply the incorporation of high-resolution SST to improve the spatial resolution of the reconstruction. We also propose a correction scheme for density reconstruction within the mixed layer to offset the absence of mixed layer dynamics in the SQG framework. Comparison against the in situ observations demonstrates a satisfactory reconstructability for subsurface velocity and density anomalies. Statistically, the zonal velocity reconstructability outperforms its meridional counterpart and the corresponding velocity phase. The reconstructed shallow-layer velocity exhibits a superior skill in eddy-active regions, when compared with drogued drifter observations. Reconstructed subsurface velocities reproduce the spatial structures of eddy-induced velocity anomaly along the GO-SHIP observation transect, although present smaller magnitudes. Results demonstrate the potential applicability of the isQG method for reconstructing mesoscale eddies, particularly in the ocean at mid-to-high latitudes, where subsurface dynamics are strongly influenced by barotropic and the first baroclinic modes. With the upcoming high-advanced satellite observations, the isQG framework is expected to achieve better subsurface estimations.

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: Limited constraints on the variability of the deep-water production in the Labrador Sea complicate reconstructions of the strength of the Atlantic Meridional Overturning Circulation (AMOC) during the Late Quaternary. Large volumes of detrital carbonates were repeatedly deposited in the Labrador Sea during the last 32 kyr, potentially affecting radiogenic Nd isotope signatures. To investigate this the Nd isotope compositions of deep and intermediate waters were extracted from the authigenic Fe-Mn oxyhydroxide fraction, foraminiferal coatings, the residual silicates and leachates of dolostone grains. We provide a first order estimation of Nd release via dissolution of detrital carbonates and its contribution to the authigenic ԑNd signatures in the Labrador Sea. During the Last Glacial Maximum the Nd isotope signatures in the Labrador Sea would allow active water mass mixing with more radiogenic ɛNd values (-12.6 and -14) prevailing in its eastern part whereas less radiogenic values (ɛNd ∼ -18.4) were found on the western Labrador slope. The deposition of detrital carbonates during Heinrich stadials (2,1) was accompanied by negative detrital and authigenic Nd isotope excursions (ɛNd ∼ -31) that were likely controlled by dissolution of dolostone or dolostone associated mineral inclusions. This highly unradiogenic signal dominated the authigenic phases and individual water masses in the Labrador Sea, serving as potential source of highly unradiogenic Nd to the North Atlantic region, while exported southward. The Holocene authigenic ɛNd signatures of the coatings and leachates significantly differed from those of the detrital silicates, approaching modern bottom water mass signatures during the Late Holocene. Key Points - Estimation of Nd release via dissolution of detrital carbonates and its contribution to the authigenic ԑNd signatures in the Labrador Sea - Dissolution of detrital dolostones in the water column during Heinrich stadials at least partially controlled ɛNd signatures - During the LGM generally more radiogenic signatures possibly indicate active water mass advection and mixing in the Labrador Sea

Description: Mapping and sampling four sections of the slow-spreading Reykjanes Ridge provide insight into how tectonic and volcanic activity varies with distance from the Iceland plume. The studied areas are characterized by significant variations in water depth, lava chemistry, crustal thickness, thermal structure, and ridge morphology. For each study area, fault pattern and dimension, tectonic strain, seamount morphology, and density are inferred from 15 m-resolution bathymetry. These observations are combined with geochemical analysis from glass samples and sediment thickness estimations along Remotely Operated Vehicle-dive videos. They reveal that (a) tectonic and volcanic activity along the Reykjanes Ridge, do not systematically vary with distance from the plume center. (b) The tectonic geometry appears directly related to the deepening of the brittle/ductile transition and the maximum change in tectonic strain related to the rapid change in crustal thickness and the transition between axial-high and axial valley (∼59.5°N). (c) Across-axis variations in the fault density and sediment thickness provide similar widths for the neo-volcanic zone except in regions of increased seamount emplacement. (d) The variations in seamount density (especially strong for flat-topped seamounts) are not related to the distance from the plume but appear to be correlated with the interaction between the V-shape ridges (VSR) flanking the ridge and the ridge axis. These observations are more compatible with the buoyant upwelling melting instability hypothesis for VSR formation and suggest that buoyant melting instabilities create many small magma batches which by-pass the normal subaxial magmatic plumbing system, erupting over a wider-than-normal area. Key Points The distance from the plume center is not the only factor controlling tectonic and volcanic activity along the Reykjanes Ridge Fault dimensions are primarily controlled by the variation of crustal thermal structure with distance from the hotspot Flat-topped seamount abundances peak where a V-shaped ridge intersects the axis, consistent with a buoyant upwelling melting instability

Description: Three volcanic arcs have been the source of New Zealand's volcanic activity since the Neogene: Northland arc, Coromandel Volcanic Zone (CVZ) and Taupō Volcanic Zone (TVZ). The eruption chronology for the Quaternary, sourced by the TVZ, is well studied and established, whereas the volcanic evolution of the precursor arc systems, like the CVZ (central activity c. 18 to 2 Ma), is poorly known due to limited accessibility to, or identification of, onshore volcanic deposits and their sources. Here, we investigate the marine tephra record of the Neogene, mostly sourced by the CVZ, of cores from IODP Exp. 375 (Sites U1520 and U1526), ODP Leg 181 (Sites 1123, 1124 and 1125), IODP Leg 329 (Site U1371) and DSDP Leg 90 (Site 594) offshore of New Zealand. In total, we identify 306 primary tephra layers in the marine sediments. Multi-approach age models (e.g. biostratigraphy, zircon ages) are used in combination with geochemical fingerprinting (major and trace element compositions) and the stratigraphic context of each marine tephra layer to establish 168 tie-lines between marine tephra layers from different holes and sites. Following this approach, we identify 208 explosive volcanic events in the Neogene between c. 17.5 and 2.6 Ma. This is the first comprehensive study of New Zealand's Neogene explosive volcanism established from tephrochronostratigraphic studies, which reveals continuous volcanic activity between c. 12 and 2.6 Ma with an abrupt compositional change at c. 4.5 Ma, potentially associated with the transition from CVZ to TVZ. Key Points New Zealand's Neogene explosive volcanism based on the marine tephra record Geochemical fingerprinting of marine tephra layers across the study area to establish volcanic events Insights into geochemical variations with time, repose times and spatiotemporal distribution

Description: The Canadian Arctic Southern Beaufort Sea is characterized by prominent relict submarine permafrost and gas hydrate occurrences formed by subaerial exposure during extensive glaciations in Pliocene and Pleistocene. Submarine permafrost is still responding to the thermal change as a consequence of the marine transgression that followed the last glaciation. Submarine permafrost is still underexplored and is currently the focus of several research projects as its degradation releases greenhouse gases that contribute to climate change. In this study, seismic reflection indicators are used to investigate the presence of submarine permafrost and gas hydrates on the outer continental shelf where the base of permafrost is expected to cross-cut geological layers. To address the challenges of marine seismic data collected in shallow water environments, we utilize a representative synthetic model to assess the data processing and the detection of submarine permafrost and gas hydrate by seismic data. The synthetic model allows us to minimize the misinterpretation of acquisition and processing artifacts. In the field data, we identify features along with characteristics arising from the top and base of submarine permafrost and the base of the gas hydrate stability zone. This work shows the distribution of the present submarine permafrost along the southern Canadian Beaufort Sea region and confirms its extension to the outer continental shelf. It supports the general shape suggested by previous works and previously published numerical models. Key Points Seismic reflection data reveal occurrences and extent of submarine permafrost and associated gas hydrates at the Canadian Beaufort Shelf Synthetic modeling of permafrost and gas hydrate is required to assess seismic processing minimizing the potential for misinterpretation Indicators of top and base of permafrost and the base of gas hydrate stability support previously published numerical models

Description: Key Points: - Novel micro-analytical techniques allow seasonally resolved climate proxy data from varved marine sediments - Potential to generate seasonal and inter annual resolution sea surface temperature proxy time series spanning 〉1,000 years - Thorough assessment of processes that influence the climate signal recovered from proxies, validated with careful replication, is required Three recently published papers including Napier et al. (2022, https://doi.org/10.1029/2021PA004355) utilize novel microanalytical approaches with varved marine sediments to demonstrate the potential to reconstruct seasonal and inter-annual climate variability. Obtaining paleoclimate data at a resolution akin to the observational record is vitally important for improving our understanding of climate phenomena such as monsoons and modes of variability such as the El Niño Southern Oscillation, for which appraisals of past inter-annual variability is critical. The ability to generate seasonal and inter annual resolution sea surface temperature proxy time series spanning a thousand years or more is revolutionary and has the potential to fill gaps in our knowledge of climate variability. Although generally limited to sediments from regions with oxygen depleted bottom waters, there is great potential to integrate shorter seasonal resolution climate “snap shots” from other archives such as annually banded corals into composite time series. But as paleoceanographic data are used more by the observational and modeling fields, we make the case for conducting a thorough case-by-case assessment of the processes that influence the climate signal recovered from proxies, using careful replication to validate new approaches. Understanding or exploring the potential influence of processes which effectively filter the climate signal will lead to more quantitative paleoceanographic data that will better serve the broader climate science community.

Description: We assess the Southern Ocean CO2 uptake (1985–2018) using data sets gathered in the REgional Carbon Cycle Assessment and Processes Project Phase 2. The Southern Ocean acted as a sink for CO2 with close agreement between simulation results from global ocean biogeochemistry models (GOBMs, 0.75 ± 0.28 PgC yr−1) and pCO2-observation-based products (0.73 ± 0.07 PgC yr−1). This sink is only half that reported by RECCAP1 for the same region and timeframe. The present-day net uptake is to first order a response to rising atmospheric CO2, driving large amounts of anthropogenic CO2 (Cant) into the ocean, thereby overcompensating the loss of natural CO2 to the atmosphere. An apparent knowledge gap is the increase of the sink since 2000, with pCO2-products suggesting a growth that is more than twice as strong and uncertain as that of GOBMs (0.26 ± 0.06 and 0.11 ± 0.03 Pg C yr−1 decade−1, respectively). This is despite nearly identical pCO2 trends in GOBMs and pCO2-products when both products are compared only at the locations where pCO2 was measured. Seasonal analyses revealed agreement in driving processes in winter with uncertainty in the magnitude of outgassing, whereas discrepancies are more fundamental in summer, when GOBMs exhibit difficulties in simulating the effects of the non-thermal processes of biology and mixing/circulation. Ocean interior accumulation of Cant points to an underestimate of Cant uptake and storage in GOBMs. Future work needs to link surface fluxes and interior ocean transport, build long overdue systematic observation networks and push toward better process understanding of drivers of the carbon cycle. Key Points: - Ocean models and machine learning estimates agree on the mean Southern Ocean CO2 sink, but the trend since 2000 differs by a factor of two - REgional Carbon Cycle Assessment and Processes Project Phase 2 estimates a 50% smaller Southern Ocean CO2 sink for the same region and timeframe as RECCAP1 - Large model spread in summer and winter indicates that sustained efforts are required to understand driving processes in all seasons

Description: The coupling at the interface between tectonic plates is a key geophysical parameter to capture the frictional locking across plate boundaries and provides a means to estimate where tectonic strain is accumulating through time. Here, we use both interferometric radar (InSAR) and Global Navigation Satellite System (GNSS) data to investigate the plate coupling of the Hikurangi subduction zone beneath the North Island of New Zealand, where multiple slow slip cycles are superimposed on the long‐term loading. We estimate the plate coupling across the subduction zone over three multi‐year observational periods targeting different stages of the slow slip cycle. Our results highlight the importance of the observational time period when interpreting coupling maps, emphasizing the temporal variability of plate coupling. Leveraging multiple geodetic data sets, we demonstrate how InSAR provides powerful constraints on the spatial resolution of both plate coupling and slow fault slip, even in a region where a dense GNSS network exists. Plain Language Summary Plate coupling as a concept describes to what degree the boundaries between tectonic plates are locked and building up stress. Such accumulated stress (over hundreds to thousands of years) will eventually be released in earthquakes, and therefore provides important information about the potential for future earthquakes. Our study uses satellite data to investigate how coupling between tectonic plates along the Hikurangi subduction zone (New Zealand's largest and most dangerous plate boundary fault) changes with time. We analyzed Interferometric Synthetic Aperture Radar and Global Navigation Satellite System data to map the areas where the plates are stuck together (coupled) and where they move past each other (uncoupled). We show that plate coupling varies significantly in space over 2, 4, and 10‐year time scales, highlighting the importance of carefully considering the observational time period when interpreting coupling maps. Key Points Integration of high‐resolution displacement maps from radar imagery captures plate coupling at fine scales Estimates of plate coupling depend strongly on the time period over which surface velocities are measured Temporal variations in plate coupling highlight when and where slow slip dominates the slip budget

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: 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: 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: 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: 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: 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: 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: 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: 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: This study investigates the mesoscale dynamics involved in the 8–11 October 2008 unseasonably strong African dust episode, during which dust was transported to the Iberian Peninsula (IP). We employ observational datasets and a high-resolution Weather Research and Forecasting model coupled with Chemistry simulations. The analysis shows that during 0900–1200 UTC 9 October, a mesoscale convective system developed over the Atlas Mountains and resulted in a southwestward propagating convective cold pool outflow on the southern foothills of the Anti-Atlas, which lifted dust from the source region. Between 1200 and 1800 UTC 9 October, new moist convection was enhanced over the Atlas Mountains due to intensifying confluence among a heat low, moist southwesterly Atlantic sea-breeze front, and northeasterly flow associated with the convective cold pool near western Algeria. This new moist convection intensified the strength of the convective cold pool outflow and haboob, both of which continued propagating southwestward. At 1200 UTC 10 October, the low-pressure system migrated poleward on the southern slopes of the Anti-Atlas Mountains in association with a mountain-plains solenoidal circulation due to the daytime differential heating between the southern slopes of the Anti-Atlas and nearby atmosphere. The deepening low-pressure and strengthening Atlantic sea-breeze redirected an equatorward advancing dust plume into the poleward direction. The dust plume ultimately crossed the Saharan Atlas Mountains on 11 October and finally impacted the IP. Key Points: - WRF-Chem simulation of an unseasonably strong haboob on the southern slopes of the Atlas Mountains - The equatorward-advancing dust plume was recirculated in the poleward direction by an Atlantic sea-breeze front - The Atlantic sea-breeze front and an intensified upper-level cutoff vortex are instrumental for dust transport over the Iberian Peninsula

Description: Despite the implication of aerosols for the radiation budget, there are persistent differences in data for the aerosol optical depth (τ) for 1998–2019. This study presents a comprehensive evaluation of the large-scale spatio-temporal patterns of mid-visible τ from modern data sets. In total, we assessed 94 different global data sets from eight satellite retrievals, four aerosol-climate model ensembles, one operational ensemble product, two reanalyses, one climatology and one merged satellite product. We include the new satellite data SLSTR and aerosol-climate simulations from the Coupled Model Intercomparison Project Phase 6 (CMIP6) and the Aerosol Comparisons between Observations and Models Phase 3 (AeroCom-III). Our intercomparison highlights model differences and observational uncertainty. Spatial mean τ for 60°N – 60°S ranges from 0.124 to 0.164 for individual satellites, with a mean of 0.14. Averaged τ from aerosol-climate model ensembles fall within this satellite range, but individual models do not. Our assessment suggests no systematic improvement compared to CMIP5 and AeroCom-I. Although some regional biases have been reduced, τ from both CMIP6 and AeroCom-III are for instance substantially larger along extra-tropical storm tracks compared to the satellite products. The considerable uncertainty in observed τ implies that a model evaluation based on a single satellite product might draw biased conclusions. This underlines the need for continued efforts to improve both model and satellite estimates of τ, for example, through measurement campaigns in areas of particularly uncertain satellite estimates identified in this study, to facilitate a better understanding of aerosol effects in the Earth system. Key Points: - Present-day patterns in aerosol optical depth differ substantially between 94 modern global data sets - The range in spatial means from individual satellites is −11% to +17% of the multi-satellite mean - Spatial means from climate model intercomparison projects fall within the satellite range but strong regional differences are identified

Description: Solar radiation received at the Earth's surface (Rs) is comprised of two components, the direct radiation (Rd) and the diffuse radiation (Rf). Rd, the direct beam from the sun, is essential for concentrated solar power generation. Rf, scattered by atmospheric molecules, aerosols, or cloud droplets, has a fertilization effect on plant photosynthesis. But how Rd and Rf change diurnally is largely unknown owing to the lack of long-term measurements. Taking advantage of 22 years of homogeneous hourly surface observations over China, this study documents the climatological means and evolutions in the diurnal cycles of Rd and Rf since 1993, with an emphasis on their implications for solar power and agricultural production. Over the solar energy resource region, we observe a loss of Rd which is relatively large near sunrise and sunset at low solar elevation angles when the sunrays pass through the atmosphere on a longer pathway. However, the concentrated Rd energy covering an average 10-hr period around noon during a day is relatively unaffected. Over the agricultural crop resource region, the large amounts of clouds and aerosols scattering more of the incoming light result in Rf taking the main proportion of Rs during the whole day. Rf resources and their fertilization effect in the main crop region of China further enhances since 1993 over almost all hours of the day. Key Points: - The loss of direct radiation over China since 1993 is relatively large at sunrise and sunset with little effect on solar power generation - The diffuse component dominates solar radiation normally near sunrise and sunset, but for the whole day over the main sown area of China - The diffuse fraction is further enhanced in the main sown area of China over almost all hours of the day since 1993

Description: Key Points High-resolution reflection seismic data reveals that the internal architecture of the Kolumbo Volcanic Chain The Kolumbo Volcanic Chain evolved during two episodes along NE-SW striking normal faults A prominent volcanic ridge connects the Kolumbo Volcanic Chain with Santorini highlighting a former connection between both systems Abstract The Christiana-Santorini-Kolumbo volcanic field in the southern Aegean Sea is one of the most hazardous volcanic regions in the world. Forming the northeastern part of this volcanic field, the Kolumbo Volcanic Chain (KVC) comprises more than submarine volcanic cones. However, due to their inaccessibility, little is known about the spatio-temporal evolution and tectonic control of these submarine volcanoes and their link to the volcanic plumbing system of Santorini. In this study, we use multichannel reflection seismic imaging to study the internal architecture of the KVC and its link to Santorini. We show that the KVC evolved during two episodes, which initiated at ~1 Ma with the formation of mainly effusive volcanic edifices along a NE-SW trending zone. The cones of the second episode were formed mainly by submarine explosive eruptions between 0.7 and 0.3 Ma and partly developed on top of volcanic edifices from the first episode. We identify two prominent normal faults that underlie and continue the two main trends of the KVC, indicating a direct link between tectonics and volcanism. In addition, we reveal several buried volcanic centers and a distinct volcanic ridge connecting the KVC with Santorini, suggesting a connection between the two volcanic centers in the past. This connection was interrupted by a major tectonic event and, as a result, the two volcanic systems now have separate, largely independent plumbing systems despite their proximity

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: 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.