ALBERT

Description: Marine methane hydrate is an ice-like substance that is stable in sediment around marine continental margins where water depths are greater than ~450–700 m. The release of methane due to melting of hydrates is considered to be a mechanism for past global carbon-cycle perturbations and could exacerbate ongoing anthropogenic climate change. Increases in bottom-water temperature at the landward limit of marine hydrate around continental margins, where vulnerable hydrate exists at or below the seabed, cause methane to vent into the ocean. However, this setting represents only ~3.5% of the global hydrate reservoir. The potential for methane from hydrate in deeper water to reach the atmosphere was considered negligible. Here we use three-dimensional (3D) seismic imagery to show that, on the Mauritanian margin, methane migrated at least 40 km below the base of the hydrate stability zone and vented through 23 pockmarks at the shelf break, probably during warmer Quaternary interglacials. We demonstrate that, under suitable circumstances, some of the 96.5% of methane bound in deeper water distal hydrates can reach the seafloor and vent into the ocean beyond the landward limit of marine hydrate. This reservoir should therefore be considered for estimating climate change-induced methane release during a warming world.

Description: Methylmercury is a potent toxin threatening the global population mainly through the consumption of marine fish. Hydrothermal venting directly delivers natural mercury to the ocean, yet its global flux remains poorly constrained. To determine the extent to which anthropogenic inputs have increased oceanic mercury levels, it is crucial to estimate natural mercury levels. Here we combine observations of vent fluids, plume waters, seawater and rock samples to quantify the release of mercury from the Trans-Atlantic Geotraverse hydrothermal vent at the Mid-Atlantic Ridge. The majority (67–95%) of the mercury enriched in the vent fluids (4,966 ± 497 pmol l −1 ) is rapidly diluted to reach background seawater levels (0.80 pmol l −1 ). A small Hg fraction (2.6–10%) is scavenged to the Trans-Atlantic Geotraverse mound rocks. Scaling up our findings and previous work, we propose a mercury flux estimate of 1.5–64.7 t per year from mid-ocean ridges. This hydrothermal flux is small in comparison to anthropogenic inputs. This suggests that most of the mercury present in the ocean must be of anthropogenic origin and that the implementation of emissions reduction measures outlined in the Minamata Convention could effectively reduce mercury levels in the global ocean and subsequently in marine fish.

Description: The Mediterranean Sea has been sampled irregularly by research vessels in the past, mostly by national expeditions in regional waters. To monitor the hydrographic, biogeochemical and circulation changes in the Mediterranean Sea, a systematic repeat oceanographic survey programme called Med-SHIP was recommended by the Mediterranean Science Commission (CIESM) in 2011, as part of the Global Ocean Ship-based Hydrographic Investigations Program (GO-SHIP). Med-SHIP consists of zonal and meridional surveys with different frequencies, where comprehensive physical and biogeochemical properties are measured with the highest international standards. The first zonal survey was done in 2011 and repeated in 2018. In addition, a network of meridional (and other key) hydrographic sections were designed: the first cycle of these sections was completed in 2016, with three cruises funded by the EU project EUROFLEETS2. This paper presents the physical and chemical data of the meridional and key transects in the Western and Eastern Mediterranean Sea collected during those cruises.

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: Large explosive volcanic eruptions from island arcs pour pyroclastic currents into marine basins, impacting ecosystems and generating tsunamis that threaten coastal communities and infrastructures. Risk assessments require robust records of such highly hazardous events, which is challenging as most of the products lie buried under the sea. Here we report the discovery by IODP Expedition 398 of a giant rhyolitic pumice deposit emplaced 520 ± 10 ky ago at water depths of 200 to 1000 m during a high-intensity, shallow submarine eruption of ancestral Santorini Volcano. Pyroclastic currents discharged into the sea transformed into turbidity currents and slurries, forming a 〉89 ± 8 km 3 volcaniclastic megaturbidite up to 150 m thick in the surrounding marine basins, while breaching of the sea surface by the eruption column laid down veneers of ignimbrite on three islands. The eruption is one of the largest recorded on the South Aegean Volcanic Arc, and highlights the hazards from submarine explosive eruptions.

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: 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: Photosynthesis fuels primary production at the base of marine food webs. Yet, in many surface ocean ecosystems, diel-driven primary production is tightly coupled to daily loss. This tight coupling raises the question: which top-down drivers predominate in maintaining persistently stable picocyanobacterial populations over longer time scales? Motivated by high-frequency surface water measurements taken in the North Pacific Subtropical Gyre (NPSG), we developed multitrophic models to investigate bottom-up and top-down mechanisms underlying the balanced control of Prochlorococcus populations. We find that incorporating photosynthetic growth with viral- and predator-induced mortality is sufficient to recapitulate daily oscillations of Prochlorococcus abundances with baseline community abundances. In doing so, we infer that grazers in this environment function as the predominant top-down factor despite high standing viral particle densities. The model-data fits also reveal the ecological relevance of light-dependent viral traits and non-canonical factors to cellular loss. Finally, we leverage sensitivity analyses to demonstrate how variation in life history traits across distinct oceanic contexts, including variation in viral adsorption and grazer clearance rates, can transform the quantitative and even qualitative importance of top-down controls in shaping Prochlorococcus population dynamics.

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: Climate change is opening the Arctic Ocean to increasing human impact and ecosystem changes. Arctic fjords, the region’s most productive ecosystems, are sustained by a diverse microbial community at the base of the food web. Here we show that Arctic fjords become more prokaryotic in the picoplankton (0.2–3 µm) with increasing water temperatures. Across 21 fjords, we found that Arctic fjords had proportionally more trophically diverse (autotrophic, mixotrophic, and heterotrophic) picoeukaryotes, while subarctic and temperate fjords had relatively more diverse prokaryotic trophic groups. Modeled oceanographic connectivity between fjords suggested that transport alone would create a smooth gradient in beta diversity largely following the North Atlantic Current and East Greenland Current. Deviations from this suggested that picoeukaryotes had some strong regional patterns in beta diversity that reduced the effect of oceanographic connectivity, while prokaryotes were mainly stopped in their dispersal if strong temperature differences between sites were present. Fjords located in high Arctic regions also generally had very low prokaryotic alpha diversity. Ultimately, warming of Arctic fjords could induce a fundamental shift from more trophic diverse eukaryotic- to prokaryotic-dominated communities, with profound implications for Arctic ecosystem dynamics including their productivity patterns.

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 sailboat Seaexplorer collected underway sea surface partial pressure of CO2 (pCO2) data for 129 days (2018–2021), including an Antarctic circumnavigation. By comparing ensembles of data-driven air-sea CO2 fluxes computed with and without sailboat data and applying a detection algorithm, we show that these sailboat observations significantly increase the regional carbon uptake in the North Atlantic and decrease it in the Southern Ocean. While compensating changes in both basins limit the global effect, the Southern Ocean–particularly frontal regions (40°S–60°S) during summertime—exhibited the largest air-sea CO2 flux changes, averaging 20% of the regional mean. Assessing the sensitivity of the air-sea CO2 flux to measurement uncertainty, the results stay robust within the expected random measurement uncertainty (± 5 μatm) but remain undetectable with a measurement offset of 5 µatm. We thus conclude that sailboats fill essential measurement gaps in remote ocean regions.

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

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

Description: We present chromosome-level genome assemblies from representative species of three independently evolved seagrass lineages: Posidonia oceanica, Cymodocea nodosa, Thalassia testudinum and Zostera marina. We also include a draft genome of Potamogeton acutifolius, belonging to a freshwater sister lineage to Zosteraceae. All seagrass species share an ancient whole-genome triplication, while additional whole-genome duplications were uncovered for C. nodosa, Z. marina and P. acutifolius. Comparative analysis of selected gene families suggests that the transition from submerged-freshwater to submerged-marine environments mainly involved fine-tuning of multiple processes (such as osmoregulation, salinity, light capture, carbon acquisition and temperature) that all had to happen in parallel, probably explaining why adaptation to a marine lifestyle has been exceedingly rare. Major gene losses related to stomata, volatiles, defence and lignification are probably a consequence of the return to the sea rather than the cause of it. These new genomes will accelerate functional studies and solutions, as continuing losses of the ‘savannahs of the sea’ are of major concern in times of climate change and loss of biodiversity.

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

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

Description: The relationship between energy reserves of cold-water corals (CWCs) and their physiological performance remains largely unknown. In addition, it is poorly understood how the energy allocation to different metabolic processes might change with projected decreasing food supply to the deep sea in the future. This study explores the temporal and spatial variations of total energy reserves (proteins, carbohydrates and lipids) of the CWC Desmophyllum dianthus and their correlation with its calcification rate. We took advantage of distinct horizontal and vertical physico-chemical gradients in Comau Fjord (Chile) and examined the changes in energy reserves over one year in an in situ reciprocal transplantation experiment (20 m vs. 300 m and fjord head vs. mouth). Total energy reserves correlated positively with calcification rates. The fast-growing deep corals had higher and less variable energy reserves, while the slower-growing shallow corals showed pronounced seasonal changes in energy reserves. Novel deep corals (transplanted from shallow) were able to quickly increase both their calcification rates and energy reserves to similar levels as native deep corals. Our study shows the importance of energy reserves in sustaining CWC growth in spite of aragonite undersaturated conditions (deep corals) in the present, and potentially also future ocean.

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: 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: 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: 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: Weather causes extremes in photovoltaic and wind power production. Here we present a comprehensive climatology of anomalies in photovoltaic and wind power production associated with weather patterns in Europe considering the 2019 and potential 2050 installations, and hourly to ten-day events. To that end, we performed kilometer-scale numerical simulations of hourly power production for 23 years and paired the output with a weather classification which allows a detailed assessment of weather-driven spatio-temporal production anomalies. Our results highlight the dependency of low-power production events on the installed capacities and the event duration. South-shifted Westerlies (Anticyclonic South-Easterlies) are associated with the lowest hourly (ten-day) extremes for the 2050 (both) installations. Regional power production anomalies can differ from the ones in the European mean. Our findings suggest that weather patterns can serve as indicators for expected photovoltaic and wind power production anomalies and may be useful for early warnings in the energy sector.

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

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

Description: Enhancing ocean productivity by artificial upwelling is evaluated as a nature-based solution for food security and climate change mitigation. Fish production is intended through diatom-based plankton food webs as these are assumed to be short and efficient. However, our findings from mesocosm experiments on artificial upwelling in the oligotrophic ocean disagree with this classical food web model. Here, diatoms did not reduce trophic length and instead impaired the transfer of primary production to crustacean grazers and small pelagic fish. The diatom-driven decrease in trophic efficiency was likely mediated by changes in nutritional value for the copepod grazers. Whilst diatoms benefitted the availability of essential fatty acids, they also caused unfavorable elemental compositions via high carbon-to-nitrogen ratios (i.e. low protein content) to which the grazers were unable to adapt. This nutritional imbalance for grazers was most pronounced in systems optimized for CO 2 uptake through carbon-to-nitrogen ratios well beyond Redfield. A simultaneous enhancement of fisheries production and carbon sequestration via artificial upwelling may thus be difficult to achieve given their opposing stoichiometric constraints. Our study suggest that food quality can be more critical than quantity to maximize food web productivity during shorter-term fertilization of the oligotrophic ocean.

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

Description: The evolution of the northern hemispheric climate during the last glacial period was beset by quasi-episodic iceberg discharge events from the Laurentide ice sheet, known as Heinrich events (HEs). The paleo record places most HEs into the cold stadial of the Dansgaard-Oeschger cycle. However, not every Dansgaard-Oeschger cycle is associated with a HE, revealing a complex interplay between the two modes of glacial variability. Here, using a coupled ice sheet-solid earth model, we introduce a mechanism that explains the synchronicity of HEs and Dansgaard-Oeschger cycles. Unlike earlier studies, our mechanism does not require a trigger during the stadial. Instead, the atmospheric warming signal during the interstadial of the Dansgaard-Oeschger cycle causes enhanced ice stream thickening that leads to the HE during the late interstadial. We demonstrate that this mechanism reproduces the key HE characteristics and provides an explanation for synchronous HEs from different regions of the Laurentide ice sheet.

Description: There has been extensive research into the nonlinear responses of the Earth system to astronomical forcing during the last glacial cycle. However, the speed and spatial geometry of ice sheet expansion to its largest extent at the Last Glacial Maximum 21 thousand years ago remains uncertain. Here we use an Earth system model with interactive ice sheets to show that distinct initial North American (Laurentide) ice sheets at 38 thousand years ago converge towards a configuration consistent with the Last Glacial Maximum due to feedbacks between atmospheric circulation and ice sheet geometry. Notably, ice advance speed and spatial pattern in our model are controlled by the amount of summer snowfall, which is dependent on moisture transport pathways from the North Atlantic warm pool linked to ice sheet geometry. The consequence of increased summer snowfall on the surface mass balance of the ice sheet is not only the direct increase in accumulation but the indirect reduction in melt through the snow/ice–albedo feedback. These feedbacks provide an effective mechanism for ice growth for a range of initial ice sheet states and may explain the rapid North American ice volume increase during the last ice age and potentially driving growth during previous glacial periods.

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

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

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

Description: Due to climate change the pathogenic bacterium Vibrio vulnificus proliferates along brackish coastlines, posing risks to public health, tourism, and aquaculture. Here we investigated previously suggested regulation measures to reduce the prevalence of V. vulnificus, locally through seagrass and regionally through the reduction of eutrophication and consequential formation of algal blooms. Field samples collected in the summer of 2021 covered the salinity and eutrophication gradients of the Baltic Sea, one of the largest brackish areas worldwide. Physico-, biological- and hydrochemical parameters were measured and variables explaining V. vulnificus occurrence were identified by machine learning. The best V. vulnificus predictors were eutrophication-related features, such as particulate organic carbon and nitrogen, as well as occurrence of potential phytoplankton blooms and associated species. V. vulnificus abundance did not vary significantly between vegetated and non-vegetated areas. Thus, reducing nutrient inputs could be an effective method to control V. vulnificus populations in eutrophied brackish coasts.

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

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

Description: Abstract Robust chronologies and time equivalent tephra markers are essential to better understand spatial palaeoenvironmental response to past abrupt climatic changes. Identification of well-dated and widely dispersed volcanic ash by tephra and cryptotephra (microscopic volcanic ash) provide time synchronous tie-points and strongly reduce chronological uncertainties. Here, we present the major, minor and trace element analyses of cryptotephra shards in the Dead Sea Deep Drilling sedimentary record (DSDDP 5017-1A) matching the Campanian Ignimbrite (CI). This geochemical identification expands the know dispersal range of the CI to the southeastern Mediterranean, over 2,300 km from the volcanic source. Due to the CI eruption occurring near-synchronous with North Atlantic ice surge of Heinrich Event 4 (HE4), this tephra provides insights into regional responses to large-scale climatic change in the Mediterranean. In the Dead Sea, the CI layer is associated with wetter climatic conditions. This contrasts with the contemporaneous occurrence of the CI deposition and dry conditions in the northern and western Mediterranean suggesting a possible climate time-transgressive expansion of HE4. Our finding underscores the temporal and spatial complexity of regional climate responses and emphasises the importance of tephra as a time marker for studying large-scale climatic changes verses regional variations.

Description: As a consequence of on-going global warming, the ocean is losing oxygen, which has implications not only in terms of marine resources management and food supply but also in terms of the potentially important feedback on the global carbon cycle and climate. Of particular scrutiny are the extended zones of already low levels of oxygen called the oxygen minimum zones (OMZs) embedded in the subsurface waters of the productive Eastern Boundary Upwelling Systems (EBUS). These OMZs are currently diversely simulated by state-of-the-art Earth System Models (ESM) hampering a reliable projection of ocean deoxygenation on marine ecosystem services in these regions. Here we focus on the most emblematic EBUS OMZs of the planet, that of the South Eastern Pacific (SEP), which is under the direct influence of the El Ni & ntilde;o Southern Oscillation (ENSO), the main climate mode on interannual timescales at global scale. We show that, despite the low consensus among ESM long-term projections of oxygen levels, the sensitivity of the depth of the upper margin (oxycline) of the SEP OMZ to El Ni & ntilde;o events in an ensemble of ESMs can be used as a predictor of its long-term trend, which establishes an emergent constraint for the SEP OMZ. Because the oxycline along the coast of Peru and Chile deepens during El Ni & ntilde;o events, the upper bound of the SEP OMZ is thus likely to deepen in the future climate, therefore oxygenating the SEP OMZ. This has implications not only for understanding the nitrogen and carbon cycles at global scale but also for designing adaptation strategies for regional upper-ocean ecosystem services. The upper bound of the southeast Pacific oxygen minimum zone deepens during El Ni & ntilde;o events across an ensemble of Earth system models and is therefore projected to likely contract in the future climate, according to an analysis of multiple Earth System Models.

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

Description: Reef-building corals live in highly hydrodynamic environments, where water flow largely controls the complex chemical microenvironments surrounding them—the concentration boundary layer (CBL). The CBL may be key to alleviate ocean acidification (OA) effects on coral colonies by partially isolating them. However, OA effects on coral CBL remain poorly understood, particularly under different flow velocities. Here, we investigated these effects on the reef-building corals Acropora cytherea, Pocillopora verrucosa, and Porites cylindrica. We preconditioned corals to a control (pH 8.0) and OA (pH 7.8) treatment for four months and tested how low flow (2 cm s−1) and moderate flow (6 cm s−1) affected O2 and H+ CBL traits (thickness, surface concentrations, and flux) inside a unidirectional-flow chamber. We found that CBL traits differed between species and flow velocities. Under OA, traits remained generally stable across flows, except surface pH. In all species, the H+ CBL was thin and led to lower surface pH. Still, low flow thickened H+ CBLs and increased light elevation of surface pH. In general, our findings reveal a weak to null OA modulation of the CBL. Moreover, the OA-buffering capacity by the H+ CBL may be limited in coral species, though low flow could enhance CBL sheltering.

Description: Turbulent mixing in the ocean, lakes and reservoirs facilitates the transport of momentum, heat, nutrients, and other passive tracers. Turbulent fluxes are proportional to the rate of turbulent kinetic energy dissipation per unit mass, ε. A common method for ε measurements is using microstructure profilers with shear probes. Such measurements are now widespread, and a non-expert practitioner will benefit from best practice guidelines and benchmark datasets. As a part of the Scientific Committee on Oceanographic Research (SCOR) working group on “Analysing ocean turbulence observations to quantify mixing” (ATOMIX), we compiled a collection of five benchmark data of ε from measurements of turbulence shear using shear probes. The datasets are processed using the ATOMIX recommendations for best practices documented separately. Here, we describe and validate the datasets. The benchmark collection is from different types of instruments and covers a wide range of environmental conditions. These datasets serve to guide the users to test their ε estimation methods and quality-assurance metrics, and to standardize their data for archiving.

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

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

Description: Age and longevity are key parameters for demography and life-history evolution of organisms. In clonal species, a widespread life history among animals, plants, macroalgae and fungi, the sexually produced offspring (genet) grows indeterminately by producing iterative modules, or ramets, and so obscure their age. Here we present a novel molecular clock based on the accumulation of fixed somatic genetic variation that segregates among ramets. Using a stochastic model, we demonstrate that the accumulation of fixed somatic genetic variation will approach linearity after a lag phase, and is determined by the mitotic mutation rate, without direct dependence on asexual generation time. The lag phase decreased with lower stem cell population size, number of founder cells for the formation of new modules, and the ratio of symmetric versus asymmetric cell divisions. We calibrated the somatic genetic clock on cultivated eelgrass Zostera marina genets (4 and 17 years respectively). In a global data set of 20 eelgrass populations, genet ages were up to 1,403 years. The somatic genetic clock is applicable to any multicellular clonal species where the number of founder cells is small, opening novel research avenues to study longevity and, hence, demography and population dynamics of clonal species.

Description: Observations of spatio-temporal variability of the deep ocean are rare and little is known about occurrence of deep ocean mesoscale dynamics. Here, we make use of 2.5 years of time series data from three distributed sensor arrays, which acquired high-resolution temperature, pressure and sound speed data of the bottom layer offshore northern Chile. Estimating salinity and density from the direct observations enable access to the full spectrum of hydrographic variability from a multi-hourly to annual time scale and with average inter-station distances of less than 1 km. Analyses revealed interannual warming over the continental slope of 0.002 °C yr−1–0.003 °C yr−1, and could trace periodic hydrographic anomalies, likely related to coastal-trapped waves, as far as to the lower continental slope. A concurrent change in the shape of the warm anomalies and the rate of deep-sea warming that occurs with the crossing of the deep-sea trench suggests that the abyssal part of the eastern boundary current system off Chile does not extend past the deep sea trench. Furthermore, the comparison of anomaly timing and shape in between stations implies southwards flow over the mid to lower continental slope, centred closer to the trench.

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

Description: Caldera-forming eruptions of silicic volcanic systems are among the most devastating events on Earth. By contrast, post-collapse volcanic activity initiating new caldera cycles is generally considered less hazardous. Formed after Santorini’s latest caldera-forming eruption of ~1600 bce , the Kameni Volcano in the southern Aegean Sea enables the eruptive evolution of a recharging multi-cyclic caldera to be reconstructed. Kameni’s eruptive record has been documented by onshore products and historical descriptions of mainly effusive eruptions dating back to 197 bce . Here we combine high-resolution seismic reflection data with cored lithologies from International Ocean Discovery Program Expedition 398 at four sites to determine the submarine architecture and volcanic history of intra-caldera deposits from Kameni. Our shore-crossing analysis reveals the deposits of a submarine explosive eruption that produced up to 3.1 km 3 of pumice and ash, which we relate to a historical eruption in 726 ce . The estimated volcanic explosivity index of magnitude 5 exceeds previously considered worst-case eruptive scenarios for Santorini. Our finding that the Santorini caldera is capable of producing large explosive eruptions at an early stage in the caldera cycle implies an elevated hazard potential for the eastern Mediterranean region, and potentially for other recharging silicic calderas.

Description: Accurate reconstruction of seasonal atmospheric patterns of the past is essential for reliable prediction of how climate will evolve due to anthropogenic CO2 forcing. The Eocene ‘hot house’ climate, as the warmest epoch during the Cenozoic, is considered as a potential analogue for ‘high-CO2’ future climate scenarios. In this context, the reconstruction of variations in seasonality are as important as changes in mean annual conditions. Here we combine stable oxygen (δ18O) and dual clumped isotope (Δ47 + Δ48) measurements of a bivalve shell to determine sub-annual variations in sea surface temperatures and oceanic freshening in the Paris Basin during the Mid-Eocene Climate Optimum, 40 million years ago. Our reconstruction indicates to high mean annual temperatures with a small seasonal amplitude (33.3 °C ± 4.4 °C) and an enhanced fresh water input during the summer period. Our results implying a substantially warmer climate state with different hydrological conditions for Western Europe during the Eocene than previously suggested by proxy data or climate modelling.

Description: Thriving in both epipelagic and mesopelagic layers, Rhizaria are biomineralizing protists, mixotrophs or flux-feeders, often reaching gigantic sizes. In situ imaging showed their contribution to oceanic carbon stock, but left their contribution to element cycling unquantified. Here, we compile a global dataset of 167,551 Underwater Vision Profiler 5 Rhizaria images, and apply machine learning models to predict their organic carbon and biogenic silica biomasses in the uppermost 1000 m. We estimate that Rhizaria represent up to 1.7% of mesozooplankton carbon biomass in the top 500 m. Rhizaria biomass, dominated by Phaeodaria, is more than twice as high in the mesopelagic than in the epipelagic layer. Globally, the carbon demand of mesopelagic, flux-feeding Phaeodaria reaches 0.46 Pg C y −1 , representing 3.8 to 9.2% of gravitational carbon export. Furthermore, we show that Rhizaria are a unique source of biogenic silica production in the mesopelagic layer, where no other silicifiers are present. Our global census further highlights the importance of Rhizaria for ocean biogeochemistry.

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

Description: The Gulf of St. Lawrence is increasingly affected by bottom water hypoxia; however, the timescales and pathways of deep water transport remain unclear. Here, we present results from the Deep Tracer Release eXperiment (TReX Deep), during which an inert SF 5 CF 3 tracer was released inshore of Cabot Strait at 279 m depth to investigate deep inflow transport and mixing rates. Dispersion was also assessed via neutrally-buoyant Swish floats. Our findings indicate that the tracer moves inland at 0.5 cm s −1 , with an effective lateral diffusivity of 2 × 10 2 m 2 s −1 over 1 year. Simplified 1D simulations suggest inflow water should reach the estuary head in 1.7 years, with the bulk arriving after 4.7 years. Basin-wide effective vertical diffusivity is around 10 −5 m 2 s −1 over 1 year; however, vertical diffusivity increases near the basin slopes, suggesting that turbulent boundary processes influence mixing. These results are compared to Lagrangian simulations in a regional 3D model to evaluate the capacity to model dispersion in the Gulf.

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: 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: While basaltic volcanism is dominate during rifting and continental breakup, felsic magmatism may also comprise important components of some rift margins. During International Ocean Discovery Program (IODP) Expedition 396 on the continental margin of Norway, a graphite-garnet-cordierite bearing dacitic, pyroclastic unit was recovered within early Eocene sediments on Mimir High (Site U1570), a marginal high on the Vøring transform margin. Here, we present a comprehensive textural, mineralogical, and petrological study of the dacite in order to assess its melting origin and emplacement. The major mineral phases (garnet, cordierite, quartz, plagioclase, alkali feldspar) are hosted in a fresh rhyolitic, highly vesicular, glassy matrix, locally mingled with sediments. The xenocrystic major element chemistry of garnet and cordierite, the presence of zircon inclusions with inherited cores, and thermobarometric calculations all support a crustal metapelite origin. While most magma-rich margin models favor crustal anatexis in the lower crust, thermobarometric calculations performed here show that the dacite was produced at upper-crustal depths (〈 5 kbar) and high temperature (750–800 °C) with up to 3 wt% water content. In situ U-Pb analyses on zircon inclusions give a magmatic age of 54.6 ± 1.1 Ma, revealing the emplacement of the dacite post-dates the Paleocene-Eocene Thermal Maximum (PETM). Our results suggest that the opening of the North Atlantic was associated with a phase of low-pressure, high-temperature crustal melting at the onset of the main phase of magmatism.

Description: Volcanic flank collapses, especially those in island settings, have generated some of the most voluminous mass transport deposits on Earth and can trigger devastating tsunamis. Reliable tsunami hazard assessments for flank collapse-driven tsunamis require an understanding of the complex emplacement processes involved. The seafloor sequence southeast of Montserrat (Lesser Antilles) is a key site for the study of volcanic flank collapse emplacement processes that span subaerial to submarine environments. Here, we present new 2D and 3D seismic data as well as MeBo drill core data from one of the most extensive mass transport deposits offshore Montserrat, which exemplifies multi-phase landslide deposition from volcanic islands. The deposits reveal emplacement in multiple stages including two blocky volcanic debris avalanches, secondary seafloor failure and a late-stage erosive density current that carved channel-like incisions into the hummocky surface of the deposit about 15 km from the source region. The highly erosive density current potentially originated from downslope-acceleration of fine-grained material that was suspended in the water column earlier during the slide. Late-stage erosive turbidity currents may be a more common process following volcanic sector collapse than has been previously recognized, exerting a potentially important control on the observed deposit morphology as well as on the runout and the overall shape of the deposit. Key Points Landslide emplacement offshore Montserrat included volcanic flank collapses, sediment incorporation, and a late-stage erosive flow Highly erosive flows are likely to be common processes during volcanic flank collapse deposition Pre-existing topography plays a major role in shaping flank collapse-associated mass transport deposits Plain Language Summary Disintegration of volcanic islands can cause very large landslides and destructive tsunamis. To assess the tsunami hazard of such events, it is crucial to understand the processes that are involved in their formation. We present new insights from seismic data and drill cores from a landslide deposit offshore Montserrat, a volcanic island in the Lesser Antilles Arc in the Caribbean. Our analysis reveals the emplacement of landslide material in several stages, including multiple volcanic flank collapses, incorporation of seafloor sediments and an erosive flow that carved channels into the top of the deposit right after its emplacement. We suggest that highly erosive flows are a common process during volcanic flank collapse deposition and that they play a significant role in the shaping of the deposit's appearance.

Description: The Arctic Ocean plays an important role in the regulation of the earth's climate system, for instance by storing large amounts of carbon dioxide within its interior. It also plays a critical role in the global thermohaline circulation, transporting water entering from the Atlantic Ocean to the interior and initializing the southward transport of deep waters. Currently, the Arctic Ocean is undergoing rapid changes due to climate warming. The resulting consequences on ventilation patterns, however, are scarce. In this study we present transient tracer (CFC-12 and SF6) measurements, in conjunction with dissolved oxygen concentrations, to asses ventilation and circulation changes in the Eurasian Arctic Ocean over three decades (1991–2021). We constrained transit time distributions of water masses in different areas and quantified temporal variability in ventilation. Specifically, mean ages of intermediate water layers in the Eurasian Arctic Ocean were evaluated, revealing a decrease in ventilation in each of the designated areas from 2005 to 2021. This intermediate layer (250–1,500 m) is dominated by Atlantic Water entering from the Nordic Seas. We also identify a variability in ventilation during the observation period in most regions, as the data from 1991 shows mean ages comparable to those from 2021. Only in the northern Amundsen Basin, where the Arctic Ocean Boundary Current is present at intermediate depths, the ventilation in 1991 is congruent to the one in 2005, increasing thereafter until 2021. This suggests a reduced ventilation and decrease in the strength of the Boundary Current during the last 16 years. Key Points Temporal variability of ventilation in the Eurasian Arctic Ocean during the past 30 years is estimated by observations of transient tracers We found a slow down of the ventilation between 2005 and 2021 in the intermediate waters Evidence of multidecadal variability of ventilation in the intermediate waters of the Eurasian Arctic Ocean is present Plain Language Summary The Eurasian Arctic Ocean, the region of the Arctic Ocean connected to the European and Asian continents, is an important pathway for recently ventilated water from the Nordic Seas. These waters are exported back to the North Atlantic following their travel through the Arctic Ocean. Ventilation describes the process of surface waters being transported into the interior ocean due to increasing density, which affects the underlying water masses. In this study we investigate how the ventilation patterns have evolved in the Eurasian Arctic Ocean over the past three decades, using transient tracer (CFC-12 and SF6) measurements. We observed a significant change in the intermediate layer (250–1,500 m) with older waters found in measurements in 1991 and 2021 compared to 2005 and 2015. Moreover, our data suggest a slowdown in ventilation throughout the three decades in the northern Amundsen Basin, implying a decrease in the circulation time-scale of the Arctic Ocean Boundary Current over the past 16 years. This has potentially important implications for the transport of, for example, heat, salt or oxygen from the Atlantic Ocean around the Arctic Ocean, and back.

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

Description: 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: The South Atlantic Meridional Overturning Circulation initiative began as a grassroots effort to study the South Atlantic Ocean and its impact on climate. In striving towards this goal, it has also become a platform for the empowerment of women and international scientists.

Description: Phytoplankton forms the base of aquatic food webs and element cycling in diverse aquatic systems. The fate of phytoplankton-derived organic matter, however, often remains unresolved as it is controlled by complex, interlinked remineralization and sedimentation processes. We here investigate a rarely considered control mechanism on sinking organic matter fluxes: fungal parasites infecting phytoplankton. We demonstrate that bacterial colonization is promoted 3.5-fold on fungal-infected phytoplankton cells in comparison to non-infected cells in a cultured model pathosystem (diatom Synedra, fungal microparasite Zygophlyctis, and co-growing bacteria), and even ≥17-fold in field-sampled populations (Planktothrix, Synedra, and Fragilaria). Additional data obtained using the Synedra–Zygophlyctis model system reveals that fungal infections reduce the formation of aggregates. Moreover, carbon respiration is 2-fold higher and settling velocities are 11–48% lower for similar-sized fungal-infected vs. non-infected aggregates. Our data imply that parasites can effectively control the fate of phytoplankton-derived organic matter on a single-cell to single-aggregate scale, potentially enhancing remineralization and reducing sedimentation in freshwater and coastal systems.

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: In search for critical elements, polymetallic nodules at the deep abyssal seafloor are targeted for mining operations. Nodules efficiently scavenge and retain several naturally occurring uranium-series radioisotopes, which predominantly emit alpha radiation during decay. Here, we present new data on the activity concentrations of thorium-230, radium-226, and protactinium-231, as well as on the release of radon-222 in and from nodules from the NE Pacific Ocean. In line with abundantly published data from historic studies, we demonstrate that the activity concentrations for several alpha emitters are often higher than 5 Bq g −1 at the surface of the nodules. These observed values can exceed current exemption levels by up to a factor of 1000, and even entire nodules commonly exceed these limits. Exemption levels are in place for naturally occurring radioactive materials (NORM) such as ores and slags, to protect the public and to ensure occupational health and radiation safety. In this context, we discuss three ways of radiation exposure from nodules, including the inhalation or ingestion of nodule fines, the inhalation of radon gas in enclosed spaces and the potential concentration of some radioisotopes during nodule processing. Seen in this light, inappropriate handling of polymetallic nodules poses serious health risks.

Description: Recent studies indicate that mantle plumes, which transfer material and heat from the earth’s interior to its surface, represent multifaceted upwellings. The Tristan-Gough hotspot track (South Atlantic), which formed above a mantle plume, documents spatial geochemical zonation in two distinct sub-tracks since ~70 Ma. The origin and the sudden appearance of two distinct geochemical flavors is enigmatic, but could provide insights into the structural evolution of mantle plumes. Sr–Nd–Pb–Hf isotope data from the Late Cretaceous Rio Grande Rise and adjacent Jean Charcot Seamount Chain (South American Plate), which represent the counterpart of the older Tristan-Gough volcanic track (African Plate), extends the bilateral-zonation to ~100 Ma. Our results support recent numerical models, demonstrating that mantle plumes can split into distinct upper mantle conduits, and provide evidence that these plumelets formed at the plume head-to-plume tail transition. We attribute the plume zonation to sampling the geochemically-graded margin of the African Large Low-Shear-Velocity Province.

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: Changes in the Atlantic Meridional Overturning Circulation (AMOC) represent a crucial component of Northern Hemisphere climate variability. In modelling studies decadal overturning variability has been attributed to the intensity of deep winter convection in the Labrador Sea. This linkage is challenged by transport observations at sections across the subpolar gyre. Here we report simulations with an eddy-rich ocean model which captures the observed concentration of downwelling in the northeastern Atlantic and the negligible impact of interannual variations in Labrador Sea convection during the last decade. However, the exceptionally cold winters in the Labrador Sea during the first half of the 1990s induced a positive AMOC anomaly of more than 20%, mainly by augmenting the downwelling in the northeastern North Atlantic. The remote effect of excessive Labrador Sea buoyancy forcing is related to rapid spreading of mid-depth density anomalies into the Irminger Sea and their entrainment into the deep boundary current off Greenland.

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: 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: Global warming may alter the dynamics of infectious diseases by affecting important steps in the transmission of pathogens and parasites. In trematode parasites, the emergence of cercarial stages from their hosts is temperature-dependent, being highest around a thermal optimum. If environmental temperatures exceed this optimum as a consequence of global warming, this may affect cercarial transmission. However, our knowledge of cercarial emergence patterns of species from high temperature environments is currently very limited. Here, we investigated the effect of temperature on the emergence of two common trematode species from an abundant mud snail Pirenella cingulata in the Persian Gulf, the warmest sea on Earth. Infected snails were incubated in the laboratory at 6 temperatures from 10 to 40°C for 3 days. We found an optimal temperature for cercarial emergence of 32.0°C and 33.5°C for Acanthotrema tridactyla and Cyathocotylidae gen. sp., respectively, which are the warmest recorded thermal optima for any aquatic trematode species. Emergence of both species dropped at 40°C, suggesting upper thermal limits to emergence. Overall, Persian Gulf trematodes may be among the most heat-tolerant marine trematode species, indicating a potential for dispersing to regions that will continue to warm in the future.

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: Knowledge of the characteristics of natural climate variability is vital when assessing the range of plausible future climate trajectories in the next decades to centuries. The reliable detection of climate fluctuations on multidecadal to centennial timescales depends on proxy reconstructions and model simulations, as the instrumental record extends back only a few decades in most parts of the world. Systematic comparisons between model-simulated and proxy-based inferences of natural variability, however, often seem contradictory. Locally, simulated temperature variability is consistently smaller on multidecadal and longer timescales than is indicated by proxy-based reconstructions, implying that climate models or proxy interpretations might have deficiencies. In contrast, at global scales, studies found agreement between simulated and proxy reconstructed temperature variations. Here we review the evidence regarding the scale of natural temperature variability during recent millennia. We identify systematic reconstruction deficiencies that may contribute to differing local and global model–proxy agreement but conclude that they are probably insufficient to resolve such discrepancies. Instead, we argue that regional climate variations persisted for longer timescales than climate models simulating past climate states are able to reproduce. This would imply an underestimation of the regional variability on multidecadal and longer timescales and would bias climate projections and attribution studies. Thus, efforts are needed to improve the simulation of natural variability in climate models accompanied by further refining proxy-based inferences of variability.

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: Abyssal seafloor communities cover more than 60% of Earth’s surface. Despite their great size, abyssal plains extend across modest environmental gradients compared to other marine ecosystems. However, little is known about the patterns and processes regulating biodiversity or potentially delimiting biogeographical boundaries at regional scales in the abyss. Improved macroecological understanding of remote abyssal environments is urgent as threats of widespread anthropogenic disturbance grow in the deep ocean. Here, we use a new, basin-scale dataset to show the existence of clear regional zonation in abyssal communities across the 5,000 km span of the Clarion–Clipperton Zone (northeast Pacific), an area targeted for deep-sea mining. We found two pronounced biogeographic provinces, deep and shallow-abyssal, separated by a transition zone between 4,300 and 4,800 m depth. Surprisingly, species richness was maintained across this boundary by phylum-level taxonomic replacements. These regional transitions are probably related to calcium carbonate saturation boundaries as taxa dependent on calcium carbonate structures, such as shelled molluscs, appear restricted to the shallower province. Our results suggest geochemical and climatic forcing on distributions of abyssal populations over large spatial scales and provide a potential paradigm for deep-sea macroecology, opening a new basis for regional-scale biodiversity research and conservation strategies in Earth’s largest biome.

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: The Palaeocene–Eocene Thermal Maximum (PETM) was a global warming event of 5–6 °C around 56 million years ago caused by input of carbon into the ocean and atmosphere. Hydrothermal venting of greenhouse gases produced in contact aureoles surrounding magmatic intrusions in the North Atlantic Igneous Province have been proposed to play a key role in the PETM carbon-cycle perturbation, but the precise timing, magnitude and climatic impact of such venting remains uncertain. Here we present seismic data and the results of a five-borehole transect sampling the crater of a hydrothermal vent complex in the Northeast Atlantic. Stable carbon isotope stratigraphy and dinoflagellate cyst biostratigraphy reveal a negative carbon isotope excursion coincident with the appearance of the index taxon Apectodinium augustum in the vent crater, firmly tying the infill to the PETM. The shape of the crater and stratified sediments suggests large-scale explosive gas release during the initial phase of vent formation followed by rapid, but largely undisturbed, diatomite-rich infill. Moreover, we show that these vents erupted in very shallow water across the North Atlantic Igneous Province, such that volatile emissions would have entered the atmosphere almost directly without oxidation to CO 2 and at the onset of the PETM.

Description: Using population genomics, we reveal the global dispersal history of eelgrass (Zostera marina) from its origin in the Northwest Pacific Ocean, across to the East Pacific, Atlantic and the Mediterranean Sea. Striking differences in genetic diversity among the locations reflect past glaciations and repeated bottlenecks during Z. marina’s worldwide colonization.

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: Armed conflicts have, in addition to severe impacts on human lives and infrastructure, also impacts on the environment, which needs to be assessed and documented. On September the 26th 2022, unknown perpetrators deliberately ruptured the two gas pipelines Nord Stream 1 and 2 with four coordinated explosions near a major chemical munition dump site near the Danish island of Bornholm in the Baltic Sea. While the massive release of natural gas into atmosphere raised serious concerns concerning the contribution to climate change—this paper assesses the overlooked direct impact of the explosions on the marine ecosystem. Seals and porpoises within a radius of four km would be at high risk of being killed by the shockwave, while temporary impact on hearing would be expected up to 50 km away. As the Baltic Proper population of harbour porpoises ( Phocoena phocoena ) is critically endangered, the loss or serious injury of even a single individual is considered a significant impact on the population. The rupture moreover resulted in the resuspension of 250000 metric tons of heavily contaminated sediment from deep-sea sedimentary basin for over a week, resulting in unacceptable toxicological risks towards fish and other biota in 11 km 3 water in the area for more than a month.

Description: The presence of anti-nutritive compounds like glucosinolates (GSLs) in the rapeseed meal severely restricts its utilization as animal feed. Therefore, reducing the GSL content to 〈 18 µmol/g dry weight in the seeds is a major breeding target. While candidate genes involved in the biosynthesis of GSLs have been described in rapeseed, comprehensive functional analyses are missing. By knocking out the aliphatic GSL biosynthesis genes BnMYB28 and BnCYP79F1 encoding an R2R3 MYB transcription factor and a cytochrome P450 enzyme, respectively, we aimed to reduce the seed GSL content in rapeseed. After expression analyses on single paralogs, we used an ethyl methanesulfonate (EMS) treated population of the inbred winter rapeseed 'Express617' to detect functional mutations in the two gene families. Our results provide the first functional analysis by knock-out for the two GSL biosynthesis genes in winter rapeseed. We demonstrate that independent knock-out mutants of the two genes possessed significantly reduced seed aliphatic GSLs, primarily progoitrin. Compared to the wildtype Express617 control plants (36.3 µmol/g DW), progoitrin levels were decreased by 55.3% and 32.4% in functional mutants of BnMYB28 (16.20 µmol/g DW) and BnCYP79F1 (24.5 µmol/g DW), respectively. Our study provides a strong basis for breeding rapeseed with improved meal quality in the future.

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