Age determination and Heinrich Events in sediment cores SU90-08 and KET80-03

The Mediterranean Sea hydrology at the time of the Heinrich formation in the North Atlantic Ocean was analyzed by comparing sea surface temperatures (SSTs) and oxygen isotope composition of seawater (dw) changes during the past 75 kyr in two marine cores. These were compared to the palynological variations derived in the Mediterranean Sea core. During the last glacial the two oceanic SST records show similar and synchronous patterns, with several long-term cooling periods, ending by abrupt SST increases. At the time of the Heinrich events, cold SSTs and low salinity prevailed in the Mediterranean Sea. The freshwater budget was similar to the modern one, permitting the presence of a mixed forest on the Mediterranean borderlands. The post-Heinrich periods are marked by a freshwater budget decrease, limiting oak and fir tree growth in the Mediterranean region. Increase of precipitation or reduction of evaporation is observed before the Heinrich episode, and is associated with a well-developed mixed Mediterranean forest.

A time-slice study in Heinrich Stadial 1 and Mid Holocene sediment cores off northern Senegal

Sediment dynamics in limnic, fluvial and marine environments can be assessed by granulometric and rock-magnetic methodologies. While classical grain-size analysis by sieving or settling mainly bears information on composition and transport, the magnetic mineral assemblages reflect to a larger extent the petrology and weathering conditions in the sediment source areas. Here, we combine both methods to investigate Late Quaternary marine sediments from five cores along a transect across the continental slope off Senegal. This region near the modern summer Intertropical Convergence Zone is particularly sensitive to climate change and receives sediments from several aeolian, fluvial and marine sources. From each of the investigated five GeoB sediment cores (494-2956 m water depth) two time slices were processed which represent contrasting climatic conditions: the arid Heinrich Stadial 1 (~ 15 kyr BP) and the humid Mid Holocene (~ 6 kyr BP). Each sediment sample was split into 16 grain-size fractions ranging from 1.6 to 500 µm. Concentration and grain-size indicative magnetic parameters (susceptibility, SIRM, HIRM, ARM and ARM/IRM) were determined at room temperature for each of these fractions. The joint consideration of whole sediment and magnetic mineral grain-size distributions allows to address several important issues: (i) distinction of two aeolian sediment fractions, one carried by the north-easterly trade winds (40-63 µm) and the other by the overlying easterly Harmattan wind (10-20 µm) as well as a fluvial fraction assigned to the Senegal River (< 10 µm); (ii) identification of three terrigenous sediment source areas: southern Sahara and Sahel dust (low fine-grained magnetite amounts and a comparatively high haematite content), dust from Senegalese coastal dunes (intermediate fine-grained magnetite and haematite contents) and soils from the upper reaches of the Senegal River (high fine-grained magnetite content); (iii) detection of partial diagenetic dissolution of fine magnetite particles as a function of organic input and shore distance; (iv) analysis of magnetic properties of marine carbonates dominating the grain-size fractions 63-500 µm.

Abrupt shifts of the Sahara-Sahel boundary during Heinrich stadials demonstrated on a sediment core transect

Relict dune fields that are found as far south as 14° N in the modern-day African Sahel are testament to equatorward expansions of the Sahara desert during the Late Pleistocene. However, the discontinuous nature of dune records means that abrupt millennial-timescale climate events are not always resolved. High-resolution marine core studies have identified Heinrich stadials as the dustiest periods of the last glacial in West Africa although the spatial evolution of dust export on millennial timescales has so far not been investigated. We use the major-element composition of four high-resolution marine sediment cores to reconstruct the spatial extent of Saharan-dust versus river-sediment input to the continental margin from West Africa over the last 60 ka. This allows us to map the position of the sediment composition corresponding to the Sahara-Sahel boundary. Our records indicate that the Sahara-Sahel boundary reached its most southerly position (13° N) during Heinrich stadials and hence suggest that these were the periods when the sand dunes formed at 14° N on the continent. Heinrich stadials are associated with cold North Atlantic sea surface temperatures which appear to have triggered abrupt increases of aridity and wind strength in the Sahel. Our study illustrates the influence of the Atlantic meridional overturning circulation on the position of the Sahara-Sahel boundary and on global atmospheric dust loading.

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.

Results of Heinrich event 1 simulation from the University of Victoria Earth System-Climate Model (UVic ESCM) and the Community Climate System Model (CCSM3)

We investigated changes in tropical climate and vegetation cover associated with abrupt climate change during Heinrich Event 1 (HE1, ca. 17.5 ka BP) using two different global climate models: the University of Victoria Earth System-Climate Model (UVic ESCM) and the Community Climate System Model version 3 (CCSM3). Tropical South American and African pollen records suggest that the cooling of the North Atlantic Ocean during HE1 influenced the tropics through a southward shift of the rain belt. In this study, we simulated the HE1 by applying a freshwater perturbation to the North Atlantic Ocean. The resulting slowdown of the Atlantic Meridional Overturning Circulation was followed by a temperature seesaw between the Northern and Southern Hemispheres, as well as a southward shift of the tropical rain belt. The shift and the response pattern of the tropical vegetation around the Atlantic Ocean were more pronounced in the CCSM3 than in the UVic ESCM simulation. For tropical South America, opposite changes in tree and grass cover were modeled around 10° S in the CCSM3 but not in the UVic ESCM. In tropical Africa, the grass cover increased and the tree cover decreased around 15° N in the UVic ESCM and around 10° N in the CCSM3. In the CCSM3 model, the tree and grass cover in tropical Southeast Asia responded to the abrupt climate change during the HE1, which could not be found in the UVic ESCM. The biome distributions derived from both models corroborate findings from pollen records in southwestern and equatorial western Africa as well as northeastern Brazil.

Biomarker in Heinrich event layers

There are controversies regarding the origin of Heinrich layer 3 (H3), the massive ice-rafting and meltwater event in the North Atlantic during the last glacial cycle spanning a time window between 29 and 30 kyr B.P. Some argue in favor of a Laurentide Ice Sheet source similar to other Heinrich layers, while a contending view argues for the European ice sheet source. Existing geochemical proxies such as 40Ar/39Ar, 206Pb/204Pb, or epsilon-Nd, etc., could not be used to distinguish among various sources of ice-rafted debris in H3 because of their low abundances, suggesting a background glacial sediment signal. In order to circumvent this problem a biomarker-based approach is used to characterize the provenance of H layers 2, 3, and 4 and other non-Heinrich layers. The presence of hopanes and steranes and their aromatic counterparts in the H layers is incompatible with Recent sediments and is attributed to the transportation of organic matter because of the glacial erosion of source rocks. The most diagnostic and useful signatures of this ancient organic matter in the H layers are the dominance of C34 hopanoids over C33 and the occurrence of isorenieratane along with palaerenieratane. Biomarkers signatures in H layers 2 and 3 of the Labrador Sea suggest no difference in their source. Hydrocarbon distributions suggest that these sediments were derived from the Middle to Late Ordovician and Silurian source rocks of the Hudson Bay of eastern Canada. Biomarker data of the H layer 4 from the northwest Atlantic reveal that the sediments of this layer have a similar source to the H layers in the Labrador Sea.