ALBERT

Description: An investigation was carried out on larvae of the oceanic tropical squid Sthenoteuthis pteropus in the equatorial Atlantic (2°30′N–7°S;12°W –8°30′E) The age of the larvae was calculated from the statolith microstructure of 20 larvae; mortality was estimated from the size structure of 1128 larvae. The larval stage lasts 32–38 days. At ages ranging from 14 to 38 days. the daily relative growth rates of mantle length decrease from 7.5 to 2.8% day −1 and from 14–16 to 5.8% of body weight day −1 At age 12–24 days, mortality rates were estimated using both raw catch data and corrected data accounting for net avoidance. The mean value of raw mortality rates was 0.189, the corrected value was 0.158. During the proboscis division (transformation of the larva into juvenile) at age 25–35 days, a sharp decrease in larval growth rates and a simultaneous increase in mortality rates (raw 0.443, corrected 0.379) were observed.

Description: A number of reproductive indices were compared with a subjective maturity scale for assessment of Loligo forbesi maturity. The ratio between nidamental gland length and mantle length corresponded well with female maturation, as did the ovary mass-soma mass and nidamental gland mass-soma mass ratios. For males, the ratio between spermatophoric complex mass and somatic mass was found to be the most suitable for maturity assessment. The timing of recruitment and maturation of L. forbesi in Irish waters was described from the size and maturity of squid in commercial landings in the south of Ireland during the years 1991–1993. Immature squid first appeared in commercial catches in July and August, and this represented the main period of recruitment. A second period of recruitment was apparent in December 1991, but was not identified in the 1992–1993 season. Mature females were present in the commercially exploited population between November and April, with a small number also found in the summer. The abundance of egg masses was used to indicate timing of spawning. Egg masses recovered from the Cork coast indicated that peak spawning occurred during the winter months, but continued on a small scale for much of the year.

Description: The transition from benthos to plankton requires multiple adaptations, yet so far it remains unclear how these are acquired in the course of the transition. To investigate this process, we analyzed the genetic diversity and distribution patterns of a group of foraminifera of the genus Bolivina with a tychopelagic mode of life (same species occurring both in benthos and plankton). We assembled a global sequence data set for this group from single-cell DNA extractions and occurrences in metabarcodes from pelagic environmental samples. The pelagic sequences all cluster within a single monophyletic clade within Bolivina. This clade harbors three distinct genetic lineages, which are associated with incipient morphological differentiation. All lineages occur in the plankton and benthos, but only one lineage exhibits no limit to offshore dispersal and has been shown to grow in the plankton. These observations indicate that the emergence of buoyancy regulation within the clade preceded the evolution of pelagic feeding and that the evolution of both traits was not channeled into a full transition into the plankton. We infer that in foraminifera, colonization of the planktonic niche may occur by sequential cooptation of independently acquired traits, with holoplanktonic species being recruited from tychopelagic ancestors

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: The sea surface microlayer (SML) is the uppermost layer of the water column that links the ocean and atmosphere. It accumulates a variety of biogenic surface-active and buoyant substances, including gelatinous material, such as transparent exopolymer particles (TEP) and Coomassie stainable particles (CSP), potentially affecting air–sea exchange processes. Here, we studied the influence of the annual cycle of phytoplankton production on organic matter (OM) accumulation in the SML relative to the subsurface water (SSW). Sampling was performed monthly from April 2012 to November 2013 at the Boknis Eck Time Series Station (Baltic Sea). For SML sampling, we used the Garrett screen, while SSW samples were collected by Niskin bottles at 1 m depth. Samples were analyzed for carbohydrates, amino acids, TEP, CSP, chlorophyll a (SSW only) and bacterial abundance. Our data showed that the SML reflected the SSW during most parts of the year, with changes mainly responding to bloom formation and decay. OM composition during phytoplankton blooms clearly differed from periods of higher bacterial abundance. Of all components investigated, only the enrichment of total carbohydrates in the SML was inversely related to the wind speed indicating that wind-driven mixing also affected the accumulation of OM in the SML during our study.

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: On 25 December 2016, a Mw 7.6 earthquake broke a portion of the Southern Chilean subduction zone south of Chiloé Island, located in the central part of the Mw 9.5 1960 Valdivia earthquake. This region is characterized by repeated earthquakes in 1960 and historical times with very sparse interseismic activity due to the subduction of a young (~15 Ma), and therefore hot, oceanic plate. We estimate the co-seismic slip distribution based on a kinematic finite fault source model, and through joint inversion of teleseismic body waves and strong motion data. The coseismic slip model yields a total seismic moment of 3.94×1020 Nm that occurred over ~30 s, with the rupture propagating mainly downdip, reaching a peak-slip of ~4.2 m. Regional moment tensor inversion of stronger aftershocks reveals thrust type faulting at depths of the plate interface. The fore- and aftershock seismicity is mostly related to the subduction interface with sparse seismicity in the overriding crust. The 2016 Chiloé event broke a region with increased locking and most likely broke an asperity of the 1960 earthquake. The updip limit of the main event, aftershocks, foreshocks and interseismic activity are spatially similar, located ~15 km offshore and parallel to Chiloé Islands west coast. The coseismic slip model of the 2016 Chiloé earthquake suggests a peak slip of 4.2 m that locally exceeds the 3.38 m slip deficit that has accumulated since 1960. Therefore, the 2016 Chiloé earthquake possibly released strain that has built up prior to the 1960 Valdivia earthquake.

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: We used a molecular phylogenetic approach to investigate species delimitation and diversification in the northeastern Atlantic and Mediterranean musseldrills of the Ocinebrina aciculata complex, based on molecular data from topotypical material of many of the nominal taxa. The complex is shown to consist of at least five species: Ocinebrina aciculata (Lamarck, 1822) from the Atlantic and western Mediterranean; O. cf. corallina (Scacchi, 1836) from the central Mediterranean Sea; O. reinai Bonomolo & Crocetta, 2012 from the Tyrrhenian Sea; O. corallinoides Pallary, 1912 from the Gulf of Gabès; and O. aegeensis n. sp. currently known from the Aegean Sea only. The new species is differentiated from the other taxa by very subtle morphological diagnostic features, although it is clearly identified by genetic distance and apomorphic DNA-sequence characters. The identity of Murex corallinus Scacchi, 1836 (type species of Ocinebrina Jousseaume, 1880) could not be defined with certainty, pending genetic comparison of specimens of the â €? large Tyrrhenian morphotype' (corresponding to the neotype, but not assayed herein) with the assayed â €? small morphotype'.

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: In 1988 and 1989 the occurrence of external diseases in fish from the German part of the Wadden Sea was quantified. Thirty-two stations, located along seven transects including four estuaries, were sampled monthly or at 3-month intervals with shrimp trawlers: more than 124 000 fish were studied. Six percent of the seven dominant fish species of 〉 12 cm length were diseased. In fish species which complete their life cycle in the Wadden Sea (gobies, hooknose, eelpout, sea scorpion) the total disease incidence was below 0.4%. In whiting, plaice and sole it was between 0.5 and 2.5%, while 6–8% of smelt, cod, eel, dab and flounder suffered from external lesions. The incidence of most diseases increased with increasing fish length. Similar geographical patterns in prevalence were observed for: (1) two types of skeletal deformities, (2) lymphocystis in dab and flounder, (3) papillomatosis in smelt and dab and (4) several infectious ulcerative diseases. Most of these ulcerative diseases of cod and flounder occurred on central estuarine stations, suggesting an impact of the relatively low, but greatly variable salinity on disease development. Buccal granulomatosis of smelt and papillomatosis of eels showed the highest incidence in the Elbe estuary as the most heavily polluted region of the Wadden Sea. A causal relationship between disease development and pollution, however, is not yet clear and requires experimental evidence. The disease types found are ranked into three priority groups regarding their indicator value for a pollution monitoring on the basis of fish diseases. There is evidence that the total disease load of fish in the Wadden Sea is higher than in other shallow coastal regions outside the North Sea. Little is known of the effects of these diseases on single fish and fish populations. “Yellow pest” of cod, which was recorded for the first time during this survey, and which occurred with an incidence of up to 14% in single estuarine samples, causes the most serious lesions and is supposed to be lethal. In total, however, the effects of a number of pathogenic parasites on fish survival is considered to be more serious than that of deformities and infectious diseases.

Description: Marine sponges are early-branched metazoans known to harbor dense and diverse microbial communities. Yet the role of the so far uncultivable alphaproteobacterial lineages that populate these sessile invertebrates remains unclear. We applied a sequence composition-dependent binning approach to assemble one Rhodospirillaceae genome from the Spongia officinalis microbial metagenome and contrast its functional features with those of closely related sponge-associated and free-living genomes. Both symbiotic and free-living Rhodospirillaceae shared a suite of common features, possessing versatile carbon, nitrogen, sulfur and phosphorus metabolisms. Symbiotic genomes could be distinguished from their free-living counterparts by the lack of chemotaxis and motility traits, enrichment of genes required for the uptake and utilization of organic sulfur compounds—particularly taurine—, higher diversity and abundance of ABC transporters, and a distinct repertoire of genes involved in natural product biosynthesis, plasmid stability, cell detoxification and oxidative stress remediation. These sessile symbionts may more effectively contribute to host fitness via nutrient exchange, and also host detoxification and chemical defense. Considering the worldwide occurrence and high diversity of sponge-associated Rhodospirillaceae verified here using a tailored in silico approach, we suggest that these organisms are not only relevant to holobiont homeostasis but also to nutrient cycling in benthic ecosystems.

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: Earlier work by the second author on the growth of oceanic squids, based on sizefrequency distributions of beaks sampled from sperm whale stomachs and on structural marks on those beaks, showed that these squids apparently had growth rates far in excess of those reported for the fastest-growing fishes, e.g. bluefin tuna. The application of recently developed methods for analysis of length-frequency distributions to some of these earlier data, and new approaches for assessing and comparing the growth performance offish and aquatic invertebrates, suggest the need for a downward revision of these high growth estimates. This is illustrated here with data on Kondakovia longimana (Cephalopoda, Onychoteuthidae) sampled off Durban and Donkergat, South Africa, in the early to mid-1960s.

Description: The antarctic krill Euphausia superba forms abundant, well-organized schools in the waters off the Antarctic Peninsula. Mean avian density is 2.6 times greater in waters where krill schools are present than in waters without krill schools. Seabird density is a good predictor of the presence of krill. Seabird density did not correlate with krill density or krill school depth. Disoriented krill routinely were observed swimming near the surface above submerged schools, providing potential prey for surface-feeding birds. Responses of seabird species to the distribution of krill schools varied. The small to medium-size procellariiform species were the best indicators of krill schools; large procellariiforms and coastal species were poor indicators. Pygoscelis penguins occurred at high densities only in the presence of krill schools. These responses are consistent with the constraints imposed by the metabolic requirements and reproductive strategies of each of these groups. Krill schools were detected near the sea surface throughout the day. Correlations between seabird density and the presence of krill during daylight hours suggest that diurnal foraging is important to the seabirds of this region.

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: Artificial selection processes associated with harvesting may operate over relatively short time scales in short-lived semelparous species. The ommastrephid squid Illex argentinus on the Patagonian Shelf is the target of a major fishery. Recent work has provided new information on the biological characteristics of squid in this fishery. That information has been utilized in the development of a model of the cohort dynamics and some of the within-seasons selective effects considered. The model results are consistent with earlier data from the shelf fisheries, supporting suggestions that the males nature and migrate earlier towards the spawning grounds than the females. Earlier maturation results in a smaller mean size in the spawning stock, while later maturation results in greater exposure to the fishery and a reduced numbers of individuals surviving to spawn. Under the current fishing regime greater egg production and a larger spermatophoric complex mass for the whole cohort is achieved by relatively late maturation. In general, however, the earlier maturation occurs, the earlier is the peak in total egg production. The within-season pattern of effort expenditure in the fishery can affect not only yield from the fishery but also the reproductive potential of the spawning stock. The management policy adopted for this fishery is likely to be conservative in terms of maintenance of a spawning stock, however, the potential for selective effects is larger and this may affect both yield and reproductive potential. The results are discussed in relation to short- and long-term effects in the fishery and the implications for future research requirements.

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: Statolith morphology and microstructure were studied in two common species of panktonic cranchiid squids, Belonella borealis [four juveniles with mantle length (ML) 375–450 mm] and Galiteuthis phyllura (13 paralarvae and juveniles, ML 9–235mm), caught near the bottom and in pelagic layers over the continental slope of Siberia in the northwest Bering Sea. The total number of growth increments within the statoliths ranged from 277 to 294 in B.borealis and from 10 to 209 in G.phyllura. Assuming that these increments were produced daily, both species grow rapidly in length (daily growth rate = 1.13mm day−1 during the first 8–10 months of their juvenile phase in the mesopelagic layers, prior to migration into deeper waters for maturation.

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: Continental hyperextension during magma-poor rifting at the Deep Galicia Margin is characterised by a complex pattern of faulting, thin continental fault blocks, and the serpentinisation, with local exhumation, of mantle peridotites along the S-reflector, interpreted as a detachment surface. In order to understand fully the evolution of these features, it is important to image seismically the structure and to model the velocity structure to the greatest resolution possible. Travel-time tomography models have revealed the long-wavelength velocity structure of this hyperextended domain, but are often insufficient to match accurately the short-wavelength structure observed in reflection seismic imaging. Here we demonstrate the application of two-dimensional (2D) time-domain acoustic full-waveform inversion to deep water seismic data collected at the Deep Galicia Margin, in order to attain a high resolution velocity model of continental hyperextension. We have used several quality assurance procedures to assess the velocity model, including comparison of the observed and modelled waveforms, checkerboard tests, testing of parameter and inversion strategy, and comparison with the migrated reflection image. Our final model exhibits an increase in the resolution of subsurface velocities, with particular improvement observed in the westernmost continental fault blocks, with a clear rotation of the velocity field to match steeply dipping reflectors. Across the S-reflector there is a sharpening in the velocity contrast, with lower velocities beneath S indicative of preferential mantle serpentinisation. This study supports the hypothesis that normal faulting acts to hydrate the upper mantle peridotite, observed as a systematic decrease in seismic velocities, consistent with increased serpentinisation. Our results confirm the feasibility of applying the full-waveform inversion method to sparse, deep water crustal datasets.

Description: Processes linked with the genesis, evolution and emplacement of silicic complexes in arcs are still poorly constrained. Of particular interest are the depth of magma production, the relative contribution of crystal fractionation versus crustal partial melting and the timescales involved. The Soufrière Volcanic Complex (SVC) on St Lucia is one of the largest silicic centres in the Lesser Antilles arc. Here we present the results of a detailed mineralogical study, including in situ Sr isotopes in plagioclase and in situ δ18O in dated zircons, of both SVC and Pre-SVC volcanic rocks to place constraints on the processes intrinsic to the development and evolution of the silicic complex. These data suggest that the production of silicic magma in the SVC occurs in two stages. The first stage involves differentiation of mafic magma by crustal assimilation and mineral fractionation in the middle–lower crust of the arc to produce magmas with intermediate compositions. These intermediate magmas are water-rich (∼7 wt %) and have high 87Sr/86Sr, Ba, Sr and La/Sm (∼5) compared with Pre-SVC lavas. Near-constant trace element and isotopic compositions throughout the SVC lifespan indicate that the same process was persistent over the last 600 kyr. In the second stage, the intermediate magmas are transferred to a shallower and more differentiated chamber (∼6 km depth). During ascent, any crystals or xenocrysts residual from stage one in the deeper chamber become fully resorbed and the magma crystallizes calcic amphibole microphenocrysts, followed by anorthite-rich plagioclase close to or at the water saturation depth. During mixing upon recharge within the shallow chamber, anorthite-rich plagioclase from the recharging magma is partially resorbed; so are the crystals in equilibrium with the resident differentiated magma. The recharge event probably causes chamber-wide convection. Mixing is thought to trigger eruption of the silicic complex magmas.

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: Eukaryotic phytoplankton exhibit an enormous species richness, displaying a range of phylogenetic, morphological and physiological diversity. Yet, until recently, very little was known about the diversity, genetic variation and evolutionary processes within species and populations. An approach to explore this diversity and to understand evolution of phytoplankton is to use population genetics as a conceptual framework and methodology. Here, we discuss the patterns, processes and questions that population genetic studies have revealed in eukaryotic phytoplankton. First, we describe the main biological processes generating genetic variation. We specifically discuss the importance of life-cycle complexity for genetic and phenotypic diversity and consider how such diversity can be maintained during blooms when rapid asexual proliferation dominates. Next, we explore how genetic diversity is partitioned over time and space, with a focus on the processes shaping this structure, in particular selection and genetic exchange. Our aim is also to show how population genetics can be used to make inferences about realized dispersal and sexual recombination, as these processes are so difficult to study directly. Finally, we highlight important open questions and suggest promising avenues for future studies that will be made possible by new sequencing technologies

Description: Marine spatial planning (MSP) is considered a valuable tool in the ecosystem-based management of marine areas. Predictive modelling may be applied in the MSP framework to obtain spatially explicit information about biodiversity patterns. The growing number of statistical approaches used for this purpose implies the urgent need for comparisons between different predictive techniques. In this study, we evaluated the performance of selected machine learning and regression-based methods that were applied for modelling fish community indices. We hypothesized that habitat features can influence fish assemblage and investigated the effect of environmental gradients on demersal fish diversity (species richness and Shannon–Weaver Index). We used fish data from the Baltic International Trawl Surveys (2001–2014) and maps of six potential predictors: bottom salinity, depth, seabed slope, growth season bottom temperature, seabed sediments and annual mean bottom current velocity. We compared the performance of six alternative modelling approaches: generalized linear models, generalized additive models, multivariate adaptive regression splines, support vector machines, boosted regression trees and random forests. We applied repeated 10-fold cross-validation, using accuracy as the measure of model quality. Finally, we selected random forest as the best performing algorithm and implemented it for the spatial prediction of fish diversity from the Baltic Proper to the Kattegat. To obtain information on the data reliability and confidence of the developed models, which are essential for MSP, we estimated the uncertainty of predictions with standard deviation of predictions obtained from all the trees in the ensemble random forest method. We showed how state-of-the-art predictive techniques, based on easily available data and simple Geographic Information System tools, can be used to obtain reliable spatial information about fish diversity. Our comparative work highlighted the potential of machine learning method to reduce prediction error in modelling of demersal fish diversity in the framework of MSP.

Description: The Mozambique Ridge, a prominent basement high in the southwestern Indian Ocean, consists of four major geomorphological segments associated with numerous phases of volcanic activity in the Lower Cretaceous. The nature and origin of the Mozambique Ridge have been intensely debated with one hypothesis suggesting a Large Igneous Province origin. High-resolution seismic reflection data reveal a large number of extrusion centres with a random distribution throughout the southern Mozambique Ridge and the nearby Transkei Rise. Intrabasement reflections emerge from the extrusion centres and are interpreted to represent massive lava flow sequences. Such lava flow sequences are characteristic of eruptions leading to the formation of continental and oceanic flood basalt provinces, hence supporting a Large Igneous Province origin of the Mozambique Ridge. We observe evidence for widespread post-sedimentary magmatic activity that we correlate with a southward propagation of the East African Rift System. Based on our volumetric analysis of the southern Mozambique Ridge we infer a rapid sequential emplacement between ∼131 and ∼125 Ma, which is similar to the short formation periods of other Large Igneous Provinces like the Agulhas Plateau.

Description: Gray whales (Eschrichtius robustus) feeding in the northern Bering Sea produce prey-rich mud plumes that provide ephemeral foraging opportunities for seabirds. Approximately 67% of all gray whales were attended by birds. In four whale-associating bird species (Northern Fulmar, Fulmarus glacialis; Red Phalarope, Phalaropus fulicaria; Black-legged Kitti-wake, Rissa tridactyla; and Thick-billed Murre, Uria lomvia), from 17 to 87% of all individuals that we observed on the water or foraging were in the whales' mud plumes. The combined density of these same four species was strongly correlated with whale density over a broad range of spatial scales. The whale-associating seabirds exhibited species-specific patterns of foraging behavior at plumes, including differences in mean group size, mean residence time, and patterns of movement between plumes. Birds tended to form larger groups and to form more mixed-species flocks in association with whales. The association of marine birds with gray whales in the Bering Sea provides a model system for examining seabird interactions at fine-scale oceanographic patches and demonstrates the importance of these patches in shaping patterns of seabird distribution and behavior.

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 Nifonea submarine volcano rises 1000 m above the seafloor of the Vate Trough back-arc basin behind the New Hebrides island arc. This large volcanic edifice has a caldera of ∼8 km diameter and is connected to two ∼20 km long volcanic rift zones in the back-arc basin. We present new chemical and isotope data for volcanic glasses and whole-rocks from both the volcano and its rift zones. Lavas from Nifonea volcano show an evolution towards more incompatible element enrichment, with the most enriched lavas being the youngest eruption products on the caldera floor. These are products of significant fractional crystallization, show minor contamination by hydrothermal fluids (〈0·3%) and reflect mixing of melts derived from depleted upper mantle and melts from an enriched source similar to those occurring in the North Fiji Basin. The enrichment in Nb of these lavas is comparable with that of some lavas from the New Hebrides island arc (e.g. Mota Lava island), where these coexist with typical island arc basalts. The lavas erupted along the rift zones in the Vate Trough back-arc basin are relatively depleted in incompatible elements, indicating melting of depleted upper mantle with a minor addition of a sediment-derived fluid. Our observations suggest that the mantle beneath Vate Trough is heterogeneous on a small scale (〈20 km) and that the occurrence of these enriched and fertile mantle portions has a stronger control on melting processes than the influx from the subducting slab, as all samples were recovered at a similar distance from the trench.

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 influence of winter on the selection of dominant taxa for the phytoplankton spring bloom was studied in batch culture experiments. Different natural phytoplankton assemblages from different phases of the temperate zone winter were exposed to varying periods of darkness (0, 6/7, 13 and 19 weeks) followed by a re-exposure to saturating light intensity for 14 days to experimentally simulate the onset of spring. The results showed that dark incubation has a strong effect on shaping the phytoplankton community composition. Many taxa disappeared in the absolute darkness. Dark survival ability might be an important contributing factor for the success of diatoms in spring. Different phytoplankton starting assemblages were dominated by the same bloom-forming diatoms, Skeletonema marinoi and Thalassosira spp., after dark incubation for only 6 weeks, irrespective of the high dissimilarities between phytoplankton communities. The growth capacity of surviving phytoplankton is almost unimpaired by darkness. Similar growth rates as that before darkness could be resumed for the surviving taxa with a potential lag time of 1–7 days dependent on taxon and the duration of darkness.

Description: The general positive effect of warmer winters on the abundance of small-sized zooplankton in the following spring and early summer has been reported from different parts of the Baltic Sea, but the mechanism of this link is not clear. Although causal links cannot be deduced with confidence from observational data, sufficiently detailed analyses can nevertheless provide insights to the potential mechanisms. We present an example of such an analysis, scrutinizing the effects of winter and spring hydroclimate on the abundance of small-sized dominant calanoid copepods (Eurytemora affinis and Acartia spp.), using data from 2080 zooplankton samples collected over 55 years (1957–2012) from a shallow coastal habitat (Pärnu Bay, Gulf of Riga) in the Baltic Sea. Our results indicated that the milder winters brought about higher abundances, and reduced seasonality of small-sized copepods, whereas ambient sea surface temperature (SST) mostly affected the relative abundance of adult stages. The sliding window correlation tests revealed temporal shifts in the effects of controlling variables: with the continuous increase in SST, the effect of winter temperature on the abundance of Acartia spp. weakened. In contrast, E. affinis was consistently affected by SST, but the effect of winter temperature was more pronounced during the period of on average colder winters.

Description: Overfishing and rapid environmental shifts pose severe challenges to the resilience and viability of marine fish populations. To develop and implement measures that enhance species’ adaptive potential to cope with those pressures while, at the same time, ensuring sustainable exploitation rates is part of the central goal of fisheries management. Here, we argue that a combination of biophysical modelling and population genomic assessments offer ideal management tools to define stocks, their physical connectivity and ultimately, their short-term adaptive potential. To date, biophysical modelling has often been confined to fisheries ecology whereas evolutionary hypotheses remain rarely considered. When identified, connectivity patterns are seldom explored to understand the evolution and distribution of adaptive genetic variation, a proxy for species’ evolutionary potential. Here, we describe a framework that expands on the conventional seascape genetics approach by using biophysical modelling and population genomics. The goals are to identify connectivity patterns and selective pressures, as well as putative adaptive variants directly responding to the selective pressures and, ultimately, link both to define testable hypotheses over species response to shifting ecological conditions and overexploitation.

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: Mortality of two dogs, after exposure to water of a brackish lake on the German North Sea coast in 1990, is considered to be caused by a toxic Nodularia spumigena Mertens bloom.

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: This study presents a new method (LBB) for the analysis of length frequency data from commercial catches. LBB works for species that grow throughout their lives, such as most commercially-important fish and invertebrates, and requires no input in addition to length frequency data. It estimates asymptotic length, length at first capture, relative natural mortality, and relative fishing mortality. Standard fisheries equations can then be used to approximate current exploited biomass relative to unexploited biomass. In addition, these parameters allow the estimation of length at first capture that would maximize catch and biomass for a given fishing effort, and estimation of a proxy for the relative biomass capable of producing maximum sustainable yields. Relative biomass estimates of LBB were not significantly different from the “true” values in simulated data and were similar to independent estimates from full stock assessments. LBB also presents a new indicator for assessing whether an observed size structure is indicative of a healthy stock. LBB results will obviously be misleading if the length frequency data do not represent the size composition of the exploited size range of the stock or if length frequencies resulting from the interplay of growth and mortality are masked by strong recruitment pulses.

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: The Rhynchoteuthion larvae of Illex argentinus resulting from summer spawning in North Patagonic shelf waters, its distribution and abundance, are described in this paper. The material was collected in the Argentine Sea (35–55°S) by means of plankton nets. The research cruises were made by the R/V Shinkai Maru and the R/V Walther Herwig during the period April 1978 to April 1979. The most important spawning ground of the summer spawning subpopulation is found in continental-shelf waters (between 43 and 46°S) during the period December-February. This area was established on the basis of both ripe (December) and spent females (February). The larvae which were caught during the same period, especially in March, confirmed the spawning area of this demographic unit. The larvae showed the length of the mantle (ML) to be from 1.2 to 6.5 mm. Tentacles were splitting in specimens from 5.0 to 6.5 mm ML (transition stage). When 7.0 mm ML or more, all specimens were juveniles and had their tentacles completely separated. Larvae were characterized as type C, following the proposal of Sato (1973) and Sato and Sawada (1974) in the Bulletin of the Shizuoka Prefectural Fisheries Experimental Station.

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: A joint research cruise (Japan-Argentina-Uruguay) was carried out in the South-western Atlantic during August-September 1989 in order to study the winter-spawning and hatchery areas of Illex argentinus, and also the migration pattern of the juveniles towards the continental shelf. A few Rhynchoteuthion larvae were found in subtropical waters of the Brazil Current, next to the Brazil-Malvinas confluence, and in the frontal zone with shelf water, but never at temperatures below 14°C. Large numbers of juveniles found in subantarctic waters (6–10°C) on the shelf were probably migrating southward from their hatchery grounds following the zooplankton concentrations on which they were feeding.