Millenial-Scale stable oxygen isotope record of ODP Site 172-1058

We use the oxygen isotopic composition of planktonic foraminifera Globigerinoides ruber (white) from Ocean Drilling Program Site 1058 in the subtropical northwestern Atlantic to construct a high-resolution (~800 year) climate record spanning the mid-Pleistocene climate transition (~410 ka to 1350 ka). We investigate whether or not millennial-scale instabilities in the proxy record are associated with the extent of continental glaciation. G. ruber d18O values display high-frequency fluctuations throughout the record, but the amplitude about mean glacial and interglacial d18O values increases at marine isotope stage (MIS) 22 (880 ka) and is highest during MIS 12. These observations support that millennial-scale climate instabilities are associated with ice sheet size. Time series analysis illustrates that these variations have significant concentration of spectral power centered on periods of ~10-12 ka and ~5 ka. The timing of these fluctuations agrees well, or coincides with, the periodicities of the second and fourth harmonics, respectively, of precessional forcing at the equator. An insolation-based origin of the millennial-scale instabilities would be independent of ice volume and explains the presence of these fluctuations before the mid-Pleistocene climate transition as well as during interglacial intervals (e.g., MIS 37 and 17). Because the amplitude of the millennial-scale variations increases during the mid-Pleistocene transition, feedback mechanisms associated with the growth of large, 100-ka-paced, polar ice sheets may be important amplifiers of regional surface water hydrographic changes.

Micropaleontological and geochemical proxy studies of sediment core GeoB7926-2

[1] The low-latitude upwelling regime off the Mauritanian coast in the subtropical NE Atlantic accounts for a significant part of global export production. Although productivity variations in coastal upwelling areas are usually attributed to changes in wind stress and upwelling intensity, productivity dynamics off Mauritania are less straightforward because of the complex atmospheric and hydrographic setting. Here we integrate micropaleontological (diatoms) and geochemical (bulk biogenic sediment components, X-ray fluorescence, and alkenones) proxies to examine on submillennial-to-millennial changes in diatom production that occurred off Mauritania, NW Africa, for the last 25 ka. During the Last Glacial Maximum (LGM, 19.0-23.0 ka B.P.), moderate silicate content of upwelled waters coupled with weakened NE trade winds determined moderate diatom productivity. No significant cooling is observed during the LGM, suggesting that our alkenone-based SST reconstruction represents a local, upwelling-related signal rather than a global insolation related one. Extraordinary increases in diatom and opal concentrations during Heinrich event 1 (H1, 15.5-18.0 ka B.P.) and the Younger Dryas (YD, 13.5-11.5 ka B.P.) are attributed to enhanced upwelling of silica-rich waters and an enlarged upwelling filament, due to more intense NE trade winds. The synchronous increase of CaCO3 and K intensity and the decreased opal and diatoms values mark the occurrence of the Bølling/Allerød (BA, 13.5-15.5 ka B.P.) due to weakened eolian input and more humid conditions on land. Although the high export of diatoms is inextricably linked to upwelling intensity off Mauritania, variability in the nutrient content of the thermocline also plays a decisive role.

Dust fluxes and iron fertilization in Holocene and Last Glacial Maximum climates

Mineral dust aerosols play a major role in present and past climates. To date, we rely on climate models for estimates of dust fluxes to calculate the impact of airborne micronutrients on biogeochemical cycles. Here we provide a new global dust flux data set for Holocene and Last Glacial Maximum (LGM) conditions based on observational data. A comparison with dust flux simulations highlights regional differences between observations and models. By forcing a biogeochemical model with our new data set and using this model's results to guide a millennial-scale Earth System Model simulation, we calculate the impact of enhanced glacial oceanic iron deposition on the LGM-Holocene carbon cycle. On centennial timescales, the higher LGM dust deposition results in a weak reduction of <10?ppm in atmospheric CO2 due to enhanced efficiency of the biological pump. This is followed by a further ~10?ppm reduction over millennial timescales due to greater carbon burial and carbonate compensation.

Hydroclimate records and radiocarbon dating of sediment cores SO189/2 _39KL, SO189/2 _119KL, and SO189/2 _144KL from the eastern tropical Indian Ocean

The response of the tropical climate in the Indian Ocean realm to abrupt climate change events in the North Atlantic Ocean is contentious. Repositioning of the intertropical convergence zone is thought to have been responsible for changes in tropical hydroclimate during North Atlantic cold spells1, 2, 3, 4, 5, but the dearth of high-resolution records outside the monsoon realm in the Indian Ocean precludes a full understanding of this remote relationship and its underlying mechanisms. Here we show that slowdowns of the Atlantic meridional overturning circulation during Heinrich stadials and the Younger Dryas stadial affected the tropical Indian Ocean hydroclimate through changes to the Hadley circulation including a southward shift in the rising branch (the intertropical convergence zone) and an overall weakening over the southern Indian Ocean. Our results are based on new, high-resolution sea surface temperature and seawater oxygen isotope records of well-dated sedimentary archives from the tropical eastern Indian Ocean for the past 45,000 years, combined with climate model simulations of Atlantic circulation slowdown under Marine Isotope Stages 2 and 3 boundary conditions. Similar conditions in the east and west of the basin rule out a zonal dipole structure as the dominant forcing of the tropical Indian Ocean hydroclimate of millennial-scale events. Results from our simulations and proxy data suggest dry conditions in the northern Indian Ocean realm and wet and warm conditions in the southern realm during North Atlantic cold spells.

Nd isotopes of sediments, foraminifera and fish debris of the western North Atlantic

Understanding changes in ocean circulation during the last deglaciation is crucial to unraveling the dynamics of glacial-interglacial and millennial climate shifts. We used neodymium isotope measurements on postdepositional iron-manganese oxide coatings precipitated on planktonic foraminifera to reconstruct changes in the bottom water source of the deep western North Atlantic at the Bermuda Rise. Comparison of our deep water source record with overturning strength proxies shows that both the deep water mass source and the overturning rate shifted rapidly and synchronously during the last deglacial transition. In contrast, any freshwater perturbation caused by Heinrich event 1 could have only affected shallow overturning. These findings show how changes in upper-ocean overturning associated with millennial-scale events differ from those associated with whole-ocean deglacial climate events.

Stable oxygen isotope data of Globorotalia inflata and radiocarbon dates for sediment cores GeoB6211-2, GeoB13862-1, and GeoB6308-3

The Southern Westerly Winds (SWW) exert a crucial influence over the world ocean and climate. Nevertheless, a comprehensive understanding of the Holocene temporal and spatial evolution of the SWW remains a significant challenge due to the sparsity of high-resolution marine archives and appropriate SWW proxies. Here, we present a north-south transect of high-resolution planktonic foraminiferal oxygen isotope records from the western South Atlantic. Our proxy records reveal Holocene migrations of the Brazil- Malvinas Confluence (BMC), a highly sensitive feature for changes in the position and strength of the northern portion of the SWW. Through the tight coupling of the BMC position to the large-scale wind field, the records allow a quantitative reconstruction of Holocene latitudinal displacements of the SWW across the South Atlantic. Our data reveal a gradual poleward movement of the SWW by about 1-1.5° from the early to the mid-Holocene. Afterwards variability in the SWW is dominated by millennial-scale displacements in the order of 1° in latitude with no recognizable longer-term trend. These findings are confronted with results from a state-of-the-art transient Holocene climate simulation using a comprehensive coupled atmosphere-ocean general circulation model. Proxy-inferred and modeled SWW shifts compare qualitatively, but the model underestimates both orbitally forced multi-millennial and internal millennial SWW variability by almost an order of magnitude. The underestimated natural variability implies a substantial uncertainty in model projections of future SWW shifts.

Nitrogen isotope variations along the northeast Pacific margin

Glacial-interglacial changes in sedimentary d15N over the last 120 kyr display a remarkably similar pattern in timing and amplitude in core records extending from the denitrification zone in the eastern tropical North Pacific (ETNP), where subsurface denitrification is active, to the Oregon margin, where no denitrification occurs today. Low d15N values (4-6 per mil) generally characterize glacial stages 2 and 4, and higher d15N values (7-10 per mil) are representative of the Holocene, millennial-scale periods within stage 3, and stage 5. The inferred synchroneity of d15N variations along the entire margin implies that the nitrate isotopic signal produced in the oxygen-poor subsurface waters in the ETNP is rapidly advected northward and recorded at sites far beyond the boundaries of the modern denitrification zone. Similar to d15N, primary production indicators (percent Corg, Ba/Al, and percent opal) show glacial-interglacial as well as millennial-scale variations along the NE Pacific margin, with higher primary production during warm periods. However, the relative phasing between d15N and paleoproduction tracers within individual records changes latitudinally. Whereas d15N and primary production vary approximately synchronously in the midlatitudes, production lags d15N in the ETNP by several kiloyears. This lag calls for a new understanding of the processes driving denitrification in the ETNP. We suggest that oxygen input by the Equatorial Undercurrent as well as local organic matter flux controls denitrification rates in the ETNP. Moreover, the differences in relative timing point to a time-transgressive development of upwelling-favorable winds along the NE Pacific margin after the last glaciation, with those in the north developing several kiloyears earlier.

Multi-proxy analyses of sediment core GeoB13801-2

Over the Uruguayan shelf and uppermost slope the coalescence of northward flowing Subantarctic Shelf Water and southward flowing Subtropical Shelf Water forms a distinct thermohaline front termed the Subtropical Shelf Front (STSF). Running in a SW direction diagonally across the shelf from the coastal waters at 32°S towards the shelf break at ca. 36°S, the STSF represents the shelf-ward extension of the Brazil-Malvinas Confluence zone. This study reconstructs latitudinal STSF shifts during the Holocene based on benthic foraminifera d18O and d13C, total organic carbon, carbonate contents, Ti/Ca, and grain-size distribution from a high-accumulation sedimentary record located at an uppermost continental-slope terrace. Our data provide direct evidence for: (1) a southern STSF position (to the South of the core site) at the beginning of the early Holocene (>9.4 cal ka BP) linked to a more southerly position of the Southern Westerly Winds in combination with restricted shelf circulation intensity due to lower sea level; (2) a gradual STSF northward migration (bypassing the core site towards the North) primarily forced by the northward migration of the Southern Westerly Winds from 9.4 cal ka BP onwards; (3) a relatively stable position of the front in the interval between 7.2 and 4.0 cal ka BP; (4) millennial-scale latitudinal oscillations close to 36°S of the STSF after 4.0 cal ka BP probably linked to the intensification in El Niño Southern Oscillation; and (5) a southward migration of the STSF during the last 200 years possibly linked to anthropogenic influences on the atmosphere.

Stable oxygen isotope, UK'37, productivity, CaCO3 and XRF records from sediment core MD06-3075

Proxy records of hydrologic variability in the West Pacific Warm Pool (WPWP) have revealed wide-scale changes in past convective activity in response to orbital and sub-orbital climate forcings. However, attributing proxy responses to regional changes in WPWP hydrology versus local variations in precipitation requires independent records linking the terrestrial and marine realms. We present high-resolution stable isotope, UK'37 sea-surface temperature, X-ray fluorescence (XRF) core scanning and coccolithophore-derived paleoproductivity records covering the past 120 ka from International Marine Global Change (IMAGES) Program Core MD06-3075 (6°29' N, 125°50' E, water depth 1878 m), situated in the Davao Gulf on the southern side of Mindanao. XRF-derived log(Fe/Ca) records provide a robust proxy for runoff-driven sedimentary discharge from Mindanao, whilst past changes in local productivity are associated with variable freshwater runoff and stratification of the surface layer. Significant precessional-scale variability in sedimentary discharge occurred during Marine Isotope Stage (MIS) 5, with peaks in discharge contemporaneous with Northern Hemisphere summer insolation minima. We attribute these changes to the latitudinal migration of the Intertropical Convergence Zone (ITCZ) over the WPWP together with variability in the strength of the Walker circulation acting on precessional timescales. Between 60 and 15 ka sedimentary discharge at Mindanao was muted, displaying little orbital- or millennial-scale variability, likely in response to weakened precessional insolation forcing and lower sea level driving increased subsidence of air masses over the exposed Sunda Shelf. These results highlight the high degree of local variability in the precipitation response to past climate changes in the WPWP.

Biomarker and x-ray diffraction data from Site U1313 (MIS 16 - 9)

A reconstruction of Milankovitch to millennial-scale variability of sea-surface temperature (SST) and sea-surface productivity in the Pleistocene mid-latitude North Atlantic Ocean (MIS 16-9) and its relationship to ice sheet instability was carried out on sediments from IODP Site U1313. This reconstruction is based on alkenone and n-alkane concentrations, Uk37' index, total organic carbon (TOC) and carbonate contents, X-Ray diffraction (XRD) data, magnetic susceptibility, and accumulation rates. Increased input of ice-rafted debris (IRD) occurred during MIS 16, 12, and 10, characterized by high concentrations of dolomite, quartz, and feldspars and elevated accumulation rates of terrigenous matter. Minimum input values of terrigenous matter, on the other hand, were determined for MIS 13 and 11. Peak values of dolomite, coinciding with quartz, plagioclase, and kalifeldspar peaks and maxima in long-chain n-alkanes indicative for land plants, are interpreted as Heinrich-like Events related to sudden instability of the Laurentide Ice Sheet during early and late (deglacial) phases of the glacials. The coincidence of increased TOC values with elevated absolute concentrations of alkenones suggest increased glacial productivity, probably due to a more southern position of the Polar Front. Alkenone-based SST reached absolute maxima of about 19°C during MIS 11.3 and absolute minima of <10°C during MIS 12 and 10. Within MIS 11, prominent cooling events (MIS 11.22 and 11.24) occurred. The absolute SST minima recorded directly before and after the glacial maxima MIS 10.2 and 12.2, are related to Heinrich-like Event meltwater pulses, as supported by the coincidence of SST minima and maxima in C37:4 alkenones and dolomite. These sudden meltwater pulses - especially during Terminations IV and V - probably caused a collapse of phytoplankton productivity as indicated by the distinct drop in alkenone concentrations. Ice-sheet disintegration and subsequent surges and outbursts of icebergs and meltwater discharge may have been triggered by increased insolation in the Northern High Latitudes.

Stable carbon and oxygen isotope record od IODP Hole 306-U1313

The Mid-Pleistocene transition (MPT) was the time when quasi-periodic (? 100 kyr), high-amplitude glacial variability developed in the absence of any significant change in the character of orbital forcing, leading to the establishment of the characteristic pattern of late Pleistocene climate variability. It has long been known that the interval around 900 ka stands out as a critical point of the MPT, when major glaciations started occurring most notably in the northern hemisphere. Here we examine the record of climatic conditions during this significant interval, using high-resolution stable isotope records from benthic and planktonic foraminifera from a sediment core in the North Atlantic (Integrated Ocean Drilling Program Expedition 306, Site U1313). We have considered the time interval from late in Marine Isotope Stage (MIS) 23 to MIS 20 (910 to 790 ka). Our data indicate that interglacial MIS 21 was a climatically unstable period and was broken into four interstadial periods, which have been identified and correlated across the North Atlantic region. These extra peaks tend to contradict previous studies that interpreted the MIS 21 variability as consisting essentially of a linear response to cyclical changes in orbital parameters. Cooling events in the surface record during MIS 21 were associated with low benthic carbon isotope excursions, suggesting a coupling between surface temperature changes and the strength of the Atlantic meridional overturning circulation. Time series analysis performed on the whole interval indicates that benthic and planktonic oxygen isotopes have significant concentrations of spectral power centered on periods of 10.7 kyr and 6 kyr, which is in agreement with the second and forth harmonic of precession. The excellent correspondence between the foraminifera d18O records and insolation variations at the Equator in March and September suggests that a mechanism related to low-latitude precession variations, advected to the high latitudes by tropical convective processes, might have generated such a response. This scenario accounts for the presence of oscillations at frequencies equal to precession harmonics at Site U1313, as well as the occurrence of higher amplitude oscillations between the MIS22/21 transition and most of MIS 21, times of enhanced insolation variability.

Stable carbon and oxygen isotope record of IODP Site 306-U1313, MIS 23-19

Since the seminal work by Hays et al. (1976), a plethora of studies has demonstrated a correlation between orbital variations and climatic change. However, information on how changes in orbital boundary conditions affected the frequency and amplitude of millennial-scale climate variability is still fragmentary. The Marine Isotope Stage (MIS) 19, an interglacial centred at around 785 ka, provides an opportunity to pursue this question and test the hypothesis that the long-term processes set up the boundary conditions within which the short-term processes operate. Similarly to the current interglacial, MIS 19 is characterised by a minimum of the 400-kyr eccentricity cycle, subdued amplitude of precessional changes, and small amplitude variations in insolation. Here we examine the record of climatic conditions during MIS 19 using high-resolution stable isotope records from benthic and planktonic foraminifera from a sedimentary sequence in the North Atlantic (Integrated Ocean Drilling Program Expedition 306, Site U1313) in order to assess the stability and duration of this interglacial, and evaluate the climate system's response in the millennial band to known orbitally induced insolation changes. Benthic and planktonic foraminiferal d18O values indicate relatively stable conditions during the peak warmth of MIS 19, but sea-surface and deep-water reconstructions start diverging during the transition towards the glacial MIS 18, when large, cold excursions disrupt the surface waters whereas low amplitude millennial scale fluctuations persist in the deep waters as recorded by the oxygen isotope signal. The glacial inception occurred at ~779 ka, in agreement with an increased abundance of tetra-unsaturated alkenones, reflecting the influence of icebergs and associated meltwater pulses and high-latitude waters at the study site. After having combined the new results with previous data from the same site, and using a variety of time series analysis techniques, we evaluate the evolution of millennial climate variability in response to changing orbital boundary conditions during the Early-Middle Pleistocene. Suborbital variability in both surface- and deep-water records is mainly concentrated at a period of ~11 kyr and, additionally, at ~5.8 and ~3.9 kyr in the deep ocean; these periods are equal to harmonics of precession band oscillations. The fact that the response at the 11 kyr period increased over the same interval during which the amplitude of the response to the precessional cycle increased supports the notion that most of the variance in the 11 kyr band in the sedimentary record is nonlinearly transferred from precession band oscillations. Considering that these periodicities are important features in the equatorial and intertropical insolation, these observations are in line with the view that the low-latitude regions play an important role in the response of the climate system to the astronomical forcing. We conclude that the effect of the orbitally induced insolation is of fundamental importance in regulating the timing and amplitude of millennial scale climate variability.