ALBERT

Description: Cold-water coral ecosystems are considered hot-spots of biodiversity and biomass production and may be a regionally important contributor to carbonate production. The impact of these ecosystems on biogeochemical processes and carbonate preservation in associated sediments were studied at Røst Reef and Traenadjupet Reef, two modern (post-glacial) cold-water coral reefs on the Mid-Norwegian shelf. Sulfate and iron reduction as well as carbonate dissolution and precipitation were investigated by combining pore-water geochemical profiles, steady state modeling, as well as solid phase analyses and sulfate reduction rate measurements on gravity cores of up to 3.25 m length. Low extents of sulfate depletion and dissolved inorganic carbon (DIC) production, combined with sulfate reduction rates not exceeding 3 nmol S cm−3 d−1, suggested that overall anaerobic carbon mineralization in the sediments was low. These data showed that the coral fragment-bearing siliciclastic sediments were effectively decoupled from the productive pelagic ecosystem by the complex reef surface framework. Organic matter being mineralized by sulfate reduction was calculated to consist of 57% carbon bound in CH2O groups and 43% carbon in -CH2- groups. Methane concentrations were below 1 μM, and failed to support the hypothesis of a linkage between the distribution of cold-water coral reefs and the presence of hydrocarbon seepage. Reductive iron oxide dissolution linked to microbial sulfate reduction buffered the pore-water carbonate system and inhibited acid-driven coral skeleton dissolution. A large pool of reactive iron was available leading to the formation of iron sulfide minerals. Constant pore-water Ca2+, Mg2+ and Sr2+ concentrations in most cores and decreasing Ca2+ and Sr2+ concentrations with depth in core 23–18 GC indicated diagenetic carbonate precipitation. This was consistent with the excellent preservation of buried coral fragments.

Description: In order to test the influences of ocean acidification on the ocean pelagic ecosystem, so far the largest CO2 manipulation mesocosm study (European Project on Ocean Acidification, EPOCA) was performed in Kings Bay (Kongsfjorden), Spitsbergen. During a 30 day incubation, bacterial diversity was investigated using DNA fingerprinting and clone library analysis of bacterioplankton samples. Terminal restriction fragment length polymorphism (T-RFLP) analysis of the PCR amplicons of the 16S rRNA genes revealed that general bacterial diversity, taxonomic richness and community structure were influenced by the variation of productivity during the time of incubation, but not the degree of ocean acidification. A BIOENV analysis suggested a complex control of bacterial community structure by various biological and chemical environmental parameters. The maximum apparent diversity of bacterioplankton (i.e., the number of T-RFs) in high and low pCO2 treatments differed significantly. A negative relationship between the relative abundance of Bacteroidetes and pCO2 levels was observed for samples at the end of the experiment by the combination of T-RFLP and clone library analysis. Our study suggests that ocean acidification affects the development of bacterial assemblages and potentially impacts the ecological function of the bacterioplankton in the marine ecosystem.

Description: The Atlantic meridional overturning circulation (AMOC) carries large amounts of heat into the North Atlantic influencing climate regionally as well as globally. Palaeo-records and simulations with comprehensive climate models suggest that the positive salt-advection feedback may yield a threshold behaviour of the system. That is to say that beyond a certain amount of freshwater flux into the North Atlantic, no meridional overturning circulation can be sustained. Concepts of monitoring the AMOC and identifying its vicinity to the threshold rely on the fact that the volume flux defining the AMOC will be reduced when approaching the threshold. Here we advance conceptual models that have been used in a paradigmatic way to understand the AMOC, by introducing a density-dependent parameterization for the Southern Ocean eddies. This additional degree of freedom uncovers a mechanism by which the AMOC can increase with additional freshwater flux into the North Atlantic, before it reaches the threshold and collapses: an AMOC that is mainly wind-driven will have a constant upwelling as long as the Southern Ocean winds do not change significantly. The downward transport of tracers occurs either in the northern sinking regions or through Southern Ocean eddies. If freshwater is transported, either atmospherically or via horizontal gyres, from the low to high latitudes, this would reduce the eddy transport and by continuity increase the northern sinking which defines the AMOC until a threshold is reached at which the AMOC cannot be sustained. If dominant in the real ocean this mechanism would have significant consequences for monitoring the AMOC.

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: 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: Pliensbachian and Toarcian Ostracoda first described by Exton (Geological Paper, Carleton University, Ottawa, 79: 1–104 1979) from the Lusitanian Basin, west-central Portugal have been re-examined. As a result, a greater diversity in the Ostracoda (80 species) is now recognized. Two species are newly described (Eucytherura zambujalensis sp. nov., Ektyphocythere mediodepressa sp. nov.) from the marls and calcareous shales of the Maria Pares Hill section near the village of Zambujal. Poor preservation precludes a complete taxonomic review of the present material. Five ostracod zones are proposed; Gammacythere ubiquita–Ogmoconchella gruendeli Zone, Poly cope cerasia–Polycope cincinnata Zone, Liasina lanceolata–Ogmoconcha convexa Zone, Bairdiacypris rectangularis–Kinkelinella sermoisensis Zone, and Cytherella toarcensis-Kinkelinella costata Zone. Although the ostracod assemblages possess strong similarities to those described from Northwest Europe, some of the Zambujal assemblages are dominated by the genus Polycope. A marked faunal turnover, in association with the extinction of the Metacopina occurs in the lower Subzone of the tenuicostatum Zone of Lower Toarcian age. These faunal events are discussed in relation to changing environmental conditions.

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: The late stage of the North East Atlantic (NEA) spring bloom was investigated during June 2005 along a transect section from 45 to 66° N between 15 and 20° W in order to characterize the contribution of siliceous and calcareous phytoplankton groups and describe their distribution in relation to environmental factors. We measured several biogeochemical parameters such as nutrients, surface trace metals, algal pigments, biogenic silica (BSi), particulate inorganic carbon (PIC) or calcium carbonate, particulate organic carbon, nitrogen and phosphorus (POC, PON and POP, respectively), as well as transparent exopolymer particles (TEP). Results were compared with other studies undertaken in this area since the JGOFS NABE program. Characteristics of the spring bloom generally agreed well with the accepted scenario for the development of the autotrophic community. The NEA seasonal diatom bloom was in the late stages when we sampled the area and diatoms were constrained to the northern part of our transect, over the Icelandic Basin (IB) and Icelandic Shelf (IS). Coccolithophores dominated the phytoplankton community, with a large distribution over the Rockall-Hatton Plateau (RHP) and IB. The Porcupine Abyssal Plain (PAP) region at the southern end of our transect was the region with the lowest biomass, as demonstrated by very low Chla concentrations and a community dominated by picophytoplankton. Early depletion of dissolved silicic acid (DSi) and increased stratification of the surface layer most likely triggered the end of the diatom bloom, leading to coccolithophore dominance. The chronic Si deficiency observed in the NEA could be linked to moderate Fe limitation, which increases the efficiency of the Si pump. TEP closely mirrored the distribution of both biogenic silica at depth and prymnesiophytes in the surface layer suggesting the sedimentation of the diatom bloom in the form of aggregates, but the relative contribution of diatoms and coccolithophores to carbon export in this area still needs to be resolved.

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: 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: 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: 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: 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: 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: 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: 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: 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: Here we show the use of the 210Pb-226Ra excess method to determine the growth rate of two corals from the world's largest known cold-water coral reef, Røst Reef, north of the Arctic circle off Norway. Colonies of each of the two species that build the reef, Lophelia pertusa and Madrepora oculata, were collected alive at 350 m depth using a submersible. Pb and Ra isotopes were measured along the major growth axis of both specimens using low level alpha and gamma spectrometry and trace element compositions were studied. 210Pb and 226Ra differ in the way they are incorporated into coral skeletons. Hence, to assess growth rates, we considered the exponential decrease of initially incorporated 210Pb, as well as the increase in 210Pb from the decay of 226Ra and contamination with 210Pb associated with Mn-Fe coatings that we were unable to remove completely from the oldest parts of the skeletons. 226Ra activity was similar in both coral species, so, assuming constant uptake of 210Pb through time, we used the 210Pb-226Ra chronology to calculate growth rates. The 45.5 cm long branch of M. oculata was 31 yr with an average linear growth rate of 14.4 ± 1.1 mm yr−1 (2.6 polyps per year). Despite cleaning, a correction for Mn-Fe oxide contamination was required for the oldest part of the colony; this correction corroborated our radiocarbon date of 40 yr and a mean growth rate of 2 polyps yr−1. This rate is similar to the one obtained in aquarium experiments under optimal growth conditions.

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: 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: 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: 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: 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: Cold-water coral (CWC) reefs are heterogeneous ecosystems comprising numerous microhabitats. A typical European CWC reef provides various biogenic microhabitats (within, on and surrounding colonies of coral species such as Lophelia pertusa, Paragorgia arborea and Primnoa resedaeformis, or formed by their remains after death). These microhabitats may be surrounded and intermixed with non-biogenic microhabitats (soft sediment, hard ground, gravel/pebbles, steep walls). To date, studies of distribution of sessile fauna across CWC reefs have been more numerous than those investigating mobile fauna distribution. In this study we quantified shrimp densities associated with key CWC microhabitat categories at the Røst Reef, Norway, by analysing image data collected by towed video sled in June 2007. We also investigated shrimp distribution patterns on the local scale (〈40 cm) and how these may vary with microhabitat. Shrimp abundances at the Røst Reef were on average an order of magnitude greater in biogenic reef microhabitats than in non-biogenic microhabitats. Greatest shrimp densities were observed in association with live Paragorgia arborea microhabitat (43 shrimp m−2, SD = 35.5), live Primnoa resedaeformis microhabitat (41.6 shrimp m−2, SD = 26.1) and live Lophelia pertusa microhabitat (24.4 shrimp m−2, SD = 18.6). In non-biogenic microhabitat, shrimp densities were 〈2 shrimp m−2. CWC reef microhabitats appear to support greater shrimp densities than the surrounding non-biogenic microhabitats at the Røst Reef, at least at the time of survey.

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: Submarine canyons are known as one of the seafloor morphological features where living cold-water coral (CWC) communities develop in the Mediterranean Sea. We investigated the CWC community of the two westernmost submarine canyons of the Gulf of Lions canyon system: the Cap de Creus Canyon (CCC) and Lacaze-Duthiers Canyon (LDC). Coral associations have been studied through video material recorded by means of a manned submersible and a remotely operated vehicle. Video transects have been conducted and analyzed in order to obtain information on (1) coral bathymetric distribution and density patterns, (2) size structure of coral populations, and (3) coral colony position with respect to the substrate. Madrepora oculata was the most abundant CWC in both canyons, while Lophelia pertusa and Dendrophyllia cornigera mostly occurred as isolated colonies or in small patches. An important exception was detected in a vertical cliff in LDC where a large L. pertusa framework was documented. This is the first record of such an extended L. pertusa framework in the Mediterranean Sea. In both canyons coral populations were dominated by medium and large colonies, but the frequent presence of small-sized colonies also indicate active recruitment. The predominant coral orientation (90° and 135°) is probably driven by the current regime as well as by the sediment load transported by the current flows. In general, no clear differences were observed in the abundance and in the size structure of the CWC populations between CCC and LDC, despite large differences in particulate matter between canyons.

Description: As a major source for atmospheric CO2, the Peruvian upwelling region exhibits strong variability in surface fCO2 on short spatial and temporal scales. Understanding the physical processes driving the strong variability is of fundamental importance for constraining the effect of marine emissions from upwelling regions on the global CO2 budget. In this study, a frontal decay on length scales of 𝒪(10 km) was observed off the Peruvian coast following a pronounced decrease in down-frontal (equatorward) wind speed with a time lag of 9 h. Simultaneously, the sea-to-air flux of CO2 on the inshore (cold) side of the front dropped from up to 80 to 10 mmol m−2 day−1, while the offshore (warm) side of the front was constantly outgassing at a rate of 10–20 mmol m−2 day−1. Based on repeated ship transects the decay of the front was observed to occur in two phases. The first phase was characterized by a development of coherent surface temperature anomalies which gained in amplitude over 6–9 h. The second phase was characterized by a disappearance of the surface temperature front within 6 h. Submesoscale mixed-layer instabilities were present but seem too slow to completely remove the temperature gradient in this short time period. Dynamics such as a pressure-driven gravity current appear to be a likely mechanism behind the evolution of the front.

Description: A major role in regulation of global methane fluxes has been attributed to the process of anaerobic oxidation of methane (AOM), which is performed by consortia of methanotrophic archaea and sulfate reducing bacteria. An important question remains how these energy limited, slow growing Microorganisms with generation times of 3–7 months respond to rapid natural variations in methane fluxes at cold seeps. We used an experimental flow-through column system filled with cold seep sediments naturally enriched in methanotrophic communities, to test their responses to short-term variations in methane and sulfate fluxes. At stable methane and sulfate concentrations of ∼2 mM and 28 mM, respectively, we measured constant rates of AOM and sulfate reduction (SR) for up to 160 days of incubation. When percolated with methane-free medium, the anaerobic methanotrophs ceased to produce sulfide. After a starvation phase of 40 days, the addition of methane restored former AOM and SR rates immediately. At methane concentrations between 0–2.3 mM we measured a linear correlation between methane availability, AOM and SR. At constant fluid flow velocities of 30 m yr−1, ca. 50% of the methane was consumed by the anaerobic methanotrophic (ANME) population at all concentrations tested. Reducing the sulfate concentration from 28 to 1 mM, a decrease in AOM and SR by 50% was observed, and 45% of the methane was consumed. Hence, the marine anaerobic methanotrophs(ANME) are capable of oxidizing substantial amounts of methane over a wide and variable range of fluxes of the reaction educts.

Description: The island of South Georgia is situated in the iron (Fe) depleted Antarctic Circumpolar Current of the Southern Ocean. Iron emanating from its shelf system fuels large phytoplankton blooms downstream of the island, but the actual supply mechanisms are unclear. To address this we present the first inventory of Fe, manganese (Mn) and aluminium (Al) in shelf sediments, pore waters and the water column in the vicinity of South Georgia, alongside data on zooplankton-mediated Fe cycling processes. The seafloor sediments were the main particulate Fe source to shelf bottom waters as indicated by Fe / Mn and Fe / Al ratios for shelf sediments and suspended particles in the water column. Less than 1 % of the total particulate Fe pool was leachable surface adsorbed (labile) Fe, and therefore potentially available to organisms. Pore waters formed the primary dissolved Fe (DFe) source to shelf bottom waters supplying 0.1–4 μmol DFe m−2 d−1. However, only 0.41 ± 0.26 μmol DFe m−2 d−1 was transferred to the surface mixed layer by vertical diffusive and advective mixing. Other trace metal sources to surface waters included glacial flour released by melting glaciers and zooplankton excretion processes. On average 6.5 ± 8.2 μmol m−2 d−1 of labile particulate Fe was supplied to the surface mixed layer via krill faecal pellets, with further DFe released by krill at around 1.1 ± 2.2 μmol m−2 d−1. The faecal pellets released by krill constituted of seafloor derived lithogenic material and settled algae debris, in addition to freshly ingested suspended phytoplankton specimen. The phytoplankton Fe requirement in the blooms ca. 1250 km downstream the island of South Georgia was 0.33 ± 0.11 μmol m−2 d−1, with the DFe supply by horizontal/vertical mixing, deep winter mixing and via aeolian dust estimated as ~ 0.12 μmol m−2 d−1. We suggest that additionally required DFe was provided through recycling of biogenically stored Fe following luxury Fe uptake by phytoplankton on the Fe rich shelf. This process would allow Fe to be retained in the surface mixed layer of waters downstream of South Georgia through continuous recycling and biological uptake, and facilitate the large scale blooms.

Description: In estuarine coastal systems such as the Baltic Sea, mussels suffer from low salinity which limits their distribution. Anthropogenic climate change is expected to cause further desalination which will lead to local extinctions of mussels in the low saline areas. It is commonly accepted that mussel distribution is limited by osmotic stress. However, along the salinity gradient environmental conditions for biomineralization are successively becoming more adverse as a result of reduced [Ca2+] and dissolved inorganic carbon (CT) availability. In larvae, calcification is an essential process starting during early development with formation of the prodissoconch I (PD I) shell which is completed under optimal conditions within 2 days. Experimental manipulations of seawater [Ca2+] start to impair PD I formation in Mytilus larvae at concentrations below 3 mM, which corresponds to conditions present in the Baltic at salinities below 8 g kg-1. In addition, lowering dissolved inorganic carbon to critical concentrations (〈 1 mM) similarly affected PD I size which was well correlated with calculated ΩAragonite and [Ca2+][HCO3-]/[H+] in all treatments. Comparing results for larvae from the western Baltic with a population from the central Baltic revealed significantly higher tolerance of PD I formation to lowered [Ca2+] and [Ca2+][HCO3-]/[H+] in the low saline adapted population. This may result from genetic adaptation to the more adverse environmental conditions prevailing in the low saline areas of the Baltic. The combined effects of lowered [Ca2+] and adverse carbonate chemistry represent major limiting factors for bivalve calcification and can thereby contribute to distribution limits of mussels in the Baltic Sea.

Description: The response of the marine carbon cycle to changes in atmospheric CO2 concentrations will be determined, in part, by the relative response of calcifying and non-calcifying organisms to global change. Planktonic foraminifera are responsible for a quarter or more of global carbonate production, therefore understanding the sensitivity of calcification in these organisms to environmental change is critical. Despite this, there remains little consensus as to whether, or to what extent, chemical and physical factors affect foraminiferal calcification. To address this, we directly test the effect of multiple controls on calcification in culture experiments and core-top measurements of Globigerinoides ruber. We find that two factors, body size and the carbonate system, strongly influence calcification intensity in life, but that exposure to corrosive bottom waters can overprint this signal post mortem. Using a simple model for the addition of calcite through ontogeny, we show that variable body size between and within datasets could complicate studies that examine environmental controls on foraminiferal shell weight. In addition, we suggest that size could ultimately play a role in determining whether calcification will increase or decrease with acidification. Our models highlight that knowledge of the specific morphological and physiological mechanisms driving ontogenetic change in calcification in different species will be critical in predicting the response of foraminiferal calcification to future change in atmospheric pCO2.

Description: The distribution of dissolved aluminium (dAl) in the water column of the North Atlantic and Labrador Sea was studied along GEOTRACES section GA01 to unravel the sources and sinks of this element. Surface water dAl concentrations were low (median of 2.5 nM) due to low aerosol deposition and removal by phytoplankton. However, surface water dAl concentrations were enhanced on the Iberian and Greenland shelves (up to 30.9 nM) due to continental inputs (rivers, glacial flour and ice melt). A negative correlation was observed between dAl in surface waters and primary production, phytoplankton community structure and biogenic opal production. The abundance of diatoms exerted a significant (p 〈 0.01) control on the surface particulate Al (pAl) to dAl ratios by decreasing dAl levels and increasing pAl levels. Dissolved Al concentrations generally increased with depth and correlated strongly with silicate (R2 〉 0.76) west of the Iberian Basin, suggesting net release of dAl at depth during remineralization of sinking biogenic opal containing particles. Enrichment of dAl at near-bottom depths was observed due to resuspension of sediments near the sediment-water interface. The highest dAl (up to 38.7 nM) concentrations were observed in Mediterranean Overflow Waters which act as a major source of dAl to mid depth waters of the eastern North Atlantic. This study clearly shows that the vertical and lateral distribution of dAl in the North Atlantic differs when compared to other regions of the North Atlantic and global ocean due to the large spatial differences both in the main source of Al, atmospheric deposition, and the main sink for Al, particle scavenging, between different oceanic regions.

Description: Autonomous underwater vehicles (AUVs) offer unique possibilities for exploring the deep seafloor in high resolution over large areas. We highlight the results from AUV-based multibeam echosounder (MBES) bathymetry / backscatter and digital optical imagery from the DISCOL area acquired during research cruise SO242 in 2015. AUV bathymetry reveals a morphologically complex seafloor with rough terrain in seamount areas and low-relief variations in sedimentary abyssal plains which are covered in Mn-nodules. Backscatter provides valuable information about the seafloor type and particularly about the influence of Mn-nodules on the response of the transmitted acoustic signal. Primarily, Mn-nodule abundances were determined by means of automated nodule detection on AUV seafloor imagery and nodule metrics such as nodules m−2 were calculated automatically for each image allowing further spatial analysis within GIS in conjunction with the acoustic data. AUV-based backscatter was clustered using both raw data and corrected backscatter mosaics. In total, two unsupervised methods and one machine learning approach were utilized for backscatter classification and Mn-nodule predictive mapping. Bayesian statistical analysis was applied to the raw backscatter values resulting in six acoustic classes. In addition, Iterative Self-Organizing Data Analysis (ISODATA) clustering was applied to the backscatter mosaic and its statistics (mean, mode, 10th, and 90th quantiles) suggesting an optimum of six clusters as well. Part of the nodule metrics data was combined with bathymetry, bathymetric derivatives and backscatter statistics for predictive mapping of the Mn-nodule density using a Random Forest classifier. Results indicate that acoustic classes, predictions from Random Forest model and image-based nodule metrics show very similar spatial distribution patterns with acoustic classes hence capturing most of the fine-scale Mn-nodule variability. Backscatter classes reflect areas with homogeneous nodule density. A strong influence of mean backscatter, fine scale BPI and concavity of the bathymetry on nodule prediction is seen. These observations imply that nodule densities are generally affected by local micro-bathymetry in a way that is not yet fully understood. However, it can be concluded that the spatial occurrence of Mn-covered areas can be sufficiently analysed by means of acoustic classification and multivariate predictive mapping allowing to determine the spatial nodule density in a much more robust way than previously possible.

Description: Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the "global carbon budget" – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, respectively, while emissions from land-use change (ELUC), mainly deforestation, are based on land-cover change data and bookkeeping models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2007–2016), EFF was 9.4 ± 0.5 GtC yr−1, ELUC 1.3 ± 0.7 GtC yr−1, GATM 4.7 ± 0.1 GtC yr−1, SOCEAN 2.4 ± 0.5 GtC yr−1, and SLAND 3.0 ± 0.8 GtC yr−1, with a budget imbalance BIM of 0.6 GtC yr−1 indicating overestimated emissions and/or underestimated sinks. For year 2016 alone, the growth in EFF was approximately zero and emissions remained at 9.9 ± 0.5 GtC yr−1. Also for 2016, ELUC was 1.3 ± 0.7 GtC yr−1, GATM was 6.1 ± 0.2 GtC yr−1, SOCEAN was 2.6 ± 0.5 GtC yr−1, and SLAND was 2.7 ± 1.0 GtC yr−1, with a small BIM of −0.3 GtC. GATM continued to be higher in 2016 compared to the past decade (2007–2016), reflecting in part the high fossil emissions and the small SLAND consistent with El Niño conditions. The global atmospheric CO2 concentration reached 402.8 ± 0.1 ppm averaged over 2016. For 2017, preliminary data for the first 6–9 months indicate a renewed growth in EFF of +2.0 % (range of 0.8 to 3.0 %) based on national emissions projections for China, USA, and India, and projections of gross domestic product (GDP) corrected for recent changes in the carbon intensity of the economy for the rest of the world. This living data update documents changes in the methods and data sets used in this new global carbon budget compared with previous publications of this data set (Le Quéré et al., 2016, 2015b, a, 2014, 2013). All results presented here can be downloaded from https://doi.org/10.18160/GCP-2017 (GCP, 2017).

Description: A long-term mean flow field for the subpolar North Atlantic region with a horizontal resolution of approximately 25km is created by gridding Argo-derived velocity vectors using two different topography-following interpolation schemes. The 10-day float displacements in the typical drift depths of 1000 to 1500m represent the flow in the Labrador Sea Water density range. Both mapping algorithms separate the flow field into potential vorticity (PV) conserving, i.e., topography-following contribution and a deviating part, which we define as the eddy contribution. To verify the significance of the separation, we compare the mean flow and the eddy kinetic energy (EKE), derived from both mapping algorithms, with those obtained from multiyear mooring observations. The PV-conserving mean flow is characterized by stable boundary currents along all major topographic features including shelf breaks and basin-interior topographic ridges such as the Reykjanes Ridge or the Rockall Plateau. Mid-basin northward advection pathways from the northeastern Labrador Sea into the Irminger Sea and from the Mid-Atlantic Ridge region into the Iceland Basin are well-resolved. An eastward flow is present across the southern boundary of the subpolar gyre near 52°N, the latitude of the Charlie Gibbs Fracture Zone (CGFZ). The mid-depth EKE field resembles most of the satellite-derived surface EKE field. However, noticeable differences exist along the northward advection pathways in the Irminger Sea and the Iceland Basin, where the deep EKE exceeds the surface EKE field. Further, the ratio between mean flow and the square root of the EKE, the Peclet number, reveals distinct advection-dominated regions as well as basin-interior regimes in which mixing is prevailing.

Description: Past ocean temperatures and salinities can be approximated from combined stable oxygen isotopes (δ18O) and Mg ∕ Ca measurements in fossil foraminiferal tests with varying success. To further refine this approach, we collected living planktic foraminifers by net sampling and pumping of sea surface water from the Caribbean Sea, the eastern Gulf of Mexico and the Florida Straits. Analyses of δ18O and Mg ∕ Ca in eight living planktic species (Globigerinoides sacculifer, Orbulina universa, Neogloboquadrina dutertrei, Pulleniatina obliquiloculata, Globorotalia menardii, Globorotalia ungulata, Globorotalia truncatulinoides and Globorotalia tumida) were compared to measured in situ properties of the ambient seawater (temperature, salinity and δ18Oseawater) and fossil tests of underlying surface sediments. “Vital effects” such as symbiont activity and test growth cause δ18O disequilibria with respect to the ambient seawater and a large scatter in foraminiferal Mg ∕ Ca. Overall, ocean temperature is the most prominent environmental influence on δ18Ocalcite and Mg ∕ Ca. Enrichment of the heavier 18O isotope in living specimens below the mixed layer and in fossil tests is clearly related to lowered in situ temperatures and gametogenic calcification. Mg ∕ Ca-based temperature estimates of G. sacculifer indicate seasonal maximum accumulation rates on the seafloor in early spring (March) at Caribbean stations and later in the year (May) in the Florida Straits, related to the respective mixed layer temperatures of ∼26 ∘C. Notably, G. sacculifer reveals a weak positive linear relationship between foraminiferal derived δ18Oseawater estimates and both measured in situ δ18Oseawater and salinity. Our results affirm the applicability of existing δ18O and Mg ∕ Ca calibrations for the reconstruction of past ocean temperatures and δ18Oseawater reflecting salinity due to the convincing accordance of proxy data in both living and fossil foraminifers, and in situ environmental parameters. Large vital effects and seasonally varying proxy signals, however, need to be taken into account.

Description: We propose a satellite mission that uses a near-nadir Ka-band Doppler radar to measure surface currents, ice drift and ocean waves at spatial scales of 40 km and more, with snapshots at least every day for latitudes 75 to 82°, and every few days for other latitudes. The use of incidence angles of 6 and 12° allows for measurement of the directional wave spectrum, which yields accurate corrections of the wave-induced bias in the current measurements. The instrument's design, an algorithm for current vector retrieval and the expected mission performance are presented here. The instrument proposed can reveal features of tropical ocean and marginal ice zone (MIZ) dynamics that are inaccessible to other measurement systems, and providing global monitoring of the ocean mesoscale that surpasses the capability of today's nadir altimeters. Measuring ocean wave properties has many applications, including examining wave–current interactions, air–sea fluxes, the transport and convergence of marine plastic debris and assessment of marine and coastal hazards.

Description: Gel particles, such as transparent exopolymer particles (TEP) and Coomassie stainable particles (CSP), are important organic components in the sea surface microlayer (SML). Here, we present results on the effect of different wind speeds on the accumulation and size distribution of TEP and CSP during a wind wave channel experiment in the Aeolotron. Total areas of TEP (TEPSML) and CSP (CSPSML) in the surface microlayer were exponentially related to wind speed. At wind speeds  〈  6 m s−1, accumulation of TEPSML and CSPSML occurred, decreasing at wind speeds of  〉  8 m s−1. Wind speeds  〉  8 m s−1 also significantly altered the size distribution of TEPSML in the 2–16 µm size range towards smaller sizes. The response of the CSPSML size distribution to wind speed varied through time depending on the biogenic source of gels. Wind speeds  〉  8 m s−1 decreased the slope of CSPSML size distribution significantly in the absence of autotrophic growth. For the slopes of TEP and CSP size distribution in the bulk water, no significant difference was observed between high and low wind speeds. Changes in spectral slopes between high and low wind speed were higher for TEPSML than for CSPSML, indicating that the impact of wind speed on size distribution of gel particles in the SML may be more pronounced for TEP than for CSP, and that CSPSML are less prone to aggregation during the low wind speeds. Addition of an E. huxleyi culture resulted in a higher contribution of submicron gels (0.4–1 µm) in the SML at higher wind speed ( 〉  6 m s−1), indicating that phytoplankton growth may potentially support the emission of submicron gels with sea spray aerosol.

Description: In the ocean, sinking of particulate organic matter (POM) drives carbon export from the euphotic zone and supplies nutrition to mesopelagic communities, the feeding and degradation activities of which in turn lead to export flux attenuation. Oxygen minimum zones (OMZs) with suboxic water layers (〈 5 μmol O2 kg−1) show a lower carbon flux attenuation compared to well oxygenated waters (〉 100 μmol O2 kg−1), supposedly due to reduced heterotrophic activity. This study focuses on sinking particle fluxes through hypoxic mesopelagic waters (〈 60 % μmol O2 kg−1); these represent ~ 100-times more ocean volume globally compared to suboxic waters, but have less been studied. Particle export fluxes and attenuation coefficients were determined in the Eastern Tropical North Atlantic (ETNA) using two surface tethered drifting sediment trap arrays with 7 trapping depths located between 100 and 600 m. Data on particulate matter fluxes were fitted to the normalized power function Fz = F100 (z/100)−b, with F100 being the flux at a depth (z) of 100 m and b being the attenuation coefficient. Higher b-values suggest stronger flux attenuation and are influenced by factors such as faster degradation at higher temperatures. In this study, b-values of organic carbon fluxes varied between 0.74 and 0.80 and were in the intermediate range of previous reports, but lower than expected from seawater temperatures within the upper 500 m. During this study, highest b-values were determined for fluxes of particulate hydrolysable amino acids (PHAA), followed by particulate organic phosphorus (POP), nitrogen (PN), carbon (POC), chlorophyll a, and transparent exopolymer particles (TEP), pointing to a sequential degradation of organic matter components during sinking. Our study suggests that in addition to oxygen concentration, organic matter composition co-determines transfer efficiency through the mesopelagial. The magnitude of future carbon export fluxes may therefore also depend on how organic matter quality in the surface ocean changes under influence of warming, acidification, and enhanced stratification.

Description: Well water level changes associated with magmatic unrest can be interpreted as a result of pore pressure changes in the aquifer due to crustal deformation, and so could provide constraints on the subsurface processes causing this strain. We use finite element analysis to demonstrate the response of aquifers to volumetric strain induced by pressurized magma reservoirs. Two different aquifers are invoked – an unconsolidated pyroclastic deposit and a vesicular lava flow – and embedded in an impermeable crust, overlying a magma chamber. The time-dependent, fully coupled models simulate crustal deformation accompanying chamber pressurization and the resulting hydraulic head changes as well as flow through the porous aquifer, i.e. porous flow. The simulated strain leads to centimetres (pyroclastic aquifer) to metres (lava flow aquifer) of hydraulic head changes; both strain and hydraulic head change with time due to substantial porous flow in the hydrological system. Well level changes are particularly sensitive to chamber volume, shape and pressurization strength, followed by aquifer permeability and the phase of the pore fluid. The depths of chamber and aquifer, as well as the aquifer's Young's modulus also have significant influence on the hydraulic head signal. While source characteristics, the distance between chamber and aquifer and the elastic stratigraphy determine the strain field and its partitioning, flow and coupling parameters define how the aquifer responds to this strain and how signals change with time. We find that generic analytical models can fail to capture the complex pre-eruptive subsurface mechanics leading to strain-induced well level changes, due to aquifer pressure changes being sensitive to chamber shape and lithological heterogeneities. In addition, the presence of a pore fluid and its flow have a significant influence on the strain signal in the aquifer and are commonly neglected in analytical models. These findings highlight the need for numerical models for the interpretation of observed well level signals. However, simulated water table changes do indeed mirror volumetric strain, and wells are therefore a valuable addition to monitoring systems that could provide important insights into pre-eruptive dynamics.

Description: The meridional Ekman volume, heat, and salt transport across two trans-Atlantic sections near 14.5° N and 11° S were estimated using in situ observations, wind products, and model data. A meridional ageostrophic velocity was obtained as the difference between the directly measured total velocity and the geostrophic velocity derived from observations. Interpreting the section mean ageostrophy to be the result of an Ekman balance, the meridional Ekman transport of 6. 2 ± 2. 3 Sv northward at 14.5° N and 11. 7 ± 2. 1 Sv southward at 11° S is estimated. The integration uses the top of the pycnocline as an approximation for the Ekman depth, which is on average about 20 m deeper than the mixed layer depth. The Ekman transport estimated based on the velocity observations agrees well with the predictions from in situ wind stress data of 6. 7 ± 3. 5 Sv at 14.5° N and 13. 6 ± 3. 3 Sv at 11° S. The meridional Ekman heat and salt fluxes calculated from sea surface temperature and salinity data or from high-resolution temperature and salinity profile data differ only marginally. The errors in the Ekman heat and salt flux calculation were dominated by the uncertainty of the Ekman volume transport estimates.

Description: Halogenated very short-lived substances (VSLSs) are naturally produced in the ocean and emitted to the atmosphere. When transported to the stratosphere, these compounds can have a significant influence on the ozone layer and climate. During a research cruise on RV Sonne in the subtropical and tropical west Indian Ocean in July and August 2014, we measured the VSLSs, methyl iodide (CH3I) and for the first time bromoform (CHBr3) and dibromomethane (CH2Br2), in surface seawater and the marine atmosphere to derive their emission strengths. Using the Lagrangian particle dispersion model FLEXPART with ERA-Interim meteorological fields, we calculated the direct contribution of observed VSLS emissions to the stratospheric halogen burden during the Asian summer monsoon. Furthermore, we compare the in situ calculations with the interannual variability of transport from a larger area of the west Indian Ocean surface to the stratosphere for July 2000–2015. We found that the west Indian Ocean is a strong source for CHBr3 (910 pmol m−2 h−1), very strong source for CH2Br2 (930 pmol m−2 h−1), and an average source for CH3I (460 pmol m−2 h−1). The atmospheric transport from the tropical west Indian Ocean surface to the stratosphere experiences two main pathways. On very short timescales, especially relevant for the shortest-lived compound CH3I (3.5 days lifetime), convection above the Indian Ocean lifts oceanic air masses and VSLSs towards the tropopause. On a longer timescale, the Asian summer monsoon circulation transports oceanic VSLSs towards India and the Bay of Bengal, where they are lifted with the monsoon convection and reach stratospheric levels in the southeastern part of the Asian monsoon anticyclone. This transport pathway is more important for the longer-lived brominated compounds (17 and 150 days lifetime for CHBr3 and CH2Br2). The entrainment of CHBr3 and CH3I from the west Indian Ocean to the stratosphere during the Asian summer monsoon is lower than from previous cruises in the tropical west Pacific Ocean during boreal autumn and early winter but higher than from the tropical Atlantic during boreal summer. In contrast, the projected CH2Br2 entrainment was very high because of the high emissions during the west Indian Ocean cruise. The 16-year July time series shows highest interannual variability for the shortest-lived CH3I and lowest for the longest-lived CH2Br2. During this time period, a small increase in VSLS entrainment from the west Indian Ocean through the Asian monsoon to the stratosphere is found. Overall, this study confirms that the subtropical and tropical west Indian Ocean is an important source region of halogenated VSLSs, especially CH2Br2, to the troposphere and stratosphere during the Asian summer monsoon.

Description: Treatment of the underwater light field in ocean biogeochemical models has been attracting increasing interest, with some models moving towards more complex parameterisations. We conduct a simple sensitivity study of a typical, highly simplified parameterisation. In our study, we vary the phytoplankton light attenuation parameter over a range constrained by data during both pre-industrial equilibrated and future climate scenario RCP8.5. In equilibrium, lower light attenuation parameters (weaker self-shading) shift net primary production (NPP) towards the high latitudes, while higher values of light attenuation (stronger shelf-shading) shift NPP towards the low latitudes. Climate forcing magnifies this relationship through changes in the distribution of nutrients both within and between ocean regions. Where and how NPP responds to climate forcing can determine the magnitude and sign of global NPP trends in this high CO2 future scenario. Ocean oxygen is particularly sensitive to parameter choice. Under higher CO2 concentrations, two simulations establish a strong biogeochemical feedback between the Southern Ocean and low-latitude Pacific that highlights the potential for regional teleconnection. Our simulations serve as a reminder that shifts in fundamental properties (e.g. light attenuation by phytoplankton) over deep time have the potential to alter global biogeochemistry.

Description: Occupying about 14 % of the world's surface, the Southern Ocean plays a fundamental role in ocean and atmosphere circulation, carbon cycling and Antarctic ice-sheet dynamics. Unfortunately, high interannual variability and a dearth of instrumental observations before the 1950s limits our understanding of how marine–atmosphere–ice domains interact on multi-decadal timescales and the impact of anthropogenic forcing. Here we integrate climate-sensitive tree growth with ocean and atmospheric observations on southwest Pacific subantarctic islands that lie at the boundary of polar and subtropical climates (52–54° S). Our annually resolved temperature reconstruction captures regional change since the 1870s and demonstrates a significant increase in variability from the 1940s, a phenomenon predating the observational record. Climate reanalysis and modelling show a parallel change in tropical Pacific sea surface temperatures that generate an atmospheric Rossby wave train which propagates across a large part of the Southern Hemisphere during the austral spring and summer. Our results suggest that modern observed high interannual variability was established across the mid-twentieth century, and that the influence of contemporary equatorial Pacific temperatures may now be a permanent feature across the mid- to high latitudes.

Description: The outer western Crimean shelf of the Black Sea is a natural laboratory to investigate effects of stable oxic versus varying hypoxic conditions on seafloor biogeochemical processes and benthic community structure. Bottom-water oxygen concentrations ranged from normoxic (175 μmol O2 L−1) and hypoxic (〈 63 μmol O2 L−1) or even anoxic/sulfidic conditions within a few kilometers' distance. Variations in oxygen concentrations between 160 and 10 μmol L−1 even occurred within hours close to the chemocline at 134 m water depth. Total oxygen uptake, including diffusive as well as fauna-mediated oxygen consumption, decreased from 15 mmol m−2 d−1 on average in the oxic zone, to 7 mmol m−2 d−1 on average in the hypoxic zone, correlating with changes in macrobenthos composition. Benthic diffusive oxygen uptake rates, comprising respiration of microorganisms and small meiofauna, were similar in oxic and hypoxic zones (on average 4.5 mmol m−2 d−1), but declined to 1.3 mmol m−2 d−1 in bottom waters with oxygen concentrations below 20 μmol L−1. Measurements and modeling of porewater profiles indicated that reoxidation of reduced compounds played only a minor role in diffusive oxygen uptake under the different oxygen conditions, leaving the major fraction to aerobic degradation of organic carbon. Remineralization efficiency decreased from nearly 100 % in the oxic zone, to 50 % in the oxic–hypoxic zone, to 10 % in the hypoxic–anoxic zone. Overall, the faunal remineralization rate was more important, but also more influenced by fluctuating oxygen concentrations, than microbial and geochemical oxidation processes.

Description: In this study we present a comparative quantification of CaCO3 production rates by rhodolith-forming coralline red algal communities situated in high polar latitudes and assess which environmental parameters control these production rates. The present rhodoliths act as ecosystem engineers, and their carbonate skeletons provide an important ecological niche to a variety of benthic organisms. The settings are distributed along the coasts of the Svalbard archipelago, being Floskjeret (78◦180N) in Isfjorden, Krossfjorden (79◦080N) at the eastern coast of Haakon VII Land, Mosselbukta (79◦530N) at the eastern coast of Mosselhalvøya, and Nordkappbukta (80◦310N) at the northern coast of Nordaustlandet. All sites feature Arctic climate and strong seasonality.

Description: The impact of changes in incoming solar irradiance on stratospheric ozone abundances should be included in climate simulations to aid in capturing the atmospheric response to solar cycle variability. This study presents the first systematic comparison of the representation of the 11-year solar cycle ozone response (SOR) in chemistry–climate models (CCMs) and in pre-calculated ozone databases specified in climate models that do not include chemistry, with a special focus on comparing the recommended protocols for the Coupled Model Intercomparison Project Phase 5 and Phase 6 (CMIP5 and CMIP6). We analyse the SOR in eight CCMs from the Chemistry–Climate Model Initiative (CCMI-1) and compare these with results from three ozone databases for climate models: the Bodeker Scientific ozone database, the SPARC/Atmospheric Chemistry and Climate (AC&C) ozone database for CMIP5 and the SPARC/CCMI ozone database for CMIP6. The peak amplitude of the annual mean SOR in the tropical upper stratosphere (1–5hPa) decreases by more than a factor of 2, from around 5 to 2%, between the CMIP5 and CMIP6 ozone databases. This substantial decrease can be traced to the CMIP5 ozone database being constructed from a regression model fit to satellite and ozonesonde measurements, while the CMIP6 database is constructed from CCM simulations. The SOR in the CMIP6 ozone database therefore implicitly resembles the SOR in the CCMI-1 models. The structure in latitude of the SOR in the CMIP6 ozone database and CCMI-1 models is considerably smoother than in the CMIP5 database, which shows unrealistic sharp gradients in the SOR across the middle latitudes owing to the paucity of long-term ozone measurements in polar regions. The SORs in the CMIP6 ozone database and the CCMI-1 models show a seasonal dependence with enhanced meridional gradients at mid- to high latitudes in the winter hemisphere. The CMIP5 ozone database does not account for seasonal variations in the SOR, which is unrealistic. Sensitivity experiments with a global atmospheric model without chemistry (ECHAM6.3) are performed to assess the atmospheric impacts of changes in the representation of the SOR and solar spectral irradiance (SSI) forcing between CMIP5 and CMIP6. The larger amplitude of the SOR in the CMIP5 ozone database compared to CMIP6 causes a likely overestimation of the modelled tropical stratospheric temperature response between 11-year solar cycle minimum and maximum by up to 0.55K, or around 80% of the total amplitude. This effect is substantially larger than the change in temperature response due to differences in SSI forcing between CMIP5 and CMIP6. The results emphasize the importance of adequately representing the SOR in global models to capture the impact of the 11-year solar cycle on the atmosphere. Since a number of limitations in the representation of the SOR in the CMIP5 ozone database have been identified, we recommend that CMIP6 models without chemistry use the CMIP6 ozone database and the CMIP6 SSI dataset to better capture the climate impacts of solar variability. The SOR coefficients from the CMIP6 ozone database are published with this paper.

Description: The radiative effects induced by the zonally asymmetric part of the ozone field have been shown to significantly change the temperature of the NH winter polar cap, and correspondingly the strength of the polar vortex. In this paper, we aim to understand the physical processes behind these effects using the National Center for Atmospheric Research (NCAR)'s Whole Atmosphere Community Climate Model, run with 1960s ozone-depleting substances and greenhouse gases. We find a mid-winter polar vortex influence only when considering the quasi-biennial oscillation (QBO) phases separately, since ozone waves affect the vortex in an opposite manner. Specifically, the emergence of a midlatitude QBO signal is delayed by 1–2 months when radiative ozone-wave effects are removed. The influence of ozone waves on the winter polar vortex, via their modulation of shortwave heating, is not obvious, given that shortwave heating is largest during fall, when planetary stratospheric waves are weakest. Using a novel diagnostic of wave 1 temperature amplitude tendencies and a synoptic analysis of upward planetary wave pulses, we are able to show the chain of events that lead from a direct radiative effect on weak early fall upward-propagating planetary waves to a winter polar vortex modulation. We show that an important stage of this amplification is the modulation of individual wave life cycles, which accumulate during fall and early winter, before being amplified by wave–mean flow feedbacks. We find that the evolution of these early winter upward planetary wave pulses and their induced stratospheric zonal mean flow deceleration is qualitatively different between QBO phases, providing a new mechanistic view of the extratropical QBO signal. We further show how these differences result in opposite radiative ozone-wave effects between east and west QBOs.

Description: The kinetics of particulate organic carbon (POC) mineralization in marine surface sediments is not well constrained. This creates considerable uncertainties when benthic processes are considered in global biogeochemical or Earth system circulation models to simulate climate-ocean interactions and biogeochemical tracers in the ocean. In an attempt to improve our understanding of the rate and depth distribution of organic carbon mineralization in bioturbated (0–10 cm) sediments, we parameterized a 1-D diagenetic model that simulates the reactivity of three discrete POC pools at global scale (a "multi-G" model). The rate constants of the three reactive classes (highly reactive, reactive, refractory) are fixed and determined to be 70 yr−1, 0.5 yr−1, and ~0.001 yr−1, respectively, based on the Martin curve model for pelagic POC degradation. In contrast to previous approaches, the reactivity of the organic material degraded in the seafloor is continuous with, and set by, the apparent reactivity of material sinking through the water column. The model is able to simulate a global database (185 stations) of benthic oxygen and nitrate fluxes across the sediment-water interface in addition to porewater oxygen and nitrate distributions and organic carbon burial efficiencies. It is further consistent with degradation experiments of fresh phytoplankton. We propose that an important yet mostly overlooked consideration in previous upscaling approaches is the proportion of the relative reactive POC classes reaching the seafloor in addition to their reactivity. The approach presented is applicable to both steady-state and non-steady state scenarios, and links POC degradation kinetics in sedimentary environments to water depth and the POC rain rate to the seafloor.

Description: The Atlantic Meridional Overturning Circulation (AMOC) is analyzed by applying a box inverse model to hydrographic data from transatlantic sections along 14.5°N, occupied in 1989 and 2013, and along 24.5°N, occupied in 1992 and 2015. Direct comparison of water mass properties among the different realizations at the respective latitudes shows that the Antarctic Intermediate Water (AAIW) became warmer and saltier at 14.5°N, and the densest Antarctic Bottom Water became lighter, while the North Atlantic Deep Water freshened at both latitudes. The inverse solution shows that the intermediate layer transport at 14.5°N was also markedly weaker in 2013 than in 1989, indicating that the AAIW property changes at this latitude may be related to changes in the circulation. The inverse solution was validated using the RAPID and MOVE array data, and the GECCO2 ocean state estimate. Comparison among these datasets indicates that the AMOC has not significantly weakened over the past 2 decades at both latitudes. Sensitivity tests of the inverse solution suggest that the overturning structure and heat transport across the 14.5°N section are sensitive to the Ekman transport, while freshwater transport is sensitive to the transport-weighted salinity at the western boundary.

Description: We evaluated NOVAC (Network for Observation of Volcanic and Atmospheric Change) gas emission data from the 2015 eruption of Cotopaxi volcano (Ecuador) for BrO/SO2 molar ratios. Statistical analysis of the data revealed a conspicuous periodic pattern with a periodicity of about two weeks in a three month time series. While the time series is too short to rule out a chance recurrence of transient geological or meteorological events as a possible origin for the periodic signal, we nevertheless took this observation as a motivation to examine the influence of natural forcings with periodicities of around two weeks on volcanic gas emissions. One strong aspirant with such a periodicity are the Earth tides, which are thus central in this study. We present the BrO/SO2 data, analyse the reliability of the periodic signal, discuss a possible meteorological or eruption-induced origin of this signal, and compare the signal with the theoretical ground surface displacement pattern caused by the Earth tides. Central result is the observation of a significant correlation between the BrO/SO2 molar ratios with the North-South and vertical components of the calculated tide-induced surface displacement with correlation coefficients of 47 % and 36 %, respectively. From all other investigated parameters, only the correlation between the BrO/SO2 molar ratios and the relative humidity in the local atmosphere resulted in a comparable correlation coefficient of about 33 %.

Description: Numerous studies have been conducted on the effect of ocean acidification on calcifiers inhabiting nearshore benthic habitats, such as the blue mussel Mytilus edulis. The majority of these experiments was performed under stable CO2 partial pressure (pCO2), carbonate chemistry and oxygen (O2) levels, reflecting present or expected future open ocean conditions. Consequently, levels and variations occurring in coastal habitats, due to biotic and abiotic processes, were mostly neglected, even though these variations largely override global long-term trends. To highlight this hiatus and guide future research, state-of-the-art technologies were deployed to obtain high-resolution time series of pCO2 and [O2] on a mussel patch within a Zostera marina seagrass bed, in Kiel Bay (western Baltic Sea) in August and September 2013. Combining the in situ data with results of discrete sample measurements, a full seawater carbonate chemistry was derived using statistical models. An average pCO2 more than 50 % (~ 640 µatm) higher than current atmospheric levels was found right above the mussel patch. Diel amplitudes of pCO2 were large: 765 ± 310 (mean ± SD). Corrosive conditions for calcium carbonates (Ωarag and Ωcalc 〈 1) centered on sunrise were found, but the investigated habitat never experienced hypoxia throughout the study period. It is estimated that mussels experience conditions limiting calcification for 12–15 h per day, based on a regional calcium carbonate concentration physiological threshold. Our findings call for more extensive experiments on the impact of fluctuating corrosive conditions on mussels. We also stress the complexity of the interpretation of carbonate chemistry time series data in such dynamic coastal environments.

Description: Atmospheric carbon dioxide (CO2) levels continue to rise, increasing the risk of severe impacts on the Earth system, and on the ecosystem services that it provides. Artificial ocean alkalinization (AOA) is capable of reducing atmospheric CO2 concentrations and surface warming and addressing ocean acidification. Here, we simulate global and regional responses to alkalinity (ALK) addition (0.25 PmolALK yr−1) over the period 2020–2100 using the CSIRO-Mk3L-COAL Earth System Model, under high (Representative Concentration Pathway 8.5; RCP8.5) and low (RCP2.6) emissions. While regionally there are large changes in alkalinity associated with locations of AOA, globally we see only a very weak dependence on where and when AOA is applied. On a global scale, while we see that under RCP2.6 the carbon uptake associated with AOA is only ∼ 60 % of the total, under RCP8.5 the relative changes in temperature are larger, as are the changes in pH (140 %) and aragonite saturation state (170 %). The simulations reveal AOA is more effective under lower emissions, therefore the higher the emissions the more AOA is required to achieve the same reduction in global warming and ocean acidification. Finally, our simulated AOA for 2020–2100 in the RCP2.6 scenario is capable of offsetting warming and ameliorating ocean acidification increases at the global scale, but with highly variable regional responses.

Description: This study introduces the Systematic Correlation Matrix Evaluation (SCoMaE) method, a bottom–up approach which combines expert judgment and statistical information to systematically select transparent, nonredundant indicators for a comprehensive assessment of the state of the Earth system. The methods consists of two basic steps: (1) the calculation of a correlation matrix among variables relevant for a given research question and (2) the systematic evaluation of the matrix, to identify clusters of variables with similar behavior and respective mutually independent indicators. Optional further analysis steps include (3) the interpretation of the identified clusters, enabling a learning effect from the selection of indicators, (4) testing the robustness of identified clusters with respect to changes in forcing or boundary conditions, (5) enabling a comparative assessment of varying scenarios by constructing and evaluating a common correlation matrix, and (6) the inclusion of expert judgment, for example, to prescribe indicators, to allow for considerations other than statistical consistency. The example application of the SCoMaE method to Earth system model output forced by different CO2 emission scenarios reveals the necessity of reevaluating indicators identified in a historical scenario simulation for an accurate assessment of an intermediate–high, as well as a business-as-usual, climate change scenario simulation. This necessity arises from changes in prevailing correlations in the Earth system under varying climate forcing. For a comparative assessment of the three climate change scenarios, we construct and evaluate a common correlation matrix, in which we identify robust correlations between variables across the three considered scenarios.

Description: Magmatic sill intrusions into organic-rich sediments cause the release of thermogenic CH4 and CO2. Pore fluids from the Guaymas Basin (Gulf of California), a sedimentary basin with recent magmatic activity, were investigated to constrain the link between sill intrusions and fluid seepage as well as the timing of sill-induced hydrothermal activity. Sampling sites were close to a hydrothermal vent field at the northern rift axis and at cold seeps located up to 30km away from the rift. Pore fluids close to the active hydrothermal vent field showed a slight imprint by hydrothermal fluids and indicated a shallow circulation system transporting seawater to the hydrothermal catchment area. Geochemical data of pore fluids at cold seeps showed a mainly ambient diagenetic fluid composition without any imprint related to high temperature processes at greater depth. Seep communities at the seafloor were mainly sustained by microbial methane, which rose along pathways formed earlier by hydrothermal activity, driving the anaerobic oxidation of methane (AOM) and the formation of authigenic carbonates. Overall, our data from the cold seep sites suggest that at present, sill-induced hydrothermalism is not active away from the ridge axis, and the vigorous venting of hydrothermal fluids is restricted to the ridge axis. Using the sediment thickness above extinct conduits and carbonate dating, we calculated that deep fluid and thermogenic gas flow ceased 28 to 7kyr ago. These findings imply a short lifetime of hydrothermal systems, limiting the time of unhindered carbon release as suggested in previous modeling studies. Consequently, activation and deactivation mechanisms of these systems need to be better constrained for the use in climate modeling approaches.

Description: The relative importance of potential source and sink terms for bromoform (CHBr3) in the tropical Atlantic Ocean is investigated with a coupled physical-biogeochemical water column model. Bromoform production is either assumed to be linked to primary production or to phytoplankton losses; bromoform decay is treated as light dependent (photolysis), and in addition either vertically uniform, proportional to remineralisation or to nitrification. All experiments lead to the observed subsurface maximum of bromoform, corresponding to the subsurface phytoplankton biomass maximum. In the surface mixed layer, the concentration is set by entrainment from below, photolysis in the upper few meters and the outgassing to the atmosphere. The assumed bromoform production mechanism has only minor effects on the solution, but the various loss terms lead to significantly different bromoform concentrations below 200 m depth. The best agreement with observations is obtained when the bromoform decay is coupled to nitrification (parameterised by an inverse proportionality to the light field). Our model results reveal a pronounced seasonal cycle of bromoform outgassing, with a minimum in summer and a maximum in early winter, when the deepening surface mixed layer reaches down into the bromoform production layer

Description: The Mediterranean Sea is a semi-enclosed sea characterized by high salinities, temperatures and densities. The net evaporation exceeds the precipitation, driving an anti-estuarine circulation through the Strait of Gibraltar, contributing to very low nutrient concentrations. The Mediterranean Sea has an active overturning circulation, one shallow cell that communicates directly with the Atlantic Ocean, and two deep overturning cells, one in each of the two main basins. It is surrounded by populated areas and is thus sensitive to anthropogenic forcing. Several dramatic changes in the oceanographic and biogeochemical conditions have been observed during the past several decades, emphasizing the need to better monitor and understand the changing conditions and their drivers. During 2011 three oceanographic cruises were conducted in a coordinated fashion in order to produce baseline data of important physical and biogeochemical parameters that can be compared to historic data and be used as reference for future observational campaigns. In this article we provide information on the Mediterranean Sea oceanographic situation, and present a short review that will serve as background information for the special issue in Ocean Science on "Physical, chemical and biological oceanography of the Mediterranean Sea". An important contribution of this article is the set of figures showing the large-scale distributions of physical and chemical properties along the full length

Description: Recent studies on the impacts of ocean acidification on pelagic communities have identified changes in carbon to nutrient dynamics with related shifts in elemental stoichiometry. In principle, mesocosm experiments provide the opportunity of determining temporal dynamics of all relevant carbon and nutrient pools and, thus, calculating elemental budgets. In practice, attempts to budget mesocosm enclosures are often hampered by uncertainties in some of the measured pools and fluxes, in particular due to uncertainties in constraining air–sea gas exchange, particle sinking, and wall growth. In an Arctic mesocosm study on ocean acidification applying KOSMOS (Kiel Off-Shore Mesocosms for future Ocean Simulation), all relevant element pools and fluxes of carbon, nitrogen and phosphorus were measured, using an improved experimental design intended to narrow down the mentioned uncertainties. Water-column concentrations of particulate and dissolved organic and inorganic matter were determined daily. New approaches for quantitative estimates of material sinking to the bottom of the mesocosms and gas exchange in 48 h temporal resolution as well as estimates of wall growth were developed to close the gaps in element budgets. However, losses elements from the budgets into a sum of insufficiently determined pools were detected, and are principally unavoidable in mesocosm investigation. The comparison of variability patterns of all single measured datasets revealed analytic precision to be the main issue in determination of budgets. Uncertainties in dissolved organic carbon (DOC), nitrogen (DON) and particulate organic phosphorus (POP) were much higher than the summed error in determination of the same elements in all other pools. With estimates provided for all other major elemental pools, mass balance calculations could be used to infer the temporal development of DOC, DON and POP pools. Future elevated pCO2 was found to enhance net autotrophic community carbon uptake in two of the three experimental phases but did not significantly affect particle elemental composition. Enhanced carbon consumption appears to result in accumulation of dissolved organic carbon under nutrient-recycling summer conditions. This carbon over-consumption effect becomes evident from mass balance calculations, but was too small to be resolved by direct measurements of dissolved organic matter. Faster nutrient uptake by comparatively small algae at high CO2 after nutrient addition resulted in reduced production rates under future ocean CO2 conditions at the end of the experiment. This CO2 mediated shift towards smaller phytoplankton and enhanced cycling of dissolved matter restricted the development of larger phytoplankton, thus pushing the system towards a retention type food chain with overall negative effects on export potential.

Description: Here we report on data from an oceanographic cruise on the German research vessel Meteor covering large parts of the Mediterranean Sea during spring of 2011. The main objectives of this cruise was to conduct measurements of physical, chemical and biological variables on a section across the Mediterranean Sea with the goal of producing a synoptic picture of the distribution of relevant physical and biogeochemical properties, in order to compare those to historic data sets. During the cruise, a comprehensive data set of relevant variables following the guide lines for repeat hydrography outlined by the GO-SHIP group (http://www.go-ship.org/) was collected. The measurements include; salinity and temperature (CTD), an over-determined carbonate system, inorganic nutrients, oxygen, transient tracers (CFC-12, SF6), Helium isotopes and tritium, and carbon isotopes. The cruise sampled all major basins of the Mediterranean Sea following roughly an east-to-west section from the coast of Lebanon to through the Strait of Gibraltar, and to the coast of Portugal. Also a south-to-north section from the Ionian Sea to the Adriatic Sea was carried out. Additionally, sampling in the Aegean, Adriatic and Tyrrhenian Seas were carried out. The sections roughly followed lines and positions that have been sampled previously during other programs, thus providing the opportunity for comparative investigations of the temporal development of various parameters.

Description: Based upon thermal-infrared satellite imagery in combination with ERA-Interim atmospheric reanalysis data, we derive long-term polynya characteristics such as polynya area, thin-ice thickness distribution, and ice-production rates for a 13-year investigation period (2002–2014) for the austral winter (1 April to 30 September) in the Antarctic southern Weddell Sea. All polynya parameters are derived from daily cloud-cover corrected thin-ice thickness composites. The focus lies on coastal polynyas which are important hot spots for new-ice formation, bottom-water formation, and heat/moisture release into the atmosphere. MODIS has the capability to resolve even very narrow coastal polynyas. Its major disadvantage is the sensor limitation due to cloud cover. We make use of a newly developed and adapted spatial feature reconstruction scheme to account for cloud-covered areas. We find the sea-ice areas in front of the Ronne and Brunt ice shelves to be the most active with an annual average polynya area of 3018 ± 1298 and 3516 ± 1420 km2 as well as an accumulated volume ice production of 31 ± 13 and 31 ± 12 km3, respectively. For the remaining four regions, estimates amount to 421 ± 294 km2 and 4 ± 3 km3 (Antarctic Peninsula), 1148 ± 432 km2 and 12 ± 5 km3 (iceberg A23A), 901 ± 703 km2 and 10 ± 8 km3 (Filchner Ice Shelf), as well as 499 ± 277 km2 and 5 ± 2 km3 (Coats Land). Our findings are discussed in comparison to recent studies based on coupled sea-ice/ocean models and passive-microwave satellite imagery, each investigating different parts of the southern Weddell Sea.

Description: Localized open-ocean low-oxygen “dead zones” in the eastern tropical North Atlantic are recently discovered ocean features that can develop in dynamically isolated water masses within cyclonic eddies (CE) and anticyclonic mode-water eddies (ACME). Analysis of a comprehensive oxygen dataset obtained from gliders, moorings, research vessels and Argo floats reveals that “dead-zone” eddies are found in surprisingly high numbers and in a large area from about 4 to 22° N, from the shelf at the eastern boundary to 38° W. In total, 173 profiles with oxygen concentrations below the minimum background concentration of 40 µmol kg−1 could be associated with 27 independent eddies (10 CEs; 17 ACMEs) over a period of 10 years. Lowest oxygen concentrations in CEs are less than 10 µmol kg−1 while in ACMEs even suboxic (〈 1 µmol kg−1) levels are observed. The oxygen minimum in the eddies is located at shallow depth from 50 to 150 m with a mean depth of 80 m. Compared to the surrounding waters, the mean oxygen anomaly in the core depth range (50 and 150 m) for CEs (ACMEs) is −38 (−79) µmol kg−1. North of 12° N, the oxygen-depleted eddies carry anomalously low-salinity water of South Atlantic origin from the eastern boundary upwelling region into the open ocean. Here water mass properties and satellite eddy tracking both point to an eddy generation near the eastern boundary. In contrast, the oxygen-depleted eddies south of 12° N carry weak hydrographic anomalies in their cores and seem to be generated in the open ocean away from the boundary. In both regions a decrease in oxygen from east to west is identified supporting the en-route creation of the low-oxygen core through a combination of high productivity in the eddy surface waters and an isolation of the eddy cores with respect to lateral oxygen supply. Indeed, eddies of both types feature a cold sea surface temperature anomaly and enhanced chlorophyll concentrations in their center. The low-oxygen core depth in the eddies aligns with the depth of the shallow oxygen minimum zone of the eastern tropical North Atlantic. Averaged over the whole area an oxygen reduction of 7 µmol kg−1 in the depth range of 50 to 150 m (peak reduction is 16 µmol kg−1 at 100 m depth) can be associated with the dispersion of the eddies. Thus the locally increased oxygen consumption within the eddy cores enhances the total oxygen consumption in the open eastern tropical North Atlantic Ocean and seems to be an contributor to the formation of the shallow oxygen minimum zone.

Description: The eastern Pacific benthic foraminifer Nonionella stella Cushman & Moyer, 1930 was recorded for the first time in the Skagerrak (North Sea) and its fjords. In this short note we evaluate its migration, considering both dispersal by propagules and ship ballast tanks. We suggest that the predominantly southward surface currents along the western European seaboard and Morocco would impede a wide-range dispersal of N. stella propagules and hypothesize transportation by ship ballast tanks as the possible vector of N. stella immigration into northern European seas.

Description: Oxygen minimum zones (OMZs) that impinge on continental margins favor the release of phosphorus (P) from the sediments to the water column, enhancing primary productivity and the maintenance or expansion of low-oxygen waters. A comprehensive field program in the Peruvian OMZ was undertaken to identify the sources of benthic P at six stations, including the analysis of particles from the water column, surface sediments, and pore fluids, as well as in situ benthic flux measurements. A major fraction of solid-phase P was bound as particulate inorganic P (PIP) both in the water column and in sediments. Sedimentary PIP increased with depth in the sediment at the expense of particulate organic P (POP). The ratio of particulate organic carbon (POC) to POP exceeded the Redfield ratio both in the water column (202 ± 29) and in surface sediments (303 ± 77). However, the POC to total particulate P (TPP = POP + PIP) ratio was close to Redfield in the water column (103 ± 9) and in sediment samples (102 ± 15). This suggests that the relative burial efficiencies of POC and TPP are similar under low-oxygen conditions and that the sediments underlying the anoxic waters on the Peru margin are not depleted in P compared to Redfield. Benthic fluxes of dissolved P were extremely high (up to 1.04 ± 0.31 mmol m−2 d−1), however, showing that a lack of oxygen promotes the intensified release of dissolved P from sediments, whilst preserving the POC / TPP burial ratio. Benthic dissolved P fluxes were always higher than the TPP rain rate to the seabed, which is proposed to be caused by transient P release by bacterial mats that had stored P during previous periods when bottom waters were less reducing. At one station located at the lower rim of the OMZ, dissolved P was taken up by the sediments, indicating ongoing phosphorite formation. This is further supported by decreasing porewater phosphate concentrations with sediment depth, whereas solid-phase P concentrations were comparatively high.

Description: The temporal evolution of the physical and biogeochemical structure of an oxygen-depleted anticyclonic modewater eddy is investigated over a 2-month period using high-resolution glider and ship data. A weakly stratified eddy core (squared buoyancy frequency N2  ∼  0.1  ×  10−4 s−2) at shallow depth is identified with a horizontal extent of about 70 km and bounded by maxima in N2. The upper N2 maximum (3–5  ×  10−4 s−2) coincides with the mixed layer base and the lower N2 maximum (0.4  ×  10−4 s−2) is found at about 200 m depth in the eddy centre. The eddy core shows a constant slope in temperature/salinity (T∕S) characteristic over the 2 months, but an erosion of the core progressively narrows down the T∕S range. The eddy minimal oxygen concentrations decreased by about 5 µmol kg−1 in 2 months, confirming earlier estimates of oxygen consumption rates in these eddies. Separating the mesoscale and perturbation flow components reveals oscillating velocity finestructure ( ∼  0.1 m s−1) underneath the eddy and at its flanks. The velocity finestructure is organized in layers that align with layers in properties (salinity, temperature) but mostly cross through surfaces of constant density. The largest magnitude in velocity finestructure is seen between the surface and 140 m just outside the maximum mesoscale flow but also in a layer underneath the eddy centre, between 250 and 450 m. For both regions a cyclonic rotation of the velocity finestructure with depth suggests the vertical propagation of near-inertial wave (NIW) energy. Modification of the planetary vorticity by anticyclonic (eddy core) and cyclonic (eddy periphery) relative vorticity is most likely impacting the NIW energy propagation. Below the low oxygen core salt-finger type double diffusive layers are found that align with the velocity finestructure. Apparent oxygen utilization (AOU) versus dissolved inorganic nitrate (NO3−) ratios are about twice as high (16) in the eddy core compared to surrounding waters (8.1). A large NO3− deficit of 4 to 6 µmol kg−1 is determined, rendering denitrification an unlikely explanation. Here it is hypothesized that the differences in local recycling of nitrogen and oxygen, as a result of the eddy dynamics, cause the shift in the AOU : NO3− ratio. High NO3− and low oxygen waters are eroded by mixing from the eddy core and entrain into the mixed layer. The nitrogen is reintroduced into the core by gravitational settling of particulate matter out of the euphotic zone. The low oxygen water equilibrates in the mixed layer by air–sea gas exchange and does not participate in the gravitational sinking. Finally we propose a mesoscale–submesoscale interaction concept where wind energy, mediated via NIWs, drives nutrient supply to the euphotic zone and drives extraordinary blooms in anticyclonic mode-water eddies.

Description: To understand solar cycle signals on the Earth's surface and identify the physical mechanisms responsible, surface temperature variations from observations as well as climate model data are analyzed to characterize their spatial structure. The solar signal in the annual mean surface temperature is characterized by (i) mid-latitude warming and (ii) no warming in the tropics. The mid-latitude warming during solar maxima in both hemispheres is associated with a downward penetration of zonal mean zonal wind anomalies from the upper stratosphere during late winter. During Northern Hemisphere winter this is manifested in a modulation of the polar-night jet whereas in the Southern Hemisphere the subtropical jet plays the major role. Warming signals are particularly apparent over the Eurasian continent and ocean frontal zones, including a previously reported lagged response over the North Atlantic. In the tropics, local warming occurs over the Indian and central Pacific oceans during high solar activity. However, this warming is counter balanced by cooling over the cold tongue sectors in the southeastern Pacific and the South Atlantic, and results in a very weak zonally averaged tropical mean signal. The cooling in the ocean basins is associated with stronger cross-equatorial winds resulting from a northward shift of the ascending branch of the Hadley circulation during solar maxima. To understand the complex processes involved in the solar signal transfer, results of an idealized middle atmosphere–ocean coupled model experiment on the impact of stratospheric zonal wind changes are compared with solar signals in observations. The model results suggest that both tropical and extra-tropical solar surface signals can result from circulation changes in the upper stratosphere through (i) a downward migration of wave–zonal mean flow interactions and (ii) changes in the stratospheric mean meridional circulation. These experiments support earlier evidence of an indirect solar influence from the stratosphere.

Description: The stable carbon isotope composition of dissolved inorganic carbon (δ13C-DIC) can be used to quantify fluxes within the carbon system. For example, knowing the δ13C signature of the inorganic carbon pool can help in describing the amount of anthropogenic carbon in the water column. The measurements can also be used for evaluating modeled carbon fluxes, for making basin-wide estimates of anthropogenic carbon, and for studying seasonal and interannual variability or decadal trends in interior ocean biogeochemistry. For all these purposes, it is not only important to have a sufficient amount of data, but these data must also be internally consistent and of high quality. In this study, we present a δ13C-DIC dataset for the North Atlantic which has undergone secondary quality control. The data originate from oceanographic research cruises between 1981 and 2014. During a primary quality control step based on simple range tests, obviously bad data were flagged. In a second quality control step, biases between measurements from different cruises were quantified through a crossover analysis using nearby data of the respective cruises, and values of biased cruises were adjusted in the data product. The crossover analysis was possible for 24 of the 32 cruises in our dataset, and adjustments were applied to 11 cruises. The internal accuracy of this dataset is 0.017 ‰.

Description: Southeast Asian rivers convey large amounts of organic carbon, but little is known about the fate of this terrestrial material in estuaries. Although Southeast Asia is, by area, considered a hotspot of estuarine carbon dioxide (CO2) emissions, studies in this region are very scarce. We measured dissolved and particulate organic carbon, as well as CO2 partial pressures and carbon monoxide (CO) concentrations in two tropical estuaries in Sarawak, Malaysia, whose coastal area is covered by carbon-rich peatlands. We surveyed the estuaries of the rivers Lupar and Saribas during the wet and dry season, respectively. Carbon-to-nitrogen ratios suggest that dissolved organic matter (DOM) is largely of terrestrial origin. We found evidence that a large fraction of this carbon is respired. The median pCO(2) in the estuaries ranged between 640 and 5065 mu atm with little seasonal variation. CO2 fluxes were determined with a floating chamber and estimated to amount to 14-268 mol m(-2) yr(-1), which is high compared to other studies from tropical and subtropical sites. Estimates derived from a merely wind-driven turbulent diffusivity model were considerably lower, indicating that these models might be inappropriate in estuaries, where tidal currents and river discharge make an important contribution to the turbulence driving water-air gas exchange. Although an observed diurnal variability of CO concentrations suggested that CO was photochemically produced, the overall concentrations and fluxes were relatively moderate (0.4-1.3 nmol L-1 and 0.7-1.8 mmol m(-2) yr(-1)) if compared to published data for oceanic or upwelling systems. We attributed this to the large amounts of suspended matter (4-5004 mg L-1), limiting the light penetration depth and thereby inhibiting CO photoproduction. We concluded that estuaries in this region function as an efficient filter for terrestrial organic carbon and release large amounts of CO2 to the atmosphere. The Lupar and Saribas rivers deliver 0.3 +/- 0.2 TgC yr(-1) to the South China Sea as organic carbon and their mid-estuaries release approximately 0.4 +/- 0.2 TgC yr(-1) into the atmosphere as CO2.

Description: Changes in incoming solar ultraviolet radiation over the 11-year solar cycle affect stratospheric ozone abundances. It is important to quantify the magnitude, structure, and seasonality of the associated solar-ozone response (SOR) to understand the impact of the 11-year solar cycle on climate. Part 1 of this two-part study uses multiple linear regression analysis to extract the SOR in a number of recently updated satellite ozone datasets covering different periods within the epoch 1970 to 2013. The annual mean SOR in the updated version 7.0 (v7.0) Stratospheric Aerosol and Gas Experiment (SAGE) II number density dataset (1984–2004) is very consistent with that found in the previous v6.2. In contrast, we find a substantial decrease in the magnitude of the SOR in the tropical upper stratosphere in the SAGE II v7.0 mixing ratio dataset (∼ 1 %) compared to the v6.2 (∼ 4 %). This difference is shown to be largely attributable to the change in the independent stratospheric temperature dataset used to convert SAGE II ozone number densities to mixing ratios. Since these temperature records contain substantial uncertainties, we suggest that datasets based on SAGE II number densities are currently most reliable for evaluating the SOR. We further analyse three extended ozone datasets that combine SAGE II v7.0 number densities with more recent GOMOS (Global Ozone Monitoring by Occultation of Stars) or OSIRIS (Optical Spectrograph and Infrared Imager System) measurements. The extended SAGE–OSIRIS dataset (1984–2013) shows a smaller and less statistically significant SOR across much of the tropical upper stratosphere compared to the SAGE II data alone. In contrast, the two SAGE–GOMOS datasets (1984–2011) show SORs that are in closer agreement with the original SAGE II data and therefore appear to provide a more reliable estimate of the SOR. We also analyse the SOR in the recent Solar Backscatter Ultraviolet Instrument (SBUV) Merged Ozone Dataset (SBUVMOD) version 8.6 (VN8.6) (1970–2012) and SBUV Merged Cohesive VN8.6 (1978–2012) datasets and compare them to the previous SBUVMOD VN8.0 (1970–2009). Over their full lengths, the three records generally agree in terms of the broad magnitude and structure of the annual mean SOR. The main difference is that SBUVMOD VN8.6 shows a smaller and less significant SOR in the tropical upper stratosphere and therefore more closely resembles the SAGE II v7.0 mixing ratio data than does the SBUV Merged Cohesive VN8.6, which has a more continuous SOR of ∼ 2 % in this region. The sparse spatial and temporal sampling of limb satellite instruments prohibits the extraction of sub-annual variations in the SOR from SAGE-based datasets. However, the SBUVMOD VN8.6 dataset suggests substantial month-to-month variations in the SOR, particularly in the winter extratropics, which may be important for the proposed high-latitude dynamical response to the solar cycle. Overall, the results highlight substantial uncertainties in the magnitude and structure of the observed SOR from different satellite records. The implications of these uncertainties for understanding and modelling the effects of solar variability on climate should be explored.

Description: Halocarbons are produced naturally in the oceans by biological and chemical processes. They are emitted from surface seawater into the atmosphere, where they take part in numerous chemical processes such as ozone destruction and the oxidation of mercury and dimethyl sulfide. Here we present oceanic and atmospheric halocarbon data for the Peruvian upwelling zone obtained during the M91 cruise onboard the research vessel METEOR in December 2012. Surface waters during the cruise were characterized by moderate concentrations of bromoform (CHBr3) and dibromomethane (CH2Br2) correlating with diatom biomass derived from marker pigment concentrations, which suggests this phytoplankton group is a likely source. Concentrations measured for the iodinated compounds methyl iodide (CH3I) of up to 35.4 pmol L−1, chloroiodomethane (CH2ClI) of up to 58.1 pmol L−1 and diiodomethane (CH2I2) of up to 32.4 pmol L−1 in water samples were much higher than previously reported for the tropical Atlantic upwelling systems. Iodocarbons also correlated with the diatom biomass and even more significantly with dissolved organic matter (DOM) components measured in the surface water. Our results suggest a biological source of these compounds as a significant driving factor for the observed large iodocarbon concentrations. Elevated atmospheric mixing ratios of CH3I (up to 3.2 ppt), CH2ClI (up to 2.5 ppt) and CH2I2 (3.3 ppt) above the upwelling were correlated with seawater concentrations and high sea-to-air fluxes. During the first part of the cruise, the enhanced iodocarbon production in the Peruvian upwelling contributed significantly to tropospheric iodine levels, while this contribution was considerably smaller during the second part.

Description: The stable carbon isotope composition of dissolved inorganic carbon (δ13CDIC) in seawater was measured in samples collected during June—July 2014 in the subpolar North Atlantic. Sample collection was carried out on the RRS James Clark Ross cruise JR302, part of the "Radiatively Active Gases from the North Atlantic Region and Climate Change" (RAGNARoCC) research programme. The observed δ13CDIC values for cruise JR302 fall in a range from 0.07 ‰ to +1.95 ‰, relative to the Vienna Peedee Belemnite standard. From duplicate samples collected during the cruise, the 1σ precision for the 341 results is 0.08 ‰, which is similar to our previous work and other studies of this kind. We also performed a cross-over analysis using nearby historical δ13CDIC data, which indicated that there were no significant systematic offsets between our measurements and previously published results. We also included seawater reference material (RM) produced by A. G. Dickson (Scripps Institution of Oceanography, USA) in every batch of analysis, enabling us to improve upon the calibration and quality-control procedures from a previous study. The δ13CDIC is consistent within each RM batch, although its value is not certified. We report δ13CDIC values of 1.15 ± 0.03 ‰ and 1.27 ± 0.05 ‰ for batches 141 and 144 respectively. Our JR302 δ13CDIC data can be used – along with measurements of other biogeochemical variables – to constrain the processes that control DIC in the interior ocean, in particular the oceanic uptake of anthropogenic carbon dioxide and the biological carbon pump. Our δ13CDIC results are available from the British Oceanographic Data Centre – doi:10.5285/22235f1a-b7f3-687f-e053-6c86abc0c8a6.

Description: The ocean is a substantial source of nitrous oxide (N2O) to the atmosphere, but little is known about how this flux might change in the future. Here, we investigate the potential evolution of marine N2O emissions in the 21st century in response to anthropogenic climate change using the global ocean biogeochemical model NEMO-PISCES. Assuming nitrification as the dominant N2O formation pathway, we implemented two different parameterizations of N2O production which differ primarily under low-oxygen (O2) conditions. When forced with output from a climate model simulation run under the business-as-usual high-CO2 concentration scenario (RCP8.5), our simulations suggest a decrease of 4 to 12 % in N2O emissions from 2005 to 2100, i.e., a reduction from 4.03/3.71 to 3.54/3.56 TgN yr−1 depending on the parameterization. The emissions decrease strongly in the western basins of the Pacific and Atlantic oceans, while they tend to increase above the oxygen minimum zones (OMZs), i.e., in the eastern tropical Pacific and in the northern Indian Ocean. The reduction in N2O emissions is caused on the one hand by weakened nitrification as a consequence of reduced primary and export production, and on the other hand by stronger vertical stratification, which reduces the transport of N2O from the ocean interior to the ocean surface. The higher emissions over the OMZ are linked to an expansion of these zones under global warming, which leads to increased N2O production, associated primarily with denitrification. While there are many uncertainties in the relative contribution and changes in the N2O production pathways, the increasing storage seems unequivocal and determines largely the decrease in N2O emissions in the future. From the perspective of a global climate system, the averaged feedback strength associated with the projected decrease in oceanic N2O emissions amounts to around −0.009 W m−2 K−1, which is comparable to the potential increase from terrestrial N2O sources. However, the assessment for a potential balance between the terrestrial and marine feedbacks calls for an improved representation of N2O production terms in fully coupled next-generation Earth system models.

Description: In this study we look beyond the previously studied effects of oceanic CO2 injections on atmospheric and oceanic reservoirs and also account for carbon cycle and climate feedbacks between the atmosphere and the terrestrial biosphere. Considering these additional feedbacks is important since backfluxes from the terrestrial biosphere to the atmosphere in response to reducing atmospheric CO2 can further offset the targeted reduction. To quantify these dynamics we use an Earth system model of intermediate complexity to simulate direct injection of CO2 into the deep ocean as a means of emissions mitigation during a high CO2 emission scenario. In three sets of experiments with different injection depths, we simulate a 100-year injection period of a total of 70 GtC and follow global carbon cycle dynamics over another 900 years. In additional parameter perturbation runs, we varied the default terrestrial photosynthesis CO2 fertilization parameterization by ±50 % in order to test the sensitivity of this uncertain carbon cycle feedback to the targeted atmospheric carbon reduction through direct CO2 injections. Simulated seawater chemistry changes and marine carbon storage effectiveness are similar to previous studies. As expected, by the end of the injection period avoided emissions fall short of the targeted 70 GtC by 16–30 % as a result of carbon cycle feedbacks and backfluxes in both land and ocean reservoirs. The target emissions reduction in the parameter perturbation simulations is about 0.2 and 2 % more at the end of the injection period and about 9 % less to 1 % more at the end of the simulations when compared to the unperturbed injection runs. An unexpected feature is the effect of the model's internal variability of deep-water formation in the Southern Ocean, which, in some model runs, causes additional oceanic carbon uptake after injection termination relative to a control run without injection and therefore with slightly different atmospheric CO2 and climate. These results of a model that has very low internal climate variability illustrate that the attribution of carbon fluxes and accounting for injected CO2 may be very challenging in the real climate system with its much larger internal variability.

Description: Deoxygenation in the Baltic Sea endangers fish yields and favours noxious algal blooms. Yet, vertical transport processes ventilating the oxygen-deprived waters at depth and replenishing nutrient-deprived surface waters (thereby fuelling export of organic matter to depth), are not comprehensively understood. Here, we investigate the effects of the interaction between surface currents and winds (also referred to as eddy/wind effects) on upwelling in an eddy-rich general ocean circulation model of the Baltic Sea. Contrary to expectations we find that accounting for current/wind effects does inhibit the overall vertical exchange between oxygenated surface waters and oxygen-deprived water at depth. At major upwelling sites, however, as e.g. off the south coast of Sweden and Finland, the reverse holds: the interaction between topographically steered surface currents with winds blowing over the sea results in a climatological sea surface temperature cooling of 0.5 K. This implies that current/wind effects drive substantial local upwelling of cold and nutrient-replete waters.

Description: The biogeochemistry of the river-sea interface was studied in the Kem' River (the largest river flowing to the White Sea from Karelian coast) estuary and adjacent area of the White Sea onboard the RV "Ekolog" in summer 2001, 2002 and 2003. The study area can be divided into 3 zones: I - the estuary itself, with water depth from 1 to 5m and low salinity in the surface layer (salinity is lower than 0.2psu in the Kem' River and varies from 15 to 20psu in outer part of this zone); II - the intermediate zone with depths from 5 to 10m and salinity at the surface from 16 to 22psu; III - the marine zone with depths from 10 to 29 m and salinity 21-24.5psu. Highest concentrations of the suspended particulate matter (SPM) were registered in the Kem' mouth (5-7mg/l). They sharply decreased to values 〈1mg/l towards the sea. At beginning of July 2001, particulate organic carbon (POC) concentration in the river mouth was 404µg/l and POC content in total SPM was 5.64%. In the marine part of the studied area the POC concentration varied from 132 to 274µg/l and the POC contents in suspended matter increased to 19-52.6%. These studies show, that the majority of riverborne suspended matter in the Kem' estuary deposits near the river mouth within the 20psu isohaline, where sedimentation of the suspended matter takes place. The role of fresh-water phytoplankton species decreases and the role of marine species increases from the river to sea and the percentage of green algae decreases and the role of diatoms increases. The organic carbon (Corg) to nitrogen (N) ratio (Corg/N) in both suspended matter and bottom sediments decreases from the river to the marine part of the mixing zone (from 8.5 to 6.1 in the suspended matter and from 14.6 to 7.5 in the bottom sediments), demonstrating that content of terrestrial-derived organic matter decreases and content of marine organic matter increases from the river mouth to the sea. The Kem' estuary exhibits a similar character of biogeochemial processes as in the large Arctic estuaries, but the scale of these processes (amount of river input of SPM, POC, area of estuaries) is different.

Description: We use numerical climate simulations, paleoclimate data, and modern observations to study the effect of growing ice melt from Antarctica and Greenland. Meltwater tends to stabilize the ocean column, inducing amplifying feedbacks that increase subsurface ocean warming and ice shelf melting. Cold meltwater and induced dynamical effects cause ocean surface cooling in the Southern Ocean and North Atlantic, thus increasing Earth's energy imbalance and heat flux into most of the global ocean's surface. Southern Ocean surface cooling, while lower latitudes are warming, increases precipitation on the Southern Ocean, increasing ocean stratification, slowing deepwater formation, and increasing ice sheet mass loss. These feedbacks make ice sheets in contact with the ocean vulnerable to accelerating disintegration. We hypothesize that ice mass loss from the most vulnerable ice, sufficient to raise sea level several meters, is better approximated as exponential than by a more linear response. Doubling times of 10, 20 or 40 years yield multi-meter sea level rise in about 50, 100 or 200 years. Recent ice melt doubling times are near the lower end of the 10–40-year range, but the record is too short to confirm the nature of the response. The feedbacks, including subsurface ocean warming, help explain paleoclimate data and point to a dominant Southern Ocean role in controlling atmospheric CO2, which in turn exercised tight control on global temperature and sea level. The millennial (500–2000-year) timescale of deep-ocean ventilation affects the timescale for natural CO2 change and thus the timescale for paleo-global climate, ice sheet, and sea level changes, but this paleo-millennial timescale should not be misinterpreted as the timescale for ice sheet response to a rapid, large, human-made climate forcing. These climate feedbacks aid interpretation of events late in the prior interglacial, when sea level rose to +6–9 m with evidence of extreme storms while Earth was less than 1 °C warmer than today. Ice melt cooling of the North Atlantic and Southern oceans increases atmospheric temperature gradients, eddy kinetic energy and baroclinicity, thus driving more powerful storms. The modeling, paleoclimate evidence, and ongoing observations together imply that 2 °C global warming above the preindustrial level could be dangerous. Continued high fossil fuel emissions this century are predicted to yield (1) cooling of the Southern Ocean, especially in the Western Hemisphere; (2) slowing of the Southern Ocean overturning circulation, warming of the ice shelves, and growing ice sheet mass loss; (3) slowdown and eventual shutdown of the Atlantic overturning circulation with cooling of the North Atlantic region; (4) increasingly powerful storms; and (5) nonlinearly growing sea level rise, reaching several meters over a timescale of 50–150 years. These predictions, especially the cooling in the Southern Ocean and North Atlantic with markedly reduced warming or even cooling in Europe, differ fundamentally from existing climate change assessments. We discuss observations and modeling studies needed to refute or clarify these assertions.

Description: Coastal zones are important source regions for a variety of trace gases including halocarbons and sulphur-bearing species. While salt-marshes, macroalgae and phytoplankton communities have been intensively studied, little is known about trace gas fluxes in seagrass meadows. Here we report results of a newly developed dynamic flux chamber system that can be deployed in intertidal areas over full tidal cycles allowing for high time resolved measurements. The trace gases measured in this study included carbon dioxide (CO2), methane (CH4) and a variety of hydrocarbons, halocarbons and sulphur-bearing compounds. The high time resolved CO2 and CH4 flux measurements revealed a complex dynamic mediated by tide and light. In contrast to most previous studies our data indicate significantly enhanced fluxes during tidal immersion relative to periods of air exposure. Short emission peaks occured with onset of the feeder current at the sampling site. We suggest an overall strong effect of advective transport processes to explain the elevated fluxes during tidal immersion. Many emission estimates from tidally influenced coastal areas still rely on measurements carried out during low tide only. Hence, our results may have significant implications for budgeting trace gases in coastal areas. This dynamic flux chamber system provides intensive time series data of community respiration (at night) and net community production (during the day) of shallow coastal systems.