ALBERT

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

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

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

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

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

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

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

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

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

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