Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Oxygen isotopes and Sr/Ca ratios of fossil corals from Bonaire

Explanations for the demise of the Classic Maya civilization on the Yucatán Peninsula during the Terminal Classic Period (TCP; CE 750-1050) are controversial. Multiyear droughts are one likely cause, but the role of the Caribbean Sea, the dominant moisture source for Mesoamerica, remains largely unknown. Here we present bimonthly resolved snapshots of reconstructed sea surface temperature (SST) and salinity (SSS) variability in the southern Caribbean from precisely dated fossil corals. Our fossil coral results from Bonaire indicate strong interannual to decadal SST and SSS variability in the southern Caribbean Sea during the TCP with multiyear extremes of high SSS and high SST that coincide with droughts on the Yucatán Peninsula. The results are best explained by changed Caribbean SST gradients affecting the Caribbean low-level atmospheric jet with consequences for Mesoamerican precipitation, possibly linked to changes in Atlantic Meridional Overturning Circulation strength. Our findings provide a new perspective on the anomalous hydrological changes on the Yucatán Peninsula during the TCP that complement the often-suggested southward displacement of the Intertropical Convergence Zone. We advocate for a strong role of ocean-atmosphere interactions in the Caribbean Sea related to the multiyear variations in Caribbean Sea surface conditions as an important driver of the spatially complex pattern of hydrological anomalies during the TCP.

Temperature proxies of sediment cores from the Sea of Marmara

Reconstructing ocean temperature values is a major target in paleoceanography and climate research. However, most temperature proxies are organism-based and thus suffer from an "ecological bias". Multiproxy approaches can potentially overcome this bias but typically require more investment in time and resources, while being susceptible to errors induced by sample preparation steps necessary before analysis. Three lipid-based temperature proxies are widely used: UK'37 (based on long chain alkenones from phytoplanktonic haptophytes), TEX86 [based on glycerol dialkyl glycerol tetraethers (GDGTs) from pelagic archaea] and LDI (based on long chain diols from phytoplanktonic eustigmatophytes). So far, separate analytical methods, including gas chromatography (GC) and liquid chromatography (LC), have been used to determine these proxies. Here we present a sensitive method for determining all three in a single normal phase high performance LC-atmospheric pressure chemical ionization mass spectrometry (NP-HPLC-APCI-MS) analysis. Each of the long chain alkenones and long chain diols was separated and unambiguously identified from the accurate masses and characteristic fragmentation during multiple stage MS analysis (MS2). Comparison of conventional GC and HPLC-MS methods showed similar results for UK'37 and LDI, respectively, using diverse environmental samples and an Emiliania huxleyi culture. Including the three sea surface temperature (SST) proxies; the NP-HPLC-APCI-MS method in fact allows simultaneous determination of nine paleoenvironmental proxies. The extent to which the ecology of the source organisms (ecological bias) influences lipid composition and thereby the reconstructed temperature values was demonstrated by applying the new method to a sediment core from the Sea of Marmara, covering the last 21 kyr BP. Reconstructed SST values differed considerably between the proxies for the Last Glacial Maximum (LGM) and the period of Sapropel S1 formation at ca. 10 kyr BP, whereas the trends during the late Holocene were similar. Changes in the composition of alkenone-producing species at the transition from the LGM to the Bølling/Allerød (B/A) were inferred from unreasonably high UK'37-derived SST values (ca. 20 °C) during the LGM. We ascribe discrepancies between the reconstructed temperature records during S1 deposition to habitat change, e.g. a different depth due to changes in nutrient availability.

Mid-Pleistocene planktonic foraminiferal composition and IRD content in ODP Sites 162-980 and 162-984

We investigated surface and deep ocean variability in the subpolar North Atlantic from 1000 to 500 thousand years ago (ka) based on two Ocean Drilling Program (ODP) sites, Feni drift site 980 (55°29'N, 14°42'W) and Bjorn drift site 984 (61°25'N, 24°04'W). Benthic foraminiferal stable isotope data, planktic foraminiferal faunas, ice-rafted debris data, and faunally based sea-surface temperature estimates help test the hypothesis that oceanographic changes in the North Atlantic region were associated with the onset of the 100-kyr world during the mid-Pleistocene revolution. Based on percentage of Neogloboquadrina pachyderma (s) records from both sites, surface waters during interglacials and glacials were cooler in the mid-Pleistocene than during marine isotope stages (MIS) 5 and 6. In particular, interglaciations at Bjorn drift site 984 were significantly cooler. Faunal evidence suggests that the interglacial Arctic front shifted from a position between the two sites to a position northwest of Bjorn drift site 984 after ca. 610 ka. As during the late Pleistocene, we find faunal evidence for lagging surface warmth at most of the glacial initiations during the mid-Pleistocene. Each initiation is associated with high benthic d13C values that are maintained into the succeeding glaciation, which we term "lagging NADW production." These findings indicate that lagging warmth and lagging NADW production are robust features of the regional climate system that persist in the middle to late Pleistocene.

Stable carbon isotopes, Methane Index, and distributions of different GDGTs

Gas hydrates represent one of the largest pools of readily exchangeable carbon on Earth's surface. Releases of the greenhouse gas methane from hydrates are proposed to be responsible for climate change at numerous events in geological history. Many of these inferred events, however, were based on carbonate carbon isotopes which are susceptible to diagenetic alterations. Here we propose a molecular fossil proxy, i.e., the "Methane Index (MI)", to detect and document the destabilization and dissociation of marine gas hydrates. MI consists of the relative distribution of glycerol dibiphytanyl glycerol tetraethers (GDGTs), the core membrane lipids of archaea. The rational behind MI is that in hydrate-impacted environments, the pool of archaeal tetraether lipids is dominated by GDGT-1, -2 and -3 due to the large contribution of signals from the methanotrophic archaeal community. Our study in the Gulf of Mexico cold-seep sediments demonstrates a correlation between MI and the compound-specific carbon isotope of GDGTs, which is strong evidence supporting the MI-methane consumption relationship. Preliminary applications of MI in a number of hydrate-impacted and/or methane-rich environments show diagnostic MI values, corroborating the idea that MI may serve as a robust indicator for hydrate dissociation that is useful for studies of global carbon cycling and paleoclimate change.

Calibration and application of B/Ca, Cd/Ca, and d11B in Neogloboquadrina pachyderma (sinistral)

The North Atlantic and Norwegian Sea are prominent sinks of atmospheric CO2 today, but their roles in the past remain poorly constrained. In this study, we attempt to use B/Ca and d11B ratios in the planktonic foraminifera Neogloboquadrina pachyderma (sinistral variety) to reconstruct subsurface water pH and pCO2 changes in the polar North Atlantic during the last deglaciation. Comparison of core-top results with nearby hydrographic data shows that B/Ca in N. pachyderma (s) is mainly controlled by seawater [B(OH)4]?/[HCO3]? with a roughly constant partition coefficient (KD =([B/Ca]of CaCO3)/([B(OH)4]-/[HCO3]-)of seawater) of 1.48 ± 0.15 * 10**-3 (2sigma), and d11B in this species is offset below d11B of the borate in seawater by 3.38 ± 0.71 per mil (2sigma). These values represent our best estimates with the sparse available hydrographic data close to our core-tops. More culturing and sediment trap work is needed to improve our understanding of boron incorporation into N. pachyderma (s). Application of a constant KD of 1.48 * 10**-3 to high resolution N. pachyderma (s) B/Ca records from two adjacent cores off Iceland shows that subsurface pCO2 at the habitat depth of N. pachyderma (s) (~50 m) generally followed the atmospheric CO2 trend but with negative offsets of ~10-50 ppmv during 19-10 ka. These B/Ca-based reconstructions are supported by independent estimates from low-resolution d11B measurements in the same cores. We also calibrate and apply Cd/Ca in N. pachyderma (s) to reconstruct nutrient levels for the same down cores. Like today's North Atlantic, past subsurface pCO2 variability off Iceland was significantly correlated with nutrient changes that might be linked to surface nutrient utilization and mixing within the upper water column. Because surface pCO2 (at 0 m water depth) is always lower than at deeper depths and if the application of a constant KD is valid, our results suggest that the polar North Atlantic has remained a CO2 sink during the calcification seasons of N. pachyderma (s) over the last deglaciation.

Sapropel cycles in the Mediterranean Sea

Water column stratification increased at climatic transitions from cold to warm periods during the late Quaternary and led to anoxic conditions and sapropel formation in the deep eastern Mediterranean basins. High-resolution data sets on sea-surface temperatures (SST) (estimated from UK'37 indices) and d18O of planktonic foraminifer calcite (d18Ofc) across late Pleistocene sapropel intervals show that d18Ofc decreased (between 1 and 4.6 per mil) and SST increased (between 0.7° and 6.7°C). Maximal d18Oseawater depletion of eastern Mediterranean surface waters at the transition is between 0.5 and 3.0 per mil, and in all but one case exceeded the depletion seen in a western Mediterranean core. The depletion in d18Oseawater is most pronounced at sapropel bases, in agreement with an initial sudden input of monsoon-derived freshwater. Most sapropels coincide with warming trends of SST. The density decrease by initial freshwater input and continued warming of the sea surface pooled fresh water in the surface layer and prohibited deep convection down to ageing deep water emplaced during cold and arid glacial conditions. An exception to this pattern is "glacial" sapropel S6; its largest d18Oseawater depletion (3 per mil) is almost matched by the depletion in the western Mediterranean Sea, and it is accompanied by surface water cooling following an initially rapid warming phase. A second period of significant isotopic depletion is in isotope stage 6 at the 150 kyr insolation maximum. While not expressed as a sapropel due to cold SST, it is in accord with a strengthened monsoon in the southern catchment.

Stable isotope ratios of foraminifera and sea surface temperature estimation of sediment core W2804A-14

Carbon isotopically based estimates of CO2 levels have been generated from a record of the photosynthetic fractionation of 13C (epsilon p) in a central equatorial Pacific sediment core that spans the last ~255 ka. Contents of 13C in phytoplanktonic biomass were determined by analysis of C37 alkadienones. These compounds are exclusive products of Prymnesiophyte algae which at present grow most abundantly at depths of 70-90 m in the central equatorial Pacific. A record of the isotopic compostion of dissolved CO2 was constructed from isotopic analyses of the planktonic foraminifera Neogloboquadrina dutertrei, which calcifies at 70-90 m in the same region. Values of epsilon p, derived by comparison of the organic and inorganic delta values, were transformed to yield concentrations of dissolved CO2 (c e) based on a new, site-specific calibration of the relationship between epsilon p and c e. The calibration was based on reassessment of existing epsilon p versus c e data, which support a physiologically based model in which epsilon p is inversely related to c e. Values of PCO2, the partial pressure of CO2 that would be in equilibrium with the estimated concentrations of dissolved CO2, were calculated using Henry's law and the temperature determined from the alkenone-unsaturation index UK 37. Uncertainties in these values arise mainly from uncertainties about the appropriateness (particularly over time) of the site-specific relationship between epsilon p and 1/c e. These are discussed in detail and it is concluded that the observed record of epsilon p most probably reflects significant variations in Delta pCO2, the ocean-atmosphere disequilibrium, which appears to have ranged from ~110 µatm during glacial intervals (ocean > atmosphere) to ~60 µatm during interglacials. Fluxes of CO2 to the atmosphere would thus have been significantly larger during glacial intervals. If this were characteristic of large areas of the equatorial Pacific, then greater glacial sinks for the equatorially evaded CO2 must have existed elsewhere. Statistical analysis of air-sea pCO2 differences and other parameters revealed significant (p < 0.01) inverse correlations of Delta pCO2 with sea surface temperature and with the mass accumulation rate of opal. The former suggests response to the strength of upwelling, the latter may indicate either drawdown of CO2 by siliceous phytoplankton or variation of [CO2]/[Si(OH)4] ratios in upwelling waters.

Sable oxygen isotope and Mg/Ca ratios of planktonic foraminifera from DSDP Hole 31-292

This study presents new evidence of when and how the Western Pacific Warm Pool (WPWP) was established in its present form. We analyzed planktic foraminifera, oxygen isotopes, and Mg/Ca ratios in upper Miocene through Pleistocene sediments collected at Deep Sea Drilling Program (DSDP) Site 292. These data were then compared with those reported from Ocean Drilling Program (ODP) Site 806. Both drilling sites are located in the western Pacific Ocean. DSDP Site 292 is located in the northern margin of the modern WPWP and ODP Site 806 near the center of the WPWP. Three stages of development in surface-water conditions are identified in the region using planktic foraminferal data. During the initial stage, from 8.5 to 4.4 Ma, Site 806 was overlain by warm surface water but Site 292 was not, as indicated by the differences in faunal compositions and sea-surface temperature (SST) between the two sites. In addition, the vertical thermal gradient at Site 292 was weak during this period, as indicated by the small differences in the delta18O values between Globigerinoides sacculifer and Pulleniatina spp. During stage two, from 4.4 to 3.6 Ma, the SST at Site 292 rapidly increased to 27 °C, but the vertical thermal gradient had not yet be strengthened, as shown by Mg/Ca ratios and the presence of both mixed-layer dwellers and thermocline dwellers. Finally, a warm mixed layer with a high SST ca. 28 °C and a strong vertical thermal gradient were established at Site 292 by 3.6 Ma. This event is marked by the dominance of mixed-layer dwellers, a high and stable SST, and a larger differences in the delta18O values between G. sacculifer and Pulleniatina spp. Thus, evidence of surface-water evolution in the western Pacific suggests that Site 292 came under the influence of the WPWP at 3.6 Ma. The northward expansion of the WPWP from 4.4 to 3.6 Ma and the establishment of the modern WPWP by 3.6 Ma appear to be closely related to the closure of the Indonesian and Central American seaways.

Distribution of deep-sea foraminifera in surface sediments east of New Zealand

Paleobathymetric assessments of fossil foraminiferal faunas play a significant role in the analysis of the paleogeographic, sedimentary, and tectonic histories of New Zealand's Neogene marine sedimentary basins. At depths >100 m, these assessments often have large uncertainties. This study, aimed at improving the precision of paleodepth assessments, documents the present-day distribution of deep-sea foraminifera (>63 µm) in 66 samples of seafloor sediment at 90-700 m water depth (outer shelf to mid-abyssal), east of New Zealand. One hundred and thirty-nine of the 465 recorded species of benthic foraminifera are new records for the New Zealand region. Characters of the foraminiferal faunas which appear to provide the most useful information for estimating paleobathymetry are, in decreasing order of reliability: relative abundance of common benthic species; benthic species associations; upper depth limits of key benthic species; and relative abundance of planktic foraminifera. R mode cluster analysis on the quantitative census data of the 58 most abundant species of benthic foraminifera produced six species associations within three higher level clusters: (1) calcareous species most abundant at mid-bathyal to outer shelf depths (<1000 m); (2) calcareous species most abundant at mid-bathyal and greater depths (>600 m); (3) agglutinated species mostly occurring at deep abyssal depths (>3000 m). A detrended correspondence analysis ordination plot exhibits a strong relationship between these species associations and bathymetry. This is manifest in the bathymetric ranges of the relative abundance peaks of many of the common benthic species (e.g., Abditodentrix pseudothalmanni 500-2800 m, Bolivina robusta 200-650 m, Bulimina marginata f. marginata 20-600 m, B. marginata f. aculeata 400-3000 m, Cassidulina norvangi 1000-4500 m, Epistominella exigua 1000-4700 m, and Trifarina angulosa 10-650 m), which should prove useful in paleobathymetric estimates. The upper depth limits of 28 benthic foraminiferal species (e.g., Fursenkoina complanata 200 m, Bulimina truncana 450 m, Melonis affinis 550 m, Eggerella bradyi 750 m, and Cassidulina norvangi 1000 m) have potential to improve the precision of paleobathymetric estimates based initially on the total faunal composition. The planktic percentage of foraminiferal tests increases from outer shelf to upper abyssal depths followed by a rapid decline within the foraminiferal lysocline (below c. 3600 m). A planktic percentage <50% is suggestive of shelf depths, and >50% is suggestive of bathyal or abyssal depths above the CCD. In the abyssal zone there is dramatic taphonomic loss of most agglutinated tests (except some textulariids) at burial depths of 0.1-0.2 m, which negates the potential usefulness of these taxa in paleobathymetric assessments.

Dinoflagellate cyst assemblages, Mg/Ca based reconstructed temperatures, stable carbon isotopes, and geochemical parameters of sediment core GeoB9508-5

In order to investigate a possible connection between tropical northeast (NE) Atlantic primary productivity, Atlantic meridional overturning circulation (AMOC), and drought in the Sahel region during Heinrich Stadial 1 (HS1), we used dinoflagellate cyst (dinocyst) assemblages, Mg/Ca based reconstructed temperatures, stable carbon isotopes (d13C) and geochemical parameters of a marine sediment core (GeoB 9508-5) from the continental slope offshore Senegal. Our results show a two-phase productivity pattern within HS1 that progressed from an interval of low marine productivity between ~ 19 and 16 kyr BP to a phase with an abrupt and large productivity increase from ~ 16 to 15 kyr BP. The second phase is characterized by distinct heavy planktonic d13C values and high concentrations of heterotrophic dinocysts in addition to a significant cooling signal based on reconstructions of past sea surface temperatures (SST). We conclude that productivity variations within HS1 can be attributed to a substantial shift of West African atmospheric processes. Taken together our results indicate a significant intensification of the North East (NE) trade winds over West Africa leading to more intense upwelling during the last millennium of HS1 between ~ 16 and 15 kyr BP, thus leaving a strong imprint on the dinocyst assemblages and sea surface conditions. Therefore, the two-phase productivity pattern indicates a complex hydrographic setting suggesting that HS1 cannot be regarded as uniform as previously thought.

Sea surafe temperture calculation for the Arabian Sea

Sea surface temperature (SST) and seawater d18O (d18Ow) were reconstructed in a suite of sediment cores from throughout the Arabian Sea for four distinct time intervals (0 ka, 8 ka, 15 ka, and 20 ka) with the aim of understanding the history of the Indian Monsoon and the climate of the Arabian Sea region. This was accomplished through the use of paired Mg/Ca and d18O measurements of the planktonic foraminifer Globigerinoides ruber. By analyzing basin-wide changes and changes in cross-basinal gradients, we assess both monsoonal and regional-scale climate changes. SST was colder than present for the majority of sites within all three paleotime slices. Furthermore, both the Indian Monsoon and the regional Arabian Sea mean climate have varied substantially over the past 20 kyr. The 20 ka and 15 ka time slices exhibit average negative temperature anomalies of 2.5°-3.5°C attributable, in part, to the influences of glacial atmospheric CO2 concentrations and large continental ice sheets. The elimination of the cross-basinal SST gradient during these two time slices likely reflects a decrease in summer monsoon and an increase in winter monsoon strength. Changes in d18Ow that are smaller than the d18O signal due to global ice volume reflect decreased evaporation and increased winter monsoon mixing. SSTs throughout the Arabian Sea were still cooler than present by an average of 1.4°C in the 8 ka time slice. These cool SSTs, along with lower d18Ow throughout the basin, are attributed to stronger than modern summer and winter monsoons and increased runoff and precipitation. The results of this study underscore the importance of taking a spatial approach to the reconstruction of processes such as monsoon upwelling.

40-year long monthly-resolved Sr/Ca record from 2.35 ka Bonaire coral, sample BON-6-A

We present a 40-year long monthly resolved Sr/Ca record from a fossil Diploria strigosa coral from Bonaire (Southern Caribbean Sea) dated with U/Th at 2.35 ka before present (BP). Secondary modifiers of this sea surface temperature (SST) proxy in annually-banded corals such as diagenetic alteration of the skeleton and skeletal growth-rate are investigated. Extensive diagenetic investigations reveal that this fossil coral skeleton is pristine which is further supported by clear annual cycles in the coral Sr/Ca record. No significant correlation between annual growth rate and Sr/Ca is observed, suggesting that the Sr/Ca record is not affected by coral growth. Therefore, we conclude that the observed interannual Sr/Ca variability was influenced by ambient SST variability. Spectral analysis of the annual mean Sr/Ca record reveals a dominant frequency centred at 6-7 years that is not associated with changes of the annual growth rate. The first monthly resolved coral Sr/Ca record from the Southern Caribbean Sea for preindustrial time suggests that fossil corals from Bonaire are suitable tools for reconstructing past SST variability. Coastal deposits on Bonaire provide abundant fossil D. strigosa colonies of Holocene age that can be accurately dated and used to reconstruct climate variability. Comparisons of long monthly resolved Sr/Ca records from multiple fossil corals will provide a mean to estimate seasonality and interannual to interdecadal SST variability of the Southern Caribbean Sea during the Holocene.

Benthic foraminiferal assemblages of middle to late Pliocene sediments in the North Atlantic Ocean

Records of benthic foraminifera from North Atlantic DSDP Site 607 and Hole 610A indicate changes in deep water conditions through the middle to late Pliocene (3.15 to 2.85 Ma). Quantitative analyses of modem associations in the North Atlantic indicate that seven species, Fontbotia wuellerstorfi, Cibicidoides kullenbergi, Uvigerina peregrina, Nuttallides umboniferus, Melonis pompilioides, Globocassidulina subglobosa and Epistominella exigua are useful for paleoenvironmental interpretation. The western North Atlantic basin (Site 607) was occupied by North Atlantic Deep Water (NADW) until c. 2.88 Ma. At that time, N. umboniferus increased, indicating an influx of Southern Ocean Water (SOW). The eastern North Atlantic basin (Hole 610A) was occupied by a relatively warm water mass, possibly Northeastern Atlantic Deep Water (NEADW), through c. 2.94 Ma when SOW more strongly influenced the site. These interpretations are consistent with benthic delta18O and delta13C records from 607 and 610A (Raymo et al., 1992). The results presented in this paper suggest that the North Atlantic was strongly influenced by northern component deep water circulation until 2.90-2.95 Ma. After that there was a transition toward a glacially driven North Atlantic circulation more strongly influenced by SOW associated with the onset of Northern Hemisphere glaciation. The circulation change follows the last significant SST and atmospheric warming prior to c. 2.6 Ma.

Planktic foraminiferal abundance and stable oxygen isotope ratios of sediment cores from the Southern California Borderlands

A variety of evidence suggests that average sea surface temperatures (SSTs) during the last glacial maximum in the California Borderlands region were significantly colder than during the Holocene. Planktonic foraminiferal delta18O evidence and average SST estimates derived by the modern analog technique indicate that temperatures were 6°-10°C cooler during the last glacial relative to the present. The glacial plankton assemblage is dominated by the planktonic foraminifer Neogloboquadrina pachyderma (sinistral coiling) and the coccolith Coccolithus pelagicus, both of which are currently restricted to subpolar regions of the North Pacific. The glacial-interglacial average SST change determined in this study is considerably larger than the 2°C change estimated by Climate: Long-Range Investigation, Mapping, and Prediction (CLIMAP) [1981]. We propose that a strengthened California Current flow was associated with the advance of subpolar surface waters into the Borderlands region during the last glacial.