ALBERT

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

Description: An abrupt lithofacies change between calcareous shale and noncalcareous shale occurs in strata deposited in the mid-Cretaceous Greenhorn Seaway in the southeastern corner of Montana. The facies were correlated lithostratigraphically using bentonites and calcarenites. The lithocorrelations were then refined using ammonites, foraminifera, and calcareous nannofossils. Twenty-five time slices were defined within the upper middle and lower upper Cenomanian strata. Biofacies analysis indicate that lithofacies changes record the boundary or oceanic front between two water masses with distinctly different paleoceanographic conditions. One water mass entered the seaway from the Arctic and the other from the Gulf of Mexico/Tethys. The microfauna and microflora permit interpretation of the environmental conditions in each water mass. At times when the front was near vertical, the two water masses were of the same density but of different temperatures and salinities.

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

Description: Alteration patterns in the lavas and dykes of the Troodos Ophiolite, Cyprus, record a complex history of axial hydrothermal alteration, crustal aging, and subsequent uplift and emplacement of the ophiolite. Field mapping shows that distribution of five alteration zones, each with distinct mineralogical, geochemical, and hydrologie characteristics, is influenced by igneous stratigraphy, structure, and the nature and thickness of the overlying sediments. Paragenetic sequences of secondary minerals indicate that alteration conditions changed progressively as the crust cooled and moved off-axis. Along spreading axes, low temperatures (≤50°C) were maintained by the rapid flow of seawater in and out of the lavas, and only minimal alteration took place. In contrast, lower water/rock ratios and higher temperatures (〉200°C) in the dykes promoted extensive seawater-rock interaction. Although the sharp rise in temperature between the two regimes generally coincides with the lava-dyke transition, late-stage intrusions or hydrothermal upwelling zones locally cause high-temperature alteration to extend upward into the lavas. As a segment of crust moved off-axis, temperatures remained low in the lavas and progressively decreased, from 〉250° to 〈80°C, in the dykes. High permeability in the uppermost lavas led to the downward migration of an oxidative alteration front whose thickness and spatial distribution was dependent upon the rate and nature of sedimentation and, thus, the original seafloor morphology. Although field relations show that alteration has a consistent vertical pattern in Troodos, the alteration zones are not laterally continuous, and the stratigraphie depth of their boundaries varies considerably.

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

Description: The morphology and structure of the submarine flanks of the Canary Islands were mapped using the GLORIA long-range side-scan sonar system, bathymetric multibeam systems, and sediment echosounders. Twelve young (〈2 Ma) giant landslides have been identified on the submarine flanks of the Canary Islands up to now. Older landslide events are long buried under a thick sediment cover due to high sedimentation rates around the Canary Islands. Most slides were found on the flanks of the youngest and most active islands of La Palma, El Hierro, and Tenerife, but young giant landslides were also identified on the flanks of the older (15–20 Ma) but still active eastern islands. Large-scale mass wasting is an important process during all periods of major magmatic activity. The long-lived volcanic constructive history of the islands of the Canary Archipelago is balanced by a correspondingly long history of destruction, resulting in a higher landslide frequency for the Canary Islands compared to the Hawaiian Islands, where giant landslides only occur late in the period of active shield growth. The lower stability of the flanks of the Canaries is probably due to the much steeper slopes of the islands, a result of the abundance of highly evolved intrusive and extrusive rocks. Another reason for the enhanced slope instability is the abundance of pyroclastic deposits on Canary Islands resulting from frequent explosive eruptions due to the elevated volatile contents in the highly alkalic magmas. Dike-induced rifting is most likely the main trigger mechanism for destabilization of the flanks. Flank collapses are a major geological hazard for the Canary Islands due to the sector collapses themselves as well as triggering of tsunamis. In at least one case, a giant lateral blast occurred when an active magmatic or hydrothermal system became unroofed during flank collapse.

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

Description: In the western equatorial Pacific, the El Niño/Southern Oscillation (ENSO) phenomenon is characterized by precipitation variability associated with the migration of the Indonesian low pressure cell to the region of the date line and the equator. During ENSO events, Tarawa Atoll (1°N, 172°E) experiences heavy rainfall which has an estimated δ18O of about −8 to −10‰ δ18OSMOW. At Tarawa, sufficient precipitation of this composition falls during ENSO events to alter the δ18O and the salinity of the surface waters. Oxygen isotope records from two corals collected off the reef crest of Tarawa reflect rainfall variations associated with both weak and strong ENSO conditions, with approximately monthly resolution. Coral skeletal δ18O variations due to small sea surface temperature (SST) changes are secondary. These records demonstrate the remarkable ability of this technique to reconstruct variations in the position of the Indonesian Low from coral δ18O records in the western equatorial Pacific, a region which has few paleoclimatic records. The coral isotopic data correctly resolve the relative magnitudes of recent variations in the Southern Oscillation Index. Combining the Tarawa record with an oxygen isotopic history from a Galápagos Islands coral demonstrates the ability to distinguish the meteorologic (precipitation) and oceanographic (SST) anomalies that characterize ENSO events across the Pacific Basin over the period of common record (1960–1979). Comparison of the intensity of climatic anomalies at these two sites yields insight into the spatial variability of ENSO events. Isotope records from older corals can provide high-resolution, Pacific-wide reconstructions of ENSO behavior during periods of different climate boundary conditions.

Description: As shown by the work of Dansgaard and his colleagues, climate oscillations of one or so millennia duration punctuate much of glacial section of the Greenland ice cores. These oscillations are characterized by 5°C air temperature changes, severalfold dust content changes and 50 ppm CO2 changes. Both the temperature and CO2 change are best explained by changes in the mode of operation of the ocean. In this paper we provide evidence which suggests that oscillations in surface water conditions of similar duration are present in the record from a deep sea core at 50°N. Based on this finding, we suggest that the Greenland climate changes are driven by oscillations in the salinity of the Atlantic Ocean which modulate the strength of the Atlantic's conveyor circulation.

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

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

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

Description: AGU considers only original scientific contributions that have not been accepted or published elsewhere and are not under consideration by another publisher. A contribution is considered previously published if its data and conclusions are offered for sale or are generally available in other ways to the public. Regardless of the original publication medium, including print, magnetic tape, or microform, such contributions are not eligible for republication in AGU journals or books.

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 seasonal variation of the intrusion of the Philippine Sea Water into the South China Sea was studied by analyzing the historical hydrographic station data in the northern South China Sea and the Philippine Sea. Water masses at 150, 200, and 250 m were classified by discriminant analysis according to their temperature-salinity characteristics. At each depth, most water in the study region was classified into two groups representing the Philippine Sea Water and the South China Sea Water, respectively. The geographic distribution of water masses in the South China Sea shows that the Philippine Sea Water was present along the continental margin south of China between October and January. A westward current in the northern South China Sea in winter was inferred from the distribution of the intrusion water.

Description: Characteristics of water masses were analyzed to study the Kuroshio intrusion into the sea southwest of Taiwan. Hydrographic data were obtained from CTD (conductivity, temperature, and depth) casts during two cruises in May and August 1986. In May, remnants of water intruding from the Kuroshio were found on the continental slope south of the Penghu Channel. By August, these were replaced by water from the South China Sea. During this period, water from the Kuroshio also appeared near the southern tip of Taiwan. The intrusion current reached a depth of at least 500 m and was probably part of a cyclonic circulation in the northern South China Sea. The results support the hypothesis of a seasonal pattern of the intrusion process: intrusion of water from the Kuroshio begins in late summer, intensifies in winter, and ceases by late spring when South China Sea waters again enter this region.

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

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

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

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

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

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

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

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

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

Description: 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: Multichannel seismic reflection data were used to determine the evolutionary history of the forearc region of the central Aleutian Ridge. Since at least late Miocene time this sector of the ridge has been obliquely underthrust 30° west of orthogonal convergence by the northwestward converging Pacific plate at a rate of 80–90 km/m.y. Our data indicate that prior to late Eocene time the forearc region was composed of rocks of the arc massif thinly mantled by slope deposits; the forearc region probably lacked both major depositional basins and a tectonically attached accretionary prism of offscraped oceanic deposits. Beginning in latest Miocene or earliest Pliocene time, a zone of outer-arc structural highs and a forearc basin began to form. Formation of these companion intraarc structures may be linked to the late Neogene growth of an accretionary wedge that formed as the result of the deposition of a thick turbidite wedge in the Aleutian Trench. Initial structures of the zone of outer-arc highs formed as the thickening wedge underran, compressively deformed, and uplifted the seaward edge of the arc massif above a landward dipping backstop thrust. Forearc basin strata ponded arcward of the elevating zone of outer-arc highs. However, most younger structures of the zone of outer-arc highs cannot be ascribed simply to the orthogonal effects of an underrunning wedge. Oblique convergence created a major right-lateral shear zone (the Hawley Ridge shear zone) that longitudinally disrupted the zone of outer-arc highs, truncating the seaward flank of the forearc basin and shearing the southern limb of Hawley Ridge, an exceptionally large antiformal outer-arc high structure. Slivers of forearc basement rocks and overlying strata have been transported along the shear zone that is flanked by differentially elevated structures attributed to localized transpressive and transtensional processes. Uplift of Hawley Ridge may be related to the thickening of the arc massif by westward directed basement duplexes. In addition, the forearc is disrupted by structures transverse to the margin that occur where unusually high-stress accumulations have resulted in the rupture of repeated great earthquakes. It is likely that many ancient active margins evolved in tectonic and depositional settings similar to those of the central Aleutian Ridge. Great structural complexity, including the close juxtaposition of coeval structures recording compression, extension, differential vertical movements, and strike-slip displacement, should be expected, even within areas of generally kindred tectonostratigraphic terranes.

Description: Several geochemical anomalies were observed before the Haichen, Longling, Tangshan, and Songpan earthquakes and their strong aftershocks. They included changes in groundwater radon levels; chemical composition of the groundwater (concentration of Ca++, Mg++, Cl−, SO4= and HCO3− ions); conductivity; and dissolved gases such as H2, CO2, etc. In addition, anomalous changes in water color and quality were observed before these large earthquakes. Before some events gases escaped from the surface, and there were reports of "ground odors" being smelled by local residents. The large amount of radon data can be grouped into long-term and short-term anomalies. The long-term anomalies have a radon emission build up time of from a few months to more than a year. The short-term anomalies have durations from a few hours or less to a few months.

Description: Current measurements from two consecutive yearlong deployments of three moored stations at the western end of the equator in the Atlantic, along 44°W, are used to determine the northwestward flow of warm water in the upper several 100 m and of the southeastward counterflow of North Atlantic Deep Water (NADW). Measurements from three acoustic Doppler current profilers (ADCPs) looking upward from 300 m toward the surface allowed calculation of a time series of upper layer transports over 1 year. Mean transport through the array for the upper 300 m is 23.8 Sv with an annual cycle of only ±3 Sv that has its maximum in June-August and minimum in northern spring. Estimated additional mean northwestward transport in the range 300–600 m is 6.7 Sv, based on moored data and shipboard Pegasus and lowered ADCP profiling. In the depth range 1400–3100 m a current core with maximum annual mean southeastward speed of 30 cm s−1 is found along the continental slope that carries an estimated upper NADW transport of 14.2–17.3 Sv, depending on the extrapolation used between the mooring in the core and the continental slope. This transport is higher than off-equatorial estimates and suggests near-equatorial recirculation at the upper NADW level, in agreement with northwestward mean flow found about 140 km offshore. Below 3100 m and above the 1.8°C isotherm, only a small core of lower NADW flow with speeds of 10–15 cm s−1 is found over the flat part of the basin near 1.5°N, clearly separated from the continental slope by a zone of near-zero mean speeds. Estimated transport of that small current core is about 4.5 Sv, which is significantly below other estimates of near-equatorial transport of lower NADW and suggests that a major fraction of lower NADW may cross the 44°W meridian north of the Ceara Rise. Intraseasonal variability is large, although smaller than observed at 8°N near the western boundary. It occurs at a period of about 1 month when it is dominant in the near-surface records and corresponds to earlier observations in the equatorial zones of all oceans and at a period of about 2 months when it is dominant at the NADW level and could be imported either from the north along the boundary or from the east along the equator. The existence of an annual cycle in the deep currents of a few centimeters per second amplitude, as suggested by high-resolution numerical model results, could neither be proven nor disproven because of the high amount of shorter-period variability.

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

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

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

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

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

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

Description: A non-Fickian physico-chemical model for electrolyte transport in high-ionic strength systems is developed and tested with laboratory experiments with copper sulfate as an example electrolyte. The new model is based on irreversible thermodynamics and uses measured mutual diffusion coefficients, varying with concentration. Compared to a traditional Fickian model, the new model predicts less diffusion and asymmetric diffusion profiles. Laboratory experiments show diffusion rates even smaller than those predicted by our non-Fickian model, suggesting that there are additional, unaccounted for processes retarding diffusion. Ionic diffusion rates may be a limiting factor in transporting salts whose effect on fluid density will in turn significantly affect the flow regime. These findings have important implications for understanding and predicting solute transport in geologic settings where dense, saline solutions occur.

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

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

Description: The supply of limiting nutrients to the low latitude ocean is controlled by physical processes linked to climate variations, but methods for reconstructing past nutrient concentrations in the surface ocean are few and indirect. Here, we present laser ablation mass spectrometry results that reveal annual cycles of P/Ca in a 4-year record from the scleractinian coral Pavona gigantea (mean P/Ca = 118 μmol mol−1). The P/Ca cycles track variations in past seawater phosphate concentration synchronously with skeletal Sr/Ca-derived temperature variations associated with seasonal upwelling in the Gulf of Panamá. Skeletal P/Ca varies seasonally by 2–3 fold, reflecting the timing and magnitude of dissolved phosphate variations. Solution cleaning experiments on drilled coral powders show that over 60% of skeletal P occurs in intracrystalline organic phases. Coral skeleton P/Ca holds promise as a proxy record of nutrient availability on time scales of decades to millennia.

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 carbon isotopic composition of methane emitted by the Alaskan emergent aquatic plants Arctophila fulva, a tundra mid-lake macrophyte, and Carex rostrata, a tundra lake margin macrophyte, was −58.6 ± 0.5 (n=2) and −66.6±2.5 (n= 6) ‰ respectively. The methane emitted by these species was found to be depleted in 13C by 12‰ and 18‰, relative to methane withdrawn from plant stems 1 to 2 cm below the waterline. As the macrophyte-mediated methane flux represented approximately 97% of the flux from these sites, these results suggest the more rapid transport of 12CH4 relative to 13CH4 through plants to the atmosphere. This preferential release of the light isotope of methane, possibly combined with CH4 oxidation, caused the buildup of the heavy isotope within plant stems. Plant stem methane concentrations ranged from 0.2 to 4.0% ( math formula, 1.4; standard deviation (sd), 0.9; n=28) in Arctophila, with an isotopic composition of −46.1±4.3 ‰ (n = 8). Carex stem methane concentrations were lower, ranging from 150 to 1200 ppm ( math formula, 500; standard deviation, 360; n = 8), with an isotopic composition of −48.3±1.4‰ (n=3). Comparisons of the observed isotopic fractionations with those predicted from gas phase effusion and diffusion coefficients suggest a combination of one or both of these gas transport mechanisms with bulk (non-fractionationating) flow.

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

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

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

Description: Until reliable procedures have been developed to preserve the phosphorus contained in particulate matter captured by in situ pumps and sediment traps and until these procedures are applied over a wide range of locations and depths in the sea, indirect methods will have to be used to determine the C/P ratio in marine detritus. We have taken two such approaches: (1) the use of C/N ratios for particulates captured in the upper thermocline in conjunction with 02/P and N/P ratios obtained from deconvolutions of ocean chemical data and (2) regression along isopycnals in the deep‐sea waters free of fossil fuel CO2. While neither approach yields a definitive answer, both suggest that a value of 127 carbon atoms per phosphorus atom would be a more appropriate interim value than that of 106 adopted long ago by A. C. Redfield and his associates.

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

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

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

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

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

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

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

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

Description: The Multitracers Experiment studied a transect of water column, sediment trap, and sediment data taken across the California Current to develop quantitative methods for hindcasting paleoproductivity. The experiment used three sediment trap moorings located 120 km, 270 km, and 630 km from shore at the Oregon/California border in North America. We report here about the sedimentation and burial of particulate organic carbon (Corg) and CaCO3. In order to observe how the integrated CaCO3 and Corg burial across the transect has changed since the last glacial maximum, we have correlated core from the three sites using time scales constrained by both radiocarbon and oxygen isotopes. By comparing surface sediments to a two-and-a-half year sediment trap record, we have also defined the modern preservation rates for many of the labile sedimentary materials. Our analysis of the Corg data indicates that significant amounts (20–40%) of the total Corg being buried today in surface sediments is terrestrial. At the last glacial maximum, the terrestrial Corg fraction within 300 km of the coast was about twice as large. Such large fluxes of terrestrial Corg obscure the marine Corg record, which can be interpreted as productivity. When we corrected for the terrestrial organic matter, we found that the mass accumulation rate of marine Corg roughly doubled from the glacial maximum to the present. Because preservation rates of organic carbon are high in the high sedimentation rate cores, corrections for degradation are straightforward and we can be confident that organic carbon rain rate (new productivity) also doubled. As confirmation, the highest burial fluxes of other biogenic components (opal and Ba) also occur in the Holocene. Productivity off Oregon has thus increased dramatically since the last glacial maximum. CaCO3 fluxes also changed radically through the deglaciation; however, they are linked not to CaCO3 production but rather to changes in deepwater carbonate chemistry between 18 Ka and now.

Description: Axial volcanic ridges (AVRs) are found on most slow-spreading mid-ocean ridges and are thought to be the main locus of volcanism there. In this study we present high-resolution mapping of a typical, well-defined AVR on the Mid-Atlantic Ridge at 45°N. The AVR is characterized by “hummocky terrain,” composed typically of hummocks with pillowed or elongate pillowed flanks with pillowed or lobate lava flow summits, often with small haystacks sitting on their highest points. The AVR is surrounded by several areas of “flat seafloor,” composed of lobate and sheet lava flows. The spatial and morphological differences between these areas indicate different eruption processes operating on and off the AVR. Volcanic fissures are found all around and on the AVR, although those with the greatest horizontal displacement are found on the ridge crest and flat seafloor. Clusters of fissures may represent volcanic vents. Extremely detailed comparisons of sediment coverage and examination of contact relations around the AVR suggest that many of the areas of flat seafloor are of a similar age or younger than the hummocky terrain of the AVR. Additionally, all the lavas surveyed have similar degrees of sediment cover, suggesting that the AVR was either built or resurfaced in the same 50 ka time frame as the flat seafloor.

Description: This review is intended to cover the principal developments that have occurred within the last six years in the paleomagnetic study of marine sediments. Recent work utilizing the reflecting-light microscope indicates that detrital high-temperature Fe-Ti oxides are probably responsible for most of the magnetic remanence in marine sediments. These minerals possess a spectrum of coercivities that makes it necessary to use alternating-field—demagnetization techniques to isolate stable components. It is possible to use the standard magnetic stratigraphy for the last 4 m.y. of earth history derived from terrestrial lavas. Using the ages of the magnetic boundaries from this time scale it is possible by extrapolation and interpolation to better determine the ages of the major events. The ages of these events in increasing age are Jaramillo, 0.87 to 0.92 m.y.; Olduvai, 1.71 to 1.86 m.y.; Kaena, 2.82 to 2.90 m.y.; Mammoth, 3.0 to 3.085 m.y.; Cochiti, 3.72 to 3.82 m.y.; Nunivak, 3.97 to 4.14 m.y.; ‘c’ event of the Gilbert series, 4.33 to 4.65 m.y. Through the use of long cores from the central Pacific and through correlation using fossil datums, it has been possible to extend the magnetic stratigraphy back to the upper middle Miocene to magnetic epoch 5. It is concluded that very short magnetic events are probably short-term excursions of the field and not true magnetic events. It is shown that the field of the earth averages to an axial-dipole field within a period of 27,000 years and that the field over the last two million years has acted as a geocentric axial dipole. The evidence shows that when reversals of the dipole occur, the values of the reversed inclination are not significantly different from the normal values. The use of magnetic stratigraphy in marine geology has opened up a new era in study of sedimentary processes and evolution of marine organisms.

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

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

Description: Digital hydrographic data combined with satellite thermal infrared and visible band remote sensing provide a synoptic climatological view of the shallow planktonic environment. This paper uses wind, hydrographic, and ocean remote sensing data to examine southwest monsoon controls on the foraminiferal faunal composition of Recent seafloor sediments of the northwestern Arabian Sea. Ekman pumping resulting in open-ocean upwelling and coastal upwelling create two distinctly different mixed layer plankton environments in the northwestern Arabian Sea during the summer monsoon. Open-sea upwelling to the northwest of the mean July position of the Findlater Jet axis yields a mixed layer environment with temperatures of less than 25°C to about 26.5°C, phytoplankton pigment concentrations between 1.5 and 5.0 mg/m³, and mixed layer depths less than 50 m. Convergence in the Ekman layer in the central Arabian Sea drives the formation of a mixed layer that is greater than 50 m thick, warmer than about 26.5°C, and has phytoplankton pigment concentrations generally below 2.0 mg/m³. Coastal upwelling creates an extremely eutrophic plankton environment that persists over and immediately adjacent to the Omani shelf and undergoes significant offshore transport only within topographically induced coastal squirts. The foraminiferal faunal composition of upper Pleistocene deep-sea sediments of the northwestern Arabian Sea are mainly controlled by vertical nutrient fluxes caused by Ekman pumping, not coastal upwelling. Transfer functions for late Pleistocene mixed layer depth, temperature, and chlorophyll have been obtained through factor analysis and nonlinear multiple regression between late summer mixed layer environment and Recent sediment faunal observations.

Description: This paper provides a detailed hydrographic climatology for the shallow northwestern Arabian Sea prior to and during the southwest monsoon, presented as multiple-year composite vertical hydrographic sections based on National Oceanographic Data Center historical ocean station data. Temperature and salinity measurements are used to infer the water masses present in the upper 500 m. The hydrographic evolution depicted on bimonthly sections is inferred to result from wind-driven physical processes. In the northwestern Arabian Sea the water mass in the upper 50 m is the Arabian Sea Surface Water. Waters from 50 to 500 m are formed by mixing of Arabian Sea Surface Water with Antarctic and Indonesian intermediate waters. The inflow of Persian Gulf Water does not significantly influence the hydrography of the northwestern Arabian Sea along the Omani coast. Nitrate has a high inverse correlation with temperature and oxygen in the premonsoon thermocline in the depth interval 50–150 m. During the southwest monsoon, coastal upwelling off Oman and adjacent offshore upward Ekman pumping alter the shallow hydrography.

Description: The Norwegian Atlantic Current (NwAC) and its eddy field are examined using data from surface drifters. The data set used spans nearly 20 years, from June 1991 to December 2009. The results are largely consistent with previous estimates, which were based on data from the first decade only. With our new data set, statistical analysis of the mean fields can be calculated with larger confidence. The two branches of the NwAC, one over the continental slope and a second further offshore, are clearly captured. The Norwegian Coastal Current is also resolved. In addition, we observe a semipermanent anticylonic eddy in the Lofoten Basin, a feature seen previously in hydrography and in models. The eddy kinetic energy (EKE) is intensified along the path of the NwAC, with the largest values occurring in the Lofoten Basin. The strongest currents, exceeding 100 cm s−1, occur west of Lofoten. Lateral diffusivities were computed in five domains and ranged from 1–5 × 107 cm2 s−1. The Lagrangian integral time and space scales are 1–2 days and 7–23 km, respectively. The data set allows studies of seasonal and interannual variations as well. The strongest seasonal signal is in the NwAC itself, as the mean flow strengthens by approximately 20% in winter. The EKE and diffusivities on the other hand do not exhibit consistent seasonality in the sampled regions. There are no consistent indications of changes in either the mean or fluctuating surface velocities between the 1990s and 2000s.

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

Description: The biological variability of the northwestern Arabian Sea during the 1979 southwest monsoon has been investigated by the synthesis of satellite ocean color remote sensing with analysis of in situ hydrographic and meteorological data sets and the results of wind-driven modeling of upper ocean circulation. The phytoplankton bloom in the northwestern Arabian Sea peaked during August-September, extended from the Oman coast to about 65°E, and lagged the development of open-sea upwelling by at least 1 month. In total, the pigment distributions, hydrographic data, and model results all suggest that the bloom was driven by spatially distinct upward nutrient fluxes to the euphotic zone forced by the physical processes of coastal upwelling and offshore Ekman pumping. Coastal upwelling was evident from May through September, yielded the most extreme concentrations of phytoplankton biomass, and along the Arabian coast was limited to the continental shelf in the promotion of high concentrations of phytoplankton. Upward Ekman pumping to the northwest of the Somali Jet axis stimulated the development of a broad open-sea phytoplankton bloom oceanward of the Oman shelf. Vertical mixing during the 1979 southwest monsoon was apparently not a primary cause of the regional-scale phytoplankton bloom.

Description: Samples collected at hourly intervals on May 18–19, 1980, at three sites 200 km downwind from Mount St. Helens, have made possible a detailed reconstruction of the conditions that contribute to the compositional heterogeneity of mineral and glass components observed in distal tephra layers. The air fall tephra deposited at the sites during the first 7 hours of the May 18 eruption is mostly coarse grained, microlite-rich, nonjuvenile glass and feldspar. Grain-size maxima in this initial tephra can be related to the cataclysmic blast at 0832 and a subsequent pulse of the eruption at 1200. Juvenile, microlite-free glass increases in relative abundance at the sampling sites beginning at about 1900. Such a change between nonjuvenile and juvenile tephra can be related to a 5-km increase in column height associated with the last major pulse of the eruption which occurred at 1700 at the volcano. Electron microprobe study of both microlite-rich and microlite-free pumice in the time series samples reveals significant compositional differences. Interstitial glass in nonjuvenile pumice deposited during the first few hours at the sampling sites is enriched in SiO2 and K2O and depleted in TiO2, FeO*, and MgO relative to juvenile glass. By comparison, major element composition of the least evolved juvenile glass sampled during the last several hours of the eruption displays a slight trend toward less evolved composition. Least squares calculations suggest that the more evolved character of the nonjuvenile glass can be explained by greater fractional crystallization brought about by enhanced cooling in a cryptodome prior to eruption, whereas the temporal changes observed in juvenile glass composition during the last several hours of the eruption suggest the presence of a small, slightly zoned magma chamber at depth. Electron microprobe study of glass-coated ilmenites, magnetites, and plagioclases provides the following estimates of the physical conditions in this reservoir: 865°±50°C, PH2O = 2.2 kbar and -log ƒO2 = 11.7. Analyses of bulk pumice, glass and selected mineral phases from May 25, June 12, July 22, and October 16–18 pumices erupted from Mount St. Helens indicate that the bulk pumice (magma) compositions have become slightly more andesitic with time, while mineral and co-existing glass compositions have changed significantly in post-May 18 eruptions with both being more highly evolved than those associated with the May 18 eruption. An application of the magnetite-ilmenite geothermometer to June 12 and July 22 samples indicates temperatures of 919°±30°C and 930°±50°C, respectively. Least squares calculations suggest that such evolved post-May 18 glass and mineral phases can be derived by fractional crystallization of a magma composition like bulk May 18 pumice into approximately 50% crystals and 50% residual liquid. Such partitioning between crystals and residual liquid appears to have occurred on the scale of centimeters and is interpreted as a consequence of accelerated crystallization under reduced water pressure.

Description: What is the relationship between volcanic eruptions and climate change? More than 200 years after the connection was first proposed, it remains a thorny question. This article provides a brief historical overview of the problem and a review of the various data bases used in evaluating volcanic events and associated climatic change. We use the term “climate” to describe changes in the atmosphere over wide regions for periods of several months and longer. We use “weather” to describe shorter-term, variable atmospheric fluctuations experienced over more restricted areas. We appraise the present state of knowledge and highlight some pitfalls involved in using available information. Cautiously, we suggest future avenues for study, including the possibility of “volcanic winters,” or severe eruption-induced coolings.

Description: The Stratigraphie record from both deep-sea and shallow-water depositional environments Indicates that during late Aptian through Cenomanian time (1) global climates were considerably warmer than at present; (2) latitudinal gradients of atmospheric and oceanic temperatures were considerably less than at present; (3) rates of accumulation of organic matter of both marine and terrestrial origin were as high as or higher than during any other interval in the Mesozoic or Cenozoic; (4) the rate and volume of accumulation of CaC02 in the deep sea were reduced in response to a marked shoaling of the carbonate compensation depth; (5) seafloor spreading rates were somewhat more rapid than at any other time in the Cretaceous or Cenozoic; (6) off-ridge volcanism was intense and widespread, particularly in the ancestral Pacific Ocean basin; and (7) sea level was relatively high, forming widespread areas of shallow shelf seas. A marked increase in the rate of C02 outgassing due to volcanic activity between about 110 and 70 m.y. ago may have resulted in a buildup of atmospheric C02. A significant fraction of this atmospheric C02 may have been reduced by an increase in the production and burial of terrestrial organic carbon. Some excess C02 may have been consumed by marine algal photosynthesis, but marine productivity apparently was low during the Aptian-Albian relative to terrestrial productivity. Terrestrial productivity also may have been stimulated by increased rainfall that resulted from a warm global climate and increased marine transgression as well as by the higher C02.

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

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

Description: In accretionary wedges, often morphologically similar sedimentary intrusions, when observed by remote geophysical means, may have one of two quite different driving mechanisms and a highly variable significance for the regional hydrogeologic picture. For example, mud diapirs are driven by buoyancy forces that arise from bulk density contrasts. In them, mud and pore fluid upwell en masse and fluid migration is a related (sometimes important) but generally subsidiary process. In contrast, diatremes contain sediments fluidized during rapid fluid advection and are forcibly and directly driven by the hydrogeologic system. The nature of fluid input from local and exotic source regions can, therefore, strongly affect sedimentary intrusive processes and vice versa. This complicates the process of defining the main features of the hydrogeological systems operating in accretionary wedges. Focused vertical advection through steep sided (piercement) mud diapirs requires conduit systems, otherwise flow will be diffuse and directed more horizontally out of the low-permeability mud mass. However, where the permeability of the overburden is less than that of the diapir, the whole diapir may act as a conduit. Apart from this special case, conduits will be associated with highly anisotropic scaly fabrics that can sometimes develop in the marginal shear zone of diapirs. Scaly fabrics form during deformation and compaction of a mud matrix under conditions of constant or increasing effective stress. However, the effective stress path can be complex as it is both controlled by the relative rates of upward intrusion and burial (by sedimentation and/or structural thickening) and the hydrogeologic system. Due to this, it appears likely that even in a geographically related group of diapirs, effective stress histories will vary widely between intrusions so that some can form advective pathways for fluids and some cannot. Mixed systems of behavior may also be present with local diatremes developing within diapirs above the terminations of conduit systems and rapidly expanding methane gas pockets. The potentially heterogeneous near-surface behavior may be why the surface manifestations of sedimentary intrusions are so variable when observed in the field. Diatremes can also form separately as large primary features above any structural or stratigraphic conduit that rapidly expels water or gas into the base of an unlithified sediment column. When active, large diatremes require enormous quantities of fluid (water or gas) to drive them, particularly if they are long lived features and hence are a direct indication of at least an episodically vigorous hydrogeologic system.

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

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

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

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

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

Description: The well-known caldera of Thira (Santorini), Greece, was not formed during a single eruption but is composed of two overlapping calderas superimposed upon a complex volcanic field that developed along a NE trending line of vents. Before the Minoan eruption of 1400 B.C., Thira consisted of three Java shields in the northern half of the island and a flooded depression surrounded by tuff deposits in the southern half. Andesitic lavas formed the overlapping shields of the north and were contemporaneous with and, in many places, interbedded with the southern tuff deposits. Although there appears to be little difference between the composition of magmas erupted, differences in eruption style indicate that most of the activity in the northern half of the volcanic field was subaerial, producing lava flows, whereas in the south, eruptions within a flooded depression produced a sequence of mostly phreatomagmatic tuffs. Many of these tuffs are plastered onto the walls of what appears to have been an older caldera, most probably associated with an eruption of rhyodacitic tephra 100,000 years ago. The Minoan eruption of about 1400 B.C. had four distinct phases, each reflecting a different vent geometry and eruption mechanism. The Minoan activity was preceded by minor eruptions of fine ash. (1) The eruption began with a Plinian phase, from subaerial vent(s) located on the easternmost of the lava shields. (2) Vent(s) grew toward the SW into the flooded depression. Subsequent activity deposited large-scale base surge deposits during vent widening by phreatomagmatic activity. (3) The third eruptive phase was also phreatomagmatic and produced 60% of the volume of the Minoan Tuff. This activity was nearly continuous and formed a large featureless tuff ring with poorly defined bedding. This deposit contains 5–40% lithic fragments that are typical of the westernmost lava shield and appears to have been erupted when caldera collapse began. (4) The last phase consisted of eruption of ignimbrites from vent(s) on the eastern shield, not yet involved in collapse. Collapse continued after eruption of the ignimbrites with foundering of the eastern half of the caldera. Total volume of the collapse was about 19 km3, overlapping the older caldera to form the caldera complex visible today. Intracaldera eruptions have formed the Kameni Islands along linear vents concomitant with vents that may have been sources for the Minoan Tuff.

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: Magnetic lineation mapping in the western central Pacific has revealed a pair of opposite-sensed, fanned lineation patterns that define the accretionary boundaries of the fossil Magellan microplate. This tectonic synthesis results from extensive magnetic mapping of two new lineation patterns over a large area and extended mapping of previously identified lineations. The entire evolutionary history of the Magellan microplate is well constrained to a 9-m.y. period in the Early Cretaceous by synchronous spreading patterns and associated geologic data. During this period the microplate grew and evolved as a generally rectangular structure to a final size of 700 km×600 km with spreading centers on two opposing sides and transform faults on the other two sides. The lifetime and size of the Magellan microplate are somewhat longer and larger, respectively, than presently active microplates on the East Pacific Rise. However, these modern structures are still evolving and growing, and the tectonic behavior of the modern and Cretaceous systems appears to be similar. Study of both active and fossilized microplates should provide additional insights on their common tectonic histories. In particular, we show that the Magellan Trough spreading center behaved as an asymmetric accretionary plate boundary that can be described with two separate poles of motion very close to this spreading center during much of its history. The Magellan Trough spreading center then failed as a result of a larger ridge reorganization at the triple junction of the Pacific, Farallon, and Phoenix plates at Ml0N time. Microplate activity ceased when the microplate became welded to the Pacific plate at M9 time.

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

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

Description: We present a plate kinematic evolution of the South Atlantic which is based largely on the determination of the equatorial fracture zone trends between the African and South American continental margins. Four main opening phases are dated by oceanic magnetic anomalies, notably MO, A34, and A13, and are correlated with volcanism and tectonic events on land around the South Atlantic Ocean. The Ceara and Sierra Leone rises are probably of oceanic origin and were created 80 m.y. ago or later in their present-day positions with respect to South America and Africa.

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

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

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

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

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

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

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

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

Description: [1] While recent studies have confirmed the ecological importance of vitamin B12, it is unclear whether the production of this vitamin could be limited by dissolved Co, a trace metal required for B12 biosynthesis, but found at only subnanomolar concentrations in the open ocean. Herein, we demonstrate that the spatial distribution of dissolved B12 (range: 0.13–5 pmol L−1) in the North Atlantic Ocean follows the abundance of total dissolved Co (range: 15–81 pmol L−1). Similar patterns were observed for bacterial productivity (range: 20–103 pmol 3H leucine L−1 hr−1) and algal biomass (range: 0.4–3.9 μg L−1). In contrast, vitamin B1 concentrations (range: 0.7–30 pM) were decoupled from both Co and B12 concentrations. Cobalt amendment experiments carried out in low-dissolved Co waters (∼20 pmol L−1) enhanced B12 production two-fold over unamended controls. This study provides evidence that B12 synthesis could be limited by the availability of Co in some regions of the world ocean.

Description: We present the first comprehensive investigation of the concentrations, fluxes and sources of aerosol trace elements over the Gulf of Aqaba. We found that the mean atmospheric concentrations of crustally derived elements such as Al, Fe and Mn (1081, 683, and 16.7 ng m�3) are about 2–3 times higher than those reported for the neighboring Mediterranean area. This is indicative of the dominance of the mineral dust component in aerosols over the Gulf. Anthropogenic impact was lower in comparison to the more heavily populated areas of the Mediterranean. During the majority of time (69%) the air masses over the Gulf originated from Europe or Mediterranean Sea areas delivering anthropogenic components such as Cu, Cd, Ni, Zn, and P. Airflows derived from North Africa in contrast contained the highest concentrations of Al, Fe, and Sr but generally lower Cu, Cd, Ni, Zn, and P. Relatively high Pb, Ni, and V were found in the local and Arabian airflows suggesting a greater influence of local emission of fuel burning. We used the data and the measured trace metal seawater concentrations to calculate residence times of dissolved trace elements in the upper 50 m surface water of the Gulf (with respect to atmospheric input) and found that the residence times for most elements are in the range of 5–37 years while Cd and V residence times are longer.

Description: Seismic techniques provide the highest-resolution measurements of the structure of the crust and have been conducted on a worldwide basis. We summarize the structure of the continental crust based on the results of seismic refraction profiles and infer crustal composition as a function of depth by comparing these results with high-pressure laboratory measurements of seismic velocity for a wide range of rocks that are commonly found in the crust. The thickness and velocity structure of the crust are well correlated with tectonic province, with extended crust showing an average thickness of 30.5 km and orogens an average of 46.3 km. Shields and platforms have an average crustal thickness nearly equal to the global average. We have corrected for the nonuniform geographical distribution of seismic refraction profiles by estimating the global area of each major crustal type. The weighted average crustal thickness based on these values is 41.1 km. This value is 10% to 20% greater than previous estimates which underrepresented shields, platforms, and orogens. The average compressional wave velocity of the crust is 6.45 km/s, and the average velocity of the uppermost mantle (Pn velocity) is 8.09 km/s. We summarize the velocity structure of the crust at 5-km depth intervals, both in the form of histograms and as an average velocity-depth curve, and compare these determinations with new measurements of compressional wave velocities and densities of over 3000 igneous and metamorphic rock cores made to confining pressures of 1 GPa. On the basis of petrographic studies and chemical analyses, the rocks have been classified into 29 groups. Average velocities, densities, and standard deviations are presented for each group at 5-km depth intervals to crustal depths of 50 km along three different geotherms. This allows us to develop a model for the composition of the continental crust. Velocities in the upper continental crust are matched by velocities of a large number of lithologies, including many low-grade metamorphic rocks and relatively silicic gneisses of amphibolite facies grade. In midcrustal regions, velocity gradients appear to originate from an increase in metamorphic grade, as well as a decrease in silica content. Tonalitic gneiss, granitic gneiss, and amphibolite are abundant midcrustal lithologies. Anisotropy due to preferred mineral orientation is likely to be significant in upper and midcrustal regions. The bulk of the lower continental crust is chemically equivalent to gabbro, with velocities in agreement with laboratory measurements of mafic granulite. Garnet becomes increasingly abundant with depth, and mafic garnet granulite is the dominant rock type immediately above the Mohorovicic discontinuity. Average compressional wave velocities of common crustal rock types show excellent correlations with density. The mean crustal density calculated from our model is 2830 kg/m3, and the average SiO2 content is 61.8%.

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