ALBERT

Description: Ice calved from the Antarctic and Greenland Ice Sheets or tidewater glaciers ultimately melts in the ocean contributing to sea-level rise. Icebergs have also been described as biological hotspots due to their potential roles as platforms for marine mammals and birds, and as micronutrient fertilizing agents. Icebergs may be especially important in the Southern Ocean where availability of the micronutrients iron and manganese extensively limits marine primary production. Whilst icebergs have long been described as a source of iron to the ocean, their nutrient signature is poorly constrained and it is unclear if there are regional differences. Here we show that 589 ice fragments collected from floating ice in contrasting regions spanning the Antarctic Peninsula, Greenland, and smaller tidewater systems in Svalbard, Patagonia and Iceland have similar characteristic (micro)nutrient signatures with limited or no significant differences between regions. Icebergs are a minor or negligible source of macronutrients to the ocean with low concentrations of NOx (NO3 + NO2, median 0.51 µM), PO4 (median 0.04 µM), and dissolved Si (dSi, median 0.02 µM). In contrast, icebergs deliver elevated concentrations of dissolved Fe (dFe; mean 82 nM, median 12 nM) and Mn (dMn; mean 26 nM, median 2.6 nM). A tight correlation between total dissolvable Fe and Mn (R2 = 0.95) and a Mn:Fe ratio of 0.024 suggested a lithogenic origin for the majority of sediment present in ice. Total dissolvable Fe and Mn retained a strong relationship with sediment load (both R2 = 0.43, p〈0.001), whereas weaker relationships were observed for dFe, dMn and dSi. Sediment load for Antarctic ice (median 9 mg L-1, n=144) was low compared to prior reported values for the Arctic. A particularly curious incidental finding was that melting samples of ice were observed to rapidly lose their sediment load, even when sediment layers were embedded within the ice and stored in the dark. Our results demonstrated that the nutrient signature of icebergs is consistent with an atmospheric source of NOx and PO4. Conversely, high Fe and Mn, and modest dSi concentrations, are associated with englacial sediment, which experiences limited biogeochemical processing prior to release into the ocean.

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

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

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

Description: Riverine nutrient export is an important process in marine coastal biogeochemistry and also impacts global marine biology. The nitrogen cycle is a key player here. Internal feedbacks regulate not only nitrogen distribution, but also primary production and thereby oxygen concentrations. Phosphorus is another essential nutrient and interacts with the nitrogen cycle via different feedback mechanisms. After a previous study of the marine nitrogen cycle response to riverine nitrogen supply, we here additionally include phosphorus from river export with different phosphorus burial scenarios and study the impact of phosphorus alone and in combination with nitrogen in a global 3-D ocean biogeochemistry model. Again, we analyse the effects on near coastal and open ocean biogeochemistry. We find that the addition of bio-available riverine phosphorus alone or together with nitrogen affects marine biology on millennial timescales more than riverine nitrogen alone. Biogeochemical feedbacks in the marine nitrogen cycle are strongly influenced by the additional phosphorus. Where bio-available phosphorus is increased by river input, nitrogen concentrations increase as well, except for regions with high denitrification rates. High phosphorus burial rates decrease biological production significantly. Globally, riverine phosphorus leads to elevated primary production rates in the coastal and open oceans.

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

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

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

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

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

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

Description: Ocean alkalinity enhancement (OAE) stands as a promising carbon dioxide removal technology. Yet, this solution to climate change entails shifts in water chemistry with unknown consequences for marine fish that are critical to ecosystem health and food security. With a laboratory and mesocosm experiment, we show that early life stages of fish can be resistant to OAE. We examined metabolic rate, swimming behavior, growth and survival in Atlantic herring (Clupea harengus) and other temperate coastal fish species. Neither direct physiological nor indirect food web-mediated impacts of OAE were apparent. This was despite non-CO2-equilibrated OAE (ΔTA = +600 µmol kg-1) that induces strong perturbations (ΔpH = +0.7, pCO2 = 75 µatm) compared to alternative deployment scenarios. Whilst our results give cause for optimism regarding the large-scale application of OAE, other life history stages (embryos) and habitats (open ocean) may prove more vulnerable. Still, our study across ecological scales (organism to community) and exposure times (short- to long-term) suggests that some fish populations, including key fisheries species, may be resilient to the carbonate chemistry changes under OAE.

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: Ocean alkalinity enhancement (OAE) is considered for the long-term removal of gigatons of carbon dioxide (CO2) from the atmosphere to achieve our climate goals. Little is known, however, about the ecosystem-level changes in biogeochemical functioning that may result from the chemical sequestration of CO2 in seawater, and how stable the sequestration is. We studied these two aspects in natural plankton communities under carbonate-based, CO2-equilibrated OAE in the nutrient-poor North Atlantic. During a month-long mesocosm experiment, the majority of biogeochemical pools, including inorganic nutrients, particulate organic carbon and phosphorus as well as biogenic silica, remained unaltered across all OAE levels of up to a doubling of ambient alkalinity (+2400 µeq kg-1). Noticeable exceptions were a minor decrease in particulate organic nitrogen and an increase in the carbon to nitrogen ratio (C:N) of particulate organic matter in response to OAE. Thus, in our nitrogen limited system, nitrogen turnover processes appear more susceptible than those of other elements leading to decreased food quality and increased organic carbon storage. However, alkalinity and chemical CO2 sequestration were not stable at all levels of OAE. Two weeks after alkalinity addition, we measured a loss of added alkalinity and of the initially stored CO2 in the mesocosm where alkalinity was highest (+2400 µeq kg-1, Ωaragonite ~10). The loss rate accelerated over time. Additional tests showed that such secondary precipitation can be initiated by particles acting as precipitation nuclei and that this process can occur even at lower levels of OAE. In conclusion, on the one hand, our study under carbonate-based OAE where the carbon is already sequestered, the risk of major and sustained impacts on biogeochemical functioning may be low in the nutrient-poor ocean. On the other hand, the durability of carbon sequestration using OAE could be constrained by alkalinity loss in supersaturated waters with precipitation nuclei present. Our study provides evaluation of ecosystem impacts of an idealised OAE deployment for monitoring, reporting and verification (MRV) in an oligotrophic system. Whether biogeochemical functioning is resilient to more technically simple and economically more viable approaches that induce stronger water chemistry perturbations remains to be seen.

Description: Circulation anomalies accompanying Sudden Stratospheric Warmings (SSWs) can have a significant impact on the troposphere. This surface response is observed for some but not all SSWs, and their downward coupling is not fully understood. We use an existing classification method to separate downward- and non-propagating SSWs (d/nSSWs) in ERA5 reanalysis data for the years 1979–2019. The differences in SSW downward propagation in composites of spatial patterns clearly show that dSSWs dominate the surface regional impacts following SSWs. During dSSWs, the upper-tropospheric jet stream is significantly displaced equatorward. Wave activity analysis shows remarkable differences between d/nSSWs for planetary and synoptic-scale waves. Enhanced stratospheric planetary eddy kinetic energy (EKE) and heat fluxes around the central date of dSSWs are followed by increased synoptic-scale wave activity and even surface coupling for synoptic-scale EKE. An observed significant reduction in upper-tropospheric synoptic-scale momentum fluxes following dSSWs confirms the important role of tropospheric eddy feedbacks for coupling to the surface. Our findings emphasize the role of the lower stratosphere and synoptic-scale waves in coupling the SSW signal to the surface and agree with mechanisms suggested in earlier modeling studies.

Description: Ocean alkalinity enhancement (OAE) is considered one of the most promising approaches to actively remove carbon dioxide (CO2) from the atmosphere by accelerating the natural process of rock weathering. This approach involves introducing alkaline substances sourced from natural mineral deposits such as olivine, basalt, and carbonates or obtained from industrial waste products such as steel slags, into seawater and dispersing them over coastal areas. Some of these natural and industrial substances contain trace metals, which would be released into the oceans along with the alkalinity enhancement. The trace metals could serve as micronutrients for marine organisms at low concentrations, but could potentially become toxic at high concentrations, adversely affecting marine biota. To comprehensively assess the feasibility of OAE, it is crucial to understand how the phytoplankton, which forms the base of marine food webs, responds to ocean alkalinization and associated trace metal perturbations. In this study, we investigated the toxicity of nickel on three representative phytoplankton species across a range of Ni concentrations (from 0 to 100 µmol L-1 with 12 µmol L-1 synthetic organic ligand). The results showed that the growth of the tested species was impacted differently. The low growth inhibition and high IC50 (concentration to inhibit growth rate by 50 %) revealed that both the coccolithophore Emiliania huxleyi and the dinoflagellate Amphidinium carterae were mildly impacted by the increase in Ni concentrations while the rapid response to exposure of Ni, high growth rate inhibition, and low IC50 of Thalassiosira weissflogii indicate low tolerance to Ni in this species. In conclusion, the variability in phytoplankton sensitivity to Ni suggests that for OAE applications with Ni-rich materials caution is required and critical toxic thresholds for Ni must be avoided.

Description: The central Arctic Ocean (CAO) plays an important role in the global carbon cycle, but the current and future exchange of the climate-forcing trace gases methane (CH4) and carbon dioxide (CO2) between the CAO and the atmosphere is highly uncertain. In particular, there are very few observations of near-surface gas concentrations or direct air–sea CO2 flux estimates and no previously reported direct air–sea CH4 flux estimates from the CAO. Furthermore, the effect of sea ice on the exchange is not well understood. We present direct measurements of the air–sea flux of CH4 and CO2, as well as air–snow fluxes of CO2 in the summertime CAO north of 82.5∘ N from the Synoptic Arctic Survey (SAS) expedition carried out on the Swedish icebreaker Oden in 2021. Measurements of air–sea CH4 and CO2 flux were made using floating chambers deployed in leads accessed from sea ice and from the side of Oden, and air–snow fluxes were determined from chambers deployed on sea ice. Gas transfer velocities determined from fluxes and surface-water-dissolved gas concentrations exhibited a weaker wind speed dependence than existing parameterisations, with a median sea-ice lead gas transfer rate of 2.5 cm h−1 applicable over the observed 10 m wind speed range (1–11 m s−1). The average observed air–sea CO2 flux was −7.6 ..., and the average air–snow CO2 flux was −1.1 . Extrapolating these fluxes and the corresponding sea-ice concentrations gives an August and September flux for the CAO of −1.75 ... , within the range of previous indirect estimates. The average observed air–sea CH4 flux of 3.5 ..., accounting for sea-ice concentration, equates to an August and September CAO flux of 0.35 , lower than previous estimates and implying that the CAO is a very small (≪ 1 %) contributor to the Arctic flux of CH4 to the atmosphere.

Description: Identification of seismically active fault zones and the definition of sufficiently large respect distances from these faults which enable avoiding the damaged rock zone surrounding the ruptured ground commonly are amongst the first steps to take in the geoscientific evaluation of sites suitable for nuclear waste disposal. In this work we present a GIS-based approach, using the earthquake-epicentre locations from the instrumental earthquake record of South-Korea to identify potentially active fault zones in the country, and compare different strategies for fault zone buffer creation as originally developed for site search in the high seismicity country Japan, and the low-to-moderate seismicity countries Germany and Sweden. In order to characterize the hazard potential of the Korean fault zones, we moreover conducted slip tendency analysis, here for the first time covering the fault zones of the entire Korean Peninsula. For our analyses we used the geo-spatial information from a new version of the Geological map of South-Korea, containing the outlines of 11 rock units, which we simplified to distinguish between 4 different rock types (granites, metamorphic rocks, sedimentary rocks and igneous rocks) and the surface traces of 1,528 fault zones and 6,654 lineaments identified through years of field work and data processing, a rich geo-dataset which we will publish along with this manuscript. Our approach for identification of active fault zones was developed without prior knowledge of already known seismically active fault zones, and as a proof of concept the results later were compared to a map containing already identified active fault zones. The comparison revealed that our approach identified 16 of the 21 known seismically active faults and added 472 previously unknown potentially active faults. The 5 seismically active fault zones which were not identified by our approach are located in the NE- and SW-sectors of the Korean Peninsula, which haven’t seen much recent seismic activity, and thus are not sufficiently well covered by the seismic record. The strike directions of fault zones identified as active are in good agreement with the orientation of the current stress field of the peninsula and slip tendency analysis provided first insights into subsurface geometry such as the dip angles of both active and inactive fault zones. The results of our work are of major importance for the early-stage seismic hazard assessment that has to be conducted in support of the nuclear waste disposal siting in South-Korea. Moreover, the GIS-based methods for identification of active fault zones and buffering of respect areas around fault zone traces presented here, are applicable also elsewhere.

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: Total alkalinity (AT) and dissolved inorganic carbon (CT) in the oceans are important properties with respect to understanding the ocean carbon cycle and its link to global change (ocean carbon sinks and sources, ocean acidification) and ultimately finding carbon-based solutions or mitigation procedures (marine carbon removal). We present a database of more than 44 400 AT and CT observations along with basic ancillary data (spatiotemporal location, depth, temperature and salinity) from various ocean regions obtained, mainly in the framework of French projects, since 1993. This includes both surface and water column data acquired in the open ocean, coastal zones and in the Mediterranean Sea and either from time series or dedicated one-off cruises. Most AT and CT data in this synthesis were measured from discrete samples using the same closed-cell potentiometric titration calibrated with Certified Reference Material, with an overall accuracy of ±4 µmol kg−1 for both AT and CT. The data are provided in two separate datasets – for the Global Ocean and the Mediterranean Sea (https://doi.org/10.17882/95414, Metzl et al., 2023), respectively – that offer a direct use for regional or global purposes, e.g., AT–salinity relationships, long-term CT estimates, and constraint and validation of diagnostic CT and AT reconstructed fields or ocean carbon and coupled climate–carbon models simulations as well as data derived from Biogeochemical-Argo (BGC-Argo) floats. When associated with other properties, these data can also be used to calculate pH, the fugacity of CO2 (fCO2) and other carbon system properties to derive ocean acidification rates or air–sea CO2 fluxes.

Description: The upper wind-driven circulation in the tropical Atlantic Ocean plays a key role in the basin-wide distribution of water mass properties and affects the transport of heat, freshwater, and biogeochemical tracers such as oxygen or nutrients. It is crucial to improve our understanding of its long-term behaviour, which largely relies on model simulations and applied forcing due to sparse observational data coverage, especially before the mid-2000s. Here, we apply two different forcing products, the Coordinated Ocean-ice Reference Experiments (CORE) v2 and the Japanese 55-year Reanalysis (JRA55-do) surface dataset, to a high-resolution ocean model. Where possible, we compare the simulated results to long-term observations. We find large discrepancies between the two simulations regarding the wind and current field. In the CORE simulation, strong, large-scale wind stress curl amplitudes above the upwelling regions of the eastern tropical North Atlantic seem to cause an overestimation of the mean and seasonal variability in the eastward subsurface current just north of the Equator. The wind stress curl of JRA55-do forcing shows much finer structures, and the JRA55-do simulation is in better agreement with the mean and intraseasonal fluctuations in the subsurface current found in observations. The northern branch of the South Equatorial Current flows westward at the surface just north of the Equator. On interannual to decadal timescales, it shows a high correlation of R=0.9 with the zonal wind stress in the CORE simulation but only a weak correlation of R=0.35 in the JRA55-do simulation. We also identify similarities between the two simulations. The strength of the eastward-flowing North Equatorial Counter Current located between 3 and 10° N covaries with the strength of the meridional wind stress just north of the Equator on interannual to decadal timescales in the two simulations. Both simulations present a comparable mean, seasonal cycle and trend of the eastward off-equatorial subsurface current south of the Equator but underestimate the current strength by half compared to observations. In both simulations, the eastward-flowing Equatorial Undercurrent weakened between 1990 and 2009. In the JRA simulation, which covers the modern period of observations, the Equatorial Undercurrent strengthened again between 2008 to 2018, which agrees with observations, although the simulation underestimates the strengthening by over a third. We propose that long-term observations, once they have reached a critical length, need to be used to test the quality of wind-driven simulations. This study presents one step in this direction.

Description: We conducted extensive sediment trap experiments in the Benguela Upwelling System (BUS) in the south-eastern Atlantic Ocean to study the influence of zooplankton on the flux of particulate organic carbon (POC) through the water column and its sedimentation. Two long term moored and sixteen short term free-floating sediment trap systems were deployed. The mooring experiments were conducted for several years and the sixteen drifters were deployed on three different research cruises between 2019 and 2021. Zooplankton was separated from the trapped material and divided into 8 different zooplankton groups. In contrast to zooplankton which actively carries POC into the traps in the form of biomass (active POC flux), the remaining fraction of the trapped material was assumed to fall passively into the traps along with sinking particles (passive POC flux). The results show, in line with other studies, that copepods dominate the active POC flux, with the active POC flux in the southern BUS (sBUS) being about three times higher than in the northern BUS (nBUS). In contrast, the differences between the passive POC fluxes in the nBUS and sBUS were small. Despite large variations, which reflected the variability within the two subsystems, the mean passive POC fluxes from the drifters and the moored traps could be described using a common POC flux attenuation equation. However, the almost equal passive POC flux, on the one hand, and large variations in the POC concentration in the surface sediments between the nBUS and sBUS, on the other hand, imply that factors others than the POC supply exert the main control on POC sedimentation in the BUS. The varying intensity of the near-bottom oxygen minimum zone (OMZ), which is more pronounced in the nBUS than in the sBUS, could in turn explain the differences in the sediments, as the lack of oxygen reduces the POC degradation. Hence, globally expanding OMZs might favour POC sedimentation in regions formerly exposed to oxygenated bottom water but bear the risk of increasing the frequency of anoxic events in the oxygen-poor upwelling systems. Apart from associated release of CH4, which is a much more potent greenhouse gas than CO2, such events pose a major threat to the pelagic ecosystem and fisheries.

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: Here we present a confocal Fe K-edge μ-XANES method (where XANES stands for X-ray absorption near-edge spectroscopy) for the analysis of Fe oxidation state in heterogeneous and one-side-polished samples. The new technique allows for an analysis of small volumes with high spatial 3D resolution of 〈100 µm3. The probed volume is restricted to that just beneath the surface of the exposed object. This protocol avoids contamination of the signal by the host material and minimizes self-absorption effects. This technique has been tested on a set of experimental glasses with a wide range of Fe3+  ΣFe ratios. The method was applied to the analysis of natural melt inclusions trapped in forsteritic to fayalitic olivine crystals of the Hekla volcano, Iceland. Our measurements reveal changes in Fe3+  ΣFe from 0.17 in basaltic up to 0.45 in dacitic melts, whereas the magnetite–ilmenite equilibrium shows redox conditions with Fe3+  ΣFe ≤0.20 (close to FMQ, fayalite–magnetite–quartz redox equilibrium) along the entire range of Hekla melt compositions. This discrepancy indicates that the oxidized nature of glasses in the melt inclusions could be related to the post-entrapment process of diffusive hydrogen loss from inclusions and associated oxidation of Fe in the melt. The Fe3+  ΣFe ratio in silicic melts is particularly susceptible to this process due to their low FeO content, and it should be critically evaluated before petrological interpretation.

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: We examine the impact of horizontal resolution and model time step on the climate of the OpenIFS version 43r3 atmospheric general circulation model. A series of simulations for the period 1979–2019 are conducted with various horizontal resolutions (i.e. ∼100, ∼50, and ∼25 km) while maintaining the same time step (i.e. 15 min) and using different time steps (i.e. 60, 30, and 15 min) at 100 km horizontal resolution. We find that the surface zonal wind bias is significantly reduced over certain regions such as the Southern Ocean and the Northern Hemisphere mid-latitudes and in tropical and subtropical regions at a high horizontal resolution (i.e. ∼25 km). Similar improvement is evident too when using a coarse-resolution model (∼100 km) with a smaller time step (i.e. 30 and 15 min). We also find improvements in Rossby wave amplitude and phase speed, as well as in weather regime patterns, when a smaller time step or higher horizontal resolution is used. The improvement in the wind bias when using the shorter time step is mostly due to an increase in shallow and mid-level convection that enhances vertical mixing in the lower troposphere. The enhanced mixing allows frictional effects to influence a deeper layer and reduces wind and wind speed throughout the troposphere. However, precipitation biases generally increase with higher horizontal resolutions or smaller time steps, whereas the surface air temperature bias exhibits a small improvement over North America and the eastern Eurasian continent. We argue that the bias improvement in the highest-horizontal-resolution (i.e. ∼25 km) configuration benefits from a combination of both the enhanced horizontal resolution and the shorter time step. In summary, we demonstrate that, by reducing the time step in the coarse-resolution (∼100 km) OpenIFS model, one can alleviate some climate biases at a lower cost than by increasing the horizontal resolution.

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

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

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

Description: Phytoplankton forms the base of the marine food web by transforming CO2 into organic carbon via photosynthesis. Some of the organic carbon is then transferred through the food web and exported into the deep ocean, a process known as the biological carbon pump. Despite the importance of phytoplankton for marine ecosystems and the global carbon cycle, projections of phytoplankton biomass in response to climate change differ strongly across Earth system models, illustrating uncertainty in our understanding of the underlying processes. Differences are especially large in the Southern Ocean, a region that is notoriously difficult to represent in models. Here, we argue that water column-integrated phytoplankton biomass in the Southern Ocean is projected to largely remain unchanged under climate change by the CMIP6 multi-model ensemble because of a shifting balance of bottom-up and top-down processes driven by a shoaling mixed layer depth. A shallower mixed layer is projected to improve growth conditions and consequently weaken bottom-up control. In addition to enhanced phytoplankton growth, the shoaling of the mixed layer also compresses phytoplankton closer to the surface and promotes zooplankton grazing efficiency, thus intensifying top-down control. Overall, our results suggest that while changes in bottom-up conditions stimulate enhanced growth, intensified top-down control opposes an increase in phytoplankton and becomes increasingly important for phytoplankton response under climate change in the Southern Ocean.

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

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

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

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

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

Description: The climate science community aims to improve our understanding of climate change due to anthropogenic influences on atmospheric composition and the Earth's surface. Yet not all climate interactions are fully understood and diversity in climate model experiments persists as assessed in the latest Intergovernmental Panel on Climate Change (IPCC) assessment report. This article synthesizes current challenges and emphasizes opportunities for advancing our understanding of climate change and model diversity. The perspective of this article is based on expert views from three multi-model intercomparison projects (MIPs) – the Precipitation Driver Response MIP (PDRMIP), the Aerosol and Chemistry MIP (AerChemMIP), and the Radiative Forcing MIP (RFMIP). While there are many shared interests and specialisms across the MIPs, they have their own scientific foci and specific approaches. The partial overlap between the MIPs proved useful for advancing the understanding of the perturbation-response paradigm through multi-model ensembles of Earth System Models of varying complexity. It specifically facilitated contributions to the research field through sharing knowledge on best practices for the design of model diagnostics and experimental strategies across MIP boundaries, e.g., for estimating effective radiative forcing. We discuss the challenges of gaining insights from highly complex models that have specific biases and provide guidance from our lessons learned. Promising ideas to overcome some long-standing challenges in the near future are kilometer-scale experiments to better simulate circulation-dependent processes where it is possible, and machine learning approaches for faster and better sub-grid scale parameterizations where they are needed. Both would improve our ability to adopt a smart experimental design with an optimal tradeoff between resolution, complexity and simulation length. Future experiments can be evaluated and improved with sophisticated methods that leverage multiple observational datasets, and thereby, help to advance the understanding of climate change and its impacts.

Description: Nitrogen (N) is a crucial limiting nutrient for phytoplankton growth in the ocean. The main source of bioavailable N in the ocean is delivered by N2-fixing diazotrophs in the surface layer. Since field observation of N2 fixation are spatially and temporally sparse, the fundamental processes and mechanisms controlling N2 fixation are not well understood and constrained. Here, we implement benthic denitrification in an Earth System Model of intermediate complexity (UVic-ESCM 2.9) coupled to an optimality-based plankton ecosystem model (OPEM v1.1). Benthic denitrification occurs mostly in coastal upwelling regions and on shallow continental shelves, and is the largest N-loss process in the global ocean. We calibrate our model against three different combinations of observed Chl, NO3-, PO43-, O2 and N* = NO3- −16PO43- +2.9. The inclusion of N* provides a powerful constraint on biogeochemical model behavior. Our new model version including benthic denitrification simulates higher global rates of N2 fixation with a more realistic distribution extending to higher latitudes that are supported by independent estimates based on geochemical data. Oxygen deficient zone volume and water column denitrification rates are reduced in the new version, indicating that including benthic denitrification may improve global biogeochemical models that commonly overestimate anoxic zones. With the improved representation of the ocean N cycle, our new model configuration also yields better global net primary production (NPP) when compared to the independent datasets not included in the calibration. Benthic denitrification plays an important role shaping N2 fixation and NPP throughout the global ocean in our model, and should be considered when evaluating and predicting their response to environmental change.

Description: Since a pH sensor has become available that is principally suitable for use on demanding autonomous measurement platforms, the marine CO2 system can be observed independently and continuously by Biogeochemical Argo floats. This opens the potential to detect variability and long-term changes in interior ocean inorganic carbon storage and quantify the ocean sink for atmospheric CO2. In combination with a second parameter of the marine CO2 system, pH can be a useful tool to derive the surface ocean CO2 partial pressure (pCO2). The large spatiotemporal variability in the marine CO2 system requires sustained observations to decipher trends and study the impacts of short-term events (e.g., eddies, storms, phytoplankton blooms) but also puts a high emphasis on the quality control of float-based pH measurements. In consequence, a consistent and rigorous quality control procedure is being established to correct sensor offsets or drifts as the interpretation of changes depends on accurate data. By applying current standardized routines of the Argo data management to pH measurements from a pH / O2 float pilot array in the subpolar North Atlantic Ocean, we assess the uncertainties and lack of objective criteria associated with the standardized routines, notably the choice of the reference method for the pH correction (CANYON-B, LIR-pH, ESPER-NN, and ESPER-LIR) and the reference depth for this adjustment. For the studied float array, significant differences ranging between ca. 0.003 pH units and ca. 0.04 pH units are observed between the four reference methods which have been proposed to correct float pH data. Through comparison against discrete and underway pH data from other platforms, an assessment of the adjusted float pH data quality is presented. The results point out noticeable discrepancies near the surface of 〉 0.004 pH units. In the context of converting surface ocean pH measurements into pCO2 data for the purpose of deriving air–sea CO2 fluxes, we conclude that an accuracy requirement of 0.01 pH units (equivalent to a pCO2 accuracy of 10 µatm as a minimum requirement for potential future inclusion in the Surface Ocean CO2 Atlas, SOCAT, database) is not systematically achieved in the upper ocean. While the limited dataset and regional focus of our study do not allow for firm conclusions, the evidence presented still calls for the inclusion of an additional independent pH reference in the surface ocean in the quality control routines. We therefore propose a way forward to enhance the float pH quality control procedure. In our analysis, the current philosophy of pH data correction against climatological reference data at one single depth in the deep ocean appears insufficient to assure adequate data quality in the surface ocean. Ideally, an additional reference point should be taken at or near the surface where the resulting pCO2 data are of the highest importance to monitor the air–sea exchange of CO2 and would have the potential to very significantly augment the impact of the current observation network.

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: A natural plankton community from oligotrophic subtropical waters of the Atlantic near Gran Canaria, Spain, was subjected to varying degrees of ocean alkalinity enhancement (OAE) to assess the potential physiological effects, in the context of the application of ocean carbon dioxide removal (CDR) techniques. We employed 8.3 m3 mesocosms with a sediment trap attached to the bottom, creating a gradient in total alkalinity (TA). The lowest point on this gradient was 2400 μmol · L-1, which corresponded to the natural alkalinity of the environment, and the highest point was 4800 μmol · L-1. Over the course of the 33-day experiment, the plankton community exhibited two distinct phases. In phase-I (days 5–20), a notable decline in the photosynthetic efficiency (Fv/Fm) was observed. This change was accompanied by substantial reductions in the abundances of picoeukaryotes, small size nanoeukaryotes (nanoeukaryotes-1), and microplankton. The cell viability of picoeukaryotes, as indicated by fluorescein-di-acetate hydrolysis by cellular esterases (FDA- green fluorescence), slightly increased by the end of phase-I whilst the viability of nanoeukaryotes 1 and Synechococcus spp . did not change. Reactive oxygen species levels (ROS-green fluorescence) showed no significant changes for any of the functional groups. In contrast, in phase-II (days 21–33), a pronounced community response was observed. Increases in Fv/Fm in the intermediate OAE treatments of ∆900 to ∆1800 μmol · L-1 and in chlorophyll-a (Chl-a), chlorophyll-c2 (Chl-c2) , fucoxanthin and divinyl-Chl-a were attributed to the emergence of blooms of large size nanoeukaryotes (nanoeukaryotes-2) from the genera Chrysochromulina, as well as picoeukaryotes. Synechococcus spp. also flourished towards the end of this phase. In parallel, we observed a total 20 % significant change in the metaproteome of the phytoplankton community. This is considered a significant alteration in protein expression, having substantial impacts on cellular functions and the physiology of the organisms. Medium levels of ∆TA showed more upregulated and less downregulated proteins than higher ∆TA treatments. Under these conditions, cell viability significantly increased in pico and nanoeukaryotes-1 in intermediate alkalinity levels, while in Synechococcus spp., nanoeukaryotes-2 and microplankton remained stable. ROS levels did not significantly change in any functional group. The pigment ratios DD+DT : FUCO, and DD+DT : Chl-a increased in medium ∆TA treatments, supporting the idea of nutrient deficiency alleviation and the absence of physiological stress. Taken all data together, this study shows that there is minimal evidence indicating a harmful impact of high alkalinity on the plankton community. The OAE treatments did not result in physiological fitness impairment, thus OAE did not cause cellular stress in the phytoplankton community studied.

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: The presented pilot for the Synthesis Product for Ocean Time Series (SPOTS) includes data from 12 fixed ship-based time-series programs. The related stations represent unique open-ocean and coastal marine environments within the Atlantic Ocean, Pacific Ocean, Mediterranean Sea, Nordic Seas, and Caribbean Sea. The focus of the pilot has been placed on biogeochemical essential ocean variables: dissolved oxygen, dissolved inorganic nutrients, inorganic carbon (pH, total alkalinity, dissolved inorganic carbon, and partial pressure of CO2), particulate matter, and dissolved organic carbon. The time series used include a variety of temporal res- olutions (monthly, seasonal, or irregular), time ranges (10–36 years), and bottom depths (80–6000 m), with the oldest samples dating back to 1983 and the most recent one corresponding to 2021. Besides having been harmo- nized into the same format (semantics, ancillary data, units), the data were subjected to a qualitative assessment in which the applied methods were evaluated and categorized. The most recently applied methods of the time- series programs usually follow the recommendations outlined by the Bermuda Time Series Workshop report (Lorenzoni and Benway, 2013), which is used as the main reference for “method recommendations by prevalent initiatives in the field”. However, measurements of dissolved oxygen and pH, in particular, still show room for improvement. Additional data quality descriptors include precision and accuracy estimates, indicators for data variability, and offsets compared to a reference and widely recognized data product for the global ocean: the GLobal Ocean Data Analysis Project (GLODAP). Generally, these descriptors indicate a high level of continuity in measurement quality within time-series programs and a good consistency with the GLODAP data product, even though robust comparisons to the latter are limited. The data are available as (i) a merged comma-separated file that is compliant with the World Ocean Circulation Experiment (WOCE) exchange format and (ii) a format dependent on user queries via the Environmental Research Division’s Data Access Program (ERDDAP) server of the Global Ocean Observing System (GOOS). The pilot increases the data utility, findability, accessibility, interoperability, and reusability following the FAIR philosophy, enhancing the readiness of biogeochemical time series. It facilitates a variety of applications that benefit from the collective value of biogeochemical time-series observations and forms the basis for a sustained time-series living data product, SPOTS, complementing relevant products for the global interior ocean carbon data (GLobal Ocean Data Analysis Project), global surface ocean carbon data (Surface Ocean CO2 Atlas; SOCAT), and global interior and surface methane and nitrous oxide data (MarinE MethanE and NiTrous Oxide product). Aside from the actual data compilation, the pilot project produced suggestions for reporting metadata, im- plementing quality control measures, and making estimations about uncertainty. These recommendations aim to encourage the community to adopt more consistent and uniform practices for analysis and reporting and to update these practices regularly. The detailed recommendations, links to the original time-series programs, the original data, their documentation, and related efforts are available on the SPOTS website. This site also pro- vides access to the data product (DOI: https://doi.org/10.26008/1912/bco-dmo.896862.2, Lange et al., 2024) and ancillary data.

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: In this paper we describe the implementation of the carbon isotopes 13C and 14C (radiocarbon) into the marine biogeochemistry model REcoM3. The implementation is tested in long-term equilibrium simulations where REcoM3 is coupled with the ocean general circulation model FESOM2.1, applying a low-resolution configuration and idealized climate forcing. Focusing on the carbon-isotopic composition of dissolved inorganic carbon (δ13CDIC and Δ14CDIC), our model results are largely consistent with reconstructions for the pre-anthropogenic period. Our simulations also exhibit discrepancies, e.g. in upwelling regions and the interior of the North Pacific. Some of these differences are due to the limitations of our ocean circulation model setup, which results in a rather shallow meridional overturning circulation. We additionally study the accuracy of two simplified modelling approaches for dissolved inorganic 14C, which are faster (15 % and about a factor of five, respectively) than the complete consideration of the marine radiocarbon cycle. The accuracy of both simplified approaches is better than 5 %, which should be sufficient for most studies of Δ14CDIC.

Description: We compare Holocene tree-cover changes in Europe derived from a transient MPI-ESM1.2 simulation with high spatial resolution LPJ-GUESS time-slice simulations and pollen-based quantitative reconstructions of tree cover based on the REVEALS model. The dynamic vegetation models and REVEALS agree with respect to the general temporal trends in tree cover for most parts of Europe, with a large tree cover during the mid-Holocene and a substantially smaller tree cover closer to the present time. However, the decrease in tree cover in REVEALS starts much earlier than in the models indicating much earlier anthropogenic deforestation than the prescribed land-use in the models. While LPJ-GUESS generally overestimates tree cover compared to the reconstructions, MPI-ESM indicates lower percentages of tree cover than REVEALS, particularly in Central Europe and the British Isles. A comparison of the simulated climate with chironomid-based climate reconstructions reveals that model-data mismatches in tree cover are in most cases not driven by biases in the climate. Instead, sensitivity experiments indicate that the model results strongly depend on the tuning of the models regarding natural disturbance regimes (e.g. fire and wind throw). The frequency and strength of disturbances are – like most of the parameters in the vegetation models – static and calibrated to modern conditions. However, these parameter values may not be valid during climate and vegetation states totally different from todays. In particular, the mid-Holocene natural forests were probably more stable and less sensitive to disturbances than present day forests that are heavily altered by human interventions. Our analysis highlights the fact that such model settings are inappropriate for palaeo-simulations and complicate model-data comparisons with additional challenges. Moreover, our study suggests that land-use is the main driver of forest decline in Europe during the mid- and late-Holocene.

Description: Nitrous oxide (N2O) is a long-lived potent greenhouse gas and stratospheric ozone-depleting substance, which has been accumulating in the atmosphere since the pre-industrial period. The mole fraction of atmospheric N2O has increased by nearly 25 % from 270 parts per billion (ppb) in 1750 to 336 ppb in 2022, with the fastest annual growth rate since 1980 of more than 1.3 ppb yr-1 in both 2020 and 2021. As a core component of our global greenhouse gas assessments coordinated by the Global Carbon Project (GCP), we present a global N2O budget that incorporates both natural and anthropogenic sources and sinks, and accounts for the interactions between nitrogen additions and the biochemical processes that control N2O emissions. We use Bottom-Up (BU: inventory, statistical extrapolation of flux measurements, process-based land and ocean modelling) and Top-Down (TD: atmospheric measurement-based inversion) approaches. We provide a comprehensive quantification of global N2O sources and sinks in 21 natural and anthropogenic categories in 18 regions between 1980 and 2020. We estimate that total annual anthropogenic N2O emissions increased 40 % (or 1.9 Tg N yr-1) in the past four decades (1980–2020). Direct agricultural emissions in 2020, 3.9 Tg N yr−1 (best estimate) represent the large majority of anthropogenic emissions, followed by other direct anthropogenic sources (including ‘Fossil fuel and industry’, ‘Waste and wastewater’, and ‘Biomass burning’ (2.1 Tg N yr−1), and indirect anthropogenic sources (1.3 Tg N yr−1). For the year 2020, our best estimate of total BU emissions for natural and anthropogenic sources was 18.3 (lower-upper bounds: 10.5–27.0) Tg N yr-1, close to our TD estimate of 17.0 (16.6–17.4) Tg N yr-1. For the period 2010–2019, the annual BU decadal-average emissions for natural plus anthropogenic sources were 18.1 (10.4–25.9) Tg N yr-1 and TD emissions were 17.4 (15.8–19.20 Tg N yr-1. The once top emitter Europe has reduced its emissions since the 1980s by 31 % while those of emerging economies have grown, making China the top emitter since the 2010s. The observed atmospheric N2O concentrations in recent years have exceeded projected levels under all scenarios in the Coupled Model Intercomparison Project Phase 6 (CMIP6), underscoring the urgency to reduce anthropogenic N2O emissions. To evaluate mitigation efforts and contribute to the Global Stocktake of the United Nations Framework Convention on Climate Change, we propose establishing a global network for monitoring and modeling N2O from the surface through the stratosphere. The data presented in this work can be downloaded from https://doi.org/10.18160/RQ8P-2Z4R (Tian et al. 2023).

Description: Black shale sediments from the Barremian to Aptian South Atlantic document intense and widespread burial of marine organic carbon during the initial stages of seafloor spreading between Africa and South America. The enhanced sequestration of atmospheric CO2 makes these young ocean basins potential drivers of the Early Cretaceous carbon cycle and climate perturbations. The opening of marine gateways between initially restricted basins and related circulation and ventilation changes are a commonly invoked explanation for the transient formation and disappearance of these regional carbon sinks. However, large uncertainties in paleogeographic reconstructions limit the interpretation of available paleoceanographic data and prevent any robust model-based quantifications of the proposed circulation and carbon burial changes. Here, we present a new approach to assess the principal controls on the Early Cretaceous South Atlantic and Southern Ocean circulation changes under full consideration of the uncertainties in available boundary conditions. Specifically, we use a large ensemble of 36 climate model experiments to simulate the Barremian to Albian progressive opening of the Falkland Plateau and Georgia Basin gateways with different configurations of the proto-Drake Passage, the Walvis Ridge, and atmospheric CO2 concentrations. The experiments are designed to complement available geochemical data across the regions and to test circulation scenarios derived from them. All simulations show increased evaporation and intermediate water formation at subtropical latitudes that drive a meridional overturning circulation whose vertical extent is determined by the sill depth of the Falkland Plateau. Densest water masses formed in the southern Angola Basin and potentially reached the deep Cape Basin as Walvis Ridge Overflow Water. Paleogeographic uncertainties are as important as the lack of precise knowledge of atmospheric CO2 levels for the simulated temperature and salinity spread in large parts of the South Atlantic. Overall temperature uncertainties are up to 15 °C and increase significantly with water depth. The ensemble approach reveals temporal changes in the relative importance of geographic and radiative forcings for the simulated oceanographic conditions and, importantly, nonlinear interactions between them. Progressive northward opening of the highly restricted Angola Basin increased the sensitivity of local overturning and upper ocean stratification to atmospheric CO2 concentrations due to large-scale changes in the hydrological cycle, while the chosen proto-Drake Passage depth is critical for the ocean dynamics and CO2 response in the southern South Atlantic. Finally, the simulated processes are integrated into a recent carbon burial framework to document the principal control of the regional gateway evolution on the progressive shift from the prevailing saline and oxygen-depleted subtropical water masses to the dominance of ventilated high-latitude deep waters.

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: Ocean alkalinity enhancement (OAE) has been proposed as a carbon dioxide removal technology (CDR) allowing for long term storage of carbon dioxide in the ocean. By changing the carbonate speciation in seawater, OAE may potentially alter marine ecosystems with implications for the biological carbon pump. Using mesocosmsthe subtropical North Atlantic, we provide first empirical insights into impacts of carbonate-based OAE on the vertical flux and attenuation of sinking particles in an oligotrophic plankton community. We enhanced total alkalinity (TA) in increments of 300 μmol kg-1, reaching up to ΔTA = 2400 µmol kg-1 compared to ambient TA. We applied a pCO2-equilibrated OAE approach, i.e. dissolved inorganic carbon (DIC) was raised simultaneously with TA to maintain seawater pCO2 in equilibrium with the atmosphere, thereby keeping perturbations of seawater carbonate chemistry moderate. The vertical flux of major elements including carbon, nitrogen, phosphorus and silicon, as well as their stoichiometric ratios (e.g. carbon-to-nitrogen) remained unaffected over 29 days of OAE. The particle properties controlling the flux attenuationinking velocities and remineralization rates also remained unaffected by OAE. However, we observed abiotic mineral precipitation at high OAE levels (ΔTA = 1800 μmol kg-1 and higher) that resulted in a substantial increase in PIC formation. The associated consumption of alkalinity reduces the efficiency of CO2 removal and emphasizes the importance of maintaining OAE within a carefully defined operating range. Our findings suggest that carbon export by oligotrophic plankton communities is insensitive to OAE perturbations using a CO2 pre-equilibrated approach. The integrity of ecosystem services is a prerequisite for large-scale application and should be further tested across a variety of nutrient-regimes and for less idealized OAE approaches.

Description: Nitric oxide (NO) is an intermediate of various microbial nitrogen cycle processes and the open ocean and coastal areas are generally a source of NO in the atmosphere. However, our knowledge about its distribution and the main production processes in coastal areas and estuaries is rudimentary at best. To this end, dissolved NO concentrations were measured for the first time in surface waters along the lower Elbe Estuary and Hamburg Port area in July 2021. The discrete surface water samples were analyzed using a chemiluminescence detection method. The NO concentrations ranged from below the limit of detection (9.1 × 10−12 mol L−1) to 17.7 × 10−12 mol L−1, averaging at 12.5 × 10−12 mol L−1 and were supersaturated in the surface layer of both the lower Elbe Estuary and the Hamburg Port area, indicating that the study site was a source of NO to the atmosphere during the study period. On the basis of a comprehensive comparison of NO concentrations with parallel nutrient, oxygen, and nitrous oxide concentration measurements, we conclude that the observed distribution of dissolved NO was most likely resulting from microbial nitrogen transformation processes, particularly nitrification in the coastal-brackish and limnic zones of the lower Elbe Estuary and nitrifier-denitrification in the Hamburg Port area.

Description: Ocean alkalinity enhancement (OAE) is an emerging strategy that aims to mitigate climate change by increasing the alkalinity of seawater. This approach involves increasing the alkalinity of the ocean to enhance its capacity to absorb and store carbon dioxide (CO2) from the atmosphere. This chapter presents an overview of the technical aspects associated with the full range of OAE methods being pursued and discusses implications for undertaking research on these approaches. Various methods have been developed to implement OAE, including the direct injection of alkaline liquid into the surface ocean; dispersal of alkaline particles from ships, platforms, or pipes; the addition of minerals to coastal environments; and the electrochemical removal of acid from seawater. Each method has its advantages and challenges, such as scalability, cost effectiveness, and potential environmental impacts. The choice of technique may depend on factors such as regional oceanographic conditions, alkalinity source availability, and engineering feasibility. This chapter considers electrochemical methods, the accelerated weathering of limestone, ocean liming, the creation of hydrated carbonates, and the addition of minerals to coastal environments. In each case, the technical aspects of the technologies are considered, and implications for best-practice research are drawn. The environmental and social impacts of OAE will likely depend on the specific technology and the local context in which it is deployed. Therefore, it is essential that the technical feasibility of OAE is undertaken in parallel with, and informed by, wider impact assessments. While OAE shows promise as a potential climate change mitigation strategy, it is essential to acknowledge its limitations and uncertainties. Further research and development are needed to understand the long-term effects, optimize techniques, and address potential unintended consequences. OAE should be viewed as complementary to extensive emission reductions, and its feasibility may be improved if it is operated using energy and supply chains with minimal CO2 emissions.

Description: The Paris Agreement to limit global warming to well below 2 °C requires ambitious emission reduction and the balancing of remaining emissions through carbon sinks, i.e. the deployment of carbon dioxide removal (CDR). While ambitious climate protection scenarios until now consider primarily land-based CDR methods, there is growing concern about their potential to deliver sufficient CDR, and marine CDR options receive more and more interest. Based on idealized theoretical studies, Ocean Alkalinity Enhancement (OAE) appears as a promising marine CDR method. However, the knowledge base is insufficient for a robust assessment of its practical feasibility, of its side effects, social and governance aspects as well as monitoring, reporting and verification issues. A number of research efforts aim to improve this in a timely manner. We provide an overview on the current situation of developing OAE as marine CDR method, and describe the history that has led to the creation of the OAE research Best Practices Guide.

Description: Ocean alkalinity enhancement (OAE) is a proposed marine carbon dioxide removal (mCDR) approach that has the potential for large-scale uptake of significant amounts of atmospheric carbon dioxide (CO2). Removing anthropogenic legacy CO2 will be required to stabilise global surface temperatures below the 1.5–2 ∘C Paris Agreement target of 2015. In this chapter we describe the impacts of various OAE feedstocks on seawater carbonate chemistry, as well as pitfalls that need to be avoided during sampling, storage, and measurement of the four main carbonate chemistry parameters, i.e. dissolved inorganic carbon (DIC), total alkalinity (TA), pH, and CO2 fugacity (fCO2). Finally, we also discuss considerations in regard to calculating carbonate chemistry speciation from two measured parameters. Key findings are that (1) theoretical CO2 uptake potential (global mean of 0.84 mol of CO2 per mole of TA added) based on carbonate chemistry calculations is probably secondary in determining the oceanic region in which OAE would be best; (2) carbonate chemistry sampling is recommended to involve gentle pressure filtration to remove calcium carbonate (CaCO3) that might have been precipitated upon TA increase as it would otherwise interfere with a number of analyses; (3) samples for DIC and TA can be stabilised to avoid the risk of secondary CaCO3 precipitation during sample storage; and (4) some OAE feedstocks require additional adjustments to carbonate chemistry speciation calculations using available programs and routines, for instance if seawater magnesium or calcium concentrations are modified.

Description: The deliberate increase in ocean alkalinity (referred to as ocean alkalinity enhancement, or OAE) has been proposed as a method for removing CO2 from the atmosphere. Before OAE can be implemented safely, efficiently, and at scale several research questions have to be addressed, including (1) which alkaline feedstocks are best suited and the doses in which they can be added safely, (2) how net carbon uptake can be measured and verified, and (3) what the potential ecosystem impacts are. These research questions cannot be addressed by direct observation alone but will require skilful and fit-for-purpose models. This article provides an overview of the most relevant modelling tools, including turbulence-, regional-, and global-scale biogeochemical models and techniques including approaches for model validation, data assimilation, and uncertainty estimation. Typical bio- geochemical model assumptions and their limitations are discussed in the context of OAE research, which leads to an identification of further development needs to make models more applicable to OAE research questions. A description of typical steps in model validation is followed by proposed minimum criteria for what constitutes a model that is fit for its intended purpose. After providing an overview of approaches for sound integration of models and observations via data assimilation, the application of observing system simulation experiments (OSSEs) for observing system design is described within the context of OAE research. Criteria for model val- idation and intercomparison studies are presented. The article concludes with a summary of recommendations and potential pitfalls to be avoided.

Description: An essential prerequisite for the implementation of ocean alkalinity enhancement (OAE) applications is their environmental safety. Only if it can be ensured that ecosystem health and ecosystem services are not at risk will the implementation of OAE move forward. Public opinion on OAEs will depend first and foremost on reliable evidence that no harm will be done to marine ecosystems and licensing authorities will demand measurable criteria against which environmental sustainability can be determined. In this context mesocosm experiments represent a highly valuable tool in determining the safe operating space of OAE applications. By combining realism and biological complexity with controllability and replication they provide an ideal OAE test bed and a critical stepping stone towards field applications. Mesocosm approaches can also be helpful in testing the efficacy, efficiency and permanence of OAE applications. This chapter outlines strengths and weaknesses of mesocosm approaches, illustrates mesocosm facilities and suitable experimental designs presently employed in OAE research, describes critical steps in mesocosm operation, and discusses possible approaches for alkalinity manipulation and monitoring. Building on a general treatise on each of these aspects, the chapter describes pelagic and benthic mesocosm approaches separately, given their inherent differences. The chapter concludes with recommendations for best practices in OAE-related mesocosm research.

Description: Carbon monoxide (CO) is an atmospheric trace gas that plays a crucial role in the oxidizing capacity of the Earth’s atmosphere. Moreover, it functions as an indirect greenhouse gas, influencing the lifetimes of potent greenhouse gases such as methane. Albeit being an overall source of atmospheric CO, the role of coastal regions in the marine cycling of CO and how its budget can be affected by anthropogenic activities, remain uncertain. Here, we present the first measurements of dissolved CO in the Ria Formosa Lagoon, an anthropogenically influenced system in southern Portugal. The dissolved CO concentrations in the surface layer ranged from 0.16 to 3.1 nmol L−1 with an average concentration of 0.75 ± 0.57 nmol L−1. The CO saturation ratio ranged from 1.7 to 32.2, indicating that the lagoon acted as a source of CO to the atmosphere in May 2021. The estimated average sea-to-air flux density was 1.53 μmol m−2 d−1, mainly fueled by CO photochemical production. Microbial consumption accounted for 83 % of the CO production, suggesting that the resulting CO emissions to the atmosphere were modulated by microbial consumption in the surface waters of the Ria Formosa Lagoon. The results from an irradiation experiment with aquaculture effluent water indicated that aquaculture facilities in the Ria Formosa Lagoon seem to be a negligible source of atmospheric CO.

Description: Ocean alkalinity enhancement (OAE) is a marine carbon dioxide removal (CDR) approach. Publicly funded research projects have begun, and philanthropic funding and start-ups are collectively pushing the field forward. This rapid progress in research activities has created an urgent need to learn if and how OAE can work at scale. The Best Practices Guide to OAE research contains 7 topics broken down into 13 chapters that compare and synthesise previously published methods and offer guidance for future research.

Description: As one of Earth's most productive marine ecosystems, the Peruvian upwelling system transports large amounts of biogenic matter from the surface to the deep ocean. Whilst particle sinking velocity is a key factor controlling the biological pump, thereby affecting carbon sequestration and O2-depletion, it has not yet been measured in this system. During a 50 d mesocosm experiment in the surface waters off the coast of Peru, we assessed particle sinking velocities and their biogeochemical and physical drivers. We further characterized the general properties of exported particles under different phytoplankton communities and nutritional states. Average sinking velocities varied between size classes and ranged from 12.8 ± 0.7 m d−1 (particles 40–100 µm) to 19.4 ± 0.7 m d−1 (particles 100–250 µm) and 34.2 ± 1.5 m d−1 (particles 250–1000 µm) (± 95 % CI). Despite a distinct plankton succession from diatoms to dinoflagellates with concomitant 5-fold drop in opal ballasting, substantial changes in sinking velocity were not observed. This illustrates the complexity of counteracting factors driving the settling behaviour of marine particles. In contrast, we found higher sinking velocities with increasing particle size and roundness and decreasing porosity. Size had by far the strongest influence among these physical particle properties, despite a high amount of unexplained variability. Our study provides a detailed analysis of the drivers of particle sinking velocity in the Peruvian upwelling system, which allows modellers to optimize local particle flux parameterization. This will help to better project oxygen concentrations and carbon sequestration in a region that is subject to substantial climate-driven changes.

Description: Global biogeochemical ocean models help to investigate the present and potential future state of the ocean, its productivity and cascading effects on higher trophic levels such as fish. They are often subjectively tuned against data sets of inorganic tracers and surface chlorophyll and only very rarely against organic components such as particulate organic carbon or zooplankton. The resulting uncertainty in biogeochemical model parameters (and parameterisations) associated with these components can explain some of the large spread of global model solutions with regard to the cycling of organic matter and its impacts on biogeochemical tracer distributions, such as oxygen minimum zones (OMZs). A second source of uncertainty arises from differences in the model spin-up length as, so far, there seems to be no agreement on the required simulation time that should elapse before a global model is assessed against observations. We investigated these two sources of uncertainty by optimising a global biogeochemical ocean model against the root-mean-squared error (RMSE) of six different combinations of data sets and different spin-up times. Besides nutrients and oxygen, the observational data sets also included phyto- and zooplankton, as well as dissolved and particulate organic phosphorus (DOP and POP, respectively). We further analysed the optimised model performance with regard to global biogeochemical fluxes, oxygen inventory and OMZ volume. Following the optimisation procedure, we evaluated the RMSE for all tracers located in the upper 100 m (except for POP, for which we considered the entire vertical domain), regardless of their consideration during optimisation. For the different optimal model solutions, we find a narrow range of the RMSE, between 14 % of the average RMSE after 10 years and 24 % after 3000 years of simulation. Global biogeochemical fluxes, global oxygen bias and OMZ volume showed a much stronger divergence among the models and over time than RMSE, indicating that even models that are similar with regard to local surface tracer concentrations can perform very differently when assessed against the global diagnostics for oxygen. Considering organic tracers in the optimisation had a strong impact on the particle flux exponent (Martin b) and may reduce much of the uncertainty in this parameter and the resulting deep particle flux. Independent of the optimisation setup, the OMZ volume showed a particularly sensitive response with strong trends over time, even after 3000 years of simulation time (despite the constant physical forcing); a high sensitivity to simulation time; and the highest sensitivity to model parameters arising from the tuning strategy setup (variation of almost 80 % of the ensemble mean). In conclusion, calibration against observations of organic tracers can help to improve global biogeochemical models even after short spin-up times; here especially, observations of deep particle flux could provide a powerful constraint. However, a large uncertainty remains with regard to global OMZ volume and its evolution over time, which can show very dynamic behaviour during the model spin-up, which renders temporal extrapolation to a final equilibrium state difficult if not impossible. Given that the real ocean shows variations on many timescales, the assumption of observations representing a steady-state ocean may require some reconsideration.

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: The Peruvian upwelling system is a highly productive ecosystem with a large oxygen minimum zone (OMZ) close to the surface. In this work, we carried out a mesocosm experiment off Callao, Peru, with the addition of water masses from the regional OMZ collected at two different sites simulating two different upwelling scenarios. Here, we focus on the pelagic remineralization of organic matter by the extracellular enzyme activity of leucine aminopeptidase (LAP) and alkaline phosphatase activity (APA). After the addition of the OMZ water, dissolved inorganic nitrogen (N) was depleted, but the standing stock of phytoplankton was relatively high, even after N depletion (mostly 〉 4 µg chlorophyll a L−1). During the initial phase of the experiment, APA was 0.6 nmol L−1 h−1 even though the PO concentration was 〉 0.5 µmol L−1. Initially, the dissolved organic phosphorus (DOP) decreased, coinciding with an increase in the PO concentration that was probably linked to the APA. The LAP activity was very high, with most of the measurements in the range of 200–800 nmol L−1 h−1. This enzyme hydrolyzes terminal amino acids from larger molecules (e.g., peptides or proteins), and these high values are probably linked to the highly productive but N-limited coastal ecosystem. Moreover, the experiment took place during a rare coastal El Niño event with higher than normal surface temperatures, which could have affected enzyme activity. Using a nonparametric multidimensional scaling analysis (NMDS) with a generalized additive model (GAM), we found that biogeochemical variables (e.g., nutrient and chlorophyll-a concentrations) and phytoplankton and bacterial communities explained up to 64 % of the variability in APA. The bacterial community best explained the variability (34 %) in LAP. The high hydrolysis rates for this enzyme suggest that pelagic N remineralization, likely driven by the bacterial community, supported the high standing stock of primary producers in the mesocosms after N depletion.

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

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

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

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

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

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

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

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

Description: As Earth's atmospheric temperatures and human populations increase, more people are becoming vulnerable to natural and human-induced disasters. This is particularly true in Central America, where the growing human population is experiencing climate extremes (droughts and floods), and the region is susceptible to geological hazards, such as earthquakes and volcanic eruptions, and environmental deterioration in many forms (soil erosion, lake eutrophication, heavy metal contamination, etc.). Instrumental and historical data from the region are insufficient to understand and document past hazards, a necessary first step for mitigating future risks. Long, continuous, well-resolved geological records can, however, provide a window into past climate and environmental changes that can be used to better predict future conditions in the region. The Lake Izabal Basin (LIB), in eastern Guatemala, contains the longest known continental records of tectonics, climate, and environmental change in the northern Neotropics. The basin is a pull-apart depression that developed along the North American and Caribbean plate boundary ∼ 12 Myr ago and contains 〉 4 km of sediment. The sedimentological archive in the LIB records the interplay among several Earth System processes. Consequently, exploration of sediments in the basin can provide key information concerning: (1) tectonic deformation and earthquake history along the plate boundary; (2) the timing and causes of volcanism from the Central American Volcanic Arc; and (3) hydroclimatic, ecologic, and geomicrobiological responses to different climate and environmental states. To evaluate the LIB as a potential site for scientific drilling, 65 scientists from 13 countries and 33 institutions met in Antigua, Guatemala, in August 2022 under the auspices of the International Continental Scientific Drilling Program (ICDP) and the US National Science Foundation (NSF). Several working groups developed scientific questions and overarching hypotheses that could be addressed by drilling the LIB and identified optimal coring sites and instrumentation needed to achieve the project goals. The group also discussed logistical challenges and outreach opportunities. The project is not only an outstanding opportunity to improve our scientific understanding of seismotectonic, volcanic, paleoclimatic, paleoecologic, and paleobiologic processes that operate in the tropics of Central America, but it is also an opportunity to improve understanding of multiple geological hazards and communicate that knowledge to help increase the resilience of at-risk Central American communities.

Description: For millennia, humans have gravitated towards coastlines for their resource potential and as geopolitical centres for global trade. A basic requirement ensuring water security for coastal communities relies on a delicate balance between the supply and demand of potable water. The interaction between freshwater and saltwater in coastal settings is, therefore, complicated by both natural and human-driven environmental changes at the land-sea interface. In particular, ongoing sea level rise, warming and deoxygenation might exacerbate such perturbations. In this context, an improved understanding of the nature and variability of groundwater fluxes across the land-sea continuum is timely, yet remains out of reach. The flow of terrestrial groundwater across the coastal transition zone as well as the extent of freshened groundwater below the present-day seafloor are receiving increased attention in marine and coastal sciences because they likely represent a significant, yet highly uncertain component of (bio)geochemical budgets, and because of the emerging interest in the potential use of offshore freshened groundwater as a resource. At the same time, “reverse” groundwater flux from offshore to onshore is of prevalent socio-economic interest as terrestrial groundwater resources are continuously pressured by overpumping and seawater intrusion in many coastal regions worldwide. An accurate assessment of the land-ocean connectivity through groundwater and its potential responses to future anthropogenic activities and climate change will require a multidisciplinary approach combining the expertise of geophysicists, hydrogeologists, (bio)geochemists and modellers. Such joint activities will lay the scientific basis for better understanding the role of groundwater in societal-relevant issues such as climate change, pollution and the environmental status of the coastal oceans within the framework of the United Nations Sustainable Development Goals. Here, we present our perspectives on future research directions to better understand land-ocean connectivity through groundwater, including the spatial distributions of the essential hydrogeological parameters, highlighting technical and scientific developments, and briefly discussing its societal relevance in rapidly changing coastal oceans.

Description: The Humboldt Current Upwelling System (HCS) is the most productive eastern boundary upwelling system (EBUS) in terms of fishery yield on the planet. EBUSs are considered hotspots of climate change with predicted expansion of mesopelagic oxygen minimum zones (OMZs) and related changes in the frequency and intensity of upwelling of nutrient-rich, low-oxygen deep water. To increase our mechanistic understanding of how upwelling impacts plankton communities and trophic links, we investigated mesozooplankton community succession and gut fluorescence, fatty acid and elemental compositions (C, N, O, P), and stable isotope (δ13C, δ15N) ratios of dominant mesozooplankton and microzooplankton representatives in a mesocosm setup off Callao (Peru) after simulated upwelling with OMZ water from two different locations and different N:P signatures (moderate and extreme treatments). An oxycline between 5 and 15 m with hypoxic conditions (〈50 µmol L−1) below ∼10 m persisted in the mesocosms throughout the experiment. No treatment effects were determined for the measured parameters, but differences in nutrient concentrations established through OMZ water additions were only minor. Copepods and polychaete larvae dominated in terms of abundance and biomass. Development and reproduction of the dominant copepod genera Paracalanus sp., Hemicyclops sp., Acartia sp., and Oncaea sp. were hindered as evident from accumulation of adult copepodids but largely missing nauplii. Failed hatching of nauplii in the hypoxic bottom layer of the mesocosms and poor nutritional condition of copepods suggested from very low gut fluorescence and fatty acid compositions most likely explain the retarded copepod development. Correlation analysis revealed no particular trophic relations between dominant copepods and phytoplankton groups. Possibly, particulate organic matter with a relatively high C:N ratio was a major diet of copepods. C:N ratios of copepods and polychaetes ranged 4.8–5.8 and 4.2–4.3, respectively. δ15N was comparatively high (∼13 ‰–17 ‰), potentially because the injected OMZ source water was enriched in δ15N as a result of anoxic conditions. Elemental ratios of dinoflagellates deviated strongly from the Redfield ratio. We conclude that opportunistic feeding of copepods may have played an important role in the pelagic food web. Overall, projected changes in the frequency and intensity of upwelling hypoxic waters may make a huge difference for copepod reproduction and may be further enhanced by varying N:P ratios of upwelled OMZ water masses.

Description: When interpreting geophysical models, we need to establish a link between the models’ physical parameters and geological units. To define these connections, it is crucial to consider and compare geophysical models with multiple, independent parameters. Particularly in complex geological scenarios, such as the rifted passive margin offshore Namibia, multi-parameter analysis and joint inversion are key techniques for comprehensive geological inferences. The models resulting from joint inversion enable the definition of specific parameter combinations, which can then be ascribed to geological units. Here we perform a user-unbiased clustering analysis of the parameters electrical resistivity and density from two models derived in a joint inversion along the Namibian passive margin. We link the resulting parameter combinations to break-up related lithology, and infer the history of margin formation. This analysis enables us to clearly differentiate two types of sediment cover. Namely, one of near-shore, thick, clastic sediments, and a second one of further offshore located, more biogenic, marine sediments. Furthermore, we clearly identify areas of interlayered massive, and weathered volcanic flows, which are usually only identified in reflection seismic studies as seaward dipping reflectors. Lastly, we find a distinct difference in the signature of the transitional crust south of- and along the supposed hot-spot track Walvis Ridge. We ascribe this contrast to an increase in magmatic activity above the volcanic centre along Walvis Ridge, and potentially a change in melt sources or depth of melting. This characterizes a rift-related southern complex, and a plume-driven Walvis Ridge regime. All of these observations demonstrate the importance of multi-parameter geophysical analysis for large-scale geological interpretations. Furthermore, our results may improve future joint inversions using direct parameter coupling, by providing a guideline for the complex passive margins parameter correlations.

Description: In this paper, we review observational and modelling results on the upwelling in the tropical Atlantic between 10∘ N and 20∘ S. We focus on the physical processes that drive the seasonal variability of surface cooling and the upward nutrient flux required to explain the seasonality of biological productivity. We separately consider the equatorial upwelling system, the coastal upwelling system of the Gulf of Guinea and the tropical Angolan upwelling system. All three tropical Atlantic upwelling systems have in common a strong seasonal cycle, with peak biological productivity during boreal summer. However, the physical processes driving the upwelling vary between the three systems. For the equatorial regime, we discuss the wind forcing of upwelling velocity and turbulent mixing, as well as the underlying dynamics responsible for thermocline movements and current structure. The coastal upwelling system in the Gulf of Guinea is located along its northern boundary and is driven by both local and remote forcing. Particular emphasis is placed on the Guinea Current, its separation from the coast and the shape of the coastline. For the tropical Angolan upwelling, we show that this system is not driven by local winds but instead results from the combined effect of coastally trapped waves, surface heat and freshwater fluxes, and turbulent mixing. Finally, we review recent changes in the upwelling systems associated with climate variability and global warming and address possible responses of upwelling systems in future scenarios.

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: Increasing Greenland Ice Sheet–melting is anticipated to impact watermass transformation in the subpolar North Atlantic and ultimately the meridional overturning circulation. Complex ocean and climate models are widely applied to predict magnitude and timing of related impacts under projected future climate. We discuss the role of the ocean mean state, subpolar gyre circulation, mesoscale eddies and atmospheric coupling in shaping the response of the subpolar North Atlantic Ocean to enhanced Greenland runoff. In a suite of eight dedicated 60 to 100-year long model experiments with and without atmospheric coupling, with eddy processes parameterized and explicitly simulated, with regular and significantly enlarged Greenland runoff, we find (1) a major impact by the interactive atmosphere in enabling a compensating temperature feedback, (2) a non-negligible influence by the ocean mean state biased towards greater stability in the coupled simulations, both of which making the Atlantic Merdional Overturning Circulation less susceptible to the freshwater perturbation applied, and (3) a more even spreading of the runoff tracer in the subpolar North Atlantic and enhanced inter-gyre exchange with the subtropics in the strongly eddying simulations. Overall, our experiments demonstrate the important role of mesoscale ocean dynamics and atmosphere feedbacks in projections of the climate system response to enhanced Greenland Ice Sheet–melting and hence underline the necessity to advance scale-aware eddy parameterizations for next-generation climate models.

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: Carbon monoxide (CO) influences the radiative budget and oxidative capacity of the atmosphere over the Arctic Ocean, which is a source of atmospheric CO. Yet, oceanic CO cycling is understudied in this area, particu- larly in light of the ongoing rapid environmental changes. We present results from incubation experiments conducted in the Fram Strait in August–September 2019 under different environmental conditions: while lower pH did not affect CO production (GPCO) or consumption (kCO) rates, enhanced GPCO and kCO were positively correlated with coloured dis- solved organic matter (CDOM) and dissolved nitrate concen- trations, respectively, suggesting microbial CO uptake under oligotrophic conditions to be a driving factor for variability in CO surface concentrations. Both production and consump- tion of CO will likely increase in the future, but it is un- known which process will dominate. Our results will help to improve models predicting future CO concentrations and emissions and their effects on the radiative budget and the oxidative capacity of the Arctic atmosphere.

Description: The carbon cycle component of the newly developed Earth System Model of intermediate complexity CLIMBER-X is presented. The model represents the cycling of carbon through atmosphere, vegetation, soils, seawater and marine sediments. Exchanges of carbon with geological reservoirs occur through sediment burial, rock weathering and volcanic degassing. The state-of-the-art HAMOCC6 model is employed to simulate ocean biogeochemistry and marine sediments processes. The land model PALADYN simulates the processes related to vegetation and soil carbon dynamics, including permafrost and peatlands. The dust cycle in the model allows for an interactive determination of the input of the micro-nutrient iron into the ocean. A rock weathering scheme is implemented into the model, with the weathering rate depending on lithology, runoff and soil temperature. CLIMBER-X includes a simple representation of the methane cycle, with explicitly modelled natural emissions from land and the assumption of a constant residence time of CH4 in the atmosphere. Carbon isotopes 13C and 14C are tracked through all model compartments and provide a useful diagnostic for model-data comparison. A comprehensive evaluation of the model performance for present–day and the historical period shows that CLIMBER-X is capable of realistically reproducing the historical evolution of atmospheric CO2 and CH4, but also the spatial distribution of carbon on land and the 3D structure of biogeochemical ocean tracers. The analysis of model performance is complemented by an assessment of carbon cycle feedbacks and model sensitivities compared to state-of-the-art CMIP6 models. Enabling interactive carbon cycle in CLIMBER-X results in a relatively minor slow-down of model computational performance by ~20 %, compared to a throughput of ~10,000 simulation years per day on a single node with 16 CPUs on a high performance computer in a climate–only model setup. CLIMBER-X is therefore well suited to investigate the feedbacks between climate and the carbon cycle on temporal scales ranging from decades to 〉100,000 years.

Description: Nitrous oxide (N2O) is a greenhouse gas, with a global warming potential 298 times that of carbon dioxide. Estuaries can be sources of N2O, but their emission estimates have significant uncertainties due to limited data availability and high spatiotemporal variability. We investigated the spatial and seasonal variability of dissolved N2O and its emissions along the Elbe Estuary (Germany), a well-mixed temperate estuary with high nutrient loading from agriculture. During nine research cruises performed between 2017 and 2022, we measured dissolved N2O concentrations, as well as dissolved nutrient and oxygen concentrations along the estuary, and calculated N2O saturations, flux densities, and emissions. We found that the estuary was a year-round source of N2O, with the highest emissions in winter when dissolved inorganic nitrogen (DIN) loads and wind speeds are high. However, in spring and summer, N2O saturations and emissions did not decrease alongside lower riverine nitrogen loads, suggesting that estuarine in situ N2O production is an important source of N2O. We identified two hotspot areas of N2O production: the Port of Hamburg, a major port region, and the mesohaline estuary near the maximum turbidity zone (MTZ). N2O production was fueled by the decomposition of riverine organic matter in the Hamburg Port and by marine organic matter in the MTZ. A comparison with previous measurements in the Elbe Estuary revealed that N2O saturation did not decrease alongside the decrease in DIN concentrations after a significant improvement of water quality in the 1990s that allowed for phytoplankton growth to re-establish in the river and estuary. The overarching control of phytoplankton growth on organic matter and, subsequently, on N2O production highlights the fact that eutrophication and elevated agricultural nutrient input can increase N2O emissions in estuaries.

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

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

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

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

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

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

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

Description: Heinrich-type ice-sheet surges are one of the prominent signals of glacial climate variability. They are characterised as abrupt, quasi-periodic episodes of ice-sheet instabilities during which large numbers of icebergs are released from the Laurentide ice sheet. The mechanisms controlling the timing and occurrence of Heinrich-type ice-sheet surges remain poorly constrained to this day. Here, we use a coupled ice sheet–solid Earth model to identify and quantify the importance of boundary forcing for the surge cycle length of Heinrich-type ice-sheet surges for two prominent ice streams of the Laurentide ice sheet – the land-terminating Mackenzie ice stream and the marine-terminating Hudson ice stream. Both ice streams show responses of similar magnitude to surface mass balance and geothermal heat flux perturbations, but Mackenzie ice stream is more sensitive to ice surface temperature perturbations, a fact likely caused by the warmer climate in this region. Ocean and sea-level forcing as well as different frequencies of the same forcing have a negligible effect on the surge cycle length. The simulations also highlight the fact that only a certain parameter space exists under which ice-sheet oscillations can be maintained. Transitioning from an oscillatory state to a persistent ice streaming state can result in an ice volume loss of up to 30 % for the respective ice stream drainage basin under otherwise constant climate conditions. We show that Mackenzie ice stream is susceptible to undergoing such a transition in response to all tested positive climate perturbations. This underlines the potential of the Mackenzie region to have contributed to prominent abrupt climate change events of the last deglaciation.

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

Description: Understanding the relationship between surface marine ecosystems and the export of carbon to depth by sinking organic particles is key to representing the effect of ecosystem dynamics and diversity, and their evolution under multiple stressors, on the carbon cycle and climate in models. Recent observational technologies have greatly increased the amount of data available, both for the abundance of diverse plankton groups and for the concentration and properties of particulate organic carbon in the ocean interior. Here we use synthetic model data to test the potential of using machine learning (ML) to reproduce concentrations of particulate organic carbon within the ocean interior based on surface ecosystem and environmental data. We test two machine learning methods that differ in their approaches to data-fitting, the random forest and XGBoost methods. The synthetic data are sampled from the PlankTOM12 global biogeochemical model using the time and coordinates of existing observations. We test 27 different combinations of possible drivers to reconstruct small (POCS) and large (POCL) particulate organic carbon concentrations. We show that ML can successfully be used to reproduce modelled particulate organic carbon over most of the ocean based on ecosystem and modelled environmental drivers. XGBoost showed better results compared to random forest thanks to its gradient boosting trees' architecture. The inclusion of plankton functional types (PFTs) in driver sets improved the accuracy of the model reconstruction by 58 % on average for POCS and by 22 % for POCL. Results were less robust over the equatorial Pacific and some parts of the high latitudes. For POCS reconstruction, the most important drivers were the depth level, temperature, microzooplankton and PO4, while for POCL it was the depth level, temperature, mixed-layer depth, microzooplankton, phaeocystis, PO4 and chlorophyll a averaged over the mixed-layer depth. These results suggest that it will be possible to identify linkages between surface environmental and ecosystem structure and particulate organic carbon distribution within the ocean interior using real observations and to use this knowledge to improve both our understanding of ecosystem dynamics and of their functional representation within models.