ALBERT

Description: A new 21.3m firn core was drilled in 2015 at a coastal Antarctic high-accumulation site in Adélie Land (66.78◦ S; 139.56◦ E, 602 m a.s.l.), named Terre Adélie 192A (TA192A). The mean isotopic values (−19.3 ‰ ± 3.1 ‰ for δ18O and 5.4 ‰±2.2 ‰ for deuterium excess) are consistent with other coastal Antarctic values. No significant isotope–temperature relationship can be evidenced at any timescale. This rules out a simple interpretation in terms of local temperature. An observed asymmetry in the δ18O seasonal cycle may be explained by the precipitation of air masses coming from the eastern and western sectors in autumn and winter, recorded in the d-excess signal showing outstanding values in austral spring versus autumn. Significant positive trends are observed in the annual d-excess record and local sea ice extent (135–145◦ E) over the period 1998–2014. However, process studies focusing on resulting isotopic compositions and particularly the deuterium excess–δ18O relationship, evidenced as a potential fingerprint of moisture origins, as well as the collection of more isotopic measurements in Adélie Land are needed for an accurate interpretation of our signals.

Description: The timing and intensity of snowmelt processes on sea ice are key drivers determining the seasonal sea-ice energy and mass budgets. In the Arctic, satellite passive microwave and radar observations have revealed a trend towards an earlier snowmelt onset during the last decades, which is an important aspect of Arctic amplification and sea ice decline. Around Antarctica, snowmelt on perennial ice is weak and very different than in the Arctic, with most snow surviving the summer. Here we compile time series of snowmelt-onset dates on seasonal and perennial Antarctic sea ice from 1992 to 2014/15 using active microwave observations from European Remote Sensing Satellite (ERS-1/2), Quick Scatterometer (QSCAT) and Advanced Scatterometer (ASCAT) radar scatterometers. We define two snowmelt transition stages: A weak backscatter rise indicating the initial warming and destructive metamorphism of the snowpack (pre-melt), followed by a rapid backscatter rise indicating the onset of thaw-freeze cycles (snowmelt). Results show large interannual variability with an average pre-melt onset date of 29 November and melt onset of 10 December, respectively, on perennial ice, without any significant trends over the study period, consistent with the small trends of Antarctic sea ice extent. There was a latitudinal gradient from early snowmelt onsets in mid-November in the northern Weddell Sea to late (end-December) or even absent snowmelt conditions in the southern Weddell Sea. We show that QSCAT Ku-band (13.4 GHz signal frequency) derived pre-melt and snowmelt onset dates are earlier by 20 and 18 days, respectively, than ERS and ASCAT C-band (5.6 GHz) derived dates. This offset has been considered when constructing the time series. Snowmelt onset dates from passive microwave observations (37 GHz) are later by 14 and 6 days than those from the scatterometers, respectively. Based on these characteristic differences between melt onset dates observed by different microwave wavelengths, we developed a conceptual model which illustrates how the seasonal evolution of snow temperature profiles may affect different microwave bands with different penetration depths. These suggest that future multi-frequency active/passive microwave satellite missions could be used to resolve melt processes throughout the vertical snow column of thick snow on perennial Antarctic sea ice.

Description: Although quantitative isotope data from speleothems has been used to evaluate isotope-enabled model simulations, currently no consensus exists regarding the most appropriate methodology through which to achieve this. A number of modelling groups will be running isotope-enabled palaeoclimate simulations in the framework of the Coupled Model Intercomparison Project Phase 6, so it is timely to evaluate different approaches to using the speleothem data for data–model comparisons. Here, we illustrate this using 456 globally distributed speleothem δ18O records from an updated version of the Speleothem Isotopes Synthesis and Analysis (SISAL) database and palaeoclimate simulations generated using the ECHAM5-wiso isotope-enabled atmospheric circulation model. We show that the SISAL records reproduce the first-order spatial patterns of isotopic variability in the modern day, strongly supporting the application of this dataset for evaluating model-derived isotope variability into the past. However, the discontinuous nature of many speleothem records complicates the process of procuring large numbers of records if data–model comparisons are made using the traditional approach of comparing anomalies between a control period and a given palaeoclimate experiment. To circumvent this issue, we illustrate techniques through which the absolute isotope values during any time period could be used for model evaluation. Specifically, we show that speleothem isotope records allow an assessment of a model’s ability to simulate spatial isotopic trends. Our analyses provide a protocol for using speleothem isotope data for model evaluation, including screening the observations to take into account the impact of speleothem mineralogy on δ18O values, the optimum period for the modern observational baseline and the selection of an appropriate time window for creating means of the isotope data for palaeo-time-slices.

Description: We present here the first results, for the preindustrial and mid-Holocene climatological periods, of the newly developed isotope-enhanced version of the fully coupled Earth system model MPI-ESM, called hereafter MPI-ESM-wiso. The water stable isotopes H16O, H18O and HDO have been implemented into all components of the coupled model setup. The mid-Holocene provides the opportunity to evaluate the model response to changes in the seasonal and latitudinal distribution of insolation induced by different orbital forcing conditions. The results of our equilibrium simulations allow us to evaluate the performance of the isotopic model in simulating the spatial and temporal variations of water isotopes in the different compartments of the hydrological system for warm climates. For the preindustrial climate, MPI-ESM-wiso reproduces very well the observed spatial distribution of the isotopic content in precipitation linked to the spatial variations in temperature and precipitation rate. We also find a good model–data agreement with the observed distribution of isotopic composition in surface seawater but a bias with the presence of surface seawater that is too 18O-depleted in the Arctic Ocean. All these results are improved compared to the previous model version ECHAM5/MPIOM. The spatial relationships of water isotopic composition with temperature, precipitation rate and salinity are consistent with observational data. For the preindustrial climate, the interannual relationships of water isotopes with temperature and salinity are globally lower than the spatial ones, consistent with previous studies. Simulated results under mid-Holocene conditions are in fair agreement with the isotopic measurements from ice cores and continental speleothems. MPI-ESM-wiso simulates a decrease in the isotopic composition of precipitation from North Africa to the Tibetan Plateau via India due to the enhanced monsoons during the mid-Holocene. Over Greenland, our simulation indicates a higher isotopic composition of precipitation linked to higher summer temperature and a reduction in sea ice, shown by positive isotope–temperature gradient. For the Antarctic continent, the model simulates lower isotopic values over the East Antarctic plateau, linked to the lower temperatures during the mid-Holocene period, while similar or higher isotopic values are modeled over the rest of the continent. While variations of isotopic contents in precipitation over West Antarctica between mid-Holocene and preindustrial periods are partly controlled by changes in temperature, the transport of relatively 18O-rich water vapor near the coast to the western ice core sites could play a role in the final isotopic composition. So, more caution has to be taken about the reconstruction of past temperature variations during warm periods over this area. The coupling of such a model with an ice sheet model or the use of a zoomed grid centered on this region could help to better describe the role of the water vapor transport and sea ice around West Antarctica. The reconstruction of past salinity through isotopic content in sea surface waters can be complicated for regions with strong ocean dynamics, variations in sea ice regimes or significant changes in freshwater budget, giving an extremely variable relationship between the isotopic content and salinity of ocean surface waters over small spatial scales. These complicating factors demonstrate the complexity of interpreting water isotopes as past climate signals of warm periods like the mid-Holocene. A systematic isotope model intercomparison study for further insights on the model dependency of these results would be beneficial.

Description: In this study, we present results obtained from modelling the mid-Pliocene warm period using the Community Earth System Models (COSMOS, version: COSMOS-landveg r2413, 2009) with the two different sets of boundary conditions prescribed for the two phases of the Pliocene Model Intercomparison Project (PlioMIP). Boundary conditions, model forcing, and modelling methodology of the two phases of PlioMIP, tagged PlioMIP1 and PlioMIP2,differ considerably in palaeogeography, in particular with regards to the state of ocean gateways, ice-masks, vegetation and topography. Further differences between model setups as suggested for PlioMIP1 and PlioMIP2 consider updates to the concentration of atmospheric carbon dioxide (CO2), that is specified as 405 and 400 parts per million by volume (ppmv) for PlioMIP1 and PlioMIP2, respectively, as well as minor differences in the concentrations of methane (CH4) and nitrous oxide (N2O) due to changes in the protocol of the Paleoclimate Model Intercomparison Project (PMIP) from phase 3 to phase 4. With this manuscript, we bridge the gap between our contributions to PlioMIP1 (Stepanek and Lohmann, 2012) and PlioMIP2 (Stepanek et al., 2019). We highlight some of the effects that differences in the chosen Mid-Pliocene model setup (PlioMIP2 vs. PlioMIP1) have on the climate state as derived with the COSMOS, as this information will be valuable in the framework of the model-model and model-data-comparison within PlioMIP2. We evaluate the model sensitivity to improved mid-Pliocene boundary conditions using PlioMIP’s core mid-Pliocene experiments for PlioMIP1 and PlioMIP2, and present further simulations where we test model sensitivity to variations in palaeogeography, orbit and concentration of CO2. Firstly,we highlight major changes in boundary conditions from PlioMIP1 to PlioMIP2 and also the limitations recorded from the initial effort. The results derived from of our simulations show that COSMOS simulates a mid-Pliocene climate state that is 0.08 K colder in PlioMIP2, if compared to PlioMIP1. On one hand, high-latitude warming,which is supported by proxy evidence of the mid-Pliocene, is underestimated in simulations of both PlioMIP1 andPlioMIP2. On the other hand, spatial variations in surface air temperature (SAT), sea surface temperature (SST) as well as the distribution of sea ice suggest improvement of simulated SAT and SST in PlioMIP2 if employing the updated palaeogeography. The PlioMIP2 Mid-Pliocene simulation produces warmer SSTs in the Arctic and North Atlantic Ocean than derived from the respective PlioMIP1 climate state. The difference in prescribed CO2accountsfor 1.1 K of warming in the Arctic, leading to an ice-free summer in the PlioMIP1 simulation, and a quasi-ice-free summer in PlioMIP2. Furthermore, employing different orbital forcings in simulating the Mid-Pliocene lead to pronounced annual and seasonal variations, which is not accounted for by marine and terrestrial reconstruction of the time-slice.

Description: The Antarctic temperature changes over the past millennia remain more uncertain than in many other continental regions. This has several origins: (1) the number of high-resolution ice cores is small, in particular on the East Antarctic plateau and in some coastal areas in East Antarctica; (2) the short and spatially sparse instrumental records limit the calibration period for reconstructions and the assessment of the methodologies; (3) the link between isotope records from ice cores and local climate is usually complex and dependent on the spatial scales and timescales investigated. Here, we use climate model results, pseudo-proxy experiments and data assimilation experiments to assess the potential for reconstructing the Antarctic temperature over the last 2 millennia based on a new database of stable oxygen isotopes in ice cores compiled in the frame- work of Antarctica2k (Stenni et al., 2017). The well-known covariance between δ18O and temperature is reproduced in the two isotope-enabled models used (ECHAM5/MPI-OM and ECHAM5-wiso), but is generally weak over the different Antarctic regions, limiting the skill of the reconstructions. Furthermore, the strength of the link displays large variations over the past millennium, further affecting the potential skill of temperature reconstructions based on statistical methods which rely on the assumption that the last decades are a good estimate for longer temperature reconstructions. Using a data assimilation technique allows, in theory, for changes in the δ18O–temperature link through time and space to be taken into account. Pseudoproxy experiments confirm the benefits of using data assimilation methods instead of statistical methods that provide reconstructions with unrealistic variances in some Antarctic subregions. They also confirm that the relatively weak link between both variables leads to a limited potential for reconstructing temperature based on δ18O. However, the reconstruction skill is higher and more uniform among reconstruction methods when the reconstruction target is the Antarctic as a whole rather than smaller Antarctic subregions. This consistency between the methods at the large scale is also observed when reconstructing temperature based on the real δ18O regional composites of Stenni et al. (2017). In this case, temperature reconstructions based on data assimilation confirm the long-term cooling over Antarctica during the last millennium, and the later onset of anthropogenic warming compared with the simulations without data assimilation, which is especially visible in West Antarctica. Data assimilation also allows for models and direct observations to be reconciled by reproducing the east–west contrast in the recent temperature trends. This recent warming pattern is likely mostly driven by internal variability given the large spread of individual Paleoclimate Modelling Intercomparison Project (PMIP)/Coupled Model Intercomparison Project (CMIP) model realizations in simulating it. As in the pseudoproxy framework, the reconstruction methods perform differently at the subregional scale, especially in terms of the variance of the time series produced. While the potential benefits of using a data assimilation method instead of a statistical method have been highlighted in a pseudoproxy framework, the instrumental series are too short to confirm this in a realistic setup.

Description: Ice-wedge polygons are common features of lowland tundra in the continuous permafrost zone and prone to rapid degradation through melting of ground ice. There are many interrelated processes involved in ice-wedge thermokarst and it is a major challenge to quantify their influence on the stability of the permafrost underlying the landscape. In this study we used a numerical modelling approach to investigate the degradation of ice wedges with a focus on the influence of hydrological conditions. Our study area was Samoylov Island in the Lena River delta of northern Siberia, for which we had in situ measurements to evaluate the model. The tailored version of the CryoGrid 3 land surface model was capable of simulating the changing microtopography of polygonal tundra and also regarded lateral fluxes of heat, water, and snow. We demonstrated that the approach is capable of simulating ice-wedge degradation and the associated transition from a low-centred to a high-centred polygonal microtopography. The model simulations showed ice-wedge degradation under recent climatic conditions of the study area, irrespective of hydrological conditions. However, we found that wetter conditions lead to an earlier onset of degradation and cause more rapid ground subsidence. We set our findings in correspondence to observed types of ice-wedge polygons in the study area and hypothesized on remaining discrepancies between modelled and observed ice-wedge thermokarst activity. Our quantitative approach provides a valuable complement to previous, more qualitative and conceptual, descriptions of the possible pathways of ice-wedge polygon evolution. We concluded that our study is a blueprint for investigating thermokarst landforms and marks a step forward in understanding the complex interrelationships between various processes shaping ice-rich permafrost landscapes.

Description: Warming of the Arctic led to an increase in permafrost temperatures by about 0.3 �C during the last decade. Permafrost warming is associated with increasing sediment water content, permeability, and diffusivity and could in the long term alter microbial community composition and abundance even before permafrost thaws. We studied the long-term effect (up to 2500 years) of submarine permafrost warming on microbial communities along an onshore–offshore transect on the Siberian Arctic Shelf displaying a natural temperature gradient of more than 10 �C. We analysed the in situ development of bacterial abundance and community composition through total cell counts (TCCs), quantitative PCR of bacterial gene abundance, and amplicon sequencing and correlated the microbial community data with temperature, pore water chemistry, and sediment physicochemical parameters. On timescales of centuries, permafrost warming coincided with an overall decreasing microbial abundance, whereas millennia after warming microbial abundance was similar to cold onshore permafrost. In addition, the dissolved organic carbon content of all cores was lowest in submarine permafrost after millennial-scale warming. Based on correlation analysis, TCC, unlike bacterial gene abundance, showed a significant rank-based negative correlation with increasing temperature, while bacterial gene copy numbers showed a strong negative correlation with salinity. Bacterial community composition correlated only weakly with temperature but strongly with the pore water stable isotopes �18O and �D, as well as with depth. The bacterial community showed substantial spatial variation and an overall dominance of Actinobacteria, Chloroflexi, Firmicutes, Gemmatimonadetes, and Proteobacteria, which are amongst the microbial taxa that were also found to be active in other frozen permafrost environments. We suggest that, millennia after permafrost warming by over 10 �C, microbial community composition and abundance show some indications for proliferation but mainly reflect the sedimentation history and paleoenvironment and not a direct effect through warming.

Description: Author Posting. © The Author(s), 2018. This is the author's version of the work. It is posted here by permission of National Academy of Sciences for personal use, not for redistribution. The definitive version was published in Proceedings of the National Academy of Sciences of the United States of America 115(52), (2018): E12275-E12284. doi: 10.1073/pnas.1805243115.

Description: Diatoms are prominent eukaryotic phytoplankton despite being limited by the micronutrient iron in vast expanses of the ocean. As iron inputs are often sporadic, diatoms have evolved mechanisms such as the ability to store iron that enable them to bloom when iron is resupplied and then persist when low iron levels are reinstated. Two iron storage mechanisms have been previously described: the protein ferritin and vacuolar storage. To investigate the ecological role of these mechanisms among diatoms, iron addition and removal incubations were conducted using natural phytoplankton communities from varying iron environments. We show that among the predominant diatoms, Pseudo-nitzschia were favored by iron removal and displayed unique ferritin expression consistent with a long-term storage function. Meanwhile, Chaetoceros and Thalassiosira gene expression aligned with vacuolar storage mechanisms. Pseudo-nitzschia also showed exceptionally high iron storage under steady-state high and low iron conditions, as well as following iron resupply to iron-limited cells. We propose that bloom-forming diatoms use different iron storage mechanisms and that ferritin utilization may provide an advantage in areas of prolonged iron limitation with pulsed iron inputs. As iron distributions and availability change, this speculated ferritin-linked advantage may result in shifts in diatom community composition that can alter marine ecosystems and biogeochemical cycles.

Description: We thank the captain and crew of the R/V Melville and the CCGS J. P. Tully as well as the participants of the IRNBRU (MV1405) cruise for the California-based data, particularly K. Ellis [University of North Carolina (UNC)], T. Coale (University of California, San Diego), F. Kuzminov (Rutgers), H. McNair [University of California, Santa Barbara (UCSB)], and J. Jones (UCSB). W. Burns (UNC), S. Haines (UNC), and S. Bargu (Louisiana State University) assisted with sample processing and analysis. This work was funded by the National Science Foundation Grants OCE-1334935 (to A.M.), OCE-1334632 (to B.S.T.), OCE-1333929 (to K.T.), OCE-1334387 (to M.A.B.), OCE-1259776 (to K.W.B), and DGE-1650116 (Graduate Research Fellowship to R.H.L).

Description: 2019-06-11

Description: The Global Ocean Data Analysis Project (GLODAP) is a synthesis effort providing regular compilations of surface to bottom ocean biogeochemical data, with an emphasis on seawater inorganic carbon chemistry and related variables determined through chemical analysis of water samples. This update of GLODAPv2, v2.2019, adds data from 116 cruises to the previous version, extending its coverage in time from 2013 to 2017, while also adding some data from prior years. GLODAPv2.2019 includes measurements from more than 1.1 million water samples from the global oceans collected on 840 cruises. The data for the 12 GLODAP core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, CFC-11, CFC-12, CFC-113, and CCl4) have undergone extensive quality control, especially systematic evaluation of bias. The data are available in two formats: (i) as submitted by the data originator but updated to WOCE exchange format and (ii) as a merged data product with adjustments applied to minimize bias. These adjustments were derived by comparing the data from the 116 new cruises with the data from the 724 quality-controlled cruises of the GLODAPv2 data product. They correct for errors related to measurement, calibration, and data handling practices, taking into account any known or likely time trends or variations. The compiled and adjusted data product is believed to be consistent to better than 0.005 in salinity, 1 % in oxygen, 2 % in nitrate, 2 % in silicate, 2 % in phosphate, 4 µmol kg−1 in dissolved inorganic carbon, 4 µmol kg−1 in total alkalinity, 0.01–0.02 in pH, and 5 % in the halogenated transient tracers. The compilation also includes data for several other variables, such as isotopic tracers. These were not subjected to bias comparison or adjustments. The original data, their documentation and DOI codes are available in the Ocean Carbon Data System of NOAA NCEI (https://www.nodc.noaa.gov/ocads/oceans/GLODAPv2_2019/, last access: 17 September 2019). This site also provides access to the merged data product, which is provided as a single global file and as four regional ones – the Arctic, Atlantic, Indian, and Pacific oceans – under https://doi.org/10.25921/xnme-wr20 (Olsen et al., 2019). The product files also include significant ancillary and approximated data. These were obtained by interpolation of, or calculation from, measured data. This paper documents the GLODAPv2.2019 methods and provides a broad overview of the secondary quality control procedures and results.

Description: We present a Lagrangian convective transport scheme developed for global chemistry and transport models, which considers the variable residence time that an air parcel spends in convection. This is particularly important for accurately simulating the tropospheric chemistry of short-lived species, e.g., for determining the time available for heterogeneous chemical processes on the surface of cloud droplets. In current Lagrangian convective transport schemes air parcels are stochastically redistributed within a fixed time step according to estimated probabilities for convective entrainment as well as the altitude of detrainment. We introduce a new scheme that extends this approach by modeling the variable time that an air parcel spends in convection by estimating vertical updraft velocities. Vertical updraft velocities are obtained by combining convective mass fluxes from meteorological analysis data with a parameterization of convective area fraction profiles. We implement two different parameterizations: a parameterization using an observed constant convective area fraction profile and a parameterization that uses randomly drawn profiles to allow for variability. Our scheme is driven by convective mass fluxes and detrainment rates that originate from an external convective parameterization, which can be obtained from meteorological analysis data or from general circulation models. We study the effect of allowing for a variable time that an air parcel spends in convection by performing simulations in which our scheme is implemented into the trajectory module of the ATLAS chemistry and transport model and is driven by the ECMWF ERA-Interim reanalysis data. In particular, we show that the redistribution of air parcels in our scheme conserves the vertical mass distribution and that the scheme is able to reproduce the convective mass fluxes and detrainment rates of ERA-Interim. We further show that the estimated vertical updraft velocities of our scheme are able to reproduce wind profiler measurements performed in Darwin, Australia, for velocities larger than 0.6 m s−1. SO2 is used as an example to show that there is a significant effect on species mixing ratios when modeling the time spent in convective updrafts compared to a redistribution of air parcels in a fixed time step. Furthermore, we perform long-time global trajectory simulations of radon-222 and compare with aircraft measurements of radon activity.

Description: Holocene permafrost from ice wedge polygons in the vicinity of large seabird breeding colonies in the Thule District, NW Greenland, was drilled to explore the relation between permafrost aggradation and seabird presence. The latter is reliant on the presence of the North Water Polynya (NOW) in the northern Baffin Bay. The onset of peat accumulation associated with the arrival of little auks (Alle alle) in a breeding colony at Annikitisoq, north of Cape York, is radiocarbon-dated to 4400 cal BP. A thick-billed murre (Uria lomvia) colony on Appat (Saunders Island) in the mouth of the Wolstenholme Fjord started 5650 cal BP. Both species provide marine-derived nutrients (MDNs) that fertilize vegetation and promote peat growth. The geochemical signature of organic matter left by the birds is traceable in the frozen Holocene peat. The peat accumulation rates at both sites are highest after the onset, decrease over time, and were about 2-times faster at the little auk site than at the thick-billed murre site. High accumulation rates induce shorter periods of organic matter (OM) decomposition before it enters the perennially frozen state. This is seen in comparably high C=N ratios and less depleted 13C, pointing to a lower degree of OM decomposition at the little auk site, while the opposite pattern can be discerned at the thick-billed murre site. Peat accumulation rates correspond to 15N trends, where decreasing accumulation led to increasing depletion in 15N as seen in the little-auk-related data. In contrast, the more decomposed OM of the thick-billed murre site shows almost stable 15N. Late Holocene wedge ice fed by cold season precipitation was studied at the little auk site and provides the first stable-water isotopic record from Greenland with mean 18O of 8:00:8, mean D of 36:25:7, mean d excess of 7:70:7, and a 18O-D slope of 7.27, which is close to those of the modern Thule meteoric water line. The syngenetic ice wedge polygon development is mirrored in testacean records of the little auk site and delineates polygon low-center, dry-out, and polygon-high-center stages. The syngenetic permafrost formation directly depending on peat growth (controlled by bird activity) falls within the period of neoglacial cooling and the establishment of the NOW, thus indirectly following the Holocene climate trends.

Description: © The Author(s), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Proceedings of the National Academy of Sciences.of the United States of America 116(36), (2019): 17666-17672. doi:10.1073/pnas.1907871116.

Description: The conditions of methane (CH4) formation in olivine-hosted secondary fluid inclusions and their prevalence in peridotite and gabbroic rocks from a wide range of geological settings were assessed using confocal Raman spectroscopy, optical and scanning electron microscopy, electron microprobe analysis, and thermodynamic modeling. Detailed examination of 160 samples from ultraslow- to fast-spreading midocean ridges, subduction zones, and ophiolites revealed that hydrogen (H2) and CH4 formation linked to serpentinization within olivine-hosted secondary fluid inclusions is a widespread process. Fluid inclusion contents are dominated by serpentine, brucite, and magnetite, as well as CH4(g) and H2(g) in varying proportions, consistent with serpentinization under strongly reducing, closed-system conditions. Thermodynamic constraints indicate that aqueous fluids entering the upper mantle or lower oceanic crust are trapped in olivine as secondary fluid inclusions at temperatures higher than ∼400 °C. When temperatures decrease below ∼340 °C, serpentinization of olivine lining the walls of the fluid inclusions leads to a near-quantitative consumption of trapped liquid H2O. The generation of molecular H2 through precipitation of Fe(III)-rich daughter minerals results in conditions that are conducive to the reduction of inorganic carbon and the formation of CH4. Once formed, CH4(g) and H2(g) can be stored over geological timescales until extracted by dissolution or fracturing of the olivine host. Fluid inclusions represent a widespread and significant source of abiotic CH4 and H2 in submarine and subaerial vent systems on Earth, and possibly elsewhere in the solar system.

Description: We are indebted to J. Eckert for his support with FE-EMPA; to K. Aquinho and E. Codillo for providing samples from Zambales; to K. Aquinho for Raman analysis of some of the samples from Zambales and Mt. Dent; to H. Dick for providing access to his thin section collection; to the curators of the IODP core repositories for providing access to Ocean Drilling Program (ODP) and Integrated Ocean Drilling Program (IODP) samples; and to the captains and crews of the many cruises without whom the collection of these samples would not have been possible. Reviews by Peter Kelemen and an anonymous referee greatly improved this manuscript. This study is supported with funds provided by the National Science Foundation (NSF-OCE Award 1634032 to F.K. and J.S.S.).

Description: 2020-02-19

Description: Author Posting. © The Author(s), 2019. This is the author's version of the work. It is posted here by permission of National Academy of Sciences for personal use, not for redistribution. The definitive version was published in Proceedings of the National Academy of Sciences of the United States of America 116(20), (2019):9925-9930, doi:10.1073/pnas.1818349116.

Description: Microbial capacity to metabolize arsenic is ancient, arising in response to its pervasive presence in the environment, which was largely in the form of As(III) in the early anoxic ocean. Many biological arsenic transformations are aimed at mitigating toxicity; however, some microorganisms can respire compounds of this redox-sensitive element to reap energetic gains. In several modern anoxic marine systems concentrations of As(V) are higher relative to As(III) than what would be expected from the thermodynamic equilibrium, but the mechanism for this discrepancy has remained unknown. Here we present evidence of a complete respiratory arsenic cycle, consisting of dissimilatory As(V) reduction and chemoautotrophic As(III) oxidation, in the pelagic ocean. We identified the presence of genes encoding both subunits of the respiratory arsenite oxidase AioA and the dissimilatory arsenate reductase ArrA in the Eastern Tropical North Pacific (ETNP) oxygen-deficient zone (ODZ). The presence of the dissimilatory arsenate reductase gene arrA was enriched on large particles (〉30 um), similar to the forward bacterial dsrA gene of sulfate-reducing bacteria, which is involved in the cryptic cycling of sulfur in ODZs. Arsenic respiratory genes were expressed in metatranscriptomic libraries from the ETNP and the Eastern Tropical South Pacific (ETSP) ODZ, indicating arsenotrophy is a metabolic pathway actively utilized in anoxic marine water columns. Together these results suggest arsenic-based metabolisms support organic matter production and impact nitrogen biogeochemical cycling in modern oceans. In early anoxic oceans, especially during periods of high marine arsenic concentrations, they may have played a much larger role.

Description: We thank John Baross and Rika Anderson for helpful discussions and feedback on this project. We also thank the chief scientists of the research cruise, Al Devol and Bess Ward, as well as the captain and crew of the R/V Thomas G. Thompson. This work was supported through a NASA Earth and Space Sciences Graduate Research Fellowship to J.K.S. and National Science Foundation Grant OCE-1138368 (to G.R.).

Description: 2019-10-29

Description: Author Posting. © National Academy of Sciences, 2019. This article is posted here by permission of National Academy of Sciences for personal use, not for redistribution. The definitive version was published in Proceedings of the National Academy of Sciences 116 (24), (2019):11646-11651, doi:10.1073/pnas.1900371116.

Description: Measurements show large decadal variability in the rate of CO2 accumulation in the atmosphere that is not driven by CO2 emissions. The decade of the 1990s experienced enhanced carbon accumulation in the atmosphere relative to emissions, while in the 2000s, the atmospheric growth rate slowed, even though emissions grew rapidly. These variations are driven by natural sources and sinks of CO2 due to the ocean and the terrestrial biosphere. In this study, we compare three independent methods for estimating oceanic CO2 uptake and find that the ocean carbon sink could be responsible for up to 40% of the observed decadal variability in atmospheric CO2 accumulation. Data-based estimates of the ocean carbon sink from pCO2 mapping methods and decadal ocean inverse models generally agree on the magnitude and sign of decadal variability in the ocean CO2 sink at both global and regional scales. Simulations with ocean biogeochemical models confirm that climate variability drove the observed decadal trends in ocean CO2 uptake, but also demonstrate that the sensitivity of ocean CO2 uptake to climate variability may be too weak in models. Furthermore, all estimates point toward coherent decadal variability in the oceanic and terrestrial CO2 sinks, and this variability is not well-matched by current global vegetation models. Reconciling these differences will help to constrain the sensitivity of oceanic and terrestrial CO2 uptake to climate variability and lead to improved climate projections and decadal climate predictions.

Description: We thank Rebecca Wright and Erik Buitenhuis at University of East Anglia, Norwich, for providing updated runs from the NEMO-PlankTOM5 model. T.D. was supported by NSF Grant OCE-1658392. C.L.Q. thanks the UK Natural Environment Research Council for supporting the SONATA Project (Grant NE/P021417/1). P.L. was supported by the Max Planck Society for the Advancement of Science. J.H. was supported under Helmholtz Young Investigator Group Marine Carbon and Ecosystem Feedbacks in the Earth System (MarESys) Grant VH-NG-1301. S.B. and R.S. were supported by the H2020 project CRESCENDO “Coordinated Research in Earth Systems and Climate: Experiments, Knowledge, Dissemination and Outreach,” which received funding from the European Union’s Horizon 2020 research and innovation program under Grant No 641816. SOCAT is an international effort, endorsed by the International Ocean Carbon Coordination Project, the Surface Ocean-Lower Atmosphere Study, and the Integrated Marine Biosphere Research program, to deliver a uniformly quality-controlled surface ocean CO2 database. The many researchers and funding agencies responsible for the collection of data and quality control are thanked for their contributions to SOCAT.

Description: 2019-11-28

Description: © The Author(s), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Proceedings of the National Academy of Sciences.of the United States of America 116(36), (2019): 17934-17942, doi:10.1073/pnas.1910121116.

Description: Plastid endosymbiosis has been a major force in the evolution of eukaryotic cellular complexity, but how endosymbionts are integrated is still poorly understood at a mechanistic level. Dinoflagellates, an ecologically important protist lineage, represent a unique model to study this process because dinoflagellate plastids have repeatedly been reduced, lost, and replaced by new plastids, leading to a spectrum of ages and integration levels. Here we describe deep-transcriptomic analyses of the Antarctic Ross Sea dinoflagellate (RSD), which harbors long-term but temporary kleptoplasts stolen from haptophyte prey, and is closely related to dinoflagellates with fully integrated plastids derived from different haptophytes. In some members of this lineage, called the Kareniaceae, their tertiary haptophyte plastids have crossed a tipping point to stable integration, but RSD has not, and may therefore reveal the order of events leading up to endosymbiotic integration. We show that RSD has retained its ancestral secondary plastid and has partitioned functions between this plastid and the kleptoplast. It has also obtained genes for kleptoplast-targeted proteins via horizontal gene transfer (HGT) that are not derived from the kleptoplast lineage. Importantly, many of these HGTs are also found in the related species with fully integrated plastids, which provides direct evidence that genetic integration preceded organelle fixation. Finally, we find that expression of kleptoplast-targeted genes is unaffected by environmental parameters, unlike prey-encoded homologs, suggesting that kleptoplast-targeted HGTs have adapted to posttranscriptional regulation mechanisms of the host.

Description: We are grateful to Martin Kolisko and Fabien Burki for helpful discussion about and comments on the phylogenetic analysis; and Filip Husnik and Vittorio Boscaro for valuable comments on the manuscript. This work was supported by a grant from the National Science Foundation to R.J.G. and P.J.K. (PLR-1341362) and from the Natural Sciences and Engineering Research Council of Canada to P.J.K. (RGPIN-2014-03994).

Description: 2020-02-19

Description: © The Author(s), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Diaz, J. M., Plummer, S., Hansel, C. M., Andeer, P. F., Saito, M. A., & McIlvin, M. R. NADPH-dependent extracellular superoxide production is vital to photophysiology in the marine diatom Thalassiosira oceanica. Proceedings of the National Academy of Sciences of the United States of America, 116 (33), (2019): 16448-16453, doi: 10.1073/pnas.1821233116.

Description: Reactive oxygen species (ROS) like superoxide drive rapid transformations of carbon and metals in aquatic systems and play dynamic roles in biological health, signaling, and defense across a diversity of cell types. In phytoplankton, however, the ecophysiological role(s) of extracellular superoxide production has remained elusive. Here, the mechanism and function of extracellular superoxide production by the marine diatom Thalassiosira oceanica are described. Extracellular superoxide production in T. oceanica exudates was coupled to the oxidation of NADPH. A putative NADPH-oxidizing flavoenzyme with predicted transmembrane domains and high sequence similarity to glutathione reductase (GR) was implicated in this process. GR was also linked to extracellular superoxide production by whole cells via quenching by the flavoenzyme inhibitor diphenylene iodonium (DPI) and oxidized glutathione, the preferred electron acceptor of GR. Extracellular superoxide production followed a typical photosynthesis-irradiance curve and increased by 30% above the saturation irradiance of photosynthesis, while DPI significantly impaired the efficiency of photosystem II under a wide range of light levels. Together, these results suggest that extracellular superoxide production is a byproduct of a transplasma membrane electron transport system that serves to balance the cellular redox state through the recycling of photosynthetic NADPH. This photoprotective function may be widespread, consistent with the presence of putative homologs to T. oceanica GR in other representative marine phytoplankton and ocean metagenomes. Given predicted climate-driven shifts in global surface ocean light regimes and phytoplankton community-level photoacclimation, these results provide implications for future ocean redox balance, ecological functioning, and coupled biogeochemical transformations of carbon and metals.

Description: This work was supported by a postdoctoral fellowship from the Ford Foundation (to J.M.D.), the National Science Foundation (NSF) under grants OCE 1225801 (to J.M.D.) and OCE 1246174 (to C.M.H.), a Junior Faculty Seed Grant from the University of Georgia Research Foundation (to J.M.D.), and a National Science Foundation Graduate Research Fellowship (to S.P.). The FIRe was purchased through a NSF equipment improvement grant (1624593).The authors thank Melissa Soule for assistance with LC/MS/MS analysis of peptide samples.

Description: Author Posting. © National Academy of Sciences, 2019. This article is posted here by permission of National Academy of Sciences for personal use, not for redistribution. The definitive version was published in Proceedings of the National Academy of Sciences 116(27), (2019): 13233-13238, doi: 10.1073/pnas.1904087116.

Description: The overturning circulation of the global ocean is critically shaped by deep-ocean mixing, which transforms cold waters sinking at high latitudes into warmer, shallower waters. The effectiveness of mixing in driving this transformation is jointly set by two factors: the intensity of turbulence near topography and the rate at which well-mixed boundary waters are exchanged with the stratified ocean interior. Here, we use innovative observations of a major branch of the overturning circulation—an abyssal boundary current in the Southern Ocean—to identify a previously undocumented mixing mechanism, by which deep-ocean waters are efficiently laundered through intensified near-boundary turbulence and boundary–interior exchange. The linchpin of the mechanism is the generation of submesoscale dynamical instabilities by the flow of deep-ocean waters along a steep topographic boundary. As the conditions conducive to this mode of mixing are common to many abyssal boundary currents, our findings highlight an imperative for its representation in models of oceanic overturning.

Description: The DynOPO project is supported by the UK Natural Environment Research Council (grants NE/K013181/1 and NE/K012843/1) and the US National Science Foundation (grants OCE-1536453 and OCE-1536779). A.C.N.G. acknowledges the support of the Royal Society and the Wolfson Foundation. S.L. acknowledges the support of award NA14OAR4320106 from the National Oceanic and Atmospheric Administration, US Department of Commerce. The statements, findings, conclusions, and recommendations are those of the authors, and do not necessarily reflect the views of the National Oceanic and Atmospheric Administration, or the US Department of Commerce. We are grateful to the scientific party, crew, and technicians on the RRS James Clark Ross for their hard work during data collection.

Description: 2019-12-18

Description: Author Posting. © The Author(s), 2019. This is the author's version of the work. It is posted here by permission of National Academy of Sciences for personal use, not for redistribution. The definitive version was published in Proceedings of the National Academy of Sciences of the United States of America 116(25), (2019):12343-12352, doi:10.1073/pnas.1901080116.

Description: Genes encoding cytochrome P450 (CYP; P450) enzymes occur widely in the Archaea, Bacteria, and Eukarya, where they play important roles in metabolism of endogenous regulatory molecules and exogenous chemicals. We now report that genes for multiple and unique P450s occur commonly in giant viruses in the Mimiviridae, Pandoraviridae, and other families in the proposed order Megavirales. P450 genes were also identified in a herpesvirus (Ranid herpesvirus 3) and a phage (Mycobacterium phage Adler). The Adler phage P450 was classified as CYP102L1, and the crystal structure of the open form was solved at 2.5 Å. Genes encoding known redox partners for P450s (cytochrome P450 reductase, ferredoxin and ferredoxin reductase, and flavodoxin and flavodoxin reductase) were not found in any viral genome so far described, implying that host redox partners may drive viral P450 activities. Giant virus P450 proteins share no more than 25% identity with the P450 gene products we identified in Acanthamoeba castellanii, an amoeba host for many giant viruses. Thus, the origin of the unique P450 genes in giant viruses remains unknown. If giant virus P450 genes were acquired from a host, we suggest it could have been from an as yet unknown and possibly ancient host. These studies expand the horizon in the evolution and diversity of the enormously important P450 superfamily. Determining the origin and function of P450s in giant viruses may help to discern the origin of the giant viruses themselves.

Description: We thank Dr. David Nes (Texas Tech University) for providing sterols and Dr. Matthieu Legendre and Dr. Chantal Abergel (CNRS, Marseille) for access to the P. celtis sequences. Drs. Irina Arkhipova, Mark Hahn, Judith Luborsky, and Ann Bucklin commented on the manuscript. The research was supported by a USA-UK Fulbright Scholarship and a Royal Society grant (to D.C.L.), the Boston University Superfund Research Program [NIH Grant 5P42ES007381 (to J.J.S. and J.V.G.) and NIH Grant 5U41HG003345 (to J.V.G.)], the European Regional Development Fund and Welsh Government Project BEACON (S.L.K.), the Woods Hole Center for Oceans and Human Health [NIH Grant P01ES021923 and National Science Foundation Grant OCE-1314642 (to J.J.S.)], and NIH Grant R01GM53753 (to T.L.P.).

Description: 2019-12-05

Description: Author Posting. © National Academy of Sciences, 2019. This article is posted here by permission of National Academy of Sciences for personal use, not for redistribution. The definitive version was published in Proceedings of the National Academy of Sciences 116 (35), (2019): 17187-17192, doi:10.1073/pnas.1903067116.

Description: Mesoscale eddies are critical components of the ocean’s “internal weather” system. Mixing and stirring by eddies exerts significant control on biogeochemical fluxes in the open ocean, and eddies may trap distinctive plankton communities that remain coherent for months and can be transported hundreds to thousands of kilometers. Debate regarding how and why predators use fronts and eddies, for example as a migratory cue, enhanced forage opportunities, or preferred thermal habitat, has been ongoing since the 1950s. The influence of eddies on the behavior of large pelagic fishes, however, remains largely unexplored. Here, we reconstruct movements of a pelagic predator, the blue shark (Prionace glauca), in the Gulf Stream region using electronic tags, earth-observing satellites, and data-assimilating ocean forecasting models. Based on 〉2,000 tracking days and nearly 500,000 high-resolution time series measurements collected by 15 instrumented individuals, we show that blue sharks seek out the interiors of anticyclonic eddies where they dive deep while foraging. Our observations counter the existing paradigm that anticyclonic eddies are unproductive ocean “deserts” and suggest anomalously warm temperatures in these features connect surface-oriented predators to the most abundant fish community on the planet in the mesopelagic. These results also shed light on the ecosystem services provided by mesopelagic prey. Careful consideration will be needed before biomass extraction from the ocean twilight zone to avoid interrupting a key link between planktonic production and top predators. Moreover, robust associations between targeted fish species and oceanographic features increase the prospects for effective dynamic ocean management.

Description: We thank D. McGillicuddy, G. Lawson, and G. Flierl for helpful discussions while developing this work and 2 anonymous reviewers whose feedback significantly improved the manuscript. We also thank C. Fischer and the OCEARCH team for their support of this research. This work was funded by awards to C.D.B. from the Martin Family Society of Fellows for Sustainability Fellowship at the Massachusetts Institute of Technology; the Grassle Fellowship and Ocean Venture Fund at the Woods Hole Oceanographic Institution; and the National Aeronatics and Space Administration (NASA) Earth and Space Science Fellowship. C.D.B. and P.G. acknowledge support from the NASA New Investigator Program Award 80NSSC18K0757, and P.G. acknowledges support from NSF Award OCE-1558809. This research is partially supported by funding to S.R.T. as part of the Audacious Project, a collaborative endeavor, housed at TED. We thank donors to the Woods Hole Oceanographic Institution (WHOI) ProjectWHOI crowdfunding campaign: The Secret Lives of Sharks. Computational support was provided by the Amazon Web Services Cloud Credits for Research program. Funding for the development of HYCOM has been provided by the National Ocean Partnership Program and the Office of Naval Research.

Description: 2020-02-06

Description: Author Posting. © National Academy of Sciences, 2019. This article is posted here by permission of National Academy of Sciences for personal use, not for redistribution. The definitive version was published in Proceedings of the National Academy of Sciences 116 (35), (2019): 17207-17212, doi:10.1073/pnas.1900325116.

Description: It has been hypothesized that the overall size of—or efficiency of carbon export from—the biosphere decreased at the end of the Great Oxidation Event (GOE) (ca. 2,400 to 2,050 Ma). However, the timing, tempo, and trigger for this decrease remain poorly constrained. Here we test this hypothesis by studying the isotope geochemistry of sulfate minerals from the Belcher Group, in subarctic Canada. Using insights from sulfur and barium isotope measurements, combined with radiometric ages from bracketing strata, we infer that the sulfate minerals studied here record ambient sulfate in the immediate aftermath of the GOE (ca. 2,018 Ma). These sulfate minerals captured negative triple-oxygen isotope anomalies as low as ∼ −0.8‰. Such negative values occurring shortly after the GOE require a rapid reduction in primary productivity of 〉80%, although even larger reductions are plausible. Given that these data imply a collapse in primary productivity rather than export efficiency, the trigger for this shift in the Earth system must reflect a change in the availability of nutrients, such as phosphorus. Cumulatively, these data highlight that Earth’s GOE is a tale of feast and famine: A geologically unprecedented reduction in the size of the biosphere occurred across the end-GOE transition.

Description: Olivia M. J. Dagnaud assisted during fieldwork. S. V. Lalonde and E. A. Sperling provided helpful comments on an early version of the manuscript. We thank N. J. Planavsky and an anonymous reviewer for their constructive feedback. M.S.W.H. was supported by an NSERC PGS-D and student research grants from National Geographic, the APS Lewis and Clark Fund, Northern Science Training Program, McGill University Graduate Research Enhancement and Travel Awards, Geological Society of America, Mineralogical Association of Canada, and Stanford University. P.W.C. acknowledges support from the University of Colorado Boulder, the Agouron Institute Geobiology postdoctoral Fellowship program, a Natural Sciences and Engineering Council of Canada Postgraduate Scholarship–Doctoral Program scholarship, and the NSTP. Y.P. was supported by the Strategic Priority Research Program of CAS (XDB26000000). T.J.H. thanks Maureen E. Auro for laboratory assistance and the NSF for supporting isotope research in the NIRVANA Labs.

Description: 2020-02-12

Description: Author Posting. © National Academy of Sciences, 2019. This article is posted here by permission of National Academy of Sciences for personal use, not for redistribution. The definitive version was published in Proceedings of the National Academy of Sciences 116(45), (2019): 22518-22525, doi:10.1073/pnas.1913714116.

Description: The Ganges–Brahmaputra (G-B) River system transports over a billion tons of sediment every year from the Himalayan Mountains to the Bay of Bengal and has built the world’s largest active sedimentary deposit, the Bengal Fan. High sedimentation rates drive exceptional organic matter preservation that represents a long-term sink for atmospheric CO2. While much attention has been paid to organic-rich fine sediments, coarse sediments have generally been overlooked as a locus of organic carbon (OC) burial. However, International Ocean Discovery Program Expedition 354 recently discovered abundant woody debris (millimeter- to centimeter-sized fragments) preserved within the coarse sediment layers of turbidite beds recovered from 6 marine drill sites along a transect across the Bengal Fan (∼8°N, ∼3,700-m water depth) with recovery spanning 19 My. Analysis of bulk wood and lignin finds mostly lowland origins of wood delivered episodically. In the last 5 My, export included C4 plants, implying that coarse woody, lowland export continued after C4 grassland expansion, albeit in reduced amounts. Substantial export of coarse woody debris in the last 1 My included one wood-rich deposit (∼0.05 Ma) that encompassed coniferous wood transported from the headwaters. In coarse layers, we found on average 0.16 weight % OC, which is half the typical biospheric OC content of sediments exported by the modern G-B Rivers. Wood burial estimates are hampered by poor drilling recovery of sands. However, high-magnitude, low-frequency wood export events are shown to be a key mechanism for C burial in turbidites.

Description: This work was funded by National Science Foundation Grants OCE-1401217 and COL-T354A55 to S.J.F. and OCE-1400805 to V.G. Graduate student participation in the project received support from University of Southern California Provost’s Fellowship to H.L. Samples were provided by the International Ocean Discovery Program. We are grateful for the efforts of the Expedition 354 Science Party, Carl Johnson, and Zongguang Liu. C.F.-L. and A.G. were supported by IODP-France. We thank Colin Osborne and Maria Vorontsova for helpful discussions.

Description: 2020-04-21