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.

Treeline dynamics in Siberia under changing climates as inferred from an individual-based model for Larix, supplementary material

Siberian boreal forests are expected to expand northwards in the course of global warming. However, processes of the treeline ecotone transition, as well astiming and related climate feedbacks are still not understood. Here, we present 'Larix Vegetation Simulator' LAVESI, an individual-based spatially-explicit model that can simulate Larix gmelinii (RUPR.) RUPR. stand dynamics in an attempt to improve our understanding about past and future treeline movements under changing climates. The relevant processes (growth, seed production and dispersal, establishment and mortality) are incorporated and adjusted to observation data mainly gained from the literature. Results of a local sensitivity analysis support the robustness of the model's parameterization by giving relatively small sensitivity values. We tested the model by simulating tree stands under modern climate across the whole Taymyr Peninsula, north-central Siberia (c. 64-80° N; 92-119° E). We find tree densities similar to observed forests in the northern to mid-treeline areas, but densities are overestimated in the southern parts of the simulated region. Finally, from a temperature-forcing experiment, we detect that the responses of tree stands lag the hypothetical warming by several decades, until the end of 21st century. With our simulation experiments we demonstrate that the newly-developed model captures the dynamics of the Siberian latitudinal treeline.

Dinoflagellate cyst assemblages, oxygen isotope records of foraminifera, and alkenone and Mg/Ca-based SST estimates for the time interval 3.40-3.18 Ma (Late Pliocene) from five North Atlantic and Caribbean sediment cores

The early Late Pliocene (3.6 to ~3.0 million years ago) is the last extended interval in Earth's history when atmospheric CO2 concentrations were comparable to today's and global climate was warmer. Yet a severe global glaciation during marine isotope stage (MIS) M2 interrupted this phase of global warmth ~3.30 million years ago, and is seen as a premature attempt of the climate system to establish an ice-age world. Our geochemical and palynological records from five marine sediment cores along a Caribbean to eastern North Atlantic transect show that increased Pacific-to-Atlantic flow via the Central American Seaway weakened the North Atlantic Current (NAC) and attendant northward heat transport prior to MIS M2. The consequent cooling of the northern high latitude oceans permitted expansion of the Greenland ice sheet during MIS M2, despite near-modern atmospheric CO2 concentrations. Before and after MIS M2, heat transport via the NAC was crucial in maintaining warm climates comparable to those predicted for the end of this century.

Radiocarbon dates of sediments from ODP Site 161-976

Multiproxy paleoenvironmental records (pollen and planktonic isotope) from Ocean Drilling Program Site 976 (Alboran Sea) document rapid ocean and climate variations during the last glacial that follow the Dansgaard-Oeschger climate oscillations seen in the Greenland ice core records, thus suggesting a close link of the Mediterranean climate swings with North Atlantic climates. Continental conditions rapidly oscillated through cold-arid and warm-wet conditions in the course of stadial-interstadial climate jumps. At the time of Heinrich events, i.e., maximum meltwater flux to the North Atlantic, western Mediterranean marine microflora and microfauna show rapid cooling correlated with increasing continental dryness. Enhanced aridity conceivably points to prolonged wintertime stability of atmospheric high-pressure systems over the southwestern Mediterranean in conjunction with cooling of the North Atlantic.

Particle concentration and size distribution in the NGRIP ice core

A novel laser microparticle detector used in conjunction with continuous sample melting has provided a more than 1500 m long record of particle concentration and size distribution of the NGRIP ice core, covering continuously the period approximately from 9.5-100 kyr before present; measurements were at 1.65 m depth resolution, corresponding to approximately 35-200 yr. Particle concentration increased by a factor of 100 in the Last Glacial Maximum (LGM) compared to the Preboreal, and sharp variations of concentration occurred synchronously with rapid changes in the delta18O temperature proxy. The lognormal mode µ of the volume distribution shows clear systematic variations with smaller modes during warmer climates and coarser modes during colder periods. We find µ ~ 1.7 µm diameter during LGM and µ ~ 1.3 µm during the Preboreal. On timescales below several 100 years µ and the particle concentration exhibit a certain degree of independence present especially during warm periods, when µ generally is more variable. Using highly simplifying considerations for atmospheric transport and deposition of particles we infer that (1) the observed changes of µ in the ice largely reflect changes in the size of airborne particles above the ice sheet and (2) changes of µ are indicative of changes in long range atmospheric transport time. From the observed size changes we estimate shorter transit times by roughly 25% during LGM compared to the Preboreal. The associated particle concentration increase from more efficient long range transport is estimated to less than one order of magnitude.

Sedimentology of ODP Leg 119 samples

Cretaceous to Quaternary sediments recovered at Leg 119 Sites 738 and 744 on the southern tip of the Kerguelen Plateau were studied in order to determine the depositional environment and the paleoceanography of the southern Indian Ocean and especially the long-term glacial history of East Antarctica. Emphasis is given to bulk-sediment composition, grain-size data, and clay mineralogy. The sediment sequence at the two sites is generally of a highly pelagic character, with nannofossil oozes, chalks, and limestones dominant from the Turanian to upper Miocene and diatom oozes dominant within the uppermost Miocene to Holocene interval. The first indication of glaciation at sea level is the occurrence of isolated gravel and terrigenous sand grains, which indicate ice rafting in the middle Eocene interval of 45.0-42.3 Ma. A major intensification of glaciation, probably the onset of continental East Antarctic glaciation, is recorded in sediments of early Oligocene age (36.0 Ma). All major sediment parameters document this event. The clay mineralogy changes from smectite-dominated assemblages, typical of moderately warm and humid climatic conditions in which chemical weathering processes are prevalent, to illite- and chlorite-dominated assemblages, indicative of cooler climates and physical weathering. Ice-rafted debris of both gravel and sand size occurs in large quantities in that interval and coincides with a change in the mode of carbonate deposition. Carbonate contents are relatively high and uniform (90%-95%) in strata younger than early Oligocene; in Oligocene to upper Miocene strata they fluctuate between 65% and 95%. Oligocene and Neogene hiatuses reflect an intensification of oceanic circulation and the increased erosional force of Circumpolar Deep Water. The long-term Cenozoic cooling trend was interrupted by a phase of early Miocene warming indicated by maximum Neogene smectite concentrations. Although ice-rafted debris is present only in minor amounts and mainly in the silt fraction of early Oligocene to late Miocene age, it shows that glaciers advanced to the East Antarctic shoreline throughout that time. Ice-rafting activity drastically increased in latest Miocene time, when carbonate deposition decreased and diatom ooze sedimentation started. This suggests a pronounced intensification of Antarctic glaciation combined with a northward migration of the Polar Front.

(Table 1) Dissolved organic carbon, water temperature and conductivity in three subarctic ponds near Kilpisjärvi, North Finland

Daphnia was collected from five subarctic ponds which differed greatly in their DOC contents and, consequently, their underwater light (UV) climates. Irrespective of which Daphnia species was present, and contrary to expectations, the ponds with the lowest DOC concentrations (highest UV radiation levels) contained Daphnia with the highest eicosapentaenoic acid (EPA) concentrations. In addition, EPA concentrations in these Daphnia generally decreased in concert with seasonally increasing DOC concentrations. Daphnia from three of the ponds was also tested for its tolerance to solar ultraviolet radiation (UVR) with respect to survival. Daphnia pulex from the clear water pond showed, by far, the best UV-tolerance, followed by D. longispina from the moderately humic and D. longispina from the very humic pond. In addition, we measured sublethal parameters related to UV-damage such as the degree to which the gut of Daphnia appeared green (as a measure of their ability to digest algae), and whether their guts appeared damaged. We developed a simple, noninvasive scoring system to quantify the proportion of the gut in which digestive processes were presumably active. This method allowed repeated measurement of the same animals over the course of the experiment. We demonstrated, for the first time, that sublethal damage of the gut precedes mortality caused by exposure to UVR. In a parallel set of experiments we fed UV-exposed and non-exposed algae to UV-exposed and non-exposed daphnids. UVR pretreatment of algae enhanced the negative effects of exposure to natural solar UV-irradiation in Daphnia. These UV-related effects were generally not specific to the species of Daphnia.

Stable carbon and oxygen isotope record of foraminifera of IODP Site 303-1307

A major tipping point of Earth's history occurred during the mid-Pliocene: the onset of major Northern-Hemisphere Glaciation (NHG) and of pronounced, Quaternary-style cycles of glacial-to-interglacial climates, that contrast with more uniform climates over most of the preceding Cenozoic and continue until today (Zachos et al., 2001, doi:10.1126/science.1059412). The severe deterioration of climate occurred in three steps between 3.2 Ma (warm MIS K3) and 2.7 Ma (glacial MIS G6/4) (Lisiecki and Raymo, 2005, doi:10.1029/2004PA001071). Various models (sensu Driscoll and Haug, 1998, doi:10.1126/science.282.5388.436) and paleoceanographic records (intercalibrated using orbital age control) suggest clear linkages between the onset of NHG and the three steps in the final closure of the Central American Seaways (CAS), deduced from rising salinity differences between Caribbean and the East Pacific. Each closing event led to an enhanced North Atlantic meridional overturning circulation and this strengthened the poleward transport of salt and heat (warmings of +2-3°C) (Bartoli et al., 2005, doi:10.1016/j.epsl.2005.06.020). Also, the closing resulted in a slight rise in the poleward atmospheric moisture transport to northwestern Eurasia (Lunt et al., 2007, doi:10.1007/s00382-007-0265-6), which probably led to an enhanced precipitation and fluvial run-off, lower sea surface salinity (SSS), and an increased sea-ice cover in the Arctic Ocean, hence promoting albedo and the build-up of continental ice sheets. Most important, new evidence shows that the closing of the CAS led to greater steric height of the North Pacific and thus doubled the low-saline Arctic Throughflow from the Bering Strait to the East Greenland Current (EGC). Accordingly, Labrador Sea IODP Site 1307 displays an abrupt but irreversible EGC cooling of 6°C and freshening by ~2 psu from 3.25/3.16-3.00 Ma, right after the first but still reversible attempt of closing the CAS.

Temperature reconstruction of sediment core SO9-93KL

The Toba volcanic event, one of the largest eruptions during the Quaternary, is documented in marine sediment cores from the northeastern Arabian Sea. On the crest of the Murray Ridge and along the western Indian continental margin, we detected distinct concentration spikes and ash layers of rhyolithic volcanic shards near the marine isotope stage 5-4 boundary with the chemical composition of the "Youngest Toba Tuff". Time series of the Uk'37-alkenone index, planktic foraminiferal species, magnetic susceptibility, and sediment accumulation rates from this interval show that the Toba event occurred between two warm periods lasting a few millennia. Using Toba as an instantaneous stratigraphic marker for correlation between the marine- and ice-core chronostratigraphies, these two Arabian Sea climatic events correspond to Greenland interstadials 20 and 19, respectively. Our data sets thus depict substantial interstadial/stadial fluctuations in sea-surface temperature and surface-water productivity. We show that variable terrigenous (eolian) sediment supply played a crucial role in transferring and preserving the productivity signal in the sediment record. Within the provided stratigraphic resolution of several decades to centennials, none of these proxies shows a particular impact of the Toba eruption. However, our results are additional support that Toba, despite its exceptional magnitude, had only a minor impact on the evolution of low-latitude monsoonal climate on centennial to millennial time scales.

Age determination and grain size distribution of sediment core GeoB3375-1

The palaeoclimatic conditions during the Last Glacial Maximum (LGM) of southern South America and especially latitudinal shifts of the southern westerly wind belt are still discussed controversially. Longer palaeoclimatic records covering the Late Quaternary are rare. A particularly sensitive area to Late Quaternary climatic changes is the Norte Chico, northern Chile, because of its extreme climatic gradients. Small shifts of the present climatic zonation could cause significant variations of the terrestrial sedimentary environment which would be recorded in marine terrigenous sediments. To unveil the history of shifting climatic zones in northern Chile, we present a sedimentological study of a marine sediment core (GeoB 3375-1) from the continental slope off the Norte Chico (27.5°S). Sedimentological investigations include bulk- and silt grain-size determinations by sieving, Atterberg separation, and detailed SediGraph analyses. Additionally, clay mineralogical parameters were obtained by X-ray diffraction methods. The 14C-dated core, covering the time span from approximately 10,000 to 120,000 cal. yr B.P., consists of hemipelagic sediments. Terrigenous sedimentological parameters reveal a strong cyclicity, which is interpreted in terms of variations of sediment provenance, modifications of the terrestrial weathering regimes, and modes of sediment input to the ocean. These interpretations imply cyclic variations between comparatively arid climates and more humid conditions with seasonal precipitation for northern Chile (27.5°S) through the Late Quaternary. The cyclicity of the terrigenous sediment parameters is strongly dominated by precessional cycles. For the palaeoclimatic signal, this means that more humid conditions coincide with maxima of the precession index, as e.g. during the LGM. Higher seasonal precipitation for this part of Chile is most likely derived from frontal winter rain of the Southern Westerlies. Thus, the data presented here favour not only an equatorward shift of this atmospheric circulation system during the LGM, but also precession-controlled latitudinal movements throughout the Late Quaternary. Precessional forcing of latitudinal movements of the westerly atmospheric circulation system may be conceivable through teleconnections to the Northern Hemisphere monsoonal system in the Atlantic Ocean region.

Iron (Fe) concentrations and fluxes in TALDICE ice core 0.6 to 314 ky BP

Atmospheric fluxes of iron (Fe) over the past 200 kyr are reported for the coastal Antarctic Talos Dome ice core, based on acid leachable Fe concentrations. Fluxes of Fe to Talos Dome were consistently greater than those at Dome C, with the greatest difference observed during interglacial climates. We observe different Fe flux trends at Dome C and Talos Dome during the deglaciation and early Holocene, attributed to a combination of deglacial activation of dust sources local to Talos Dome and the reorganisation of atmospheric transport pathways with the retreat of the Ross Sea ice shelf. This supports similar findings based on dust particle sizes and fluxes and Rare Earth Element fluxes. We show that Ca and Fe should not be used as quantitative proxies for mineral dust, as they all demonstrate different deglacial trends at Talos Dome and Dome C. Considering that a 20 ppmv decrease in atmospheric CO2 at the coldest part of the last glacial maximum occurs contemporaneously with the period of greatest Fe and dust flux to Antarctica, we confirm that the maximum contribution of aeolian dust deposition to Southern Ocean sequestration of atmospheric CO2 is approximately 20 ppmv.

Northwards and eastwards velocities extracted from a historical run (for the year 2000) of the UVic Earth System Climate Model of Weaver et al. (2001)

This study provides a theoretical assessment of the potential bias due to differential lateral transport on multi-proxy studies based on a range of marine microfossils. Microfossils preserved in marine sediments are at the centre of numerous proxies for paleoenvironmental reconstructions. The precision of proxies is based on the assumption that they accurately represent the overlying watercolumn properties and faunas. Here we assess the possibility of a syn-depositional bias in sediment assemblages caused by horizontal drift in the water column, due to differential settling velocities of sedimenting particles based on their shape, size and density, and due to differences in current velocities. Specifically we calculate the post-mortem lateral transport undergone by planktic foraminifera and a range of other biological proxy carriers (diatoms, radiolaria and fecal pellets transporting coccolithophores) in several regions with high current velocities. We find that lateral transport of different planktic foraminiferal species is minimal due to high settling velocities. No significant shape- or size-dependent sorting occurs before reaching the sediment, making planktic foraminiferal ideal proxy carriers. In contrast, diatoms, radiolaria and fecal pellets can be transported up to 500km in some areas. For example in the Agulhas current, transport can lead to differences of up to 2°C in temperature reconstructions between different proxies in response to settling velocities. Therefore, sediment samples are likely to contain different proportions of local and imported particles, decreasing the precision of proxies based on these groups and the accuracy of the temperature reconstruction.

Planktonic foraminiferal oxygen isotopes of ODP Site 143-865 samples

Cool tropical sea surface temperatures (SSTs) are reported for warm Paleogene greenhouse climates based on the d18O of planktonic foraminiferal tests. These results are difficult to reconcile with models of greenhouse gas-forced climate. It has been suggested that this "cool tropics paradox" arises from postdepositional alteration of foraminiferal calcite, yielding erroneously high d18O values. Recrystallization of foraminiferal tests is cryptic and difficult to quantify, and the compilation of robust d18O records from moderately altered material remains challenging. Scanning electron microscopy of planktonic foraminiferal chamber-wall cross sections reveals that the basal area of muricae, pustular outgrowths on the chamber walls of species belonging to the genus Morozovella, contain no mural pores and may be less susceptible to postdepositional alteration. We analyzed the d18O in muricae bases of morozovellids from the central Pacific (Ocean Drilling Program Site 865) by ion microprobe using 10 ?m pits with an analytical reproducibility of ±0.34 per mil (2 standard deviations). In situ measurements of d18O in these domains yield consistently lower values than those published for conventional multispecimen analyses. Assuming that the original d18O is largely preserved in the basal areas of muricae, this new d18O record indicates Early Paleogene (~49-56 Ma) tropical SSTs in the central Pacific were 4°-8°C higher than inferred from the previously published d18O record and that SSTs reached at least ~33°C during the Paleocene-Eocene thermal maximum. This study demonstrates the utility of ion microprobe analysis for generating more reliable paleoclimate records from moderately altered foraminiferal tests preserved in deep-sea sediments.

Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Stable isotope record, Mg/Ca ratios and age model of St. Stephens Quarry, Alabama

In the largest global cooling event of the Cenozoic Era, between 33.8 and 33.5 Myr ago, warm, high-CO2 conditions gave way to the variable 'icehouse' climates that prevail today. Despite intense study, the history of cooling versus ice-sheet growth and sea-level fall reconstructed from oxygen isotope values in marine sediments at the transition has not been resolved. Here, we analyse oxygen isotopes and Mg/Ca ratios of benthic foraminifera, and integrate the results with the stratigraphic record of sea-level change across the Eocene-Oligocene transition from a continental-shelf site at Saint Stephens Quarry, Alabama. Comparisons with deep-sea (Sites 522 (South Atlantic) and 1218 (Pacific)) d18O and Mg/Ca records enable us to reconstruct temperature, ice-volume and sea-level changes across the climate transition. Our records show that the transition occurred in at least three distinct steps, with an increasing influence of ice volume on the oxygen isotope record as the transition progressed. By the early Oligocene, ice sheets were ~25% larger than present. This growth was associated with a relative sea-level decrease of approximately 105 m, which equates to a 67 m eustatic fall.