ALBERT

Description: The possibility of using the 15% excess U234 activity in oceanic uranium for dating pelagic sediments in the age range 100,000 years to more than 1 m.y. has been explored. Results from a series of analyses of bulk samples, mechanical separates, and acid leach fractions indicate that separation of authigenic uranium from detrital uranium by either mechanical or chemical means is impractical. Measurements on totally dissolved samples reveal that the sediments do not form a closed system; post-depositional migration of U234 in the sedimentary column takes place. Based on the experimental data obtained from three red-clay cores with sedimentation rates ranging from 2 to 6 mm/1000 yr, a model depicting diffusion of the U234 generated within the sediments is proposed. The diffusion equation includes three parameters: sedimentation rate, diffusion coefficient for U234, and fraction of the internally produced U234 subject to mobility. If the amount of U234 lost from these cores is typical, a sizeable part of the U234 excess in the sea must be from this source.

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: 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: 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: 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: 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: 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: 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: 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 Arctic Ocean is particularly vulnerable to ocean acidification, a process that is mainly driven by the uptake of anthropogenic carbon (Cant) from the atmosphere. Although Cant concentrations cannot be measured directly in the ocean, they have been estimated using data‐based methods such as the transient time distribution (TTD) approach, which characterizes the ventilation of water masses with inert transient tracers, such as CFC‐12. Here, we evaluate the TTD approach in the Arctic Ocean using an eddying ocean model as a test bed. When the TTD approach is applied to simulated CFC‐12 in that model, it underestimates the same model's directly simulated Cant concentrations by up to 12%, a bias that stems from its idealized assumption of gas equilibrium between atmosphere and surface water, both for CFC‐12 and anthropogenic CO2. Unlike the idealized assumption, the simulated partial pressure of CFC‐12 (p CFC‐12) in Arctic surface waters is undersaturated relative to that in the atmosphere in regions and times of deep‐water formation, while the simulated equivalent for Cant is supersaturated. After accounting for the TTD approach's negative bias, the total amount of Cant in the Arctic Ocean in 2005 increases by 8% to 3.3 ± 0.3 Pg C. By combining the adjusted TTD approach with scenarios of future atmospheric CO2, it is estimated that all Arctic waters, from surface to depth, would become corrosive to aragonite by the middle of the next century even if atmospheric CO2 could be stabilized at 540 ppm.

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: We investigate the lag between warm interannual Sea Surface Temperature (SST) events in the eastern-equatorial Atlantic, the Atlantic Niños, and the occurrence of Benguela Niños along the southwestern Angolan coast. While it is commonly agreed that both events are associated with equatorial and subsequent coastal-trapped wave propagations driven remotely by a relaxation of the trade-winds, it is surprising that SST anomalies off Angola tend to precede the ones in the eastern-equatorial sector by ~1 month. To explain this counterintuitive behavior, our methodology is based on the experimentation with a Tropical Atlantic Ocean model. Using idealized wind-stress perturbations from a composite analysis, we trigger warm equatorial and coastal events over a stationary and then, seasonally varying ocean mean-state. In agreement with the linear dynamics, our results show that when the interannual wind-stress forcing is restricted to the western-central equatorial Atlantic, the model yields equatorial events leading the coastal ones. This implies that neither the differences in the ocean stratification between the two regions (thermocline depths or modal wave contributions) nor the seasonal phasing of the events explains the observed temporal sequence. Only if wind-stress anomalies are also prescribed in the coastal fringe, the coastal warming precedes the eastern-equatorial SST anomaly peak, emphasizing the role of the local forcing in the phenology of Benguela Niños. A weaker South-Atlantic Anticyclone initiates the coastal warming before the development of eastern-equatorial SST anomalies. Then, equatorward coastal wind anomalies, driven by a convergent anomalous circulation located on the warm Atlantic Niño, stop the remotely forced coastal warming prematurely.

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: 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: 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: 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: 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: 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: 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 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: 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: Observations during the satellite era 1979–2018 only depict small sea surface temperature (SST) trends over the Equatorial Atlantic cold tongue region in boreal summer. This lack of surface warming of the cold tongue, termed warming hole here, denotes an 11% amplification of the mean SST annual cycle. The warming hole is driven by a shoaling of the equatorial thermocline, linked to increased wind stress forcing, and damped by the surface turbulent heat fluxes. The satellite era warming deficit is not unusual during the twentieth century—similar weak trends were also observed during the 1890s–1910s and 1940s–1960s. The tendency for surface cooling appears to reflect an interaction of external forcing, which controls the timing and magnitude of the cooling, with the intrinsic variability of the climate system. The hypothesis for externally forced modulation of internal variability is supported by climate model simulations forced by the observed time-varying concentrations of atmospheric greenhouse gases and natural aerosols. These show that increased greenhouse forcing warmed the cold tongue and aerosols cooled it during the satellite era. However, internal variability, as derived from control integrations with fixed, preindustrial values of greenhouse gases and aerosols, can potentially cause larger cooling than observed during the satellite era. Large uncertainties remain on the relative roles of external forcing and intrinsic variability in both observations and coupled climate models.

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.

Description: High‐temperature hydrothermal venting has been discovered on all modern mid‐ocean ridges at all spreading rates. Although significant strides have been made in understanding the underlying processes that shape such systems, several first‐order discrepancies between model predictions and observations remain. One key paradox is that numerical experiments consistently show entrainment of cold ambient seawater in shallow high permeability ocean crust causing a temperature drop that is difficult to reconcile with high vent temperatures. We investigate this conundrum using a thermo‐hydro‐chemical model that couples hydrothermal fluid flow with anhydrite‐ and pyrite‐forming reactions in the shallow subseafloor. The models show that precipitation of anhydrite in warming seawater and in cooling hydrothermal fluids during mixing results in the formation of a chimney‐like subseafloor structure around the upwelling, high‐temperature plume. The establishment of such anhydrite‐sealed zones reduces mixing between the hydrothermal fluid and seawater and results in an increase in vent temperature. Pyrite subsequently precipitates close to the seafloor within the anhydrite chimney. Although anhydrite thus formed may be dissolved when colder seawater circulates through the crust away from the spreading axis, the inside pyrite walls would be preserved as veins in present‐day metal deposits, thereby preserving the history of hydrothermal circulation through shallow oceanic crust.

Description: An optimized method is presented to determine dissolved free (DFCHO) and dissolved combined carbohydrates (DCCHO) in saline matrices, such as oceanic seawater, Arctic ice core samples or brine using a combination of a desalination with electro-dialysis (ED) and high-performance anion exchange chromatography coupled to pulsed amperometric detection (HPAEC-PAD). Free neutral sugars, such as glucose and galactose, were found with 95 %–98 % recovery rates. Free amino sugars and free uronic acids were strongly depleted during ED at pH=8, but an adjustment of the pH could result in higher recoveries (58 %–59 % for amino sugars at pH=11; 45 %–49 % for uronic acids at pH=1.5). The applicability of this method for the analysis of DCCHO was evaluated with standard solutions and seawater samples compared with another established desalination method using membrane dialysis. DFCHO in field samples from different regions on Earth ranged between 11 and 118 nM and DCCHO between 260 and 1410 nM. This novel method has the potential to contribute to a better understanding of biogeochemical processes in the oceans and sea–air transfer processes of organic matter into the atmosphere in future studies.

Description: The Labrador Sea is important for the modern global thermohaline circulation system through the formation of intermediate Labrador Sea Water (LSW) that has been hypothesized to stabilize the modern mode of North Atlantic deep-water circulation. The rate of LSW formation is controlled by the amount of winter heat loss to the atmosphere, the expanse of freshwater in the convection region and the inflow of saline waters from the Atlantic. The Labrador Sea, today, receives freshwater through the East and West Greenland currents (EGC, WGC) and the Labrador Current (LC). Several studies have suggested the WGC to be the main supplier of freshwater to the Labrador Sea, but the role of the southward flowing LC in Labrador Sea convection is still debated. At the same time, many paleoceanographic reconstructions from the Labrador Shelf focussed on late deglacial to early Holocene meltwater run-off from the Laurentide Ice Sheet (LIS), whereas little information exists about LC variability since the final melting of the LIS about 7000 years ago. In order to enable better assessment of the role of the LC in deep-water formation and its importance for Holocene climate variability in Atlantic Canada, this study presents high-resolution middle to late Holocene records of sea surface and bottom water temperatures, freshening, and sea ice cover on the Labrador Shelf during the last 6000 years. Our records reveal that the LC underwent three major oceanographic phases from the mid- to late Holocene. From 6.2 to 5.6 ka, the LC experienced a cold episode that was followed by warmer conditions between 5.6 and 2.1 ka, possibly associated with the late Holocene thermal maximum. While surface waters on the Labrador Shelf cooled gradually after 3 ka in response to the neoglaciation, Labrador Shelf subsurface or bottom waters show a shift to warmer temperatures after 2.1 ka. Although such an inverse stratification by cooling of surface and warming of subsurface waters on the Labrador Shelf would suggest a diminished convection during the last 2 millennia compared to the mid-Holocene, it remains difficult to assess whether hydrographic conditions in the LC have had a significant impact on Labrador Sea deep-water formation.

Description: Recent winters have been unique due to the rapid and extreme cooling of the subpolar North Atlantic. Here, we present a novel view on its causes and consequences. Combining in-situ observations with remote sensing and atmospheric reanalysis data, we show that increased freshening of the subpolar region gives rise to a faster surface cooling in fall and winter. Large freshwater events, in particular, result in pronounced cold anomalies with sharp temperature gradients that promote an enhanced storminess. The storms reinforce the cooling by driving stronger heat losses and modulating the surface flow. Consistent with this mechanism, past freshwater events have been followed by cold anomalies in winter of approximately −2°C and increases in the North Atlantic Oscillation index of up to ∼0.6 within 3 years. We expect that future freshwater discharges into the North Atlantic will amplify the cold anomaly and trigger an enhanced wintertime storminess with far-reaching climatic implications.

Description: Key Points: - Fluid flow is focused along Nootka Fault traces resulting in shallow bright spots - Two seafloor mounds are the result of basaltic intrusions in the Nootka Fault zone - Gas hydrates occur at the Nootka Slope and are imaged seismically as bottom- simulating reflectors suggesting a regional heat-flow of ~80 mW/m2 along the slope Abstract Geophysical and geochemical data indicate there is abundant fluid expulsion in the Nootka fault zone (NFZ) between the Juan de Fuca and Explorer plates and the Nootka continental slope. Here we combine observations from 〉 20 years of investigations to demonstrate the nature of fluid‐flow along the NFZ, which is the seismically most active region off Vancouver Island. Seismicity reaching down to the upper mantle is linked to near‐seafloor manifestation of fluid flow through a network of faults. Along the two main fault traces, seismic reflection data imaged bright spots 100 – 300 m below seafloor that lie above changes in basement topography. The bright spots are conformable to sediment layering, show opposite‐to‐seafloor reflection polarity, and are associated with frequency‐reduction and velocity push‐down indicating the presence of gas in the sediments. Two seafloor mounds ~15 km seaward of the Nootka slope are underlain by deep, non‐conformable high amplitude reflective zones. Measurements in the water column above one mound revealed a plume of warm water, and bottom‐video observations imaged hydrothermal vent system biota. Pore fluids from a core at this mound contain predominately microbial methane (C1) with a high proportion of ethane (C2) yielding C1/C2 ratios 〈 500 indicating a possible slight contribution from a deep source. We infer the reflective zones beneath the two mounds are basaltic intrusions that create hydrothermal circulation within the overlying sediments. Across the Nootka continental slope, gas hydrate related bottom‐simulating reflectors are widespread and occur at depths indicating heat‐flow values of 80 – 90 mW/m2.

Description: The distribution of dissolved silicon isotopes (δ30Si) was examined along the US GEOTRACES East Pacific Zonal Transect (GP16) extending from Peru to Tahiti (10°S and 15°S latitude). Surface waters in the subtropical gyre displayed high δ30Si due to strong utilization of silicic acid (DSi). In contrast, surface waters close to the Peruvian coast where upwelling prevailed were less depleted and only moderately fractionated. δ30Si of water masses along the transect were compared with the results of an Optimum Multiparameter Analysis that quantified the fractional contributions of endmember water masses in each sample. Strong admixture of intermediate waters obscured the expected heavy isotopic signatures of Subantarctic Mode Water and Antarctic Intermediate Water. Isotope values were nearly homogenous below 2000 m (Average: +1.3 ± 0.1 ‰, 1 s.d.) despite the 25 μmol kg‐1 range in the DSi content among water masses. This homogeneity confirms prior observations and model results that predict nearly constant δ30Si values of +1.0 to +1.2 ‰ for Pacific deep waters with [DSi] 〉 100 μmol kg‐1. Waters above the East Pacific Rise (EPR) influenced by hydrothermal activity showed a small increase in [DSi] together with dissolved iron, but overall stations close to the EPR were slightly depleted in [DSi] (3 to 6 μmol kg‐1) with no significant shift in δ30Si compared to adjacent waters. Hydrothermal [DSi] appears to precipitate within the conduit of the EPR or upon contact with cold seawater resulting in a negligible influence of hydrothermal fluids on δ30Si in deep water. Key Points Surface waters have a large range in dissolved silicon isotopes covering nutrient‐rich coastal upwelling to oligotrophic waters Deep water masses with DSi concentrations 〉 100 μmol kg‐1 show homogenous silicon isotope signatures despite up to 25 μmol kg‐1 differences in [DSi] Hydrothermal fluids have a negligible effect on Si isotope distributions in the deep Pacific

Description: We simulate the two Coupled Model Intercomparison Project scenarios RCP4.5 and RCP8.5, to assess the effects of melt‐induced fresh water on the Atlantic meridional overturning circulation (AMOC). We use a newly developed climate model with high resolution at the coasts, resolving the complex ocean dynamics. Our results show an AMOC recovery in simulations run with and without an included ice sheet model. We find that the ice sheet adds a strong decadal variability on the freshwater release, resulting in intervals in which it reduces the surface runoff by high accumulation rates. This compensating effect is missing in climate models without dynamic ice sheets. Therefore, we argue to assess those freshwater hosing experiments critically, which aim to parameterize Greenland's freshwater release. We assume the increasing net evaporation over the Atlantic and the resulting increase in ocean salinity, to be the main driver of the AMOC recovery.

Description: Nowadays various methods and sensors are available for 3D reconstruction tasks; however, it is still necessary to integrate advantages of different technologies for optimizing the quality 3D models. Computed tomography (CT) is an imaging technique which takes a large number of radiographic measurements from different angles, in order to generate slices of the object, however, without colour information. The aim of this study is to put forward a framework to extract colour information from photogrammetric images for corresponding Computed Tomography (CT) surface data with high precision. The 3D models of the same object from CT and photogrammetry methods are generated respectively, and a transformation matrix is determined to align the extracted CT surface to the photogrammetric point cloud through a coarse-to-fine registration process. The estimated pose information of images to the photogrammetric point clouds, which can be obtained from the standard image alignment procedure, also applies to the aligned CT surface data. For each camera pose, a depth image of CT data is calculated by projecting all the CT points to the image plane. The depth image is in principle should agree with the corresponding photogrammetric image. The points, which cannot be seen from the pose, but are also projected on the depth image, are excluded from the colouring process. This is realized by comparing the range values of neighbouring pixels and finding the corresponding 3D points with larger range values. The same procedure is implemented for all the image poses to obtain the coloured CT surface. Thus, by using photogrammetric images, we achieve a coloured CT dataset with high precision, which combines the advantages from both methods. Rather than simply stitching different data, we deep-dive into the photogrammetric 3D reconstruction process and optimize the CT data with colour information. This process can also provide an initial route and more options for other data fusion processes.

Description: The coastal ocean is strongly affected by ocean acidification because of its shallow water depths, low volume, and the closeness to terrestrial dynamics. Earlier observations of dissolved inorganic carbon (DIC) and total alkalinity (TA) in the southern part of the North Sea, a northwest European shelf sea, revealed lower acidification effects than expected. It has been assumed that anaerobic degradation and subsequent TA release in the adjacent back-barrier tidal areas (Wadden Sea) in summertime is responsible for this phenomenon. In this study the exchange rates of TA and DIC between the Wadden Sea tidal basins and the North Sea and the consequences for the carbonate system in the German Bight are estimated using a 3D ecosystem model. The aim of this study is to differentiate the various sources contributing to observed high summer TA in the southern North Sea. Measured TA and DIC in the Wadden Sea are considered as model boundary conditions. This procedure acknowledges the dynamic behaviour of the Wadden Sea as an area of effective production and decomposition of organic material. According to the modelling results, 39 Gmol TA yr−1 were exported from the Wadden Sea into the North Sea, which is less than a previous estimate but within a comparable range. The interannual variabilities in TA and DIC, mainly driven by hydrodynamic conditions, were examined for the years 2001–2009. Dynamics in the carbonate system are found to be related to specific weather conditions. The results suggest that the Wadden Sea is an important driver for the carbonate system in the southern North Sea. On average 41 % of TA inventory changes in the German Bight were caused by riverine input, 37 % by net transport from adjacent North Sea sectors, 16 % by Wadden Sea export, and 6 % were caused by internal net production of TA. The dominant role of river input for the TA inventory disappears when focusing on TA concentration changes due to the corresponding freshwater fluxes diluting the marine TA concentrations. The ratio of exported TA versus DIC reflects the dominant underlying biogeochemical processes in the Wadden Sea. Whereas aerobic degradation of organic matter played a key role in the North Frisian Wadden Sea during all seasons of the year, anaerobic degradation of organic matter dominated in the East Frisian Wadden Sea. Despite the scarcity of high-resolution field data, it is shown that anaerobic degradation in the Wadden Sea is one of the main contributors of elevated summer TA values in the southern North Sea.

Description: Acetone is one of the most abundant oxygenated volatile organic compounds (VOCs) in the atmosphere. The oceans impose a strong control on atmospheric acetone, yet the oceanic fluxes of acetone remain poorly constrained. In this work, the global budget of acetone is evaluated using two global models: CAM‐chem and GEOS‐Chem. CAM‐chem uses an online air‐sea exchange framework to calculate the bidirectional oceanic acetone fluxes, which is coupled to a data‐oriented machine‐learning approach. The machine‐learning algorithm is trained using a global suite of seawater acetone measurements. GEOS‐Chem uses a fixed surface seawater concentration of acetone to calculate the oceanic fluxes. Both model simulations are compared to airborne observations from a recent global‐scale, multiseasonal campaign, the NASA Atmospheric Tomography Mission (ATom). We find that both CAM‐chem and GEOS‐Chem capture the measured acetone vertical distributions in the remote atmosphere reasonably well. The combined observational and modeling analysis suggests that (i) the ocean strongly regulates the atmospheric budget of acetone. The tropical and subtropical oceans are mostly a net source of acetone, while the high‐latitude oceans are a net sink. (ii) CMIP6 anthropogenic emission inventory may underestimate acetone and/or its precursors in the Northern Hemisphere. (iii) The MEGAN biogenic emissions model may overestimate acetone and/or its precursors, and/or the biogenic oxidation mechanisms may overestimate the acetone yields. (iv) The models consistently overestimate acetone in the upper troposphere‐lower stratosphere over the Southern Ocean in austral winter. (v) Acetone contributes up to 30–40% of hydroxyl radical production in the tropical upper troposphere/lower stratosphere.

Description: Eastern boundary upwelling systems (EBUS) are among the most productive marine ecosystems on Earth. The production of organic material is fueled by upwelling of nutrient-rich deep waters and high incident light at the sea surface. However, biotic and abiotic factors can modify surface production and related biogeochemical processes. Determining these factors is important because EBUS are considered hotspots of climate change, and reliable predictions of their future functioning requires understanding of the mechanisms driving the biogeochemical cycles therein. In this field experiment, we used in situ mesocosms as tools to improve our mechanistic understanding of processes controlling organic matter cycling in the coastal Peruvian upwelling system. Eight mesocosms, each with a volume of ∼55 m3, were deployed for 50 d ∼6 km off Callao (12∘ S) during austral summer 2017, coinciding with a coastal El Niño phase. After mesocosm deployment, we collected subsurface waters at two different locations in the regional oxygen minimum zone (OMZ) and injected these into four mesocosms (mixing ratio ≈1.5 : 1 mesocosm: OMZ water). The focus of this paper is on temporal developments of organic matter production, export, and stoichiometry in the individual mesocosms. The mesocosm phytoplankton communities were initially dominated by diatoms but shifted towards a pronounced dominance of the mixotrophic dinoflagellate (Akashiwo sanguinea) when inorganic nitrogen was exhausted in surface layers. The community shift coincided with a short-term increase in production during the A. sanguinea bloom, which left a pronounced imprint on organic matter C : N : P stoichiometry. However, C, N, and P export fluxes did not increase because A. sanguinea persisted in the water column and did not sink out during the experiment. Accordingly, export fluxes during the study were decoupled from surface production and sustained by the remaining plankton community. Overall, biogeochemical pools and fluxes were surprisingly constant for most of the experiment. We explain this constancy by light limitation through self-shading by phytoplankton and by inorganic nitrogen limitation which constrained phytoplankton growth. Thus, gain and loss processes remained balanced and there were few opportunities for blooms, which represents an event where the system becomes unbalanced. Overall, our mesocosm study revealed some key links between ecological and biogeochemical processes for one of the most economically important regions in the oceans.

Description: We present a new framework for global ocean- sea-ice model simulations based on phase 2 of the Ocean Model Intercomparison Project (OMIP-2), making use of the surface dataset based on the Japanese 55-year atmospheric reanalysis for driving ocean-sea-ice models (JRA55-do).We motivate the use of OMIP-2 over the framework for the first phase of OMIP (OMIP-1), previously referred to as the Coordinated Ocean-ice Reference Experiments (COREs), via the evaluation of OMIP-1 and OMIP-2 simulations from 11 state-of-the-science global ocean-sea-ice models. In the present evaluation, multi-model ensemble means and spreads are calculated separately for the OMIP-1 and OMIP-2 simulations and overall performance is assessed considering metrics commonly used by ocean modelers. Both OMIP-1 and OMIP-2 multi-model ensemble ranges capture observations in more than 80% of the time and region for most metrics, with the multi-model ensemble spread greatly exceeding the difference between the means of the two datasets. Many features, including some climatologically relevant ocean circulation indices, are very similar between OMIP-1 and OMIP- 2 simulations, and yet we could also identify key qualitative improvements in transitioning from OMIP-1 to OMIP- 2. For example, the sea surface temperatures of the OMIP- 2 simulations reproduce the observed global warming during the 1980s and 1990s, as well as the warming slowdown in the 2000s and the more recent accelerated warming, which were absent in OMIP-1, noting that the last feature is part of the design of OMIP-2 because OMIP-1 forcing stopped in 2009. A negative bias in the sea-ice concentration in summer of both hemispheres in OMIP-1 is significantly reduced in OMIP-2. The overall reproducibility of both seasonal and interannual variations in sea surface temperature and sea surface height (dynamic sea level) is improved in OMIP-2. These improvements represent a new capability of the OMIP-2 framework for evaluating processlevel responses using simulation results. Regarding the sensitivity of individual models to the change in forcing, the models show well-ordered responses for the metrics that are directly forced, while they show less organized responses for those that require complex model adjustments. Many of the remaining common model biases may be attributed either to errors in representing important processes in ocean-sea-ice models, some of which are expected to be reduced by using finer horizontal and/or vertical resolutions, or to shared biases and limitations in the atmospheric forcing. In particular, further efforts are warranted to resolve remaining issues in OMIP-2 such as the warm bias in the upper layer, the mismatch between the observed and simulated variability of heat content and thermosteric sea level before 1990s, and the erroneous representation of deep and bottom water formations and circulations. We suggest that such problems can be resolved through collaboration between those developing models (including parameterizations) and forcing datasets. Overall, the present assessment justifies our recommendation that future model development and analysis studies use the OMIP-2 framework.

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

Description: Publications on temperate deciduous tree refugia in Europe are abundant, but little is known about the patterns of temperate tree refugia in eastern Asia, an area where biodiversity survived Quaternary glaciations and which has the world's most diverse temperate flora. Our goal is to compare climate model simulations with pollen data in order to establish the location of glacial refugia during the Last Glacial Maximum (LGM). Limits in which temperate deciduous trees can survive are taken from the literature. The model outputs are first tested for the present by comparing climate models with published modern pollen data. As this method turned out to be satisfactory for the present, the same approach was used for the LGM. Climate model simulations (ECHAM5 T106), statistically further downscaled, are used to infer the temperate deciduous tree distribution during the LGM. These were compared with available fossil temperate tree pollen occurrences. The impact of the LGM on the eastern Asian climate was much weaker than on the European climate. The area of possible tree growth shifts only by about 2∘ to the south between the present and the LGM. This contributes to explaining the greater biodiversity of forests in eastern Asia compared to Europe. Climate simulations and the available, although fractional, fossil pollen data agree. Therefore, climate estimations can safely be used to fill areas without pollen data by mapping potential refugia distributions. The results show two important areas with population connectivity: the Yellow Sea emerged shelf and the southern Himalayas. These two areas were suitable for temperate deciduous tree growth, providing corridors for population migration and connectivity (i.e. less population fragmentation) in glacial periods. Many tree populations live in interglacial refugia, not glacial ones. The fact that the model simulation for the LGM fits so well with observed pollen distribution is another indication that the model used is good enough to also simulate the LGM period.

Description: Key Points: • Observational transport time series of the Atlantic Subtropical Cells reveals dominant seasonal variability for horizontal branches • On time scales longer than ~5 years, interior thermocline layer transport convergence modulates equatorial sea surface temperature anomalies • Western boundary current and interior transport anomalies are partly compensating each other at thermocline level on all time scales The shallow meridional overturning cells of the Atlantic Ocean, the subtropical cells (STCs), consist of poleward Ekman transport at the surface, subduction in the subtropics, equatorward flow at thermocline level and upwelling along the equator and at the eastern boundary. In this study, we provide the first observational estimate of transport variability associated with the horizontal branches of the Atlantic STCs in both hemispheres based on Argo float data and supplemented by reanalysis products. Thermocline layer transport convergence and surface layer transport divergence between 10°N and 10°S are dominated by seasonal variability. Meridional thermocline layer transport anomalies at the western boundary and in the interior basin are anti‐correlated and partially compensate each other at all resolved time scales. It is suggested that the seesaw‐like relation is forced by the large‐scale off‐equatorial wind stress changes through low‐baroclinic‐mode Rossby wave adjustment. We further show that anomalies of the thermocline layer interior transport convergence modulate sea surface temperature (SST) variability in the upwelling regions along the equator and at the eastern boundary at time scales longer than 5 years. Phases of weaker (stronger) interior transport are associated with phases of higher (lower) equatorial SST. At these time scales, STC transport variability is forced by off‐equatorial wind stress changes, especially by those in the southern hemisphere. At shorter time scales, equatorial SST anomalies are, instead, mainly forced by local changes of zonal wind stress.

Description: Stable oxygen isotope records from central Greenland suggest disproportionally large long‐term surface air temperature (SAT) variability during the Last Glacial Maximum (LGM) relative to preindustrial times. Large perturbations in mean atmospheric circulation and its variability forced by extensive Northern Hemisphere ice sheet coverage have been suggested as cause for the enhanced Greenland SAT variability. Here, we assess the factors driving Greenland SAT variability during the LGM by means of dedicated climate model simulations and find remote forcing from the Pacific of critical importance. Atmospheric teleconnections from the Interdecadal Pacific Oscillation (IPO), a multidecadal oscillation of sea surface temperature in the Pacific Ocean, strongly intensify under LGM conditions, driving enhanced surface wind variability over Greenland, which in turn amplifies SAT variability by anomalous atmospheric heat transport. A major role of the IPO in forcing Greenland SAT variability also is supported by a number of models from the Paleoclimate Modeling Intercomparison Project Phase III

Description: Lipids, in their function as trophic markers in food webs and organic matter source indicators in the water column and sediments, provide a tool for reconstructing the complexity of global change effects on aquatic ecosystems. It remains unclear how ongoing changes in multiple environmental drivers affect the production of key lipid biomarkers in marine phytoplankton. Here, we tested the responses of sterols, alkenones and fatty acids (FAs) in the diatom Phaeodactylum tricornutum, the cryptophyte Rhodomonas sp. and the haptophyte Emiliania huxleyi under a full-factorial combination of three temperatures (12, 18 and 24 ∘C), three N : P supply ratios (molar ratios 10 : 1, 24 : 1 and 63 : 1) and two pCO2 levels (560 and 2400 µatm) in semicontinuous culturing experiments. Overall, N and P deficiency had a stronger effect on per-cell contents of sterols, alkenones and FAs than warming and enhanced pCO2. Specifically, P deficiency caused an overall increase in biomarker production in most cases, while N deficiency, warming and high pCO2 caused nonsystematic changes. Under future ocean scenarios, we predict an overall decrease in carbon-normalized contents of sterols and polyunsaturated fatty acids (PUFAs) in E. huxleyi and P. tricornutum and a decrease in sterols but an increase in PUFAs in Rhodomonas sp. Variable contents of lipid biomarkers indicate a diverse carbon allocation between marine phytoplankton species in response to changing environments. Thus, it is necessary to consider the changes in key lipids and their consequences for food-web dynamics and biogeochemical cycles, when predicting the influence of global change on marine ecosystems.

Description: The gridded sea ice thickness (SIT) climate data record (CDR) produced by the European Space Agency (ESA) Sea Ice Climate Change Initiative Phase 2 (CCI-2) is the longest available, Arctic-wide SIT record covering the period from 2002 to 2017. SIT data are based on radar altimetry measurements of sea ice freeboard from the Environmental Satellite (ENVISAT) and CryoSat-2 (CS2). The CCI-2 SIT has previously been validated with in situ observations from drilling, airborne remote sensing, electromagnetic (EM) measurements and upward-looking sonars (ULSs) from multiple ice-covered regions of the Arctic. Here we present the Laptev Sea CCI-2 SIT record from 2002 to 2017 and use newly acquired ULS and upward-looking acoustic Doppler current profiler (ADCP) sea ice draft (VAL) data for validation of the gridded CCI-2 and additional satellite SIT products. The ULS and ADCP time series provide the first long-term satellite SIT validation data set from this important source region of sea ice in the Transpolar Drift. The comparison of VAL sea ice draft data with gridded monthly mean and orbit trajectory CCI-2 data, as well as merged CryoSat-2–SMOS (CS2SMOS) sea ice draft, shows that the agreement between the satellite and VAL draft data strongly depends on the thickness of the sampled ice. Rather than providing mean sea ice draft, the considered satellite products provide modal sea ice draft in the Laptev Sea. Ice drafts thinner than 0.7 m are overestimated, while drafts thicker than approximately 1.3 m are increasingly underestimated by all satellite products investigated for this study. The tendency of the satellite SIT products to better agree with modal sea ice draft and underestimate thicker ice needs to be considered for all past and future investigations into SIT changes in this important region. The performance of the CCI-2 SIT CDR is considered stable over time; however, observed trends in gridded CCI-2 SIT are strongly influenced by the uncertainties of ENVISAT and CS2 and the comparably short investigation period.

Description: There has been a steady increase in interest in mining of deep-sea minerals in the Clarion–Clipperton Zone (CCZ) in the eastern Pacific Ocean during the last decade. This region is known to be one of the most eddy-rich regions in the world ocean. Typically, mesoscale eddies are generated by intense wind bursts channeled through gaps in the Sierra Madre mountains in Central America. Here, we use a combination of satellite and in situ observations to evaluate the relationship between deep-sea current variability in the region of potential future mining and eddy kinetic energy (EKE) in the vicinity of gap winds. A geometry-based eddy detection algorithm has been applied to altimetry sea surface height data for a period of 24 years, from 1993 to 2016, in order to analyze the main characteristic parameters and the spatiotemporal variability of mesoscale eddies in the northeast tropical Pacific Ocean (NETP). Significant differences between the characteristics of eddies with different polarity (cyclonic vs. anticyclonic) were found. For eddies with lifetimes longer than 1 d, cyclonic polarity is more common than anticyclonic rotation. However, anticyclonic eddies are larger in size, show stronger vorticity, and survive longer in the ocean than cyclonic eddies (often 90 d or more). Besides the polarity of eddies, the location of eddy formation should be taken into consideration when investigating the impacted deep-ocean region as we found eddies originating from the Tehuantepec (TT) gap winds lasting longer in the ocean and traveling farther distances in a different direction compared to eddies produced by the Papagayo (PP) gap winds. Long-lived anticyclonic eddies generated by the TT gap winds are observed to travel distances up to 4500 km offshore, i.e., as far as west of 110∘ W. EKE anomalies observed in the surface of the central ocean at distances of ca. 2500 km from the coast correlate with the seasonal variability of EKE in the region of the TT gap winds with a time lag of 5–6 months. A significant seasonal variability of deep-ocean current velocities at water depths of 4100 m was observed in multiple-year time series data, likely reflecting the energy transfer of the surface EKE generated by the gap winds to the deep ocean. Furthermore, the influence of mesoscale eddies on deep-ocean currents is examined by analyzing the deep-ocean current measurements when an anticyclonic eddy crosses the study region. Our findings suggest that despite the significant modulation of dominant current directions driven by the bottom-reaching eddy, the current magnitude intensification was not strong enough to trigger local sediment resuspension in this region. A better insight into the annual variability of ocean surface mesoscale activity in the CCZ and its effects on deep-ocean current variability can be of great help to mitigate the impact of future potential deep-sea mining activities on the benthic ecosystem. On an interannual scale, a significant relationship between cyclonic eddy characteristics and El Niño–Southern Oscillation (ENSO) was found, whereas a weaker correlation was detected for anticyclonic eddies.

Description: Gravitationally consistent solutions of the Sea Level Equation from leakage‐corrected monthly‐mean GFZ RL06 Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow‐On (GRACE‐FO) Stokes coefficients reveal that barystatic sea level averaged over the whole global ocean was rising by 1.72 mm a−1 during the period April 2002 until August 2016. This rate refers to a truely global ocean averaging domain that includes all polar and semienclosed seas. The result corresponds to 2.02 mm a−1 mean barystatic sea level rise in the open ocean with a 1,000 km coastal buffer zone as obtained from a direct spatial integration of monthly GRACE data. The bias of +0.3 mm a−1 is caused by below‐average barystatic sea level rise in close proximity to coastal mass losses induced by the smaller gravitational attraction of the remaining continental ice and water masses. Alternative spherical harmonics solutions from CSR, JPL, and TU Graz reveal open‐ocean rates between 1.94 and 2.08 mm a−1, thereby demonstrating that systematic differences among the processing centers are much reduced in the latest release. We introduce in this paper a new method to approximate spatial leakage from the differences of two differently filtered global gravity fields. A globally constant and time‐invariant scale factor required to obtain full leakage from those filter differences is found to be 3.9 for GFZ RL06 when filtered with DDK3, and lies between 3.9 and 4.4 for other processing centers. Spatial leakage is estimated for every month in terms of global grids, thereby providing also valuable information of intrabasin leakage that is potentially relevant for hydrologic and hydrometeorological applications.

Description: Understanding and quantifying the global methane (CH4) budget is important for assessing realistic pathways to mitigate climate change. Atmospheric emissions and concentrations of CH4 continue to increase, making CH4 the second most important human-influenced greenhouse gas in terms of climate forcing, after carbon dioxide (CO2). The relative importance of CH4 compared to CO2 depends on its shorter atmospheric lifetime, stronger warming potential, and variations in atmospheric growth rate over the past decade, the causes of which are still debated. Two major challenges in reducing uncertainties in the atmospheric growth rate arise from the variety of geographically overlapping CH4 sources and from the destruction of CH4 by short-lived hydroxyl radicals (OH). To address these challenges, we have established a consortium of multidisciplinary scientists under the umbrella of the Global Carbon Project to synthesize and stimulate new research aimed at improving and regularly updating the global methane budget. Following Saunois et al. (2016), we present here the second version of the living review paper dedicated to the decadal methane budget, integrating results of top-down studies (atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up estimates (including process-based models for estimating land surface emissions and atmospheric chemistry, inventories of anthropogenic emissions, and data-driven extrapolations). For the 2008–2017 decade, global methane emissions are estimated by atmospheric inversions (a top-down approach) to be 576 Tg CH4 yr−1 (range 550–594, corresponding to the minimum and maximum estimates of the model ensemble). Of this total, 359 Tg CH4 yr−1 or ∼ 60 % is attributed to anthropogenic sources, that is emissions caused by direct human activity (i.e. anthropogenic emissions; range 336–376 Tg CH4 yr−1 or 50 %–65 %). The mean annual total emission for the new decade (2008–2017) is 29 Tg CH4 yr−1 larger than our estimate for the previous decade (2000–2009), and 24 Tg CH4 yr−1 larger than the one reported in the previous budget for 2003–2012 (Saunois et al., 2016). Since 2012, global CH4 emissions have been tracking the warmest scenarios assessed by the Intergovernmental Panel on Climate Change. Bottom-up methods suggest almost 30 % larger global emissions (737 Tg CH4 yr−1, range 594–881) than top-down inversion methods. Indeed, bottom-up estimates for natural sources such as natural wetlands, other inland water systems, and geological sources are higher than top-down estimates. The atmospheric constraints on the top-down budget suggest that at least some of these bottom-up emissions are overestimated. The latitudinal distribution of atmospheric observation-based emissions indicates a predominance of tropical emissions (∼ 65 % of the global budget, 〈 30∘ N) compared to mid-latitudes (∼ 30 %, 30–60∘ N) and high northern latitudes (∼ 4 %, 60–90∘ N). The most important source of uncertainty in the methane budget is attributable to natural emissions, especially those from wetlands and other inland waters. Some of our global source estimates are smaller than those in previously published budgets (Saunois et al., 2016; Kirschke et al., 2013). In particular wetland emissions are about 35 Tg CH4 yr−1 lower due to improved partition wetlands and other inland waters. Emissions from geological sources and wild animals are also found to be smaller by 7 Tg CH4 yr−1 by 8 Tg CH4 yr−1, respectively. However, the overall discrepancy between bottom-up and top-down estimates has been reduced by only 5 % compared to Saunois et al. (2016), due to a higher estimate of emissions from inland waters, highlighting the need for more detailed research on emissions factors. Priorities for improving the methane budget include (i) a global, high-resolution map of water-saturated soils and inundated areas emitting methane based on a robust classification of different types of emitting habitats; (ii) further development of process-based models for inland-water emissions; (iii) intensification of methane observations at local scales (e.g., FLUXNET-CH4 measurements) and urban-scale monitoring to constrain bottom-up land surface models, and at regional scales (surface networks and satellites) to constrain atmospheric inversions; (iv) improvements of transport models and the representation of photochemical sinks in top-down inversions; and (v) development of a 3D variational inversion system using isotopic and/or co-emitted species such as ethane to improve source partitioning.

Description: Reconstructions of global hydroclimate during the Common Era (CE; the past ∼2000 years) are important for providing context for current and future global environmental change. Stable isotope ratios in water are quantitative indicators of hydroclimate on regional to global scales, and these signals are encoded in a wide range of natural geologic archives. Here we present the Iso2k database, a global compilation of previously published datasets from a variety of natural archives that record the stable oxygen (δ18O) or hydrogen (δ2H) isotopic compositions of environmental waters, which reflect hydroclimate changes over the CE. The Iso2k database contains 759 isotope records from the terrestrial and marine realms, including glacier and ground ice (210); speleothems (68); corals, sclerosponges, and mollusks (143); wood (81); lake sediments and other terrestrial sediments (e.g., loess) (158); and marine sediments (99). Individual datasets have temporal resolutions ranging from sub-annual to centennial and include chronological data where available. A fundamental feature of the database is its comprehensive metadata, which will assist both experts and nonexperts in the interpretation of each record and in data synthesis. Key metadata fields have standardized vocabularies to facilitate comparisons across diverse archives and with climate-model-simulated fields. This is the first global-scale collection of water isotope proxy records from multiple types of geological and biological archives. It is suitable for evaluating hydroclimate processes through time and space using large-scale synthesis, model–data intercomparison and (paleo)data assimilation. The Iso2k database is available for download at https://doi.org/10.25921/57j8-vs18 (Konecky and McKay, 2020) and is also accessible via the NOAA/WDS Paleo Data landing page: https://www.ncdc.noaa.gov/paleo/study/29593 (last access: 30 July 2020).

Description: We use observations from novel biogeochemical profiling floats deployed by the Southern Ocean Carbon and Climate Observations and Modeling program to estimate annual net community production (ANCP; associated with carbon export) from the seasonal drawdown of mesopelagic oxygen and surface nitrate in the Southern Ocean. Our estimates agree with previous observations in showing an increase in ANCP in the vicinity of the polar front (∼3 mol C m−2 y−1), compared to lower rates in the subtropical zone (≤ 1 mol C m−2 y−1) and the seasonal ice zone (〈2 mol C m−2 y−1). Paradoxically, the increase in ANCP south of the subtropical front is associated with elevated surface nitrate and silicate concentrations, but decreasing surface iron. We hypothesize that iron limitation promotes silicification in diatoms, which is evidenced by the low silicate to nitrate ratio of surface waters around the Antarctic polar front. High diatom silicification increases the ballasting effect of particulate organic carbon and overall ANCP in this region. A model-based assessment of our methods shows a good agreement between ANCP estimates based on oxygen and nitrate drawdown and the modeled downward organic carbon flux at 100 m. This agreement supports the presumption that net biological consumption is the dominant process affecting the drawdown of these chemical tracers and that, given sufficient data, ANCP can be inferred from observations of oxygen and/or nitrate drawdown in the Southern Ocean.

Description: With the mining of polymetallic nodules from the deep-sea seafloor once more evoking commercial interest, decisions must be taken on how to most efficiently regulate and monitor physical and community disturbance in these remote ecosystems. Image-based approaches allow non-destructive assessment of the abundance of larger fauna to be derived from survey data, with repeat surveys of areas possible to allow time series data collection. At the time of writing, key underwater imaging platforms commonly used to map seafloor fauna abundances are autonomous underwater vehicles (AUVs), remotely operated vehicles (ROVs) and towed camera “ocean floor observation systems” (OFOSs). These systems are highly customisable, with cameras, illumination sources and deployment protocols changing rapidly, even during a survey cruise. In this study, eight image datasets were collected from a discrete area of polymetallic-nodule-rich seafloor by an AUV and several OFOSs deployed at various altitudes above the seafloor. A fauna identification catalogue was used by five annotators to estimate the abundances of 20 fauna categories from the different datasets. Results show that, for many categories of megafauna, differences in image resolution greatly influenced the estimations of fauna abundance determined by the annotators. This is an important finding for the development of future monitoring legislation for these areas. When and if commercial exploitation of these marine resources commences, robust and verifiable standards which incorporate developing technological advances in camera-based monitoring surveys should be key to developing appropriate management regulations for these regions.

Description: Sediments in oxygen-depleted marine environments can be an important sink or source of bio-essential trace metals in the ocean. However, the key mechanisms controlling the release from or burial of trace metals in sediments are not exactly understood. Here, we investigate the benthic biogeochemical cycling of Fe and Cd in the oxygen minimum zone off Peru. We combine bottom water profiles, pore water profiles, as well as benthic fluxes determined from pore water profiles and in-situ from benthic chamber incubations along a depth transect at 12° S. In agreement with previous studies, both concentration-depth profiles and in-situ benthic fluxes indicate a Fe release from sediments into bottom waters. Diffusive Fe fluxes and Fe fluxes from benthic chamber incubations are roughly consistent (0.3–17.1 mmol m−2 y−1), indicating that diffusion is the main transport mechanism of dissolved Fe across the sediment-water interface. The occurrence of mats of sulfur oxidizing bacteria on the seafloor represents an important control on the spatial distribution of Fe fluxes by regulating hydrogen sulfide (H2S) concentrations and, potentially, Fe sulfide precipitation within the surface sediment. Removal of dissolved Fe after its release to anoxic bottom waters is rapid in the first 4 m away from the seafloor (half-life 〈 3 min) which hints to oxidative removal by nitrite or interaction with particles in the benthic boundary layer. Benthic flux estimates of Cd are indicative of a flux into the sediment within the oxygen minimum zone. Fluxes from benthic chamber incubations (up to 22.6 µmol m−2 y−1) exceed the diffusive fluxes (〈 1 µmol m−2 y−1) by a factor 〉 25, indicating that downward diffusion of Cd across the sediment-water interface is of subordinate importance for Cd removal from benthic chambers. As Cd removal in benthic chambers co-varies with H2S concentrations in the pore water of surface sediments, we argue that Cd removal is mediated by precipitation of CdS within the chamber. A mass balance approach, taking into account the contributions of diffusive fluxes and fluxes measured in benthic chambers as well as Cd delivery with organic material suggests that CdS precipitation in the near-bottom water could make an important contribution to the overall Cd mass accumulation in the sediment solid phase. According to our results, the solubility of trace metal sulfide minerals (Cd 〈〈 Fe) is a key-factor controlling trace metal removal and consequently the magnitude as well as the temporal and spatial heterogeneity of sedimentary fluxes. We argue that depending on their sulfide solubility, sedimentary source or sink fluxes of trace metals will change differentially as a result of declining oxygen concentrations and an associated expansion of sulfidic surface sediments. Such a trend could cause a change in the trace metal stoichiometry of upwelling water masses with potential consequences for marine ecosystems in the surface ocean.

Description: Oceanic Anoxic Events (OAEs) document major perturbations of the global carbon cycle with repercussions on the Earth’s climate and ocean circulation that are relevant to understand future climate trends. Here, we compare sedimentation patterns, nutrient cycling, organic carbon accumulation and carbon isotope variability across Cretaceous Oceanic Anoxic Events OAE1a and OAE2 in two drill cores with unusually high sedimentation rates from the Vocontian Basin (southern France) and Tarfaya Basin (southern Morocco). OAE1a and OAE2 exhibit remarkable similarities in the evolution of their δ13C excursion with long-lasting negative carbon isotope excursions preceding the onset of both anoxic events, supporting the view that OAEs were triggered by massive emissions of volcanic CO2 into the atmosphere. Based on analysis of cyclic sediment variations, we estimated the duration of the individual phases within the carbon isotope excursions. For both events, we identify: (1) a precursor phase lasting ~ 430 kyr and ~ 130 kyr, (2) an onset phase of ~ 390 and ~ 70 kyr, (3) a peak phase of ~ 600 and ~ 90 kyr, (4) a plateau phase of ~ 1400 and ~ 200 kyr and (5) a recovery phase of ~ 630 and ~ 440 kyr, respectively. The total duration of the positive carbon isotope excursion is estimated as 3400 kyr and 790 kyr and that of the main carbon accumulation phase as 980 kyr and 180 kyr, for OAE1a and OAE 2 respectively. The extended duration of the peak, plateau and recovery phases requires fundamental changes in global nutrient cycles either (1) through excess nutrient inputs to the oceans by increasing continental weathering and river discharge or (2) through nutrient-recycling from the marine sediment reservoir. We investigated the role of phosphorus on the development of carbon accumulation by analysing phosphorus speciation across OAE2 and the mid-Cenomanian Event (MCE) in the Tarfaya Basin. The ratios of organic carbon and total nitrogen to reactive phosphorus (Corg/Preact and Ntotal/Preact) prior to OAE2 and the MCE hover close to or below the Redfield ratio characteristic of marine organic matter. Decreases in reactive phosphorus resulting in Corg/Preact and Ntotal/Preact above the Redfield ratio during the later phase of OAE2 and the MCE indicate leakage from the sedimentary column into the water column under the influence of intensified and expanded oxygen minimum zones. These results suggest that a positive feedback loop, rooted in the benthic phosphorus cycle, contributed to increased marine productivity and carbon burial over an extended period of time during OAEs.

Description: Controlled manipulation of environmental conditions within large enclosures in the ocean, so-called pelagic mesocosms, has become a standard method to explore potential responses of marine plankton communities to anthropogenic change. Among the challenges of interpreting mesocosm data is the often uncertain role of vertical mixing, which usually is not observed directly. To account for mixing nonetheless, two pragmatic assumptions are common: either that the water column is homogeneously mixed or that it is divided into two water bodies with a horizontal barrier inhibiting turbulent exchange. In this study, we present a model-based reanalysis of vertical turbulent diffusion in the mesocosm experiments PeECE III and KOSMOS 2013. Our diffusivity estimates indicate intermittent mixing events along with stagnating periods and yield simulated temperature and salinity profiles that are consistent with the observations. Here, we provide the respective diffusivities as a comprehensive data product in the Network Common Data Format (NetCDF). This data product will help to guide forthcoming model studies that aim at deepening our understanding of biogeochemical processes in the PeECE III and KOSMOS 2013 mesocosms, such as the CO2-related changes in marine carbon export. In addition, we make our model code available, providing an adjustable tool to simulate vertical mixing in any other pelagic mesocosm. The data product and the model code are available at https://doi.org/10.1594/PANGAEA.905311 (Mathesius et al., 2019).