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.

Gridded climate data from 5 Global Climate Models (GCM) of the Last Glacial Maximum (LGM) downscaled to 30 arc seconds for Europe, with links to NetCDF files

Studies on the impact of historical, current and future global change require very high-resolution climate data (less or equal 1km) as a basis for modelled responses, meaning that data from digital climate models generally require substantial rescaling. Another shortcoming of available datasets on past climate is that the effects of sea level rise and fall are not considered. Without such information, the study of glacial refugia or early Holocene plant and animal migration are incomplete if not impossible. Sea level at the last glacial maximum (LGM) was approximately 125m lower, creating substantial additional terrestrial area for which no current baseline data exist. Here, we introduce the development of a novel, gridded climate dataset for LGM that is both very high resolution (1km) and extends to the LGM sea and land mask. We developed two methods to extend current terrestrial precipitation and temperature data to areas between the current and LGM coastlines. The absolute interpolation error is less than 1°C and 0.5 °C for 98.9% and 87.8% of all pixels for the first two 1 arc degree distance zones. We use the change factor method with these newly assembled baseline data to downscale five global circulation models of LGM climate to a resolution of 1km for Europe. As additional variables we calculate 19 'bioclimatic' variables, which are often used in climate change impact studies on biological diversity. The new LGM climate maps are well suited for analysing refugia and migration during Holocene warming following the LGM.

(Appendix table DR1) Foraminiferal distribution data for ODP Hole 178-1098B

Benthic foraminiferal data from Ocean Drilling Program Site 1098 indicate significant changes in deep-water conditions of the Palmer Deep, western Antarctic Peninsula margin, throughout the Holocene (13 ka to present). The earliest Holocene represents a period of transition from the Last Glacial Maximum (LGM). Cold bottom waters, similar to saline shelf water (SSW), dominated the middle Holocene. The late Holocene in the Palmer Deep has been characterized by alternating dominance of circumpolar deep water (CDW) and saline shelf water. These changes have global oceanographic and climatic implications. We suggest that the middle Holocene bottom-water record, in the absence of circumpolar deep water on the western Antarctic Peninsula shelf, indicates high saline shelf water production and/or weakened circumpolar deep water production during the middle Holocene climatic optimum. The late Holocene benthic foraminiferal record indicates rapidly fluctuating sea-ice conditions and may indicate a teleconnection between the South Pacific and Southern Ocean, thus having implications related to the Southern Oscillation Index.

Clay mineralogical assemblages, Sm-Nd concentrations, Sm and Nd isotopic ratios and oxygen isotope ratios, OIS number and the age model of ODP Site 105-646

We present 40 Sm-Nd isotope measurements of the clay-size (<2 µm) fractions of sediments from the Southern Greenland rise (ODP-646) that span the last 365 kyr. These data track changes in the relative supply of fine particles carried into the deep Labrador Sea by the Western Boundary Under Current (WBUC) back to the fourth glacial-interglacial cycles. Earlier studies revealed three general sources of particles to the core site: (i) Precambrian crustal material from Canada, Greenland, and/or Scandinavia (North American Shield - NAS), (ii) Palaeozoic or younger crustal material from East Greenland, NW Europe, and/or western Scandinavia (Young Crust - YC) and (iii) volcanic material from Iceland and the Mid-Atlantic Ridge (MAR). Clay-size fractions from glacial sediments have the lowest Nd isotopic ratios. Supplies of young crustal particles were similar during glacial oxygen isotope stages (OIS) 2, 6, and 10. In contrast the mean volcanic contributions decreased relative to old craton material from OIS 10 to OIS 6 and then from OIS 6 to OIS 2. The glacial OIS 8 interval displays a mean Sm/Nd ratio similar to those of interglacials OIS 1, 5, and 9. Compared with other interglacials, OIS 7 was marked by a higher YC contribution but a similar ~30% MAR supply. The overall NAS contribution dropped by a factor of 2 during each glacial/interglacial transition, with the MAR contribution broadly replacing it during interglacials. To decipher between higher supplies and/or dilution, particle fluxes from each end member were estimated. Glacial NAS fluxes were systematically higher than interglacial fluxes. During the time interval examined, fine particle supplies to the Labrador Sea were strongly controlled by proximal ice-margin erosion and thus echoed the glacial stage intensity. In contrast, the WBUC-carried MAR supplies from the eastern basins did not change significantly throughout the last 365 kyr, except for a marked increase in surface-sediments that suggests unique modern conditions. Distal WBUC-controlled inputs from the Northern and NE North Atlantic seem to have been less variable than proximal supplies linked with glacial erosion rate.

Silt grain-size analysis of ODP holes from the western equatorial Atlantic

Grain-size records of the terrigenous and calcareous silt fraction, preservation of planktic foraminifera, and benthic foraminiferal stable-isotope data (delta13C, delta18O values of C. wuellerstorfi) at ODP Site 927 on the Ceará Rise (5°27.7'N, 44°28.8'W), are used to reconstruct variations in the history of bottom current strength, ventilation, and carbonate corrosiveness of deep waters during the time interval from 0.8 to 0.3 Ma. Glacial periods are characterized by generally smaller mean sizes of the terrigenous sortable silt fraction (mean(SS)), lower delta13C values, and poorer preservation of planktic foraminifera compared to interglacials. This indicates lower bottom current speeds, larger nutrient contents and more corrosive deep water. By contrast, larger mean(SS) sizes, higher delta13C values, and well preserved planktic foraminifera indicate strong circulation and a well ventilated deep-water mass during interglacials. The observed changes are most likely related to the weakening and strengthening of circulation of Lower North Atlantic Deep Water (LNADW). Cross-spectral analysis between the mean(SS) and benthic delta18O records reveals that minima in mean(SS) occur about 7.6 k.y. after the maximum in ice volume. This indicates a considerable lag time between ice-shield induced changes in LNADW production and subsequent changes in the velocity of LNADW flow in the western equatorial Atlantic. Striking changes in bottom current speed occur regularly during glacial to interglacial transitions. Extremely fine mean(SS) minima point to an almost complete shutdown of bottom current vigor in response to a cessation of LNADW production caused by an enhanced melt water release during the initial phases of deglaciation. However, each of the fine minima extremes is followed by a rapid shift to very high mean(SS) values that indicate strong bottom currents, and hence, vigorous LNADW flow during the early interglacials. After the onset of glacial Stage 12, generally poorer carbonate preservation and higher variability is registered. This coincides with a global decrease in carbonate preservation during the mid-Brunhes (mid-Brunhes dissolution event). Detailed grain-size analysis of the calcareous fine fraction (<63 µm) revealed a considerable reduction of particles in the fraction from 7 to 63 µm during periods of enhanced dissolution. This indicates a preferential dissolution of larger planktic foraminiferal fragments which leads to an enrichment of coccoliths in the calcareous fine fraction.

Beryllium-10 and aluminium-26 concentration and production rates in glacial deposits from the Quartermain Mountains

Many glacial deposits in the Quartermain Mountains, Antarctica present two apparent contradictions regarding the degradation of unconsolidated deposits. The glacial deposits are up to millions of years old, yet they have maintained their meter-scale morphology despite the fact that bedrock and regolith erosion rates in the Quartermain Mountains have been measured at 0.1-4.0 m/Ma. Additionally, ground ice persists in some Miocene-aged soils in the Quartermain Mountains even though modeled and measured sublimation rates of ice in Antarctic soils suggest that without any recharge mechanisms ground ice should sublimate in the upper few meters of soil on the order of 10**3 to 10**5 years. This paper presents results from using the concentration of cosmogenic nuclides beryllium-10 (10Be) and aluminum-26 (26Al) in bulk sediment samples from depth profiles of three glacial deposits in the Quartermain Mountains. The measured nuclide concentrations are lower than expected for the known ages of the deposits, erosion alone does not always explain these concentrations, and deflation of the tills by the sublimation of ice coupled with erosion of the overlying till can explain some of the nuclide concentration profiles. The degradation rates that best match the data range 0.7-12 m/Ma for sublimation of ice with initial debris concentrations ranging 12-45% and erosion of the overlying till at rates of 0.4-1.2 m/Ma. Overturning of the tills by cryoturbation, vertical mixing, or soil creep is not indicated by the cosmogenic nuclide profiles, and degradation appears to be limited to within a few centimeters of the surface. Erosion of these tills without vertical mixing may partially explain how some glacial deposits in the Quartermain Mountains maintain their morphology and contain ground ice close to the surface for millions of years.

Paleoenvironmental changes in the western Mediterranean since the last glacial maximum: High resolution multiproxy record from the Algero-Balearic basin

The present study uses a multiproxy approach in order to further understand the evolution of climate responses in the western Mediterranean as of the Last Glacial Maximum. Sediments from ODP Site 975 in the Algero-Balearic basin have been analysed at high resolution, both geochemically andmineralogicallly. The resulting data have been used as proxies to establish a sedimentary regime, primary marine productivity, the preservation of the proxies and oxygen conditions. Fluctuations in detrital element concentrations were mainly the consequence of wet/arid oscillations. Productivity has been established using Ba excess, according to which marine productivity appears to have been greatest during cold events Heinrich 1 and Younger Dryas. The S1 time interval was not as marked by increases in productivity as was the eastern Mediterranean. In contrast, the S1 interval was first characterized by a decreasing trend and then by a fall in productivity after the 8.2 ky BP dry-cold event. Since then productivity has remained low. Here we report that there was an important redox event in this basin, probably a consequence of the major oceanographic circulation change occurring in the western Mediterranean at 7.7 ky BP. This circulation change led to reventilation as well as to diagenetic remobilization of redox-sensitive elements and organic matter oxidation. Comparisons between our paleoceanographic reconstruction for this basin and those regarding other Mediterranean basins support the hypothesis that across the Mediterranean there were different types of responses to climate forcing mechanism. The Algero-Balearic basin is likely to be a key area for further understanding of the relationships between the North Atlantic and the eastern Mediterranean basins.

Late Quaternary siliciclastic/carbonate sedimentation in the Cabricorn Channel

A model is presented for hemipelagic siliciclastic and carbonate sedimentation during the last glacial-interglacial cycle in the Capricorn Channel, southern Great Barrier Reef (GBR). Stable isotope ratios, grainsize, carbonate content and mineralogy were analysed for seven cores in a depth transect from 166 to 2892 m below sea level (mbsl). Results show variations in the flux of terrigenous, neritic and pelagic sediments to the continental slope over the last sea level cycle. During the glacial lowstand terrigenous sediment influenced all the cores down to 2000 mbsl. The percentages of quartz and feldspar in the cores decreased with water depth, while the percentage of clay increased. X-ray diffraction analysis of the glacial lowstand clay mineralogy suggests that the siliciclastic sediment was primarily sourced from the Fitzroy River, which debouched directly into the northwest sector of the Capricorn Channel at this time. The cores also show a decrease in pelagic calcite and an increase in aragonite and high magnesium calcite (HMC) during the glacial. The influx of HMC and aragonite is most likely from reworking of coral reefs exposed on the continental shelf during the glacial, and also from HMC ooids precipitated at the head of the Capricorn Channel at this time. Mass accumulation rates (MARs) are high (13.5 g/cm**/kyr) during the glacial and peak at ~20 g/cm** 3/kyr in the early transgression (16-14 ka BP). MARs then decline with further sea level rise as the Fitzroy River mouth retreats from the edge of the continental shelf after 13.5 ka BP. MARs remain low (4 g/cm**3/kyr) throughout the Holocene highstand. Data for the Holocene highstand indicate there is a reduction in siliciclastic influx to the Capricorn Channel with little quartz and feldspar below 350 mbsl. However, fine-grained fluvial sediments, presumably from the Fitzroy River, were still accumulating on the mid slope down to 2000 mbsl. The proportion of pelagic calcite in the core tops increases with water depth, while HMC decreases, and is present only in trace amounts in cores below 1500 mbsl. The difference in the percentage of HMC in the deeper cores between the glacial and Holocene may reflect differences in supply or deepening of the HMC lysocline during the glacial. Sediment accumulation rates also vary between cores in the Capricorn Channel and do not show the expected exponential decrease with depth. This may be due to intermediate or deep water currents reworking the sediments. It is also possible that present bathymetry data are too sparse to detect the potential role that submarine channels may play in the distribution and accumulation of sediments. Comparison of the Capricorn Channel MARs with those for other mixed carbonate/siliciclastic provinces from the northeast margin of Australia indicates that peak MARs in the early transgression in the Capricorn Channel precede those from the central GBR and south of Fraser Island. The difference in the timing of the carbonate and siliciclastic MAR peaks along the northeast margin is primarily related to differences in the physiography and climate of the provinces. The only common trend in the MARs from the northeast margin of Australia is the near synchronicity of the carbonate and siliciclastic MAR peaks in individual sediment cores, which supports a coeval sedimentation model.

(Table 3) Aeolian component of North Pacific pelagic sediments

Isolation and analysis of the eolian component of late Cenozoic pelagic sediments from the North Pacific provides direct information concerning changes in atmospheric circulation. A 50% increase in intensity of both the prevailing westerlies and the tradewinds coincides with increasing pole-to-equator temperature gradients resulting from the onset of northern hemisphere glaciation. At the same time, the mass flux of dust from continents to the North Pacific increased by a factor of 4.5, apparently reflecting significantly increased continental aridity associated with the late Cenozoic glacial ages.

Sedimentation rates and foraminiferal assemblages from ODP Site 162-184

Sedimentological and faunal records from the transitional period marking the onset of widespread northern hemisphere glaciation have been investigated at Ocean Drilling Program Site 984. The late Pliocene interglacial sediments of the northeast Atlantic are carbonate rich and show evidence of vigorous bottom water circulation at intermediate water depths. Contrasting this, the late Pliocene glacial sediments are characterised by carbonate dissolution and slower bottom current velocities. Weak or "leaky" Norwegian Sea overflows, undersaturated with respect to carbonate, influenced this region during the late Pliocene glacials. The early Pleistocene pattern of intermediate water circulation appears to have changed radically in the northeast Atlantic. At this time, interglacial carbonate values and inferred bottom current velocities are low. This suggests slow-flowing, undersaturated Norwegian Sea water bathing the site. The overflow increased during the early Pleistocene interglacials as the exchange between the Atlantic and Norwegian-Greenland Seas improved. The most significant feature of the early Pleistocene glacials is the increase in inferred bottom current velocity. These changes reflect a switch in deep North Atlantic convection to shallower depths during glacial periods, possibly in a manner similar to the increasing contribution of glacial intermediate water to the North Atlantic during the late Pleistocene glacials. Our results suggest that the late Pleistocene climate variability of the North Atlantic is a pervasive feature of the late Pliocene-early Pleistocene record.

Sea surface temperature calculated from UK'37 and IP25 sea-ice biomarker of sediment core SO201-2-85KL

During the past decades, remarkable changes in sea-surface temperature (SST) and sea-ice extent have been observed in the marginal seas of the subarctic Pacific. However, little is known about natural climate variability at millennial time scales far beyond instrumental observations. Geological proxy records, such as those derived from marine sediments, offer a unique opportunity to investigate millennial-scale natural climate variability of the Artic and subarctic environments during past glacial-interglacial cycles. Here we provide reconstructions of sea-ice variability inferred from IP25 (Ice Proxy with 25 carbon atoms) sea-ice biomarker and SST fluctuations based on alkenone unsaturation index (UK'37) of the subarctic Pacific realm between 138 and 70 ka. Warmest sea-surface conditions were found during the early Eemian interglacial (128 to 126 ka), exceeding modern SSTs by ~2 °C. The further North Pacific climate evolu- tion is marked by pronounced oscillations in SST and sea-ice extent on millennial time scales, which correspond remarkably well to short-term temperature oscillations known from Green- land and the North Atlantic. These results imply a common forcing, which seems to be closely coupled to dynamics of the Atlantic meridional overturning circulation. However, immediate propagation of such climate fluctuations far beyond the North Atlantic basin suggests a rapid circumpolar coupling mechanism probably acting through the atmosphere, a prerequisite to explain the apparent synchronicity of remote climatic reorganizations in the subarctic Pacific.

Glacial-interglacial sediment record from the North Atlantic

Studies from the subtropical western and eastern Atlantic Ocean, using the 231Pa/230Th ratio as a kinematic proxy for deep water circulation, provided compelling evidence for a strong link between climate and the rate of meridional overturning circulation (MOC) over the last deglaciation. In this study, we present a compilation of existing and new sedimentary 231Pa/230Th records from North Atlantic cores between 1710 and 4550 m water depth. Comparing sedimentary 231Pa/230Th from different depths provides new insights into the evolution of the geometry and rate of deep water formation in the North Atlantic during the last 20,000 years. The 231Pa/230Th ratio measured in upper Holocene sediments indicates slow water renewal above ?2500 m and rapid flushing below, consistent with our understanding of modern circulation. In contrast, during the Last Glacial Maximum (LGM), Glacial North Atlantic Intermediate Water (GNAIW) drove a rapid overturning circulation to a depth of at least ?3000 m depth. Below ~4000 m, water renewal was much slower than today. At the onset of Heinrich event 1, transport by the overturning circulation declined at all depths. GNAIW shoaled above 3000 m and significantly weakened but did not totally shut down. During the Bølling-Allerød (BA) that followed, water renewal rates further decreased above 2000 m but increased below. Our results suggest for the first time that ocean circulation during that period was quite distinct from the modern circulation mode, with a comparatively higher renewal rate above 3000 m and a lower renewal rate below in a pattern similar to the LGM but less accentuated. MOC during the Younger Dryas appears very similar to BA down to 2000 m and slightly slower below.