ALBERT

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: The climactic Los Chocoyos (LCY) eruption from Atitlán caldera (Guatemala) is a key chronostratigraphic marker for the Quaternary period given the extensive distribution of its deposits that reached both the Pacific and Atlantic Oceans. Despite LCY tephra being an important marker horizon, a radioisotopic age for this eruption has remained elusive. Using zircon (U–Th)/He geochronology, we present the first radioisotopically determined eruption age for the LCY of 75 ± 2 ka. Additionally, the youngest zircon crystallization 238U–230Th rim ages in their respective samples constrain eruption age maxima for two other tephra units that erupted from Atitlán caldera, W-Fall (130 +16/−14 ka) and I-Fall eruptions (56 +8.2/−7.7 ka), which under- and overlie LCY tephra, respectively. Moreover, rim and interior zircon dating and glass chemistry suggest that before eruption silicic magma was stored for 〉80 kyr, with magma accumulation peaking within ca. 35 kyr before the LCY eruption during which the system may have developed into a vertically zoned magma chamber. Based on an updated distribution of LCY pyroclastic deposits, a new conservatively estimated volume of ~1220 ± 150 km3 is obtained (volcanic explosivity index VEI 〉 8), which confirms the LCY eruption as the first-ever recognized supereruption in Central America.

Description: A new autotrophic hydrogen‐oxidizing Chromatiaceae bacterium, namely bacterium CTD079, was enriched from a water column sample at 1500 m water depth in the southern Pacific Ocean. Based on the phylogeny of 16S rRNA genes it was closely related to a scaly snail endosymbiont (99.2% DNA sequence identity) whose host so far is only known to colonize hydrothermal vents along the Indian ridge. The average nucleotide identity between the genomes of CTD079 and the snail endosymbiont was 91%. The observed differences likely reflect adaptations to their specific habitats. For example, CTD079 encodes additional enzymes like the formate dehydrogenase increasing the organism's spectrum of energy generation pathways. Other additional physiological features of CTD079 included the increase of viral defense strategies, secretion systems and specific transporters for essential elements. These important genome characteristics suggest an adaptation to life in the open ocean.

Description: Idealized models or emulators of volcanic aerosol forcing have been widely used to reconstruct the spatiotemporal evolution of past volcanic forcing. However, existing models, including the most recently developed Easy Volcanic Aerosol (EVA; Toohey et al., doi: 10.5194/gmd‐2016‐83), (i) do not account for the height of injection of volcanic SO urn:x-wiley:jgrd:media:jgrd55987:jgrd55987-math-0001; (ii) prescribe a vertical structure for the forcing; and (iii) are often calibrated against a single eruption. We present a new idealized model, EVA_H, that addresses these limitations. Compared to EVA, EVA_H makes predictions of the global mean stratospheric aerosol optical depth that are (i) similar for the 1979–1998 period characterized by the large and high‐altitude tropical SO urn:x-wiley:jgrd:media:jgrd55987:jgrd55987-math-0002 injections of El Chichón (1982) and Mount Pinatubo (1991); (ii) significantly improved for the 1998–2015 period characterized by smaller eruptions with a large variety of injection latitudes and heights. Compared to EVA, the sensitivity of volcanic forcing to injection latitude and height in EVA_H is much more consistent with results from climate models that include interactive aerosol chemistry and microphysics, even though EVA_H remains less sensitive to eruption latitude than the latter models. We apply EVA_H to investigate potential biases and uncertainties in EVA‐based volcanic forcing data sets from phase 6 of the Coupled Model Intercomparison Project (CMIP6). EVA and EVA_H forcing reconstructions do not significantly differ for tropical high‐altitude volcanic injections. However, for high‐latitude or low‐altitude injections, our reconstructed forcing is significantly lower. This suggests that volcanic forcing in CMIP6 last millenium experiments may be overestimated for such eruptions.

Description: Three volcanic arcs have been the source of New Zealand's volcanic activity since the Neogene: Northland arc, Coromandel Volcanic Zone (CVZ) and Taupō Volcanic Zone (TVZ). The eruption chronology for the Quaternary, sourced by the TVZ, is well studied and established, whereas the volcanic evolution of the precursor arc systems, like the CVZ (central activity c. 18 to 2 Ma), is poorly known due to limited accessibility to, or identification of, onshore volcanic deposits and their sources. Here, we investigate the marine tephra record of the Neogene, mostly sourced by the CVZ, of cores from IODP Exp. 375 (Sites U1520 and U1526), ODP Leg 181 (Sites 1123, 1124 and 1125), IODP Leg 329 (Site U1371) and DSDP Leg 90 (Site 594) offshore of New Zealand. In total, we identify 306 primary tephra layers in the marine sediments. Multi-approach age models (e.g. biostratigraphy, zircon ages) are used in combination with geochemical fingerprinting (major and trace element compositions) and the stratigraphic context of each marine tephra layer to establish 168 tie-lines between marine tephra layers from different holes and sites. Following this approach, we identify 208 explosive volcanic events in the Neogene between c. 17.5 and 2.6 Ma. This is the first comprehensive study of New Zealand's Neogene explosive volcanism established from tephrochronostratigraphic studies, which reveals continuous volcanic activity between c. 12 and 2.6 Ma with an abrupt compositional change at c. 4.5 Ma, potentially associated with the transition from CVZ to TVZ. Key Points New Zealand's Neogene explosive volcanism based on the marine tephra record Geochemical fingerprinting of marine tephra layers across the study area to establish volcanic events Insights into geochemical variations with time, repose times and spatiotemporal distribution

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: 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: 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: A site at the gas hydrate stability limit was investigated offshore northwestern Svalbard to study methane transport in sediment. The site was characterized by chemosynthetic communities (sulfur bacteria mats, tubeworms) and gas venting. Sediments were sampled with in‐situ porewater collectors and by gravity coring followed by analyses of porewater constituents, sediment and carbonate geochemistry, and microbial activity, taxonomy, and lipid biomarkers. Sulfide and alkalinity concentrations showed concentration maxima in near‐surface sediments at the bacterial mat and deeper maxima at the gas vent site. Sediments at the periphery of the chemosynthetic field were characterized by two sulfate‐methane transition zones (SMTZ) at ~204 and 45 cm depth, where activity maxima of microbial anaerobic oxidation of methane (AOM) with sulfate were found. Amplicon sequencing and lipid biomarker indicate that AOM at the SMTZs was mediated by ANME‐1 archaea. A 1D numerical transport reaction model suggests that the deeper SMTZ‐1 formed on centennial scale by vertical advection of methane, while the shallower SMTZ‐2 could only be reproduced by non‐vertical methane injections starting on decadal scale. Model results were supported by age distribution of authigenic carbonates, showing youngest carbonates within SMTZ‐2. We propose that non‐vertical methane injection was induced by increasing blockage of vertical transport or formation of sediment fractures. Our study further suggests that the methanotrophic response to the non‐vertical methane injection was commensurate with new methane supply. This finding provides new information about for the response time and efficiency of the benthic methane filter in environments with fluctuating methane transport.

Description: In recent decades, the increase in terrestrial inputs to freshwater and coastal ecosystems, especially occurring at northern latitudes, has led to a process of water color darkening known as “brownification.” To assess how brownification affects plankton community composition and functioning in northern coastal areas, an in situ mesocosm experiment using a highly colored humic substance to simulate a brownification event was performed in a North Atlantic bay (Hopavågen, Norway) in August 2019. Manual sampling for analyses of nutrient concentrations, phytoplankton pigments and zooplankton abundances was combined with high-frequency (every 15 min) monitoring of key environmental variables to investigate the response of the plankton community in terms of oxygen metabolism and community composition. In response to brownification, the oxygen gross primary production (GPP) and community respiration (R) slowed down significantly, by almost one-third. However, GPP and R both decreased to the same extent; thus, the oxygen metabolic balance was not affected. Moreover, the chlorophyll-a concentration significantly decreased under brownification, by 9% on average, and the chemotaxonomic pigment composition of the phytoplankton changed, indicating their acclimation to the reduced light availability. In addition, brownification seemed to favor appendicularians, the dominant mesozooplankton group in the mesocosms, which potentially contributed to lowering the phytoplankton biomass. In conclusion, the results of this in situ mesocosm experiment suggest that brownification could induce significant changes in phytoplankton and zooplankton community composition and significantly alter the overall oxygen metabolism of plankton communities in a northern Atlantic bay.

Description: Submarine groundwater discharge (SGD) into coastal areas is a common global phenomenon and is rapidly gaining scientific interest due to its influence on marine ecology, the coastal sedimentary environment and its potential as a future freshwater resource. We conducted an integrated study of hydroacoustic surveys combined with geochemical porewater and water column investigations at a well‐known groundwater seep site in Eckernförde Bay (Germany). We aim to better constrain the effects of shallow gas and SGD on high frequency multibeam backscatter data and to present acoustic indications for submarine groundwater discharge. Our high‐quality hydroacoustic data reveal hitherto unknown internal structures within the pockmarks in Eckernförde Bay. Using precisely positioned sediment core samples, our hydroacoustic‐geochemical approach can differentiate intra‐pockmark regimes that were formerly assigned to pockmarks of a different nature. We demonstrate that high‐frequency multibeam data, in particular the backscatter signals, can be used to detect shallow free gas in areas of enhanced groundwater advection in muddy sediments. Intriguingly, our data reveal relatively small (typically 〈15 m across) pockmarks within the much larger, previously mapped, pockmarks. The small pockmarks, which we refer to as “intra‐pockmarks”, have formed due to the localized ascent of gas and groundwater; they manifest themselves as a new type of ‘eyed’ pockmarks, revealed by their acoustic backscatter pattern. Our data suggest that, in organic‐rich muddy sediments, morphological lows combined with a strong multibeam backscatter signal can be indicative of free shallow gas and subsequent advective groundwater flow.

Description: Males and females follow distinct life-history strategies that have co-evolved with several sex-specific traits. Higher investment into parental investment (PI) demands an increased lifespan. Thus, resource allocation toward an efficient immune system is mandatory. In contrast, resources allocated toward secondary sexual signals (ornamentation) may negatively correlate with investment into immunity and ultimately result in a shorter lifespan. Previous studies have addressed how resource allocation toward single sex-specific traits impacts lifetime reproductive success (LRS). However, the trade-offs between diverse sex-specific characteristics and their impact on LRS remain largely unassessed impeding our understanding of life-history evolution. We have designed a theoretical framework (informed by experimental data and evolutionary genetics) that explores the effects of multiple sex-specific traits and assessed how they influence LRS. From the individual sex-specific traits, we inferred the consequences at the population level by evaluating adult sex ratios (ASR). Our theory implies that sex-specific resource allocation toward the assessed traits resulted in a biased ASR. Our model focuses on the impact of PI, ornamentation, and immunity as causal to biased ASR. The framework developed herein can be employed to understand the combined impact of diverse sex-specific traits on the LRS and the eventual population dynamics of particular model systems.

Description: The Los Chocoyos (14.6°N, 91.2°W) supereruption happened ∼75,000 years ago in Guatemala and was one of the largest eruptions of the past 100,000 years. It emitted enormous amounts of sulfur, chlorine, and bromine, with multi‐decadal consequences for the global climate and environment. Here, we simulate the impact of a Los Chocoyos‐like eruption on the quasi‐biennial oscillation (QBO), an oscillation of zonal winds in the tropical stratosphere, with a comprehensive aerosol chemistry Earth System Model. We find a ∼10‐year disruption of the QBO starting 4 months post eruption, with anomalous easterly winds lasting ∼5 years, followed by westerlies, before returning to QBO conditions with a slightly prolonged periodicity. Volcanic aerosol heating and ozone depletion cooling leads to the QBO disruption and anomalous wind regimes through radiative changes and wave‐mean flow interactions. Different model ensembles, volcanic forcing scenarios and results of a second model back up the robustness of our results.

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: Dinoflagellates are a diverse group of eukaryotic microbes that are ubiquitous in aquatic environments. Largely photosynthetic, they encompass symbiotic, parasitic, and free-living lineages with a broad spectrum of trophism. Many free-living taxa can produce bioactive secondary metabolites such as biotox- ins, some of which cause harmful algal blooms. In contrast, most symbiotic species are crucial for sustaining coral reef health. The year 2023 marked a decade since the first genome data of dinoflagellates became available. The growing genome-scale resources for these taxa are highlighting their remarkable evolutionary and genomic complexities. Here, we discuss the prospect of developing dinoflagellate models using the criteria of accessibil- ity, tractability, resources, research support, and promise. Moving forward in the post-genomic era, we argue for the development of fit-to-purpose models that tailor to specific biological contexts, and that a one-size-fits-all model is inadequate for encapsulating the complex biology, ecology, and evolutionary history of dinoflagellates.

Description: Mesoscale eddies are frequently observed in the Eastern Tropical North Atlantic (ETNA), yet their effects on the transport and distribution of biogeochemical solutes, and specifically on the production and remineralization of dissolved organic matter (DOM) remain difficult to elucidate. Here, we investigated the submesoscale variability of chromophoric DOM (CDOM) and fluorescent DOM (FDOM) together with microbial production and remineralization processes in two cyclonic eddies (CEs) in the ETNA during summer and winter 2019. One CE, formed near the coast off Mauritania during the post-upwelling season, was sampled along a ∼900 km zonal corridor between Mauritania and the Cape Verde Islands. The other CE, formed nearby Brava Island, was out of coastal influence. Four fluorescent components were identified with parallel factor analysis, two humic-like, and two protein-like components. Humic-like FDOM components correlated to optode-based community respiration and were also good indicators of upwelling associated with the Brava Island CE as they correlated to physical parameters (e.g., temperature) and to dissolved inorganic nitrogen. The tryptophan-like FDOM components correlated with the carbon and nitrogen content of semi-labile DOM, phytoplankton biomass, community respiration, and bacterial production. Overall, our study revealed that DOM optical properties are suitable for tracing freshly produced organic matter and the transport of remineralized DOM within offshore eddies. Key Points: - Four fluorescent dissolved organic matter (FDOM) components were studied in two cyclonic eddies (CEs) in the Eastern Tropical North Atlantic - Tryptophan-like FDOM was an indicator of the CEs' productivity as it correlated with semi-labile dissolved organic matter and microbial metabolic activities - Humic-like FDOM was a by-product of microbial respiration; its distribution within an offshore CE covaried with nutrient upwelling

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: Whereas fungal symbionts of terrestrial plants are among the most widespread and well-studied symbioses, relatively little is known about fungi that are associated with macroalgae. To fill the gap in marine fungal taxonomy, we combined simple culture methods with amplicon sequencing to characterize the fungal communities associated with three brown (Sargassum muticum, Pelvetia canaliculata, Himanthalia elongata) and two red (Mastocarpus stellatus, Chondrus crispus) macroalgae from one intertidal zone. In addition to characterizing novel fungal diversity, we tested three hypotheses: fungal diversity and community composition vary (i) among species distributed at different tidal heights, (ii) among tissue types (apices, mid-thallus, and stipe), and (iii) among ‘isomorphic’ C. crispus life cycle stages. Almost 70% of our reads were classified as Ascomycota, 29% as Basidiomycota, and 1% that could not be classified to phylum. Thirty fungal isolates were obtained, 18 of which were also detected with amplicon sequencing. Fungal communities differed by host and tissue type. Interestingly, P. canaliculata, a fucoid at the extreme high intertidal, did not show differences in fungal diversity across the thallus. As found in filamentous algal endophytes, fungal diversity varied among the three life cycle stages in C. crispus. Female gametophytes were also compositionally more dispersed as compared to the less variable tetrasporophytes and male gametophytes. We demonstrate the power of combining relatively simple cultivation and sequencing approaches to characterize and study macroalgal-fungal associations and highlight the need to understand the role of fungi in near-shore marine ecosystems.

Description: Communities are shaped by scale dependent processes. To study the diversity and variation of microbial communities across scales, the invasive and widespread seaweed Agarophyton vermiculophyllum presents a unique opportunity. We characterized pro‐ and eukaryotic communities associated with this holobiont across its known distribution range, which stretches over the northern hemisphere. Our data reveal that community composition and diversity in the holobiont vary at local but also larger geographic scales. While processes acting at the local scale (i.e., within population) are the main structuring drivers of associated microbial communities, changes in community composition also depend on processes acting at larger geographic scales. Interestingly, the largest analysed scale (i.e., native and non‐native ranges) explained variation in the prevalence of predicted functional groups, which could suggest a functional shift in microbiota occurred over the course of the invasion process. While high variability in microbiota at the local scale supports A. vermiculophyllum to be a generalist host, we also identified a number of core taxa. These geographically independent holobiont members imply that cointroduction of specific microbiota may have additionally promoted the invasion process.

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 unique male pregnancy in pipefishes and seahorses ranges from basic attachment (pouch-less species: Nerophinae) of maternal eggs to specialized internal gestation in pouched species (e.g. Syngnathus and Hippocampus) with many transitions in between. Due to this diversity, male pregnancy offers a unique platform for assessing physiological and molecular adaptations in pregnancy evolution. These insights will contribute to answering long-standing questions of why and how pregnancy evolved convergently in so many vertebrate systems. To understand the molecular congruencies and disparities in male pregnancy evolution, we compared transcriptome-wide differentially expressed genes in four syngnathid species, at four pregnancy stages (nonpregnant, early, late and parturition). Across all species and pregnancy forms, metabolic processes and immune dynamics defined pregnancy stages, especially pouched species shared expression features akin to female pregnancy. The observed downregulation of adaptive immune genes in early-stage pregnancy and its reversed upregulation during late/parturition in pouched species, most notably in Hippocampus, combined with directionless expression in the pouch-less species, suggests immune modulation to be restricted to pouched species that evolved placenta-like systems. We propose that increased foeto-paternal intimacy in pouched syngnathids commands immune suppression processes in early gestation, and that the elevated immune response during parturition coincides with pouch opening and reduced progeny reliance. Immune response regulation in pouched species supports the recently described functional MHC II pathway loss as critical in male pregnancy evolution. The independent co-option of similar genes and pathways both in male and female pregnancy highlights immune modulation as crucial for the evolutionary establishment of pregnancy.

Description: Forecasting and early warning systems are important investments to protect lives, properties and livelihood. While early warning systems are frequently used to predict the magnitude, location and timing of potentially damaging events, these systems rarely provide impact estimates, such as the expected amount and distribution of physical damage, human consequences, disruption of services or financial loss. Complementing early warning systems with impact forecasts has a two‐fold advantage: it would provide decision makers with richer information to take informed decisions about emergency measures, and focus the attention of different disciplines on a common target. This would allow capitalizing on synergies between different disciplines and boosting the development of multi‐hazard early warning systems. This review discusses the state‐of‐the‐art in impact forecasting for a wide range of natural hazards. We outline the added value of impact‐based warnings compared to hazard forecasting for the emergency phase, indicate challenges and pitfalls, and synthesize the review results across hazard types most relevant for Europe. Plain language summary Forecasting and early warning systems are important investments to protect lives, properties and livelihood. While such systems are frequently used to predict the magnitude, location and timing of potentially damaging events, they rarely provide impact estimates, such as the expected physical damage, human consequences, disruption of services or financial loss. Extending hazard forecast systems to include impact estimates promises many benefits for the emergency phase, for instance, for organising evacuations. We review and compare the state‐of‐the‐art of impact forcasting across a wide range of natural hazards, and outline opportunities and key challenges for research and development of impact forecasting.

Description: Predation is a pervasive force that structures food webs and directly influences ecosystem functioning. The relative body sizes of predators and prey may be an important determinant of interaction strengths. However, studies quantifying the combined influence of intra‐ and interspecific variation in predator–prey body size ratios are lacking. We use a comparative functional response approach to examine interaction strengths between three size classes of invasive bluegill and largemouth bass toward three scaled size classes of their tilapia prey. We then quantify the influence of intra‐ and interspecific predator–prey body mass ratios on the scaling of attack rates and handling times. Type II functional responses were displayed by both predators across all predator and prey size classes. Largemouth bass consumed more than bluegill at small and intermediate predator size classes, while large predators of both species were more similar. Small prey were most vulnerable overall; however, differential attack rates among prey were emergent across predator sizes. For both bluegill and largemouth bass, small predators exhibited higher attack rates toward small and intermediate prey sizes, while larger predators exhibited greater attack rates toward large prey. Conversely, handling times increased with prey size, with small bluegill exhibiting particularly low feeding rates toward medium–large prey types. Attack rates for both predators peaked unimodally at intermediate predator–prey body mass ratios, while handling times generally shortened across increasing body mass ratios. We thus demonstrate effects of body size ratios on predator–prey interaction strengths between key fish species, with attack rates and handling times dependent on the relative sizes of predator–prey participants. Considerations for intra‐ and interspecific body size ratio effects are critical for predicting the strengths of interactions within ecosystems and may drive differential ecological impacts among invasive species as size ratios shift.

Description: Along ultraslow spreading ridges melt is distributed unequally, but melt focusing guides melt away from amagmatic segments toward volcanic centers. An interplay of tectonism and magmatism is thought to control melt ascent, but the detailed process of melt extraction is not yet understood. We present a detailed image of the seismic velocity structure of the Logachev volcanic center and adjacent region along the Knipovich Ridge. With travel times of P- and S-waves of 3,959 earthquakes we performed a local earthquake tomography. We simultaneously inverted for source locations, velocity structure and the Vp/Vs-ratio. An extensive low velocity anomaly coincident with high Vp/Vs-ratios 〉1.9 lies underneath the volcanic center at depths of 10 km below sea level in an aseismic area. More shallow, tightly clustered earthquake swarms connect the anomaly to a shallow anomaly with high Vp/Vs-ratio beneath the basaltic seafloor. We consider the deep low-velocity anomaly to represent an area of partial melt from which melts ascent vertically to the surface and northwards into the adjacent segment. By comparing tomographic studies of the Logachev and Southwest Indian Ridge Segment-8 volcano we conclude that volcanic centers of ultraslow spreading ridges host spatially confined, circular partial melt areas below 10 km depth, in contrast to the shallow extended melt lenses along fast spreading ridges. Lateral feeding over distances of 35 km is possible at orthogonal spreading segments, but limited at the obliquely spreading Knipovich Ridge. Key Points - Active volcanic centers at ultraslow spreading ridges host deeper and more confined partial melt areas than faster spreading ridges - Earthquake swarms delineate melt ascent paths from the partial melt area to the surface - Lateral feeding at shallow depths into subordinate segments is prevented by ridge obliquity

Description: Spreading centers in the proximity of back‐rolling subduction zones constitute an ideal natural laboratory to investigate the interaction of magmatism and tectonism during the early evolution of back‐arc basins. Using 32 days of ocean bottom seismometer data, we located 697 microeathquakes at the southern Fonualei Rift and Spreading Center (S‐FRSC). The majority of epicenters concentrate along the central region of the axial valley, marking the active ridge axis. Only odd events were associated with the prominent faults bounding the axial valley. About 450 events are spatially clustered around 17°42’S and their waveforms show a pronounced similarity. Most of these events are associated with a 138 hours lasting earthquake swarm. The tectonic structure of the ridge axis in the S‐FRSC resembles a series of left‐stepping en echelon segments, expressed at the seafloor by numerous volcanic ridges. The recorded earthquake swarm is located at the stepover of two en echelon segments suggesting that the earthquake swarm is mainly tectonically driven. The events directly beneath our seismic network indicate a maximum depth of brittle faulting down to about 14 km below the seafloor. This is within the maximum depth range of brittle faulting at ultraslow mid‐ocean ridges. Since the thickness of the brittle lithosphere is mainly controlled by temperature, our results suggest a sub‐axial thermal structure similar to that of ultraslow mid‐ocean ridges of similar opening rates.

Description: Seaweeds are colonized by a microbial community, which can be directly linked to their performance. This community is shaped by an interplay of stochastic and deterministic processes, including mechanisms which the holobiont host deploys to manipulate its associated microbiota. The Anna Karenina principle predicts that when a holobiont is exposed to suboptimal or stressful conditions, these host mechanisms may be compromised. This leads to a relative increase of stochastic processes that may potentially result in the succession of a microbial community harmful to the host. Based on this principle, we used the variability in microbial communities (i.e., beta diversity) as a proxy for stability within the invasive holobiont Gracilaria vermiculophylla during a simulated invasion in a common garden experiment. Independent of host range, host performance declined at elevated temperature (22°C) and disease incidence and beta diversity increased. Under thermally stressful conditions, beta diversity increased more in epibiota from native populations, suggesting that epibiota from non-native holobionts are thermally more stable. This pattern reflects an increase in deterministic processes acting on epibiota associated with non-native hosts, which in the setting of a common garden can be assumed to originate from the host itself. Therefore, these experimental data suggest that the invasion process may have selected for hosts better able to maintain stable microbiota during stress. Future studies are needed to identify the underlying host mechanisms.

Description: Carbon cycling by Antarctic microbial plankton is poorly understood but it plays a major role in CO2 sequestration in the Southern Ocean. We investigated the summer bacterioplankton community in the largely understudied Weddell Sea, applying Illumina amplicon sequencing, measurements of bacterial production and chemical analyses of organic matter. The results revealed that the patchy distribution of productive coastal polynyas and less productive, mostly ice-covered sites was the major driver of the spatial changes in the taxonomic composition and activity of bacterioplankton. Gradients in organic matter availability induced by phytoplankton blooms were reflected in the concentrations and composition of dissolved carbohydrates and proteins. Bacterial production at bloom stations was, on average, 2.7 times higher than at less productive sites. Abundant bloom-responsive lineages were predominately affiliated with ubiquitous marine taxa, including Polaribacter, Yoonia-Loktanella, Sulfitobacter, the SAR92 clade, and Ulvibacter, suggesting a widespread genetic potential for adaptation to sub-zero seawater temperatures. A co-occurrence network analysis showed that dominant taxa at stations with low phytoplankton productivity were highly connected, indicating beneficial interactions. Overall, our study demonstrates that heterotrophic bacterial communities along Weddell Sea ice shelves were primarily constrained by the availability of labile organic matter rather than low seawater temperature.

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: 1. Geolocating aquatic animals with acoustic tags has been ongoing for decades, relying on the detection of acoustic signals at multiple receivers with known positions to calculate a 2D or 3D position, and ultimately recreate the path of an aquatic animal from detections at fixed stations. 2. This method of underwater geolocation is evolving with new software and hardware options available to help investigators design studies and calculate positions using solvers based predominantly on time-difference-of-arrival and time-of-arrival. 3. We provide an overview of the considerations necessary to implement positioning in aquatic acoustic telemetry studies, including how to design arrays of receivers, test performance, synchronize receiver clocks and calculate positions from the detection data. We additionally present some common positioning algorithms, including both the free open-source solvers and the ‘black-box’ methods provided by some manufacturers for calculating positions. 4. This paper is the first to provide a comprehensive overview of methods and considerations for designing and implementing better positioning studies that will support users, and encourage further knowledge advances in aquatic systems.

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: What controls subduction megathrust seismogenesis downdip of the mantle wedge corner (MWC)? We propose that, in the region of the 2010 Mw=8.8 Maule, Chile, earthquake, serpentine minerals derived from the base of the hydrated mantle wedge exert a dominant control. Based on modeling, we predict that the megathrust fault zone near the MWC contains abundant lizardite/chrysotile‐rich serpentinite that transforms to antigorite‐rich serpentinite at greater depths. From the MWC at 32–40 km depth to at least 55 km, the predominantly velocity‐strengthening megathrust accommodated dynamic propagation of the 2010 rupture but with small slip and negative stress drop. The downdip distribution of interplate aftershocks exhibits a gap around the MWC that can be explained by the velocity‐strengthening behavior of lizardite/chrysotile. Interspersed velocity‐weakening and dynamic weakening antigorite‐rich patches farther downdip may be responsible for increased abundance of aftershocks and possibly for some of the high‐frequency energy radiation during the 2010 rupture.

Description: The green seaweed Ulva is a model system to study seaweed–bacteria interactions, but the impact of environmental drivers on the dynamics of these interactions is little understood. In this study, we investigated the stability and variability of the seaweed-associated bacteria across the Atlantic–Baltic Sea salinity gradient. We characterized the bacterial communities of 15 Ulva sensu lato species along 2,000 km of coastline in a total of 481 samples. Our results demonstrate that the Ulva-associated bacterial composition was strongly structured by both salinity and host species (together explaining between 34% and 91% of the variation in the abundance of the different bacterial genera). The largest shift in the bacterial consortia coincided with the horohalinicum (5–8 PSU, known as the transition zone from freshwater to marine conditions). Low-salinity communities especially contained high relative abundances of Luteolibacter, Cyanobium, Pirellula, Lacihabitans and an uncultured Spirosomaceae, whereas high-salinity communities were predominantly enriched in Litorimonas, Leucothrix, Sulfurovum, Algibacter and Dokdonia. We identified a small taxonomic core community (consisting of Paracoccus, Sulfitobacter and an uncultured Rhodobacteraceae), which together contributed to 14% of the reads per sample, on average. Additional core taxa followed a gradient model, as more core taxa were shared between neighbouring salinity ranges than between ranges at opposite ends of the Atlantic–Baltic Sea gradient. Our results contradict earlier statements that Ulva-associated bacterial communities are taxonomically highly variable across individuals and largely stochastically defined. Characteristic bacterial communities associated with distinct salinity regions may therefore facilitate the host's adaptation across the environmental gradient.

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: 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: 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: Key Points: • Surface seawater dimethylsulfide (DMS) concentrations remain unchanged after sea ice retreat in the western Canada Basin • Increased wind speed is a critical factor driving the enhancement of dimethylsulfide (DMS) flux after sea ice retreat at high latitudes in the Arctic Ocean • Nutrient supply is hypothesized to significantly impact dimethylsulfide (DMS) distribution in the western Arctic Ocean Abstract: The receding of the seasonal ice cover in the Arctic due to climate change has been predicted by models to increase climate-active biogenic trace gas emissions, specifically those of dimethylsulfide (DMS). However, insufficient DMS measurements are currently available to either support or refute this hypothesis and to fully understand the various responses of oceanic DMS in a rapidly changing Arctic Ocean environment. Here, we present high-resolution surface water DMS data collected in the summer of 2014 in combination with a suite of ancillary variables including sea ice cover, salinity, and nutrients. We show that surface seawater DMS concentrations, generally below 0.5 nmol L−1, remained unchanged in the Canada Basin after sea ice retreat probably due to insufficient nutrients supply to the upper mixed layer and resulting low primary production. Moreover, in the Chukchi shelf region, DMS concentrations decreased following a phytoplankton bloom due to the rapid depletion and slow resupply of nutrients. Although the DMS sea-to-air fluxes were not high from a global perspective, they increased by a factor of 4-fold after sea ice retreat in the Arctic Ocean high latitudes. This increase in DMS flux was mainly driven by increased wind speed. This work provides unique observations and insights on how surface seawater DMS and flux to the atmosphere may change in the future Arctic Ocean.

Description: The Tierra Blanca (TB) eruptive suite comprises the last four major eruptions of Ilopango caldera in El Salvador (≤45 ka), including the youngest Tierra Blanca Joven eruption (TBJ; ∼106 km3): the most voluminous event during the Holocene in Central America. Despite the protracted and productive history of explosive silicic eruptions at Ilopango caldera, many aspects regarding the longevity and the prevailing physicochemical conditions of the underlying magmatic system remain unknown. Zircon 238U-230Th geochronology of the TB suite (TBJ, TB2, TB3, and TB4) reveals a continuous and overlapping crystallization history among individual eruptions, suggesting persistent melt presence in thermally and compositionally distinct magma reservoirs over the last ca. 80 kyr. The longevity of zircon is in contrast to previously determined crystallization timescales of 〈10 kyr for major mineral phases in TBJ. This dichotomy is explained by a process of rhyolitic melt segregation from a crystal-rich refractory residue that incorporates zircon, whereas a new generation of major mineral phases crystallized shortly before eruption. Ti-in-zircon temperatures and amphibole geothermobarometry suggest that rhyolitic melt was extracted from different storage zones of the magma reservoir as indicated by distinct but synchronous thermochemical zircon histories among the TB suite eruptions. Zircon from TBJ and TB2 suggests magma differentiation within deeper and hotter parts of the reservoir, whereas zircon from TB3 and TB4 instead hints at crystallization in comparatively shallower and cooler domains. The assembly of the voluminous TBJ magma reservoir was also likely enhanced by cannibalization of hydrothermally altered components as suggested by low-δ18O values in zircon (+4.5 ± 0.3‰).