Elevation lines of modelled ice sheet extension during last glacial maximum

A numerical ice-sheet model was used to reconstruct the Late Weichselian glaciation of the Eurasian High Arctic, between Franz Josef Land and Severnaya Zemlya. An ice sheet was developed over the entire Eurasian High Arctic so that ice flow from the central Barents and Kara seas toward the northern Russian Arctic could be accounted for. An inverse approach to modeling was utilized, where ice-sheet results were forced to be compatible with geological information indicating ice-free conditions over the Taymyr Peninsula during the Late Weichselian. The model indicates complete glaciation of the Barents and Kara seas and predicts a "maximum-sized" ice sheet for the Late Weichselian Russian High Arctic. In this scenario, full-glacial conditions are characterized by a 1500-m-thick ice mass over the Barents Sea, from which ice flowed to the north and west within several bathymetric troughs as large ice streams. In contrast to this reconstruction, a "minimum" model of glaciation involves restricted glaciation in the Kara Sea, where the ice thickness is only 300 m in the south and which is free of ice in the north across Severnaya Zemlya. Our maximum reconstruction is compatible with geological information that indicates complete glaciation of the Barents Sea. However, geological data from Severnaya Zemlya suggest our minimum model is more relevant further east. This, in turn, implies a strong paleoclimatic gradient to colder and drier conditions eastward across the Eurasian Arctic during the Late Weichselian.

(Table 3) Sedimentation rates of stage 1 and last glacial maximum from different Holes during cruises of RRS Discovery in 1984-1985 and 1987-1988, and RRS James Clark Ross in 1992-1993

Within the Scotia Sea, the axis of the Antarctic Circumpolar Current (ACC) is geographically confined, and sediments therefore contain a record of palaeo-flow speed uncomplicated by ACC axis migration. We outline Holocene and Last Glacial Maximum (LGM) current-controlled sedimentation using data from 3.5-kHz profiles, cores and current meter moorings. Geophysical surveys show areas of erosion and deposition controlled by Neogene basement topography. Deposition occurs in mounded sediment drifts or flatter areas, where 500-1000 m of sediment overlies acoustic basement. 3.5-kHz profiles show parallel, continuous sub-bottom reflectors with highest sedimentation rates in the centre of the drifts, and reflectors converging towards marginal zones of non-deposition. Locally, on the flanks of continental blocks (e.g. South Georgia), downslope processes are dominant. The absence of mudwaves on the sediment drifts may result from the unsteadiness of ACC flow. A core transect from the ACC axis south to the boundary with the Weddell Gyre shows a southward decrease in biogenic content, controlled by the Polar Front and the spring sea-ice edge. Both these features lay farther north at LGM. The cores have been dated by relative abundance of the radiolarian Cycladophora davisiana, and by changes in the biogenic Ba content, a palaeoproductivity indicator. Sedimentation rates range from 3 to 17 cm/ka. The grain size of Holocene sediments shows a coarsening trend from south to north, consistent with strongest bottom-current flow near the ACC axis, though interpretation is complicated by the presence of biogenic grains. Year-long current meter records indicate mean speeds from 7 cm/s in the south to 12 cm/s in the north, with benthic storm frequency increasing northwards. LGM sediments are predominantly terrigenous and show a clearer northward-coarsening trend, with well-sorted silts in the northern Scotia Sea. Assuming a constant terrigenous source, this implies stronger ACC flow at the LGM, contrasting with weaker Weddell Gyre flow deduced from earlier work.

Stable oxygen isotope and Mg/Ca ratios of planktonic foraminifera from the Holocene and the Last Glacial Maximum

We measured oxygen isotopes and Mg/Ca ratios in the surface-dwelling planktonic foraminifer Globigerinoides ruber (white s.s.) and the thermocline dweller Pulleniatina obliquiloculata to investigate upper ocean spatial variability in the Indo-Pacific Warm Pool (IPWP). We focused on three critical time intervals: the Last Glacial Maximum (LGM; 18-21.5 ka), the early Holocene (8-9 ka), and the late Holocene (0-2 ka). Our records from 24 stations in the South China Sea, Timor Sea, Indonesian seas, and western Pacific indicate overall dry and cool conditions in the IPWP during the LGM with a low thermal gradient between surface and thermocline waters. During the early Holocene, sea surface temperatures increased by ~3°C over the entire region, indicating intensification of the IPWP. However, in the eastern Indian Ocean (Timor Sea), the thermocline gradually shoaled from the LGM to early Holocene, reflecting intensification of the subsurface Indonesian Throughflow (ITF). Increased surface salinity in the South China Sea during the Holocene appears related to northward displacement of the monsoonal rain belt over the Asian continent together with enhanced influx of saltier Pacific surface water through the Luzon Strait and freshwater export through the Java Sea. Opening of the freshwater portal through the Java Sea in the early Holocene led to a change in the vertical structure of the ITF from surface- to thermocline-dominated flow and to substantial freshening of Timor Sea thermocline waters.

Constraints on the magnitude and patterns of ocean cooling at the Last Glacial Maximum - Supplementary material

Observation-based reconstructions of sea surface temperature from relatively stable periods in the past, such as the Last Glacial Maximum, represent an important means of constraining climate sensitivity and evaluating model simulations. The first quantitative global reconstruction of sea surface temperatures during the Last Glacial Maximum was developed by the Climate Long-Range Investigation, Mapping and Prediction (CLIMAP) project in the 1970s and 1980s. Since that time, several shortcomings of that earlier effort have become apparent. Here we present an updated synthesis of sea surface temperatures during the Last Glacial Maximum, rigorously defined as the period between 23 and 19 thousand years before present, from the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface (MARGO) project. We integrate microfossil and geochemical reconstructions of surface temperatures and include assessments of the reliability of individual records. Our reconstruction reveals the presence of large longitudinal gradients in sea surface temperature in all of the ocean basins, in contrast to the simulations of the Last Glacial Maximum climate available at present.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Last Glacial Maximum carbon isotopic composition

We present Holocene and last glacial maximum (LGM) oxygen and carbon isotope measurements on Planulina wuellerstorfi in six southeast Pacific cores. Sedimentation rates are low in this part of the ocean, and measurements were made on individual foraminiferal shells in order to identify the Holocene and glacial individuals on the basis of their extreme d18O. The new d13C data were combined with previous P. wuellerstorfi data for interpretation of global thermohaline circulation. Data from the Southern Ocean were examined closely for regional coherency and a few anomalous d13C values suspected of having productivity overprint were removed. The resulting global d13C distributions and gradients indicate that the deep water circulation was similar during the Holocene and LGM. This interpretation brings d13C data to a better agreement with Cd/Ca data and marks a sharp contrast with a widely held view based on d13C measurements that the glacial Southern Ocean was the terminus of the thermohaline circulation. The proposed presence of glacial North Atlantic Deep Water does not necessarily contradict the postulated presence of Glacial North Atlantic Intermediate Water.

Global compilation of last glacial maximum oxygen isotopic ratios for planktonic and benthic foraminifera

We use the fully coupled atmosphere-ocean three-dimensional model of intermediate complexity iLOVECLIM to simulate the climate and oxygen stable isotopic signal during the Last Glacial Maximum (LGM, 21 000 yr). By using a model that is able to explicitly simulate the sensor (d18O), results can be directly compared with data from climatic archives in the different realms. Our results indicate that iLOVECLIM reproduces well the main feature of the LGM climate in the atmospheric and oceanic components. The annual mean d18O in precipitation shows more depleted values in the northern and southern high latitudes during the LGM. The model reproduces very well the spatial gradient observed in ice core records over the Greenland ice-sheet. We observe a general pattern toward more enriched values for continental calcite d18O in the model at the LGM, in agreement with speleothem data. This can be explained by both a general atmospheric cooling in the tropical and subtropical regions and a reduction in precipitation as confirmed by reconstruction derived from pollens and plant macrofossils. Data-model comparison for sea surface temperature indicates that iLOVECLIM is capable to satisfyingly simulate the change in oceanic surface conditions between the LGM and present. Our data-model comparison for calcite d18O allows investigating the large discrepancies with respect to glacial temperatures recorded by different microfossil proxies in the North Atlantic region. The results argue for a trong mean annual cooling between the LGM and present (>6°C), supporting the foraminifera transfer function reconstruction but in disagreement with alkenones and dinocyst reconstructions. The data-model comparison also reveals that large positive calcite d18O anomaly in the Southern Ocean may be explained by an important cooling, although the driver of this pattern is unclear. We deduce a large positive d18Osw anomaly for the north Indian Ocean that contrasts with a large negative d18Osw anomaly in the China Sea between the LGM and present. This pattern may be linked to changes in the hydrological cycle over these regions. Our simulation of the deep ocean suggests that changes in d18Osw between the LGM and present are not spatially homogenous. This is supported by reconstructions derived from pore fluids in deep-sea sediments. The model underestimates the deep ocean cooling thus biasing the comparison with benthic calcite d18O data. Nonetheless, our data-model comparison support a heterogeneous cooling of few degrees (2-4°C) in the LGM Ocean.

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.

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.