skip to main content


Title: Global Core Top Calibration of δ 18 O in Planktic Foraminifera to Sea Surface Temperature
Abstract

The oxygen isotopic composition of planktic foraminiferal calcite () is one of the most prevalent proxies used in the paleoceanographic community. The relationship between, temperature, and seawater oxygen isotopic composition () is firmly rooted in thermodynamics, and experimental constraints are commonly used for sea surface temperature (SST) reconstructions. However, in marine sedimentary applications, additional sources of uncertainty emerge, and these uncertainty constraints have not as of yet been included in global calibration models. Here, we compile a global data set of over 2,600 marine sediment core top samples for five planktic species:Globigerinoides ruber,Trilobatus sacculifer,Globigerina bulloides,Neogloboquadrina incompta, andNeogloboquadrina pachyderma. We developed a suite of Bayesian regression models to calibrate the relationship betweenand SST. Spanning SSTs from 0.0 to 29.5 °C, our annual model with species pooled together has a mean standard error of approximately 0.54‰. Accounting for seasonality and species‐specific differences improves model validation, reducing the mean standard error to 0.47‰. Example applications spanning the Late Quaternary show good agreement with independent alkenone‐based estimates. Our pooled calibration model may also be used for reconstruction in the deeper geological past, using modern planktic foraminifera as an analog for non‐extant species. Our core top‐based models provide a robust assessment of uncertainty in thepaleothermometer that can be used in statistical assessments of interproxy and model‐proxy comparisons. The suite of models is publicly available as the Open Source software librarybayfox, for Python, R, and MATLAB/Octave.

 
more » « less
NSF-PAR ID:
10459864
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Paleoceanography and Paleoclimatology
Volume:
34
Issue:
8
ISSN:
2572-4517
Page Range / eLocation ID:
p. 1292-1315
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Flows into and out of the Gulf of Mexico (GoM) are integral to North Atlantic Ocean circulation and help facilitate poleward heat transport in the Western Hemisphere. The GoM also serves as a key source of moisture for most of North America. Modern patterns of sea‐surface temperature (SST) and salinity in the GoM are influenced by the Loop Current, its eddy‐shedding dynamics, and the ensuing interplay with coastal processes. Here, we present sub‐centennial‐scale records of SST and stable oxygen isotope composition of seawater (18Osw; a proxy for salinity) over the past 11,700 years using planktic foraminiferal geochemistry in sediments from the Garrison Basin, northwestern GoM. We measuredO and magnesium‐to‐calcium ratios in tests ofGlobigerinoides ruber(white) to generate quantitative estimates of past sea‐surface conditions. Our results replicate and extend late Holocene reconstructions from the Garrison Basin, using which we then create composites of SST and18Osw. We find considerable centennial and millennial‐scale variability in both SST and18Osw, although their evolution over the Holocene is distinct. Whereas mean‐annual SSTs display pronounced millennial‐scale variability,18Oswexhibits a secular trend spanning multiple millennia and points to increasing northwestern GoM surface salinity since the early Holocene. We then synthesize the available Holocene records from across the GoM and alongside the Garrison Basin composite uncover substantial, yet regionally consistent, spatiotemporal variability. Finally, we discuss the role of the Loop Current and coastal influx of freshwater in imposing these heterogeneities. We conclude that dynamic surface‐ocean changes occurred across the GoM over the Holocene.

     
    more » « less
  2. Abstract

    The genetic effective population size,Ne, can be estimated from the average gametic disequilibrium () between pairs of loci, but such estimates require evaluation of assumptions and currently have few methods to estimate confidence intervals.speed‐neis a suite ofmatlabcomputer code functions to estimatefromwith a graphical user interface and a rich set of outputs that aid in understanding data patterns and comparing multiple estimators.speed‐neincludes functions to either generate or input simulated genotype data to facilitate comparative studies ofestimators under various population genetic scenarios.speed‐newas validated with data simulated under both time‐forward and time‐backward coalescent models of genetic drift. Three classes of estimators were compared with simulated data to examine several general questions: what are the impacts of microsatellite null alleles on,how should missing data be treated, and does disequilibrium contributed by reduced recombination among some loci in a sample impact. Estimators differed greatly in precision in the scenarios examined, and a widely employedestimator exhibited the largest variances among replicate data sets.speed‐neimplements several jackknife approaches to estimate confidence intervals, and simulated data showed that jackknifing over loci and jackknifing over individuals provided ~95% confidence interval coverage for some estimators and should be useful for empirical studies.speed‐neprovides an open‐source extensible tool for estimation offrom empirical genotype data and to conduct simulations of both microsatellite and single nucleotide polymorphism (SNP) data types to develop expectations and to compareestimators.

     
    more » « less
  3. Summary

    Chemotrophic microorganisms gain energy for cellular functions by catalyzing oxidation–reduction (redox) reactions that are out of equilibrium. Calculations of the Gibbs energy (ΔGr) can identify whether a reaction is thermodynamically favourable and quantify the accompanying energy yield at the temperature, pressure and chemical composition in the system of interest. Based on carefully calculated values ofΔGr, we predict a novel microbial metabolism – sulfur comproportionation (3H2S ++ 2H+4S0+ 4H2O). We show that at elevated concentrations of sulfide and sulfate in acidic environments over a broad temperature range, this putative metabolism can be exergonic (ΔGr<0), yielding ~30–50 kJ mol−1. We suggest that this may be sufficient energy to support a chemolithotrophic metabolism currently missing from the literature. Other versions of this metabolism, comproportionation to thiosulfate (H2S ++ H2O) and to sulfite (H2S + 34+ 2H+), are only moderately exergonic or endergonic even at ideal geochemical conditions. Natural and impacted environments, including sulfidic karst systems, shallow‐sea hydrothermal vents, sites of acid mine drainage, and acid–sulfate crater lakes, may be ideal hunting grounds for finding microbial sulfur comproportionators.

     
    more » « less
  4. Abstract

    The mineral apatite, Ca5(PO4)3(F,Cl,OH), is a ubiquitous accessory mineral, with its volatile content and isotopic compositions used to interpret the evolution of H2O on planetary bodies. During hypervelocity impact, extreme pressures shock target rocks resulting in deformation of minerals; however, relatively few microstructural studies of apatite have been undertaken. Given its widespread distribution in the solar system, it is important to understand how apatite responds to progressive shock metamorphism. Here, we present detailed microstructural analyses of shock deformation in ~560 apatite grains throughout ~550 m of shocked granitoid rock from the peak ring of the Chicxulub impact structure, Mexico. A combination of high‐resolution backscattered electron (BSE) imaging, electron backscatter diffraction mapping, transmission Kikuchi diffraction mapping, and transmission electron microscopy is used to characterize deformation within apatite grains. Systematic, crystallographically controlled deformation bands are present within apatite, consistent with tilt boundaries that contain the <c> (axis) and result from slip in <> (direction) on(plane) during shock deformation. Deformation bands contain complex subgrain domains, isolated dislocations, and low‐angle boundaries of ~1° to 2°. Planar fractures within apatite form conjugate sets that are oriented within either {, {, {, or. Complementary electron microprobe analyses (EPMA) of a subset of recrystallized and partially recrystallized apatite grains show that there is an apparent change in MgO content in shock‐recrystallized apatite compositions. This study shows that the response of apatite to shock deformation can be highly variable, and that application of a combined microstructural and chemical analysis workflow can reveal complex deformation histories in apatite grains, some of which result in changes to crystal structure and composition, which are important for understanding the genesis of apatite in both terrestrial and extraterrestrial environments.

     
    more » « less
  5. Abstract

    Debate continues on the amount and distribution of radioactive heat producing elements (i.e., U, Th, and K) in the Earth, with estimates for mantle heat production varying by an order of magnitude. Constraints on the bulk‐silicate Earth's (BSE) radiogenic power also places constraints on overall BSE composition. Geoneutrino detection is a direct measure of the Earth's decay rate of Th and U. The geoneutrino signal has contributions from the local (40%) and global (35%) continental lithosphere and the underlying inaccessible mantle (25%). Geophysical models are combined with geochemical data sets to predict the geoneutrino signal at current and future geoneutrino detectors. We propagated uncertainties, both chemical and physical, through Monte Carlo methods. Estimated total signal uncertainties are on the order of20%, proportionally with geophysical and geochemical inputs contributing30% and70%, respectively. We find that estimated signals, calculated using CRUST2.0, CRUST1.0, and LITHO1.0, are within physical uncertainty of each other, suggesting that the choice of underlying geophysical model will not change results significantly, but will shift the central value by up to15%. Similarly, we see no significant difference between calculated layer abundances and bulk crustal heat production when using these geophysical models. The bulk crustal heat production is calculated as 7  2 TW, which includes an increase of 1 TW in uncertainty relative to previous studies. Combination of our predicted lithospheric signal with measured signals yield an estimated BSE heat production of 21.5  10.4 TW. Future improvements, including uncertainty attribution and near‐field modeling, are discussed.

     
    more » « less