North Atlantic IODP Site 303-U1308 Middle Pleistocene Heinrich Events record

Heinrich events are well documented for the last glaciation, but little is known about their occurrence in older glacial periods of the Pleistocene. Here we report scanning XRF and bulk carbonate d18O results from Integrated Ocean Drilling Program Site U1308 (reoccupation of Deep Sea Drilling Project Site 609) that are used to develop proxy records of ice-rafted detritus (IRD) for the last ~1.4 Ma. Ca/Sr is used as an indicator of IRD layers that are rich in detrital carbonate (i.e., Heinrich layers), whereas Si/Sr reflects layers that are poor in biogenic carbonate and relatively rich in detrital silicate minerals. A pronounced change occurred in the composition and frequency of IRD at ~640 ka during marine isotope stage (MIS) 16, coinciding with the end of the middle Pleistocene transition. At this time, "Hudson Strait" Heinrich layers suddenly appeared in the sedimentary record of Site U1308, and the dominant period of the Si/Sr proxy shifted from 41 ka prior to 640 ka to 100 ka afterward. The onset of Heinrich layers during MIS 16 represents either the initiation of surging of the Laurentide Ice Sheet (LIS) off Hudson Strait or the first time icebergs produced by this process survived the transport to Site U1308. We speculate that ice volume (i.e., thickness) and duration surpassed a critical threshold during MIS 16 and activated the dynamical processes responsible for LIS instability in the region of Hudson Strait. We also observe a strong coupling between IRD proxies and benthic d13C variation at Site U1308 throughout the Pleistocene, supporting a link between iceberg discharge and weakening of thermohaline circulation in the North Atlantic.

Investigation of the seasonality of N. pachyderma in the North Atlantic Ocean during Heinrich Stadial 1 using an ecosystem modeling approach

Fossil shells of planktonic foraminifera serve as the prime source of information on past changes in surface ocean conditions. Because the population size of planktonic foraminifera species changes throughout the year, the signal preserved in fossil shells is biased towards the conditions when species production was at its maximum. The amplitude of the potential seasonal bias is a function of the magnitude of the seasonal cycle in production. Here we use a planktonic foraminifera model coupled to an ecosystem model to investigate to what degree seasonal variations in production of the species Neogloboquadrina pachyderma may affect paleoceanographic reconstructions during Heinrich Stadial 1 (~18-15 cal. ka B.P.) in the North Atlantic Ocean. The model implies that during Heinrich Stadial 1 the maximum seasonal production occurred later in the year compared to the Last Glacial Maximum (~21-19 cal. ka B.P.) and the pre-industrial era north of 30 ºN. A diagnosis of the model output indicates that this change reflects the sensitivity of the species to the seasonal cycle of sea-ice cover and food supply, which collectively lead to shifts in the modeled maximum production from the Last Glacial Maximum to Heinrich Stadial 1 by up to six months. Assuming equilibrium oxygen isotopic incorporation in the shells of N. pachyderma, the modeled changes in seasonality would result in an underestimation of the actual magnitude of the meltwater isotopic signal recorded by fossil assemblages of N. pachyderma wherever calcification is likely to take place.

Tropical climate and vegetation simulations during the Heinrich event 1 using an Earth System Model of Intermediate Complexity (EMIC) - the University of Victoria Earth System-Climate Model (UVic ESCM)

Abrupt climate changes from 18 to 15 thousand years before present (kyr BP) associated with Heinrich Event 1 (HE1) had a strong impact on vegetation patterns not only at high latitudes of the Northern Hemisphere, but also in the tropical regions around the Atlantic Ocean. To gain a better understanding of the linkage between high and low latitudes, we used the University of Victoria (UVic) Earth System-Climate Model (ESCM) with dynamical vegetation and land surface components to simulate four scenarios of climate-vegetation interaction: the pre-industrial era, the Last Glacial Maximum (LGM), and a Heinrich-like event with two different climate backgrounds (interglacial and glacial). We calculated mega-biomes from the plant-functional types (PFTs) generated by the model to allow for a direct comparison between model results and palynological vegetation reconstructions. Our calculated mega-biomes for the pre-industrial period and the LGM corresponded well with biome reconstructions of the modern and LGM time slices, respectively, except that our pre-industrial simulation predicted the dominance of grassland in southern Europe and our LGM simulation resulted in more forest cover in tropical and sub-tropical South America. The HE1-like simulation with a glacial climate background produced sea-surface temperature patterns and enhanced inter-hemispheric thermal gradients in accordance with the "bipolar seesaw" hypothesis. We found that the cooling of the Northern Hemisphere caused a southward shift of those PFTs that are indicative of an increased desertification and a retreat of broadleaf forests in West Africa and northern South America. The mega-biomes from our HE1 simulation agreed well with paleovegetation data from tropical Africa and northern South America. Thus, according to our model-data comparison, the reconstructed vegetation changes for the tropical regions around the Atlantic Ocean were physically consistent with the remote effects of a Heinrich event under a glacial climate background.