Table 2. Lead and Ba concentrations and Pb isotopic compositions measured in Law Dome ice cores

Lead isotopic compositions and Pb and Ba concentrations have been measured in ice cores from Law Dome, East Antarctica, covering the past 6500 years. 'Natural' background concentrations of Pb (ab. 0.4 pg/g) and Ba (ab. 1.3 pg/g) are observed until 1884 AD, after which increased Pb concentrations and lowered 206Pb/207Pb ratios indicate the influence of anthropogenic Pb. The isotopic composition of 'natural' Pb varies within the range 206Pb/207Pb=1.20-1.25 and 208Pb/207Pb=2.46-2.50, with an average rock and soil dust Pb contribution of 8-12%. A major pollution event is observed at Law Dome between 1884 and 1908 AD, elevating the Pb concentration four-fold and changing 206Pb/207Pb ratios in the ice to ab. 1.12. Based on Pb isotopic systematics and Pb emission statistics, this is attributed to Pb mined at Broken Hill and smelted at Broken Hill and Port Pirie, Australia. Anthropogenic Pb inputs are at their greatest from 1900 to 1910 and from ab. 1960 to ab. 1980. During the 20th century, Ba concentrations are consistently higher than 'natural' levels and are attributed to increased dust production, suggesting the influence of climate change and/or changes in land coverage with vegetation.

Code and readme of a MATLAB-based graphical user interface program for computing functionals of the geopotential

We present a novel graphical user interface program GrafLab (GRAvity Field LABoratory) for spherical harmonic synthesis (SHS) created in MATLAB®. This program allows to comfortably compute 38 various functionals of the geopotential up to ultra-high degrees and orders of spherical harmonic expansion. For the most difficult part of the SHS, namely the evaluation of the fully normalized associated Legendre functions (fnALFs), we used three different approaches according to required maximum degree: (i) the standard forward column method (up to maximum degree 1800, in some cases up to degree 2190); (ii) the modified forward column method combined with Horner's scheme (up to maximum degree 2700); (iii) the extended-range arithmetic (up to an arbitrary maximum degree). For the maximum degree 2190, the SHS with fnALFs evaluated using the extended-range arithmetic approach takes only approximately 2-3 times longer than its standard arithmetic counterpart, i.e. the standard forward column method. In the GrafLab, the functionals of the geopotential can be evaluated on a regular grid or point-wise, while the input coordinates can either be read from a data file or entered manually. For the computation on a regular grid we decided to apply the lumped coefficients approach due to significant time-efficiency of this method. Furthermore, if a full variance-covariance matrix of spherical harmonic coefficients is available, it is possible to compute the commission errors of the functionals. When computing on a regular grid, the output functionals or their commission errors may be depicted on a map using automatically selected cartographic projection.

(Table S1) Globigerinoides sacculifer Boron isotope ratios of ODP Hole 165-999A

everal hypotheses have been put forward to explain the onset of intensive glaciations on Greenland, Scandinavia, and North America during the Pliocene epoch between 3.6 and 2.7 million years ago (Ma). A decrease in atmospheric CO2 may have played a role during the onset of glaciations, but other tectonic and oceanic events occurring at the same time may have played a part as well. Here we present detailed atmospheric CO2 estimates from boron isotopes in planktic foraminifer shells spanning 4.6-2.0 Ma. Maximal Pliocene atmospheric CO2 estimates gradually declined from values around 410 µatm to early Pleistocene values of 300 ?atm at 2.0 Ma. After the onset of large-scale ice sheets in the Northern Hemisphere, maximal pCO2 estimates were still at 2.5 Ma +90 µatm higher than values characteristic of the early Pleistocene interglacials. By contrast, Pliocene minimal atmospheric CO2 gradually decreased from 310 to 245 µatm at 3.2 Ma, coinciding with the start of transient glaciations on Greenland. Values characteristic of early Pleistocene glacial atmospheric CO2 of 200 ?atm were abruptly reached after 2.7 Ma during the late Pliocene transition. This trend is consistent with the suggestion that ocean stratification and iron fertilization increased after 2.7 Ma in the North Pacific and Southern Ocean and may have led to increased glacial CO2 storage in the oceanic abyss after 2.7 Ma onward.