ALBERT

Description: Beryllium isotopes have emerged as a quantitative tracer of continental weathering, but accurate and precise determination of the cosmogenic 10 Be and stable 9 Be in seawater is challenging, because seawater contains high concentrations of matrix elements but extremely low concentrations of 9 Be and 10 Be. In this study, we develop a new, time‐efficient procedure for the simultaneous preconcentration of 9 Be and 10 Be from (coastal) seawater based on the iron co‐precipitation method. The concentrations of 9 Be, 10 Be, and the resulting 10 Be/ 9 Be ratio for Changjiang Estuary water derived from the new procedure agree well with those obtained from the conventional procedure requiring separate preconcentration for 9 Be and 10 Be determinations. By avoiding the separate preconcentration, our newly developed procedure contributes toward more time‐efficient handling of samples, less sample cross‐contamination, and a more reliable 10 Be/ 9 Be ratio. Prior to this, we validated the iron co‐precipitation method using artificial seawater and natural water samples from the Amazon Estuary regarding: (1) the “matrix effect” for Be analysis, (2) its extraction efficiency for pg g −1 levels Be in the presence and absence of organic matter, and (3) the data comparability with another preconcentration method. We calculated that for the determination of 9 Be and 10 Be in most open ocean seawater with typical 10 Be concentrations of 〉 500 atoms g −1 , good precisions (〈 5%) can be achieved using less than 3 liters of seawater compared to more than 20 liters routinely used previously. Even for coastal seawater with extremely low 10 Be concentration (e.g., 100 atoms g −1 ), we estimate a maximum amount of 10 liters to be adequate.

Description: Aim: Human activities have introduced numerous non-native species (NNS) worldwide. Understanding and predicting large-scale NNS establishment patterns remain fundamental scientific challenges. Here, we evaluate if NNS composition represents a proportional subset of the total species pool available to invade (i.e. total global biodiversity), or, conversely, certain taxa are disproportionately pre-disposed to establish in non-native areas. Location: Global. Time period: Present day. Major taxa studied: Global diversity. Methods: We compiled one of the most comprehensive global databases of NNS (36,822 established species) to determine if NNS diversity is a representative proportional subset of global biodiversity. Results: Our study revealed that, while NNS diversity mirrors global biodiversity to a certain extent, due to significant deviance from the null model it is not always a representative proportional subset of global biodiversity. The strength of global biodiversity as a predictor depended on the taxonomic scale, with successive lower taxonomic levels less predictive than the one above it. Consequently, on average, 58%, 42% and 28% of variability in NNS numbers were explained by global biodiversity for phylum, class and family respectively. Moreover, global biodiversity was a similarly strong explanatory variable for NNS diversity among regions, but not habitats (i.e. terrestrial, freshwater and marine), where it better predicted NNS diversity for terrestrial than for freshwater and marine habitats. Freshwater and marine habitats were also greatly understudied relative to invasions in the terrestrial habitats. Over-represented NNS relative to global biodiversity tended to be those intentionally introduced and/or ‘hitchhikers’ associated with deliberate introductions. Finally, randomness is likely an important factor in the establishment success of NNS. Main conclusions: Besides global biodiversity, other important explanatory variables for large-scale patterns of NNS diversity likely include propagule and colonization pressures, environmental similarity between native and non-native regions, biased selection of intentionally introduced species and disparate research efforts of habitats and taxa.

Description: Severe storm flooding poses a major hazard to the coasts of north‐western Europe. However, the long‐term recurrence patterns of extreme coastal flooding and their governing factors are poorly understood. Therefore, high‐resolution sedimentary records of past North Atlantic storm flooding are required. This multi‐proxy study reconstructs storm‐induced overwash processes from coastal lake sediments on the Shetland Islands using grain‐size and geochemical data, and the re‐analysis of historical data. The chronostratigraphy is based on Bayesian age–depth modelling using accelerator mass spectrometry 14 C and 137 Cs data. A high XRF‐based Si/Ti ratio and the unimodal grain‐size distribution link the sand layers to the beach and thus storm‐induced overwash events. Periods with more frequent storm flooding occurred 980–1050, 1150–1300, 1450–1550, 1820–1900 and 1950–2000 ce, which is largely consistent with a positive North Atlantic Oscillation mode. The Little Ice Age (1400–1850 ce ) shows a gap of major sand layers suggesting a southward shift of storm tracks and a seasonal variance with more storm floods in spring and autumn. Warmer phases shifted winter storm tracks towards the north‐east Atlantic, indicating a possible trend for future storm‐track changes and increased storm flooding in the northern North Sea region.

Description: The seafloor morphology reflects both past and on‐going sedimentary, oceanographic and tectonic processes. Vertical movement is one of the drivers responsible for reshaping the seafloor through forming steep flanks that decrease slope stability, favour landslides, change current paths, form minibasins and control the sediment deposition, distribution and geometry. Here, we make use of these interactions to derive vertical movements and constrain the active tectonic processes at the western termination of the upper Calabrian accretionary wedge from the integrated analysis of bathymetric, backscatter, surface attributes and high‐resolution reflection seismic data. Within this area, we identify two types of deformational features and mechanisms that affect the depositional, erosional and tectonic processes at different scales. These include the deviation of channels, landslide scars, mass transport deposits (MTDs), separated drifts, sediment waves, lineaments and offset seafloor structures. The first type (long‐wavelength uplift) is an uplifted 22‐km‐wide region, in which seismic onlap relationships and the dip of deep reflectors suggest long‐lasting but slow tectonic uplift affecting sedimentation, and the second type (short‐wavelength uplift) includes three narrow elongated structures and one circular dome encircling the first region of uplift. We interpret that the first type of uplift feature was caused by tectonic deformation, while the second type is interpreted as formed by the fast uplift, tilting and faulting of modern sediments caused by diapirism due to rapid sedimentation in response to the first tectonically driven uplift. The study provides insight into the complex interaction of tectonic and sedimentary processes in the upper Calabrian accretionary wedge.

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: The marine biological carbon pump (BCP) stores carbon in the ocean interior, isolating it from exchange with the atmosphere and thereby coregulating atmospheric carbon dioxide (CO 2 ). As the BCP commonly is equated with the flux of organic material to the ocean interior, termed “export flux,” a change in export flux is perceived to directly impact atmospheric CO 2 , and thus climate. Here, we recap how this perception contrasts with current understanding of the BCP, emphasizing the lack of a direct relationship between global export flux and atmospheric CO 2 . We argue for the use of the storage of carbon of biological origin in the ocean interior as a diagnostic that directly relates to atmospheric CO 2 , as a way forward to quantify the changes in the BCP in a changing climate. The diagnostic is conveniently applicable to both climate model data and increasingly available observational data. It can explain a seemingly paradoxical response under anthropogenic climate change: Despite a decrease in export flux, the BCP intensifies due to a longer reemergence time of biogenically stored carbon back to the ocean surface and thereby provides a negative feedback to increasing atmospheric CO 2 . This feedback is notably small compared with anthropogenic CO 2 emissions and other carbon‐climate feedbacks. In this Opinion paper, we advocate for a comprehensive view of the BCP's impact on atmospheric CO 2 , providing a prerequisite for assessing the effectiveness of marine CO 2 removal approaches that target marine biology.

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: This study was the first to investigate the key reproductive traits of the electric lantern fish Electrona risso (Myctophidae, n = 918) and the bigscale fishes (Melamphaidae) Melamphaes polylepis (n = 260) and Scopelogadus mizolepis (n = 649). Specimens of these mesopelagic species were collected in March and April 2015 in the eastern Central Atlantic (0–24° N, 20–26° W). Sex ratio was not significantly different from 1:1 in E. risso and M. polylepis but significantly skewed toward female dominance in S. mizolepis. Reproductive phases were determined macroscopically and by histological analyses on selected individuals. Female length at 50% maturity (L50) was 55.1 mm standard length (LS) in E. risso, with an observed female maximum length (Lmax) of 81.2 mm LS. M. polylepis females had an L50 of 40.2 mm LS and an Lmax of 86.7 mm LS. S. mizolepis had an L50 of 46 mm LS and an Lmax of 97.9 mm LS. The three species show histological features of iteroparity, but the E. risso population appears to occur in two year-classes and experience only one spawning season per lifetime in the study region. All three species are batch-spawners. A batch fecundity of 2668 eggs was estimated from one E. risso individual, with a relative batch fecundity of 369 eggs g−1 gonad-free body mass. M. polylepis had a batch fecundity of 1027 eggs and a relative batch fecundity of 149 eggs g−1 (n = 3). S. polylepis had a batch fecundity of 1545 eggs and a relative batch fecundity of 215 eggs g−1 (n = 21). The median gonado-somatic index during the actively spawning phase of E. risso was 4.5, significantly lower than that of M. polylepis (7.5) and S. mizolepis (7.1). No regressing or regenerating phases were observed in this study. Batch-spawning in all three species is suggested to be advantageous to cope with intra-annual variability in food supply and other risks for offspring survival. With what appears to be in effect a (facultative) semelparous strategy in combination with a short life span in E. risso, interannual differences would have a great effect on population dynamics of this species. Knowledge is still lacking on temporal aspects of reproduction such as the duration of the spawning season and the frequency of spawning, as well as age and growth.

Description: Gelatinous zooplankton are increasingly recognized to play a key role in the ocean's biological carbon pump. Appendicularians, a class of pelagic tunicates, are among the most abundant gelatinous plankton in the ocean, but it is an open question how their contribution to carbon export might change in the future. Here, we conducted an experiment with large volume in situ mesocosms (~55–60 m3 and 21 m depth) to investigate how ocean acidification (OA) extreme events affect food web structure and carbon export in a natural plankton community, particularly focusing on the keystone species Oikopleura dioica, a globally abundant appendicularian. We found a profound influence of O. dioica on vertical carbon fluxes, particularly during a short but intense bloom period in the high CO2 treatment, during which carbon export was 42%–64% higher than under ambient conditions. This elevated flux was mostly driven by an almost twofold increase in O. dioica biomass under high CO2. This rapid population increase was linked to enhanced fecundity (+20%) that likely resulted from physiological benefits of low pH conditions. The resulting competitive advantage of O. dioica resulted in enhanced grazing on phytoplankton and transfer of this consumed biomass into sinking particles. Using a simple carbon flux model for O. dioica, we estimate that high CO2 doubled the carbon flux of discarded mucous houses and fecal pellets, accounting for up to 39% of total carbon export from the ecosystem during the bloom. Considering the wide geographic distribution of O. dioica, our findings suggest that appendicularians may become an increasingly important vector of carbon export with ongoing OA.

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: Three boreholes drilled during the International Ocean Discovery Program (IODP) Expedition 396 have yielded unexpected findings of altered granitic rocks covered by basalt flows, interbedded sediments and glacial mud near the continent‐ocean transition of the mid‐Norwegian margin. U‐Pb and K‐Ar geochronological analyses were conducted on both protolithic and authigenically formed K‐bearing minerals to determine the age of granite crystallisation and subsequent alteration episodes. The granite's crystallisation age based on 104 zircons is 56.3 ± 0.2 Ma, and subsequent exhumation along with alteration/weathering events took place between 54.7 ± 1 and 37.1 ± 1 Ma. This intrusion represents the youngest granite discovered in Norway and intruded at an extremely shallow crustal level before a rapid rift‐to‐drift transition. The shallow emplacement of granitic rock and its fast exhumation before and during the onset of volcanism holds significant implications for the syn‐ and post‐breakup tectonic evolution of volcanic margins.

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: Key Points: - Freshwater input has significantly contributed to the surface warming at the peak of the 1995 Benguela Niño - Anomalously high river discharge and precipitation increased stratification and reduced turbulent heat loss by creating barrier layers - Combination of high freshwater input and strong poleward surface current might play a role in temperature variability off Angola Benguela Nino events are characterized by strong warm sea surface temperature (SST) anomalies off the Angolan and Namibian coasts. In 1995, the strongest event in the satellite era took place, impacting fish availability in both Angolan and Namibian waters. In this study, we use direct observations, satellite data, and reanalysis products to investigate the impact that the up-until-now unnoticed mechanism of freshwater input from Congo River discharge (CRD) and precipitation had on the evolution of the 1995 Benguela Nino. In the onset phase of the event, anomalous rainfall in November/December 1994 at around 6 degrees S, combined with a high CRD, generated a low salinity plume. The plume was advected into the Angola-Namibia region in the following February/March 1995 by an anomalously strong poleward surface current generated by the relaxation of the southerly winds and shifts in the coastal wind stress curl. The presence of this low surface salinity anomaly of about -2 psu increased ocean stability by generating barrier layers, thereby reducing the turbulent heat loss, since turbulent mixing acted on a weak vertical temperature gradient. A mixed layer heat budget analysis demonstrates that southward advection of Angolan waters drove the warming at the onset, while reduced mixing played the main role at the event's peak. We conclude that a freshwater input contributed to the SST increase in this exceptionally strong event and suggest that this input can influence the SST variability in Angola-Namibia waters through a combination of high CRD, precipitation, and the presence of a strong poleward surface current. Benguela Nino events are characterized by excessive warming of the sea surface temperature off the Angolan and Namibian coasts. One of the strongest-ever recorded warm events dates back to 1995, impacting fish availability in both Angolan and Namibian waters. In our research, we investigate if freshwater from rain and from the Congo River could have impacted the evolution of this 1995 Benguela Nino. In the event's early stage, high precipitation and river discharge generated a low salinity pool at the Congo River mouth, which in February/March 1995 was taken to the south by an exceptionally strong surface current, generated by changes in wind strength and direction at the African coast. This low sea surface salinity in a shallow layer in the upper meters of the ocean increased the ocean's stability. As the stabilized waters diminished the usual mixing from the depths below which cools down the surface waters, it contributed to an increase in warming in the surface layer of the ocean. We conclude that the warming of the surface waters in the region was indeed influenced by the combination of high precipitation and high Congo River discharge with a strong surface current toward the south. Freshwater input has significantly contributed to the surface warming at the peak of the 1995 Benguela Nino Anomalously high river discharge and precipitation increased stratification and reduced turbulent heat loss by creating barrier layers Combination of high freshwater input and strong poleward surface current might play a role in temperature variability off Angola

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: Mixed turbidite–contourite depositional systems result from interactions between down‐slope turbidity currents and along‐slope bottom currents, comprising excellent records of past oceanographic currents. Modern and ancient systems have been widely documented along the continental margins of the Atlantic Ocean. Yet, few examples have so far been identified on the North‐west African continental margin, limiting understanding of the sedimentary and palaeoceanographic evolution in this area. This work uses two‐dimensional seismic reflection profiles to report, for the first time, the presence of three giant sediment mounds beneath the headwall region of the Sahara Slide Complex. The sediment mounds are elongated and separated by two broad canyons, showing a north‐west/south‐east orientation that is roughly perpendicular to the continental margin. These mounds are 24 to 37 km long and 12 to 17 km wide, reaching a maximum height of ca 1000 m. Numerous slide scarps are observed within and along the flanks of the mounds, hinting at the occurrence of submarine landslides during their development. Based on their geometries, external shapes, internal seismic architecture and stratigraphic stacking patterns, it is proposed that these sediment mounds comprise down‐slope elongated mounded drifts formed in a mixed turbidite–contourite system during four evolutionary stages: onset, growth, maintenance and burial. The significance of this work is that it demonstrates the gradual transition from a turbidite system to a full mixed turbidite–contourite system to be associated, in the study area, with the establishment of strong ocean currents along north‐west Africa.

Description: “Flip‐flop” detachment mode represents an endmember type of lithosphere‐scale faulting observed at almost amagmatic sections of ultraslow‐spreading mid‐ocean ridges. Recent numerical experiments using an imposed steady temperature structure show that an axial temperature maximum is essential to trigger flip‐flop faults by focusing flexural strain in the footwall of the active fault. However, ridge segments without significant melt budget are more likely to be in a transient thermal state controlled, at least partly, by the faulting dynamics themselves. Therefore, we investigate which processes control the thermal structure of the lithosphere and how feedbacks with the deformation mechanisms can explain observed faulting patterns. We present results of 2‐D thermo‐mechanical numerical modeling including serpentinization reactions and dynamic grain size evolution. The model features a novel form of parametrized hydrothermal cooling along fault zones as well as the thermal and rheological effects of periodic sill intrusions. We find that the interplay of hydrothermal fault zone cooling and periodic sill intrusions in the footwall facilitates the flip‐flop detachment mode. Hydrothermal cooling of the fault zone pushes the temperature maximum into the footwall, while intrusions near the temperature maximum further weaken the rock and promote the formation of new faults with opposite polarity. Our model allows us to put constraints on the magnitude of two processes, and we obtain most reasonable melt budgets and hydrothermal heat fluxes if both are considered. Furthermore, we frequently observe two other faulting modes in our experiments complementing flip‐flop faulting to yield a potentially more robust alternative interpretation for existing observations. Plain Language Summary At mid‐ocean ridges, two plates diverge and new seafloor is created. The nature and appearance of this new seafloor strongly depend on spreading velocity and the availability of magmatic melts. At one of the melt‐poorest and slowest‐spreading ridges, a special form of large‐scale tectonic faults, so‐called flip‐flop detachments, can be observed. Tectonic faults can act as pathways for fluids circulating through the seafloor, which provides a significant cooling effect for the young plate. The interplay of magmatic activity, faulting and fluid circulation is evident at many different ridges with different magmatic activity and spreading rates. Flip‐flop faulting is restricted to only a few ridge sections worldwide, and we here investigate the prerequisites for this special spreading mode. To do so, we set up a computer model of an ultraslow‐spreading mid‐ocean ridge including the effects of sparse magmatism as well as the cooling effect associated with fluid circulation. We find that feedbacks between faulting dynamics, hydrothermal cooling and magmatic activity control the magnitude and spatial location of each individual process. Seafloor and subsurface observations are best explained by calculations with moderate melt input and hydrothermal circulation acting together. Key Points We implemented hydrothermal cooling and magmatic intrusion in a thermo‐mechanical model to explain detachment faulting at ultraslow ridges Stable flip‐flop detachment faulting is observed for setups considering both melt input and hydrothermal heat fluxes at realistic magnitudes Two other faulting modes frequently observed in our model offer potential alternative interpretations for existing seafloor observations

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

Description: Nutrient transfer into the sunlit surface ocean by cyclonic eddies is potentially crucial for sustaining primary productivity in the stratified subtropical gyres. However, the nature of productivity enhancements, including the flow of matter to higher trophic levels and its impact on carbon fluxes, remain poorly resolved. Here, we report a detailed assessment of the biogeochemical response to a cyclonic eddy in the subtropical Northwest Pacific via a combination of ship‐based and autonomous platforms. Primary production was enhanced twofold within the eddy core relative to reference sites outside, whereas phytoplankton biomass even decreased. Pico‐phytoplankton (〈 2 μ m) dominated (〉 80%) both phytoplankton biomass and primary production inside and outside the eddy. The stimulated primary production in the eddy core was accompanied by an approximately twofold increase in mesozooplankton abundance, an approximately threefold increase in particle formation in the deep chlorophyll maximum layer, as well as significantly enhanced surface oceanic CO 2 uptake and net community production. We suggest these observations carry important implications for understanding carbon export in the subtropical ocean and highlight the need to include such subtropical eddy features in ocean carbon budget analyses.

Description: Concentrations of bioavailable inorganic nitrogen (N) and phosphorus (P) are simultaneously depleted in the (sub)tropical North Atlantic Ocean, but it remains unclear if phytoplankton growth rates are N limited or N–P co‐limited. Here we present findings from three bottle‐scale experiments using a four‐by‐four matrix of low‐level N and P additions, conducted at one site in the subtropical North Atlantic Ocean. Phytoplankton responses were assessed both in terms of bulk chlorophyll a (Chl a ) concentrations and intracellular Chl a of dominant Prochlorococcus and Synechococcus groups. Two matrix experiments suggested that N was independently limiting in situ growth, with no co‐limiting role for P, while the third showed co‐limitation by both N and P in this region. This switch from N limitation to N–P co‐limitation was attributed to an episodic wet deposition event that supplied N, thereby stimulating phytoplankton growth and consuming available P. Such rapid transitions in nutrient limitation in response to environmental forcing might be common in oceanic systems with multiple depleted nutrients.

Description: Rationale Potassium (K) is a major component of several silicate minerals and seawater, and, therefore, constraining past changes in the potassium cycle is a promising way of tracing large‐scale geological processes on Earth. However, [K] measurement using inductively coupled plasma mass spectrometry (ICP‐MS) is challenging due to an ArH + interference, which may be of a similar magnitude to the K + ion beam in samples with 〈0.1% m/m [K]. Methods In this work, we investigated the effect of the ArH + interference on K/Ca data quality by comparing results from laser‐ablation (LA)‐ICP‐MS measured in medium and high mass resolution modes and validating our LA results via solution ICP‐optical emission spectroscopy (OES) and solution ICP‐MS measurements. To do so, we used a wide range of geological reference materials, with a particular focus on marine carbonates, which are potential archives of past changes in the K cycle but are typically characterised by [K] 〈 200 μg/g. In addition, we examine the degree to which trace‐element data quality is driven by downhole fractionation during LA‐ICP‐MS measurements. Results Our results show that medium mass resolution (MR) mode is sufficiently capable of minimising the effect of the ArH + interference on K + . However, the rate of downhole fractionation for Na and K varies between different samples as a result of their differing bulk composition, resulting in matrix‐specific inaccuracy. We show how this can be accounted for via downhole fractionation corrections, resulting in an accuracy of better than 1% and a long‐term reproducibility (intermediate precision) of 〈6% (relative standard deviation) in JCp‐1NP using LA‐ICP‐MS in MR mode. Conclusion Our [K] measurement protocol is demonstrably precise and accurate and applicable to a wide range of materials. The measurement of K/Ca in relatively low‐[K] marine carbonates is presented here as a key example of a new application opened up by these advances.

Description: The Black Sea is a permanently anoxic, marine basin serving as model system for the deposition of organic-rich sediments in a highly stratified ocean. In such systems, archaeal lipids are widely used as paleoceanographic and biogeochemical proxies; however, the diverse planktonic and benthic sources as well as their potentially distinct diagenetic fate may complicate their application. To track the flux of archaeal lipids and to constrain their sources and turnover, we quantitatively examined the distributions and stable carbon isotopic compositions (delta 13C) of intact polar lipids (IPLs) and core lipids (CLs) from the upper oxic water column into the underlying sediments, reaching deposits from the last glacial. The distribution of IPLs responded more sensitively to the geochemical zonation than the CLs, with the latter being governed by the deposition from the chemocline. The isotopic composition of archaeal lipids indicates CLs and IPLs in the deep anoxic water column have negligible influence on the sedimentary pool. Archaeol substitutes tetraether lipids as the most abundant IPL in the deep anoxic water column and the lacustrine methanic zone. Its elevated IPL/CL ratios and negative delta 13C values indicate active methane metabolism. Sedimentary CL- and IPL-crenarchaeol were exclusively derived from the water column, as indicated by non-variable delta 13C values that are identical to those in the chemocline and by the low BIT (branched isoprenoid tetraether index). By contrast, in situ production accounts on average for 22% of the sedimentary IPL-GDGT-0 (glycerol dibiphytanyl glycerol tetraether) based on isotopic mass balance using the fermentation product lactate as an endmember for the dissolved substrate pool. Despite the structural similarity, glycosidic crenarchaeol appears to be more recalcitrant in comparison to its non-cycloalkylated counterpart GDGT-0, as indicated by its consistently higher IPL/CL ratio in sediments. The higher TEX86, CCaT, and GDGT-2/-3 values in glacial sediments could plausibly result from selective turnover of archaeal lipids and/or an archaeal ecology shift during the transition from the glacial lacustrine to the Holocene marine setting. Our in-depth molecular-isotopic examination of archaeal core and intact polar lipids provided new constraints on the sources and fate of archaeal lipids and their applicability in paleoceanographic and biogeochemical studies.

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: Retrieving the physical properties and water content of marine aerosols requires understanding the links between the particles' optical and microphysical properties. By using a morphologically realistic model with varying salt mass fractions fm, describing the transition from irregularly shaped, dry salt crystals to brine-coated geometries, optical properties relevant to polarimetric remote sensing are computed at wavelengths of 532 and 1,064 nm. The extinction cross section and its color ratio depend on particle size, but are insensitive to changes in fm; thus, measured extinction coefficients at two wavelengths contain information on both particle number and size. The lidar ratio's dependence on both size and wavelength has implications for inverting the lidar equation. The results suggest that active observations of the backscattering cross section's color ratio and the depolarization ratio, as well as, passive observations of the degree of linear polarization offer avenues to obtain the water content of marine aerosols. Key Points: - Modeled extinction coefficient of marine aerosol depends on particle radius and wavelength, but not on water content - The depolarization ratio and the color ratio of the backscattering cross section generally decrease with growing aerosol water content - The linear polarization peak near backscattering angles at NIR wavelength could be used in passive polarimetry to retrieve water content

Description: Current earthquake forecasting approaches are mainly based on probabilistic assumptions, as earthquakes seem to occur randomly. Such apparent randomness can however be caused by deterministic chaos, rendering deterministic short‐term forecasts possible. Due to the short historical and instrumental record of earthquakes, chaos detection has proven challenging, but more frequently occurring slow slip events (SSE) are promising candidates to probe for determinism. Here, we characterize the SSE signatures obtained from GNSS position time series in the Hikurangi Subduction Zone (New Zealand) to investigate whether the seemingly random SSE occurrence is governed by chaotic determinism. We find evidence for deterministic chaos for stations recording shallow SSEs, suggesting that short‐term deterministic forecasting of SSEs, similar to weather forecasts, might indeed be possible over timescales of a few weeks. We anticipate that our findings could open the door for next‐generation SSE forecasting, adding new tools to existing probabilistic approaches. Plain Language Summary Since earthquakes appear to occur randomly, the currently available probabilistic predictions are based on past earthquake records. These predictions estimate the likelihood of an earthquake of a given magnitude occurring within a defined time period. In contrast to such probabilistic approaches, deterministic systems are fully predictable, albeit often confined to short time scales due to their potential chaotic behavior. Probing for deterministic predictability in the earthquake cycle is intractable due to the limited historical instrumental record. However, frequently occurring slow slip events ‐ captured by transient GNSS displacements that can last several weeks ‐ provide a unique opportunity to explore deterministic predictability in these types of slow earthquakes. By studying GNSS time series from various stations on New Zealand’s North Island, we have discovered evidence suggesting that these irregularly occurring slow slip events might be governed by chaotic determinism. This implies the potential to forecast both timing and magnitude of slow slip events a few weeks in advance using deterministic methods, much like we predict weather patterns. Consequently, our theoretical findings could therefore pave the way for innovative approaches to short‐term slow slip forecasting. Key Points Nonlinear analysis of GNSS displacement time series unveils evidence for deterministic chaos in slow slip events in New Zealand Our theoretical findings imply that irregularly occurring slow slip events could potentially be forecasted a few weeks in advance

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

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

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

Description: The Bay of Bengal (BoB) spans 〉2.2 million km(2) in the northeastern Indian Ocean and is bordered by dense populations that depend upon its resources. Over recent decades, a shift from larger phytoplankton to picoplankton has been reported, yet the abundance, activity, and composition of primary producer communities are not well-characterized. We analysed the BoB regions during the summer monsoon. Prochlorococcus ranged up to 3.14 x 10(5) cells mL(-1 )in the surface mixed layer, averaging 1.74 +/- 0.46 x 10(5) in the upper 10 m and consistently higher than Synechococcus and eukaryotic phytoplankton. V1-V2 rRNA gene amplicon analyses showed the High Light II (HLII) ecotype formed 98 +/- 1% of Prochlorococcus amplicons in surface waters, comprising six oligotypes, with the dominant oligotype accounting for 65 +/- 4% of HLII. Diel sampling of a coherent water mass demonstrated evening onset of cell division and rapid Prochlorococcus growth between 1.5 and 3.1 div day-1, based on cell cycle analysis, as confirmed by abundance-based estimates of 2.1 div day(-1). Accumulation of Prochlorococcus produced by ultradian growth was restricted by high loss rates. Alongside prior Arabian Sea and tropical Atlantic rates, our results indicate Prochlorococcus growth rates should be reevaluated with greater attention to latitudinal zones and influences on contributions to global primary production.

Description: Seismic data from the North Sea commonly show vertical acoustic blanking (VAB) often interpreted as fluid conduits with implications for Quaternary development. The robustness of this interpretation has long been controversial as the infill of tunnel valleys can also cause vertical blanking. Using 2D and 3D seismic data and sediment echosounder data from the German North Sea, we investigate VAB to determine a geological or imaging origin of these anomalies. We detected multiple VAB occurrences throughout the North Sea. 3D data from the Ducks Beak (‘Entenschnabel’) reveal a correlation of VAB with bright spots in incised channels directly below the seafloor. Large source–receiver distances allow imaging the subsurface below the channel without signal penetrating through it (undershooting). This method removes the blanking. Energy absorption by shallow biogenic gas trapped within the channels explains the observed VAB. Hence, the blanking represents an imaging artifact, highlighting the need for careful seismic processing with sufficient offset before interpreting such anomalies as fluid pathways. The channels belong to a postglacial channel system related to the now submerged lowlands of Doggerland. This work demonstrates the usability of mapping VAB to detect shallow features for paleo‐landscape reconstruction and identification of shallow gas for hazard assessments, for example.

Description: 1. The expansion of scientific image data holds great promise to quantify individuals, size distributions and traits. Computer vision tools are especially powerful to automate data mining of images and thus have been applied widely across studies in aquatic and terrestrial ecology. Yet marine benthic communities, especially infauna, remain understudied despite their dominance of marine biomass, biodiversity and playing critical roles in ecosystem functioning. 2. Here, we disaggregated infauna from sediment cores taken throughout the spring transition (April-June) from a near-natural mesocosm setup under experimental warming (Ambient, +1.5 degrees C, +3.0 degrees C). Numerically abundant mudsnails were imaged in batches under stereomicroscopy, from which we automatically counted and sized individuals using a superpixel-based segmentation algorithm. Our segmentation approach was based on clustering superpixels, which naturally partition images by low-level properties (e.g., colour, shape and edges) and allow instance-based segmentation to extract all individuals from each image. 3. We demonstrate high accuracy and precision for counting and sizing individuals, through a procedure that is robust to the number of individuals per image (5-65) and to size ranges spanning an order of magnitude (〈750 mu m to 7.4 mm). The segmentation routine provided at least a fivefold increase in efficiency compared with manual measurements. Scaling this approach to a larger dataset tallied 〉40k individuals and revealed overall growth in response to springtime warming. 4. We illustrate that image processing and segmentation workflows can be built upon existing open-access R packages, underlining the potential for wider adoption of computer vision tools among ecologists. The image-based approach also generated reproducible data products that, alongside our scripts, we have made freely available. This work reinforces the need for next-generation monitoring of benthic communities, especially infauna, which can display differential responses to average warming.

Description: The long-term dynamics of microbial communities across geographic, hydrographic, and biogeochemical gradients in the Arctic Ocean are largely unknown. To address this, we annually sampled polar, mixed, and Atlantic water masses of the Fram Strait (2015–2019; 5–100 m depth) to assess microbiome composition, substrate concentrations, and oceanographic parameters. Longitude and water depth were the major determinants (~30%) of microbial community variability. Bacterial alpha diversity was highest in lower-photic polar waters. Community composition shifted from west to east, with the prevalence of, for example, Dadabacteriales and Thiotrichales in Arctic- and Atlantic-influenced waters, respectively. Concentrations of dissolved organic carbon peaked in the western, compared to carbohydrates in the chlorophyll-maximum of eastern Fram Strait. Interannual differences due to the time of sampling, which varied between early (June 2016/2018) and late (September 2019) phytoplankton bloom stages, illustrated that phytoplankton composition and resulting availability of labile substrates influence bacterial dynamics. We identified 10 species clusters with stable environmental correlations, representing signature populations of distinct ecosystem states. In context with published metagenomic evidence, our microbial-biogeochemical inventory of a key Arctic region establishes a benchmark to assess ecosystem dynamics and the imprint of climate change.

Description: Standardised terminology in science is important for clarity of interpretation and communication. In invasion science - a dynamic and rapidly evolving discipline - the proliferation of technical terminology has lacked a standardised framework for its development. The result is a convoluted and inconsistent usage of terminology, with various discrepancies in descriptions of damage and interventions. A standardised framework is therefore needed for a clear, universally applicable, and consistent terminology to promote more effective communication across researchers, stakeholders, and policymakers. Inconsistencies in terminology stem from the exponential increase in scientific publications on the patterns and processes of biological invasions authored by experts from various disciplines and countries since the 1990s, as well as publications by legislators and policymakers focusing on practical applications, regulations, and management of resources. Aligning and standardising terminology across stakeholders remains a challenge in invasion science. Here, we review and evaluate the multiple terms used in invasion science (e.g. 'non-native', 'alien', 'invasive' or 'invader', 'exotic', 'non-indigenous', 'naturalised', 'pest') to propose a more simplified and standardised terminology. The streamlined framework we propose and translate into 28 other languages is based on the terms (i) 'non-native', denoting species transported beyond their natural biogeographic range, (ii) 'established non-native', i.e. those non-native species that have established self-sustaining populations in their new location(s) in the wild, and (iii) 'invasive non-native' - populations of established non-native species that have recently spread or are spreading rapidly in their invaded range actively or passively with or without human mediation. We also highlight the importance of conceptualising 'spread' for classifying invasiveness and 'impact' for management. Finally, we propose a protocol for classifying populations based on (i) dispersal mechanism, (ii) species origin, (iii) population status, and (iv) impact. Collectively and without introducing new terminology, the framework that we present aims to facilitate effective communication and collaboration in invasion science and management of non-native species.

Description: The zebra mussel Dreissena polymorpha is one of the most successful, notorious, and detrimental aquatic invasive non-native species worldwide, having invaded Europe and North America while causing substantial ecological and socio-economic impacts. Here, we investigated the spatiotemporal trends in this species' invasion success using 178 macroinvertebrate abundance time series, containing 1451 records of D. polymorpha collected across nine European countries between 1972–2019. Using these raw (absolute) abundance data, we examined trends and drivers of occurrences and relative abundances of D. polymorpha within invaded communities. Meta-regression models revealed non-significant trends both at the European level and for the majority of the invaded countries, except for France (significant decreasing trend) and Hungary (marginally positive trend). At the European level, the number of D. polymorpha occurrences over time followed a flat-top bell-shaped distribution, with a steep increase between 1973–1989 followed by a plateau phase prior to significantly declining post-1998. Using a series of climatic and hydromorphological site-specific characteristics of invaded and uninvaded sites from two periods (1998–2002; 2011–2015), we found that native richness, non-native abundance, distance to the next barrier, and elevation were associated with the occurrence of D. polymorpha. We also found that higher native richness and lower latitude were related to lower relative abundances. Using Cohen's D as a measure of D. polymorpha impact, we found that biodiversity within the invaded sites was initially higher than in uninvaded ones, but then declined, suggesting differences in biodiversity trends across invaded and uninvaded sites. While our results emphasise the high invasion success of D. polymorpha, increasing stressors within the context of global change – particularly ongoing climate change – are likely to enhance invasion rates and the impact of D. polymorpha in the near future, exacerbated by the lack of timely and effective management actions.

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: Promoting effects of aluminum addition on chlorophyll biosynthesis and growth of two cultured iron‐limited marine diatoms Linbin Zhou CAS Key Laboratory of Tropical Marine Bio‐resources and Ecology, Guangdong Provincial Key Laboratory of Applied Marine Biology, South China Sea Institute of Oceanology Chinese Academy of Sciences Guangzhou China Marine Biogeochemistry Division GEOMAR Helmholtz Centre for Ocean Research Kiel Germany University of Chinese Academy of Sciences Beijing China https://orcid.org/0000-0001-7230-4116 Fengjie Liu Marine Biogeochemistry Division GEOMAR Helmholtz Centre for Ocean Research Kiel Germany Grantham Institute—Climate Change and the Environment, Department of Life Sciences Imperial College London London UK Eric P. Achterberg Marine Biogeochemistry Division GEOMAR Helmholtz Centre for Ocean Research Kiel Germany Anja Engel Marine Biogeochemistry Division GEOMAR Helmholtz Centre for Ocean Research Kiel Germany https://orcid.org/0000-0002-1042-1955 Peter G.C. Campbell Institut National de la Recherche Scientifique Centre Eau Terre Environnement Quebec Canada https://orcid.org/0000-0001-7160-4571 Claude Fortin Institut National de la Recherche Scientifique Centre Eau Terre Environnement Quebec Canada https://orcid.org/0000-0002-2479-1869 Liangmin Huang CAS Key Laboratory of Tropical Marine Bio‐resources and Ecology, Guangdong Provincial Key Laboratory of Applied Marine Biology, South China Sea Institute of Oceanology Chinese Academy of Sciences Guangzhou China University of Chinese Academy of Sciences Beijing China Yehui Tan CAS Key Laboratory of Tropical Marine Bio‐resources and Ecology, Guangdong Provincial Key Laboratory of Applied Marine Biology, South China Sea Institute of Oceanology Chinese Academy of Sciences Guangzhou China University of Chinese Academy of Sciences Beijing China Abstract Aluminum (Al) may play a role in the ocean's capacity for absorbing atmospheric CO 2 via influencing carbon fixation, export, and sequestration. Aluminum fertilization, especially in iron (Fe)‐limited high‐nutrient, low‐chlorophyll ocean regions, has been proposed as a potential CO 2 removal strategy to mitigate global warming. However, how Al addition would influence the solubility and bioavailability of Fe as well as the physiology of Fe‐limited phytoplankton has not yet been examined. Here, we show that Al addition (20 and 100 nM) had little influence on the Fe solubility in surface seawater and decreased the Fe bio‐uptake by 11–22% in Fe‐limited diatom Thalassiosira weissflogii in Fe‐buffered media. On the other hand, the Al addition significantly increased the rate of chlorophyll biosynthesis by 45–60% for Fe‐limited T. weissflogii and 81–102% for Fe‐limited Thalassiosira pseudonana , as well as their cell size, cellular chlorophyll content, photosynthetic quantum efficiency ( F v / F m ) and growth rate. Under Fe‐sufficient conditions, the Al addition still led to an increased growth rate, though the beneficial effects of Al addition on chlorophyll biosynthesis were no longer apparent. These results suggest that Al may facilitate chlorophyll biosynthesis and benefit the photosynthetic efficiency and growth of Fe‐limited diatoms. We speculate that Al addition may enhance intracellular Fe use efficiency for chlorophyll biosynthesis by facilitating the superoxide‐mediated intracellular reduction of Fe(III) to Fe(II). Our study provides new evidence and support for the iron–aluminum hypothesis.

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 simulation of deep-sea conditions in laboratories is technically challenging but necessary for experiments that aim at a deeper understanding of physiological mechanisms or host-symbiont interactions of deep-sea organisms. In a proof-of-concept study, we designed a recirculating system for long-term culture (〉2 yr) of deep-sea mussels Gigantidas childressi (previously Bathymodiolus childressi). Mussels were automatically (and safely) supplied with a maximum stable level of ~60 μmol L−1 methane in seawater using a novel methane–air mixing system. Experimental animals also received daily doses of live microalgae. Condition indices of cultured G. childressi remained high over the years, and low shell growth rates could be detected, too, which is indicative of positive energy budgets. Using stable isotope data, we demonstrate that G. childressi in our culture system gained energy, both, from the digestion of methane-oxidizing endosymbionts and from digesting particulate food (microalgae). Limitations of the system, as well as opportunities for future experimental approaches involving deep-sea mussels, are discussed.

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: Infections by filamentous phages, which are usually nonlethal to the bacterial cells, influence bacterial fitness in various ways. While phage-encoded accessory genes, for example virulence genes, can be highly beneficial, the production of viral particles is energetically costly and often reduces bacterial growth. Consequently, if costs outweigh benefits, bacteria evolve resistance, which can shorten phage epidemics. Abiotic conditions are known to influence the net-fitness effect for infected bacteria. Their impact on the dynamics and trajectories of host resistance evolution, however, remains yet unknown. To address this, we experimentally evolved the bacterium Vibrio alginolyticus in the presence of a filamentous phage at three different salinity levels, that is (1) ambient, (2) 50% reduction and (3) fluctuations between reduced and ambient. In all three salinities, bacteria rapidly acquired resistance through super infection exclusion (SIE), whereby phage-infected cells acquired immunity at the cost of reduced growth. Over time, SIE was gradually replaced by evolutionary fitter surface receptor mutants (SRM). This replacement was significantly faster at ambient and fluctuating conditions compared with the low saline environment. Our experimentally parameterized mathematical model explains that suboptimal environmental conditions, in which bacterial growth is slower, slow down phage resistance evolution ultimately prolonging phage epidemics. Our results may explain the high prevalence of filamentous phages in natural environments where bacteria are frequently exposed to suboptimal conditions and constantly shifting selections regimes. Thus, our future ocean may favour the emergence of phage-born pathogenic bacteria and impose a greater risk for disease outbreaks, impacting not only marine animals but also humans.

Description: A sense of kuleana (personal responsibility) in caring for the land and sea. An appreciation for laulima (many hands cooperating). An understanding of aloha 'āina (love of the land). The University of Hawai'i at Manoa hosted the 2023 Ecological Dissertations in Aquatic Sciences (Eco-DAS) program, which fostered each of these intentions by bringing together a team of early career aquatic ecologists for a week of networking and collaborative, interdisciplinary project development

Description: The spectral composition of light is an important factor for the metabolism of photosynthetic organisms. Several blue light-regulated metabolic processes have already been identified in the industrially relevant microalga Monoraphidium braunii. However, little is known about the spectral impact on this species' growth, fatty acid (FA), and pigment composition. In this study, M. braunii was cultivated under different light spectra (white light: 400–700 nm, blue light: 400–550 nm, green light: 450–600 nm, and red light: 580–700 nm) at 25°C for 96 h. The growth was monitored daily. Additionally, the FA composition, and pigment concentration was analyzed after 96 h. The highest biomass production was observed upon white light and red light irradiation. However, green light also led to comparably high biomass production, fueling the scientific debate about the contribution of weakly absorbed light wavelengths to microalgal biomass production. All light spectra (white, blue, and green) that comprised blue-green light (450–550 nm) led to a higher degree of FA unsaturation and a greater concentration of all identified pigments than red light. These results further contribute to the growing understanding that blue-green light is an essential trigger for maximized pigment concentration and FA unsaturation in green microalgae.

Description: This article presents risk factors that are associated with the handling of unexploded ordnance (UXO) during explosive ordnance disposal (EOD) operations in German waters. The construction of offshore wind parks and the German immediate action program are expected to increase the number of EOD operations. Existing literature and guidelines do not offer a structured and reproducible framework for assessing EOD risk. To fill this gap, a network of EOD risk factors was developed by means of a literature review and validation via expert consultation. The study was scoped to “personnel and equipment at the EOD location” as the risk receptor and “undesired detonation” as the undesired event under investigation. Factors are subdivided into UXO factors that depend on the object that should be handled and factors that describe the object's surrounding environment. While the former can be researched by an EOD expert, the latter must be measured on site or acquired from a model. Each of these factors contributes to risk, some directly and others indirectly via other factors. The complexity of the resulting network, with its 33 factors, demonstrates the need for a reliable and reproducible model to quantify EOD risk. Its purpose is not to replace EOD experts but to aid them in their decision‐making process. Such a tool can provide valuable support for the high‐cost and high‐risk EOD operations.

Description: The collection of zooplankton swimmers and sinkers in time-series sediment traps provides unique insight into year-round and interannual trends in zooplankton population dynamics. These samples are particularly valuable in remote and difficult to access areas such as the Arctic Ocean, where samples from the ice-covered season are rare. In the present study, we investigated zooplankton composition based on swimmers and sinkers collected by sediment traps at water depths of 180–280, 800–1320, and 2320–2550 m, over a period of 16 yr (2000–2016) at the Long-Term Ecological Research observatory HAUSGARTEN located in the eastern Fram Strait (79°N, 4°E). The time-series data showed seasonal and interannual trends within the dominant zooplankton groups including copepoda, foraminifera, ostracoda, amphipoda, pteropoda, and chaetognatha. Amphipoda and copepoda dominated the abundance of swimmers while pteropoda and foraminifera were the most important sinkers. Although the seasonal occurrence of these groups was relatively consistent between years, there were notable interannual variations in abundance, suggesting the influence of various environmental conditions such as sea-ice dynamic and lateral advection of water masses, for example, meltwater and Atlantic water. Statistical analyses revealed a correlation between the Arctic dipole climatic index and sea-ice dynamics (i.e., ice coverage and concentration), as well as the importance of the distance from the ice edge on swimmer composition patterns and carbon export.

Description: Communities and their functioning are jointly shaped by ecological and evolutionary processes that manifest in diversity shifts of their component species and genotypes. How both processes contribute to community functional change over time is rarely studied. We here repeatedly quantified eco-evolutionary contributions to CO2-driven total abundance and mean cell size changes after short-, mid-, and longer-term (80, 168, and 〉 168 d, respectively) in experimental phytoplankton communities. While the CO2-driven changes in total abundance and mean size in the short- and mid-term could be predominantly attributed to ecological shifts, the relative contribution of evolution increased. Over the longer-term, the CO2-effect and underlying eco-evolutionary changes disappeared, while total abundance increased, and mean size decreased significantly independently of CO2. The latter could be presumably attributed to CO2-independent genotype selection which fed back to species composition. In conclusion, ecological changes largely dominated the regulation of environmentally driven phytoplankton functional shifts at first. However, evolutionary changes gained importance with time, and can ultimately feedback on species composition, and thus must be considered when predicting phytoplankton change.

Description: Massive swarms of the red crab Pleuroncodes planipes (Stimpson, 1860), a species of squat lobster, are a dominant functional component of the upwelling ecosystem in the eastern Pacific Ocean (Boyd, 1967; Smith et al., 1975). These swarms can wash ashore on the coast, creating mass depositions of crustacean carcasses, a striking phenomenon that has been long documented in Baja California and California (Aurioles-Gamboa et al., 1994; Boyd, 1967). However, little is known about the fate of crab swarms transported offshore by oceanic currents. In May 2015, using an autonomous deep-sea robot, we discovered an unexpectedly large fall of red crab carcasses (〉1000 carcasses ha−1) at a depth of 4050 m on the abyssal Pacific seafloor (Figure 1), almost 1500 km from their spawning areas off the northwest American coast. Several questions arise from this unexpected finding that may help unveil additional close linkages in nutritional transport between processes at the sea surface and the remote abyssal seafloor.

Description: End-member modelling of bulk grain-size distributions allows the unravelling of natural and anthropogenic depositional processes in salt marshes and quantification of their respective contribution to marsh accretion. The sedimentology of two marshes is presented: (1) a sheltered back-barrier marsh; and (2) an exposed, reinstated foreland marsh. Sedimentological data are supplemented by an age model based on lead-210 decay and caesium-137, as well as geochemical data. End-member modelling of grain-size data shows that marsh growth in back-barrier settings is primarily controlled by the settling of fines from suspension during marsh inundation. In addition, nearby active dunes deliver aeolian sediment (up to 77% of the total sediment accretion), potentially enhancing the capability of salt marshes to adapt to sea-level rise. Growth of exposed marshes, by contrast, primarily results from high-energy inundation and is attributed to two sediment-transport processes. On the seaward edge of the marsh, sedimentation is dominated by coarser-grained traction load, whereas further inland, settling of fine-grained suspension load prevails. In addition, a third, coarse-grained sediment sub-population is interpreted to derive from anthropogenic land-reclamation measures, that is material from drainage channels relocated onto the marsh surface. This process contributed up to 34% to the total marsh accretion and terminated synchronously with the end of land reclamation measures. Data suggest that natural sediment supply to marshes alone is sufficient to outpace contemporary sea-level rise in the study area. This underlines the resilience potential of salt marshes in times of rising sea levels. The comparison of grain-size sub-populations with observed climate variability implies that even managed marshes allow for the extraction of environmental signals if natural and anthropogenic sedimentary processes are determined and their relative contribution to bulk sediment composition is quantified. Data series based solely on bulk sediments, however, seem to be of limited use because it is difficult to exclude bias of natural signals by anthropogenic measures.

Description: Driven by climate change, marine biodiversity is undergoing a phase of rapid change that has proven to be even faster than changes observed in terrestrial ecosystems. Understanding how these changes in species composition will affect future marine life is crucial for conservation management, especially due to increasing demands for marine natural resources. Here, we analyse predictions of a multiparameter habitat suitability model covering the global projected ranges of 〉33,500 marine species from climate model projections under three CO2 emission scenarios (RCP2.6, RCP4.5, RCP8.5) up to the year 2100. Our results show that the core habitat area will decline for many species, resulting in a net loss of 50% of the core habitat area for almost half of all marine species in 2100 under the high-emission scenario RCP8.5. As an additional consequence of the continuing distributional reorganization of marine life, gaps around the equator will appear for 8% (RCP2.6), 24% (RCP4.5), and 88% (RCP8.5) of marine species with cross-equatorial ranges. For many more species, continuous distributional ranges will be disrupted, thus reducing effective population size. In addition, high invasion rates in higher latitudes and polar regions will lead to substantial changes in the ecosystem and food web structure, particularly regarding the introduction of new predators. Overall, our study highlights that the degree of spatial and structural reorganization of marine life with ensued consequences for ecosystem functionality and conservation efforts will critically depend on the realized greenhouse gas emission pathway.

Description: 1. In a warming world, both the metabolic rates of ectotherm predators and the phenology of their prey organisms is subject to change. Knowledge on how intrinsic and extrinsic factors govern predator-prey interactions is essential in order to understand how the environment regulates the vital rates of consumers. Controlled experiments, however, simultaneously testing behavioural and growth responses of the larvae of fish and other ectotherm organisms in different feeding regimes are scarce. 2. Prey size (PS) selection was determined for young Atlantic herring (Clupea harengus L.) larvae offered 100- to 850-μm copepods (Acartia tonsa) at five different concentrations. In separate, 4- (13°C) or 7-day (7°C) trials, the effect of prey size on larval foraging behaviour, specific growth rate (SGR) and biochemical condition (RNA:DNA, RD, a proxy for individual instantaneous growth) was tested. 3. Preferred (selected) PS was similar at all prey concentrations but increased from 3 to 5% predator length with increasing larval size. At various temperatures, dome-shaped relationships existed between PS and larval RD (and accordingly SGR). Compensatory changes in foraging behaviour (pause and feeding strike frequencies) existed but were not adequate to maintain positive SGR when available prey were substantially smaller than those preferred by larvae. 5. A physiology-based model predicted that larvae depended more heavily on optimal prey sizes at the colder versus warmer temperature to grow well and that the profitable prey niche breadth (the range in prey sizes in which growth was positive) increased at warmer temperatures. 6. Seemingly subtle match-mismatch dynamics between ectotherm predators and their preferred prey size can have large, temperature-dependent consequences for rates of growth and likely survival of the predator. To the best knowledge, this was the first study to directly quantify the “costs and trade-offs” of optimal foraging in marine fish larvae.

Description: Dykes are magma-filled fractures propagating through the brittle crust. Understanding the physics of dyking process is essential to mitigate the volcanic hazard associated with the opening of new eruptive fissures at the surface. Often, physics-based models view either fracturing of the host rock or viscous flow of the magma as the dominating energy sink during dyke propagation. Here, we provide a numerical model that captures the coupling of fracturing at the crack tip and the transport of a viscous fluid. Built with the boundary element technique, our model allows for computation of the shape and velocity of a growing fluid-filled crack accounting for the viscosity of the fluid: The fluid flow induces a viscous pressure drop acting at the crack walls, and modifies the shape of the crack. The energy conservation equation provides the constraints to solve for the crack growth velocity, assuming that brittle fracturing and viscous flow are the main processes that dissipate energy. Using a parameter range that represents typical magmatic intrusions, we obtain crack shapes displaying some typical characteristics, including a tear-drop head and an open tail that depend on rock rigidity, magma viscosity and buoyancy. We show that viscous forces significantly contribute to the energy dissipated during the propagation of magmatic dykes. Applied to the 1998 intrusion at Piton de la Fournaise (La Réunion Island), we provide ranges of dyke lengths and openings by adjusting the numerical velocity to the one deduced from the migration of volcano-tectonic events. Key Points We present a new modeling scheme to compute the shape and velocity of a growing fluid-filled crack Our magmatic dykes show a tear drop head and open tail, on a wide range of propagation velocities We reproduce the velocity and fit important parameters for the 1998 Piton de la Fournaise intrusion

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: Laser ablation multi-collector-inductively coupled plasma-mass spectrometry (LA-MC-ICP-MS) has become a valuable tool for the in situ measurement of the boron isotope composition of geological samples at high (tens to hundreds of μm) spatial resolution. That said, this application suffers from significant analytical challenges. We focus in this study on the underlying processes of two of the main causes for inaccuracies using this technique. We provide empirical evidence that not only Ca ions (Sadekov et al. 2019, Standish et al. 2019, Evans et al. 2019) but also Ar ions, that are reflected within the flight tube of the mass spectrometer, are the source for previously reported issues with spectral baselines. We also address the impact of plasma conditions on the instrumental mass fractionation as a source for matrix- and mass-load-related analytical biases. Comparing experimental data with the results of a dedicated release and diffusion model (RDM) we estimate that a close to complete (~ 97%) release of boron from the sample aerosol is needed to allow for consistently accurate LA boron isotope measurement results without the need for corrections.

Description: Range expansions can lead to increased contact of divergent populations, thus increasing the potential of hybridization events. Whether viable hybrids are produced will most likely depend on the level of genomic divergence and associated genomic incompatibilities between the different entities as well as environmental conditions. By taking advantage of historical Baltic cod (Gadus morhua) otolith samples combined with genotyping and whole genome sequencing, we here investigate the genetic impact of the increased spawning stock biomass of the eastern Baltic cod stock in the mid 1980s. The eastern Baltic cod is genetically highly differentiated from the adjacent western Baltic cod and locally adapted to the brackish environmental conditions in the deeper Eastern basins of the Baltic Sea unsuitable for its marine counterparts. Our genotyping results show an increased proportion of eastern Baltic cod in western Baltic areas (Mecklenburg Bay and Arkona Basin)—indicative of a range expansion westwards—during the peak population abundance in the 1980s. Additionally, we detect high frequencies of potential hybrids (including F1, F2 and backcrosses), verified by whole genome sequencing data for a subset of individuals. Analysis of mitochondrial genomes further indicates directional gene flow from eastern Baltic cod males to western Baltic cod females. Our findings unravel that increased overlap in distribution can promote hybridization between highly divergent populations and that the hybrids can be viable and survive under specific and favourable environmental conditions. However, the observed hybridization had seemingly no long-lasting impact on the continuous separation and genetic differentiation between the unique Baltic cod stocks.

Description: Genome sequencing enables answering fundamental questions about the genetic basis of adaptation, population structure and epigenetic mechanisms. Yet, we usually need a suitable reference genome for mapping population-level resequencing data. In some model systems, multiple reference genomes are available, giving the challenging task of determining which reference genome best suits the data. Here, we compared the use of two different reference genomes for the three-spined stickleback (Gasterosteus aculeatus), one novel genome derived from a European gynogenetic individual and the published reference genome of a North American individual. Specifically, we investigated the impact of using a local reference versus one generated from a distinct lineage on several common population genomics analyses. Through mapping genome resequencing data of 60 sticklebacks from across Europe and North America, we demonstrate that genetic distance among samples and the reference genomes impacts downstream analyses. Using a local reference genome increased mapping efficiency and genotyping accuracy, effectively retaining more and better data. Despite comparable distributions of the metrics generated across the genome using SNP data (i.e. π, Tajima's D and FST), window-based statistics using different references resulted in different outlier genes and enriched gene functions. A marker-based analysis of DNA methylation distributions had a comparably high overlap in outlier genes and functions, yet with distinct differences depending on the reference genome. Overall, our results highlight how using a local reference genome decreases reference bias to increase confidence in downstream analyses of the data. Such results have significant implications in all reference-genome-based population genomic analyses.

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: An increasingly globalised world has facilitated the movement of non-native species (NNS) via the poorly regulated international pet trade. While focus is increasingly being placed on preventative action to combat invasive NNS—often cheaper and less difficult than the management of established populations—successful prevention requires controlling potential pathways and obtaining baseline knowledge of species' availability. Here we performed an in-depth analysis of the freshwater pet trade as one major vector of NNS, compiling its species inventory and deriving threats of NNS release and establishment in the wild. With Germany as our study region, we surveyed pet stores, websites and the country's largest online classified portal, eBay Kleinanzeigen, recording the taxa encountered. For each species, we determined the likelihood of release based on availability and price (cheaper and/or more readily available species have been shown to be of greater risk), and the likelihood of establishment based on ecological niche breadth and niche overlap with environmental conditions in Germany. The survey revealed 669 species, of which 651 were non-native to Germany. Looking at release likelihood, more readily available species in pet stores and on websites proved to be cheaper. For websites, there was a significant effect of occurrence status (i.e. released, not released, native) on price, with released and native species being significantly cheaper. Species previously released in Germany and elsewhere demonstrated greater niche breadths and greater niche overlaps between their source regions and Germany; and for species released in Germany, there was a significantly positive relationship between the magnitude of niche overlap and the number of documented occurrences. Finally, we combined our release and establishment likelihood findings under ‘Release Risk’ metrics to highlight the species most worthy of prioritisation. We propose these metrics as proactive methods for screening species in the trade, which can inform future policy direction and intervention. Read the free Plain Language Summary for this article on the Journal blog.

Description: Compound‐specific isotope analysis has opened up a new realm for resolving the sources and transformation processes of marine organic matter. However, the stable carbon isotope patterns of amino sugars remain unknown. We examined δ 13 C of amino sugars in marine phytoplankton and heterotrophic bacteria, and the variations in amino sugar δ 13 C during 66‐d planktonic organic matter degradation experiments, to investigate the metabolic sources and transformations of amino sugars by bacterial reworking. The δ 13 C values of glucosamine (GlcN) and galactosamine (GalN) were comparable in heterotrophic bacteria (difference Δδ 13 C GlcN–GalN = 0.4–4.0‰) but pronouncedly different in phytoplankton (Δδ 13 C GlcN–GalN = 4.3–16.6‰), suggesting similar synthesis pathways of GlcN and GalN in bacteria that differed from phytoplankton. Compared to GlcN and GalN, bacteria preferentially use isotopically light organic compounds for muramic acid (MurA) synthesis. During simulated microbial degradation of organic matter, the δ 13 C difference between GlcN and GalN decreased from 5.8‰ on the initial day to 1‰ at a late stage in the experiment, but the difference between GlcN and MurA remained at 5.3‰. This difference is consistent with the pattern in cultured phytoplankton (average Δδ 13 C GlcN–GalN = 5.9‰ ± 1.4‰) and heterotrophic bacteria (average Δδ 13 C GlcN–MurA = 4.6‰ ± 3.4‰), indicating enhanced bacterial resynthesis as degradation proceeded. Based on the difference in δ 13 C among GlcN, GalN, and MurA, we propose a novel index of variation in amino sugar δ 13 C, representing amino sugar resynthesis, to describe the diagenetic state of organic matter. Together, these findings suggest that amino sugar δ 13 C can be used as a new tool to track heterotrophic processes of marine organic matter.

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: Understanding individual growth in commercially exploited fish populations is key to successful stock assessment and informed ecosystem-based fisheries management. Traditionally, growth rates in marine fish are estimated using otolith age-readings in combination with age-length relationships from field samples, or tag-recapture field experiments. However, for some species, otolith-based approaches have been proven unreliable and tag-recapture experiments suffer from high working effort and costs as well as low recapture rates. An important alternative approach for estimating fish growth is represented by bioenergetic modelling which in addition to pure growth estimation can provide valuable insights into the processes leading to temporal growth changes resulting from environmental and related behavioural changes. We here developed an individual-based bioenergetic model for Western Baltic cod (Gadus morhua), traditionally a commercially important fish species that however collapsed recently and likely suffers from climate change effects. Western Baltic cod is an ideal case study for bioenergetic modelling because of recently gained in-situ process knowledge on spatial distribution and feeding behaviour based on highly resolved data on stomachs and fish distribution. Additionally, physiological processes such as gastric evacuation, consumption, net-conversion efficiency and metabolic rates have been well studied for cod in laboratory experiments. Our model reliably reproduced seasonal growth patterns observed in the field. Importantly, our bioenergetic modelling approach implementing depth-use patterns and food intake allowed us to explain the potentially detrimental effect summer heat periods have on the growth of Western Baltic cod that likely will increasingly occur in the future. Hence, our model simulations highlighted a potential mechanism on how warming due to climate change affects the growth of a key species that may apply for similar environments elsewhere.

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: Pocillopora damicornis (Linnaeus, 1758), a species complex, consists of several genetic lineages, some of which likely represent reproductively isolated species, including the species Pocillopora acuta Lamarck, 1816. Pocillopora acuta can exhibit similar morphological characteristics as P. damicornis, thus making it difficult to identify species-level taxonomic units. To determine whether the P. damicornis-like colonies on the reefs in the Andaman Sea (previously often identified as P. damicornis) consist of different species, we sampled individual colonies at five sites along a 50 km coastal stretch at Phuket Island and four island sites towards Krabi Province, Thailand. We sequenced 210 coral samples for the mitochondrial open reading frame and identified six distinct haplotypes, all belonging to P. acuta according to the literature. Recently, P. acuta was observed to efficiently recolonize heat-damaged reefs in Thailand as well as globally, making it a potentially important coral species in future reefs. Specifically in the light of global change, this study underscores the importance of high-resolution molecular species recognition, since taxonomic units are important factors for population genetic studies, and the latter are crucial for management and conservation efforts.

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: The relevance of considering environmental variability for understanding and predicting biological responses to environmental changes has resulted in a recent surge in variability-focused ecological research. However, integration of findings that emerge across studies and identification of remaining knowledge gaps in aquatic ecosystems remain critical. Here, we address these aspects by: (1) summarizing relevant terms of variability research including the components (characteristics) of variability and key interactions when considering multiple environmental factors; (2) identifying conceptual frameworks for understanding the consequences of environmental variability in single and multifactorial scenarios; (3) highlighting challenges for bridging theoretical and experimental studies involving transitioning from simple to more complex scenarios; (4) proposing improved approaches to overcome current mismatches between theoretical predictions and experimental observations; and (5) providing a guide for designing integrated experiments across multiple scales, degrees of control, and complexity in light of their specific strengths and limitations.

Description: Aim: Invasive alien species are a growing problem worldwide due to their ecological, economic and human health impacts. The “killer shrimp” Dikerogammarus villosus is a notorious invasive alien amphipod from the Ponto-Caspian region that has invaded many fresh and brackish waters across Europe. Understandings of large-scale population dynamics of highly impactful invaders such as D. villosus are lacking, inhibiting predictions of impact and efficient timing of management strategies. Hence, our aim was to assess trends and dynamics of D. villosus as well as its impacts in freshwater rivers and streams. Location: Europe. Methods: We analysed 96 European time series between 1994 and 2019 and identified trends in the relative abundance (i.e. dominance %) of D. villosus in invaded time series, as well as a set of site-specific characteristics to identify drivers and determinants of population changes and invasion dynamics using meta-regression modelling. We also looked at the spread over space and time to estimate the invasion speed (km/year) of D. villosus in Europe. We investigated the impact of D. villosus abundance on recipient community metrics (i.e. abundance, taxa richness, temporal turnover, Shannon diversity and Pielou evenness) using generalized linear models. Results: Population trends varied across the time series. Nevertheless, community dominance of D. villosus increased over time across all time series. The frequency of occurrences (used as a proxy for invader spread) was well described by a Pareto distribution, whereby we estimated a lag phase (i.e. the time between introduction and spatial expansion) of approximately 28 years, followed by a gradual increase before new occurrences declined rapidly in the long term. D. villosus population change was associated with decreased taxa richness, community turnover and Shannon diversity. Main Conclusion: Our results show that D. villosus is well-established in European waters and its abundance significantly alters ecological communities. However, the multidecadal lag phase prior to observed spatial expansion suggests that initial introductions by D. villosus are cryptic, thus signalling the need for more effective early detection methods.

Description: Background: Syngnathids are a highly derived and diverse fish clade comprising the pipefishes, pipe-horses, and seahorses. They are characterized by a plethora of iconic traits that increasingly capture the attention of biologists, including geneticists, ecologists, and developmental biologists. The current understanding of the origins of their derived body plan is, however, hampered by incomplete and limited descriptions of the early syngnathid ontogeny. Results: We provide a comprehensive description of the development of Nerophis ophidion, Syngnathus typhle, and Hippocampus erectus from early cleavage stages to release from the male brooding organ and beyond, including juvenile development. We comparatively describe skeletogenesis with a particular focus on dermal bony plates, the snout-like jaw morphology, and appendages. Conclusions: This most comprehensive and detailed account of syngnathid development to date suggests that convergent phenotypes (e.g., reduction and loss of the caudal fins), likely arose by distinct ontogenetic means in pipefishes and seahorses. Comparison of the ontogenetic trajectories of S. typhle and H. erectus provides indications that characteristic features of the seahorse body plan result from developmental truncation. Altogether, this work provides a valuable resource and framework for future research to understand the evolution of the outlandish syngnathid morphology from a developmental perspective.

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

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

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

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

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

Description: Aquaculture has been one of the fastest-growing food production systems sectors for over three decades. With its growth, the demand for alternative, cheaper and high-quality feed ingredients is also increasing. Innovation investments on providing new functional feed alternatives have yielded several viable alternative raw materials. Considering all the current feed ingredients, their circular adaption in the aquafeed manufacturing industry is clearly of the utmost importance to achieve sustainable aquaculture in the near future. The use of terrestrial plant materials and animal by-products predominantly used in aquafeed ingredients puts a heavily reliance on terrestrial agroecosystems, which also has its own sustainability concerns. Therefore, the aquafeed industry needs to progress with functional and sustainable alternative raw materials for feed that must be more resilient and consistent, considering a circular perspective. In this review, we assess the current trends in using various marine organisms, ranging from microorganisms (including fungi, thraustochytrids, microalgae and bacteria) to macroalgae and macroinvertebrates as viable biological feed resources. This review focuses on the trend of circular use of resources and the development of new value chains. In this, we present a perspective of promoting novel circular economy value chains that promote the re-use of biological resources as valuable feed ingredients. Thus, we highlight some potentially important marine-derived resources that deserve further investigations for improving or addressing circular aquaculture.

Description: Invasive species can successfully and rapidly colonize new niches and expand ranges via founder effects and enhanced tolerance towards environmental stresses. However, the underpinning molecular mechanisms (i.e., gene expression changes) facilitating rapid adaptation to harsh environments are still poorly understood. The red seaweed Gracilaria vermiculophylla, which is native to the northwest Pacific but invaded North American and European coastal habitats over the last 100 years, provides an excellent model to examine whether enhanced tolerance at the level of gene expression contributed to its invasion success. We collected G. vermiculophylla from its native range in Japan and from two non-native regions along the Delmarva Peninsula (Eastern United States) and in Germany. Thalli were reared in a common garden for 4 months at which time we performed comparative transcriptome (mRNA) and microRNA (miRNA) sequencing. MRNA-expression profiling identified 59 genes that were differently expressed between native and non-native thalli. Of these genes, most were involved in metabolic pathways, including photosynthesis, abiotic stress, and biosynthesis of products and hormones in all four non-native sites. MiRNA-based target-gene correlation analysis in native/non-native pairs revealed that some target genes are positively or negatively regulated via epigenetic mechanisms. Importantly, these genes are mostly associated with metabolism and defence capability (e.g., metal transporter Nramp5, senescence-associated protein, cell wall-associated hydrolase, ycf68 protein and cytochrome P450-like TBP). Thus, our gene expression results indicate that resource reallocation to metabolic processes is most likely a predominant mechanism contributing to the range-wide persistence and adaptation of G. vermiculophylla in the invaded range. This study, therefore, provides molecular insight into the speed and nature of invasion-mediated rapid adaption.

Description: Clonal reproduction, the formation of nearly identical individuals via mitosis in the absence of genetic recombination, is a very common reproductive mode across plants, fungi and animals. To detect clonal genetic structure, genetic similarity indices based on shared alleles are widely used, such as the Jaccard index, or identity by state. Here we propose a new pairwise genetic similarity index, the SH index, based on segregating genetic marker loci (typically single nucleotide polymorphisms) that are identically heterozygous for pairs of samples (NSH). To test our method, we analyse two old seagrass clones (Posidonia australis, estimated to be around 8500 years old; Zostera marina, 〉750 years old) along with two young Z. marina clones of known age (17 years old). We show that focusing on shared heterozygosity amplifies the power to distinguish sample pairs belonging to different clones compared to methods focusing on all shared alleles. Our proposed workflow can successfully detect clonemates at a location dominated by a single clone. When the collected samples involve two or more clones, the SH index shows a clear gap between clonemate pairs and interclone sample pairs. Ideally NSH should be on the order of approximately ≥3000, a number easily achievable via restriction-site associated DNA (RAD) sequencing or whole-genome resequencing. Another potential application of the SH index is to detect possible parent–descendant pairs under selfing. Our proposed workflow takes advantage of the availability of the larger number of genetic markers in the genomic era, and improves the ability to distinguish clonemates from nonclonemates in multicellular diploid clonal species.

Description: Evolutionary novelties—derived traits without clear homology found in the ancestors of a lineage—may promote ecological specialization and facilitate adaptive radiations. Examples for such novelties include the wings of bats, pharyngeal jaws of cichlids and flowers of angiosperms. Belonoid fishes (flying fishes, halfbeaks and needlefishes) feature an astonishing diversity of extremely elongated jaw phenotypes with undetermined evolutionary origins. We investigate the development of elongated jaws in a halfbeak (Dermogenys pusilla) and a needlefish (Xenentodon cancila) using morphometrics, transcriptomics and in situ hybridization. We confirm that these fishes' elongated jaws are composed of distinct base and novel ‘extension’ portions. These extensions are morphologically unique to belonoids, and we describe the growth dynamics of both bases and extensions throughout early development in both studied species. From transcriptomic profiling, we deduce that jaw extension outgrowth is guided by populations of multipotent cells originating from the anterior tip of the dentary. These cells are shielded from differentiation, but proliferate and migrate anteriorly during the extension's allometric growth phase. Cells left behind at the tip leave the shielded zone and undergo differentiation into osteoblast-like cells, which deposit extracellular matrix with both bone and cartilage characteristics that mineralizes and thereby provides rigidity. Such bone has characteristics akin to histological observations on the elongated ‘kype’ process on lower jaws of male salmon, which may hint at common conserved regulatory underpinnings. Future studies will evaluate the molecular pathways that govern the anterior migration and proliferation of these multipotent cells underlying the belonoids' evolutionary novel jaw extensions.

Description: Foraging is a behavioural process and, therefore, individual behaviour and diet are theorized to covary. However, few comparisons of individual behaviour type and diet exist in the wild. We tested whether behaviour type and diet covary in a protected population of Atlantic cod, Gadus morhua. Working in a no-take marine reserve, we could collect data on natural behavioural variation and diet choice with minimal anthropogenic disturbance. We inferred behaviour using acoustic telemetry and diet from stable isotope compositions (expressed as δ13C and δ15N values). We further investigated whether behaviour and diet could have survival costs. We found cod with shorter diel vertical migration distances fed at higher trophic levels. Cod δ13C and δ15N values scaled positively with body size. Neither behaviour nor diet predicted survival, indicating phenotypic diversity is maintained without survival costs for cod in a protected ecosystem. The links between diet and diel vertical migration highlight that future work is needed to understand whether the shifts in this behaviour during environmental change (e.g. fishing or climate), could lead to trophic cascades.

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: In regions of the nitrogen limited low latitude ocean, phosphate can also be depleted to levels initiating stress responses in marine microbes. Here, we associate a broad region of phosphate depletion in the subtropical North Pacific with different levels of phosphorus stress. Nutrient and aerosol addition experiments demonstrated primary nitrogen limitation of the bulk phytoplankton community, with supply of aerosols relieving this limitation. At northern sites with depleted phosphate, alkaline phosphatase activities were enhanced, indicating elevated phosphorus stress. Analysis of satellite- and model-derived aerosol loading showed that aerosol deposition was elevated in these regions. Surface rate measurements suggested that the regional enhancement in phosphate depletion was predominantly driven by elevated nitrogen fixation, likely stimulated by the coincident supply of aerosol iron. Such observations are important for predicting future biogeochemical responses in the subtropical North Pacific to changing aerosol supply.

Description: Methane (CH4) is a climate-relevant trace gas that is emitted from the open and coastal oceans in considerable amounts. However, its distribution in remote oceanic areas is largely unknown. To fill this knowledge gap, dissolved CH4 was measured at nine stations at 75°S in the Ross Sea during austral summer in January 2020. CH4 undersaturation (mean: 82 ± 20%) was found throughout the water column. In subsurface waters, the distribution of CH4 mainly resulted from mixing of water masses and in situ consumption, whereas the CH4 concentrations in the surface mixed layer were mainly driven by air–sea exchange and diapycnal diffusion between the surface and subsurface layers, as well as consumption of CH4. With a mean air–sea CH4 flux density of −0.44 ± 0.34 μmol m−2 d−1, the Ross Sea was a substantial sink for atmospheric CH4 during austral summer, which is in contrast with most oceanic regions, which are known sources.

Description: Selection along environmental gradients can drive reproductive isolation and speciation. Among fishes, salinity is a major factor limiting species distributions, and despite its importance in generating species diversity, speciation events between marine and freshwater are rare. Here, we tested for mechanisms of reproductive isolation between locally adapted freshwater and brackish water-native populations of killifish, Fundulus heteroclitus, from either side of a hybrid zone along a salinity gradient. There was evidence for pre-zygotic endogenous reproductive isolation with reduced fertilization success between crosses of freshwater-native males and brackish water-native females. Exogenous pre-zygotic isolation was also present where females had highest fertilization in their native salinity. We used a replicated mass spawning design to test for mate choice in both brackish and fresh water. After genotyping 187 parents and 2523 offspring at 2347 SNPs across the genome, 85% of offspring were successfully assign to their parents. However, no reinforcing mate choice was observed. These results therefore demonstrate emerging, yet limited, reproductive isolation and incipient speciation across a marine to freshwater salinity gradient and suggest that both endogenous and exogenous mechanisms, but not assortative mating, contribute to divergence.

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