(Table S3) Benthic foraminiferal Mg/Ca, Li/Ca and stable isotope ratios of ODP Site 122-761

The sensitivities of benthic foraminiferal Mg/Ca and Li/Ca to bottom water temperature and carbonate saturation state have recently been assessed. Here we present a new approach that uses paired Mg/Ca and Li/Ca records to calculate simultaneous changes in temperature and saturation state. Using previously published records, we first use this approach to document a cooling of deep ocean waters associated with the establishment of the Antarctic ice sheet at the Eocene-Oligocene climate transition. We then apply this approach to new records of the Middle Miocene Climate Transition from ODP Site 761 to estimate variations in bottom water temperature and the oxygen isotopic composition of seawater. We estimate that the oxygen isotopic composition of seawater varied by ~1 per mil between the deglacial extreme of the Miocene Climatic Optimum and the glacial maximum following the Middle Miocene Climate Transition, indicating large amplitude variations in ice volume. However, the longer-term change between 15.3 and 12.5 Ma is marked by a ~1°C cooling of deep waters, and an increase in the oxygen isotopic composition of seawater of ~0.6 per mil. We find that bottom water saturation state increased in the lead up to the Middle Miocene Climate Transition and decreased shortly after. This supports decreasing pCO2 as a driver for global cooling and ice sheet expansion, in agreement with existing boron isotope and leaf stomatal index CO2 records but in contrast to the published alkenone CO2 records.

Ba/Ca ratios in Neogloboquadrina pachyderma of sedimet core 94B17

Down-core samples of the planktonic foraminifer Neogloboquadrina pachyderma sinistral from the Mendeleyev Ridge in the western Arctic Ocean have been analyzed for Ba/Ca and d18O. The apparent distribution coefficient for N. pachyderma sin. is estimated at DBa = 0.22 ± 0.02. A meltwater event is identified at around 11.8 14C kyr BP and is coincident with elevated Ba/Ca ratios. The barium enrichment is believed to be the result of enhanced weathering and erosion following glaciation. Additionally, barium may have desorbed from shelf sediments as sea level rose. Changes in Ba/Ca correlate with the retreat of the Laurentide Ice Sheet and the evolution of the Mackenzie River drainage basin. Therefore maximum Ba/Ca in Arctic surface waters at 11.8 ka may be indicative of an increase in the export of freshwater from the Arctic to the North Atlantic, potentially contributing to the onset of the Younger Dryas. This work suggests that Ba/Ca in planktonic foraminifera may be a useful indicator of the timing and processes associated with deglaciation.

Separating ITCZ- and ENSO-related rainfall changes in the Galápagos over the last 3 kyr using D/H ratios of multiple lipid biomarkers

We present a 3000-yr rainfall reconstruction from the Galápagos Islands that is based on paired biomarker records from the sediment of El Junco Lake. Located in the eastern equatorial Pacific, the climate of the Galápagos Islands is governed by movements of the Intertropical Convergence Zone (ITCZ) and the El Niño-Southern Oscillation (ENSO). We use a novel method for reconstructing past ENSO- and ITCZ-related rainfall changes through analysis of molecular and isotopic biomarker records representing several types of plants and algae that grow under differing climatic conditions. We propose that ?D values of dinosterol, a sterol produced by dinoflagellates, record changes in mean rainfall in El Junco Lake, while dD values of C34 botryococcene, a hydrocarbon unique to the green alga Botryococcus braunii, record changes in rainfall associated with moderate-to-strong El Niño events. We use these proxies to infer changes in mean rainfall and El Niño-related rainfall over the past 3000 yr. During periods in which the inferred change in El Niño-related rainfall opposed the change in mean rainfall, we infer changes in the amount of ITCZ-related rainfall. Simulations with an idealized isotope hydrology model of El Junco Lake help illustrate the interpretation of these proxy reconstructions. Opposing changes in El Niño- and ITCZ-related rainfall appear to account for several of the largest inferred hydrologic changes in El Junco Lake. We propose that these reconstructions can be used to infer changes in frequency and/or intensity of El Niño events and changes in the position of the ITCZ in the eastern equatorial Pacific over the past 3000 yr. Comparison with El Junco Lake sediment grain size records indicates general agreement of inferred rainfall changes over the late Holocene.

Age model and Plio-Pleistocene benthic stable oxygen isotope ratios of ODP Site 184-1143 in the Sourh China Sea

Based on benthic foraminiferal delta18O from ODP Site 1143, a 5-Myr astronomical timescale for the West Pacific Plio-Pleistocene was established using an automatic orbital tuning method. The tuned Brunhes/Matuyama paleomagnetic polarity reversal age agrees well with the previously published age of 0.78 Ma. The tuned ages for several planktonic foraminifer bio-events also agree well with published dates, and new ages for some other bio-events in the South China Sea were also estimated. The benthic delta18O from Site 1143 is highly coherent with the Earth's orbit (ETP) both at the obliquity and precession bands for the last 5 Myr, and at the eccentricity band for the last 2 Myr. In general, the 41-kyr cycle was dominant through the Plio-Pleistocene although the 23-kyr cycle was also very strong. The 100-kyr cycle became dominant only during the last 1 Myr. A comparison of the benthic delta18O between the Atlantic (ODP 659) and the East and West Pacific (846 and 1143) reveals that the Atlantic-Pacific benthic oxygen isotope difference ratio (Delta delta18OAtl-Pac) displays an increasing trend in three time intervals: 3.6-2.7 Ma, 2.7-2.1 Ma and 1.5-0.25 Ma. Each of the intervals begins with a rapid negative shift in Delta delta18OAtl-Pac, followed by a long period with an increasing trend, corresponding to the growth of the Northern Hemisphere ice sheet. This means that all three intervals of ice sheet growth in the Northern Hemisphere were accompanied at the beginning by a rapid relative warming of deep water in the Atlantic as compared to that of the Pacific, followed by its gradual relative cooling. This general trend, superimposed on the frequent fluctuations with glacial cycles, should yield insights into the processes leading to the boreal glaciation. Cross-spectral analyses of the Delta delta18OAtl-Pac with the Earth's orbit suggests that after the initiation of Northern Hemisphere glaciation at about 2.5 Ma, obliquity rather than precession had become the dominant force controlling the vertical structure or thermohaline circulation in the paleo-ocean.

Planktic foraminiferal abundance and stable oxygen isotope ratios of sediment cores from the Southern California Borderlands

A variety of evidence suggests that average sea surface temperatures (SSTs) during the last glacial maximum in the California Borderlands region were significantly colder than during the Holocene. Planktonic foraminiferal delta18O evidence and average SST estimates derived by the modern analog technique indicate that temperatures were 6°-10°C cooler during the last glacial relative to the present. The glacial plankton assemblage is dominated by the planktonic foraminifer Neogloboquadrina pachyderma (sinistral coiling) and the coccolith Coccolithus pelagicus, both of which are currently restricted to subpolar regions of the North Pacific. The glacial-interglacial average SST change determined in this study is considerably larger than the 2°C change estimated by Climate: Long-Range Investigation, Mapping, and Prediction (CLIMAP) [1981]. We propose that a strengthened California Current flow was associated with the advance of subpolar surface waters into the Borderlands region during the last glacial.

Mg/Ca ratios in planktonic foraminifera

The primary Mg/Ca ratio of foraminiferal shells is a potentially valuable paleoproxy for sea surface temperature (SST) reconstructions. However, the reliable extraction of this ratio from sedimentary calcite assumes that we can overcome artifacts related to foraminiferal ecology and partial dissolution, as well as contamination by secondary calcite and clay. The standard batch method for Mg/Ca analysis involves cracking, sonicating, and rinsing the tests to remove clay, followed by chemical cleaning, and finally acid-digestion and single-point measurement. This laborious procedure often results in substantial loss of sample (typically 30-60%). We find that even the earliest steps of this procedure can fractionate Mg from Ca, thus biasing the result toward a more variable and often anomalously low Mg/Ca ratio. Moreover, the more rigorous the cleaning, the more calcite is lost, and the more likely it becomes that any residual clay that has not been removed by physical cleaning will increase the ratio. These potentially significant sources of error can be overcome with a flow-through (FT) sequential leaching method that makes time- and labor-intensive pretreatments unnecessary. When combined with time-resolved analysis (FT-TRA) flow-through, performed with a gradually increasing and highly regulated acid strength, produces continuous records of Mg, Sr, Al, and Ca concentrations in the leachate sorted by dissolution susceptibility of the reacting material. Flow-through separates secondary calcite from less susceptible biogenic calcite and clay, and further resolves the biogenic component into primary and more resistant fractions. FT-TRA reliably separates secondary calcite (which is not representative of original life habitats) from the more resistant biogenic calcite (the desired signal) and clay (a contaminant of high Mg/Ca, which also contains Al), and further resolves the biogenic component into primary and more resistant fractions that may reflect habitat or other changes during ontogeny. We find that the most susceptible fraction of biogenic calcite in surface dwelling foraminifera gives the most accurate value for SST and therefore best represents primary calcite. Sequential dissolution curves can be used to correct the primary Mg/Ca ratio for clay, if necessary. However, the temporal separation of calcite from clay in FT-TRA is so complete that this correction is typically <=2%, even in clay-rich sediments. Unlike hands-on batch methods, that are difficult to reproduce exactly, flow-through lends itself to automation, providing precise replication of treatment for every sample. Our automated flow-through system can process 22 samples, two system blanks, and 48 mixed standards in <12 hours of unattended operation. FT-TRA thus represents a faster, cheaper, and better way to determine Mg/Ca ratios in foraminiferal calcite.

TanDEM-X ice thickness dataset for Wilkins Ice Shelf, Antarctic Peninsula, with link to GeoTIFF

The dataset shows the ice thickness over Wilkins Ice Shelf, Antarctic Peninsula derived from TanDEM-X Interferometry. The data has been acquired between June and August 2012. The TanDEM-X heights have been linked to CryoSAT-2 heights (V. Helm) from the respective time stamp. Elevations have been transformed from WGS84 ellipsoidal heights to the EGM2008 geoid. The ice shelf thickness was estimated assuming hydrostatic equilibrium and a mean ice density of 915 kg/m³.

Environmental variables and multivariate seabed classification maps via k means clustering of Potter Cove, Antarctica, link to input and result files in GeoTIFF format

This study subdivides the Potter Cove, King George Island, Antarctica, into seafloor regions using multivariate statistical methods. These regions are categories used for comparing, contrasting and quantifying biogeochemical processes and biodiversity between ocean regions geographically but also regions under development within the scope of global change. The division obtained is characterized by the dominating components and interpreted in terms of ruling environmental conditions. The analysis includes in total 42 different environmental variables, interpolated based on samples taken during Australian summer seasons 2010/2011 and 2011/2012. The statistical errors of several interpolation methods (e.g. IDW, Indicator, Ordinary and Co-Kriging) with changing settings have been compared and the most reasonable method has been applied. The multivariate mathematical procedures used are regionalized classification via k means cluster analysis, canonical-correlation analysis and multidimensional scaling. Canonical-correlation analysis identifies the influencing factors in the different parts of the cove. Several methods for the identification of the optimum number of clusters have been tested and 4, 7, 10 as well as 12 were identified as reasonable numbers for clustering the Potter Cove. Especially the results of 10 and 12 clusters identify marine-influenced regions which can be clearly separated from those determined by the geological catchment area and the ones dominated by river discharge.

Lithogeochemical and isotope (Rb-Sr, Sm-Nd, Pb-Pb whole rock and Lu-Hf zircon) composition of samples from the Shackleton Range, East Antarctica

Three distinct, spatially separated crustal terranes have been recognised in the Shackleton Range, East Antarctica: the Southern, Eastern and Northern Terranes. Mafic gneisses from the Southern Terrane provide geochemical evidence for a within-plate, probably back-arc origin of their protoliths. A plume-distal ridge origin in an incipient ocean basin is the favoured interpretation for the emplacement site of these rocks at c. 1850 Ma, which, together with a few ocean island basalts, were subsequently incorporated into an accretionary continental arc/supra-subduction zone tectonic setting. Magmatic underplating resulted in partial melting of the lower crust, which caused high-temperature granulite-facies metamorphism in the Southern Terrane at c. 1710-1680 Ma. Mafic and felsic gneisses there are characterised by isotopically depleted, positive Nd and Hf initials and model ages between 2100 and 2000 Ma. They may be explained as juvenile additions to the crust towards the end of the Palaeoproterozoic. These juvenile rocks occur in a narrow, c. 150 km long E-W trending belt, inferred to trace a suture that is associated with a large Palaeoproterozoic accretionary orogenic system. The Southern Terrane contains many features that are similar to the Australo-Antarctic Mawson Continent and may be its furthermost extension into East Antarctica. The Eastern Terrane is characterised by metagranitoids that formed in a continental volcanic arc setting during a late Mesoproterozoic orogeny at c. 1060 Ma. Subsequently, the rocks experienced high-temperature metamorphism during Pan-African collisional tectonics at 600 Ma. Isotopically depleted zircon grains yielded Hf model ages of 1600-1400 Ma, which are identical to Nd model ages obtained from juvenile metagranitoids. Most likely, these rocks trace the suture related to the amalgamation of the Indo-Antarctic and West Gondwana continental blocks at ~600 Ma. The Eastern Terrane is interpreted as the southernmost extension of the Pan-African Mozambique/Maud Belt in East Antarctica and, based on Hf isotope data, may also represent a link to the Ellsworth-Whitmore Mountains block in West Antarctica and the Namaqua-Natal Province of southern Africa. Geochemical evidence indicates that the majority of the protoliths of the mafic gneisses in the Northern Terrane formed as oceanic island basalts in a within-plate setting. Subsequently the rocks were incorporated into a subduction zone environment and, finally, accreted to a continental margin during Pan-African collisional tectonics. Felsic gneisses there provide evidence for a within-plate and volcanic arc/collisional origin. Emplacement of granitoids occurred at c. 530 Ma and high-temperature, high-pressure metamorphism took place at 510-500 Ma. Enriched Hf and Nd initials and Palaeoproterozoic model ages for most samples indicate that no juvenile material was added to the crust of the Northern Terrane during the Pan-African Orogeny but recycling of older crust or mixing of crustal components of different age must have occurred. Isotopically depleted mafic gneisses, which are spatially associated with eclogite-facies pyroxenites, yielded late Mesoproterozoic Nd model ages. These rocks occur in a narrow, at least 100 km long, E-W trending belt that separates alkaline ocean island metabasalts and within-plate metagranitoids from volcanic arc metabasalts and volcanic arc/syn-collisional metagranitoids in the Northern Terrane. This belt is interpreted to trace the late Neoproterozoic/early Cambrian Pan-African collisional suture between the Australo-Antarctic and the combined Indo-Antarctic/West Gondwana continental blocks that formed during the final amalgamation of Gondwana.