ALBERT

Description: Glacial isostatic adjustment is largely governed by the rheological properties of the Earth's mantle. Large mass redistributions in the ocean–cryosphere system and the subsequent response of the viscoelastic Earth have led to dramatic sea level changes in the past. This process is ongoing, and in order to understand and predict current and future sea level changes, the knowledge of mantle properties such as viscosity is essential. In this study, we present a method to obtain estimates of mantle viscosities by the assimilation of relative sea level rates of change into a viscoelastic model of the lithosphere and mantle. We set up a particle filter with probabilistic resampling. In an identical twin experiment, we show that mantle viscosities can be recovered in a glacial isostatic adjustment model of a simple three-layer Earth structure consisting of an elastic lithosphere and two mantle layers of different viscosity. We investigate the ensemble behaviour on different parameters in the following three set-ups: (1) global observations data set since last glacial maximum with different ensemble initialisations and observation uncertainties, (2) regional observations from Fennoscandia or Laurentide/Greenland only, and (3) limiting the observation period to 10 ka until the present. We show that the recovery is successful in all cases if the target parameter values are properly sampled by the initial ensemble probability distribution. This even includes cases in which the target viscosity values are located far in the tail of the initial ensemble probability distribution. Experiments show that the method is successful if enough near-field observations are available. This makes it work best for a period after substantial deglaciation until the present when the number of sea level indicators is relatively high.

Description: Accurate and reliable biodiversity estimates of marine zooplankton are a prerequisite to understand how changes in diversity can affect whole ecosystems. Species identification in the deep sea is significantly impeded by high numbers of new species and decreasing numbers of taxonomic experts, hampering any assessment of biodiversity. We used in parallel morphological, genetic, and proteomic characteristics of specimens of calanoid copepods from the abyssal South Atlantic to test if proteomic fingerprinting can accelerate estimating biodiversity. We cross-validated the respective molecular discrimination methods with morphological identifications to establish COI and proteomic reference libraries, as they are a pre-requisite to assign taxonomic information to the identified molecular species clusters. Due to the high number of new species only 37% of the individuals could be assigned to species or genus level morphologically. COI sequencing was successful for 70% of the specimens analysed, while proteomic fingerprinting was successful for all specimens examined. Predicted species richness based on morphological and molecular methods was 42 morphospecies, 56 molecular operational taxonomic units (MOTUs) and 79 proteomic operational taxonomic units (POTUs), respectively. Species diversity was predicted based on proteomic profiles using hierarchical cluster analysis followed by application of the variance ratio criterion for identification of species clusters. It was comparable to species diversity calculated based on COI sequence distances. Less than 7% of specimens were misidentified by proteomic profiles when compared with COI derived MOTUs, indicating that unsupervised machine learning using solely proteomic data could be used for quickly assessing species diversity.

Description: Based on the numerical weather prediction model COSMO of Germany's national meteorological service (Deutscher Wetterdienst, DWD), regional reanalysis datasets have been developed with grid spacing of up to 2 km. This development started as a fundamental research activity within the Hans-Ertel-Centre for Weather Research (HErZ) at the University of Bonn and the University of Cologne. Today, COSMO reanalyses are an established product of the DWD and have been widely used in applications on European and national German level. Successful applications of COSMO reanalyses include renewable energy assessments as well as meteorological risk estimates. The COSMO reanalysis datasets are now publicly available and provide spatio-temporal consistent data of atmospheric parameters covering both near-surface conditions and vertical profiles. This article reviews the status of the COSMO reanalyses, including evaluation results and applications. In many studies, evaluation of the COSMO reanalyses point to an overall good quality and often an added value compared to different contemporary global reanalysis datasets. We further outline current plans for the further development and application of regional reanalyses in the HErZ research group Cologne/Bonn in collaboration with the DWD.

Description: The effective radiative forcing, which includes the instantaneous forcing plus adjustments from the atmosphere and surface, has emerged as the key metric of evaluating human and natural influence on the climate. We evaluate effective radiative forcing and adjustments in 17 contemporary climate models that are participating in the Coupled Model Intercomparison Project (CMIP6) and have contributed to the Radiative Forcing Model Intercomparison Project (RFMIP). Present-day (2014) global-mean anthropogenic forcing relative to pre-industrial (1850) levels from climate models stands at 2.00 (±0.23) W m−2, comprised of 1.81 (±0.09) W m−2 from CO2, 1.08 (± 0.21) W m−2 from other well-mixed greenhouse gases, −1.01 (± 0.23) W m−2 from aerosols and −0.09 (±0.13) W m−2 from land use change. Quoted uncertainties are 1 standard deviation across model best estimates, and 90 % confidence in the reported forcings, due to internal variability, is typically within 0.1 W m−2. The majority of the remaining 0.21 W m−2 is likely to be from ozone. In most cases, the largest contributors to the spread in effective radiative forcing (ERF) is from the instantaneous radiative forcing (IRF) and from cloud responses, particularly aerosol–cloud interactions to aerosol forcing. As determined in previous studies, cancellation of tropospheric and surface adjustments means that the stratospherically adjusted radiative forcing is approximately equal to ERF for greenhouse gas forcing but not for aerosols, and consequentially, not for the anthropogenic total. The spread of aerosol forcing ranges from −0.63 to −1.37 W m−2, exhibiting a less negative mean and narrower range compared to 10 CMIP5 models. The spread in 4×CO2 forcing has also narrowed in CMIP6 compared to 13 CMIP5 models. Aerosol forcing is uncorrelated with climate sensitivity. Therefore, there is no evidence to suggest that the increasing spread in climate sensitivity in CMIP6 models, particularly related to high-sensitivity models, is a consequence of a stronger negative present-day aerosol forcing and little evidence that modelling groups are systematically tuning climate sensitivity or aerosol forcing to recreate observed historical warming.

Description: We present the result of the third Marine Ice Sheet Intercomparison project, MISMIP+. MISMIP+ is intended to be a test of ice flow models which include fast sliding marine ice streams and floating ice shelves and in particular a treatment of viscous stress that is sufficient for buttressing, where upstream ice flow is restrained by a downstream ice shelf. A set of idealized experiments test the models in circumstances where buttressing contributes to a stable steady state, and where a reduction in that buttressing causes ice stream acceleration, thinning, and grounding line retreat. We find that the most important distinction between models in this particular type of simulation is in the treatment of sliding at the bed, with other distinctions – notably the difference between the simpler and more complete treatments of englacial stress, but also the differences between numerical methods – taking a secondary role.

Description: Varved lake sediments provide long climatic records with high temporal resolution and low associated age uncertainty. Robust and detailed comparison of well-dated and annually laminated sediment records is crucial for reconstructing abrupt and regionally time-transgressive changes as well as validation of spatial and temporal trajectories of past climatic changes. The VARved sediments DAtabase (VARDA) presented here is the first data compilation for varve chronologies and associated palaeoclimatic proxy records. The current version 1.0 allows detailed comparison of published varve records from 95 lakes. VARDA is freely accessible and was created to assess outputs from climate models with high-resolution terrestrial palaeoclimatic proxies. VARDA additionally provides a technical environment that enables to explore the database of varved lake sediments using a connected data-model and can generate a state-of-the-art graphic representation of multi-site comparison. This allows to reassess existing chronologies and tephra events to synchronize and compare even distant varved lake records. Furthermore, the present version of VARDA permits to explore varve thickness data. In this paper, we report in detail on the data mining and compilation strategies for the identification of varved lakes and assimilation of high-resolution chronologies as well as the technical infrastructure of the database. Additional paleoclimate proxy data will be provided in forthcoming updates. The VARDA graph-database and user interface can be accessed online at https://varve.gfz-potsdam.de, all datasets of version 1.0 are available at http://doi.org/10.5880/GFZ.4.3.2019.003 (Ramisch et al., 2019).

Description: Various observational estimates indicate growing mass loss at Antarctica's margins but also heavier precipitation across the continent. In the future, heavier precipitation fallen on Antarctica will counteract any stronger iceberg discharge and increased basal melting of floating ice shelves driven by a warming ocean. Here, we use from nine CMIP5 models future projections, ranging from strong mitigation efforts to business-as-usual, to run an ensemble of ice-sheet simulations. We test, how the precipitation boundary condition determines Antarctica's sea-level contribution. The spatial and temporal varying climate forcings drive ice-sheet simulations. Hence, our ensemble inherits all spatial and temporal climate patterns, which is in contrast to a spatial mean forcing. Regardless of the applied boundary condition and forcing, some areas will lose ice in the future, such as the glaciers from the West Antarctic Ice Sheet draining into the Amundsen Sea. In general the simulated ice-sheet thickness grows in a broad marginal strip, where incoming storms deliver topographically controlled precipitation. This strip shows the largest ice thickness differences between the applied precipitation boundary conditions too. On average Antarctica's ice mass shrinks for all future scenarios if the precipitation is scaled by the spatial temperature anomalies coming from the CMIP5 models. In this approach, we use the relative precipitation increment per degree warming as invariant scaling constant. In contrast, Antarctica gains mass in our simulations if we apply the simulated precipitation anomalies of the CMIP5 models directly. Here, the scaling factors show a distinct spatial pattern across Antarctica. Furthermore, the diagnosed mean scaling across all considered climate forcings is larger than the values deduced from ice cores. In general, the scaling is higher across the East Antarctic Ice Sheet, lower across the West Antarctic Ice Sheet, and lowest around the Siple Coast. The latter is located on the east side of the Ross Ice Shelf.

Description: A finely laminated lake sediment record with a basal age of 11,619 ± 603 years BP was retrieved from Lake Chatyr Kol (Kyrgyz Republic). Microfacies analysis reveals the presence of seasonal laminae (varves) from the sediment basis to ~ 360 ± 40 years BP. The Chatvd19 floating varve chronology covers the time span from 360 ± 40 years BP to the base and relies on replicate varve counts on overlapping petrographic thin sections with an uncertainty of ± 5 %. The uppermost non-varved interval was chronologically constrained by 210Pb and 137Cs γ-spectrometry and interpolation based on varve thickness measurements of adjacent varved intervals with an assumed uncertainty of 10 %. Six varve types were distinguished, are described in detail and show a changing predominance of clastic-organic, clastic-calcitic or -aragonitic, calcitic-clastic, organic-clastic and clastic-diatom varves throughout the Holocene. Variations in varve thickness and the number and composition of seasonal sublayers are attributed to 1) changes in the amount of summer or winter/spring precipitation affecting local runoff and erosion and/or to 2) evaporative conditions during summer. Radiocarbon dating of bulk organic matter, daphnia remains, aquatic plant remains and Ruppia maritima seeds reveal reservoir ages with a clear decreasing trend up core from ~ 6,150 years in the early Holocene, to ~ 3,000 years in the mid-Holocene, to ~ 1,000 years and less in the late Holocene and modern times. In contrast, two radiocarbon dates from terrestrial plant remains are in good agreement with the varve-based chronology.

Description: Plankton is a massive and phylogenetically diverse group of thousands of prokaryotes, protists (unicellular eukaryotic organisms), and metazoans (multicellular eukaryotic organisms; Fig. 1). Plankton functional diversity is at the core of various ecological processes, including productivity, carbon cycling and sequestration, nutrient cycling (Falkowski 2012), interspecies interactions, and food web dynamics and structure (D'Alelio et al. 2016). Through these functions, plankton play a critical role in the health of the coastal and open ocean and provide essential ecosystem services. Yet, at present, our understanding of plankton dynamics is insufficient to project how climate change and other human-driven impacts affect the functional diversity of plankton. That limits our ability to predict how critical ecosystem services will change in the future and develop strategies to adapt to these changes.

Description: This study investigates the mesoscale dynamics involved in the 8–11 October 2008 unseasonably strong African dust episode, during which dust was transported to the Iberian Peninsula (IP). We employ observational datasets and a high-resolution Weather Research and Forecasting model coupled with Chemistry simulations. The analysis shows that during 0900–1200 UTC 9 October, a mesoscale convective system developed over the Atlas Mountains and resulted in a southwestward propagating convective cold pool outflow on the southern foothills of the Anti-Atlas, which lifted dust from the source region. Between 1200 and 1800 UTC 9 October, new moist convection was enhanced over the Atlas Mountains due to intensifying confluence among a heat low, moist southwesterly Atlantic sea-breeze front, and northeasterly flow associated with the convective cold pool near western Algeria. This new moist convection intensified the strength of the convective cold pool outflow and haboob, both of which continued propagating southwestward. At 1200 UTC 10 October, the low-pressure system migrated poleward on the southern slopes of the Anti-Atlas Mountains in association with a mountain-plains solenoidal circulation due to the daytime differential heating between the southern slopes of the Anti-Atlas and nearby atmosphere. The deepening low-pressure and strengthening Atlantic sea-breeze redirected an equatorward advancing dust plume into the poleward direction. The dust plume ultimately crossed the Saharan Atlas Mountains on 11 October and finally impacted the IP. Key Points: - WRF-Chem simulation of an unseasonably strong haboob on the southern slopes of the Atlas Mountains - The equatorward-advancing dust plume was recirculated in the poleward direction by an Atlantic sea-breeze front - The Atlantic sea-breeze front and an intensified upper-level cutoff vortex are instrumental for dust transport over the Iberian Peninsula

Description: Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and inform on the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimated the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes and the forcings employed. This study presents results from 18 simulations from 15 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015–2100, forced with different scenarios from the Coupled Model Intercomparison Project Phase 5 (CMIP5) representative of the spread in climate model results. The contribution of the Antarctic ice sheet in response to increased warming during this period varies between −7.8 and 30.0 cm of Sea Level Equivalent (SLE). The evolution of the West Antarctic Ice Sheet varies widely among models, with an overall mass loss up to 21.0 cm SLE in response to changes in oceanic conditions. East Antarctica mass change varies between −6.5 and 16.5 cm SLE, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional mass loss of 8 mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 AOGCMs show an overall mass loss of 10 mm SLE compared to simulations done under present-day conditions, with limited mass gain in East Antarctica.

Description: An extensive, reprocessed two‐dimensional (2D) seismic data set was utilized together with available well data to study the Tiddlybanken Basin in the southeastern Norwegian Barents Sea, which is revealed to be an excellent example of base salt rift structures, evaporite accumulations and evolution of salt structures. Late Devonian–early Carboniferous NE‐SW regional extensional stress affected the study area and gave rise to three half‐grabens that are separated by a NW‐SE to NNW‐SSE trending horst and an affiliated interference transfer zone. The arcuate nature of the horst is believed to be the effect of pre‐existing Timanian basement grain, whereas the interference zone formed due to the combined effect of a Timanian (basement) lineament and the geometrical arrangement of the opposing master faults. The interference transfer zone acted as a physical barrier, controlling the facies distribution and sedimentary thickness of three‐layered evaporitic sequences (LES). During the late Triassic, the northwestern part of a salt wall was developed due to passive diapirism and its evolution was influenced by halite lithology between the three‐LES. The central and southeastern parts of the salt wall did not progress beyond the pedestal stage due to lack of halite in the deepest evaporitic sequence. During the Triassic–Jurassic transition, far‐field stresses from the Novaya Zemlya fold‐and‐thrust belt reactivated the pre‐salt Carboniferous rift structures. The reactivation led to the development of the Signalhorn Dome, rejuvenated the northwestern part of the salt wall and affected the sedimentation rates in the southeastern broad basin. The salt wall together with the Signalhorn Dome and the Carboniferous pre‐salt structures were again reactivated during post‐Early Cretaceous, in response to regional compressional stresses. During this main tectonic inversion phase, the northwestern and southeastern parts of the salt wall were rejuvenated; however, salt reactivation was minimized towards the interference transfer zone beneath the centre of the salt wall.

Description: The Net-Zero-2050 cluster aims for a national roadmap for net zero CO2 emissions by 2050, including integrated scenario analyses and negative emission technology assessment (see fact sheet Net-Zero-2050 Structure Project 1). This national target to substantially reduce national CO2 emissions by 2050 stems from the objective to comply with the global long-term temperature goal of well below 2°C of the Paris Agreement (UNFCCC, 2015). Within the cluster it is therefore important to decide on an approach for deriving a national remaining carbon budget from global emissions trajectories in agreement with the Paris Climate Agreement’s longterm temperature goal (UNFCCC, 2015). Allocating national carbon budgets is a balance of environmental effectiveness, equity, national capacity and ability, political feasibility, economic efficiency and technical requirements (Gignac and Matthews, 2015; Höhne et al., 2003; 2014). Given Germany’s capacity and abilities, we decided to follow a sustainable growth trajectory with a convergence phase to equal-per-capita CO2 emissions by 2035, and a net zero CO2 emissions trajectory after 2050 until the end of the century. This approach leads to a remaining Germany CO2 budget of 9 GtCO2 (from 1st January 2018 to 2050 and 2100), which we propose to be used across the Net-Zero-2050 cluster. The remaining carbon budget will serve as a target to be used in all work packages in a concerted way, either qualitatively or quantitatively, and in accordance with other work packages (see also fact sheet Net-Zero-2050 Energy Scenario Approach). The calculated budget is at the lower end of the national budget if allocated by the grandfathering approach (emissions are allocated with respect to today’s emissions shares: 5.5-13.1 GtCO2), but slightly higher than the highest estimate of an equal-per-capita remaining carbon budget (emissions are allocated with respect to Germany’s share of the global population: 3.5-8.4 GtCO2) The 9 GtCO2 national remaining CO2 budget, 6.9 GtCO2 from 1st January 2021, will need to be broken down by category (e.g. energy, land use, industrial processes, and man-made sinks and sources; see Gap Analysis Report) in order to provide a consistent approach across work packages.

Description: Doppler wind lidars (DWLs) have increasingly been used over the last decade to derive the mean wind in the atmospheric boundary layer. DWLs allow the determination of wind vector profiles with high vertical resolution and provide an alternative to classic meteorological tower observations. They also receive signals from altitudes higher than a tower and can be set up flexibly in any power-supplied location. In this work, we address the question of whether and how wind gusts can be derived from DWL observations. The characterization of wind gusts is one central goal of the Field Experiment on Sub-Mesoscale Spatio-Temporal Variability in Lindenberg (FESSTVaL). Obtaining wind gusts from a DWL is not trivial because a monostatic DWL provides only a radial velocity per line of sight, i.e., only one component of a three-dimensional vector, and measurements in at least three linearly independent directions are required to derive the wind vector. Performing them sequentially limits the achievable time resolution, while wind gusts are short-lived phenomena. This study compares different DWL configurations in terms of their potential to derive wind gusts. For this purpose, we develop a new wind retrieval method that is applicable to different scanning configurations and various time resolutions. We test eight configurations with StreamLine DWL systems from HALO Photonics and evaluate gust peaks and mean wind over 10 min at 90 m a.g.l. against a sonic anemometer at the meteorological tower in Falkenberg, Germany. The best-performing configuration for retrieving wind gusts proves to be a fast continuous scanning mode (CSM) that completes a full observation cycle within 3.4 s. During this time interval, about 11 radial Doppler velocities are measured, which are then used to retrieve single gusts. The fast CSM configuration was successfully operated over a 3-month period in summer 2020. The CSM paired with our new retrieval technique provides gust peaks that compare well to classic sonic anemometer measurements from the meteorological tower.

Description: Despite the implication of aerosols for the radiation budget, there are persistent differences in data for the aerosol optical depth (τ) for 1998–2019. This study presents a comprehensive evaluation of the large-scale spatio-temporal patterns of mid-visible τ from modern data sets. In total, we assessed 94 different global data sets from eight satellite retrievals, four aerosol-climate model ensembles, one operational ensemble product, two reanalyses, one climatology and one merged satellite product. We include the new satellite data SLSTR and aerosol-climate simulations from the Coupled Model Intercomparison Project Phase 6 (CMIP6) and the Aerosol Comparisons between Observations and Models Phase 3 (AeroCom-III). Our intercomparison highlights model differences and observational uncertainty. Spatial mean τ for 60°N – 60°S ranges from 0.124 to 0.164 for individual satellites, with a mean of 0.14. Averaged τ from aerosol-climate model ensembles fall within this satellite range, but individual models do not. Our assessment suggests no systematic improvement compared to CMIP5 and AeroCom-I. Although some regional biases have been reduced, τ from both CMIP6 and AeroCom-III are for instance substantially larger along extra-tropical storm tracks compared to the satellite products. The considerable uncertainty in observed τ implies that a model evaluation based on a single satellite product might draw biased conclusions. This underlines the need for continued efforts to improve both model and satellite estimates of τ, for example, through measurement campaigns in areas of particularly uncertain satellite estimates identified in this study, to facilitate a better understanding of aerosol effects in the Earth system. Key Points: - Present-day patterns in aerosol optical depth differ substantially between 94 modern global data sets - The range in spatial means from individual satellites is −11% to +17% of the multi-satellite mean - Spatial means from climate model intercomparison projects fall within the satellite range but strong regional differences are identified

Description: Solar radiation received at the Earth's surface (Rs) is comprised of two components, the direct radiation (Rd) and the diffuse radiation (Rf). Rd, the direct beam from the sun, is essential for concentrated solar power generation. Rf, scattered by atmospheric molecules, aerosols, or cloud droplets, has a fertilization effect on plant photosynthesis. But how Rd and Rf change diurnally is largely unknown owing to the lack of long-term measurements. Taking advantage of 22 years of homogeneous hourly surface observations over China, this study documents the climatological means and evolutions in the diurnal cycles of Rd and Rf since 1993, with an emphasis on their implications for solar power and agricultural production. Over the solar energy resource region, we observe a loss of Rd which is relatively large near sunrise and sunset at low solar elevation angles when the sunrays pass through the atmosphere on a longer pathway. However, the concentrated Rd energy covering an average 10-hr period around noon during a day is relatively unaffected. Over the agricultural crop resource region, the large amounts of clouds and aerosols scattering more of the incoming light result in Rf taking the main proportion of Rs during the whole day. Rf resources and their fertilization effect in the main crop region of China further enhances since 1993 over almost all hours of the day. Key Points: - The loss of direct radiation over China since 1993 is relatively large at sunrise and sunset with little effect on solar power generation - The diffuse component dominates solar radiation normally near sunrise and sunset, but for the whole day over the main sown area of China - The diffuse fraction is further enhanced in the main sown area of China over almost all hours of the day since 1993

Description: A challenge of an energy system that nowadays more strongly depends on wind power generation is the spatial and temporal variability in winds. Nocturnal low-level jets (NLLJs) are typical wind phenomena defined as a maximum in the vertical profile of the horizontal wind speed. A NLLJ has typical core heights of 50–500 m a.g.l. (above ground level), which is in the height range of most modern wind turbines. This study presents NLLJ analyses based on new observations from Doppler wind lidars. The aim is to characterize the temporal and spatial variability in NLLJs on the mesoscale and to quantify their impacts on wind power generation. The data were collected during the Field Experiment on Submesoscale Spatio-Temporal Variability (FESSTVaL) campaign from June to August 2020 in Lindenberg and Falkenberg (Germany), located at about 6 km from each other. Both sites have seen NLLJs in about 70 % of the nights with half of them lasting for more than 3 h. Events longer than 6 h occurred more often simultaneously at both sites than shorter events, indicating the mesoscale character of very long NLLJs. Very short NLLJs of less than 1 h occurred more often in Lindenberg than Falkenberg, indicating more local influences on the wind profile. We discussed different meteorological mechanisms for NLLJ formation and linked NLLJ occurrences to synoptic weather patterns. There were positive and negative impacts of NLLJs on wind power that we quantified based on the observational data. NLLJs increased the mean power production by up to 80 % and were responsible for about 25 % of the power potential during the campaign. However, the stronger shear in the rotor layer during NLLJs can also have negative impacts. The impacts of NLLJs on wind power production depended on the relative height between the wind turbine and the core of the NLLJ. For instance, the mean increase in the estimated power production during NLLJ events was about 30 % higher for a turbine at 135 m a.g.l. compared to one at 94 m a.g.l. Our results imply that long NLLJs have an overall stronger impact on the total power production, while short events are primarily relevant as drivers for power ramps.

Description: The reaction between ozone and iodide at the sea surface is now known to be an important part of atmospheric ozone cycling, causing ozone deposition and the release of ozone-depleting reactive iodine to the atmosphere. The importance of this reaction is reflected by its inclusion in chemical transport models (CTMs). Such models depend on accurate sea surface iodide fields, but measurements are spatially and temporally limited. Hence, the ability to predict current and future sea surface iodide fields, i.e. sea surface iodide concentration on a narrow global grid, requires the development of process-based models. These models require a thorough understanding of the key processes that control sea surface iodide. The aim of this study was to explore if there are common features of iodate-to-iodide reduction amongst diverse marine phytoplankton in order to develop models that focus on sea surface iodine and iodine release to the troposphere. In order to achieve this, rates and patterns of changes in inorganic iodine speciation were determined in 10 phytoplankton cultures grown at ambient iodate concentrations. Where possible these data were analysed alongside results from previous studies. Iodate loss and some iodide production were observed in all cultures studied, confirming that this is a widespread feature amongst marine phytoplankton. We found no significant difference in log-phase, cell-normalised iodide production rates between key phytoplankton groups (diatoms, prymnesiophytes including coccolithophores and phaeocystales), suggesting that a phytoplankton functional type (PFT) approach would not be appropriate for building an ocean iodine cycling model. Iodate loss was greater than iodide formation in the majority of the cultures studied, indicating the presence of an as-yet-unidentified “missing iodine” fraction. Iodide yield at the end of the experiment was significantly greater in cultures that had reached a later senescence stage. This suggests that models should incorporate a lag between peak phytoplankton biomass and maximum iodide production and that cell mortality terms in biogeochemical models could be used to parameterise iodide production.

Description: Marine heatwaves have been observed worldwide and are expected to increase in both frequency and intensity due to climate change. Such events may cause ecosystem reconfigurations arising from species range contraction or redistribution, with ecological, economic and social implications. Macrophytes such as the brown seaweed Fucus vesiculosus and the seagrass Zostera marina are foundation species in many coastal ecosystems of the temperate northern hemisphere. Hence, their response to extreme events can potentially determine the fate of associated ecosystems. Macrophyte functioning is intimately linked to the maintenance of photosynthesis, growth and reproduction, and resistance against pathogens, epibionts and grazers. We investigated morphological, physiological, pathological and chemical defence responses of western Baltic Sea F. vesiculosus and Z. marina populations to simulated near‐natural marine heatwaves. Along with (a) the control, which constituted no heatwave but natural stochastic temperature variability (0HW), two treatments were applied: (b) two late‐spring heatwaves (June, July) followed by a summer heatwave (August; 3HW) and (c) a summer heatwave only (1HW). The 3HW treatment was applied to test whether preconditioning events can modulate the potential sensitivity to the summer heatwave. Despite the variety of responses measured in both species, only Z. marina growth was impaired by the accumulative heat stress imposed by the 3HW treatment. Photosynthetic rate, however, remained high after the last heatwave indicating potential for recovery. Only epibacterial abundance was significantly affected in F. vesiculosus. Hence both macrophytes, and in particular F. vesiculosus, seem to be fairly tolerant to short‐term marine heatwaves at least at the intensities applied in this experiment (up to 5°C above mean temperature over a period of 9 days). This may partly be due to the fact that F. vesiculosus grows in a highly variable environment, and may have a high phenotypic plasticity.

Description: Nitric oxide (NO) is a short-lived intermediate of the oceanic nitrogen cycle. However, our knowledge about its production and consumption pathways in oceanic environments is rudimentary. In order to decipher the major factors affecting NO photochemical production, we irradiated several artificial seawater samples as well as 31 natural surface seawater samples in laboratory experiments. The seawater samples were collected during a cruise to the western tropical North Pacific Ocean (WTNP, a N-S section from 36 to 2 degrees N along 146 to 143 degrees E with 6 and 12 stations, respectively, and a W-E section from 137 to 161 degrees E along the Equator with 13 stations) from November 2015 to January 2016. NO photoproduction rates from dissolved nitrite in artificial seawater showed increasing trends with decreasing pH, increasing temperature, and increasing salinity. In contrast, NO photoproduction rates (average: 0.5 +/- 0.2 x 10(-12) mol L-1 s(-1)) in the natural seawater samples from the WTNP did not show any correlations with pH, water temperature, salinity, or dissolved inorganic nitrite concentrations. The flux induced by NO photoproduction in the WTNP (average: 13 x 10(-12) mol M-2 S-1) was significantly larger than the NO air-sea flux density (average: 1.8 x 10(-12) Mol M-2 S-1), indicating a further NO loss process in the surface layer.

Description: Variations of the solar spectral irradiance (SSI) with the 11-year sunspot cycle have been shown to have a significant impact on temperatures and the mixing ratios of atmospheric constituents in the stratosphere and mesosphere. Uncertainties in modelling the effects of SSI variations arise from uncertainties in the empirical models reconstructing the prescribed SSI data set as well as from uncertainties in the chemistry-climate model (CCM) formulation. In this study CCM simulations with the ECHAM MESSy Atmospheric Chemistry (EMAC) model and the Community Earth System Model 1 (CESM1) – Whole Atmosphere Chemistry Climate Model (WACCM) have been performed to quantify the uncertainties of the solar responses in chemistry and dynamics that are due to the usage of five different SSI data sets or the two CCMs. We apply a two-way analysis of variance (ANOVA) to separate the influence of the SSI data sets and the CCMs on the variability of the solar response in shortwave heating rates, temperature and ozone. The ANOVA identifies the SSI data set with the strongest influence on the variability of the solar signal in shortwave heating rates in the upper mesosphere and in the upper stratosphere/lower mesosphere. The strongest influence on the variability of the solar signal in ozone and temperature is identified in the upper stratosphere/lower mesosphere. The largest influence of the CCMs on variability of the solar responses can be identified in the upper mesosphere. The solar response in the lower stratosphere also depends on the CCM used, especially in the tropics and northern hemispheric subtropics and mid latitudes, where the model dynamics modulate the solar responses.

Description: Drill sites in the southern Bay of Bengal at 3°N 91°E (International Ocean Discovery Program Expedition 362) have sampled for the first time a complete section of the Nicobar Fan and below to the oceanic crust. This generally overlooked part of the Bengal–Nicobar Fan System may provide new insights into uplift and denudation rates of the Himalayas and Tibetan Plateau. The Nicobar Fan comprises sediment gravity‐flow deposits, mostly turbidites, that alternate with hemipelagite drapes and pelagite intervals of varying thicknesses. The decimetre‐thick to metre‐thick oldest pre‐fan sediments (limestones/chalks) dated at 69 Ma are overlain by volcanic material and slowly accumulated pelagites (0.5 g cm−2 kyr−1). At Expedition 362 Site U1480, terrigenous input began in the early Miocene at ca 22.5 Ma as muds, overlain by very thin‐bedded and thin‐bedded muddy turbidites at ca 19.5 Ma. From 9.5 Ma, sand content and sediment supply sharply increase (from 1–5 to 10–50 g cm−2 kyr−1). Despite the abundant normal faulting in the Nicobar Fan compared with the Bengal Fan, it offers a better‐preserved and more homogeneous sedimentary record with fewer unconformities. The persistent connection between the two fans ceased at 0.28 Ma when the Nicobar Fan became inactive. The Nicobar Fan is a major sink for Himalaya‐derived material. This study presents integrated results of International Ocean Discovery Program Expedition 362 with older Deep Sea Drilling Project/Ocean Drilling Program/International Ocean Discovery Program sites that show that the Bengal–Nicobar Fan System experienced successive large‐scale avulsion processes that switched sediment supply between the Bengal Fan (middle Miocene and late Pleistocene) and the Nicobar Fan (late Miocene to early Pleistocene). A quantitative analysis of the submarine channels of the Nicobar Fan is also presented, including their stratigraphic frequency, showing that channel size/area and abundance peaked at ca 2 to 3 Ma, but with a distinct low at 3 to 7 Ma: the intervening stratigraphic [sub]unit was a time of reduced sediment accumulation rates.

Description: The generation of reliable age models for palaeoenvironmental and archaeological records in the Eurasian Arctic is often problematic when using conventional dating techniques. Tephrochronology can potentially improve the chronologies of such records and synchronise disparate sedimentary archives. However, to date, systematic tephra studies are lacking for this region. This paper presents the first cryptotephra data from the White Sea region (northwestern Russia) based on a peat core spanning the past ~1800 years. We identify seven geochemical glass populations that derive from six Icelandic volcanoes and correlate four of them to north European tephra isochrons; these include Askja ad 1875, the basaltic component of the ad 877 Landnám tephra, and tephras BTD‐15 (c. ad 1750–1650) and SL‐2/SB‐2 (ad 803–767) from unknown eruptions of Katla and Snæfellsjökull, respectively. The remaining three populations originate from Grímsvötn, Hekla and Katla; however, their attribution to individual eruptions remains ambiguous. These findings highlight the potential to extend the Late Holocene tephrochronological framework of northern Europe to the west Eurasian Arctic. The detection of at least three basaltic tephras in the core suggests that basaltic shards can be transported over larger distances than previously known and that peatlands are well suited to preserve such components.

Description: Climate engineering (CE) measures are increasingly discussed when dealing with the adverse impacts of climate change. While much research has focused on individual methods, few studies attempt to compare and rank the effectiveness of these measures. Furthermore, model uncertainties are seldom acknowledged and lesser still, estimated when CE scenarios are assessed. In this work, we quantify the variance in outcomes due to poorly constrained model parameters under several idealized CE scenarios. The four scenarios considered are (1) warming under the high emission scenario Representative Concentration Pathway 8.5 without CE applied and the same emission scenario with (2) afforestation,(3) solar radiation management, and (4) artificial ocean alkalinization. By considering the parametric uncertainty in model outputs, we demonstrate the problems with comparing these scenarios using a single parameter setting. Using statistical emulation, we estimate the probability distributions of several model outcomes. Based on such distributions, we suggest an approach to ranking the effectiveness of the scenarios considered according to their probability of avoiding climate thresholds.

Description: A new estimate of Agulhas leakage transport is calculated using profiling floats and drifters. Since Richardson's seminal estimate of 15 Sv in 2007, the number of floats and drifters passing through the Agulhas Current has quadrupled. Within uncertainties we find the same leakage percentages as Richardson, with 34% of drifters leaking at the surface and 21% of floats leaking at 1,000 m depth. We find that the drifters tend to follow a northward leakage pathway via the Benguela Current compared to the northwestward leakage pathway of the floats along the Agulhas Ring corridor. We simulate the isobaric and profiling behavior of the floats and drifters using two high resolution models and two offline Lagrangian tracking tools, quantifying for the first time the sampling biases associated with the observations. We find that the isobaric bias cannot be robustly simulated but likely causes an underestimate of observed leakage by one or two Sverdrups. The profiling behavior of the floats causes no significant bias in the leakage. Fitting a simulated vertical leakage profile to the observed leakage percentages from the floats and drifters and using the mean Agulhas transport observed by a moored array at 34°S we find an improved Agulhas leakage transport of 21.3 Sv, with an estimated error of 4.7 Sv. Our new leakage transport is higher primarily because we account for leakage at depths down to 2,000 m, while Richardson considered only the top 1,000 m of the water column.

Description: Upwelling ocean currents associated with oxygen minimum zones (OMZs) supply nutrients fuelling intense marine productivity. Perturbations in the extent and intensity of OMZs are projected in the future, but it is currently uncertain how this will impact fluxes of redox‐sensitive trace metal micronutrients to the surface ocean. Here we report seawater concentrations of Fe, Mn, Co, Cd, and Ni alongside the redox indicator iodide/iodate in the Peruvian OMZ during the 2015 El Niño event. The El Niño drove atypical upwelling of oxygen‐enriched water over the Peruvian Shelf, resulting in oxidized iodine and strongly depleted Fe (II), total dissolved Fe, and reactive particulate Fe concentrations relative to non‐El Niño conditions. Observations of Fe were matched by the redox‐sensitive micronutrients Co and Mn, but not by non‐redox‐sensitive Cd and Ni. These observations demonstrate that oxygenation of OMZs significantly reduces water column inventories of redox‐sensitive micronutrients, with potential impacts on ocean productivity. Plain Language Summary Some trace metals, including iron, are essential micronutrients for phytoplankton growth. However, the solubility of iron is very low under oxygenated conditions. Consequently, restricted iron availability in oxygen‐rich seawater can limit phytoplankton growth in the ocean, including in the Eastern Tropical South Pacific. Under typical conditions, depleted oxygen on the South American continental shelf is generally thought to enhance iron supply to the ocean, fuelling phytoplankton productivity in overlying waters. However, the impact of changes in oxygenation, which are predicted to occur in the future, are not known. The 2015 El Niño event led to unusually high oxygen on the Peruvian shelf, offering a system‐scale test on how oxygen influences seawater iron concentrations. We show that El Niño‐driven oxygenation resulted in marked decreases in iron and other metals sensitive to oxygen (cobalt and manganese), whilst metals not sensitive to oxygen (cadmium and nickel) were unaffected. The measured reductions in iron may have led to decreased phytoplankton productivity.

Description: The causes of the seasonal cycle of vertical turbulent cooling at the base of the mixed layer are assessed using observations from moored buoys in the tropical Atlantic Intertropical Convergence Zone (ITCZ) (4°N, 23°W) and trade wind (15°N, 38°W) regions together with mixing parameterizations and a one-dimensional model. At 4°N the parameterized turbulent cooling rates during 2017–2018 and 2019 agree with indirect estimates from the climatological mooring heat budget residual: both show mean cooling of 25–30 W m (Formula presented.) during November–July, when winds are weakest and the mixed layer is thinnest, and 0–10 W m (Formula presented.) during August–October. Mixing during November–July is driven by variability on multiple time scales, including subdiurnal, near-inertial, and intraseasonal. Shear associated with tropical instability waves (TIWs) is found to generate mixing and monthly mean cooling of 15–30 W m (Formula presented.) during May–July in 2017 and 2019. At 15°N the seasonal cycle of turbulent cooling is out of phase compared to 4°N, with largest cooling of up to 60 W m (Formula presented.) during boreal fall. However, the relationships between wind speed, mixed layer depth, and turbulent mixing are similar: weaker mean winds and a thinner mixed layer in the fall are associated with stronger mixing and turbulent cooling of SST. These results emphasize the importance of seasonal modulations of mixed layer depth at both locations and shear from TIWs at 4°N.

Description: Coastal areas contribute significantly to the emissions of methane (CH4) from the ocean. In order to decipher its temporal variability in the whole water column, dissolved CH4 was measured on a monthly basis at the Boknis Eck Time-series Station (BE) located in the Eckernförde Bay (SW Baltic Sea) from 2006 to 2017. BE has a water depth of about 28 m and dissolved CH4 was measured at six water depths ranging from 0 to 25 m. In general CH4 concentrations increased with depth, indicating a sedimentary release of CH4. Pronounced enhancement of the CH4 concentrations in the bottom layer (15–25 m) was found during February, May–June and October. CH4 was not correlated with Chlorophyll a or O2 over the measurement period. Unusually high CH4 concentrations (of up to 696 nM) were sporadically observed in the upper layer (0–10 m) (e.g. in November 2013 and December 2014) and were coinciding with Major Baltic Inflow (MBI) events. Surface CH4 concentrations were always supersaturated throughout the monitoring period, indicating that the Eckernförde Bay is an intense but highly variable source of atmospheric CH4. We did not detect significant temporal trends in CH4 concentrations or emissions, despite of ongoing environmental changes such as warming and deoxygenation in the Eckernförde Bay. Overall, the CH4 variability at BE is driven by a complex interplay of various biological and physical processes.

Description: The extracellular concentration of H2O2 in surface aquatic environments is controlled by a balance between photochemical production and the microbial synthesis of catalase and peroxidase enzymes to remove H2O2 from solution. In any kind of incubation experiment, the formation rates and equilibrium concentrations of reactive oxygen species (ROSs) such as H2O2 may be sensitive to both the experiment design, particularly to the regulation of incident light, and the abundance of different microbial groups, as both cellular H2O2 production and catalase–peroxidase enzyme production rates differ between species. Whilst there are extensive measurements of photochemical H2O2 formation rates and the distribution of H2O2 in the marine environment, it is poorly constrained how different microbial groups affect extracellular H2O2 concentrations, how comparable extracellular H2O2 concentrations within large-scale incubation experiments are to those observed in the surface-mixed layer, and to what extent a mismatch with environmentally relevant concentrations of ROS in incubations could influence biological processes differently to what would be observed in nature. Here we show that both experiment design and bacterial abundance consistently exert control on extracellular H2O2 concentrations across a range of incubation experiments in diverse marine environments. During four large-scale (〉1000 L) mesocosm experiments (in Gran Canaria, the Mediterranean, Patagonia and Svalbard) most experimental factors appeared to exert only minor, or no, direct effect on H2O2 concentrations. For example, in three of four experiments where pH was manipulated to 0.4–0.5 below ambient pH, no significant change was evident in extracellular H2O2 concentrations relative to controls. An influence was sometimes inferred from zooplankton density, but not consistently between different incubation experiments, and no change in H2O2 was evident in controlled experiments using different densities of the copepod Calanus finmarchicus grazing on the diatom Skeletonema costatum (〈1 % change in [H2O2] comparing copepod densities from 1 to 10 L−1). Instead, the changes in H2O2 concentration contrasting high- and low-zooplankton incubations appeared to arise from the resulting changes in bacterial activity. The correlation between bacterial abundance and extracellular H2O2 was stronger in some incubations than others (R2 range 0.09 to 0.55), yet high bacterial densities were consistently associated with low H2O2. Nonetheless, the main control on H2O2 concentrations during incubation experiments relative to those in ambient, unenclosed waters was the regulation of incident light. In an open (lidless) mesocosm experiment in Gran Canaria, H2O2 was persistently elevated (2–6-fold) above ambient concentrations; whereas using closed high-density polyethylene mesocosms in Crete, Svalbard and Patagonia H2O2 within incubations was always reduced (median 10 %–90 %) relative to ambient waters.

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: A major surface circulation feature of the Arctic Ocean is the Transpolar Drift (TPD), a current that transports river‐influenced shelf water from the Laptev and East Siberian Seas toward the center of the basin and Fram Strait. In 2015, the international GEOTRACES program included a high‐resolution pan‐Arctic survey of carbon, nutrients, and a suite of trace elements and isotopes (TEIs). The cruises bisected the TPD at two locations in the central basin, which were defined by maxima in meteoric water and dissolved organic carbon concentrations that spanned 600 km horizontally and ~25‐50 m vertically. Dissolved TEIs such as Fe, Co, Ni, Cu, Hg, Nd, and Th, which are generally particle‐reactive but can be complexed by organic matter, were observed at concentrations much higher than expected for the open ocean setting. Other trace element concentrations such as Al, V, Ga, and Pb were lower than expected due to scavenging over the productive East Siberian and Laptev shelf seas. Using a combination of radionuclide tracers and ice drift modeling, the transport rate for the core of the TPD was estimated at 0.9 ± 0.4 Sv (106 m3 s‐1). This rate was used to derive the mass flux for TEIs that were enriched in the TPD, revealing the importance of lateral transport in supplying materials beneath the ice to the central Arctic Ocean and potentially to the North Atlantic Ocean via Fram Strait. Continued intensification of the Arctic hydrologic cycle and permafrost degradation will likely lead to an increase in the flux of TEIs into the Arctic Ocean.

Description: Cichlid fishes’ famous diversity in body coloration is accompanied by a highly diverse and complex visual system. Although cichlids possess an unusually high number of seven cone opsin genes, they express only a subset of these during their ontogeny, accounting for their astonishing interspecific variation in visual sensitivities. Much of this diversity is thought to have been shaped by natural selection as cichlids inhabit a variety of habitats with distinct light environments. Also, sexual selection might have contributed to the observed visual diversity, and sexual dimorphism in coloration potentially co‐evolved with sexual dimorphism in opsin expression. We investigated sex‐specific opsin expression of several cichlids from Africa and the Neotropics and collected and integrated datasets on sex‐specific body coloration, species‐specific visual sensitivities, lens transmission and habitat light properties for some of them. We comparatively analyzed this wide range of molecular and ecological data, illustrating how integrative approaches can address specific questions on the factors and mechanisms driving diversification, and the evolution of cichlid vision in particular. We found that both sexes expressed opsins at the same levels ‐ even in sexually dimorphic cichlid species – which argues against coevolution of sexual dichromatism and differences in sex‐specific visual sensitivity. Rather, a combination of environmental light properties and body coloration shaped the diversity in spectral sensitivities among cichlids. We conclude that although cichlids are particularly colorful and diverse and often sexually dimorphic, it would appear that natural rather than sexual selection is a more powerful force driving visual diversity in this hyper‐diverse lineage.

Description: Key Points: - The Kamchatka arc lavas show across-arc variations in chalcophile elements, suggesting that the amount of fluid decreases with depth - Slab-derived fluids have a negligible contribution to the Li budget of the Kamchatka arc lavas - CKD lavas have high U/Th, Li/Y, La/Sm and B/Nb ratios, indicating that lawsonite breakdown reaction dominates the water release Chalcophile elements and lithium (Li) isotopes were measured on lavas from a 220 km transect across the Kamchatka arc in order to investigate the fluid variations below arc volcanoes and to trace the geochemical behaviour of Li in convergent plate margins. From the Eastern Volcanic Front (EVF), through the Central Kamchatka Depression (CKD), into the Sredinny Range (SR) volcanic zones, chalcophile element ratios (e.g., As/Ce and Sb/Ce) show clear across‐arc variations, decreasing (e.g., As/Ce: 0.20 to 0.03 and Sb/Ce: 0.013 to 0.002) with increasing depth above the slab (110 to 400 km). This clearly indicates a gradually decreasing influx of slab‐derived fluids added to the mantle wedge as the slab subducts below Kamchatka. In addition, the anomalously high U/Th, La/Sm and B/Nb ratios in the CKD lavas suggest lawsonite breakdown reaction dominates the fluid release in this area. However, Li/Y (0.07 to 1.78) and δ7Li (+1.8 to +5.4‰, with an exception of +8.6‰ in CKD) show limited variations and values similar to the MORB mantle. A dehydration model suggests that slab‐derived fluids, which are characterized by high Li concentration and high δ7Li, do not control the Li budget in Kamchatka arc lavas. Therefore, the isotopic heavy Li from slab‐derived fluids likely equilibrates in the sub‐arc mantle, which acts as a buffer for Li systematics. In addition, based on the Li isotopic signatures of Klyuchevskoy volcano, our study demonstrates insignificant Li isotopic fractionation during mantle melting and subsequent differentiation.

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: The Endeavour Segment of the Juan de Fuca Ridge is well known for its abundance of hydrothermal vents and chimneys. One-meter scale multibeam mapping data collected by an autonomous undersea vehicle revealed 572 chimneys along the central 14 km of the segment, although only 47 are named and known to be active. Hydrothermal deposits are restricted to the axial graben and the near-rims of the graben above a seismically mapped axial magma lens. The sparse eruptive activity on the segment during the last 4,300 years has not buried inactive chimneys, as occurs at more magmatically robust mid-ocean ridges.

Description: Cold-water corals (CWCs) constitute important deep-water ecosystems that are under increasing environmental pressure due to ocean acidification and global warming. The sensitivity of these deep-water ecosystems to environmental change is demonstrated by abundant paleorecords drilled through CWC mounds that reveal characteristic alterations between rapid formation and dormant or erosive phases. Previous studies have identified several central parameters for driving or inhibiting CWC growth such as food supply, oxygenation, and the carbon saturation state of bottom water, yet there are still large uncertainties about the relative importance of the different environmental parameters. To advance this debate we have performed a multiproxy study on a sediment core retrieved from the 25 m high Bowie Mound, located at 866 m water depth on the continental slope off southeastern Brazil, a structure built up mainly by the CWC Solenosmilia variabilis. Our results indicate a multifactorial control on CWC growth at Bowie Mound during the past ∼ 160 kyr, which reveals distinct formation pulses during northern high-latitude glacial cold events (Heinrich stadials, HSs) largely associated with anomalously strong monsoonal rainfall over the continent. The ensuing enhanced runoff elevated the terrigenous nutrient and organic-matter supply to the continental margin and likely boosted marine productivity. The dispersal of food particles towards the CWC colonies during HSs was facilitated by the highly dynamic hydraulic conditions along the continental slope that prevailed throughout glacial periods. These conditions caused the emplacement of a pronounced nepheloid layer above Bowie Mound, thereby aiding the concentration and along-slope dispersal of organic matter. Our study thus emphasizes the impact of continental climate variability on a highly vulnerable deep-marine ecosystem.

Description: Variations in the solar spectral irradiance (SSI) with the 11-year sunspot cycle have been shown to have a significant impact on temperatures and the mixing ratios of atmospheric constituents in the stratosphere and mesosphere. Uncertainties in modelling the effects of SSI variations arise from uncertainties in the empirical models reconstructing the prescribed SSI data set as well as from uncertainties in the chemistry–climate model (CCM) formulation. In this study CCM simulations with the ECHAM/MESSy Atmospheric Chemistry (EMAC) model and the Community Earth System Model 1 (CESM1)–Whole Atmosphere Chemistry Climate Model (WACCM) have been performed to quantify the uncertainties of the solar responses in chemistry and dynamics that are due to the usage of five different SSI data sets or the two CCMs. We apply a two-way analysis of variance (ANOVA) to separate the influence of the SSI data sets and the CCMs on the variability of the solar response in shortwave heating rates, temperature, and ozone. The solar response is derived from climatological differences of time slice simulations prescribing SSI for the solar maximum in 1989 and near the solar minimum in 1994. The SSI values for the solar maximum of each SSI data set are created by adding the SSI differences between November 1994 and November 1989 to a common SSI reference spectrum for near-solar-minimum conditions based on ATLAS-3 (Atmospheric Laboratory of Applications and Science-3). The ANOVA identifies the SSI data set with the strongest influence on the variability of the solar response in shortwave heating rates in the upper mesosphere and in the upper stratosphere–lower mesosphere. The strongest influence on the variability of the solar response in ozone and temperature is identified in the upper stratosphere–lower mesosphere. However, in the region of the largest ozone mixing ratio, in the stratosphere from 50 to 10 hPa, the SSI data sets do not contribute much to the variability of the solar response when the Spectral And Total Irradiance REconstructions-T (SATIRE-T) SSI data set is omitted. The largest influence of the CCMs on variability of the solar responses can be identified in the upper mesosphere. The solar response in the lower stratosphere also depends on the CCM used, especially in the tropics and northern hemispheric subtropics and mid-latitudes, where the model dynamics modulate the solar responses. Apart from the upper mesosphere, there are also regions where the largest fraction of the variability of the solar response is explained by randomness, especially for the solar response in temperature.

Description: Although the core velocity of the Atlantic North Equatorial Undercurrent (NEUC) is low (0.1−0.3 m s−1), it has been suggested to act as an important oxygen supply route towards the oxygen minimum zone in the eastern tropical North Atlantic. For the first time, the intraseasonal to interannual NEUC variability and its impact on oxygen are investigated based on shipboard and moored velocity observations around 5°N, 23°W. In contrast to previous studies that were mainly based on models or hydrographic data, we find hardly any seasonal cycle of NEUC transports in the central Atlantic. The NEUC transport variability is instead dominated by sporadic intraseasonal events. Only some of these events are associated with high oxygen levels suggesting an occasional eastward oxygen supply by NEUC transport events. Nevertheless, they likely contribute to the local oxygen maximum in the mean shipboard section along 23°W at the NEUC core position.

Description: In Kiel, in the north of Germany, marine research is rooted in a lively research community hosted mainly at Kiel University and the GEOMAR Helmholtz Centre. While the ratio of women and men is more or less balanced on all qualification levels with mainly nonpermanent junior positions, women are generally underrepresented in leading research positions. The problem of gender imbalance and inequality has been well-known for a long time. Especially in the last decade, however, manifold efforts were initiated to improve gender equality on a political and institutional level as well as within the research community itself. In our article we focus on the gender equality activities of the two large externally funded marine sciences research alliances: the Cluster of Excellence “The Future Ocean” and the Collaborative Research Centre 754 “Climate–Biogeochemistry Interactions in the Tropical Ocean”. For about a decade they offered both financial provisions and a structural framework to tackle the problem of women's underrepresentation in science and came up with innovative measures. In the following case study, we not only introduce the situation of women in marine sciences in Kiel and the structural arrangement to improve gender equality in general, but we also discuss three specific measures developed within the two collaborative research projects in detail: (i) the mentoring program via:mento_ocean for female postdocs, (ii) hiring policies integrating a gender quota for recruiting postdoctoral researchers and (iii) a code of conduct. Based on these best-practice examples we can show that progress towards gender equality has been made despite some obstacles faced when implementing the measures. This was especially the case for attracting female researchers to work in Kiel marine sciences and bringing the relevance of the topic to the surface of debates within the community. Looking at gender equality activities from a managerial point of view, we conclude from the situation in Kiel, where external funding for both research alliances ended in 2019, that even time-bound activities can initiate change. Initiatives developed by the marine sciences community were taken up by other research groups and inspired new activities at the level of the institutions

Description: The Zero Emissions Commitment (ZEC) is the change in global mean temperature expected to occur following the cessation of net CO2 emissions and as such is a critical parameter for calculating the remaining carbon budget. The Zero Emissions Commitment Model Intercomparison Project (ZECMIP) was established to gain a better understanding of the potential magnitude and sign of ZEC, in addition to the processes that underlie this metric. A total of 18 Earth system models of both full and intermediate complexity participated in ZECMIP. All models conducted an experiment where atmospheric CO2 concentration increases exponentially until 1000 PgC has been emitted. Thereafter emissions are set to zero and models are configured to allow free evolution of atmospheric CO2 concentration. Many models conducted additional second-priority simulations with different cumulative emission totals and an alternative idealized emissions pathway with a gradual transition to zero emissions. The inter-model range of ZEC 50 years after emissions cease for the 1000 PgC experiment is −0.36 to 0.29 ∘C, with a model ensemble mean of −0.07 ∘C, median of −0.05 ∘C, and standard deviation of 0.19 ∘C. Models exhibit a wide variety of behaviours after emissions cease, with some models continuing to warm for decades to millennia and others cooling substantially. Analysis shows that both the carbon uptake by the ocean and the terrestrial biosphere are important for counteracting the warming effect from the reduction in ocean heat uptake in the decades after emissions cease. This warming effect is difficult to constrain due to high uncertainty in the efficacy of ocean heat uptake. Overall, the most likely value of ZEC on multi-decadal timescales is close to zero, consistent with previous model experiments and simple theory.

Description: The development and application of geochemical techniques to identify redox conditions in modern and ancient aquatic environments has intensified over recent years. Iron (Fe) speciation has emerged as one of the most widely used procedures to distinguish different redox regimes in both the water column and sediments, and is the main technique used to identify oxic, ferruginous (anoxic, Fe(II) containing) and euxinic (anoxic, sulfidic) water column conditions. However, an international sediment reference material has never been developed. This has led to concern over the consistency of results published by the many laboratories that now utilise the technique. Here, we report an interlaboratory comparison of four Fe speciation reference materials for palaeoredox analysis, which span a range of compositions and reflect deposition under different redox conditions. We provide an update of extraction techniques used in Fe speciation and assess the effects of both test portion mass, and the use of different analytical procedures, on the quantification of different Fe fractions in sedimentary rocks. While atomic absorption spectroscopy and inductively coupled plasma‐optical emission spectrometry produced comparable Fe measurements for all extraction stages, the use of ferrozine consistently underestimated Fe in the extraction step targeting mixed ferrous–ferric minerals such as magnetite. We therefore suggest that the use of ferrozine is discontinued for this Fe pool. Finally, we report the combined data of four independent Fe speciation laboratories to characterise the Fe speciation composition of the reference materials. These reference materials are available to the community to provide an essential validation of in‐house Fe speciation measurements.

Description: Resting egg production is considered the most common form of dormancy in aquatic invertebrates. Given that many taxa at least partially terminate resting egg state using environmental cues, knowledge on environmental drivers of hatching success is important, particularly within the context of climate change and environmental degradation. Fairy shrimp (anostracans) are temporary wetland specialists that are reliant on resting egg production for population persistence. Temporary wetlands are common in many arid regions projected to experience increases in temperature, and in areas often compromised by human‐mediated activities. In this study, we assessed the combined effects of light and temperature on the hatching success of Streptocephalus cafer (Anostraca) dormant eggs from temporary wetlands in an arid environment. Both temperature and light altered hatching success, with emergent effects evident. Light caused a significant threefold increase in hatching success overall, while temperature effects were non‐linear, with hatching optimised at 27°C, and especially under light conditions. These results are discussed within the context of shifting climates and disturbances to temporary wetland ecosystems.

Description: Most molluscs possess shells, constructed from a vast array of microstructures and architectures. The fully formed shell is composed of calcite or aragonite. These CaCO3 crystals form complex biocomposites with proteins, which although typically less than 5% of total shell mass, play significant roles in determining shell microstructure. Despite much research effort, large knowledge gaps remain in how molluscs construct and maintain their shells, and how they produce such a great diversity of forms. Here we synthesize results on how shell shape, microstructure, composition and organic content vary among, and within, species in response to numerous biotic and abiotic factors. At the local level, temperature, food supply and predation cues significantly affect shell morphology, whilst salinity has a much stronger influence across latitudes. Moreover, we emphasize how advances in genomic technologies [e.g. restriction site-associated DNA sequencing (RAD-Seq) and epigenetics] allow detailed examinations of whether morphological changes result from phenotypic plasticity or genetic adaptation, or a combination of these. RAD-Seq has already identified single nucleotide polymorphisms associated with temperature and aquaculture practices, whilst epigenetic processes have been shown significantly to modify shell construction to local conditions in, for example, Antarctica and New Zealand. We also synthesize results on the costs of shell construction and explore how these affect energetic trade-offs in animal metabolism. The cellular costs are still debated, with CaCO3 precipitation estimates ranging from 1–2 J/mg to 17–55 J/mg depending on experimental and environmental conditions. However, organic components are more expensive (~29 J/mg) and recent data indicate transmembrane calcium ion transporters can involve considerable costs. This review emphasizes the role that molecular analyses have played in demonstrating multiple evolutionary origins of biomineralization genes. Although these are characterized by lineage-specific proteins and unique combinations of co-opted genes, a small set of protein domains have been identified as a conserved biomineralization tool box. We further highlight the use of sequence data sets in providing candidate genes for in situ localization and protein function studies. The former has elucidated gene expression modularity in mantle tissue, improving understanding of the diversity of shell morphology synthesis. RNA interference (RNAi) and clustered regularly interspersed short palindromic repeats - CRISPR-associated protein 9 (CRISPR-Cas9) experiments have provided proof of concept for use in the functional investigation of mollusc gene sequences, showing for example that Pif (aragonite-binding) protein plays a significant role in structured nacre crystal growth and that the Lsdia1 gene sets shell chirality in Lymnaea stagnalis. Much research has focused on the impacts of ocean acidification on molluscs. Initial studies were predominantly pessimistic for future molluscan biodiversity. However, more sophisticated experiments incorporating selective breeding and multiple generations are identifying subtle effects and that variability within mollusc genomes has potential for adaption to future conditions. Furthermore, we highlight recent historical studies based on museum collections that demonstrate a greater resilience of molluscs to climate change compared with experimental data. The future of mollusc research lies not solely with ecological investigations into biodiversity, and this review synthesizes knowledge across disciplines to understand biomineralization. It spans research ranging from evolution and development, through predictions of biodiversity prospects and future-proofing of aquaculture to identifying new biomimetic opportunities and societal benefits from recycling shell products.

Description: In continental settings, seismic failure is generally restricted to crustal depth. Crustal structure is therefore an important proxy to evaluate seismic hazard of continental fault systems. Here we present a seismic velocity model across the Gibraltar Arc System, from the Eurasian Betics Range (South Iberian margin), across offshore East Alboran and Pytheas (African margin) basins, and ending onshore in North Morocco. Our results reveal the nature and configuration of the crust supporting the coexistence of three different crustal domains: the continental crust of the Betics, the continental crust of the Pytheas Basin (south Alboran Basin) and onshore Morocco, and a distinct domain formed of magmatic arc crust under the East Alboran Basin. The magmatic arc under the East Alboran Basin is characterized by a velocity structure containing a relatively high‐velocity lower crust (~7 km/s) bounded at the top and base by reflections. The lateral extension of this crust is mapped integrating a second perpendicular wide‐angle seismic profile along the Eastern Alboran basin, together with basement samples, multibeam bathymetry, and a grid of deep‐penetrating multichannel seismic profiles. The transition between crustal domains is currently unrelated to extensional and magmatic processes that formed the basin. The abrupt transition zones between the different crustal domains support that they are bounded by crustal‐scale active fault systems that reactivate inherited structures. Seismicity in the area is constrained to upper‐middle crust depths, and most earthquakes nucleate outside of the magmatic arc domain. Key Points New velocity model reveals the lithospheric structure under the Betics (South Iberia), the Alboran Basin and the North African margin The East Alboran Basin is floored by magmatic arc crust, while the southern area of the Alboran Basin is floored by continental crust Seismic activity is constrained to the upper‐middle continental crust. Crustal domains are likely bounded by active faults

Description: We constrain the lithospheric mantle density of the North China Craton (NCC) at both in situ and standard temperature‐pressure (STP) conditions from gravity data. The lithosphere‐asthenosphere boundary (LAB) depth is constrained by our new thermal model, which is based on a new regional heat flow data set and a recent regional crustal model NCcrust. The new thermal model shows that the thermal lithosphere thickness is 〈120 km in most of the NCC, except for the northern and southern parts with the maximum depth of 170 km. The gravity calculations reveal a highly heterogeneous density structure of the lithospheric mantle with in situ and STP values of 3.22–3.29 and 3.32–3.40 g/cm3, respectively. Thick and reduced‐density cratonic‐type lithosphere is preserved mostly in the southern NCC. Most of the Eastern Block has a thin (90–140 km) and high‐density lithospheric mantle. Most of the Western Block has a high‐density lithospheric mantle and a thin (80–110 km) lithosphere typical of Phanerozoic regions, which suggests that the Archean lithosphere is no longer present there. We conclude that in almost the entire NCC the lithosphere has lost its cratonic characteristics by geodynamic processes that include, but are not limited to, the Paleozoic closure of the Paleo‐Asian Ocean in the north, the Mesozoic Yangtze Craton flat subduction in the south, the Mesozoic Pacific subduction in the east, the Cenozoic remote response to the Indian‐Eurasian collision in the west, and the Cenozoic extensional tectonics (possibly associated with the slab roll‐back) in the center.

Description: Most animals, including sponges (Porifera), have species-specific microbiomes. Which genetic or environmental factors play major roles structuring the microbial community at the intraspecific level in sponges is, however, largely unknown. In this study, we tested whether geographic location or genetic structure of conspecific sponges influences their microbial assembly. For that, we used three sponge species with different rates of gene flow, and collected samples along their entire distribution range (two from the Mediterranean and one from the Southern Ocean) yielding a total of 393 samples. These three sponge species have been previously analysed by microsatellites or single nucleotide polymorphisms, and here we investigate their microbiomes by amplicon sequencing of the microbial 16S rRNA gene. The sponge Petrosia ficiformis, with highly isolated populations (low gene flow), showed a stronger influence of the host genetic distance on the microbial composition than the spatial distance. Host-specificity was therefore detected at the genotypic level, with individuals belonging to the same host genetic cluster harbouring more similar microbiomes than distant ones. On the contrary, the microbiome of Ircinia fasciculata and Dendrilla antarctica - both with weak population structure (high gene flow) - seemed influenced by location rather than by host genetic distance. Our results suggest that in sponge species with high population structure, the host genetic cluster influence the microbial community more than the geographic location.

Description: Many coastal oceans experience not only increased loads of nutrients but also changes in the stoichiometry of nutrient supply. Excess supply of nitrogen and stable or decreased supply of silicon lower silicon to nitrogen (Si:N) ratios, which may decrease diatom proportion in phytoplankton. To examine how Si:N ratios affect plankton community composition and food web structure, we performed a mesocosm experiment where we manipulated Si:N ratios and copepod abundance in a Baltic Sea plankton community. In high Si:N treatments, diatoms dominated. Some of them were likely spared from grazing unexpectedly resulting in higher diatom biomass under high copepod grazing. With declining Si:N ratios, dinoflagellates became more abundant under low and picoplankton under high copepod grazing. This altered plankton food web structure: under high Si:N ratios, edible diatoms were directly accessible food for copepods, while under low Si:N ratios, microzooplankton and phago-mixotrophs (mixoplankton) were a more important food source for mesograzers. The response of copepods to changes in the phytoplankton community was complex and copepod density-dependent. We suggest that declining Si:N ratios favor microzoo- and mixoplankton leading to increased complexity of planktonic food webs. Consequences on higher trophic levels will, however, likely be moderated by edibility, nutritional value or toxicity of dominant phytoplankton species.

Description: Observations and reanalysis products are used to investigate the substantial weakening in the southeastern tropical Atlantic sea‐surface temperature (SST) variability since 2000. Relative to 1982‐1999, the March‐April‐May SST variability in the Angola‐Benguela area (ABA) has decreased by more than 30 %. Both equatorial remote forcing and local forcing are known to play an important role in driving SST variability in the ABA. Compared to 1982‐1999, since 2000 equatorial remote forcing had less influence on ABA SSTs whereas local forcing has become more important. In particular, the robust correlation that existed between the equatorial zonal wind stress and the ABA SSTs has substantially weakened, suggesting less influence of Kelvin waves on ABA SSTs. Moreover, the strong correlation linking the South Atlantic Anticyclone and the ABA SSTs has reduced. Finally, multidecadal surface warming of the ABA could also have played a role in the weakening of the interannual SST variability.

Description: Quantification and attribution of the food web changes associated with the invasion of non-indigenous species in the marine realm often remain a challenge. One of the pelagic non-indigenous species of concern in the recent history of aquatic bioinvasions is the predatory cladoceran Cercopagis pengoi, which invaded the Baltic Sea in the early 1990s. While several studies have reported immediate declines in abundances of its potential prey, the long-term effects of C. pengoi on the food webs remain to be examined. Based on the long-term time series (1968–2018) in the Gulf of Riga (Baltic Sea), we found significant declines in abundance of the cladoceran Pleopis spp. and copepod Eurytemora affinis by 90 and 80%, respectively, are associated with the invasion of C. pengoi as well as significant alterations in seasonal abundance patterns of Pleopis spp., E. affinis and cladoceran Bosmina spp. The invasion of the non-indigenous predator has led to the changed prey abundance–temperature relationships. Special caution was taken in data preprocessing, to minimize the likelihood that observed changes in the zooplankton prey could be associated with factors other than the invasion of C. pengoi, such as temperature and storminess.

Description: The quasi‐biennial oscillation (QBO) of the equatorial zonal wind leads to zonally symmetric temperature variations in the stratosphere that descend downward. Here we investigate the QBO‐induced temperature anomalies in the tropical tropopause layer (TTL) and detect pronounced longitudinal variations of the signal. In addition, the QBO temperature anomalies show a strong seasonal variability. The magnitude of these seasonal and longitudinal QBO variations is comparable to the magnitude of the well‐known zonal mean QBO signal in the TTL. At the cold point tropopause, the strongest QBO variations of around ±1.6 K are found over regions of active convection such as the West Pacific and Africa during boreal winter. The weakest QBO variations of ±0.25 K are detected over the East Pacific during boreal summer, while the zonal mean signal ranges around ±0.7 K. The longitudinal variations are associated with enhanced convective activity that occurs during QBO cold phases and locally enhances the cold anomalies.

Description: The last few decades have seen dramatic changes in the hydrography and biogeochemistry of the Mediterranean Sea. The complex bathymetry and highly variable spatial and temporal scales of atmospheric forcing, convective and ventilation processes contribute to generate complex and unsteady circulation patterns and significant variability in biogeochemical systems. Part of the variability of this system can be influenced by anthropogenic contributions. Consequently, it is necessary to document details and to understand trends in place to better relate the observed processes and to possibly predict the consequences of these changes. In this context we report data from an oceanographic cruise in the Mediterranean Sea on the German research vessel Maria S. Merian (MSM72) in March 2018. The main objective of the cruise was to contribute to the understanding of long-term changes and trends in physical and biogeochemical parameters, such as the anthropogenic carbon uptake and to further assess the hydrographical situation after the major climatological shifts in the eastern and western part of the basin, known as the Eastern and Western Mediterranean Transients. During the cruise, multidisciplinary measurements were conducted on a predominantly zonal section throughout the Mediterranean Sea, contributing to the Med-SHIP and GO-SHIP long-term repeat cruise section that is conducted at regular intervals in the Mediterranean Sea to observe changes and impacts on physical and biogeochemical variables.

Description: We present a mechanism for self‐sustained ocean circulation changes that cause abrupt temperature changes over Greenland in a multimillennial climate model simulation with glacial CO2 concentrations representative of Marine Isotope Stage 3. The Atlantic meridional overturning circulation (AMOC) and the subpolar gyre (SPG) oscillate on millennial time scales. When the AMOC is strong, the SPG is weak and contracted; when the AMOC is weak, the SPG is strong and extensive. The coupling between the two systems via wind‐driven and density‐driven feedbacks is key to maintaining the oscillations. The SPG controls the transport of heat and salt into the deep‐water formation sites and thus controls the AMOC strength. The strength and location of the deep‐water formation affect the density‐driven part of the SPG and thus control the mean strength and extent of the SPG. This mechanism supports the hypothesis that coupled ocean‐ice‐atmosphere interactions could have triggered abrupt glacial climate change.

Description: Decreasing concentrations of dissolved oxygen in the ocean are considered one of the main threats to marine ecosystems as they jeopardize the growth of higher organisms. They also alter the marine nitrogen cycle, which is strongly bound to the carbon cycle and climate. While higher organisms in general start to suffer from oxygen concentrations 〈 ∼ 63 µM (hypoxia), the marine nitrogen cycle responds to oxygen concentration below a threshold of about 20 µM (microbial hypoxia), whereas anoxic processes dominate the nitrogen cycle at oxygen concentrations of 〈 ∼ 0.05 µM (functional anoxia). The Arabian Sea and the Bay of Bengal are home to approximately 21 % of the total volume of ocean waters revealing microbial hypoxia. While in the Arabian Sea this oxygen minimum zone (OMZ) is also functionally anoxic, the Bay of Bengal OMZ seems to be on the verge of becoming so. Even though there are a few isolated reports on the occurrence of anoxia prior to 1960, anoxic events have so far not been reported from the open northern Indian Ocean (i.e., other than on shelves) during the last 60 years. Maintenance of functional anoxia in the Arabian Sea OMZ with oxygen concentrations ranging between 〉 0 and ∼ 0.05 µM is highly extraordinary considering that the monsoon reverses the surface ocean circulation twice a year and turns vast areas of the Arabian Sea from an oligotrophic oceanic desert into one of the most productive regions of the oceans within a few weeks. Thus, the comparably low variability of oxygen concentration in the OMZ implies stable balances between the physical oxygen supply and the biological oxygen consumption, which includes negative feedback mechanisms such as reducing oxygen consumption at decreasing oxygen concentrations (e.g., reduced respiration). Lower biological oxygen consumption is also assumed to be responsible for a less intense OMZ in the Bay of Bengal. According to numerical model results, a decreasing physical oxygen supply via the inflow of water masses from the south intensified the Arabian Sea OMZ during the last 6000 years, whereas a reduced oxygen supply via the inflow of Persian Gulf Water from the north intensifies the OMZ today in response to global warming. The first is supported by data derived from the sedimentary records, and the latter concurs with observations of decreasing oxygen concentrations and a spreading of functional anoxia during the last decades in the Arabian Sea. In the Arabian Sea decreasing oxygen concentrations seem to have initiated a regime shift within the pelagic ecosystem structure, and this trend is also seen in benthic ecosystems. Consequences for biogeochemical cycles are as yet unknown, which, in addition to the poor representation of mesoscale features in global Earth system models, reduces the reliability of estimates of the future OMZ development in the northern Indian Ocean.

Description: In the current era of rapid climate change, accurate characterization of climate-relevant gas dynamics-namely production, consumption, and net emissions-is required for all biomes, especially those ecosystems most susceptible to the impact of change. Marine environments include regions that act as net sources or sinks for numerous climateactive trace gases including methane (CH4) and nitrous oxide (N2O). The temporal and spatial distributions of CH4 and N2O are controlled by the interaction of complex biogeochemical and physical processes. To evaluate and quantify how these mechanisms affect marine CH4 and N2O cycling requires a combination of traditional scientific disciplines including oceanography, microbiology, and numerical modeling. Fundamental to these efforts is ensuring that the datasets produced by independent scientists are comparable and interoperable. Equally critical is transparent communication within the research community about the technical improvements required to increase our collective understanding of marine CH4 and N2O. A workshop sponsored by Ocean Carbon and Biogeochemistry (OCB) was organized to enhance dialogue and collaborations pertaining to marine CH4 and N2O. Here, we summarize the outcomes from the workshop to describe the challenges and opportunities for near-future CH4 and N2O research in the marine environment.

Description: It is widely accepted that orbital variations are responsible for the generation of glacial cycles during the late Pleistocene. However, the relative contributions of the orbital forcing compared to CO2 variations and other feedback mechanisms causing the waxing and waning of ice sheets have not been fully understood. Testing theories of ice ages beyond statistical inferences, requires numerical modeling experiments that capture key features of glacial transitions. Here, we focus on the glacial buildup from Marine Isotope Stage (MIS) 7 to 6 covering the period from 240 to 170 ka (ka: thousand years before present). This transition from interglacial to glacial conditions includes one of the fastest Pleistocene glaciation–deglaciation events, which occurred during MIS 7e–7d–7c (236–218 ka). Using a newly developed three-dimensional coupled atmosphere–ocean–vegetation–ice sheet model (LOVECLIP), we simulate the transient evolution of Northern Hemisphere and Southern Hemisphere ice sheets during the MIS 7–6 period in response to orbital and greenhouse gas forcing. For a range of model parameters, the simulations capture the evolution of global ice volume well within the range of reconstructions. Over the MIS 7–6 period, it is demonstrated that glacial inceptions are more sensitive to orbital variations, whereas terminations from deep glacial conditions need both orbital and greenhouse gas forcings to work in unison. For some parameter values, the coupled model also exhibits a critical North American ice sheet configuration, beyond which a stationary-wave–ice-sheet topography feedback can trigger an unabated and unrealistic ice sheet growth. The strong parameter sensitivity found in this study originates from the fact that delicate mass imbalances, as well as errors, are integrated during a transient simulation for thousands of years. This poses a general challenge for transient coupled climate–ice sheet modeling, with such coupled paleo-simulations providing opportunities to constrain such parameters.

Description: In 2015, we have collected more than 60,000 scavenging amphipod specimens during two expeditions to the Clarion-Clipperton fracture Zone (CCZ), in the Northeast (NE) Pacific and to the DISturbance and re-COLonisation (DisCOL) Experimental Area (DEA), a simulated mining impact disturbance proxy in the Peru basin, Southeast (SE) Pacific. Here, we compare biodiversity patterns of the larger specimens (〉15mm) within and between these two oceanic basins. Nine scavenging amphipod species are shared between these two areas, thus indicating connectivity. We further provide evidence that disturbance proxies seem to negatively affect scavenging amphipod biodiversity, as illustrated by a reduced alpha biodiversity in the DEA (Simpson Index (D)=0.62), when compared to the CCZ (D=0.73) and particularly of the disturbance site in the DEA and the site geographically closest to it. Community compositions of the two basins differs, as evidenced by a Non-Metric Dimensional Scaling (NMDS) analysis of beta biodiversity. The NMDS also shows a further separation of the disturbance site (D1) from its neighbouring, undisturbed reference areas (D2, D3, D4 and D5) in the DEA. A single species, Abyssorchomene gerulicorbis, dominates the DEA with 60% of all individuals.

Description: It is fundamentally important for many animal ecologists to quantify the costs of animal activities, although it is not straightforward to do so. The recording of triaxial acceleration by animal-attached devices has been proposed as a way forward for this, with the specific suggestion that dynamic body acceleration (DBA) be used as a proxy for movement-based power. Dynamic body acceleration has now been validated frequently, both in the laboratory and in the field, although the literature still shows that some aspects of DBA theory and practice are misunderstood. Here, we examine the theory behind DBA and employ modelling approaches to assess factors that affect the link between DBA and energy expenditure, from the deployment of the tag, through to the calibration of DBA with energy use in laboratory and field settings. Using data from a range of species and movement modes, we illustrate that vectorial and additive DBA metrics are proportional to each other. Either can be used as a proxy for energy and summed to estimate total energy expended over a given period, or divided by time to give a proxy for movement-related metabolic power. Nonetheless, we highlight how the ability of DBA to predict metabolic rate declines as the contribution of non-movement-related factors, such as heat production, increases. Overall, DBA seems to be a substantive proxy for movement-based power but consideration of other movement-related metrics, such as the static body acceleration and the rate of change of body pitch and roll, may enable researchers to refine movement-based metabolic costs, particularly in animals where movement is not characterized by marked changes in body acceleration.

Description: The speciation of dissolved iron (DFe) in the ocean is widely assumed to consist almost exclusively of Fe(III)-ligand complexes. Yet in most aqueous environments a poorly defined fraction of DFe also exists as Fe(II), the speciation of which is uncertain. Here we deploy flow injection analysis to measure in situ Fe(II) concentrations during a series of mesocosm/microcosm/multistressor experiments in coastal environments in addition to the decay rate of this Fe(II) when moved into the dark. During five mesocosm/microcosm/multistressor experiments in Svalbard and Patagonia, where dissolved (0.2 µm) Fe and Fe(II) were quantified simultaneously, Fe(II) constituted 24 %–65 % of DFe, suggesting that Fe(II) was a large fraction of the DFe pool. When this Fe(II) was allowed to decay in the dark, the vast majority of measured oxidation rate constants were less than calculated constants derived from ambient temperature, salinity, pH, and dissolved O2. The oxidation rates of Fe(II) spikes added to Atlantic seawater more closely matched calculated rate constants. The difference between observed and theoretical decay rates in Svalbard and Patagonia was most pronounced at Fe(II) concentrations 〈2 nM, suggesting that the effect may have arisen from organic Fe(II) ligands. This apparent enhancement of Fe(II) stability under post-bloom conditions and the existence of such a high fraction of DFe as Fe(II) challenge the assumption that DFe speciation in coastal seawater is dominated by ligand bound-Fe(III) species.

Description: A strong oxygen-deficient layer is located in the upper layers of the tropical Pacific Ocean and deeper in the North Pacific. Processes related to climate change (upper-ocean warming, reduced ventilation) are expected to change ocean oxygen and nutrient inventories. In most ocean basins, a decrease in oxygen (“deoxygenation”) and an increase in nutrients have been observed in subsurface layers. Deoxygenation trends are not linear and there could be multiple influences on oxygen and nutrient trends and variability. Here oxygen and nutrient time series since 1950 in the Pacific Ocean were investigated at 50 to 300 m depth, as this layer provides critical pelagic habitat for biological communities. In addition to trends related to ocean warming the oxygen and nutrient trends show a strong influence of the Pacific Decadal Oscillation (PDO) in the tropical and the eastern Pacific, and the North Pacific Gyre Oscillation (NPGO) in particular in the North Pacific. In the Oyashio Region the PDO, the NPGO, the North Pacific Index (NPI) and an 18.6-year nodal tidal cycle overlay the long-term trend. In most eastern Pacific regions oxygen increases and nutrients decrease in the 50 to 300 m layer during the negative PDO phase, with opposite trends during the positive PDO phase. The PDO index encapsulates the major mode of sea surface temperature variability in the Pacific, and oxygen and nutrients trends throughout the basin can be described in the context of the PDO phases. El Niño and La Niña years often influence the oxygen and nutrient distribution during the event in the eastern tropical Pacific but do not have a multi-year influence on the trends.

Description: Deciphering the dynamics of dissolved oxygen in the mid-depth ocean during the last deglaciation is essential to understand the influence of climate change on modern oxygen minimum zones (OMZs). Many paleo-proxy records from the Eastern Pacific Ocean indicate an extension of oxygen depleted conditions during the deglaciation but the degree of deoxygenation has not been quantified to date. The Peruvian OMZ, one of the largest OMZs in the world, is a key area to monitor such changes in near-bottom water oxygenation in relation to changing climatic conditions. Here, we analysed the potential to use the composition of foraminiferal assemblages from the Peruvian OMZ as a quantitative redox-proxy. A multiple regression analysis was applied to a joint dataset of living (rose Bengal stained, fossilizable calcareous species) benthic foraminiferal distributions from the Peruvian continental margin. Bottom-water oxygen concentrations ([O2]BW) during sampling were used as dependant variable. The correlation was significant (R2 = 0.82; p 〈 0.05) indicating that the foraminiferal assemblages are rather governed by oxygen availability than by the deposition of particulate organic matter (R2 = 0.53; p = 0.31). We applied the regression formula to four sediment cores from the northern part of the Peruvian OMZ between 3° S and 8° S and 600 m to 1250 m water depths; thereby recording oxygenation changes at the lower boundary of the Peruvian OMZ. Each core displayed a similar trend of decreasing oxygen levels since the Last Glacial Maximum (LGM). The overall [O2]BW change from the Last Glacial Maximum and the Holocene was constrained to 30 μmol/kg at the lower boundary of the OMZ, whereas at shallower depths [O2]BW was relatively stable along the deglaciation. The deoxygenation trend was time-transgressive. It commenced at the southern core, and gradually spread to deeper waters and to the northernmost core location. This pattern indicates a gradual expansion of the OMZ during the last deglaciation, as a result of increasing surface productivity in the Eastern Equatorial Pacific and decreasing advective oxygen supply to intermediate waters off Peru.

Description: Oxygen minimum zones (OMZs) show distinct biogeochemical processes that relate to microorganisms being able to thrive under low or even absent oxygen. Microbial degradation of organic matter is expected to be reduced in OMZs, although quantitative evidence is low. Here, we present heterotrophic bacterial production (3H leucine incorporation), extracellular enzyme rates (leucine aminopeptidase/β-glucosidase) and bacterial cell abundance for various in situ oxygen concentrations in the water column, including the upper and lower oxycline, of the eastern tropical South Pacific off Peru. Bacterial heterotrophic activity in the suboxic core of the OMZ (at in situ ≤ 5 µmol O2 kg−1) ranged from 0.3 to 281 µmol C m−3 d−1 and was not significantly lower than in waters of 5–60 µmol O2 kg−1. Moreover, bacterial abundance in the OMZ and leucine aminopeptidase activity were significantly higher in suboxic waters compared to waters of 5–60 µmol O2 kg−1, suggesting no impairment of bacterial organic-matter degradation in the core of the OMZ. Nevertheless, high cell-specific bacterial production was observed in samples from oxyclines, and cell-specific extracellular enzyme rates were especially high at the lower oxycline, corroborating earlier findings of highly active and distinct micro-aerobic bacterial communities. To assess the impact of bacterial degradation of dissolved organic matter (DOM) for oxygen loss in the Peruvian OMZ, we compared diapycnal fluxes of oxygen and dissolved organic carbon (DOC) and their microbial uptake within the upper 60 m of the water column. Our data indicate low bacterial growth efficiencies of 1 %–21 % at the upper oxycline, resulting in a high bacterial oxygen demand that can explain up to 33 % of the observed average oxygen loss over depth. Our study therewith shows that microbial degradation of DOM has a considerable share in sustaining the OMZ off Peru.

Description: Small steam-driven volcanic explosions are common at volcanoes worldwide but are rarely documented or monitored; therefore, these events still put residents and tourists at risk every year. Steam-driven explosions also occur frequently (once every 2–5 years on average) at Lascar volcano, Chile, where they are often spontaneous and lack any identifiable precursor activity. Here, for the first time at Lascar, we describe the processes culminating in such a sudden volcanic explosion that occurred on October 30, 2015, which was thoroughly monitored by cameras, a seismic network, and gas (SO2 and CO2) and temperature sensors. Prior to the eruption, we retrospectively identified unrest manifesting as a gradual increase in the number of long-period (LP) seismic events in 2014, indicating an augmented level of activity at the volcano. Additionally, SO2 flux and thermal anomalies were detected before the eruption. Then, our weather station reported a precipitation event, followed by changes in the brightness of the permanent volcanic plume and (10 days later) by the sudden volcanic explosion. The multidisciplinary data exhibited short-term variations associated with the explosion, including (1) an abrupt eruption onset that was seismically identified in the 1–10 Hz frequency band, (2) the detection of a 1.7 km high white-grey eruption column in camera images, and (3) a pronounced spike in sulfur dioxide (SO2) emission rates reaching 55 kg sec−1 during the main pulse of the eruption as measured by a mini-DOAS scanner. Continuous CO2 gas and temperature measurements conducted at a fumarole on the southern rim of the Lascar crater revealed a pronounced change in the trend of the relationship between the carbon dioxide (CO2) mixing ratio and the gas outlet temperature; we believe that this change was associated with the prior precipitation event. An increased thermal anomaly inside the active crater observed through Sentinel-2 images and drone overflights performed after the steam-driven explosion revealed the presence of a fracture ~ 50 metres in diameter truncating the dome and located deep inside the active crater, which coincides well with the location of the thermal anomaly. Altogether, these observations lead us to infer that a lava dome was present and subjected to cooling and inhibited degassing. We conjecture that a precipitation event led to the short-term build-up of pressure inside the shallow dome that eventually triggered a vent-clearing phreatic explosion. This study shows the chronology of events culminating in a steam-driven explosion but also demonstrates that phreatic explosions are difficult to forecast, even if the volcano is thoroughly monitored; these findings also emphasize why ascending to the summits of Lascar and similar volcanoes is hazardous, particularly after considerable rainfall.

Description: The intraseasonal evolution of physical and biogeochemical properties during a coastal trapped wave event off central Peru is analysed using data from an extensive shipboard observational programme conducted between April and June 2017, and remote sensing data. The poleward velocities in the Peru–Chile Undercurrent were highly variable and strongly intensified to above 0.5 m s−1 between the middle and end of May. This intensification was likely caused by a first-baroclinic-mode downwelling coastal trapped wave, excited by a westerly wind anomaly at the Equator and originating at about 95∘ W. Local winds along the South American coast did not impact the wave. Although there is general agreement between the observed cross-shore-depth velocity structure of the coastal trapped wave and the velocity structure of first vertical mode solution of a linear wave model, there are differences in the details of the two flow distributions. The enhanced poleward flow increased water mass advection from the equatorial current system to the study site. The resulting shorter alongshore transit times between the Equator and the coast off central Peru led to a strong increase in nitrate concentrations, less anoxic water, likely less fixed nitrogen loss to N2 and a decrease of the nitrogen deficit compared to the situation before the poleward flow intensification. This study highlights the role of changes in the alongshore advection due to coastal trapped waves for the nutrient budget and the cumulative strength of N cycling in the Peruvian oxygen minimum zone. Enhanced availability of nitrate may impact a range of pelagic and benthic elemental cycles, as it represents a major electron acceptor for organic carbon degradation during denitrification and is involved in sulfide oxidation in sediments.

Description: We conduct the seismic signal analysis on vintage and recently collected multichannel seismic reflection profiles from the Ionian Basin to characterize the deep basin Messinian evaporites. These evaporites were deposited in deep and marginal Mediterranean sedimentary basins as a consequence of the “salinity crisis” between 5.97 and 5.33 Ma, a basin‐wide oceanographic and ecological crisis whose origin remains poorly understood. The seismic markers of the Messinian evaporites in the deep Mediterranean basins can be divided in two end‐members, one of which is the typical “trilogy” of gypsum and clastics (Lower Unit – LU), halite (Mobile Unit – MU) and upper anhydrite and marl layers (Upper Unit – UU) traced in the Western Mediterranean Basins. The other end‐member is a single MU unit subdivided in seven sub‐units by clastic interlayers located in the Levant Basin. The causes of these different seismic expressions of the Messinian salinity crisis (MSC) appear to be related to a morphological separation between the two basins by the structural regional sill of the Sicily Channel. With the aid of velocity analyses and seismic imaging via prestack migration in time and depth domains, we define for the first time the seismic signature of the Messinian evaporites in the deep Ionian Basin, which differs from the known end‐members. In addition, we identify different evaporitic depositional settings suggesting a laterally discontinuous deposition. With the information gathered we quantify the volume of evaporitic deposits in the deep Ionian Basin as 500,000 km3 ± 10%. This figure allows us to speculate that the total volume of salts in the Mediterranean basin is larger than commonly assumed. Different depositional units in the Ionian Basin suggest that during the MSC it was separated from the Western Mediterranean by physical thresholds, from the Po Plain/Northern Adriatic Basin, and the Levant Basin, likely reflecting different hydrological and climatic conditions. Finally, the evidence of erosional surfaces and V‐shaped valleys at the top of the MSC unit, together with sharp evaporites pinch out on evaporite‐free pre‐Messinian structural highs, suggest an extreme Messinian Stage 3 base level draw down in the Ionian Basin. Such evidence should be carefully evaluated in the light of Messinian and post‐Messinian vertical crustal movements in the area. The results of this study demonstrates the importance of extracting from seismic data the Messinian paleotopography, the paleomorphology and the detailed stratal architecture in the in order to advance in the understanding of the deep basins Messinian depositional environments. Highlights First description of a new type of deepwater Messinian salt giant in the Ionian Sea. First quantification of the Messinian salt volume in the Ionian Sea. New seismic evidence of erosional surfces and Lago Mare deposits in the deep Ionian Basin. Further evidence of sea level lowering during the Messinian Salinity Crisis. Evidence for a different, physically separated deepwater Messinian salt basins in the Mediterranean.

Description: Carbonyl sulfide (OCS) and carbon disulfide (CS2) are volatile sulfur gases that are naturally formed in seawater and exchanged with the atmosphere. OCS is the most abundant sulfur gas in the atmosphere, and CS2 is its most important precursor. They have gained interest due to their direct (OCS) or indirect (CS2 via oxidation to OCS) contribution to the stratospheric sulfate aerosol layer. Furthermore, OCS serves as a proxy to constrain terrestrial CO2 uptake by vegetation. Oceanic emissions of both gases contribute a major part to their atmospheric concentration. Here we present a database of previously published and unpublished, mainly ship-borne measurements in seawater and the marine boundary layer for both gases, available at https://doi.pangaea.de/10.1594/PANGAEA.905430 (Lennartz et al., 2019). The database contains original measurements as well as data digitalized from figures in publications from 42 measurement campaigns, i.e. cruises or time series stations, ranging from 1982 to 2019. OCS data cover all ocean basins except for the Arctic Ocean, as well as all months of the year, while the CS2 dataset shows large gaps in spatial and temporal coverage. Concentrations are consistent across different sampling and analysis techniques for OCS. The database is intended to support the identification of global spatial and temporal patterns and to facilitate the evaluation of model simulations.

Description: Prochlorococcus and Synechococcus are the most abundant free-living photosynthetic microorganisms in the ocean. Uncultivated lineages of these picocyanobacteria also thrive in the dimly illuminated upper part of oxygen-deficient zones (ODZs), where an important portion of ocean nitrogen (N) loss takes place via denitrification and anaerobic ammonium oxidation. Recent metagenomic studies revealed that ODZ Prochlorococcus have the genetic potential for using different N forms, including nitrate and nitrite, uncommon N sources for Prochlorococcus, but common for Synechococcus. To determine which N sources ODZ picocyanobacteria are actually using in nature, the cellular N-15 natural abundance (delta N-15) and assimilation rates of different N compounds were determined using cell sorting by flow cytometry and mass spectrometry. The natural delta N-15 of the ODZ Prochlorococcus varied from -4.0 parts per thousand to 13.0 parts per thousand (n = 9), with 50% of the values in the range of -2.1-2.6 parts per thousand. While the highest values suggest nitrate use, most observations indicate the use of nitrite, ammonium, or a mixture of N sources. Meanwhile, incubation experiments revealed potential assimilation rates of ammonium and urea in the same order of magnitude as that expected for total N in several environments including ODZs, whereas rates of nitrite and nitrate assimilation were very low. Our results thus indicate that reduced forms of N and nitrite are the dominant sources for ODZ picocyanobacteria, although nitrate might be important on some occasions. ODZ picocyanobacteria might thus represent potential competitors with anammox bacteria for ammonium and nitrite, with ammonia-oxidizing archaea for ammonium, and with nitrite-oxidizing bacteria for nitrite.

Description: The immune system affects senescence (declines in probabilities of survival or reproduction with age), by shaping late age vulnerability to chronic inflammatory diseases and infections. It is also a dynamic interactive system that must balance competing demands across the life course. Thus, immune system function remains an important frontier in understanding the evolution of senescence. Here, we review our expanding mechanistic understanding of immune function over the life course, in the context of theoretical predictions from life-history evolution. We are especially interested in stage- and sex-dependent costs and benefits of investment in the immune system, given differential life-history priorities of the life stages and sexes. We introduce the costs likely to govern immune allocation across the life course. We then discuss theoretical expectations for differences between the sexes and their likely consequences in terms of how the immune system is both modulated by and may modulate senescence, building on information from life-history theory, experimental immunology and demography. We argue that sex differences in immune function provide a potentially powerful probe of selection pressures on the immune system across the life course. In particular, differences in 'competing' and 'caring' between the sexes have evolved across the tree of life, providing repeated instances of divergent selection pressures on immune function occurring within the same overall bauplan. We conclude by detailing an agenda for future research, including development of theoretical predictions of the differences between the sexes under an array of existing models for sex differences in immunity, and empirical tests of such predictions across the tree of life. A free Plain Language Summary can be found within the Supporting Information of this article.

Description: The plea for using more “realistic,” community‐level, investigations to assess the ecological impacts of global change has recently intensified. Such experiments are typically more complex, longer, more expensive, and harder to interpret than simple organism‐level benchtop experiments. Are they worth the extra effort? Using outdoor mesocosms, we investigated the effects of ocean warming (OW) and acidification (OA), their combination (OAW), and their natural fluctuations on coastal communities of the western Baltic Sea during all four seasons. These communities are dominated by the perennial and canopy‐forming macrophyte Fucus vesiculosus—an important ecosystem engineer Baltic‐wide. We, additionally, assessed the direct response of organisms to temperature and pH in benchtop experiments, and examined how well organism‐level responses can predict community‐level responses to the dominant driver, OW. OW affected the mesocosm communities substantially stronger than acidification. OW provoked structural and functional shifts in the community that differed in strength and direction among seasons. The organism‐level response to OW matched well the community‐level response of a given species only under warm and cold thermal stress, that is, in summer and winter. In other seasons, shifts in biotic interactions masked the direct OW effects. The combination of direct OW effects and OW‐driven shifts of biotic interactions is likely to jeopardize the future of the habitat‐forming macroalga F. vesiculosus in the Baltic Sea. Furthermore, we conclude that seasonal mesocosm experiments are essential for our understanding of global change impact because they take into account the important fluctuations of abiotic and biotic pressures.

Description: Over the past decades, coralline algae have increasingly been used as archives of palaeoclimate information due to their seasonal growth bands and their vast distribution from high latitudes to the tropics. Traditionally, these reconstructions have been performed mainly on high latitude species, limiting the geographical area of their potential use. Here we assess the use of temperate crustose fossil coralline algae from shallow water habitats for palaeoenvironmental reconstruction to generate records of past climate change. We determine the potential of three different species of coralline algae, Lithothamnion minervae, Lithophyllum stictaeforme and Mesophyllum philippii, with different growth patterns, as archives for pH (δ11B) and temperature (Mg/Ca) reconstruction in the Mediterranean Sea. Mg concentration is driven by temperature but modulated by growth rate, which is controlled by species‐specific and intraspecific growth patterns. L. minervae is a good temperature recorder, showing a moderate warming trend in specimens from 11.37 cal ka BP (from 14.2 ± 0.4°C to 14.9 ± 0.15°C) to today. In contrast to Mg, all genera showed consistent values of boron isotopes (δ11B) suggesting a common control on boron incorporation. The recorded δ11B in modern and fossil coralline specimens is in agreement with literature data about early Holocene pH, opening new perspectives of coralline‐based, high‐resolution pH reconstructions in deep time.

Description: Tephra layers produced by volcanic eruptions are widely used for correlation and dating of various deposits and landforms, for synchronization of disparate paleoenvironmental archives, and for reconstruction of magma origin. Here we present our original database TephraKam, which includes chemical compositions of volcanic glass in tephra and welded tuffs from the Kamchatka volcanic arc. The database contains 7049 major element analyses obtained by electron microprobe and 738 trace element analyses obtained by laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) on 487 samples collected in proximity of their volcanic sources in all volcanic zones in Kamchatka. The samples characterize about 300 explosive eruptions, which occurred in Kamchatka from the Pliocene until historic times. Precise or estimated ages for all samples are based on published 39Ar/40Ar dates of rocks and 14C dates of host sediments, statistical age modelling and geologic relationships with dated units. All data in TephraKam is supported by information about source volcanoes and analytical details. Using the data, we present an overview of geochemical variations of Kamchatka volcanic glasses and discuss application of this data for precise identification of tephra layers, their source volcanoes, temporal and spatial geochemical variations of pyroclastic rocks in Kamchatka.

Description: High-resolution optical and hydroacoustic seafloor data acquired in 2015 enabled the reconstruction of disturbance tracks of a past Benthic Impact Experiment that was conducted in 1989 in the Peru Basin in the course of former German environmental impact studies associated with manganese nodule mining. Based on this information, the disturbance level of the experiment regarding the plough impact and distribution and re-deposition of sediment from the evolving sediment plume was assessed qualitatively. Through this, the evolution over the 26 years of a number of the total 78 disturbance tracks could be analyzed which highlights the considerable difference between natural sedimentation in the deep-sea and sedimentation of a resettled sediment plume. Such plumes are seen as one of the most concerning impact associated with potential Mn-nodule mining. Problems in data processing became eminent while dealing with old data from the late 80s, at a time when GPS was just invented and underwater navigation was in an infant stage. However, even today the uncertainties of underwater navigation and the use of a variety of acoustical and optical sensors at different resolutions require detailed post-processing in terms of absolute geographic positioning to improve the overall accuracy of the data. In this study, a ship-based bathymetric map of the survey area was used as absolute geographic reference and a workflow was applied successfully resulting in the most accurate geo-referenced dataset of the DISCOL Experimental Area to date. The new field data were acquired with sensors attached to GEOMARs AUV Abyss and the 0.5 × 1° EM122 multibeam system of RV SONNE during cruise SO242 -1 while the old data first needed to be found and compiled before they could be digitized and properly georeferenced for the presented joined analyses.

Description: We report on the geochemistry of hydrocarbons and pore waters down to 62.5 mbsf, collected by drilling with the MARUM‐MeBo70 and by gravity coring at the Lunde pockmark in the Vestnesa Ridge. Our data document the origin and transformations of volatiles feeding gas emissions previously documented in this region. Gas hydrates are present where a fracture network beneath the pockmark focusses migration of thermogenic hydrocarbons characterized by their C1/C2+ and stable isotopic compositions (δ2H‐CH4, δ13C‐CH4). Measured geothermal gradients (~80°C km‐1) and known formation temperatures (〉70°C) suggest that those hydrocarbons are formed at depths 〉800 mbsf. A combined analytical/modeling approach, including concentration and isotopic mass balances, reveals that pockmark sediments experience diffuse migration of thermogenic hydrocarbons. However, at sites without channeled flow this appears to be limited to depths 〉 ~50 mbsf. At all sites we document a contribution of microbial methanogenesis to the overall carbon cycle that includes a component of secondary carbonate reduction (CR) – i.e. reduction of dissolved inorganic carbon (DIC) generated by anaerobic oxidation of methane (AOM) in the uppermost methanogenic zone. AOM and CR rates are spatially variable within the pockmark and are highest at high‐flux sites. These reactions are revealed by δ13C‐DIC depletions at the sulfate‐methane interface at all sites. However, δ13C‐CH4 depletions are only observed at the low methane flux sites because changes in the isotopic composition of the overall methane pool are masked at high‐flux sites. 13C‐depletions of TOC suggest that at seeps sites, methane‐derived carbon is incorporated into de novo synthesized biomass.

Description: Sympatric speciation occurs without geographical barriers and is thought to often be driven by ecological specialization of individuals that eventually diverge genetically and phenotypically. Distinct morphologies between sympatric populations occupying different niches have been interpreted as such differentiating adaptive phenotypes, yet differences in performance and thus likely adaptiveness between them were rarely tested. Here, we investigated if divergent body shapes of two sympatric crater lake cichlid species from Nicaragua, one being a shore‐associated (benthic) species while the other prefers the open water zones (limnetic), affect cruising (Ucrit) and sprinting (Usprint) swimming abilities ‐ performances particularly relevant to their respective lifestyles. Furthermore, we investigated species differences in oxygen consumption (MO2) across different swimming speeds and compare gene expression in gills and white muscle at rest and during exercise. We found a superior cruising ability in the limnetic Amphilophus zaliosus compared to the benthic A. astorquii, while sprinting was not different, suggesting that their distinct morphologies affect swimming performance. Increased cruising swimming ability in A. zaliosus was linked to a higher oxygen demand during activity (but not rest), indicating different metabolic rates during exercise ‐ a hypothesis supported by coinciding gene expression patterns of gill transcriptomes. We identified differentially expressed genes linked to swimming physiology, regulation of swimming behaviour and oxygen intake. A combination of physiological and morphological differences may thus underlie adaptations to these species’ distinct niches. This complex ecological specialization likely resulted in morphological and physiological trade‐offs that contributed to the rapid establishment and maintenance of divergence with gene flow.

Description: The duration and magnitude of the North Atlantic spring bloom impacts both higher trophic levels and oceanic carbon sequestration. Nutrient exhaustion offers a general explanation for bloom termination, but detail on which nutrients and their relative influence on phytoplankton productivity, community structure, and physiology is lacking. Here, we address this using nutrient addition bioassay experiments conducted across the midlatitude North Atlantic in June 2017 (late spring). In four out of six experiments, phytoplankton accumulated over 48–72 h following individual additions of either iron (Fe) or nitrogen (N). In the remaining two experiments, Fe and N were serially limiting, that is, their combined addition sequentially enhanced phytoplankton accumulation. Silicic acid (Si) added in combination with N + Fe led to further chlorophyll a (Chl a) enhancement at two sites. Conversely, addition of zinc, manganese, cobalt, vitamin B12, or phosphate in combination with N + Fe did not. At two sites, the simultaneous supply of all six nutrients, in combination with N + Fe, also led to no further Chl a enhancement, but did result in an additional 30–60% particulate carbon accumulation. This particulate carbon accumulation was not matched by a Redfield equivalent of particulate N, characteristic of high C:N organic exudates that enhance cell aggregation and sinking. Our results suggest that growth rates of larger phytoplankton were primarily limited by Fe and/or N, making the availability of these nutrients the main bottom‐up factors contributing to spring bloom termination. In addition, the simultaneous availability of other nutrients could modify bloom characteristics and carbon export efficiency.

Description: Oxygen-deficient zones (ODZs) are major sites of net natural nitrous oxide (N2O) production and emissions. In order to understand changes in the magnitude of N2O production in response to global change, knowledge on the individual contributions of the major microbial pathways (nitrification and denitrification) to N2O production and their regulation is needed. In the ODZ in the coastal area off Peru, the sensitivity of N2O production to oxygen and organic matter was investigated using 15N tracer experiments in combination with quantitative PCR (qPCR) and microarray analysis of total and active functional genes targeting archaeal amoA and nirS as marker genes for nitrification and denitrification, respectively. Denitrification was responsible for the highest N2O production with a mean of 8.7 nmol L−1 d−1 but up to 118±27.8 nmol L−1 d−1 just below the oxic–anoxic interface. The highest N2O production from ammonium oxidation (AO) of 0.16±0.003 nmol L−1 d−1 occurred in the upper oxycline at O2 concentrations of 10–30 µmol L−1 which coincided with the highest archaeal amoA transcripts/genes. Hybrid N2O formation (i.e., N2O with one N atom from NH+4 and the other from other substrates such as NO−2) was the dominant species, comprising 70 %–85 % of total produced N2O from NH+4, regardless of the ammonium oxidation rate or O2 concentrations. Oxygen responses of N2O production varied with substrate, but production and yields were generally highest below 10 µmol L−1 O2. Particulate organic matter additions increased N2O production by denitrification up to 5-fold, suggesting increased N2O production during times of high particulate organic matter export. High N2O yields of 2.1 % from AO were measured, but the overall contribution by AO to N2O production was still an order of magnitude lower than that of denitrification. Hence, these findings show that denitrification is the most important N2O production process in low-oxygen conditions fueled by organic carbon supply, which implies a positive feedback of the total oceanic N2O sources in response to increasing oceanic deoxygenation.

Description: Considerable uncertainty remains over how increasing atmospheric CO2 and anthropogenic climate changes are affecting open‐ocean marine ecosystems from phytoplankton to top predators. Biological time series data are thus urgently needed for the world's oceans. Here, we use the carbon stable isotope composition of tuna to provide a first insight into the existence of global trends in complex ecosystem dynamics and changes in the oceanic carbon cycle. From 2000 to 2015, considerable declines in δ13C values of 0.8‰–2.5‰ were observed across three tuna species sampled globally, with more substantial changes in the Pacific Ocean compared to the Atlantic and Indian Oceans. Tuna recorded not only the Suess effect, that is, fossil fuel‐derived and isotopically light carbon being incorporated into marine ecosystems, but also recorded profound changes at the base of marine food webs. We suggest a global shift in phytoplankton community structure, for example, a reduction in 13C‐rich phytoplankton such as diatoms, and/or a change in phytoplankton physiology during this period, although this does not rule out other concomitant changes at higher levels in the food webs. Our study establishes tuna δ13C values as a candidate essential ocean variable to assess complex ecosystem responses to climate change at regional to global scales and over decadal timescales. Finally, this time series will be invaluable in calibrating and validating global earth system models to project changes in marine biota.

Description: Key Points Calibration of XRF core scanning data highlights the need for careful examination of sediment properties such as porosity/water Grain size and water content in the sediment trigger systematic artifacts in the signal intensity of light elements (e.g. Si and Al) Known terrigenous flux proxies (e.g Ti/Ca, Fe/Ca) are influenced by sea level variations X‐ray fluorescence (XRF) core scanning of marine and lake sediments has been extensively used to study changes in past environmental and climatic processes over a range of timescales. The interpretation of XRF‐derived element ratios in paleoclimatic and paleoceanographic studies primarily considers differences in the relative abundances of particular elements. Here we present new XRF core scanning data from two long sediment cores in the Andaman Sea in the northern Indian Ocean and show that sea level related processes influence terrigenous inputs based proxies such as Ti/Ca, Fe/Ca, and elemental concentrations of the transition metals (e.g. Mn). Zr/Rb ratios are mainly a function of changes in median grain size of lithogenic particles and often covary with changes in Ca concentrations that reflect changes in biogenic calcium carbonate production. This suggests that a common process (i.e. sea level) influences both records. The interpretation of lighter element data (e.g. Si and Al) based on low XRF counts is complicated as variations in mean grain size and water content result in systematic artifacts and signal intensities not related to the Al or Si content of the sediments. This highlights the need for calibration of XRF core scanning data based on discrete sample analyses and careful examination of sediment properties such as porosity/water content for reliably disentangling environmental signals from other physical properties. In the case of the Andaman Sea, reliable extraction of a monsoon signal will require accounting for the sea level influence on the XRF data.

Description: Seamounts represent ideal systems to study the influence and interdependency of environmental gradients at a single geographic location. These topographic features represent a prominent habitat for various forms of life, including microbiota and macrobiota, spanning benthic as well as pelagic organisms. While it is known that seamounts are globally abundant structures, it still remains unclear how and to which extent the complexity of the sea floor is intertwined with the local oceanographic mosaic, biogeochemistry, and microbiology of a seamount ecosystem. Along these lines, the present study aimed to explore whether and to what extent seamounts can have an imprint on the microbial community composition of seawater and of sessile benthic invertebrates, sponges. For our high-resolution sampling approach of microbial diversity (16S rRNA gene amplicon sequencing) along with measurements of inorganic nutrients and other biogeochemical parameters, we focused on the Schulz Bank seamount ecosystem, a sponge ground ecosystem which is located on the Arctic Mid-Ocean Ridge. Seawater samples were collected at two sampling depths (mid-water, MW, and near-bed water, BW) from a total of 19 sampling sites. With a clustering approach we defined microbial microhabitats within the pelagic realm at Schulz Bank, which were mapped onto the seamount's topography and related to various environmental parameters (such as suspended particulate matter, SPM; dissolved inorganic carbon, DIC; silicate, SiO−4; phosphate, PO3−4; ammonia, NH+4; nitrate, NO2−3; nitrite, NO−2; depth; and dissolved oxygen, O2). The results of our study reveal a “seamount effect” (sensu stricto) on the microbial mid-water pelagic community at least 200 m above the sea floor. Further, we observed a strong spatial heterogeneity in the pelagic microbial landscape across the seamount, with planktonic microbial communities reflecting oscillatory and circulatory water movements, as well as processes of bentho-pelagic coupling. Depth, NO2−3, SiO−4, and O2 concentrations differed significantly between the determined pelagic microbial clusters close to the sea floor (BW), suggesting that these parameters were presumably linked to changes in microbial community structures. Secondly, we assessed the associated microbial community compositions of three sponge species along a depth gradient of the seamount. While sponge-associated microbial communities were found to be mainly species-specific, we also detected significant intra-specific differences between individuals, depending on the pelagic near-bed cluster they originated from. The variable microbial phyla (i.e. phyla which showed significant differences across varying depth, NO2−3, SiO−4, O2 concentrations, and different from local seawater communities) were distinct for every sponge species when considering average abundances per species. Variable microbial phyla included representatives of both those taxa traditionally counted for the variable community fraction and taxa counted traditionally for the core community fraction. Microbial co-occurrence patterns for the three examined sponge species Geodia hentscheli, Lissodendoryx complicata, and Schaudinnia rosea were distinct from each other. Over all, this study shows that topographic structures such as the Schulz Bank seamount can have an imprint (seamount effect sensu lato) on both the microbial community composition of seawater and sessile benthic invertebrates such as sponges by an interplay between the geology, physical oceanography, biogeochemistry, and microbiology of seamounts.

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: The copepod Acartia tonsa is a key component of a wide range of marine ecosystems, linking energy transfer from phytoplankton to higher trophic levels, and has a central role in productivity and biogeochemistry. The interaction of end-of-century global warming and ocean acidification scenarios with testing moderate temperature effects on a seminatural copepod community is needed to understand future community functioning. Here, we deployed a mesocosm experimental set-up with a full factorial design using two temperatures (13°C and 19°C) crossed with a pCO2 gradient ranging from ambient (550 μatm) to 3000 μatm. We used the natural bacteria, phyto- and microzooplankton species composition and biomass of the Kiel Bight and tested the response of A. tonsa development, carbon growth, mortality, size and condition. The tested traits were differently affected by the interaction of temperature and acidification. Ocean acidification increased development, carbon growth, size and mortality under the warming scenario of 19°C. At 13°C mortality rates decreased, while carbon growth, size and condition increased with acidification. We conclude from our experimental approach that a single species shows a variety of responses depending on the focal functional trait. Trait-specific mesozooplankton responses need to be further investigated and compared between geographical regions, seasons and taxonomic groups.

Description: Along the margins of continental ice sheets, lakes formed in isostatically depressed basins duringglacial retreat. Their shorelines and extent are sensitive to the ice margin and the glacial history of the region.Proglacial lakes, in turn, also impact the glacial isostatic adjustment due to loading, and ice dynamics by posing amarine‐like boundary condition at the ice margin. In this study we present a tool that efficiently identifies lake basinsand the corresponding maximum water level for a given ice sheet and topography reconstruction. This algorithm,called the LakeCC model, iteratively checks the whole map for a set of increasing water levels and fills isolated basinsuntil they overflow into the ocean. We apply it to the present‐day Great Lakes and the results show good agreement(∼1−4%) with measured lake volume and depth. We then apply it to two topography reconstructions of NorthAmerica between the Last Glacial Maximum and the present. The model successfully reconstructs glacial lakes suchas Lake Agassiz, Lake McConnell and the predecessors of the Great Lakes. LakeCC can be used to judge the quality ofice sheet reconstructions.

Description: The largest and commercially appealing mineral deposits can be found in the abyssal seafloor of the Clarion-Clipperton Zone (CCZ), a polymetallic nodule province, in the NE Pacific Ocean, where experimental mining is due to take place. In anticipation of deep-sea mining impacts, it has become essential to rapidly and accurately assess biodiversity. For this reason, ophiuroid material collected during seven scientific cruises from five exploration license areas within CCZ, one area protected from mining (APEI3, Area of Particular Environmental Interest) in the periphery of CCZ and the DIS-turbance and re-COLonisation (DISCOL) Experimental Area (DEA), in the SE Pacific Ocean, was examined. Specimens were genetically analysed using a fragment of the mitochondrial cytochrome c oxidase subunit I (COI). Maximum Likelihood and Neighbour Joining trees were constructed, while four tree-based and distance-based methods of species delineation (ABGD, BINs, GMYC, mPTP) were employed to propose Secondary Species Hypotheses (SSHs) within the ophiuroids collected. The species delimitations analyses concordant results revealed the presence of 43 deep-sea brittle stars SSHs, revealing an unexpectedly high diversity and showing that the most conspicuous invertebrates in abyssal plains have been so far considerably under-estimated. The number of SSHs found in each area varied from 5 (IFREMER area) to 24 (BGR area), while 13 SSHs were represented by singletons. None of the SSHs was found to be present in all 7 areas, while the majority of species (44.2 %) had a single-area presence (19 SSHs). The most common species were Ophioleucidae sp. (Species 29), Amphioplus daleus (Species 2) and Ophiosphalma glabrum (Species 3), present in all areas except APEI3. The biodiversity patterns could be mainly attributed to POC fluxes that could explain the highest species numbers found in BGR (German contractor area) and UKSRL (UK contractor area) areas. The five exploration contract areas belong to a mesotrophic province, while in contrary the APEI3 is located in an oligotrophic province which could explain the lowest diversity as well as very low similarity with the other six study areas. Based on these results the representativeness and the appropriateness of APEI3 to meet its purpose of preserving the biodiversity of the CCZ fauna are questioned. Finally, this study provides the foundation for biogeographic and functional analyses that will provide insight into the drivers of species diversity and its role in ecosystem function.

Description: Dissolved Fe (DFe) samples from the GEOVIDE voyage (GEOTRACES GA01, May–June 2014) in the North Atlantic Ocean were analyzed using a seaFAST-pico™ coupled to an Element XR sector field inductively coupled plasma mass spectrometer (SF-ICP-MS) and provided interesting insights into the Fe sources in this area. Overall, DFe concentrations ranged from 0.09±0.01 to 7.8±0.5 nmol L−1. Elevated DFe concentrations were observed above the Iberian, Greenland, and Newfoundland margins likely due to riverine inputs from the Tagus River, meteoric water inputs, and sedimentary inputs. Deep winter convection occurring the previous winter provided iron-to-nitrate ratios sufficient to sustain phytoplankton growth and lead to relatively elevated DFe concentrations within subsurface waters of the Irminger Sea. Increasing DFe concentrations along the flow path of the Labrador Sea Water were attributed to sedimentary inputs from the Newfoundland Margin. Bottom waters from the Irminger Sea displayed high DFe concentrations likely due to the dissolution of Fe-rich particles in the Denmark Strait Overflow Water and the Polar Intermediate Water. Finally, the nepheloid layers located in the different basins and at the Iberian Margin were found to act as either a source or a sink of DFe depending on the nature of particles, with organic particles likely releasing DFe and Mn particle scavenging DFe.

Description: In the abyssal equatorial Pacific Ocean, most of the seafloor of the Clarion-Clipperton Fracture Zone (CCFZ), a 6 million km2 polymetallic nodule province, has been preempted for future mining. In light of the large environmental footprint that mining would leave and given the diversity and the vulnerability of the abyssal fauna, the International Seabed Authority has implemented a regional management plan that includes the creation of nine Areas of Particular Environmental Interest (APEIs) located at the periphery of the CCFZ. The scientific principles for the design of the APEIs were based on the best – albeit very limited – knowledge of the area. The fauna and habitats in the APEIs are unknown, as are species' ranges and the extent of biodiversity across the CCFZ. As part of the Joint Programming Initiative Healthy and Productive Seas and Oceans (JPI Oceans) pilot action “Ecological aspects of deep-sea mining”, the SO239 cruise provided data to improve species inventories, determine species ranges, identify the drivers of beta diversity patterns and assess the representativeness of an APEI. Four exploration contract areas and an APEI (APEI no. 3) were sampled along a gradient of sea surface primary productivity that spanned a distance of 1440 km in the eastern CCFZ. Between three and eight quantitative box cores (0.25 m2; 0–10 cm) were sampled in each study area, resulting in a large collection of polychaetes that were morphologically and molecularly (cytochrome c oxidase subunit I and 16S genes) analyzed. A total of 275 polychaete morphospecies were identified. Only one morphospecies was shared among all five study areas and 49 % were singletons. The patterns in community structure and composition were mainly attributed to variations in organic carbon fluxes to the seafloor at the regional scale and nodule density at the local scale, thus supporting the main assumptions underlying the design of the APEIs. However, the APEI no. 3, which is located in an oligotrophic province and separated from the CCFZ by the Clarion Fracture Zone, showed the lowest densities, lowest diversity, and a very low and distant independent similarity in community composition compared to the contract areas, thus questioning the representativeness and the appropriateness of APEI no. 3 to meet its purpose of diversity preservation. Among the four exploration contracts, which belong to a mesotrophic province, the distance decay of similarity provided a species turnover of 0.04 species km−1, an average species range of 25 km and an extrapolated richness of up to 240 000 polychaete species in the CCFZ. By contrast, nonparametric estimators of diversity predict a regional richness of up to 498 species. Both estimates are biased by the high frequency of singletons in the dataset, which likely result from under-sampling and merely reflect our level of uncertainty. The assessment of potential risks and scales of biodiversity loss due to nodule mining thus requires an appropriate inventory of species richness in the CCFZ.

Description: Species distribution models (SDMs) correlate species occurrences with environmental predictors, and can be used to forecast distributions under future climates. SDMs have been criticized for not explicitly including the physiological processes underlying the species response to the environment. Recently, new methods have been suggested to combine SDMs with physiological estimates of performance (physiology-SDMs). In this study, we compare SDM and physiology-SDM predictions for select marine species in the Mediterranean Sea, a region subjected to exceptionally rapid climate change. We focused on six species and created physiology-SDMs that incorporate physiological thermal performance curves from experimental data with species occurrence records. We then contrasted projections of SDMs and physiology-SDMs under future climate (year 2100) for the entire Mediterranean Sea, and particularly the ‘warm’ trailing edge in the Levant region. Across the Mediterranean, we found cross-validation model performance to be similar for regular SDMs and physiology-SDMs. However, we also show that for around half the species the physiology-SDMs substantially outperform regular SDM in the warm Levant. Moreover, for all species the uncertainty associated with the coefficients estimated from the physiology-SDMs were much lower than in the regular SDMs. Under future climate, we find that both SDMs and physiology-SDMs showed similar patterns, with species predicted to shift their distribution north-west in accordance with warming sea temperatures. However, for the physiology-SDMs predicted distributional changes are more moderate than those predicted by regular SDMs. We conclude, that while physiology-SDM predictions generally agree with the regular SDMs, incorporation of the physiological data led to less extreme range shift forecasts. The results suggest that climate-induced range shifts may be less drastic than previously predicted, and thus most species are unlikely to completely disappear with warming climate. Taken together, the findings emphasize that physiological experimental data can provide valuable supplemental information to predict range shifts of marine species.

Description: Back‐arc basins open in response to subduction processes, which cause extension in the upper plate, usually along trench‐parallel spreading axes. However, global seismic databases reveal that the majority of seismic events in the Lau Basin occur along transcurrent (strike‐slip) rather than extensional faults. To better characterize active deformation in this region we compared Centroid Moment Tensors (CMTs), calculated for large (Mw 〉5), shallow (〈30 km) seismic events to the orientations of seafloor lineaments mapped throughout the Lau Basin. Ship‐based multibeam was combined with vertical gravity gradient data to provide comprehensive coverage to create the lineament map. By comparing the possible focal planes of the CMTs to the orientations of the lineaments, the most likely fault plane solutions were selected, thus classifying the faults and establishing the nature of the highly variable stress regimes in the basin. We resolved the strike, dip and dip direction of 308 faults, and classified 258 additional structures by fault type. The analysis highlights a stress regime that is dominated by a combination of left‐lateral and right‐lateral strike‐slip faults, large‐scale transcurrent motion along rigid crustal‐scale fault zones, and non‐rigid diffuse deformation along pre‐existing seafloor structures, with extension mainly limited to the tips of propagating rifts and spreading centers. By resolving many of the uncertain motions on the mapped lineaments of the Lau Basin, the CMT analysis addresses a number of questions concerning basin‐scale stress regimes and microplate development, complementing GPS measurements and providing a more complete picture of the complexities of back‐arc basin development.

Description: The Middle Eocene Climatic Optimum (MECO) was a gradual warming event and carbon cycle perturbation that occurred between 40.5 and 40.1 Ma. A number of characteristics, including greater-than-expected deep-sea carbonate dissolution, a lack of globally-coherent negative δ 13 40 C excursion in marine carbonates, a duration longer than the characteristic timescale of carbon-cycle recovery, and the absence of a clear trigger mechanism, challenge our current understanding of the Earth system and its regulatory feedbacks. This makes the MECO one of the most enigmatic events in the Cenozoic, dubbed a middle Eocene ‘carbon cycle conundrum’. Here we use boron isotopes in planktic foraminifera to better constrain pCO2 changes over the event. Over the MECO itself, we find that pCO2 rose by only 0.55-0.75 doublings, thus requiring a much more modest carbon injection than previously indicated by the alkenone δ 13 C-pCO2 proxy. In addition, this rise in pCO2 was focused around the peak of the 400 kyr warming trend. Before this, considerable global carbonate δ 18 O change was asynchronous with any coherent ocean pH (and hence pCO2) excursion. This finding suggests that middle Eocene climate (and perhaps a nascent cryosphere) was highly sensitive to small changes in radiative forcing.

Description: Climate change in Siberia is currently receiving a lot of attention because large permafrost-covered areas could provide a strong positive feedback to global warming through the release of carbon that has been sequestered there on glacial–interglacial timescales. Geological evidence and climate model experiments show that the Siberian region also played an exceptional role during glacial periods. The region that is currently known for its harsh cold climate did not experience major glaciations during the last ice age, including its severest stages around the Last Glacial Maximum (LGM). On the contrary, it is thought that glacial summer temperatures were comparable to the present day. However, evidence of glaciation has been found for several older glacial periods. We combine LGM experiments from the second and third phases of the Paleoclimate Modelling Intercomparison Project (PMIP2 and PMIP3) with sensitivity experiments using the Community Earth System Model (CESM). Together, these climate model experiments reveal that the intermodel spread in LGM summer temperatures in Siberia is much larger than in any other region of the globe and suggest that temperatures in Siberia are highly susceptible to changes in the imposed glacial boundary conditions, the included feedbacks and processes, and to the model physics of the different components of the climate model. We find that changes in the circumpolar atmospheric stationary wave pattern and associated northward heat transport drive strong local snow and vegetation feedbacks and that this combination explains the susceptibility of LGM summer temperatures in Siberia. This suggests that a small difference between two glacial periods in terms of climate, ice buildup or their respective evolution towards maximum glacial conditions can lead to strongly divergent summer temperatures in Siberia, allowing for the buildup of an ice sheet during some glacial periods, while during others, above-freezing summer temperatures preclude a multi-year snowpack from forming.

Description: In this paper we introduce a Bayesian framework, which is explicit about prior assumptions, for using model ensembles and observations together to constrain future climate change. The emergent constraint approach has seen broad application in recent years, including studies constraining the equilibrium climate sensitivity (ECS) using the Last Glacial Maximum (LGM) and the mid-Pliocene Warm Period (mPWP). Most of these studies were based on ordinary least squares (OLS) fits between a variable of the climate state, such as tropical temperature, and climate sensitivity. Using our Bayesian method, and considering the LGM and mPWP separately, we obtain values of ECS of 2.7 K (0.6–5.2, 5th–95th percentiles) using the PMIP2, PMIP3, and PMIP4 datasets for the LGM and 2.3 K (0.5–4.4) with the PlioMIP1 and PlioMIP2 datasets for the mPWP. Restricting the ensembles to include only the most recent version of each model, we obtain 2.7 K (0.7–5.2) using the LGM and 2.3 K (0.4–4.5) using the mPWP. An advantage of the Bayesian framework is that it is possible to combine the two periods assuming they are independent, whereby we obtain a tighter constraint of 2.5 K (0.8–4.0) using the restricted ensemble. We have explored the sensitivity to our assumptions in the method, including considering structural uncertainty, and in the choice of models, and this leads to 95 % probability of climate sensitivity mostly below 5 K and only exceeding 6 K in a single and most uncertain case assuming a large structural uncertainty. The approach is compared with other approaches based on OLS, a Kalman filter method, and an alternative Bayesian method. An interesting implication of this work is that OLS-based emergent constraints on ECS generate tighter uncertainty estimates, in particular at the lower end, an artefact due to a flatter regression line in the case of lack of correlation. Although some fundamental challenges related to the use of emergent constraints remain, this paper provides a step towards a better foundation for their potential use in future probabilistic estimations of climate sensitivity.

Description: Rivers are a major source of nutrients, carbon and alkalinity to the global ocean. In this study, we firstly estimate pre-industrial riverine loads of nutrients, carbon and alkalinity based on a hierarchy of weathering and terrestrial organic matter export models, while identifying regional hotspots of the riverine exports. Secondly, we implement the riverine loads into a global ocean biogeochemical model to describe their implications for oceanic nutrient concentrations, net primary production (NPP) and air–sea CO2 fluxes globally, as well as in an analysis of coastal regions. Thirdly, we quantitatively assess the terrestrial origins and the long-term fate of riverine carbon in the ocean. We quantify annual bioavailable pre-industrial riverine loads of 3.7 Tg P, 27 Tg N, 158 Tg Si and 603 Tg C delivered to the ocean globally. We thereby identify the tropical Atlantic catchments (20 % of global C), Arctic rivers (9 % of global C) and Southeast Asian rivers (15 % of global C) as dominant suppliers of carbon for the ocean. The riverine exports lead to a simulated net global oceanic CO2 source of 231 Tg C yr−1 to the atmosphere, which is mainly caused by inorganic carbon (source of 183 Tg C yr−1) and by organic carbon (source of 128 Tg C yr−1) riverine loads. Additionally, a sink of 80 Tg C yr−1 is caused by the enhancement of the biological carbon uptake from dissolved inorganic nutrient inputs from rivers and the resulting alkalinity production. While large outgassing fluxes are simulated mostly in proximity to major river mouths, substantial outgassing fluxes can be found further offshore, most prominently in the tropical Atlantic. Furthermore, we find evidence for the interhemispheric transfer of carbon in the model; we detect a larger relative outgassing flux (49 % of global riverine-induced outgassing) in the Southern Hemisphere in comparison to the hemisphere's relative riverine inputs (33 % of global C inputs), as well as an outgassing flux of 17 Tg C yr−1 in the Southern Ocean. The addition of riverine loads in the model leads to a strong NPP increase in the tropical west Atlantic, Bay of Bengal and the East China Sea (+166 %, +377 % and +71 %, respectively). On the light-limited Arctic shelves, the NPP is not strongly sensitive to riverine loads, but the CO2 flux is strongly altered regionally due to substantial dissolved inorganic and organic carbon supplies to the region. While our study confirms that the ocean circulation remains the main driver for biogeochemical distributions in the open ocean, it reveals the necessity to consider riverine inputs for the representation of heterogeneous features in the coastal ocean and to represent riverine-induced pre-industrial carbon outgassing in the ocean. It also underlines the need to consider long-term CO2 sources from volcanic and shale oxidation fluxes in order to close the framework's atmospheric carbon budget.

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: Key points:  First insights into the crustal structure of the northeastern Lau Basin, along a 290 km transect at 17°20’S.  Crust in southern Fonualei Rift and Spreading Center was created by extension of arc crust and variable amount of magmatism.  Magmatic underplating is present in some parts of the southern Niuafo’ou Microplate The northeastern Lau Basin is one of the fastest opening and magmatically most active back‐arc regions on Earth. Although the current pattern of plate boundaries and motions in this complex mosaic of microplates is reasonably understood, the internal structure and evolution of the back‐arc crust are not. We present new geophysical data from a 290 km long east‐west oriented transect crossing the Niuafo’ou Microplate (back‐arc), the Fonualei Rift and Spreading Centre (FRSC) and the Tofua Volcanic Arc at 17°20’S. Our P‐wave tomography model and density modelling suggests that past crustal accretion inside the southern FRSC was accommodated by a combination of arc crustal extension and magmatic activity. The absence of magnetic reversals inside the FRSC supports this and suggests that focused seafloor spreading has until now not contributed to crustal accretion. The back‐arc crust constituting the southern Niuafo’ou Microplate reveals a heterogeneous structure comprising several crustal blocks. Some regions of the back‐arc show a crustal structure similar to typical oceanic crust, suggesting they originate from seafloor spreading. Other crustal blocks resemble a structure that is similar to volcanic arc crust or a ‘hydrous’ type of oceanic crust that has been created at a spreading center influenced by slab‐derived water at distances 〈 50 km to the arc. Throughout the back‐arc region we observe a high‐velocity (Vp 7.2‐7.5 km s‐1) lower crust, which is an indication for magmatic underplating, which is likely sustained by elevated upper mantle temperatures in this region.

Description: Southern hemisphere lower stratospheric ozone depletion has been shown to lead to a poleward shift of the tropospheric jet stream during austral summer, influencing surface atmosphere and ocean conditions, such as surface temperatures and sea ice extent. The characteristics of stratospheric and tropospheric responses to ozone depletion, however, differ largely among climate models depending on the representation of ozone in the models. The most accurate way to represent ozone in a model is to calculate it interactively. However, due to computational costs, in particular for long-term coupled ocean-atmosphere model integrations, the more common way is to prescribe ozone from observations or calculated model fields. Here, we investigate the difference between an interactive and a specified chemistry version of the same atmospheric model in a fully-coupled setup using a 9-member chemistry-climate model ensemble. In the specified chemistry version of the model the ozone fields are prescribed using the output from the interactive chemistry model version. In contrast to earlier studies, we use daily-resolved ozone fields in the specified chemistry simulations to achieve a better comparability between the ozone forcing with and without interactive chemistry. We find that although the short-wave heating rate trend in response to ozone depletion is the same in the different chemistry settings, the interactive chemistry ensemble shows a stronger trend in polar cap stratospheric temperatures (by about 0.7 K per decade) and circumpolar stratospheric zonal mean zonal winds (by about 1.6 m/s per decade) as compared to the specified chemistry ensemble. This difference between interactive and specified chemistry in the stratospheric response to ozone depletion also affects the tropospheric response, namely the poleward shift of the tropospheric jet stream. We attribute part of these differences to the missing representation of feedbacks between chemistry and dynamics in the specified chemistry ensemble, which affect the dynamical heating rates, and part of it to the lack of spatial asymmetries in the prescribed ozone fields. This effect is investigated using a sensitivity ensemble that was forced by a three-dimensional instead of a two–dimensional ozone field. This study emphasizes the value of interactive chemistry for the representation of the southern hemisphere tropospheric jet response to ozone depletion and infers that for periods with strong ozone variability (trends) the details of the ozone forcing can be crucial for representing southern hemispheric climate variability.

Description: Morphological changes in coccoliths, tiny calcite platelets covering the outer surface of coccolithophores, can be induced by physiological responses to environmental changes. Coccoliths recovered from sedimentary successions may therefore provide information on paleo-environmental conditions prevailing at the time when the coccolithophores were alive. To calibrate the biomineralization responses of ancient coccolithophore to environmental changes, studies often compared the biological responses of living coccolithophore species with paleo-data from calcareous nannofossils. However, there is uncertainty whether the morphological responses of living coccolithophores are representative of those of the fossilized ancestors. To investigate this, we exposed four living coccolithophore species (Emiliania huxleyi, Gephyrocapsa oceanica, Coccolithus pelagicus subsp. braarudii, and Pleurochrysis carterae) that have been evolutionarily distinct for hundreds of thousands to millions of years, to a range of environmental conditions (i.e., changing light intensity, Mg∕Ca ratio, nutrient availability, temperature, and carbonate chemistry) and evaluated their responses in coccolith morphology (i.e., size, length, width, malformation). The motivation for this study was to test if there is a consistent morphological response of the four species to changes in any of the tested abiotic environmental factors. If this was the case, then this could suggest that coccolith morphology can serve as a paleo-proxy for that specific factor because this response is conserved across species that have been evolutionary distinct over geological timescales. However, we found that the four species responded differently to changing light intensity, Mg∕Ca ratio, nutrient availability, and temperature in terms of coccolith morphology. The lack of a common response reveals the difficulties in using coccolith morphology as a paleo-proxy for these environmental drivers. However, a common response was observed under changing seawater carbonate chemistry (i.e., rising CO2), which consistently induced malformations. This commonality provides some confidence that malformations found in the sedimentary record could be indicative of adverse carbonate chemistry conditions.