ALBERT

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

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

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

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

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

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