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.

Spectral reflectance library of selected biotic and abiotic coral reef features in Glovers Reef, Belize

Underwater spectral reflectance was measured for selected biotic and abiotic coral reef features of Glovers Reef, Belize from March 6 - 10, 2005. Spectral reflectance's of 63 different benthic types were obtained in-situ. An Ocean Optics USB2000 spectrometer was deployed in an custom made underwater housing with a 0.5 m fiber-optic probe mounted next to an artificial light source. Spectral readings were collected with the probe (bear fibre) about 5 cm from the target to ensure that the target would fill the field of view of the fiber optic (FOV diameter ~4.4 cm), as well as to reduce the attenuating effect of the intermediate water (Roelfsema et al., 2006). Spectral readings included for one target included: 1 reading of the covered spectral fibre to correct for instrument noise, 1 reading of spectralon panel mounted on divers wrist to measure incident ambient light, and 8 readings of the target. Spectral reflectance was calculated for each target by first subtracting the instrument noise reading from each other reading. The corrected target readings were then divided by the corrected spectralon reading resulting in spectral reflectance of each target reading. An average target spectral reflectance was calculated by averaging the eight individual spectral reflectance's of the target. If an individual target spectral reflectance was visual considered an outlier, it was not included in the average spectral reflectance calculation. See Roelfsema at al. (2006) for additional info on the methodology of underwater spectra collection.

Spectral reflectance library of selected biotic and abiotic coral reef features in Heron Reef

Underwater spectral reflectance was measured for selected biotic and abiotic coral reef features of Heron Reef from June 25-30, 2006. Spectral reflectance's of 105 different benthic types were obtained in-situ. An Ocean Optics USB2000 spectrometer was deployed in an custom made underwater housing with a 0.5 m fiber-optic probe mounted next to an artificial light source. Spectral readings were collected with the probe(bear fibre) about 5 cm from the target to ensure that the target would fill the field of view of the fiber optic (FOV diameter ~4.4 cm), as well as to reduce the attenuating effect of the intermediate water (Roelfsema et al., 2006). Spectral readings included for one target included: 1 reading of the covered spectral fibre to correct for instrument noise, 1 reading of spectralon panel mounted on divers wrist to measure incident ambient light, and 8 readings of the target. Spectral reflectance was calculated for each target by first subtracting the instrument noise reading from each other reading. The corrected target readings were then divided by the corrected spectralon reading resulting in spectral reflectance of each target reading. An average target spectral reflectance was calculated by averaging the eight individual spectral reflectance's of the target. If an individual target spectral reflectance was visual considered an outlier, it was not included in the average spectral reflectance calculation. See Roelfsema at al. (2006) for additional info on the methodology of underwater spectra collection.

Sea surface temperature and salinity reconstruction of sediment core GeoB3129-1

A sediment core from the western tropical Atlantic covering the last 21,000 yr has been analysed for centennial scale reconstruction of sea surface temperature (SST) and ice volume-corrected oxygen isotopic composition of sea water (delta18O(ivc-sw)) using Mg / Ca and delta18O of the shallow dwelling planktonic foraminifer Globigerinoides ruber (white). At a period between 15.5 and 17.5 kyr BP, the Mg / Ca SST and delta18O(ivc-sw), a proxy for sea surface salinity (SSS), reveals a warming of around 2.5 °C along with an increase in salinity. A second period of pronounced warming and SSS increase occurred between 11.6 and 13.5 kyr BP. Within age model uncertainties, both warming intervals were synchronous with air temperature increase over Antarctica and ice retreat in the southern South Atlantic and terminated with abrupt centennial scale SSS decrease and slight SST cooling in conjunction with interglacial reactivation of the meridional overturning circulation (MOC). We suggest that during these warm intervals, production of saline and warm water of the North Brazil Current resulted in pronounced heat and salt accumulation, and was associated with warming in the southern Atlantic, southward displacement of the intertropical convergence zone and weakened MOC. At the termination of the Younger Dryas and Heinrich event 1, intensification of cross-equatorial heat and salt transport caused centennial scale cooling and freshening of the western tropical Atlantic surface water. This study shows that the western tropical Atlantic served as a heat and salt reservoir during deglaciation. The sudden release of accumulated heat and salt at the end of Younger Drays and Heinrich event 1 may have contributed to the rapid reinvigoration of the Atlantic MOC.

(Table 1) Age determination of sediment core SO147_106KL

Here we present a high-resolution marine sediment record from the El Niño region off the coast of Peru spanning the last 20,000 years. Sea surface temperature, photosynthetic pigments, and a lithic proxy for El Niño flood events on the continent are used as paleo-El Niño-Southern Oscillation proxy data. The onset of stronger El Niño activity in Peru started around 17,000 calibrated years before the present, which is later than modeling experiments show but contemporaneous with the Heinrich event 1. Maximum El Niño activity occurred during the early and late Holocene, especially during the second and third millennium B.P. The recurrence period of very strong El Niño events is 60-80 years. El Niño events were weak before and during the beginning of the Younger Dryas, during the middle of the Holocene, and during medieval times. The strength of El Niño flood events during the last millennium has positive and negative relationships to global and Northern Hemisphere temperature reconstructions.

Event layer chronology from Lake Suigetsu (Japan, SG06 Project) for the last 40 kyr

Event layers in lake sediments are indicators of past extreme events, mostly the results of floods or earthquakes. Detailed characterisation of the layers allows the discrimination of the sedimentation processes involved, such as surface runoff, landslides or subaqueous slope failures. These processes can then be interpreted in terms of their triggering mechanisms. Here we present a 40 kyr event layer chronology from Lake Suigetsu, Japan. The event layers were characterised using a multi-proxy approach, employing light microscopy and µXRF for microfacies analysis. The vast majority of event layers in Lake Suigetsu was produced by flood events (362 out of 369), allowing the construction of the first long-term, quantitative (with respect to recurrence) and well dated flood chronology from the region. The flood layer frequency shows a high variability over the last 40 kyr, and it appears that extreme precipitation events were decoupled from the average long-term precipitation. For instance, the flood layer frequency is highest in the Glacial at around 25 kyr BP, at which time Japan was experiencing a generally cold and dry climate. Other cold episodes, such as Heinrich Event 1 or the Late Glacial stadial, show a low flood layer frequency. Both observations together exclude a simple, straightforward relationship with average precipitation and temperature. We argue that, especially during Glacial times, changes in typhoon genesis/typhoon tracks are the most likely control on the flood layer frequency, rather than changes in the monsoon front or snow melts. Spectral analysis of the flood chronology revealed periodic variations on centennial and millennial time scales, with 220 yr, 450 yr and a 2000 yr cyclicity most pronounced. However, the flood layer frequency appears to have not only been influenced by climate changes, but also by changes in erosion rates due to, for instance, earthquakes.

(Table 1) Age determination on monospecies planktonic foraminifera of sediment core AN25-934

Reconstruction of regional climate and the Okhotsk Sea (OS) environment for the Last Glacial Maximum (LGM), deglaciation and Holocene were performed on the basis of high-resolution records of ice rafted debris (IRD), CaCO3, opal, total organic carbon (TOC), biogenic Ba (Ba_bio) and redox sensitive element (Mn, Mo) content, and diatom and pollen results of four cores that form a north-southern transect. Age models of the studied cores were earlier established by AMS 14C data, oxygen - isotope chronostratigraphy and tephrochronology. According to received results, since 25 ka the regional climate and OS environmental conditions have changed synchronously with LGM condition, cold Heinrich event 1, Bølling -Allerød (BA) warming, Younger Dryas (YD) cooling and Pre-Boreal (PB) warming recorded in the Greenland ice core, North Atlantic sediment, and China cave stalagmites. Calculation of IRD MAR in sediment of north-south transect cores indicate an increase of sea ice formation several times in the glacial OS as compared to the Late Holocene. Accompanying ice formation, increased brine rejection and the larger potential density of surface water at the north shelf due to a drop of glacial East Asia summer monsoon precipitation and Amur River run off, led to strong enhancement of the role of the OS in glacial North Pacific Intermediate Water (NPIW) formation. The remarkable increase in OS productivity during BA and PB warming was probably related with significant reorganisation of the North Pacific deep water ventilation and nutrient input into the NPIW and OS Intermediate Water (OSIW). Seven Holocene OS millennial cold events based on the elevated values of the detrended IRD stack record over the IRD broad trend in the sediments of the studied cores have occurred synchronously with cold events recorded in the North Atlantic, Greenland ice cores and China cave stalagmites after 9 ka. Diatom production in the OS were mostly controlled by sea ice cover changes and surface water stratification induced by sea-ice melting; therefore significant opal accumulation in sediments of this basin begin from 4-6 ka ago simultaneously with a remarkable decrease of sea ice cover.

(Table 1) Age determination of sediment cores HU89038_8 and OCE326-GGC5

Models indicate that a complete shutdown of deep and intermediate water production is a possible consequence of extreme climate conditions in the northern North Atlantic, and the high ratio of 231Pa to 230Th on Bermuda Rise is evidence that this might have happened ?17 ka during Heinrich event 1 (H1). However, new radiocarbon data from bivalves that lived at ?4.6 km on the Bermuda Rise during H1 lead to a different conclusion. The bivalve data do indeed indicate ventilation of the deep western North Atlantic was suppressed during H1 but not as much as it was during the last glacial maximum. We propose that high diatom flux to the Bermuda Rise during H1 is at least in part responsible for increased 231Pa/230Th at that time. Although we cannot say for sure why opal production was so high in a gyre center location at that time, increased leakage of silica rich waters from the Southern Ocean to the North Atlantic is one possibility.

(Table 1) Oxygen isotopic composition of interstitial waters from ODP Hole 161-976

Interstitial waters recovered from Ocean Drilling Program, Leg 161, site 976 in the western Mediterranean Sea are used in conjunction with a numerical model to constrain the delta18O of seawater in the basin since the Last Glacial Maximum, including Sapropel Event 1. To resolve the oxygen isotopic composition of the deep Mediterranean, we use a model that couples fluid diffusion with advective transport, thus producing a profile of seawater delta18O variability that is unaffected by glacial-interglacial variations in marine temperature. Comparing our reconstructed seawater delta18O to recent determinations of 1.0 per mil for the mean ocean change in glacial-interglacial delta18O due to the expansion of global ice volume, we calculate an additional 0.2 per mil increase in Mediterranean delta18O caused by local evaporative enrichment. This estimate of delta18O change, due to salinity variability, is smaller than previous studies have proposed and demonstrates that Mediterranean records of foraminiferal calcite delta18O from the last glacial period include a strong temperature component. Paleotemperatures determined in combination with a stacked record of foraminiferal calcite depict almost 9°C of regional cooling for the Last Glacial Maximum. Model results suggest a decrease of ~1.1 per mil in seawater delta18O relative to the modern value caused by increased freshwater input and reduced salinity accompanying the formation of the most recent sapropel. The results additionally indicate the existence of isotopically light water circulating down to bottom water depths, at least in the western Mediterranean, supporting the existence of an 'anti-estuarine' thermohaline circulation pattern during Sapropel Event 1.

Radiocarbon, isotope, biogenic, and 230Th measurements for sediment core RC10-196

The subarctic North Pacific Ocean holds a large CO2 reservoir that is currently isolated from the atmosphere by a low-salinity layer. It has recently been hypothesized that the reorganization of these high-CO2 waters may have played a crucial role in the degassing of carbon dioxide to the atmosphere during the last deglaciation. This reorganization would leave some imprint on paleo-productivity records. Here we present 230Th-normalized biogenic fluxes from an intermediate depth sediment core in the Northwest Pacific (RC10-196, 54.7°N, 177.1°E, 1007 m) and place them within the context of a synthesis of previously-published biogenic flux data from 49 deep-sea cores north of 20°N, ranging from 420 to 3968 m water depth. The 230Th-normalized opal, carbonate, and organic carbon fluxes from RC10-196 peak approximately 13,000 calendar years BP during the Bølling/Allerød (B/A) period. Our data synthesis suggests that biogenic fluxes were in general lowest during the last glacial period, increased somewhat in the Northwest Pacific during Heinrich Event 1, and reached a maximum across the entire North Pacific during the B/A period. We evaluate several mechanisms as possible drivers of deglacial change in biogenic fluxes in the North Pacific, including changes in preservation, sediment focusing, sea ice extent, iron inputs, stratification, and circulation shifts initiated in the North Atlantic and North Pacific. Our analysis suggests that while micronutrient sources likely contributed to some of the observed changes, the heterogeneity in timing of glaciogenic retreat and sea level make these mechanisms unlikely causes of region-wide contemporaneous peaks in export production. We argue that paleo-observations are most consistent with ventilation increases in both the North Pacific (during H1) and North Atlantic (during B/A) being the primary drivers of increases in biogenic flux during the deglaciation, as respectively they were likely to bring nutrients to the surface via increased vertical mixing and shoaling of the global thermocline.

Age determination, foraminiferal abundance and stable isotope ratios of sediment core CH82-20

A core from the Mid-Atlantic Ridge at 43.5°N and ~3 km water depth shows distinct evidence of the deglacial events known as Heinrich event 1 (probably the marine equivalent of Oldest Dryas cooling in Europe) and the Younger Dryas. The Heinrich event, dated at three levels to between 14.3 and 15.0 ka, is marked by a minimum in foraminifera per gram, by maxima in rates of sedimentation, ice rafted debris per gram, and relative abundance of N. pachyderma (s.), and by a delta18O minimum in planktonic foraminifera. The Younger Dryas event is marked by peak abundance of N. pachyderma (s.) and a planktonic delta18O maximum. Benthic foraminiferal delta13C reaches minimum values during both the Heinrich event and the Younger Dryas. Our data indicate pronounced changes in surface water properties were coupled with reduced production of North Atlantic Deep Water at each of these times.

Multi-proxy analyses of sediment cores from the Pakistan continental margin

Late Holocene laminated sediments from a core transect centred in the oxygen minimum zone (OMZ) impinging at the continental slope off Pakistan indicate stable oxygen minimum conditions for the past 7000 calendar years. High SW-monsoon-controlled biological productivity and enhanced organic matter preservation during this period is reflected in high contents of total organic carbon (TOC) and redox-sensitive elements (Ni, V), as well as by a low-diversity, high-abundance benthic foraminiferal Buliminacea association and high abundance of the planktonic species Globigerina bulloides indicative of upwelling conditions. Surface-water productivity was strongest during SW monsoon maxima. Stable OMZ conditions (reflected by laminated sediments) were found also during warm interstadial events (Preboreal, Bølling-Allerød, and Dansgaard-Oeschger events), as well as during peak glacial times (17-22.5 ka, all ages in calendar years). Sediment mass accumulation rates were at a maximum during the Preboreal and Younger Dryas periods due to strong riverine input and mobilisation of fine-grained sediment coinciding with rapid deglacial sea-level rise, whereas eolian input generally decreased from glacial to interglacial times. In contrast, the occurrence of bioturbated intervals from 7 to 10.5 ka (early Holocene), in the Younger Dryas (11.7-13 ka), from 15 to 17 ka (Heinrich event 1) and from 22.5 to 25 ka (Heinrich event 2) suggests completely different conditions of oxygen-rich bottom waters, extremely low mass and organic carbon accumulation rates, a high-diversity benthic fauna, all indicating lowered surface-water productivity. During these intervals the OMZ was very poorly developed or absent and a sharp fall of the aragonite compensation depth favoured the preservation of pteropods. The abundance of lithogenic proxies suggests aridity and wind transport by northwesterly or northeasterly winds during these periods coinciding with the North Atlantic Heinrich events and dust peaks in the Tibetan Loess records. The correlation of the monsoon-driven OMZ variability in the Arabian Sea with the rapid climatic fluctuations in the high northern latitudes suggests a close coupling between the climates of the high and low latitudes at a global scale.