ALBERT

Description: The Denmark Strait Overflow (DSO) contributes roughly half to the total volume transport of the Nordic overflows. The overflow increases its volume by entraining ambient water as it descends into the subpolar North Atlantic, feeding into the deep branch of the Atlantic Meridional Overturning Circulation. In June 2012, a multiplatform experiment was carried out in the DSO plume on the continental slope off Greenland (180 km downstream of the sill in Denmark Strait), to observe the variability associated with the entrainment of ambient waters into the DSO plume. In this study, we report on two high-dissipation events captured by an autonomous underwater vehicle (AUV) by horizontal profiling in the interfacial layer between the DSO plume and the ambient water. Strong dissipation of turbulent kinetic energy of O( math formula) W kg−1 was associated with enhanced small-scale temperature variance at wavelengths between 0.05 and 500 m as deduced from a fast-response thermistor. Isotherm displacement slope spectra reveal a wave number-dependence characteristic of turbulence in the inertial-convective subrange ( math formula) at wavelengths between 0.14 and 100 m. The first event captured by the AUV was transient, and occurred near the edge of a bottom-intensified energetic eddy. Our observations imply that both horizontal advection of warm water and vertical mixing of it into the plume are eddy-driven and go hand in hand in entraining ambient water into the DSO plume. The second event was found to be a stationary feature on the upstream side of a topographic elevation located in the plume pathway. Flow-topography interaction is suggested to drive the intense mixing at this site.

Description: The ocean's potential to export carbon to depth partly depends on the fraction of primary production (PP) sinking out of the euphotic zone (i.e., the e-ratio). Measurements of PP and export flux are often performed simultaneously in the field, although there is a temporal delay between those parameters. Thus, resulting e-ratio estimates often incorrectly assume an instantaneous downward export of PP to export flux. Evaluating results from four mesocosm studies, we find that peaks in organic matter sedimentation lag chlorophyll a peaks by 2 to 15 days. We discuss the implications of these time lags (TLs) for current e-ratio estimates and evaluate potential controls of TL. Our analysis reveals a strong correlation between TL and the duration of chlorophyll a buildup, indicating a dependency of TL on plankton food web dynamics. This study is one step further toward time-corrected e-ratio estimates

Description: Stable isotope compositions can potentially be used to trace atmospheric Cd inputs to the surface ocean and anthropogenic Cd emissions to the atmosphere. Both of these applications may provide valuable insights into the effects of anthropogenic activities on the cycling of Cd in the environment. However, a lack of constraints for the Cd isotope compositions of atmospheric aerosols is currently hindering such studies. Here, we present stable Cd isotope data for aerosols collected over the Tropical Atlantic Ocean. The samples feature variable proportions of mineral dust-derived and anthropogenic Cd, yet exhibit similar isotope compositions, thus negating the distinction of these Cd sources using isotopic signatures in this region. Isotopic variability between these two atmospheric Cd sources may be identified in other areas, and thus warrants further investigation. Regardless, these data provide important initial constraints on the isotope composition of atmospheric Cd inputs to the ocean.

Description: The upstream sources and pathways of the Denmark Strait Overflow Water and their variability have been investigated using a high-resolution model hindcast. This global simulation covers the period from 1948 to 2009 and uses a fine model mesh (1/20°) to resolve mesoscale features and the complex current structure north of Iceland explicitly. The three sources of the Denmark Strait Overflow, the shelfbreak East Greenland Current (EGC), the separated EGC, and the North Icelandic Jet, have been analyzed using Eulerian and Lagrangian diagnostics. The shelfbreak EGC contributes the largest fraction in terms of volume and freshwater transport to the Denmark Strait Overflow and is the main driver of the overflow variability. The North Icelandic Jet contributes the densest water to the Denmark Strait Overflow and shows only small temporal transport variations. During summer, the net volume and freshwater transports to the south are reduced. On interannual time scales, these transports are highly correlated with the large-scale wind stress curl around Iceland and, to some extent, influenced by the North Atlantic Oscillation, with enhanced southward transports during positive phases. The Lagrangian trajectories support the existence of a hypothesized overturning loop along the shelfbreak north of Iceland, where water carried by the North Icelandic Irminger Current is transformed and feeds the North Icelandic Jet. Monitoring these two currents and the region north of the Iceland shelfbreak could provide the potential to track long-term changes in the Denmark Strait Overflow and thus also the AMOC.

Description: The North Atlantic Current (NAC) is subject to variability on multiannual to decadal time scales, influencing the transport of volume, heat, and freshwater from the subtropical to the eastern subpolar North Atlantic (NA). Current observational time series are either too short or too episodic to study the processes involved. Here we compare the observed continuous NAC transport time series at the western flank of the Mid-Atlantic Ridge (MAR) and repeat hydrographic measurements at the OVIDE line in the eastern Atlantic with the NAC transport and circulation in the high-resolution (1/20°) ocean model configuration VIKING20 (1960–2008). The modeled baroclinic NAC transport relative to 3400 m (24.5 ± 7.1 Sv) at the MAR is only slightly lower than the observed baroclinic mean of 27.4 ± 4.7 Sv from 1993 to 2008, and extends further north by about 0.5°. In the eastern Atlantic, the western NAC (WNAC) carries the bulk of the transport in the model, while transport estimates based on hydrographic measurements from five repeated sections point to a preference for the eastern NAC (ENAC). The model is able to simulate the main features of the subpolar NA, providing confidence to use the model output to analyze the influence of the North Atlantic Oscillation (NAO). Model based velocity composites reveal an enhanced NAC transport across the MAR of up to 6.7 Sv during positive NAO phases. Most of that signal (5.4 Sv) is added to the ENAC transport, while the transport of the WNAC was independent of the NAO.

Description: While the Earth's surface has considerably warmed over the past two decades, the tropical Pacific has featured a cooling of sea surface temperatures in its eastern and central parts, which went along with an unprecedented strengthening of the equatorial trade winds, the surface component of the Pacific Walker Circulation (PWC). Previous studies show that this decadal trend in the trade winds is generally beyond the range of decadal trends simulated by climate models when forced by historical radiative forcing. There is still a debate on the origin of and the potential role that internal variability may have played in the recent decadal surface wind trend. Using a number of long control (unforced) integrations of global climate models and several observational data sets, we address the question as to whether the recent decadal to multidecadal trends are robustly classified as an unusual event or the persistent response to external forcing. The observed trends in the tropical Pacific surface climate are still within the range of the long-term internal variability spanned by the models but represent an extreme realization of this variability. Thus, the recent observed decadal trends in the tropical Pacific, though highly unusual, could be of natural origin. We note that the long-term trends in the selected PWC indices exhibit a large observational uncertainty, even hindering definitive statements about the sign of the trends.

Description: We examine the mean pathways, transit timescales, and transformation of waters flowing from the Pacific and the marginal seas through the Indian Ocean (IO) on their way toward the South Atlantic within a high-resolution ocean/sea-ice model. The model fields are analyzed from a Lagrangian perspective where water volumes are tracked as they enter the IO. The IO contributes 12.6 Sv to Agulhas leakage, which within the model is 14.1 ± 2.2 Sv, the rest originates from the South Atlantic. The Indonesian Through-flow constitutes about half of the IO contribution, is surface bound, cools and salinificates as it leaves the basin within 10–30 years. Waters entering the IO south of Australia are at intermediate depths and maintain their temperature-salinity properties as they exit the basin within 15–35 years. Of these waters, the contribution from Tasman leakage is 1.4 Sv. The rest stem from recirculation from the frontal regions of the Southern Ocean. The marginal seas export 1.0 Sv into the Atlantic within 15–40 years, and the waters cool and freshen on-route. However, the model's simulation of waters from the Gulfs of Aden and Oman are too light and hence overly influenced by upper ocean circulations. In the Cape Basin, Agulhas leakage is well mixed. On-route, temperature-salinity transformations occur predominantly in the Arabian Sea and within the greater Agulhas Current region. Overall, the IO exports at least 7.9 Sv from the Pacific to the Atlantic, thereby quantifying the strength of the upper cell of the global conveyor belt.

Description: The potential of mining seafloor massive sulfide deposits for metals such as Cu, Zn, and Au is currently debated. One key challenge is to predict where the largest deposits worth mining might form, which in turn requires understanding the pattern of subseafloor hydrothermal mass and energy transport. Numerical models of heat and fluid flow are applied to illustrate the important role of fault zone properties (permeability and width) in controlling mass accumulation at hydrothermal vents at slow spreading ridges. We combine modeled mass-flow rates, vent temperatures, and vent field dimensions with the known fluid chemistry at the fault-controlled Logatchev 1 hydrothermal field of the Mid-Atlantic Ridge. We predict that the 135 kilotons of SMS at this site (estimated by other studies) can have accumulated with a minimum depositional efficiency of 5% in the known duration of hydrothermal venting (58,200 year age of the deposit). In general, the most productive faults must provide an efficient fluid pathway while at the same time limit cooling due to mixing with entrained cold seawater. This balance is best met by faults that are just wide and permeable enough to control a hydrothermal plume rising through the oceanic crust. Model runs with increased basal heat input, mimicking a heat flow contribution from along-axis, lead to higher mass fluxes and vent temperatures, capable of significantly higher SMS accumulation rates. Nonsteady state conditions, such as the influence of a cooling magmatic intrusion beneath the fault zone, also can temporarily increase the mass flux while sustaining high vent temperatures.

Description: Diazotrophic dinitrogen (N2) fixation contributes ~76% to "new" nitrogen inputs to the sunlit open ocean, but environmental factors determining N2 fixation rates are not well constrained. Excess phosphate (phosphate-nitrate/16 〉 0) and iron availability control N2 fixation rates in the eastern tropical North Atlantic (ETNA), but it remains an open question how excess phosphate is generated within or supplied to the phosphate-depleted sunlit layer. Our observations in the ETNA region (8°N-15°N, 19°W-23°W) suggest that Prochlorococcus and Synechococcus, the two ubiquitous non-diazotrophic cyanobacteria with cellular N:P ratios higher than the Redfield ratio, create an environment of excess phosphate, which cannot be explained by diapycnal mixing, atmospheric, and riverine inputs. Thus, our results unveil a new biogeochemical niche construction mechanism by non-diazotrophic cyanobacteria for their diazotrophic phylum group members (N2 fixers). Our observations may help to understand the prevalence of diazotrophy in low-phosphate, oligotrophic regions.

Description: We present bathymetry and hydrological observations collected in the summer of 2014 from two fjord systems in southeastern Greenland with a multibeam sonar system. Our results provide a detailed bathymetric map of the fjord complex around the island of Skjoldungen in Skjoldungen Fjord and the outer part of Timmiarmiut Fjord and show far greater depths compared to the International Bathymetric Chart of the Arctic Ocean. The hydrography collected shows different properties in the fjords with the bottom water masses below 240 m in Timmiarmiut Fjord being 1–2 °C warmer than in the two fjords around Skjoldungen, but data also illustrate the influence of sills on the exchange of deeper water masses within fjords. Moreover, evidence of subglacial discharge in Timmiarmiut Fjord, which is consistent with satellite observations of ice mélange set into motion, adds to our increasing understanding of the distribution of subglacial meltwater.

Description: Oceanographic observations from the Eurasian Basin north of Svalbard collected between January and June 2015 from the N-ICE2015 drifting expedition are presented. The unique winter observations are a key contribution to existing climatologies of the Arctic Ocean, and show a ∼100 m deep winter mixed layer likely due to high sea ice growth rates in local leads. Current observations for the upper ∼200 m show mostly a barotropic flow, enhanced over the shallow Yermak Plateau. The two branches of inflowing Atlantic Water are partly captured, confirming that the outer Yermak Branch follows the perimeter of the plateau, and the inner Svalbard Branch the coast. Atlantic Water observed to be warmer and shallower than in the climatology, is found directly below the mixed layer down to 800 m depth, and is warmest along the slope, while its properties inside the basin are quite homogeneous. From late May onwards, the drift was continually close to the ice edge and a thinner surface mixed layer and shallower Atlantic Water coincided with significant sea ice melt being observed.

Description: Our study followed the seasonal cycling of soluble (SFe), colloidal (CFe), dissolved (DFe), total dissolvable (TDFe), labile particulate (LPFe) and total particulate (TPFe) iron in the Celtic Sea (NE Atlantic Ocean). Preferential uptake of SFe occurred during the spring bloom, preceding the removal of CFe. Uptake and export of Fe during the spring bloom, coupled with a reduction in vertical exchange, led to Fe deplete surface waters (〈0.2 nM DFe; 0.11 nM LPFe, 0.45 nM TDFe, 1.84 nM TPFe) during summer stratification. Below the seasonal thermocline, DFe concentrations increased from spring to autumn, mirroring NO3- and consistent with supply from remineralised sinking organic material, and cycled independently of particulate Fe over seasonal timescales. These results demonstrate that summer Fe availability is comparable to the seasonally Fe limited Ross Sea shelf, and therefore is likely low enough to affect phytoplankton growth and species composition.

Description: Current climate models disagree on how much carbon dioxide land ecosystems take up for photosynthesis. Tracking the stronger carbonyl sulfide signal could help.

Description: Oceanic dimethyl sulfide (DMS) is of interest due to its critical influence on atmospheric sulfur compounds in the marine atmosphere and its hypothesized significant role in global climate. High-resolution shipboard underway measurements of surface seawater DMS and the partial pressure of carbon dioxide (pCO2) were conducted in the Atlantic Ocean and Indian Ocean sectors of the Southern Ocean (SO), the southeast Indian Ocean, and the northwest Pacific Ocean from February to April 2014 during the 30th Chinese Antarctic Research Expedition. The SO, particularly in the region south of 58°S, had the highest mean surface seawater DMS concentration of 4.1 ± 8.3 nM (ranged from 0.1 to 73.2 nM) and lowest mean seawater pCO2 level of 337 ± 50 μatm (ranged from 221 to 411 μatm) over the entire cruise. Significant variations of surface seawater DMS and pCO2 in the seasonal ice zone (SIZ) of SO were observed, which are mainly controlled by biological process and sea ice activity. We found a significant negative relationship between DMS and pCO2 in the SO SIZ using 0.1° resolution, [DMS] seawater = -0.160 [pCO2] seawater + 61.3 (r2 = 0.594, n = 924, p 〈 0.001). We anticipate that the relationship may possibly be utilized to reconstruct the surface seawater DMS climatology in the SO SIZ. Further studies are necessary to improve the universality of this approach. Key Points: • The characteristics of surface water DMS and pCO2 distributions from the Southern Ocean to northwest Pacific Ocean are investigated • The correlations between DMS, pCO2, and environmental parameters are analyzed • Anticorrelation between DMS and pCO2 is found in the seasonal ice zone of the Southern Ocean

Description: Ocean acidification resulting from the uptake of anthropogenic carbon dioxide (CO2) by the ocean is considered a major threat to marine ecosystems. Here we examined the effects of ocean acidification on microbial community dynamics in the eastern Baltic Sea during the summer of 2012 when inorganic nitrogen and phosphorus were strongly depleted. Large-volume in situ mesocosms were employed to mimic present, future and far future CO2 scenarios. All six groups of phytoplankton enumerated by flow cytometry ( 〈  20 µm cell diameter) showed distinct trends in net growth and abundance with CO2 enrichment. The picoeukaryotic phytoplankton groups Pico-I and Pico-II displayed enhanced abundances, whilst Pico-III, Synechococcus and the nanoeukaryotic phytoplankton groups were negatively affected by elevated fugacity of CO2 (fCO2). Specifically, the numerically dominant eukaryote, Pico-I, demonstrated increases in gross growth rate with increasing fCO2 sufficient to double its abundance. The dynamics of the prokaryote community closely followed trends in total algal biomass despite differential effects of fCO2 on algal groups. Similarly, viral abundances corresponded to prokaryotic host population dynamics. Viral lysis and grazing were both important in controlling microbial abundances. Overall our results point to a shift, with increasing fCO2, towards a more regenerative system with production dominated by small picoeukaryotic phytoplankton.

Description: Natural gas hydrates are considered a potential resource for gas production on industrial scales. Gas hydrates contribute to the strength and stiffness of the hydrate-bearing sediments. During gas production, the geomechanical stability of the sediment is compromised. Due to the potential geotechnical risks and process management issues, the mechanical behavior of the gas hydrate-bearing sediments needs to be carefully considered. In this study, we describe a coupling concept that simplifies the mathematical description of the complex interactions occurring during gas production by isolating the effects of sediment deformation and hydrate phase changes. Central to this coupling concept is the assumption that the soil grains form the load-bearing solid skeleton, while the gas hydrate enhances the mechanical properties of this skeleton. We focus on testing this coupling concept in capturing the overall impact of geomechanics on gas production behavior though numerical simulation of a high-pressure isotropic compression experiment combined with methane hydrate formation and dissociation. We consider a linear-elastic stress-strain relationship because it is uniquely defined and easy to calibrate. Since, in reality, the geomechanical response of the hydrate-bearing sediment is typically inelastic and is characterized by a significant shear-volumetric coupling, we control the experiment very carefully in order to keep the sample deformations small and well within the assumptions of poroelasticity. The closely coordinated experimental and numerical procedures enable us to validate the proposed simplified geomechanics-to-flow coupling, and set an important precursor toward enhancing our coupled hydro-geomechanical hydrate reservoir simulator with more suitable elastoplastic constitutive models.

Description: The Summer East Atlantic (SEA) mode is the second dominant mode of summer low-frequency variability in the Euro-Atlantic region. Using reanalysis data, we show that SEA-related circulation anomalies significantly influence temperatures and precipitation over Europe. We present evidence that part of the interannual SEA variability is forced by diabatic heating anomalies of opposing signs in the tropical Pacific and Caribbean that induce an extratropical Rossby wave train. This precipitation dipole is related to SST anomalies characteristic of the developing ENSO phases. Seasonal hindcast experiments forced with observed sea surface temperatures (SST) exhibit skill at capturing the interannual SEA variability corroborating the proposed mechanism and highlighting the possibility for improved prediction of boreal summer variability. Our results indicate that tropical forcing of the SEA likely played a role in the dynamics of the 2015 European heat wave.

Description: Biogeochemical models, capturing the major feedbacks of the pelagic ecosystem of the world ocean, are today often embedded into Earth system models which are increasingly used for decision making regarding climate policies. These models contain poorly constrained parameters (e.g., maximum phytoplankton growth rate), which are typically adjusted until the model shows reasonable behavior. Systematic approaches determine these parameters by minimizing the misfit between the model and observational data. In most common model approaches, however, the underlying functions mimicking the biogeochemical processes are nonlinear and non-convex. Thus, systematic optimization algorithms are likely to get trapped in local minima and might lead to non-optimal results. To judge the quality of an obtained parameter estimate, we propose determining a preferably large lower bound for the global optimum that is relatively easy to obtain and that will help to assess the quality of an optimum, generated by an optimization algorithm. Due to the unavoidable noise component in all observations, such a lower bound is typically larger than zero. We suggest deriving such lower bounds based on typical properties of biogeochemical models (e.g., a limited number of extremes and a bounded time derivative). We illustrate the applicability of the method with two real-world examples. The first example uses real-world observations of the Baltic Sea in a box model setup. The second example considers a three-dimensional coupled ocean circulation model in combination with satellite chlorophyll a.

Description: Very short-lived substances (VSLS) contribute as source gases significantly to the tropospheric and stratospheric bromine loading. At present, an estimated 25% of stratospheric bromine is of oceanic origin. In this study, we investigate how climate change may impact the ocean- atmosphere flux of brominated VSLS, their atmospheric transport, and chemical transformations and evaluate how these changes will affect stratospheric ozone over the 21st century. Under the assumption of fixed ocean water concentrations and RCP6.0 scenario, we find an increase of the ocean- atmosphere flux of brominated VSLS of about 8-10% by the end of the 21st century compared to present day. A decrease in the tropospheric mixing ratios of VSLS and an increase in the lower stratosphere are attributed to changes in atmospheric chemistry and transport. Our model simulations reveal that this increase is counteracted by a corresponding reduction of inorganic bromine. Therefore the total amount of bromine from VSLS in the stratosphere will not be changed by an increase in upwelling. Part of the increase of VSLS in the tropical lower stratosphere results from an increase in the corresponding tropopause height. As the depletion of stratospheric ozone due to bromine depends also on the availability of chlorine, we find the impact of bromine on stratospheric ozone at the end of the 21st century reduced compared to present day. Thus, these studies highlight the different factors influencing the role of brominated VSLS in a future climate

Description: Reanalysis data sets are widely used to understand atmospheric processes and past variability, and are often used to stand in as "observations" for comparisons with climate model output. Because of the central role of water vapor (WV) and ozone (O3) in climate change, it is important to understand how accurately and consistently these species are represented in existing global reanalyses. In this paper, we present the results of WV and O3 intercomparisons that have been performed as part of the SPARC (Stratosphere–troposphere Processes and their Role in Climate) Reanalysis Intercomparison Project (S-RIP). The comparisons cover a range of timescales and evaluate both inter-reanalysis and observation-reanalysis differences. We also provide a systematic documentation of the treatment of WV and O3 in current reanalyses to aid future research and guide the interpretation of differences amongst reanalysis fields. The assimilation of total column ozone (TCO) observations in newer reanalyses results in realistic representations of TCO in reanalyses except when data coverage is lacking, such as during polar night. The vertical distribution of ozone is also relatively well represented in the stratosphere in reanalyses, particularly given the relatively weak constraints on ozone vertical structure provided by most assimilated observations and the simplistic representations of ozone photochemical processes in most of the reanalysis forecast models. However, significant biases in the vertical distribution of ozone are found in the upper troposphere and lower stratosphere in all reanalyses. In contrast to O3, reanalysis estimates of stratospheric WV are not directly constrained by assimilated data. Observations of atmospheric humidity are typically used only in the troposphere, below a specified vertical level at or near the tropopause. The fidelity of reanalysis stratospheric WV products is therefore mainly dependent on the reanalyses' representation of the physical drivers that influence stratospheric WV, such as temperatures in the tropical tropopause layer, methane oxidation, and the stratospheric overturning circulation. The lack of assimilated observations and known deficiencies in the representation of stratospheric transport in reanalyses result in much poorer agreement amongst observational and reanalysis estimates of stratospheric WV. Hence, stratospheric WV products from the current generation of reanalyses should generally not be used in scientific studies.

Description: The Sumatran subduction zone exhibits strong seismic and tsunamogenic potential with the prominent examples of the 2004, 2005 and 2007 earthquakes. Here, we invert travel time data of local earthquakes for vp and vp/vs velocity models of the central Sumatran forearc. Data were acquired by an amphibious seismometer network consisting of 52 land stations and 10 ocean bottom seismometers located on a segment of the Sumatran subduction zone that had not ruptured in a great earthquake since 1797 but witnessed recent ruptures to the north in 2005 (Nias earthquake, Mw = 8.7) and to the south in 2007 (Bengkulu earthquake, Mw = 8.5). 2D and 3D vp velocity anomalies reveal the downgoing slab and the sedimentary basins. Although the seismicity pattern in the study area appears to be strongly influenced by the obliquely subducting Investigator Fracture Zone to at least 200 km depth, the 3D velocity model shows prevailing trench parallel structures at depths of the plate interface. The tomographic model suggests a thinned crust below the basin east of the forearc islands (Nias, Pulau Batu, Siberut) at ~ 180 km distance to the trench. Vp velocities beneath the magmatic arc and the Sumatran fault zone SFZ are around 5 km/s at 10 km depth and the vp/vs ratios in the uppermost 10 km are low, indicating the presence of felsic lithologies typical for continental crust. We find moderately elevated vp/vs values of 1.85 at ~ 150 km distance to the trench in the region of the Mentawai fault. Vp/vs ratios suggest absence of large scale alteration of the mantle wedge and might explain why the seismogenic plate interface (observed as a locked zone from geodetic data) extends below the continental forearc Moho in Sumatra. Reduced vp velocities beneath the forearc basin covering the region between Mentawai Islands and the Sumatra mainland possibly reflect a reduced thickness of the overriding crust.

Description: The analysis of high-resolution vector magnetic data acquired by deep-sea submersibles (DSSs) requires the development of specific approaches adapted to their uneven tracks. We present a method that takes advantage of (1) the varying altitude of the DSS above the seafloor and (2) high-resolution multibeam bathymetric data acquired separately, at higher altitude, by an Autonomous Underwater Vehicle, to estimate the absolute magnetization intensity and the magnetic polarity of the shallow subseafloor along the DSS path. We apply this method to data collected by DSS Nautile on a small active basalt-hosted hydrothermal site. The site is associated with a lack of magnetization, in agreement with previous findings at the same kind of sites: the contrast between nonmagnetic sulfide deposits/stockwork zone and strongly magnetized basalt is sufficient to explain the magnetic signal observed at such a low altitude. Both normal and reversed polarities are observed in the lava flows surrounding the site, suggesting complex history of accumulating volcanic flows.

Description: In this study ship- and AUV-based multibeam data from the German Mn-nodule license area in the Clarion-Clipperton Zone (CCZ; eastern Pacific) are linked to ground truth data from optical imaging. Photographs obtained by an AUV enable semi-quantitative assessments of nodule coverage at a spatial resolution in the range of meters. Together with high resolution AUV bathymetry this revealed a correlation of small-scale terrain variations (〈 5 m horizontally, 〈 1 m vertically) with nodule abundance. In the presented data set, increased nodule coverage could be correlated with slopes 〉 1.8° and concave terrain. On a more regional scale, factors such as the geological setting (existence of horst and graben structures, sediment thickness, outcropping basement) and influence of bottom currents seem to play an essential role for the spatial variation of nodule abundance and the related hard substrate habitat. AUV imagery was also successfully employed to map the distribution of re-settled sediment following a disturbance and sediment cloud generation during a sampling deployment of an Epibenthic Sledge. Data from before and after the "disturbance" allows a direct assessment of the impact. Automated image processing analyzed the nodule coverage at the seafloor, revealing nodule blanketing by resettling of suspended sediment within 16 hours after the disturbance. The visually detectable impact was spatially limited to a maximum of 100m distance from the disturbance track, downstream of the bottom water current. A correlation with high resolution AUV bathymetry reveals that the blanketing pattern varies in extent by tens of meters, strictly following the bathymetry, even in areas of only slightly undulating seafloor (〈 1 m vertical change). These results highlight the importance of detailed terrain knowledge when engaging in resource assessment studies for nodule abundance estimates and defining minable areas. At the same time, it shows the importance of high resolution mapping for detailed benthic habitat studies that show a heterogeneity at scales of 10 m to 100 m. Terrain knowledge is also needed to determine the scale of the impact by seafloor sediment blanketing during mining-operations.

Description: To maintain the chance of keeping the average global temperature increase below 2 degrees C and to limit long-term climate change, removing carbon dioxide from the atmosphere (carbon dioxide removal, CDR) is becoming increasingly necessary. We analyze optimal and cost-effective climate policies in the dynamic integrated assessment model (IAM) of climate and the economy (DICE2016R) and investigate (1) the utilization of (ocean) CDR under different climate objectives, (2) the sensitivity of policies with respect to carbon cycle feedbacks, and (3) how well carbon cycle feedbacks are captured in the carbon cycle models used in state-of-the-art IAMs. Overall, the carbon cycle model in DICE2016R shows clear improvements compared to its predecessor, DICE2013R, capturing much better long-term dynamics and also oceanic carbon outgassing due to excess oceanic storage of carbon from CDR. However, this comes at the cost of a (too) tight short-term remaining emission budget, limiting the model suitability to analyze low-emission scenarios accurately. With DICE2016R, the compliance with the 2 degrees C goal is no longer feasible without negative emissions via CDR. Overall, the optimal amount of CDR has to take into account (1) the emission substitution effect and (2) compensation for carbon cycle feedbacks.

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: The region encompassing the Kuroshio Extension (KE) in the Northwestern Pacific Ocean (25°N–45°N and 130°E–180°E) is one of the most eddy-energetic regions of the global ocean. The three-dimensional structures and transports of mesoscale eddies in this region are comprehensively investigated by combined use of satellite data and Argo profiles. With the allocation of Argo profiles inside detected eddies, the spatial variations of structures of eddy temperature and salinity anomalies are analyzed. The results show that eddies predominantly have subsurface (near-surface) intensified temperature and salinity anomalies south (north) of the KE jet, which is related to different background stratifications between these regions. A new method based on eddy trajectories and the inferred three-dimensional eddy structures is proposed to estimate heat and salt transports by eddy movements in a Lagrangian framework. Spatial distributions of eddy transports are presented over the vicinity of the KE for the first time. The magnitude of eddy-induced meridional heat (freshwater volume) transport is on the order of 0.01 PW (103 m3/s). The eddy heat transport divergence results in an oceanic heat loss south and heat gain north of the KE, thereby reinforcing and counteracting the oceanic heat loss from air-sea fluxes south and north of the KE jet, respectively. It also suggests a poleward heat transport across the KE jet due to eddy propagation.

Description: During the summer monsoon, the western tropical Indian Ocean is predicted to be a hot spot for dimethylsulfide emissions, the major marine sulfur source to the atmosphere, and an important aerosol precursor. Other aerosol relevant fluxes, such as isoprene and sea spray, should also be enhanced, due to the steady strong winds during the monsoon. Marine air masses dominate the area during the summer monsoon, excluding the influence of continentally derived pollutants. During the SO234-2/235 cruise in the western tropical Indian Ocean from July to August 2014, directly measured eddy covariance DMS fluxes confirm that the area is a large source of sulfur to the atmosphere (cruise average 9.1 μmol m−2 d−1). The directly measured fluxes, as well as computed isoprene and sea spray fluxes, were combined with FLEXPART backward and forward trajectories to track the emissions in space and time. The fluxes show a significant positive correlation with aerosol data from the Terra and Suomi-NPP satellites, indicating a local influence of marine emissions on atmospheric aerosol numbers.

Description: Climate-driven changes in environmental conditions have significant and complex effects on marine ecosystems. Variability in phytoplankton elements and biochemicals can be important for global ocean biogeochemistry and ecological functions, while there is currently limited understanding on how elements and biochemicals respond to the changing environments in key coccolithophore species such as Emiliania huxleyi. We investigated responses of elemental stoichiometry and fatty acids (FAs) in a strain of E. huxleyi under 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). Overall, C : N : P stoichiometry showed the most pronounced response to N : P supply ratios, with high ratios of particulate organic carbon vs. particulate organic nitrogen (POC : PON) and low ratios of PON vs. particulate organic phosphorus (PON : POP) in low-N media, and high POC : POP and PON : POP in low-P media. The ratio of particulate inorganic carbon vs. POC (PIC : POC) and polyunsaturated fatty acid proportions strongly responded to temperature and pCO2, both being lower under high pCO2 and higher with warming. We observed synergistic interactions between warming and nutrient deficiency (and high pCO2) on elemental cellular contents and docosahexaenoic acid (DHA) proportion in most cases, indicating the enhanced effect of warming under nutrient deficiency (and high pCO2). Our results suggest differential sensitivity of elements and FAs to the changes in temperature, nutrient availability and pCO2 in E. huxleyi, which is to some extent unique compared to non-calcifying algal classes. Thus, simultaneous changes of elements and FAs should be considered when predicting future roles of E. huxleyi in the biotic-mediated connection between biogeochemical cycles, ecological functions and climate change.

Description: The pre-industrial millennium is among the periods selected by the Paleoclimate Model Intercomparison Project (PMIP) for experiments contributing to the sixth phase of the Coupled Model Intercomparison Project (CMIP6) and the fourth phase of the PMIP (PMIP4). The past1000 transient simulations serve to investigate the response to (mainly) natural forcing under background conditions not too different from today, and to discriminate between forced and internally generated variability on interannual to centennial timescales. This paper describes the motivation and the experimental set-ups for the PMIP4-CMIP6 past1000 simulations, and discusses the forcing agents orbital, solar, volcanic, and land use/land cover changes, and variations in greenhouse gas concentrations. The past1000 simulations covering the pre-industrial millennium from 850 Common Era (CE) to 1849 CE have to be complemented by historical simulations (1850 to 2014 CE) following the CMIP6 protocol. The external forcings for the past1000 experiments have been adapted to provide a seamless transition across these time periods. Protocols for the past1000 simulations have been divided into three tiers. A default forcing data set has been defined for the Tier 1 (the CMIP6 past1000) experiment. However, the PMIP community has maintained the flexibility to conduct coordinated sensitivity experiments to explore uncertainty in forcing reconstructions as well as parameter uncertainty in dedicated Tier 2 simulations. Additional experiments (Tier 3) are defined to foster collaborative model experiments focusing on the early instrumental period and to extend the temporal range and the scope of the simulations. This paper outlines current and future research foci and common analyses for collaborative work between the PMIP and the observational communities (reconstructions, instrumental data).

Description: Although mesoscale ocean eddies are ubiquitous in the Southern Ocean, their average regional and seasonal association with phytoplankton has not been quantified systematically yet. To this end, we identify over 100000 mesoscale eddies with diameters of 50km and more in the Southern Ocean and determine the associated phytoplankton biomass anomalies using satellite-based chlorophyll-a (Chl) as a proxy. The mean Chl anomalies, δChl, associated with these eddies, comprising the upper echelon of the oceanic mesoscale, exceed ±10% over wide regions. The structure of these anomalies is largely zonal, with cyclonic, thermocline lifted, eddies having positive anomalies in the subtropical waters north of the Antarctic Circumpolar Current (ACC) and negative anomalies along its main flow path. The pattern is similar, but reversed for anticyclonic, thermocline deepened eddies. The seasonality of δChl is weak in subtropical waters, but pronounced along the ACC, featuring a seasonal sign switch. The spatial structure and seasonality of the mesoscale δChl can be explained largely by lateral advection, especially local eddy-stirring. A prominent exception is the ACC region in winter, where δChl is consistent with a modulation of phytoplankton light exposure caused by an eddy-induced modification of the mixed layer depth. The clear impact of mesoscale eddies on phytoplankton may implicate a downstream effect on Southern Ocean biogeochemical properties, such as mode water nutrient contents.

Description: Anthropogenic activities have resulted in enhanced lead (Pb) emissions to the environment over the past century, mainly through the combustion of leaded gasoline. Here, we present the first combined dissolved (DPb), labile (LpPb) and particulate (PPb) Pb dataset from the Northeast Atlantic (Celtic Sea) since the phasing out of leaded gasoline in Europe. Concentrations of DPb in surface waters have decreased by 4-fold over the last four decades. We demonstrate that anthropogenic Pb is transported from the Mediterranean Sea over long distances (〉2500 km). Benthic DPb fluxes exceeded the atmospheric Pb flux in the region, indicating the importance of sediments as a contemporary Pb source. A strong positive correlation between DPb, PPb and LpPb indicates a dynamic equilibrium between the phases and the potential for particles to ‘buffer’ the DPb pool. This study provides insights into Pb biogeochemical cycling and demonstrates the potential of Pb in constraining ocean circulation patterns.

Description: The variability of the Atlantic Meridional Overturning Circulation (AMOC) may play a role in sea surface temperature predictions on seasonal to decadal time scales. Therefore, AMOC seasonal cycles are a potential baseline for interpreting predictions. Here we present estimates for the seasonal cycle of transports of volume, temperature, and freshwater associated with the upper limb of the AMOC in the eastern subpolar North Atlantic on the Extended Ellett Line hydrographic section between Scotland and Iceland. Due to weather, ship‐based observations are primarily in summer. Recent glider observations during other seasons present an opportunity to investigate the seasonal variability in the upper layer of the AMOC. First, we document a new method to quality control and merge ship, float, and glider hydrographic observations. This method accounts for the different spatial sampling rates of the three platforms. The merged observations are used to compute seasonal cycles of volume, temperature, and freshwater transports in the Rockall Trough. These estimates are similar to the seasonal cycles in two eddy‐resolving ocean models. Volume transport appears to be the primary factor modulating other Rockall Trough transports. Finally, we show that the weakest transports occur in summer, consistent with seasonal changes in the regional‐scale wind stress curl. Although the seasonal cycle is weak compared to other variability in this region, the amplitude of the seasonal cycle in the Rockall Trough, roughly 0.5–1 Sv about a mean of 3.4 Sv, may account for up to 7–14% of the heat flux between Scotland and Greenland.

Description: Benthic foraminifera from Bottsand coastal lagoon, western Baltic Sea, have been studied since the mid-1960s. They were monitored annually in late autumn since 2003 at the terminal ditch of the lagoon. There were 12 different species recognised, of which three have not been recorded during earlier investigations. Dominant species showed strong interannual fluctuations and a steady increase in population densities over the last decade. Elphidium incertum, a stenohaline species of the Baltic deep water fauna, colonised the Bottsand lagoon in 2016, most likely during a period of salinities 〉19 units and water temperatures of 18 °C on average in early autumn. The high salinities probably triggered their germination from a propagule bank in the ditch bottom sediment. The new E. incertum population showed densities higher by an order of magnitude than those of the indigenous species. The latter did not decline, revealing that E. incertum used another food source or occupied a different microhabitat. Elphidium incertum survived transient periods of lower salinities in late autumn 2017, though with reduced abundances, and became a regular faunal constituent at the Bottsand lagoon.

Description: The quantitative reconstruction of past seawater salinity has yet to be achieved and the search for a direct and independent salinity proxy is ongoing. Recent culture and field studies show a significant positive correlation of Na/Ca with salinity in benthic and planktonic foraminiferal calcite. For accurate paleoceanographic reconstructions, consistent and reliable calibrations are necessary, which are still missing. In order to assess the reliability of foraminiferal Na/Ca as a direct proxy for seawater salinity, this study presents electron microprobe Na/Ca data, measured on cultured specimens of Trilobatus sacculifer. The culture experiments were conducted over a wide salinity range of 26 to 45, while temperature was kept constant. To further understand potential controlling factors of Na incorporation, measurements were also performed on foraminifera cultured at various temperatures in the range of 19.5 °C to 29.5 °C under constant salinity conditions. Foraminiferal Na/Ca ratios positively correlate with seawater salinity (Na/Caforam = 0.97 + 0.115 ⋅ Salinity, R = 0.97, p 〈 0.005). Temperature on the other hand exhibits no statistically significant relationship with Na/Ca ratios indicating salinity to be the dominant factor controlling Na incorporation. The culturing results are corroborated by measurements on T. sacculifer from Caribbean and Gulf of Guinea surface sediments. In conclusion, planktonic foraminiferal Na/Ca can be applied as a reliable proxy for reconstructing sea surface salinities, albeit species-specific calibrations might be necessary.

Description: We report here the results of total mercury (HgT) determinations along the 2014 Geotraces Geovide cruise (GA01 transect) in the North Atlantic Ocean (NA) from Lisbon (Portugal) to the coast of Labrador (Canada). HgT concentrations in unfiltered samples (HgTUNF) were log-normally distributed and ranged between 0.16 and 1.54 pmol L−1, with a geometric mean of 0.51 pmol L−1 for the 535 samples analysed. The dissolved fraction (〈 0.45 µm) of HgT (HgTF), determined on 141 samples, averaged 78 % of the HgTUNF for the entire data set, 84 % for open seawaters (below 100 m) and 91 % if the Labrador Sea data are excluded, where the primary production was high (with a winter convection down to 1400 m). HgTUNF concentrations increased eastwards and with depth from Greenland to Europe and from subsurface to bottom waters. The HgTUNF concentrations were similarly low in the subpolar gyre waters ( ∼  0.45 pmol L−1), whereas they exceeded 0.60 pmol L−1 in the subtropical gyre waters. The HgTUNF distribution mirrored that of dissolved oxygen concentration, with highest concentration levels associated with oxygen-depleted zones. The relationship between HgTF and the apparent oxygen utilization confirms the nutrient-like behaviour of Hg in the NA. An extended optimum multiparameter analysis allowed us to characterize HgTUNF concentrations in the different source water types (SWTs) present along the transect. The distribution pattern of HgTUNF, modelled by the mixing of SWTs, show Hg enrichment in Mediterranean waters and North East Atlantic Deep Water and low concentrations in young waters formed in the subpolar gyre and Nordic seas. The change in anthropogenic Hg concentrations in the Labrador Sea Water during its eastward journey suggests a continuous decrease in Hg content in this water mass over the last decades. Calculation of the water transport driven by the Atlantic Meridional Overturning Circulation across the Portugal–Greenland transect indicates northward Hg transport within the upper limb and southward Hg transport within the lower limb, with resulting net northward transport of about 97.2 kmol yr−1.

Description: The chemistry–climate model ECHAM-HAMMOZ contains a detailed representation of tropospheric and stratospheric reactive chemistry and state-of-the-art parameterizations of aerosols using either a modal scheme (M7) or a bin scheme (SALSA). This article describes and evaluates the model version ECHAM6.3-HAM2.3-MOZ1.0 with a focus on the tropospheric gas-phase chemistry. A 10-year model simulation was performed to test the stability of the model and provide data for its evaluation. The comparison to observations concentrates on the year 2008 and includes total column observations of ozone and CO from IASI and OMI, Aura MLS observations of temperature, HNO3, ClO, and O3 for the evaluation of polar stratospheric processes, an ozonesonde climatology, surface ozone observations from the TOAR database, and surface CO data from the Global Atmosphere Watch network. Global budgets of ozone, OH, NOx, aerosols, clouds, and radiation are analyzed and compared to the literature. ECHAM-HAMMOZ performs well in many aspects. However, in the base simulation, lightning NOx emissions are very low, and the impact of the heterogeneous reaction of HNO3 on dust and sea salt aerosol is too strong. Sensitivity simulations with increased lightning NOx or modified heterogeneous chemistry deteriorate the comparison with observations and yield excessively large ozone budget terms and too much OH. We hypothesize that this is an impact of potential issues with tropical convection in the ECHAM model.

Description: We reanalyze existing paleodata of global mean surface temperature ΔTg and radiative forcing ΔR of CO2 and land ice albedo for the last 800,000 years to show that a state‐dependency in paleoclimate sensitivity S, as previously suggested, is only found if ΔTg is based on reconstructions, and not when ΔTg is based on model simulations. Furthermore, during times of decreasing obliquity (periods of land‐ice sheet growth and sea level fall) the multi‐millennial component of reconstructed ΔTg diverges from CO2, while in simulations both variables vary more synchronously, suggesting that the differences during these times are due to relatively low rates of simulated land ice growth and associated cooling. To produce a reconstruction‐based extrapolation of S for the future we exclude intervals with strong ΔTg‐CO2 divergence and find that S is less state‐dependent, or even constant (state‐independent), yielding a mean equilibrium warming of 2–4 K for a doubling of CO2.

Description: Chlorophyll (Chl) is a distinctive component of autotrophic organisms, often used as an indicator of phytoplankton biomass in the ocean. However, assessment of phytoplankton biomass from Chl relies on the accurate estimation of the Chl:carbon(C) ratio. Here we present global patterns of Chl:C ratios in the surface ocean obtained from a phytoplankton growth model that accounts for the optimal acclimation of phytoplankton to ambient nutrient, light, and temperature conditions. The model agrees largely with observed/expected global patterns of Chl:C. Combining our Chl:C estimates with satellite Chl and particulate organic carbon (POC), we infer phytoplankton C concentration in the surface ocean and its contribution to the total POC pool. Our results suggest that the portion of POC corresponding to living phytoplankton is higher in subtropical latitudes and less productive regions (∼30–70%) and decreases to ∼10–30% toward high latitudes and productive regions. An important caveat of our model is the lack of iron limiting effects on phytoplankton physiology. Comparison of our predicted phytoplankton biomass with an independent estimate of total POC reveals a positive correlation between nitrate concentrations and nonphotosynthetic POC in the surface ocean. This correlation disappears when a constant Chl:C is applied. Our analysis is not constrained by assumptions of constant Chl:C or phytoplankton:POC ratio, providing a novel independent analysis of phytoplankton biomass in the surface ocean. These results highlight the importance of accounting for the variability in Chl:C and its application in distinguishing the autotrophic and heterotrophic components in the assemblage of the marine plankton ecosystem.

Description: Lakes and reservoirs play an important role in the carbon cycle, and therefore monitoring their metabolic rates is essential. The triple oxygen-isotope anomaly of dissolved O2 [17Δ =  ln(1+δ17O) − 0.518  ×  ln(1 + δ18O)] offers a new, in situ, perspective on primary production, yet little is known about 17Δ from freshwater systems. We investigated the 17Δ together with the oxygen : argon ratio [Δ(O2 ∕ Ar)] in the subtropical Feitsui Reservoir in Taiwan from June 2014 to July 2015. Here, we present the seasonal variations in 17Δ, GP (gross production), NP (net production) and the NP ∕ GP (net to gross ratio) in association with environmental parameters. The 17Δ varied with depth and season, with values ranging between 26 and 205 per meg. The GP rates were observed to be higher (702 ± 107 mg C m−2 d−1) in winter than those (303 ± 66 mg C m−2 d−1) recorded during the summer. The overall averaged GP was 220 g C m−2 yr−1 and NP was −3 g C m−2 yr−1, implying the reservoir was net heterotrophic on an annual basis. This is due to negative NP rates from October to February (−198 ± 78 mg C m−2 d−1). Comparisons between GP rates obtained from the isotope mass balance approach and 14C bottle incubation method (14C–GP) showed consistent values on the same order of magnitude with a GP ∕ 14C–GP ratio of 1.2 ± 1.1. Finally we noted that, although typhoon occurrences were scarce, higher than average 17Δ values and GP rates were recorded after typhoon events.

Description: The historical developments are reviewed that have led from a bottom-up responsibility initiative of concerned scientists to the emergence of a nationwide interdisciplinary Priority Program on the assessment of Climate Engineering (CE), funded by the German Research Foundation (DFG). Given the perceived lack of comprehensive and comparative appraisals of different CE methods, the Priority Program was designed to encompass both solar radiation management (SRM) and carbon dioxide removal (CDR) ideas, and to cover the atmospheric, terrestrial and oceanic realm. First key findings obtained by the ongoing Priority Program are summarized and reveal that compared to earlier assessments, such as the 2009 Royal Society report, more detailed investigations tend to indicate less efficiency, lower effectiveness and often lower safety. Emerging research trends are discussed in the context of the recent Paris agreement to limit global warming to less than two degrees and the associated increasing reliance on negative emission technologies. Our results show then when deployed at scales large enough to have a significant impact on atmospheric CO2, even CDR methods such as afforestation – often perceived as ‘benign’ – can have substantial side effects and may raise severe ethical, legal and governance issues. We suppose that before being deployed at climatically relevant scales, any negative-emission or climate engineering method will require careful analysis of efficiency, effectiveness and undesired side effects.

Description: Oxygen optode measurements on floats and gliders suffer from a slow time response and various sources of drift in the calibration coefficients. Based on two dual-O2 Argo floats, we show how to post-correct for the effect of the optode's time response and give an update on optode in situ drift stability and in-air calibration. Both floats are equipped with an unpumped Aanderaa 4330 optode and a pumped Sea-Bird SBE63 optode. Response times for the pumped SBE63 were derived following Bittig et al. (2014) and the same methods were used to correct the time response bias. Using both optodes on each float, the time response regime of the unpumped Aanderaa optode was characterized more accurately than previously possible. Response times for the pumped SBE63 on profiling floats are in the range of 25–40 s, while they are between 60 and 95 s for the unpumped 4330 optode. Our parameterization can be employed to post-correct the slow optode time response on floats and gliders. After correction, both sensors agree to within 2–3 µmol kg−1 (median difference) in the strongest gradients (120 µmol kg−1 change over 8 min or 20 dbar) and better elsewhere. However, time response correction is only possible if measurement times are known, i.e., provided by the platform as well as transmitted and stored with the data. The O2 in-air measurements show a significant in situ optode drift of −0.40 and −0.27 % yr−1 over the available 2 and 3 years of deployment, respectively. Optode in-air measurements are systematically biased high during midday surfacings compared to dusk, dawn, and nighttime. While preference can be given to nighttime surfacings to avoid this in-air calibration bias, we suggest a parameterization of the daytime effect as a function of the Sun's elevation to be able to use all data and to better constrain the result. Taking all effects into account, calibration factors have an uncertainty of 0.1 %. In addition, in-air calibration factors vary by 0.1–0.2 % when using different reanalysis models as a reference. The overall accuracy that can be achieved following the proposed correction routines is better than 1 µmol kg−1.

Description: Selecting appropriate indicators is essential to aggregate the information provided by climate model outputs into a manageable set of relevant metrics on which assessments of climate engineering (CE) can be based. From all the variables potentially available from climate models, indicators need to be selected that are able to inform scientists and society on the development of the Earth system under CE, as well as on possible impacts and side effects of various ways of deploying CE or not. However, the indicators used so far have been largely identical to those used in climate change assessments and do not visibly reflect the fact that indicators for assessing CE (and thus the metrics composed of these indicators) may be different from those used to assess global warming. Until now, there has been little dedicated effort to identifying specific indicators and metrics for assessing CE. We here propose that such an effort should be facilitated by a more decision-oriented approach and an iterative procedure in close interaction between academia, decision makers, and stakeholders. Specifically, synergies and trade-offs between social objectives reflected by individual indicators, as well as decision-relevant uncertainties should be considered in the development of metrics, so that society can take informed decisions about climate policy measures under the impression of the options available, their likely effects and side effects, and the quality of the underlying knowledge base.

Description: Conventional integration of Earth system and ocean models can accrue considerable computational expenses, particularly for marine biogeochemical applications. "Offline" numerical schemes in which only the biogeochemical tracers are time stepped and transported using a pre-computed circulation field can substantially reduce the burden and are thus an attractive alternative. One such scheme is the "transport matrix method" (TMM), which represents tracer transport as a sequence of sparse matrix–vector products that can be performed efficiently on distributed-memory computers. While the TMM has been used for a variety of geochemical and biogeochemical studies, to date the resulting solutions have not been comprehensively assessed against their "online" counterparts. Here, we present a detailed comparison of the two. It is based on simulations of the state-of-the-art biogeochemical sub-model embedded within the widely used coarse-resolution University of Victoria Earth System Climate Model (UVic ESCM). The default, non-linear advection scheme was first replaced with a linear, third-order upwind-biased advection scheme to satisfy the linearity requirement of the TMM. Transport matrices were extracted from an equilibrium run of the physical model and subsequently used to integrate the biogeochemical model offline to equilibrium. The identical biogeochemical model was also run online. Our simulations show that offline integration introduces some bias to biogeochemical quantities through the omission of the polar filtering used in UVic ESCM and in the offline application of time-dependent forcing fields, with high latitudes showing the largest differences with respect to the online model. Differences in other regions and in the seasonality of nutrients and phytoplankton distributions are found to be relatively minor, giving confidence that the TMM is a reliable tool for offline integration of complex biogeochemical models. Moreover, while UVic ESCM is a serial code, the TMM can be run on a parallel machine with no change to the underlying biogeochemical code, thus providing orders of magnitude speed-up over the online model.

Description: Mesoscale eddies are ubiquitous features in the Southern Ocean, yet their phenomenology is not well quantified. To tackle this task, we use satellite observations of sea level anomalies and sea surface temperature (SST) as well as in situ temperature and salinity measurements from profiling floats. Over the period 1997–2010, we identified over a million mesoscale eddy instances and were able to track about 105 of them over 1 month or more. The Antarctic Circumpolar Current (ACC), the boundary current systems, and the regions where they interact are hot spots of eddy presence, representing also the birth places and graveyards of most eddies. These hot spots contrast strongly to areas shallower than about 2000 m, where mesoscale eddies are essentially absent, likely due to topographical steering. Anticyclones tend to dominate the southern subtropical gyres, and cyclones the northern flank of the ACC. Major causes of regional polarity dominance are larger formation numbers and lifespans, with a contribution of differential propagation pathways of long-lived eddies. Areas of dominance of one polarity are generally congruent with the same polarity being longer-lived, bigger, of larger amplitude, and more intense. Eddies extend down to at least 2000 m. In the ACC, eddies show near surface temperature and salinity maxima, whereas eddies in the subtropical areas generally have deeper anomaly maxima, presumably inherited from their origin in the boundary currents. The temperature and salinity signatures of the average eddy suggest that their tracer anomalies are a result of both trapping in the eddy core and stirring.

Description: The assessment of the ocean biota's role in climate climate change is often carried out with global biogeochemical ocean models that contain many components, and involve a high level of parametric uncertainty. Examination the models' fit to climatologies of inorganic tracers, after the models have been spun up to steady state, is a common, but computationally expensive procedure to assess model performance and reliability. Using new tools that have become available for global model assessment and calibration in steady state, this paper examines two different model types – a complex seven-component model (MOPS), and a very simple two-component model (RetroMOPS) – for their fit to dissolved quantities. Before comparing the models, a subset of their biogeochemical parameters has been optimised against annual mean nutrients and oxygen. Both model types fit the observations almost equally well. The simple model, which contains only nutrients and dissolved organic phosphorus (DOP), is sensitive to the parameterisation of DOP production and decay. The spatio-temporal decoupling of nitrogen and oxygen, and processes involved in their uptake and release, renders oxygen and nitrate valuable tracers for model calibration. In addition, the non-conservative nature of these tracers (with respect to their upper boundary condition) introduces the global bias as a useful additional constraint on model parameters. Dissolved organic phosphorous at the surface behaves antagonistically to phosphate, and suggests that observations of this tracer – although difficult to measure – may be an important asset for model calibration

Description: Benthic microbial methanogenesis is a known source of methane in marine systems. In most sediments, the majority of methanogenesis is located below the sulfate-reducing zone, as sulfate reducers outcompete methanogens for the major substrates hydrogen and acetate. The coexistence of methanogenesis and sulfate reduction has been shown before and is possible through the usage of noncompetitive substrates by methanogens such as methanol or methylated amines. However, knowledge about the magnitude, seasonality, and environmental controls of this noncompetitive methane production is sparse. In the present study, the presence of methanogenesis within the sulfate reduction zone (SRZ methanogenesis) was investigated in sediments (0–30 cm below seafloor, cm b.s.f.) of the seasonally hypoxic Eckernförde Bay in the southwestern Baltic Sea. Water column parameters such as oxygen, temperature, and salinity together with porewater geochemistry and benthic methanogenesis rates were determined in the sampling area "Boknis Eck" quarterly from March 2013 to September 2014 to investigate the effect of seasonal environmental changes on the rate and distribution of SRZ methanogenesis, to estimate its potential contribution to benthic methane emissions, and to identify the potential methanogenic groups responsible for SRZ methanogenesis. The metabolic pathway of methanogenesis in the presence or absence of sulfate reducers, which after the addition of a noncompetitive substrate was studied in four experimental setups: (1) unaltered sediment batch incubations (net methanogenesis), (2) 14C-bicarbonate labeling experiments (hydrogenotrophic methanogenesis), (3) manipulated experiments with the addition of either molybdate (sulfate reducer inhibitor), 2-bromoethanesulfonate (methanogen inhibitor), or methanol (noncompetitive substrate, potential methanogenesis), and (4) the addition of 13C-labeled methanol (potential methylotrophic methanogenesis). After incubation with methanol, molecular analyses were conducted to identify key functional methanogenic groups during methylotrophic methanogenesis. To also compare the magnitudes of SRZ methanogenesis with methanogenesis below the sulfate reduction zone (〉 30 cm b.s.f.), hydrogenotrophic methanogenesis was determined by 14C-bicarbonate radiotracer incubation in samples collected in September 2013. SRZ methanogenesis changed seasonally in the upper 30 cm b.s.f. with rates increasing from March (0.2 nmol cm−3 d−1) to November (1.3 nmol cm−3 d−1) 2013 and March (0.2 nmol cm−3 d−1) to September (0.4 nmol cm−3 d−1) 2014. Its magnitude and distribution appeared to be controlled by organic matter availability, C / N, temperature, and oxygen in the water column, revealing higher rates in the warm, stratified, hypoxic seasons (September–November) compared to the colder, oxygenated seasons (March–June) of each year. The majority of SRZ methanogenesis was likely driven by the usage of noncompetitive substrates (e.g., methanol and methylated compounds) to avoid competition with sulfate reducers, as was indicated by the 1000–3000-fold increase in potential methanogenesis activity observed after methanol addition. Accordingly, competitive hydrogenotrophic methanogenesis increased in the sediment only below the depth of sulfate penetration (〉 30 cm b.s.f.). Members of the family Methanosarcinaceae, which are known for methylotrophic methanogenesis, were detected by PCR using Methanosarcinaceae-specific primers and are likely to be responsible for the observed SRZ methanogenesis. The present study indicates that SRZ methanogenesis is an important component of the benthic methane budget and carbon cycling in Eckernförde Bay. Although its contributions to methane emissions from the sediment into the water column are probably minor, SRZ methanogenesis could directly feed into methane oxidation above the sulfate–methane transition zone.

Description: Over the past 17 years, the western boundary current system of the Labrador Sea has been closely observed by maintaining the 53°N observatory (moorings and shipboard station data) measuring the top-to-bottom flow field offshore from the Labrador shelf break. Volume transports for the North Atlantic Deep Water (NADW) components were calculated using different methods, including gap filling procedures for deployment periods with suboptimal instrument coverage. On average the Deep Western Boundary Current (DWBC) carries 30.2 ± 6.6 Sv of NADW southward, which are almost equally partitioned between Labrador Sea Water (LSW, 14.9 ± 3.9 Sv) and Lower North Atlantic Deep Water (LNADW, 15.3 ± 3.8 Sv). The transport variability ranges from days to decades, with the most prominent multiyear fluctuations at interannual to near decadal time scales (±5 Sv) in the LNADW overflow water mass. These long-term fluctuations appear to be in phase with the NAO-modulated wind fluctuations. The boundary current system off Labrador occurs as a conglomerate of nearly independent components, namely, the shallow Labrador Current, the weakly sheared LSW range, and the deep baroclinic, bottom-intensified current core of the LNADW, all of which are part of the cyclonic Labrador Sea circulation. This structure is relatively stable over time, and the 120 km wide boundary current is constrained seaward by a weak counterflow which reduces the deep water export by 10–15%.

Description: The balance in microbial net consumption of nitrogen and phosphorus was investigated in samples collected in two mesotrophic coastal environments: the Baltic Sea (Tvärminne field station) and the North Sea (Espegrend field station). For this, we have refined a bioassay based on the response in alkaline phosphatase activity (APA) over a matrix of combinations in nitrogen and phosphorus additions. This assay not only provides information on which element (N or P) is the primary limiting nutrient, but also gives a quantitative estimate for the excess of the secondary limiting element (P+ or N+, respectively), as well as the ratio of balanced net consumption of added N and P over short timescales (days). As expected for a Baltic Sea late spring–early summer situation, the Tvärminne assays (n =  5) indicated N limitation with an average P+ =  0.30 ± 0.10 µM-P, when incubated for 4 days. For short incubations (1–2 days), the Espegrend assays indicated P limitation, but the shape of the response surface changed with incubation time, resulting in a drift in parameter estimates toward N limitation. Extrapolating back to zero incubation time gave P limitation with N+ ≈  0.9 µM-N. The N : P ratio (molar) of nutrient net consumption varied considerably between investigated locations: from 2.3 ± 0.4 in the Tvärminne samples to 13 ± 5 and 32 ± 3 in two samples from Espegrend. Our assays included samples from mesocosm acidification experiments, but statistically significant effects of ocean acidification were not found by this method.

Description: Primary productivity is limited by the availability of nitrogen (N) in most of the coastal Arctic, as a large portion of N is released by the spring freshet and completely consumed during the following summer. Thus, understanding the fate of riverine nitrogen is critical to identify the link between dissolved nitrogen dynamic and coastal primary productivity to foresee upcoming changes in the Arctic seas, such as increase riverine discharge and permafrost thaw. Here, we provide a field-based study of nitrogen dynamic over the Laptev Sea shelf based on isotope geochemistry. We demonstrate that while most of the nitrate found under the surface fresh water layer is of remineralized origin, some of the nitrate originates from atmospheric input and was probably transported at depth by the mixing of brine-enriched denser water during sea-ice formation. Moreover, our results suggest that riverine dissolved organic nitrogen (DON) represents up to 6 times the total riverine release of nitrate and that about 62 to 76% of the DON is removed within the shelf waters. This is a crucial information regarding the near-future impact of climate change on primary productivity in the Eurasian coastal Arctic.

Description: Marine controlled source electromagnetic (CSEM) data have been collected to investigate methane seep sites and associated gas hydrate deposits at Opouawe Bank on the southern tip of the Hikurangi Margin, New Zealand. The bank is located in about 1000 m water depth within the gas hydrate stability field. The seep sites are characterized by active venting and typical methane seep fauna accompanied with patchy carbonate outcrops at the seafloor. Below the seeps, gas migration pathways reach from below the bottom-simulating reflector (at around 380 m sediment depth) toward the seafloor, indicating free gas transport into the shallow hydrate stability field. The CSEM data have been acquired with a seafloor-towed, electric multi-dipole system measuring the inline component of the electric field. CSEM data from three profiles have been analyzed by using 1-D and 2-D inversion techniques. High-resolution 2-D and 3-D multichannel seismic data have been collected in the same area. The electrical resistivity models show several zones of highly anomalous resistivities (〉50 Ωm) which correlate with high amplitude reflections located on top of narrow vertical gas conduits, indicating the coexistence of free gas and gas hydrates within the hydrate stability zone. Away from the seeps the CSEM models show normal background resistivities between ~1 and 2 Ωm. Archie's law has been applied to estimate gas/gas hydrate saturations below the seeps. At intermediate depths between 50 and 200 m below seafloor, saturations are between 40 and 80% and gas hydrate may be the dominating pore filling constituent. At shallow depths from 10 m to the seafloor, free gas dominates as seismic data and gas plumes suggest.

Description: The Ocean Model Intercomparison Project (OMIP) focuses on the physics and biogeochemistry of the ocean component of Earth system models participating in the sixth phase of the Coupled Model Intercomparison Project (CMIP6). OMIP aims to provide standard protocols and diagnostics for ocean models, while offering a forum to promote their common assessment and improvement. It also offers to compare solutions of the same ocean models when forced with reanalysis data (OMIP simulations) vs. when integrated within fully coupled Earth system models (CMIP6). Here we detail simulation protocols and diagnostics for OMIP's biogeochemical and inert chemical tracers. These passive-tracer simulations will be coupled to ocean circulation models, initialized with observational data or output from a model spin-up, and forced by repeating the 1948–2009 surface fluxes of heat, fresh water, and momentum. These so-called OMIP-BGC simulations include three inert chemical tracers (CFC-11, CFC-12, SF6) and biogeochemical tracers (e.g., dissolved inorganic carbon, carbon isotopes, alkalinity, nutrients, and oxygen). Modelers will use their preferred prognostic BGC model but should follow common guidelines for gas exchange and carbonate chemistry. Simulations include both natural and total carbon tracers. The required forced simulation (omip1) will be initialized with gridded observational climatologies. An optional forced simulation (omip1-spunup) will be initialized instead with BGC fields from a long model spin-up, preferably for 2000 years or more, and forced by repeating the same 62-year meteorological forcing. That optional run will also include abiotic tracers of total dissolved inorganic carbon and radiocarbon, CTabio and 14CTabio, to assess deep-ocean ventilation and distinguish the role of physics vs. biology. These simulations will be forced by observed atmospheric histories of the three inert gases and CO2 as well as carbon isotope ratios of CO2. OMIP-BGC simulation protocols are founded on those from previous phases of the Ocean Carbon-Cycle Model Intercomparison Project. They have been merged and updated to reflect improvements concerning gas exchange, carbonate chemistry, and new data for initial conditions and atmospheric gas histories. Code is provided to facilitate their implementation.

Description: The injection of sulphur into the stratosphere by explosive volcanic eruptions is the cause of significant climate variability. Based on sulphate records from a suite of ice cores from Greenland and Antarctica, the eVolv2k database includes estimates of the magnitudes and approximate source latitudes of major volcanic stratospheric sulphur injection (VSSI) events from 500 BCE to 1900 CE, constituting an update of prior reconstructions and an extension of the record by 1000 years. The VSSI estimates incorporate improvements to the ice core records in terms of synchronization and dating, refinements to the methods used to estimate VSSI from ice core records, and includes first estimates of the random uncertainties in VSSI values. VSSI estimates for many of the largest eruptions, including Samalas (1257), Tambora (1815) and Laki (1783) are within 10% of prior estimates. A number of strong events are included in eVolv2k which are largely underestimated or not included in earlier VSSI reconstructions, including events in 540, 574, 682 and 1108 CE. The long term annual mean VSSI from major volcanic eruptions is estimated to be ∼ 0.5 Tg [S] yr−1, ∼ 50 % greater than a prior reconstruction, due to the identification of more events and an increase in the magnitude of many intermediate events. A long-term, latitudinally and monthly resolved stratospheric aerosol optical depth (SAOD) time series is reconstructed from the eVolv2k VSSI estimates, and the resulting global mean SAOD is found to be similar (within 33%) to a prior reconstruction for most of the largest eruptions. The long-term (500 BCE–900 CE) average global mean SAOD estimated from the eVolv2k VSSI estimates and including a constant "background" injection of stratospheric sulphur is ∼ 0.014, 30 % greater than a prior reconstruction. These new long-term reconstructions of past VSSI and SAOD variability give context to recent volcanic forcing, suggesting that the 20th century was a period of somewhat weaker than average volcanic forcing, with current best estimates of 20th century mean VSSI and SAOD values being 25 and 14 % less, respectively, than the mean of the 500 BCE to 1900 CE period. The reconstructed VSSI and SAOD data are available at https://doi.org/10.1594/WDCC/eVolv2k_v2〉.

Description: The Cretaceous ocean witnessed intervals of profound perturbations such as volcanic input of large amounts of CO2, anoxia, eutrophication and introduction of biologically relevant metals. Some of these extreme events were characterized by size reduction and/or morphological changes of a few calcareous nannofossil species. The correspondence between intervals of high trace metal concentrations and coccolith dwarfism suggests a negative effect of these elements on nannoplankton biocalcification processes in past oceans. In order to test this hypothesis, we explored the potential effect of a mixture of trace metals on growth and morphology of four living coccolithophore species, namely Emiliania huxleyi, Gephyrocapsa oceanica, Pleurochrysis carterae and Coccolithus pelagicus. The phylogenetic history of coccolithophores shows that the selected living species are linked to Mesozoic species showing dwarfism under excess metal concentrations. The trace metals tested were chosen to simulate the environmental stress identified in the geological record and upon known trace metal interactions with living coccolithophore algae. Our laboratory experiments demonstrated that elevated trace metal concentrations, similarly to the fossil record, affect coccolithophore algae size and/or weight. Smaller coccoliths were detected in E. huxleyi and C. pelagicus, while coccoliths of G. oceanica showed a decrease in size only at the highest trace metal concentrations. P. carterae coccolith size was unresponsive to changing trace metal concentrations. These differences among species allow discriminating the most- (P. carterae), intermediate- (E. huxleyi and G. oceanica) and least-tolerant (C. pelagicus) taxa. The fossil record and the experimental results converge on a selective response of coccolithophores to metal availability. These species-specific differences must be considered before morphological features of coccoliths are used to reconstruct paleo-chemical conditions.

Description: The Baltic Sea is a seasonally ice-covered marginal sea in northern Europe with intense wintertime ship traffic and a sensitive ecosystem. Understanding and modeling the evolution of the sea-ice pack is important for climate effect studies and forecasting purposes. Here we present and evaluate the sea-ice component of a new NEMO–LIM3.6-based ocean–sea-ice setup for the North Sea and Baltic Sea region (NEMO-Nordic). The setup includes a new depth-based fast-ice parametrization for the Baltic Sea. The evaluation focuses on long-term statistics, from a 45-year long hindcast, although short-term daily performance is also briefly evaluated. We show that NEMO-Nordic is well suited for simulating the mean sea-ice extent, concentration, and thickness as compared to the best available observational data set. The variability of the annual maximum Baltic Sea ice extent is well in line with the observations, but the 1961–2006 trend is underestimated. Capturing the correct ice thickness distribution is more challenging. Based on the simulated ice thickness distribution we estimate the undeformed and deformed ice thickness and concentration in the Baltic Sea, which compares reasonably well with observations.

Description: Receiver functions (RF) have been used for several decades to study structures beneath seismic stations. Although most available stations are deployed on-shore, the number of ocean bottom station (OBS) experiments has increased in recent years. Almost all OBSs have to deal with higher noise levels and a limited deployment time (∼1 year), resulting in a small number of usable records of teleseismic earthquakes. Here, we use OBSs deployed as mid-aperture array in the deep ocean (4.5-5.5 km water depth) of the eastern mid-Atlantic. We use evaluation criteria for OBS data and beam forming to enhance the quality of the RFs. Although some stations show reverberations caused by sedimentary cover, we are able to identify the Moho signal, indicating a normal thickness (5-8 km) of oceanic crust. Observations at single stations with thin sediments (300-400 m) indicate that a probable sharp lithosphere-asthenosphere boundary (LAB) might exist at a depth of ∼70-80 km which is in line with LAB depth estimates for similar lithospheric ages in the Pacific. The mantle discontinuities at ∼410 km and ∼660 km are clearly identifiable. Their delay times are in agreement with PREM. Overall the usage of beam formed earthquake recordings for OBS RF analysis is an excellent way to increase the signal quality and the number of usable events.

Description: We developed a new method for the determination of dissolved nitric oxide (NO) in discrete seawater samples based on the combination of a purge-and-trap setup and a fluorometric detection of NO. 2,3-diaminonaphthalene (DAN) reacts with NO in seawater to form the highly fluorescent 2,3-naphthotriazole (NAT). The fluorescence intensity was linear for NO concentrations in the range from 0.14 to 19 nmol L−1. We determined a detection limit of 0.068 nmol L−1, an average recovery coefficient of 83.8 % (80.2–90.0 %), and a relative standard deviation of ±7.2 %. With our method we determined for the first time the temporal and spatial distributions of NO surface concentrations in coastal waters of the Yellow Sea off Qingdao and in Jiaozhou Bay during a cruise in November 2009. The concentrations of NO varied from below the detection limit to 0.50 nmol L−1 with an average of 0.26 ± 0.14 nmol L−1. NO surface concentrations were generally enhanced significantly during daytime, implying that NO formation processes such as NO2− photolysis are much higher during daytime than chemical NO consumption, which, in turn, lead to a significant decrease in NO concentrations during nighttime. In general, NO surface concentrations and measured NO production rates were higher compared to previously reported measurements. This might be caused by the high NO2− surface concentrations encountered during the cruise. Moreover, additional measurements of NO production rates implied that the occurrence of particles and a temperature increase can enhance NO production rates. With the method introduced here, we have a reliable and comparably easy to use method at hand to measure oceanic NO surface concentrations, which can be used to decipher both its temporal and spatial distributions as well as its biogeochemical pathways in the oceans.

Description: There is currently no consensus on how humans are affecting the marine nitrogen (N) cycle, which limits marine biological production and CO2 uptake. Anthropogenic changes in ocean warming, deoxygenation, and atmospheric N deposition can all individually affect the marine N cycle and the oceanic production of the greenhouse gas nitrous oxide (N2O). However, the combined effect of these perturbations on marine N cycling, ocean productivity, and marine N2O production is poorly understood. Here we use an Earth system model of intermediate complexity to investigate the combined effects of estimated 21st century CO2 atmospheric forcing and atmospheric N deposition. Our simulations suggest that anthropogenic perturbations cause only a small imbalance to the N cycle relative to preindustrial conditions (∼+5 Tg N y−1 in 2100). More N loss from water column denitrification in expanded oxygen minimum zones (OMZs) is counteracted by less benthic denitrification, due to the stratification-induced reduction in organic matter export. The larger atmospheric N load is offset by reduced N inputs by marine N2 fixation. Our model predicts a decline in oceanic N2O emissions by 2100. This is induced by the decrease in organic matter export and associated N2O production and by the anthropogenically driven changes in ocean circulation and atmospheric N2O concentrations. After comprehensively accounting for a series of complex physical-biogeochemical interactions, this study suggests that N flux imbalances are limited by biogeochemical feedbacks that help stabilize the marine N inventory against anthropogenic changes. These findings support the hypothesis that strong negative feedbacks regulate the marine N inventory on centennial time scales.

Description: When volcanic mountains slide into the sea, they trigger tsunamis. How big are these waves, and how far away can they do damage? Ritter Island provides some answers.

Description: Proxy data suggest the onset of Northern Hemisphere glaciation during the Plio-Pleistocene transition from 3.2 to 2.5 Ma resulted in enhanced climate variability at the obliquity (41 kyr) frequency. Here, we investigate the influence of the expanding Greenland ice sheet (GrIS) on the mean climate and obliquity-related variability in a series of climate model simulations. These suggest that an expanding GrIS weakens the Atlantic Meridional Overturning Circulation (AMOC) by ~1 Sv, mainly due to reduced heat loss in the Greenland-Iceland-Norwegian Sea. Moreover, the growing GrIS amplifies the Hadley circulation response to obliquity forcing driving variations in freshwater export from the tropical Atlantic and in turn variations of the AMOC. The stronger AMOC response to obliquity forcing, by about a factor of two, results in a stronger global-mean near-surface temperature response. We conclude that the AMOC response to obliquity forcing is important to understand the enhanced climate variability at the obliquity frequency during the Plio-Pleistocene transition.

Description: In this pilot study we link the yield of industrial fisheries to changes in the zooplankton mortality in an idealized way accounting for different target species (planktivorous fish—decreased zooplankton mortality; large predators—increased zooplankton mortality). This indirect approach is used in a global coupled biogeochemistry circulation model to estimate the range of the potential impact of industrial fisheries on marine biogeochemistry. The simulated globally integrated response on phytoplankton and primary production is in line with expectations—a high (low) zooplankton mortality results in a decrease (increase) of zooplankton and an increase (decrease) of phytoplankton. In contrast, the local response of zooplankton and phytoplankton depends on the region under consideration: In nutrient-limited regions, an increase (decrease) in zooplankton mortality leads to a decrease (increase) in both zooplankton and phytoplankton biomass. In contrast, in nutrient-replete regions, such as upwelling regions, we find an opposing response: an increase (decrease) of the zooplankton mortality leads to an increase (decrease) in both zooplankton and phytoplankton biomass. The results are further evaluated by relating the potential fisheries-induced changes in zooplankton mortality to those driven by CO2 emissions in a business-as-usual 21st century emission scenario. In our idealized case, the potential fisheries-induced impact can be of similar size as warming-induced changes in marine biogeochemistry.

Description: Back-arc spreading centers (BASCs) form a distinct class of ocean spreading ridges distinguished by steep along-axis gradients in spreading rate and by additional magma supplied through subduction. These characteristics can affect the population and distribution of hydrothermal activity on BASCs compared to mid-ocean ridges (MORs). To investigate this hypothesis, we comprehensively explored 600 km of the southern half of the Mariana BASC. We used water column mapping and seafloor imaging to identify 19 active vent sites, an increase of 13 over the current listing in the InterRidge Database (IRDB), on the bathymetric highs of 7 of the 11 segments. We identified both high and low (i.e., characterized by a weak or negligible particle plume) temperature discharge occurring on segment types spanning dominantly magmatic to dominantly tectonic. Active sites are concentrated on the two southernmost segments, where distance to the adjacent arc is shortest (〈40 km), spreading rate is highest (〉48 mm/yr), and tectonic extension is pervasive. Re-examination of hydrothermal data from other BASCs supports the generalization that hydrothermal site density increases on segments 〈90 km from an adjacent arc. Although exploration quality varies greatly among BASCs, present data suggest that, for a given spreading rate, the mean spatial density of hydrothermal activity varies little between MORs and BASCs. The present global database, however, may be misleading. On both BASCs and MORs, the spatial density of hydrothermal sites mapped by high-quality water-column surveys is 2–7 times greater than predicted by the existing IRDB trend of site density versus spreading rate.

Description: Provenance studies of widely distributed tephras, integrated within a well-defined temporal framework, are important to deduce systematic changes in the source, scale, distribution and changes in regional explosive volcanism. Here, we establish a robust tephro-chronostratigraphy for a total of 157 marine tephra layers collected during IODP Expedition 352. We infer at least three major phases of highly explosive volcanism during Oligocene to Pleistocene time. Provenance analysis based on glass composition assigns 56 of the tephras to a Japan source, including correlations with 12 major and widespread tephra layers resulting from individual eruptions in Kyushu, Central Japan and North Japan between 115 ka and 3.5 Ma. The remaining 101 tephras are assigned to four source regions along the Izu-Bonin arc. One, of exclusively Oligocene age, is proximal to the Bonin Ridge islands; two reflect eruptions within the volcanic front and back-arc of the central Izu-Bonin arc, and a fourth region corresponds to the Northern Izu-Bonin arc source. First-order volume estimates imply eruptive magnitudes ranging from 6.3 to 7.6 for Japan-related eruptions and between 5.5 and 6.5 for IBM eruptions. Our results suggest tephras between 30 and 22 Ma that show a subtly different Izu-Bonin chemical signature compared to the recent arc. After a ∼11 m.y. gap in eruption, tephra supply from the Izu-Bonin arc predominates from 15 to 5 Ma, and finally a subequal mixture of tephra sources from the (palaeo)Honshu and Izu-Bonin arcs occurs within the last ∼5 Ma.

Description: Igneous intrusions in sedimentary basins may have a profound effect on the thermal structure and physical properties of the hosting sedimentary rocks. These include mechanical effects such as deformation and uplift of sedimentary layers, generation of overpressure, mineral reactions and porosity evolution, and fracturing and vent formation following devolatilization reactions and the generation of CO2 and CH4. The gas generation and subsequent migration and venting may have contributed to several of the past climatic changes such as the end-Permian event and the Paleocene-Eocene Thermal Maximum. Additionally, the generation and expulsion of hydrocarbons and cracking of pre-existing oil reservoirs around a hot magmatic intrusion is of significant interest to the energy industry. In this paper, we present a user-friendly 1D FEM based tool, SILLi, which calculates the thermal effects of sill intrusions on the enclosing sedimentary stratigraphy. The model is accompanied by three case studies of sills emplaced in two different sedimentary basins, the Karoo Basin in South Africa and the Vøring Basin offshore Norway. Input data for the model is the present-day well log or sedimentary column with an Excel input file and includes rock parameters such as thermal conductivity, total organic carbon (TOC) content, porosity, and latent heats. The model accounts for sedimentation and burial based on a rate calculated by the sedimentary layer thickness and age. Erosion of the sedimentary column is also included to account for realistic basin evolution. Multiple sills can be emplaced within the system with varying ages. The emplacement of a sill occurs instantaneously. The model can be applied to volcanic sedimentary basins occurring globally. The model output includes the thermal evolution of the sedimentary column through time, and the changes that take place following sill emplacement such as TOC changes, thermal maturity, and the amount of organic and carbonate-derived CO2. The TOC and vitrinite results can be readily benchmarked within the tool to present-day values measured within the sedimentary column. This allows the user to determine the conditions required to obtain results that match observables and leads to a better understanding of metamorphic processes in sedimentary basins.

Description: Past warm periods provide an opportunity to evaluate climate models under extreme forcing scenarios, in particular high ( 〉  800 ppmv) atmospheric CO2 concentrations. Although a post hoc intercomparison of Eocene ( ∼  50  Ma) climate model simulations and geological data has been carried out previously, models of past high-CO2 periods have never been evaluated in a consistent framework. Here, we present an experimental design for climate model simulations of three warm periods within the early Eocene and the latest Paleocene (the EECO, PETM, and pre-PETM). Together with the CMIP6 pre-industrial control and abrupt 4 ×  CO2 simulations, and additional sensitivity studies, these form the first phase of DeepMIP – the Deep-time Model Intercomparison Project, itself a group within the wider Paleoclimate Modelling Intercomparison Project (PMIP). The experimental design specifies and provides guidance on boundary conditions associated with palaeogeography, greenhouse gases, astronomical configuration, solar constant, land surface processes, and aerosols. Initial conditions, simulation length, and output variables are also specified. Finally, we explain how the geological data sets, which will be used to evaluate the simulations, will be developed.

Description: The Stratospheric Sulfur and its Role in Climate (SSiRC) interactive stratospheric aerosol model intercomparison project (ISA-MIP) explores uncertainties in the processes that connect volcanic emission of sulphur gas species and the radiative forcing associated with the resulting enhancement of the stratospheric aerosol layer. The central aim of ISA-MIP is to constrain and improve interactive stratospheric aerosol models and reduce uncertainties in the stratospheric aerosol forcing by comparing results of standardized model experiments with a range of observations. In this paper we present 4 co-ordinated inter-model experiments designed to investigate key processes which influence the formation and temporal development of stratospheric aerosol in different time periods of the observational record. The "Background" (BG) experiment will focus on microphysics and transport processes under volcanically quiescent conditions, when the stratospheric aerosol is controlled by the transport of aerosols and their precursors from the troposphere to the stratosphere. The "Transient Aerosol Record" (TAR) experiment will explore the role of small- to moderate-magnitude volcanic eruptions, anthropogenic sulphur emissions and transport processes over the period 1998–2012 and their role in the warming hiatus. Two further experiments will investigate the stratospheric sulphate aerosol evolution after major volcanic eruptions. The "Historical Eruptions SO2 Emission Assessment" (HErSEA) experiment will focus on the uncertainty in the initial emission of recent large-magnitude volcanic eruptions, while the "Pinatubo Emulation in Multiple models" (PoEMS) experiment will provide a comprehensive uncertainty analysis of the radiative forcing from the 1991 Mt. Pinatubo eruption.

Description: Shorelines exhibit long-range dependence (LRD) and have been shown in some environments to be described in the wavenumber domain by a power law characteristic of scale-independence. Recent evidence suggests that the geomorphology of barrier islands can, however, exhibit scale dependence as a result of systematic variations of the underlying framework geology. The LRD of framework geology, which influences island geomorphology and its response to storms and sea level rise, has not been previously examined. Electromagnetic induction (EMI) surveys conducted along Padre Island National Seashore (PAIS), Texas, USA, reveal that the EMI apparent conductivity σa signal and, by inference, the framework geology exhibits LRD at scales up to 101 to 102 km. Our study demonstrates the utility of describing EMI σa and LiDAR spatial series by a fractional auto-regressive integrated moving average process that specifically models LRD. This method offers a robust and compact way for quantifying the geological variations along a barrier island shoreline using three parameters (p,d,q). We discuss how ARIMA (0,d,0) models that use a single parameter d provide a quantitative measure for determining free and forced barrier island evolutionary behavior across different scales. Statistical analyses at regional, intermediate, and local scales suggest that the geologic framework within an area of paleo-channels exhibits a first order control on dune height. The exchange of sediment amongst nearshore, beach and dune in areas outside this region are scale-independent, implying that barrier islands like PAIS exhibit a combination of free and forced behaviors that affect the response of the island to sea level rise.

Description: The influence of mesoscale eddies on the flow field and the water masses, especially the oxygen distribution of the eastern tropical South Pacific, is investigated from a mooring, float, and satellite data set. Two anticyclonic (ACE1/2), one mode-water (MWE), and one cyclonic eddy (CE) are identified and followed in detail with satellite data on their westward transition with velocities of 3.2 to 6.0cms−1 from their generation region, the shelf of the Peruvian and Chilean upwelling regime, across the Stratus Ocean Reference Station (ORS;  ∼ 20°S, 85°W) to their decaying region far west in the oligotrophic open ocean. The ORS is located in the transition zone between the oxygen minimum zone and the well oxygenated South Pacific subtropical gyre. Velocity, hydrographic, and oxygen measurements at the mooring show the impact of eddies on the weak flow region of the eastern tropical South Pacific. Strong anomalies are related to the passage of eddies and are not associated with a seasonal signal in the open ocean. The mass transport of the four observed eddies across 85°W is between 1.1 and 1.8Sv. The eddy type-dependent available heat, salt, and oxygen anomalies are 8.1×1018J (ACE2), 1.0×1018J (MWE), and −8.9×1018J (CE) for heat; 25.2×1010kg (ACE2), −3.1×1010kg (MWE), and −41.5×1010kg (CE) for salt; and −3.6×1016µmol (ACE2), −3.5×1016µmol (MWE), and −6.5×1016µmol (CE) for oxygen showing a strong imbalance between anticyclones and cyclones for salt transports probably due to seasonal variability in water mass properties in the formation region of the eddies. Heat, salt, and oxygen fluxes out of the coastal region across the ORS region in the oligotrophic open South Pacific are estimated based on these eddy anomalies and on eddy statistics (gained out of 23 years of satellite data). Furthermore, four profiling floats were trapped in the ACE2 during its westward propagation between the formation region and the open ocean, which allows for conclusions on lateral mixing of water mass properties with time between the core of the eddy and the surrounding water. The strongest lateral mixing was found between the seasonal thermocline and the eddy core during the first half of the eddy lifetime.

Description: Subsurface coherent eddies are well-known features of ocean circulation, but the sparsity of observations prevents an assessment of their importance for biogeochemistry. Here, we use a global eddying (0.1° ) ocean-biogeochemical model to carry out a census of subsurface coherent eddies originating from eastern boundary upwelling systems (EBUS), and quantify their biogeochemical effects as they propagate westward into the subtropical gyres. While most eddies exist for a few months, moving over distances of 100s of km, a small fraction (〈 5%) of long-lived eddies propagates over distances greater than 1000km, carrying the oxygen-poor and nutrient-rich signature of EBUS into the gyre interiors. In the Pacific, transport by subsurface coherent eddies accounts for roughly 10% of the offshore transport of oxygen and nutrients in pycnocline waters. This "leakage" of subsurface waters can be a significant fraction of the transport by nutrient-rich poleward undercurrents, and may contribute to the well-known reduction of productivity by eddies in EBUS. Furthermore, at the density layer of their cores, eddies decrease climatological oxygen locally by close to 10%, thereby expanding oxygen minimum zones. Finally, eddies represent low-oxygen extreme events in otherwise oxygenated waters, increasing the area of hypoxic waters by several percent and producing dramatic short-term changes that may play an important ecological role. Capturing these non-local effects in global climate models, which typically include non-eddying oceans, would require dedicated parameterizations.

Description: Climate models depict large diversity in the strength of the El Niño/Southern Oscillation (ENSO) (ENSO amplitude). Here we investigate ENSO-amplitude diversity in the Coupled Model Intercomparison Project Phase 5 (CMIP5) by means of the linear recharge oscillator model, which reduces ENSO dynamics to a two-dimensional problem in terms of eastern equatorial Pacific sea surface temperature anomalies (T) and equatorial Pacific upper ocean heat content anomalies (h). We find that a large contribution to ENSO-amplitude diversity originates from stochastic forcing. Further, significant interactions exist between the stochastic forcing and the growth rates of T and h with competing effects on ENSO amplitude. The joint consideration of stochastic forcing and growth rates explains more than 80% of the ENSO-amplitude variance within CMIP5. Our results can readily explain the lack of correlation between the Bjerknes Stability index, a measure of the growth rate of T, and ENSO amplitude in a multimodel ensemble.

Description: Monthly mean sea level anomalies in the tropical Pacific for the period 1961-2002 are reconstructed using a linear, multi-mode model driven by monthly mean wind stress anomalies from the NCEP/NCAR and ERA-40 reanalysis products. Overall, the sea level anomalies reconstructed by both wind stress products agree well with the available tide gauge data, although with poor performance at Kanton Island in the western-central equatorial Pacific and reduced amplitude at Christmas Island. The reduced performance is related to model error in locating the pivot point in sea level variability associated with the so-called “tilt” mode. We present evidence that the pivot point was further west during the period 1993-2014 than during the period 1961-2002 and attribute this to a persistent upward trend in the zonal wind stress variance along the equator west of 160° W throughout the period 1961-2014. Experiments driven by the zonal component of the wind stress alone reproduce much of the trend in sea level found in the experiments driven by both components of the wind stress. The experiments show an upward trend in sea level in the eastern tropical Pacific over the period 1961-2002, but with a much stronger upward trend when using the NCEP/NCAR product. We argue that the latter is related to an overly strong eastward trend in zonal wind stress in the eastern-central Pacific that is believed to be a spurious feature of the NCEP/NCAR product.

Description: The eruption of Mt. Tambora in 1815 was the largest volcanic eruption of the past 500 years. The eruption had significant climatic impacts, leading to the 1816 "year without a summer", and remains a valuable event from which to understand the climatic effects of large stratospheric volcanic sulfur dioxide injections. The eruption also resulted in one of the strongest and most easily identifiable volcanic sulfate signals in polar ice cores, which are widely used to reconstruct the timing and atmospheric sulfate loading of past eruptions. As part of the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP), five state-of-the-art global aerosol models simulated this eruption. We analyse both simulated background (no Tambora) and volcanic (with Tambora) sulfate deposition to polar regions and compare to ice core records. The models simulate overall similar patterns of background sulfate deposition, al-though there are differences in regional details and magnitude. However, the volcanic sulfate deposition varies considerably between the models with differences in timing, spatial pattern and magnitude. Mean simulated deposited sulfate on Antarctica ranges from 19 to 264 kgkm-2 and on Greenland from 31 to 194 kgkm-2, as compared to the mean ice-corederived estimates of roughly 50 kgkm-2 for both Greenland and Antarctica. The ratio of the hemispheric atmospheric sulfate aerosol burden after the eruption to the average ice sheet deposited sulfate varies between models by up to a factor of 15. Sources of this inter-model variability include differences in both the formation and the transport of sulfate aerosol. Our results suggest that deriving relationships between sulfate deposited on ice sheets and atmospheric sulfate burdens from model simulations may be associated with greater uncertainties than previously thought.

Description: Oceanic very short-lived substances (VSLSs), such as bromoform (CHBr3), contribute to stratospheric halogen loading and, thus, to ozone depletion. However, the amount, timing, and region of bromine delivery to the stratosphere through one of the main entrance gates, the Indian summer monsoon circulation, are still uncertain. In this study, we created two bromoform emission inventories with monthly resolution for the tropical Indian Ocean and west Pacific based on new in situ bromoform measurements and novel ocean biogeochemistry modeling. The mass transport and atmospheric mixing ratios of bromoform were modeled for the year 2014 with the particle dispersion model FLEXPART driven by ERA-Interim reanalysis. We compare results between two emission scenarios: (1) monthly averaged and (2) annually averaged emissions. Both simulations reproduce the atmospheric distribution of bromoform from ship- and aircraft-based observations in the boundary layer and upper troposphere above the Indian Ocean reasonably well. Using monthly resolved emissions, the main oceanic source regions for the stratosphere include the Arabian Sea and Bay of Bengal in boreal summer and the tropical west Pacific Ocean in boreal winter. The main stratospheric injection in boreal summer occurs over the southern tip of India associated with the high local oceanic sources and strong convection of the summer monsoon. In boreal winter more bromoform is entrained over the west Pacific than over the Indian Ocean. The annually averaged stratospheric injection of bromoform is in the same range whether using monthly averaged or annually averaged emissions in our Lagrangian calculations. However, monthly averaged emissions result in the highest mixing ratios within the Asian monsoon anticyclone in boreal summer and above the central Indian Ocean in boreal winter, while annually averaged emissions display a maximum above the west Indian Ocean in boreal spring. In the Asian summer monsoon anticyclone bromoform atmospheric mixing ratios vary by up to 50% between using monthly averaged and annually averaged oceanic emissions. Our results underline that the seasonal and regional stratospheric bromine injection from the tropical Indian Ocean and west Pacific critically depend on the seasonality and spatial distribution of the VSLS emissions.

Description: The calcareous tubeworm Spirorbis spirorbis is a widespread serpulid species in the Baltic Sea, where it commonly grows as an epibiont on brown macroalgae (genus Fucus). It lives within a Mg-calcite shell and could be affected by ocean acidification and temperature rise induced by the predicted future atmospheric CO2 increase. However, Spirorbis tubes grow in a chemically modified boundary layer around the algae, which may mitigate acidification. In order to investigate how increasing temperature and rising pCO2 may influence S. spirorbis shell growth we carried out four seasonal experiments in the Kiel Outdoor Benthocosms at elevated pCO2 and temperature conditions. Compared to laboratory batch culture experiments the benthocosm approach provides a better representation of natural conditions for physical and biological ecosystem parameters, including seasonal variations. We find that growth rates of S. spirorbis are significantly controlled by ontogenetic and seasonal effects. The length of the newly grown tube is inversely related to the initial diameter of the shell. Our study showed no significant difference of the growth rates between ambient atmospheric and elevated (1100 ppm) pCO2 conditions. No influence of daily average CaCO3 saturation state on the growth rates of S. spirorbis was observed. We found, however, net growth of the shells even in temporarily undersaturated bulk solutions, under conditions that concurrently favoured selective shell surface dissolution. The results suggest an overall resistance of S. spirorbis growth to acidification levels predicted for the year 2100 in the Baltic Sea. In contrast, S. spirorbis did not survive at mean seasonal temperatures exceeding 24 °C during the summer experiments. In the autumn experiments at ambient pCO2, the growth rates of juvenile S. spirorbis were higher under elevated temperature conditions. The results reveal that S. spirorbis may prefer moderately warmer conditions during their early life stages but will suffer from an excessive temperature increase and from increasing shell corrosion as a consequence of progressing ocean acidification.

Description: n the Paris Agreement in 2015 countries agreed on holding global mean surface air warming to "well below 2 degrees C above pre-industrial" levels, but the emission reduction pledges under that agreement are not ambitious enough to meet this target. Therefore, the question arises of whether restoring global warming to this target after exceeding it by artificially removing CO2 from the atmosphere is possible. One important aspect is the reversibility of ocean heat uptake and associated sea level rise, which have very long (centennial to millennial) response timescales. In this study the response of sea level rise due to thermal expansion to a 1% yearly increase of atmospheric CO2 up to a quadrupling of the pre-industrial concentration followed by a 1% yearly decline back to the pre-industrial CO2 concentration is examined using the University of Victoria Earth System Climate Model (UVic ESCM). We find that global mean thermosteric sea level (GMTSL) continues to rise for several decades after atmospheric CO2 starts to decline and does not return to pre-industrial levels for over 1000 years after atmospheric CO2 is restored to the pre-industrial concentration. This finding is independent of the strength of vertical sub-grid-scale ocean mixing implemented in the model. Furthermore, GMTSL rises faster than it declines in response to a symmetric rise and decline in atmospheric CO2 concentration partly because the deep ocean continues to warm for centuries after atmospheric CO2 returns to the pre-industrial concentration. Both GMTSL rise and decline rates increase with increasing vertical ocean mixing. Exceptions from this behaviour arise if the overturning circulations in the North Atlantic and Southern Ocean intensify beyond pre-industrial levels in model versions with lower vertical mixing, which leads to rapid cooling of the deep ocean.

Description: Continental margin ecosystems in the western North Pacific Ocean are subject to strong climate forcing and anthropogenic impacts. To evaluate mechanisms controlling phytoplankton biomass and community structure variations in marginal sea‐open ocean boundary regions, brassicasterol, dinosterol and C37 alkenones were measured in suspended particles in summer and autumn from 2012 to 2013 in the northeastern East China Sea and the western Tsushima Strait (NEECS‐WTS). In summer, the concentrations of brassicasterol (40 ‐ 1535 ng L‐1) and dinosterol (4.2 ‐ 94 ng L‐1) were higher in the southwest of Cheju Island, while C37 alkenones (0 ‐ 30 ng L‐1) were higher in the south of Cheju Island. In autumn, brassicasterol (12 ‐ 106 ng L‐1), dinosterol (2.4 ‐ 21 ng L‐1) and C37 alkenones (0.7 – 7.0 ng L‐1) were higher in the southwest of Cheju Island and the WTS, and higher C37 alkenones also occurred in the Okinawa Trough. Correlation analysis of biomarkers and environmental conditions (temperature, salinity and inorganic nutrient concentrations) clearly demonstrated that phytoplankton biomass and community structure variations can be well elucidated by water masses as indexed by temperature and salinity. High nutrients from the Changjiang River were the main cause of high biomass in summer, while nutrients from subsurface water were likely the key factor regulating phytoplankton biomass in open ocean water stations in autumn. This study indicates that mechanisms controlling phytoplankton biomass in marginal sea‐open ocean boundary regions should be classified by various water masses with different nutrient concentrations, instead of by geography.

Description: A submesoscale coherent vortex (SCV) with a low oxygen core is characterized from underwater glider and mooring observations from the eastern tropical North Atlantic, north of the Cape Verde Islands. The eddy crossed the mooring with its center and a 1 month time series of the SCV's hydrographic and upper 100 m currents structure was obtained. About 45 days after, and ∼100 km west, the SCV frontal zone was surveyed in high temporal and spatial resolution using an underwater glider. Satellite altimetry showed the SCV was formed about 7 months before at the Mauritanian coast. The SCV was located at 80-100 m depth, its diameter was ∼100 km and its maximum swirl velocity ∼0.4 m s-1. A Burger number of 0.2 and a vortex Rossby number 0.15 indicate a flat lens in geostrophic balance. Mooring and glider data show in general comparable dynamical and thermohaline structures, the glider in high spatial resolution, the mooring in high temporal resolution. Surface maps of chlorophyll concentration suggest high productivity inside and around the SCV. The low potential vorticity (PV) core of the SCV is surrounded by filamentary structures, sloping down at different angles from the mixed layer base and with typical width of 10-20 km and a vertical extent of 50-100 m.

Description: We present an improved neotectonic numerical model of the complex NW Africa-SW Eurasia plate boundary segment that runs from west to east along the Gloria Fault up to the northern Algerian margin. We model the surface velocity field and the ongoing lithospheric deformation using the most recent version of the thin-shell code SHELLS and updated lithospheric model and fault map of the region. To check the presence versus the absence of an independently driven Alboran domain, we develop two alternative plate models: one does not include an Alboran plate; another includes it and determines the basal shear tractions necessary to drive it with known velocities. We also compare two alternative sets of Africa-Eurasia velocity boundary conditions, corresponding to geodetic and geological-scale averages of plate motion. Finally, we perform an extensive parametric study of fault friction coefficient, trench resistance, and velocities imposed in Alboran nodes. The final run comprises 5240 experiments, each scored to geodetic velocities (estimated for 250 stations and here provided), stress direction data, and seismic strain rates. The model with the least discrepancy to the data includes the Alboran plate driven by a basal WSW directed shear traction, slightly oblique to the westward direction of Alboran motion. We provide estimates of long-term strain rates and slip rates for the modeled faults, which can be useful for further hazard studies. Our results support that a mechanism additional to the Africa-Eurasia convergence is required to drive the Alboran domain, which can be related to subduction processes occurring within the mantle.

Description: Overriding plate topography provides constraints on subduction zone geodynamics. We investigate its evolution using fully dynamic laboratory models of subduction with techniques of stereoscopic photogrammetry and particle image velocimetry. Model results show that the topography is characterized by an area of forearc dynamic subsidence, with a magnitude scaling to 1.44–3.97 km in nature, and a local topographic high between the forearc subsided region and the trench. These topographic features rapidly develop during the slab free‐sinking phase and gradually decrease during the steady state slab rollback phase. We propose that they result from the variation of the vertical component of the trench suction force along the subduction zone interface, which gradually increases with depth and results from the gradual slab steepening during the initial transient slab sinking phase. The downward mantle flow in the nose of the mantle wedge plays a minor role in driving forearc subsidence.

Description: Earth is 180 Myr into the current supercontinent cycle, and the next supercontinent is predicted to form in 250 Myr. The continuous changes in continental configuration can move the ocean between resonant states, and the semidiurnal tides are currently large compared to the past 252 Myr due to tidal resonance in the Atlantic. This leads to the hypothesis that there is a “supertidal” cycle linked to the supercontinent cycle. Here this is tested using new tectonic predictions for the next 250 Myr as bathymetry in a numerical tidal model. The simulations support the following hypothesis: a new tidal resonance will appear 150 Myr from now, followed by a decreasing tide as the supercontinent forms 100 Myr later. This affects the dissipation of tidal energy in the oceans, with consequences for the evolution of the Earth‐Moon system, ocean circulation and climate, and implications for the ocean's capacity of hosting and evolving life.

Description: The Antarctic Circumpolar Current (ACC) has highly energetic mesoscale phenomena, but their impacts on phytoplankton biomass, productivity, and biogeochemical cycling are not understood well. We analyze satellite observations and an eddy‐rich ocean model to show that they drive chlorophyll anomalies of opposite sign in winter versus summer. In winter, deeper mixed layers in positive sea surface height (SSH) anomalies reduce light availability, leading to anomalously low chlorophyll concentrations. In summer with abundant light, however, positive SSH anomalies show elevated chlorophyll concentration due to higher iron level, and an iron budget analysis reveals that anomalously strong vertical mixing enhances iron supply to the mixed layer. Features with negative SSH anomalies exhibit the opposite tendencies: higher chlorophyll concentration in winter and lower in summer. Our results suggest that mesoscale modulation of iron supply, light availability and vertical mixing plays an important role in causing systematic variations in primary productivity over the seasonal cycle.

Description: Halogen- and sulfur-containing compounds are supersaturated in the surface ocean, which results in their emission to the atmosphere. These compounds can be transported to the stratosphere, where they impact ozone, the background aerosol layer, and climate. In this study we calculate the seasonal and interannual variability of transport from the West Indian Ocean (WIO) surface to the stratosphere for 2000-2016 with the Lagrangian transport model FLEXPART using ERA-Interim meteorological fields. We investigate the transport relevant for very short lived substances (VSLS) with tropospheric lifetimes corresponding to dimethylsulfide (1 day), methyl iodide (CH3I, 3.5 days), bromoform (CHBr3, 17 days), and dibromomethane (CH2Br2, 150 days). The stratospheric source gas injection of VSLS tracers from the WIO shows a distinct annual cycle associated with the Asian monsoon. Over the 16-year time series, a slight increase in source gas injection from the WIO to the stratosphere is found for all VSLS tracers and during all seasons. The interannual variability shows a relationship with sea surface temperatures in the WIO as well as the El Niño-Southern Oscillation. During boreal spring of El Niño, enhanced stratospheric injection of VSLS from the tropical WIO is caused by positive sea surface temperature anomalies and enhanced vertical uplift above the WIO. During boreal fall of La Niña, strong injection is related to enhanced atmospheric upward motion over the East Indian Ocean and a prolonged Indian summer monsoon season. Related physical mechanisms and uncertainties are discussed in this study

Description: The notion of a shallow northern sourced intermediate water mass is a well evidenced feature of the Atlantic circulation scheme of the Last Glacial Maximum (LGM). However, recent observations from stable carbon isotopes (δ13C) at the Corner Rise in the deep northwest Atlantic suggested a significant contribution of a Northern Component Water mass to the abyssal northwest Atlantic basin that has not been described before. Here we test the hypothesis of this northern sourced water mass underlying the southern sourced glacial Antarctic Bottom Water by measuring the authigenic neodymium (Nd) isotopic composition from the same sediments from 5,010-m water depth. Neodymium isotopes act as a semiconservative water mass tracer capable of distinguishing between Northern and Southern Component Waters at the northwest Atlantic. Our new Nd isotopic record resolves various water mass changes from the LGM to the early Holocene in agreement with existing Nd-based reconstructions from across the west Atlantic Ocean. Especially pronounced are the Younger Dryas and Bølling-Allerød with unprecedented changes in the Nd isotopic composition. For the LGM we found Nd isotopic evidence for a northern sourced water mass contributing to abyssal depths, thus being in agreement with previous δ13C data from Corner Rise. Overall, however, the deep northwest Atlantic was still dominated by southern sourced water, since we found signatures that are intermediate between northern and southern end member compositions. Furthermore, this new record indicates that C and Nd isotopes were partly decoupled, pointing to nonconservative behavior of one or more likely of both water mass proxies during the LGM.

Description: Pre-stack depth migration data across the Hikurangi margin, East Coast of the North Island, New Zealand, are used to derive subducting slab geometry, upper crustal structure and seismic velocities resolved to ∼14 km depth. We investigate the potential relationship between the crustal architecture, fluid migration and short-term geodetically determined slow-slip events. The subduction interface is a shallow dipping thrust at 〈 7 km depth near the trench and steps down to 14 km depth along an ∼18 km long ramp, beneath Porangahau Ridge. This apparent bend in the décollement is associated with splay fault branching and coincides with a zone of maximum slip (90 mm) inferred on the subduction interface during slow slip events in June and July 2011. A low-velocity zone beneath the plate interface, up-dip of the plate interface ramp, is interpreted as fluid-rich overpressured sediments capped with a low permeability condensed layer of chalk and interbedded mudstones. Fluid rich sediments have been imbricated by splay faults in a region that coincides with the step down in the décollement from the top of subducting sediments to the oceanic crust and contribute to spatial variation in frictional properties of the plate interface that may promote slow slip behavior in the region. Further, transient fluid migration along splay faults at Porangahau Ridge may signify stress changes during slow slip.

Description: We provide high-resolution foraminiferal stable carbon isotope (δ13C) records from the subarctic Pacific and Eastern Equatorial Pacific (EEP) to investigate circulation dynamics between the extra-tropical and tropical North Pacific during the past 60 kyr. We measured the δ13C composition of the epibenthic foraminiferal species Cibicides lobatulus from a shallow sediment core recovered from the western Bering Sea (SO201-2-101KL; 58°52.52’N, 170°41.45’E; 630 m water depth) to reconstruct past ventilation changes close to the source region of Glacial North Pacific Intermediate Water (GNPIW). Information regarding glacial changes in the δ13C of sub-thermocline water masses in the EEP is derived from the deep-dwelling planktonic foraminifera Globorotaloides hexagonus at ODP Site 1240 (00°01.31’N, 82°27.76’W; 2921 m water depth). Apparent similarities in the long-term evolution of δ13C between GNPIW, intermediate waters in the eastern tropical North Pacific and sub-thermocline water masses in the EEP suggest the expansion of relatively 13C-depleted, nutrient-enriched, and northern-sourced intermediate waters to the equatorial Pacific under glacial conditions. Further, it appears that additional influence of GNPIW to the tropical Pacific is consistent with changes in nutrient distribution and biological productivity in surface-waters of the glacial EEP. Our findings highlight potential links between North Pacific mid-depth circulation changes, nutrient cycling, and biological productivity in the equatorial Pacific under glacial boundary conditions.

Description: A multi-mode, linear reduced-gravity model, driven by ERA-Interim monthly mean wind stress anomalies, is used to investigate interannual variability in tropical Pacific sea level as seen in satellite altimeter data. The model output is fitted to the altimeter data along the equator, in order to derive the vertical profile for the model forcing, showing that a signature from modes higher than mode six cannot be extracted from the altimeter data. It is shown that the model has considerable skill at capturing interannual sea level variability both on and off the equator. The correlation between modelled and satellite-derived sea level data exceeds 0.8 over a wide range of longitudes along the equator and readily captures the observed ENSO events. Overall, the combination of the first, second, third and fifth modes can provide a robust estimate of the interannual sea level variability, the second mode being dominant. A remarkable feature of both the model and the altimeter data is the presence of a pivot point in the western Pacific on the equator. We show that the westward displacement of the pivot point from the centre of the basin is strongly influenced by the fact that most of the wind stress variance is found in the western part of the basin. We also show that the Sverdrup transport is not fundamental to the dynamics of the recharge/discharge mechanism in our model, although the spatial structure of the wind forcing does play a role in setting the amplitude of the “warm water volume”.

Description: El Niño-Southern Oscillation (ENSO) in the Pacific is asymmetric for warm and cold events with respect to amplitude, spatial patterns and temporal evolution. Here the symmetry of the Atlantic Niño mode, which many previous studies have argued is governed by atmosphere–ocean dynamics similar to those of ENSO, is investigated using two different ocean reanalysis products. Calculation of Bjerknes feedback terms for the Pacific reveals a pronounced asymmetry between warm and cold events, though unlike most previous studies, the largest asymmetry is found in the relationship between eastern Pacific thermocline depth and SST anomalies. For the Atlantic, cold events are effectively mirror images of warm events with Bjerknes feedbacks of similar strength. The analysis supports not only the conclusion that Atlantic Niños are more symmetric than ENSO, but the hypothesis itself that the Bjerknes feedback is operative in the Atlantic given the strength of the relationship between the key variables involved.

Description: At the eastern end of the Azores-Gloria transform fault system to the southwest of Portugal, the plate boundary between Africa and Iberia is a region where deformation is accommodated over a wide tectonically-active area. The region has unleashed large earthquakes and tsunamis, including the Mw ~ 8.5 Great Lisbon earthquake of 1755. Although the source region of the 1755 earthquake is still disputed, most proposals include a source location in the vicinity of the Horseshoe Abyssal Plain (HAP), which is bounded by the 5000 m high Gorringe Bank (GB). In this study we characterise seismic activity in the region using data recorded by two local networks of ocean-bottom seismometers (OBS). The networks were deployed in the eastern HAP and at the GB. The dataset allowed the detection of 160 local earthquakes. These earthquakes cluster around the GB, to the SW of Cabo Sao Vicente, and in the HAP. Focal depths indicate deep-seated earthquakes, with depths increasing from 20-35 km (mean of 26.1 ± 7.2 km) at the GB to 15-45 km (mean 31.5 km ± 10.5 km) under the HAP. Seismic activity thus extends down to levels that are deeper than those mapped by active seismic profiling, with the majority of events occurring within the mantle. Thermal modelling suggests that temperatures of approximately 600 °C characterise the base of the seismogenic brittle lithosphere at ~45 km depth. The large source depth and thermal structure supports previous suggestions that catastrophic seismic rupture through the lithospheric mantle may indeed occur in the area.

Description: Despite a growing literature on the climate response to solar geoengineering – proposals to cool the planet by increasing the planetary albedo – there has been little published on the impacts of solar geoengineering on natural and human systems such as agriculture, health, water resources, and ecosystems. An understanding of the impacts of different scenarios of solar geoengineering deployment will be crucial for informing decisions on whether and how to deploy it. Here we review the current state of knowledge about impacts of a solar geoengineered climate and identify major research gaps. We suggest that a thorough assessment of the climate impacts of a range of scenarios of solar geoengineering deployment is needed and can build upon existing frameworks. However, solar geoengineering poses a novel challenge for climate impacts research as the manner of deployment could be tailored to pursue different objectives making possible a wide range of climate outcomes. We present a number of ideas for approaches to extend the survey of climate impacts beyond standard scenarios of solar geoengineering deployment to address this challenge. Reducing the impacts of climate change is the fundamental motivator for emissions reductions and for considering whether and how to deploy solar geoengineering. This means that the active engagement of the climate impacts research community will be important for improving the overall understanding of the opportunities, challenges and risks presented by solar geoengineering.

Description: The supply and bioavailability of iron (Fe) controls primary productivity and N2-fixation in large parts of the global ocean. An important, yet poorly quantified, source to the ocean is particulate Fe (pFe). Here we present the first combined dataset of particulate, labile-particulate (L-pFe) and dissolved Fe (dFe) from the (sub)-tropical North Atlantic. We show a strong relationship between L-pFe and dFe, indicating a dynamic equilibrium between these two phases whereby particles ‘buffer’ dFe and maintain the elevated concentrations observed. Moreover, L-pFe can increase the overall ‘available’ (L-pFe + dFe) Fe pool by up to 55%. The lateral shelf flux of this available Fe was similar in magnitude to observed soluble aerosol-Fe deposition, a comparison that has not been previously considered. These findings demonstrate that L-pFe is integral to Fe cycling and hence plays a role in regulating carbon cycling, warranting its’ inclusion in Fe budgets and biogeochemical models.

Description: Combined seawater radiogenic hafnium (Hf) and neodymium (Nd) isotope compositions were extracted from bulk sediment leachates and foraminifera of Site 1088, ODP Leg 177, 2082 m water depth on the Agulhas Ridge. The new data provide a continuous reconstruction of long and short-term changes in ocean circulation and continental weathering inputs since the Mid-Miocene. Due to its intermediate water depth the sediments of this core sensitively recorded changes in admixture of North Atlantic Deep Water (NADW) to the Antarctic Circumpolar Current (ACC) as a function of the strength of the Atlantic Meridional Overturning Circulation (AMOC). Nd isotope compositions (εNd) range from -7 to -11 with glacial values generally 1 to 3 units more radiogenic than during the interglacials of the Quaternary. The data reveal episodes of significantly increased AMOC strength during late Miocene and Pliocene warm periods whereas peak radiogenic εNd values mark a strongly diminished AMOC during the major intensification of Northern Hemisphere Glaciation near 2.8 Ma and in the Pleistocene after 1.5 Ma. In contrast, the Hf isotope compositions (εHf) show an essentially continuous evolution from highly radiogenic values of up to +11 during the Miocene to less radiogenic present day values (+2 to +4) during the late Quaternary. The data document a long-term transition in dominant weathering inputs, where inputs from the South America are replaced by those from Southern Africa. Moreover, radiogenic peaks provide evidence for the supply of radiogenic Hf originating from Patagonian rocks to the Atlantic sector of the Southern Ocean via dust inputs.

Description: The climate research community uses atmospheric reanalysis data sets to understand a wide range of processes and variability in the atmosphere, yet different reanalyses may give very different results for the same diagnostics. The Stratosphere–troposphere Processes And their Role in Climate (SPARC) Reanalysis Intercomparison Project (S-RIP) is a coordinated activity to compare reanalysis data sets using a variety of key diagnostics. The objectives of this project are to identify differences among reanalyses and understand their underlying causes, to provide guidance on appropriate usage of various reanalysis products in scientific studies, particularly those of relevance to SPARC, and to contribute to future improvements in the reanalysis products by establishing collaborative links between reanalysis centres and data users. The project focuses predominantly on differences among reanalyses, although studies that include operational analyses and studies comparing reanalyses with observations are also included when appropriate. The emphasis is on diagnostics of the upper troposphere, stratosphere, and lower mesosphere. This paper summarizes the motivation and goals of the S-RIP activity and extensively reviews key technical aspects of the reanalysis data sets that are the focus of this activity. The special issue "The SPARC Reanalysis Intercomparison Project (S-RIP)" in this journal serves to collect research with relevance to the S-RIP in preparation for the publication of the planned two (interim and full) S-RIP reports.

Description: A high-resolution dynamical model (Nucleus for European Modelling of the Ocean, Mediterranean configuration – NEMO-MED12) was used to give the first simulation of the distribution of radiocarbon (14C) across the whole Mediterranean Sea. The simulation provides a descriptive overview of both the natural pre-bomb 14C and the entire anthropogenic radiocarbon transient generated by the atmospheric bomb tests performed in the 1950s and early 1960s. The simulation was run until 2011 to give the post-bomb distribution. The results are compared to available in situ measurements and proxy-based reconstructions. The radiocarbon simulation allows an additional and independent test of the dynamical model, NEMO-MED12, and its performance to produce the thermohaline circulation and deep-water ventilation. The model produces a generally realistic distribution of radiocarbon when compared with available in situ data. The results demonstrate the major influence of the flux of Atlantic water through the Strait of Gibraltar on the inter-basin natural radiocarbon distribution and characterize the ventilation of intermediate and deep water especially through the propagation of the anthropogenic radiocarbon signal. We explored the impact of the interannual variability on the radiocarbon distribution during the Eastern Mediterranean Transient (EMT) event. It reveals a significant increase in 14C concentration (by more than 60 ‰) in the Aegean deep water and at an intermediate level (value up to 10 ‰) in the western basin. The model shows that the EMT makes a major contribution to the accumulation of radiocarbon in the eastern Mediterranean deep waters.

Description: Reconstructed sea surface temperatures (SSTs) derived from Mg/Ca measurements in nine encrusting coralline algal skeletons from the Aleutian archipelago in the northernmost Pacific Ocean reveal an overall increase in SST from 1665 to 2007. In the Aleutian SST reconstruction, decadal-scale variability is a transient feature present during the 1700s and early 1800s and then fully emerging post-1950. SSTs vary coherently with available instrument records of cyclone variance and vacillate in and out of coherence with multicentennial Pacific Northwest drought reconstructions as a response to SST-driven alterations of storm tracks reaching North America. These results indicate that an influence of decadal-scale variability on the North Pacific storm tracks only became apparent during the midtwentieth century. Furthermore, what has been assumed as natural variability in the North Pacific, based on twentieth century instrumental data, is not consistent with the long-term natural variability evident in reconstructed SSTs predating the anthropogenic influence.

Description: Extensive biogeochemical transformation of organic matter takes place in the shallow continental shelf seas of Siberia. This, in combination with brine production from sea-ice formation, results in cold bottom waters with relatively high salinity and nutrient concentrations, as well as low oxygen and pH levels. Data from the SWERUS-C3 expedition with icebreaker Oden, from July to September 2014, show the distribution of such nutrient-rich, cold bottom waters along the continental margin from about 140 to 180° E. The water with maximum nutrient concentration, classically named the upper halocline, is absent over the Lomonosov Ridge at 140° E, while it appears in the Makarov Basin at 150° E and intensifies further eastwards. At the intercept between the Mendeleev Ridge and the East Siberian continental shelf slope, the nutrient maximum is still intense, but distributed across a larger depth interval. The nutrient-rich water is found here at salinities of up to ∼ 34.5, i.e. in the water classically named lower halocline. East of 170° E transient tracers show significantly less ventilated waters below about 150 m water depth. This likely results from a local isolation of waters over the Chukchi Abyssal Plain as the boundary current from the west is steered away from this area by the bathymetry of the Mendeleev Ridge. The water with salinities of ∼ 34.5 has high nutrients and low oxygen concentrations as well as low pH, typically indicating decay of organic matter. A deficit in nitrate relative to phosphate suggests that this process partly occurs under hypoxia. We conclude that the high nutrient water with salinity ∼ 34.5 are formed on the shelf slope in the Mendeleev Ridge region from interior basin water that is trapped for enough time to attain its signature through interaction with the sediment.

Description: Microbathymetry data, in situ observations, and sampling along the 138200N and 138200N oceanic core complexes (OCCs) reveal mechanisms of detachment fault denudation at the seafloor, links between tectonic extension and mass wasting, and expose the nature of corrugations, ubiquitous at OCCs. In the initial stages of detachment faulting and high-angle fault, scarps show extensive mass wasting that reduces their slope. Flexural rotation further lowers scarp slope, hinders mass wasting, resulting in morphologically complex chaotic terrain between the breakaway and the denuded corrugated surface. Extension and drag along the fault plane uplifts a wedge of hangingwall material (apron). The detachment surface emerges along a continuous moat that sheds rocks and covers it with unconsolidated rubble, while local slumping emplaces rubble ridges overlying corrugations. The detachment fault zone is a set of anostomosed slip planes, elongated in the alongextension direction. Slip planes bind fault rock bodies defining the corrugations observed in microbathymetry and sonar. Fault planes with extension-parallel stria are exposed along corrugation flanks, where the rubble cover is shed. Detachment fault rocks are primarily basalt fault breccia at 138200N OCC, and gabbro and peridotite at 138300N, demonstrating that brittle strain localization in shallow lithosphere form corrugations, regardless of lithologies in the detachment zone. Finally, faulting and volcanism dismember the 138300N OCC, with widespread present and past hydrothermal activity (Semenov fields), while the Irinovskoe hydrothermal field at the 138200N core complex suggests a magmatic source within the footwall. These results confirm the ubiquitous relationship between hydrothermal activity and oceanic detachment formation and evolution.