TALDICE-1 age scale of the Talos Dome deep ice core, East Antarctica

A new deep ice core drilling program, TALDICE, has been successfully handled by a European team at Talos Dome, in the Ross Sea sector of East Antarctica, down to 1620 m depth. Using stratigraphic markers and a new inverse method, we produce the first official chronology of the ice core, called TALDICE-1. We show that it notably improves an a priori chronology resulting from a one-dimensional ice flow model. It is in agreement with a posteriori controls of the resulting accumulation rate and thinning function along the core. An absolute uncertainty of only 300 yr is obtained over the course of the last deglaciation. This uncertainty remains lower than 600 yr over Marine Isotope Stage 3, back to 50 kyr BP. The phasing of the TALDICE ice core climate record with respect to the central East Antarctic plateau and Greenland records can thus be determined with a precision allowing for a discussion of the mechanisms at work at sub-millennial time scales.

Martian sub-surface ionising radiation: biosignature and geology

The surface of Mars, unshielded by thick atmosphere or global magnetic field, is exposed to high levels of cosmic radiation. This ionising radiation field is deleterious to the survival of dormant cells or spores and the persistence of molecular biomarkers in the subsurface, and so its characterisation is of prime astrobiological interest. Here, we present modelling results of the absorbed radiation dose as a function of depth through the Martian subsurface, suitable for calculation of biomarker persistence. A second major implementation of this dose accumulation rate data is in application of the optically stimulated luminescence technique for dating Martian sediments. We present calculations of the dose-depth profile in the Martian subsurface for various scenarios: variations of surface composition (dry regolith, ice, layered permafrost), solar minimum and maximum conditions, locations of different elevation (Olympus Mons, Hellas basin, datum altitude), and increasing atmospheric thickness over geological history. We also model the changing composition of the subsurface radiation field with depth compared between Martian locations with different shielding material, determine the relative dose contributions from primaries of different energies, and discuss particle deflection by the crustal magnetic fields.

How stable are 20th century calibration models? A high-resolution summer temperature reconstruction for the eastern Swiss Alps back to A.D. 1580 derived from proglacial varved sediments

We found a significant positive correlation between local summer air temperature (May-September) and the annual sediment mass accumulation rate (MAR) in Lake Silvaplana (46°N, 9°E, 1800 m a.s.l.) during the twentieth century (r = 0.69, p < 0.001 for decadal smoothed series). Sediment trap data (2001-2005) confirm this relation with exceptionally high particle yields during the hottest summer of the last 140 years in 2003. On this base we developed a decadal-scale summer temperature reconstruction back to AD 1580. Surprisingly, the comparison of our reconstruction with two other independent regional summer temperature reconstructions (based on tree-rings and documentary data) revealed a significant negative correlation for the pre-1900 data (ie, late ‘Little Ice Age’). This demonstrates that the correlation between MAR and summer temperature is not stable in time and the actualistic principle does not apply in this case. We suggest that different climatic regimes (modern/‘Little Ice Age’) lead to changing state conditions in the catchment and thus to considerably different sediment transport mechanisms. Therefore, we calibrated our MAR data with gridded early instrumental temperature series from AD 1760-1880 (r = -0.48, p < 0.01 for decadal smoothed series) to properly reconstruct the late LIA climatic conditions. We found exceptionally low temperatures between AD 1580 and 1610 (0.75°C below twentieth-century mean) and during the late Maunder Minimum from AD 1680 to 1710 (0.5°C below twentieth-century mean). In general, summer temperatures did not experience major negative departures from the twentieth-century mean during the late ‘Little Ice Age’. This compares well with the two existing independent regional reconstructions suggesting that the LIA in the Alps was mainly a phenomenon of the cold season.

Geochemical and isotopic characterization of the Bodélé Depression dust source and implications for transatlantic dust transport to the Amazon Basin

The Bodélé Depression (Chad) in the central Sahara/Sahel region of Northern Africa is the most important source of mineral dust to the atmosphere globally. The Bodélé Depression is purportedly the largest source of Saharan dust reaching the Amazon Basin by transatlantic transport. Here, we have undertaken a comprehensive study of surface sediments from the Bodélé Depression and dust deposits (Chad, Niger) in order to characterize geochemically and isotopically (Sr, Nd and Pb isotopes) this dust source, and evaluate its importance in present and past African dust records. We similarly analyzed sedimentary deposits from the Amazonian lowlands in order to assess postulated accumulation of African mineral dust in the Amazon Basin, as well as its possible impact in fertilizing the Amazon rainforest. Our results identify distinct sources of different ages and provenance in the Bodélé Depression versus the Amazon Basin, effectively ruling out an origin for the Amazonian deposits, such as the Belterra Clay Layer, by long-term deposition of Bodélé Depression material. Similarly, no evidence for contributions from other potential source areas is provided by existing isotope data (Sr, Nd) on Saharan dusts. Instead, the composition of these Amazonian deposits is entirely consistent with derivation from in-situ weathering and erosion of the Precambrian Amazonian craton, with little, if any, Andean contribution. In the Amazon Basin, the mass accumulation rate of eolian dust is only around one-third of the vertical erosion rate in shield areas, suggesting that Saharan dust is “consumed” by tropical weathering, contributing nutrients and stimulating plant growth, but never accumulates as such in the Amazon Basin. The chemical and isotope compositions found in the Bodélé Depression are varied at the local scale, and have contrasting signatures in the “silica-rich” dry lake-bed sediments and in the “calcium-rich” mixed diatomites and surrounding sand material. This unexpected finding implies that the Bodélé Depression material is not “pre-mixed” at the source to provide a homogeneous source of dust. Rather, different isotope signatures can be emitted depending on subtle vagaries of dust-producing events. Our characterization of the Bodélé Depression components indicate that the Bodélé “calcium-rich” component, identified here, is most likely released via eolian processes of sand grain saltation and abrasion and may be significant in the overall global budget of dusts carried out by the Harmattan low-level jet during the winter.

Glacial–interglacial dynamics of Antarctic firn columns: comparison between simulations and ice core air-δ15N measurements

Correct estimation of the firn lock-in depth is essential for correctly linking gas and ice chronologies in ice core studies. Here, two approaches to constrain the firn depth evolution in Antarctica are presented over the last deglaciation: outputs of a firn densification model, and measurements of δ15N of N2 in air trapped in ice core, assuming that δ15N is only affected by gravitational fractionation in the firn column. Since the firn densification process is largely governed by surface temperature and accumulation rate, we have investigated four ice cores drilled in coastal (Berkner Island, BI, and James Ross Island, JRI) and semi-coastal (TALDICE and EPICA Dronning Maud Land, EDML) Antarctic regions. Combined with available ice core air-δ15N measurements from the EPICA Dome C (EDC) site, the studied regions encompass a large range of surface accumulation rates and temperature conditions. Our δ15N profiles reveal a heterogeneous response of the firn structure to glacial–interglacial climatic changes. While firn densification simulations correctly predict TALDICE δ15N variations, they systematically fail to capture the large millennial-scale δ15N variations measured at BI and the δ15N glacial levels measured at JRI and EDML – a mismatch previously reported for central East Antarctic ice cores. New constraints of the EDML gas–ice depth offset during the Laschamp event (~41 ka) and the last deglaciation do not favour the hypothesis of a large convective zone within the firn as the explanation of the glacial firn model–δ15N data mismatch for this site. While we could not conduct an in-depth study of the influence of impurities in snow for firnification from the existing datasets, our detailed comparison between the δ15N profiles and firn model simulations under different temperature and accumulation rate scenarios suggests that the role of accumulation rate may have been underestimated in the current description of firnification models.

Local Rates of Ice-Sheet Thickness Change in Greenland

The rate of ice-sheet thickness change is calculated for 10 sites in Greenland by comparing measured values of ice vertical velocity and snow-accumulation rate. Vertical velocities are derived from repeat surveys of markers using precision global positioning system techniques, and accumulation rates are determined from stratigraphic analysis of firn cores. The results apply to time-scales covered by the firn-core records, which in most cases are a few decades. A spectrum of thickness-change rates is obtained, ranging from substantial thinning to slow thickening. The sites where ice-sheet thinning is indicated are located near the ice-sheet margin or in outlet glacier catchments. Interior and high-elevation sites are predominantly in balance or thickening slowly. Uncertainties in the rates of thickness change are dominated by errors in the determination of accumulation rates. The results of this work are broadly comparable with regional estimates of mass balance obtained from the analysis of catchment input vs discharge.

Stratigraphic Continuity in 400 Mhz Short-Pulse Radar Profiles of Firn in West Antarctica

We track dated firn horizons within 400 MHz short-pulse radar profiles to find the continuous extent over which they can be used as historical benchmarks to study past accumulation rates in West Antarctica. The 30-40 cm pulse resolution compares with the accumulation rates of most areas. We tracked a particular set that varied from 30 to 90 m in depth over a distance of 600 km. The main limitations to continuity are fading at depth, pinching associated with accumulation rate differences within hills and valleys, and artificial fading caused by stacking along dips. The latter two may be overcome through multi-kilometer distances by matching the relative amplitude and spacing of several close horizons, along with their pulse forms and phases. Modeling of reflections from thin layers suggests that the - 37 to - 50 dB range of reflectivity and the pulse waveforms we observed are caused by the numerous thin ice layers observed in core stratigraphy. Constructive interference between reflections from these close, high-density layers can explain the maintenance of reflective strength throughout the depth of the firn despite the effects of compaction. The continuity suggests that these layers formed throughout West Antarctica and possibly into East Antarctica as well.

A 700-Year Record of Atmospheric Circulation Developed from the Law Dome Ice Core, East Antarctica

A 700-year, high-resolution, multivariate ice core record from Dome Summit South (DSS) (66degrees46'S, 112degrees48'E; 1370 m), Law Dome, is used to investigate sea level pressure (SLP) variability in the region of East Antarctica. Empirical orthogonal function (EOF) analysis reveals that the first EOF (LDEOF1) of the combined glaciochemical, oxygen isotope ratio, and accumulation rate record from DSS represents most of the variability in sea salt seen in the record. LDEOF1 is positively correlated (at least 95% confidence level) to instrumental June mean SLP across most of East Antarctica. Over the last 700 years, LDEOF1 levels at Law Dome were the highest during the nineteenth century, suggesting an increase in intensification of winter circulation during this period. The Law Dome DSS oxygen isotope ratio series also indicates that the nineteenth century had the coldest winters of any century in the record. In contrast, LDEOF1 levels were the lowest at Law Dome during the eighteenth century, suggesting a significant shift in the patterns and/or intensity of East Antarctic atmospheric circulation between the eighteenth and the nineteenth centuries. The LDEOF1 sea salt record is characterized by significant decadal-scale variability with a strong 25-year periodic structure.

Influence of Equatorial Diatom Processes on Si Deposition and Atmospheric CO(2) Cycles at Glacial/Interglacial Timescales

The causes of the glacial cycle remain unknown, although the primary driver is changes in atmospheric CO(2), likely controlled by the biological pump and biogeochemical cycles. The two most important regions of the ocean for exchange of CO(2) with the atmosphere are the equatorial Pacific and the Southern Ocean ( SO), the former a net source and the latter a net sink under present conditions. The equatorial Pacific has been shown to be a Si(OH)(4)-limited ecosystem, a consequence of the low source Si(OH)(4) concentrations in upwelled water that has its origin in the SO. This teleconnection for nutrients between the two regions suggests an oscillatory relationship that may influence or control glacial cycles. Opal mass accumulation rate (MAR) data and delta(15)N measurements in equatorial cores are interpreted with predictions from a one- dimensional Si(OH)(4)-limited ecosystem model (CoSINE) for the equatorial Pacific. The results suggest that equatorial Pacific surface CO(2) processes are in opposite phase to that of the global atmosphere, providing a negative feedback to the glacial cycle. This negative feedback is implemented through the effect of the SO on the equatorial Si(OH)(4) supply. An alternative hypothesis, that the whole ocean becomes Si(OH)(4) poor during cooling periods, is suggested by low opal MAR in cores from both equatorial and Antarctic regions, perhaps as a result of low river input. terminations in this scenario would result from blooms of coccolithophorids triggered by low Si(OH)(4) concentrations.

Temperature reconstruction from 10 to 120 kyr b2k from the NGRIP ice core

In order to reconstruct the temperature of the North Greenland Ice Core Project (NGRIP) site, new measurements of δ15N have been performed covering the time period from the beginning of the Holocene to Dansgaard–Oeschger (DO) event 8. Together with previously measured and mostly published δ15N data, we present for the first time a NGRIP temperature reconstruction for the whole last glacial period from 10 to 120 kyr b2k (thousand years before 2000 AD) including every DO event based on δ15N isotope measurements combined with a firn densification and heat diffusion model. The detected temperature rises at the onset of DO events range from 5 °C (DO 25) up to 16.5 °C (DO 11) with an uncertainty of ±3 °C. To bring measured and modelled data into agreement, we had to reduce the accumulation rate given by the NGRIP ss09sea06bm timescale in some periods by 30 to 35%, especially during the last glacial maximum. A comparison between reconstructed temperature and δ18Oice data confirms that the isotopic composition of the stadial was strongly influenced by seasonality. We evidence an anticorrelation between the variations of the δ18Oice sensitivity to temperature (referred to as α) and obliquity in agreement with a simple Rayleigh distillation model. Finally, we suggest that α might be influenced by the Northern Hemisphere ice sheet volume.

Alpine lacustrine varved record reveals summer temperature as main control of glacier fluctuations over the past 2250 years

Glacier fluctuations are a key indicator of changing climate. Their reconstruction beyond historical times unravels glacier variability and its forcing factors on long time scales, which can considerably improve our understanding of the climate–glacier relationship. Here, we present a 2250-year-long reconstruction of particle-mass accumulation rates recorded in the lacustrine sediments of Lake Trüebsee (Central Swiss Alps) that are directly related to glacier extent, thus reflecting a continuous record of fluctuations of the upstream-located Titlis Glacier. Mass accumulation rate values show strong centennial to multi-centennial fluctuations and reveal 12 well-pronounced periods of enhanced values corresponding to times of maximum extent of the neighboring Lower Grindelwald Glacier. This result supports previous studies of proglacial lake sediments that documented high mass accumulation rate values during glacier advances. The strong variability in the Lake Trüebsee mass accumulation rate record thus represents a highly sensitive paleoclimatic archive, which mirrors rapid and pronounced feedbacks of Titlis Glacier to climatic changes over the past 2250years. The comparison of our data with independent paleo-temperature reconstructions from tree rings suggests that variations in mean summer temperature were the primary driving factor of fluctuations of Titlis Glacier. Also, advances of Titlis Glacier occurred during the grand solar minima (Dalton, Maunder, Spörer, Wolf) of the last millennium. This relation of glacier extent with summer temperature reveals strong evidence that the mass balance of this Alpine glacier is primarily controlled by the intensity of glacier melting during summer.

Estimation and calibration of the water isotope differential diffusion length in ice core records

Palaeoclimatic information can be retrieved from the diffusion of the stable water isotope signal during firnification of snow. The diffusion length, a measure for the amount of diffusion a layer has experienced, depends on the firn temperature and the accumulation rate. We show that the estimation of the diffusion length using power spectral densities (PSDs) of the record of a single isotope species can be biased by uncertainties in spectral properties of the isotope signal prior to diffusion. By using a second water isotope and calculating the difference in diffusion lengths between the two isotopes, this problem is circumvented. We study the PSD method applied to two isotopes in detail and additionally present a new forward diffusion method for retrieving the differential diffusion length based on the Pearson correlation between the two isotope signals. The two methods are discussed and extensively tested on synthetic data which are generated in a Monte Carlo manner. We show that calibration of the PSD method with this synthetic data is necessary to be able to objectively determine the differential diffusion length. The correlation-based method proves to be a good alternative for the PSD method as it yields precision equal to or somewhat higher than the PSD method. The use of synthetic data also allows us to estimate the accuracy and precision of the two methods and to choose the best sampling strategy to obtain past temperatures with the required precision. In addition to application to synthetic data the two methods are tested on stable-isotope records from the EPICA (European Project for Ice Coring in Antarctica) ice core drilled in Dronning Maud Land, Antarctica, showing that reliable firn temperatures can be reconstructed with a typical uncertainty of 1.5 and 2 °C for the Holocene period and 2 and 2.5 °C for the last glacial period for the correlation and PSD method, respectively.

Global pulses of organic carbon burial in deep-sea sediments during glacial maxima

The burial of organic carbon in marine sediments removes carbon dioxide from the ocean–atmosphere pool, provides energy to the deep biosphere, and on geological timescales drives the oxygenation of the atmosphere. Here we quantify natural variations in the burial of organic carbon in deep-sea sediments over the last glacial cycle. Using a new data compilation of hundreds of sediment cores, we show that the accumulation rate of organic carbon in the deep sea was consistently higher (50%) during glacial maxima than during interglacials. The spatial pattern and temporal progression of the changes suggest that enhanced nutrient supply to parts of the surface ocean contributed to the glacial burial pulses, with likely additional contributions from more efficient transfer of organic matter to the deep sea and better preservation of organic matter due to reduced oxygen exposure. These results demonstrate a pronounced climate sensitivity for this global carbon cycle sink.

Weather and Growing Season Summary

Precipitation for 2011 was less than the longterm climate average. Early in the year, precipitation lagged behind normal, but then tracked close to the normal accumulation rate from mid-April through mid-August. After that time, precipitation amounts greatly lagged behind normal, and the year ended almost 7 in. behind the long-term average. (Figure 1). Overall, 2011 will be remembered for good moisture early, but ending the season with almost no rainfall.