LPJ-LMfire simulations for Europe for the Last Glacial Maximum and preindustrial control, link to netCDF files

This data sets contains LPJ-LMfire dynamic global vegetation model output covering Europe and the Mediterranean for the Last Glacial Maximum (LGM; 21 ka) and for a preindustrial control simulation (20th century detrended climate). The netCDF data files are time averages of the final 30 years of the model simulation. Each netCDF file contains four or five variables: fractional cover of 9 plant functional types (PFTs; cover), total fractional coverage of trees (treecover), population density of hunter-gatherers (foragerPD; only for the "people" simulations), fraction of the gridcell burned on 30-year average (burnedf), and vegetation net primary productivity (NPP). The model spatial resolution is 0.5-degrees For the LGM simulations, LPJ-LMfire was driven by the PMIP3 suite of eight GCMs for which LGM climate simulations were available. Also provided in this archive is the result of an LPJ-LMfire run that was forced by the average climate of all GCMs (the "GCM-mean" files), and the average of each of the individual LPJ-LMfire runs over the eight LGM scenarios individually (the "LPJ-mean" files). The model simulations are provided that include the influence of human presence on the landscape (the "people" files), and in a "world without humans" scenario (the "natural" files). Finally this archive contains the preindustrial reference simulation with and without human influence ("PI_reference_people" and "PI_reference_nat", respectively). There are therefore 22 netCDF files in this archive: 8 each of LGM simulations with and without people (total 16) and the "GCM mean" simulation (2 files) and the "LPJ mean" aggregate (2 files), and finally the two preindustrial "control" simulations ("PI"), with and without humans (2 files). In addition to the LPJ-LMfire model output (netCDF files), this archive also contains a table of arboreal pollen percent calculated from pollen samples dated to the LGM at sites throughout (lgmAP.txt), and a table containing the location of archaeological sites dated to the LGM (LGM_archaeological_site_locations.txt).

Last glacial maximum sediment record of the southern Weddell Sea shelf

Abstract: The history of grounded ice-sheet extent on the southern Weddell Sea shelf during the Last Glacial Maximum (LGM) and the timing of post-LGM ice-sheet retreat are poorly constrained. Several glaciological models reconstructed widespread grounding and major thickening of the Antarctic Ice Sheet in the Weddell Sea sector at the LGM. In contrast, recently published onshore data and modelling results concluded only very limited LGM-thickening of glaciers and ice streams feeding into the modern Filchner and Ronne ice shelves. These studies concluded that during the LGM ice shelves rather than grounded ice covered the Filchner and Ronne troughs, two deep palaeo-ice stream troughs eroded into the southern Weddell Sea shelf. Here we review previously published and unpublished marine geophysical and geological data from the southern Weddell Sea shelf. The stratigraphy and geometry of reflectors in acoustic sub-bottom profiles are similar to those from other West Antarctic palaeo-ice stream troughs, where grounded ice had advanced to the shelf break at the LGM. Numerous cores from the southern Weddell Sea shelf recovered sequences with properties typical for subglacially deposited tills or subglacially compacted sediments. These data sets give evidence that grounded ice had advanced across the shelf during the past, thereby grounding in even the deepest parts of the Filchner and Ronne troughs. Radiocarbon dates from glaciomarine sediments overlying the subglacial deposits are limited, but indicate that the ice grounding occurred at the LGM and that ice retreat started before ~15.1 corrected 14C kyrs before present (BP) on the outer shelf and before ~7.7 corrected 14C kyrs BP on the inner shelf, which is broadly synchronous with ice retreat in other Antarctic sectors. The apparent mismatch between the ice-sheet reconstructions from marine and terrestrial data can be attributed to ice streams with very low surface profiles (similar to those of "ice plains") that had advanced through Filchner Trough and Ronne Trough at the LGM. Considering the global sea-level lowstand of ~130 metres below present, a low surface slope of the expanded LGM-ice sheet in the southern Weddell Sea can reconcile grounding-line advance to the shelf break with limited thickening of glaciers and ice streams in the hinterland. This scenario implies that ice-sheet growth in the Weddell Sea sector during the LGM and ice-sheet drawdown throughout the last deglaciation could only have made minor contributions to the major global sea-level fluctuations during these times.

Sedimentary record of the western Amundsen Sea Embayment after the last glacial maximum

Reliable dating of glaciomarine sediments deposited on the Antarctic shelf since the Last Glacial Maximum (LGM) is very challenging because of the general absence of calcareous (micro-) fossils and the recycling of fossil organic matter. As a consequence, radiocarbon (14C) ages of the acid-insoluble organic fraction (AIO) of the sediments bear uncertainties that are very difficult to quantify. In this paper we present the results of three different chronostratigraphic methods to date a sedimentary unit consisting of diatomaceous ooze and diatomaceous mud that was deposited following the last deglaciation at five core sites on the inner shelf in the western Amundsen Sea (West Antarctica). In three cores conventional 14C dating of the AIO in bulk sediment samples yielded age reversals down-core, but at all sites the AIO 14C ages obtained from diatomaceous ooze within the diatom-rich unit yielded similar uncorrected 14C ages ranging from 13,517±56 to 11,543±47 years before present (yr BP). Correction of these ages by subtracting the core-top ages, which are assumed to reflect present-day deposition (as indicated by 21044 Pb dating of the sediment surface at one core site), yielded ages between ca. 10,500 and 8,400 calibrated years before present (cal yr BP). Correction of the AIO ages of the diatomaceous ooze by only subtracting the marine reservoir effect (MRE) of 1,300 years indicated deposition of the diatom-rich sediments between 14,100 and 11,900 cal yr BP. Most of these ages are consistent with age constraints between 13.0 and 8.0 ka BP for the diatom-rich unit, which we obtained by correlating the relative palaeomagnetic intensity (RPI) records of three of the sediment cores with global and regional reference curves for palaeomagnetic intensity. As a third dating technique we applied conventional 53 radiocarbon dating of the AIO included in acid-cleaned diatom hard parts that were extracted from the diatomaceous ooze. This method yielded uncorrected 14C ages of only 5,111±38 and 5,106±38 yr BP, respectively. We reject these young ages, because they are likely to be overprinted by the adsorption of modern atmospheric carbon dioxide onto the surfaces of the extracted diatom hard parts prior to sample graphitisation and combustion for 14C dating. The deposition of the diatom-rich unit in the western Amundsen Sea suggests deglaciation of the inner shelf before ca. 13 ka BP. The deposition of diatomaceous oozes on other parts of the Antarctic shelf around the same time, however, seems to be coincidental rather than directly related.

Sea-surface temperature and sea ice distribution of the Southern Ocean at the EPILOG Last Glacial Maximum

Based on the quantitative study of diatoms and radiolarians, summer sea-surface temperature (SSST) and sea ice distribution were estimated from 122 sediment core localities in the Atlantic, Indian and Pacific sectors of the Southern Ocean to reconstruct the last glacial environment at the EPILOG (19.5-16.0 ka or 23 000-19 000 cal yr. B.P.) time-slice. The statistical methods applied include the Imbrie and Kipp Method, the Modern Analog Technique and the General Additive Model. Summer SSTs reveal greater surface-water cooling than reconstructed by CLIMAP (Geol. Soc. Am. Map Chart. Ser. MC-36 (1981) 1), reaching a maximum (4-5 °C) in the present Subantarctic Zone of the Atlantic and Indian sector. The reconstruction of maximum winter sea ice (WSI) extent is in accordance with CLIMAP, showing an expansion of the WSI field by around 100% compared to the present. Although only limited information is available, the data clearly show that CLIMAP strongly overestimated the glacial summer sea ice extent. As a result of the northward expansion of Antarctic cold waters by 5-10° in latitude and a relatively small displacement of the Subtropical Front, thermal gradients were steepened during the last glacial in the northern zone of the Southern Ocean. Such reconstruction may, however, be inapposite for the Pacific sector. The few data available indicate reduced cooling in the southern Pacific and give suggestion for a non-uniform cooling of the glacial Southern Ocean.

Holocene and Last Glacial Maximum Atlantic neodymium isotopes

Neodymium isotopes measured on chemically uncleaned planktic foraminifera from cores throughout the Atlantic Ocean. Samples are Holocene and Last Glacial Maximum in age.

Tab. 1+7: Hydrochemical analyses of firn and meltwater from different "cold" snow patches

In wide areas of Northern Siberia, glaciers have been absent since the Late Pleistocene. Therefore, ground ice and especially ice wedges are used as archives for paleoclimatic studies. In the present study, carried out on the Bykovsky Peninsula, eastern Lena Delta, we were able to distinguish ice wedges of different genetic units by means of oxygen and hydrogen isotopes. The results obtained by this study on the Ice Complex, a peculiar periglacial phenomenon, allowed the reconstruction of the climate history with a subdivision of a period of very cold winters (60-55 ka), followed by a long stable period of cold winter temperatures (50-24 ka), Between 20 ka and 11 ka, climate warming is indicated in stable isotope compositions, most probably after the Late Glacial Maximum. At that time, a change of the marine source of the precipitation from a more humid source to the present North AtIantic source region was assumed. For the Ice Complex, a continuous age-height relationship was established, indicating syngenetic vertical ice wedge growth and sediment accumulation rates of 0.7 m/ky. During the Holocene optimum, ice wedge growth was probably limited due to the extensive formation of lacustrine environments. Holocene ice wedges in thermokarst depressions (alases) and thermoerosional valleys (logs) were formed after climate deterioration from about 4.5 ka until the present. Winter temperatures were warmer at this time as compared to the cooler Pleistocene. Migration of bound water between ice wedges and segregated ice may have altered the isotopic composition of old ice wedges. The presence of ice wedges as diagnostic features for permafrost conditions since 60 ka, implies that a large glacier extending over the Laptev Sea shelf did not exist. For the remote non-glaciated areas of Northern Siberia, ice wedges were established as a powerful climate archive.

Sr/Ca ratios of Cenozoic benthic foraminifera

A Cenozoic multi-species record of benthic foraminiferal calcite Sr/Ca has been produced and is corrected for interspecific offsets (typically less than 0.3 mmol/mol) and for the linear relationship between decreasing benthic foraminiferal Sr/Ca and increasing water depth. The water depth correction, determined from Holocene, Late Glacial Maximum and Eocene paleowater-depth transects, is ~0.1 mmol/mol/km. The corrected Cenozoic benthic foraminiferal Sr/Ca record ranges from 1.2 to 2.0 mmol/mol, and has been interpreted in terms of long-term changes in seawater Sr/Ca, enabling issues related to higher-resolution variability in Sr/Ca to be ignored. We estimate that seawater Sr/Ca was ~1.5 times modern values in the late Cretaceous, but declined rapidly into the Paleogene. Following a minimum in the Eocene, seawater Sr/Ca increased gradually through to the present day with a minimum superimposed on this trend centered in the late Miocene. By assuming scenarios for changing seawater calcium concentration, and using published carbonate accumulation rate data combined with suitable values for Sr partition coefficients into carbonates, the seawater Sr/Ca record is used to estimate global average river Sr fluxes. These fluxes are used in conjunction with the seawater strontium isotope curve and estimates of hydrothermal activity/tectonic outgassing to calculate changes in global average river 87Sr/86Sr through the Cenozoic. The absolute magnitude of Sr fluxes and isotopic compositions calculated in this way are subject to relatively large uncertainties. Nevertheless, our results suggest that river Sr flux increased from 35 Ma to the present day (roughly two-fold) accompanied by an overall increase in 87Sr/86Sr (by ~0 to 0.001). Between 75 and 35 Ma, river 87Sr/86Sr also increased (by ~0.001 to 0.002) but was accompanied by a decrease (two- to three-fold) in river Sr flux.

Hydrogen and oxygen isotopes in ground lce of Bykovsky Peninsula, North Siberia

In wide areas of Northern Siberia, glaciers have been absent since the Late Pleistocene. Therefore, ground ice and especially ice wedges are used as archives for paleoclimatic studies. In the present study, carried out on the Bykovsky Peninsula, eastern Lena Delta, we were able to distinguish ice wedges of different genetic units by means of oxygen and hydrogen isotopes. The results obtained by this study on the Ice Complex, a peculiar periglacial phenomenon, allowed the reconstruction of the climate history with a subdivision of a period of very cold winters (60-55 ka), followed by a long stable period of cold winter temperatures (50-24 ka), Between 20 ka and I I ka, climate warming is indicated in stable isotope compositions, most probably after the Late Glacial Maximum. At that time, a change of the marine source of the precipitation from a more humid source to the present North Atlantic source region was assumed. For the Ice Complex, a continuous age-height relationship was established, indicating syngenetic vertical ice wedge growth and sediment accumulation rates of 0.7 m/ky. During the Holocene optimum, ice wedge growth was probably limited due to the extensive formation of lacustrine environments. Holocene ice wedges in thermokarst depressions (alases) and thermoerosional valleys (logs) were formed after climate deterioration from about 4.5 ka until the present. Winter temperatures were warmer at this time as compared to the cooler Pleistocene. Migration of bound water between ice wedges and segregated ice may have altered the isotopic composition of old ice wedges. The presence of ice wedges as diagnostic features for permafrost conditions since 60 ka, implies that a large glacier extending over the Laptev Sea shelf did not exist. For the remote non-glaciated areas of Northern Siberia, ice wedges were established as a powerful climate archive.

Stable isotope record and planktonic foraminiferal fragmentation of the late Pleistocene section of ODP Hole 134-828A

A high-resolution record of foraminiferal fragmentation (a dissolution indicator) for the last 250 k.y. (isotopic Stages 1 to 7) is identified in the upper 61.9 m of Ocean Drilling Program (ODP) Hole 828A, west Vanuatu. This record is comparable in detail to the atmospheric CO2 record and the d18O stack. Phase shifts between preservation spikes and maximum ice volumes (d18O of Globigerinoides sacculifer) are analogous to those on Ontong Java Plateau. Mass spectrometer (AMS14C) dating of a sample taken at the base of dissolution cycle B1 and the position of the last glacial maximum indicates a lag in time of ~8 k.y. in the Vanuatu region for the last glacial termination. When dissolution spikes are compared with minimum ice volumes there is no phase shift for the last two glacial terminations. The difference between Vanuatu and Ontong Java Plateau may be explained by local CO2 sinks and the interplay between intermediate and deep water masses. Terrigenous input increasingly affected sediment of Hole 828A on the North d'Entrecasteaux Ridge (NDR) as it approached Espiritu Santo Island. Mud and silt suspended in mid-water flows become important after 125 ka, while turbidites bypass the New Hebrides Trench only towards the last glacial maximum (LGM). Terrigenous supply seems to affect the lysocline profile that changed from an "open ocean" to a "near continent" type, thus favoring dissolution. Fragmentation of planktonic foraminifers is a more sensitive indicator of lysocline variations than is foraminiferal susceptibility to dissolution, the foraminiferal dissolution index, the abundance of benthic foraminifers, or CaCO3 content. A modern foraminiferal lysocline for the neighboring area (between 10°S and 30°S, and 160°E and 180°E) is found at 3.1 km below sea level, compared to west Vanuatu where it is shallower. The past lysocline level was deeper than 3086 m during intervals of dissolution minima, and ranged from ~2550 to 3000 m during intervals of dissolution maxima. The high sedimentation rates (in the order of 10 to 50 cm/k.y.) found in Hole 828A offer a great potential for future high-resolution studies either in this hole or other western localities along the NDR. Areas of high sedimentation near continental regions have been discarded for paleoceanographic and/or paleoclimatic studies. Nonetheless, conditions analogous to those found in Hole 828A are expected to occur in many trench areas around the world where mid-water flows have preserved as yet undiscovered fine high-resolution sedimentary records.

Response of deep-sea benthic foraminifera to Late Quaternary climate changes, southeast Indian Ocean, offshore Western Australia

The Late Quaternary benthic foraminifera of four deep-sea cores off Western Australia (ODP 122-760A, ODP 122-762B, BMR96GC21 and RC9-150) have been examined for evidence of increased surface productivity to explain the anomalously low sea-surface paleotemperatures inferred by planktic foraminifera for the last and penultimate glaciations. The delta13C trends of Cibicidoides wuellerstorfi, and differences between the delta13C trends of planktics (Globigerinoides sacculifer) and benthics (C. wuellerstorfi) in the four cores indicate that during stage 6 bottom waters were significantly depleted in delta13C, and strong delta13C gradients were established in the water column, while during stage 2 and the Last Glacial Maximum, delta13C trends did not differ greatly from that of the Holocene. Two main assemblages of benthic foraminifera were identified by principal component analyses: one dominated by Uvigerina peregrina, another dominated by U. proboscidea. Abundance of these Uvigerinids, and of taxa preferring an infaunal microhabitat, and of Epistominella exigua and Bulimina aculeata indicate that episodes of high influx of particulate organic matter were established in most sites during glacial episodes, and particularly so during stage 6, while evidence for upwelling during the Last Glacial Maximum is less strong. The Penultimate Glaciation upwellings were established within the areas of low sea-surface paleotemperature indicated by planktic foraminifera. During the Last Interglacial Climax, upwelling appears to have been established in an isolated region offshore from a strengthened Leeuwin Current off North West Cape. Last Glacial Maximum delta13C values of C. wuellerstorfi at waterdepths of less than 2000 m show smaller than global mean glacial-interglacial changes suggesting the development of a deep hydrological front. A similar vertical stratification/bathyal front was also established during the Penultimate Glaciation.

Late Quateranry dinoflagellate cysts of ODP Site 202-1233

The late Quaternary organic-walled dinoflagellate cyst record of Site 1233 (41°S, offshore Chile) was studied with a ?200 year resolution spanning the last 25,000 years. The study provides the first continuous record of sub-recent and recent dinoflagellate cysts in the Southeast (SE) Pacific. Major changes in the composition of the cyst association, cyst concentration and morphology of Operculodinium centrocarpum reflect changes in sea surface temperature (SST), sea surface salinity (SSS), palaeoproductivity and upwelling intensity. These changes can be associated with latitudinal shifts of the circumpolar frontal systems. The high cyst concentration, high Brigantedinium spp. abundances, low species diversity and the occurrence of certain cold water species are supportive for a 7-10° equatorward shift of the Antarctic Circumpolar Current (ACC) during the coldest phase of the Last Glacial Maximum (LGM) between 25 and 21.1 cal ka BP. Deglacial warming initiated at ~18.6 cal ka BP. Termination I (18.6-11.1 cal ka BP) is interrupted by an unstable period of extreme seasonality, rather than a cooling event, between 14.4 and 13.2 cal ka BP, synchronous with the Antarctic Cold Reversal (ACR). The Holocene Maximum is observed between 11.6 and 9.8 cal ka BP and is typified by the most southward position of the northern margin of the ACC. A cooling phase occurred during the early Holocene (until ~7 cal ka BP) and during the last ~0.8 ka. Our data indicates that the SE Pacific (41°S) climate has been influenced over the whole record by changes in the Southern Hemisphere (SH) high-latitudes, while during the mid to late Holocene, also a tropical forcing mechanism was involved, including the El Niño Southern Oscillation and the variable Hadley cell intensity. Furthermore, this study showed a relationship between the variable morphology of the spines/processes of O. centrocarpum and the combined variation of sea surface salinity and temperature (SSS/SST-ratio).

Glacial-interglacial record of ODP Site 181-1119 in the Canterbury Bright, Southwest Pacific Ocean

Ocean Drilling Program (ODP) Site 1119 is located at water depth 395 m near the subtropical front (STF; here represented by the Southland Front), just downslope from the shelf edge of eastern South Island, New Zealand. The upper 86.19 metres composite depth (mcd) of Site 1119 sediment was deposited at an average sedimentation rate of 34 cm/kyr during Marine Isotope Stages (MIS) 1-8 (0-252 ka), and is underlain across a ~25 kyr intra-MIS 8 unconformity by MIS 8.5-11 (277-367 ka) and older sediment deposited at ~14 cm/kyr. A time scale is assigned to Site 1119 using radiocarbon dates for the period back to ~39 ka, and, prior to then, by matching its climatic record with that of the Vostok ice core, which it closely resembles. Four palaeoceanographic proxy measures for surface water masses vary together with the sandy-muddy, glacial-interglacial (G/I) cyclicity at the site. Interglacial intervals are characterised by heavy delta13C, high colour reflectance (a proxy for carbonate content), low Q-ray (a proxy for clay content) and light delta18O; conversely, glacial intervals exhibit light delta13C, low reflectance, high Q-ray and heavy delta18O signatures. Early interglacial intervals are represented by silty clays with 10-105-cm-thick beds of sharp-based (Chondrites-burrowed), shelly, graded, fine sand. The sands are rich in foraminifera, and were deposited distant from the shoreline under the influence of longitudinal flow in relatively deep water. Glacial intervals comprise mostly micaceous silty clay, though with some thin (2-10 cm thick) sands present also at peak cold periods, and contain the cold-water scallop Zygochlamys delicatula. Interglacial sandy intervals are characterised by relatively low sedimentation rates of 5-32 cm/kyr; cold climate intervals MIS 10, 6 and 2 have successively higher sedimentation rates of 45, 69 and 140 cm/kyr. Counter-intuitively,and forced by the bathymetric control of a laterally-moving shoreline during G/I and I/G transitions, the 1119 core records a southeasterly (seaward) movement of the STF during early glacial periods, accompanied by the incursion of subtropical water (STW) above the site, and northwesterly (landward) movement during late glacial and interglacial times, resulting in a dominant influence then of subantarctic surface water (SAW). The history of passage of these different water masses at the site is clearly delineated by their characteristic delta13C values. The intervals of thin, graded sands-muds which occur within MIS 2-3, 6, 7.4 and 10 indicate the onset at times of peak cold of intermittent bottom currents caused by strengthened and expanded frontal flows along the STF, which at such times lay near Site 1119 in close proximity to seaward-encroaching subantarctic waters within the Bounty gyre. In common with other nearby Southern Hemisphere records, the cold period which represents the last glacial maximum lasted between ~23-18 ka at Site 1119, during which time the STF and Subantarctic Front (SAF) probably merged into a single intense frontal zone around the head of the adjacent Bounty Trough.

Sea ice diatoms from Late Quaternary sediments in the Scotia Sea

Sea-ice growth and decay in Antarctica is one of the biggest seasonal changes on Earth, expanding ice cover from 4x10**6 km**2 to a maximum of 19x10**6 km**2 during the austral winter. Analyses of six marine sediment cores from the Scotia Sea, SW Atlantic, yield records of sea-ice migration across the basin since the Lateglacial. The cores span nearly ten degrees of latitude from the modern seasonal sea-ice zone to the modern Polar Front. Surface sediments in the cores comprise predominantly diatomaceous oozes and muddy diatom oozes that reflect Holocene conditions. The cores exhibit similar down-core stratigraphies with decreasing diatom concentrations and increasing magnetic susceptibility from modern through to the Last Glacial Maximum (LGM). Sediments in all cores contain sea-ice diatoms that preserve a signal of changing sea-ice cover and permit reconstruction of past sea-ice dynamics. The sea-ice records presented here are the first to document the position of both the summer and winter sea-ice cover at the Last Glacial Maximum (LGM) in the Scotia Sea. Comparison of the LGM and Holocene sea-ice conditions shows that the average winter sea-ice extent was at least 5° further north at the LGM. Average summer sea-ice extent was south of the most southerly core site at the LGM, and suggests that sea-ice expanded from approximately 61°S to 52°S each season. Our data also suggest that the average summer sea-ice position at the LGM was not the maximum extent of summer sea-ice during the last glacial. Instead, the sediments contain evidence of a pre-LGM maximum extent of summer sea-ice between ab. 30 ka and 22 ka that extended to ab. 59°S, close to the modern average winter sea-ice limit. Based on our reconstruction we propose that the timing of the maximum extent of summer sea-ice and subsequent retreat by 22 ka, could be insolation controlled and that the strong links between sea-ice and bottom water formation provide a potential mechanism by which Southern Hemisphere regional sea-ice dynamics at the LGM could have a global impact and promote deglaciation.