Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Chemical and boron isotopic compositions of volcanics from ODP Sites 135-834 to 135-840

The influence of fluid flux on petrogenesis in the Tonga-Kermadec Arc was investigated using ion microprobe measurements of B/Be and boron isotope ratios (11B/10B) to document the source and relative volumes of the fluids released from the subducting oceanic plate. We analyzed young lavas from eight different islands along the Tonga-Kermadec Arc, as well as glass shards in volcanic sediments from Ocean Drilling Program (ODP) Site 840, which record the variations in the chemistry of Tonga magmatism since 7 Ma. B/Be is variable (5.8-122), in young Tonga-Kermadec Arc lavas. In contrast, glass shards from around 3 to 4 Ma old volcanic sediments at Site 840 have the highest B/Be values yet reported for arc lavas (18-607). These values are too high to be related simply to a sediment influence on petrogenesis. Together with very high d11B values (-11.6 to +37.5) for the same shards and lavas these data indicate that most of the B is derived from fluid escaped from the subducting altered Pacific oceanic crust, rather than from sediment. High d11B values also reflect large degrees of isotopic fractionation in this cold fast subduction zone. Lower d11B values noted in the Kermadec Arc (17 to -4.4) are related to the influence of sediment eroded from New Zealand and slower convergence. High fluid flux (B/Be) is synchronous in Tonga and the Marianas at 3 to 4 Ma and may be related to acceleration of the Pacific Plate just prior to this time. The timing of maximum B/Be at 3 to 4 Ma correlates with maximum light rare earth (LREE) and high field strength element depletion. This suggests maximum degrees of partial melting at this time. Although thinning of the arc lithosphere during rifting to form the Lau Basin is expected to influence the arc geochemistry, variable aqueous fluid flux from the subducting plate alone appears capable of explaining boron and other trace element systematics in the Tonga-Kermadec Arc with no indication of slab melting.

Water chemistry and iron speciation in samples taken during SWEDARP 97/98 campaign with S.A. Agulhas to the Southern Ocean

The distribution and speciation of iron was determined along a transect in the eastern Atlantic sector (6°E) of the Southern Ocean during a collaborative Scandinavian/South African Antarctic cruise conducted in late austral summer (December 1997/January 1998). Elevated concentrations of dissolved iron (>0.4 nM) were found at 60°S in the vicinity of the Spring Ice Edge (SIE) in tandem with a phytoplankton bloom, chiefly dominated by Phaeocystis sp. This bloom had developed rapidly after the loss of the seasonal sea ice cover. The iron that fuelled this bloom was mostly likely derived from sea ice melt. In the Winter Ice Edge (WIE), around 55°S, dissolved iron concentrations were low (<0.2 nM) and corresponded to lower biological productivity, biomass. In the Antarctic Polar Front, at approximately 50°S, a vertical profile of dissolved iron showed low concentrations (<0.2 nM); however, a surface survey showed higher concentrations (1-3 nM), and considerable patchiness in this dynamic frontal region. The chemical speciation of iron was dominated by organic complexation throughout the study region. Organic iron-complexing ligands ([L]) ranged from 0.9 to 3.0 nM Fe equivalents, with complex stability log K'(FeL) = 21.4-23.5. Estimated concentrations of inorganic iron (Fe') ranged from 0.03 to 0.79 pM, with the highest values found in the Phaeocystis bloom in the SIE. A vertical profile of iron-complexing ligands in the WIE showed a maximum consistent with a biological source for ligand production and near surface minimum possibly consistent with loss via photodecomposition. This work further confirms the role iron that has in the Southern Ocean in limiting primary productivity.

Registry of all campaigns (from port to port) from the Tara Oceans Expedition (2009-2013)

The Tara Oceans Expedition (2009-2013) sampled the world oceans on board a 36 m long schooner, collecting environmental data and organisms from viruses to planktonic metazoans for later analyses using modern sequencing and state-of-the-art imaging technologies. Tara Oceans Data are particularly suited to study the genetic, morphological and functional diversity of plankton. The present data set is a registry of all campaigns (from port to port) conducted during the Tara Oceans Expedition (2009-2013). The registry provides details about the scientific interest of each campaign, including direct links (URLs) to the corresponding (1) campaign report, (2) samples collected during the campaign, (3) environmental data published at PANGAEA, and (4) nucleotides data published at the European Nucleotides Archive (EBI-ENA).

Condensation versus air density change as a major cause of airflow: experimental evidence, link to data files and videos

The dominant model of atmospheric circulation posits that hot air rises, creating horizontal winds. A second major driver has recently been proposed by Makarieva and Gorshkov in their biotic pump theory (BPT), which suggests that evapotranspiration from natural closed-canopy forests causes intense condensation, and hence winds from ocean to land. Critics of the BPT argue that air movement to fill the partial vacuum caused by condensation is always isotropic, and therefore causes no net air movement (Bunyard, 2015, hdl:11232/397). This paper explores the physics of water condensation under mild atmospheric conditions, within a purpose-designed square-section 4.8 m-tall closed-system structure. Two enclosed vertical columns are connected at top and bottom by two horizontal tunnels, around which 19.5 m**3 of atmospheric air can circulate freely, allowing rotary airflows in either direction. This air can be cooled and/or warmed by refrigeration pipes and a heating mat, and changes in airflow, temperature, humidity and barometric pressure measured in real time. The study investigates whether the "hot-air-rises" or an implosive condensation model can better explain the results of more than 100 experiments. The data show a highly significant correlation (R2 >0.96, p value <0.001) between observed airflows and partial pressure changes from condensation. While the kinetic energy of the refrigerated air falls short of that required in bringing about observed airflows by a factor of at least 30, less than a tenth of the potential kinetic energy from condensation is shown to be sufficient. The assumption that condensation of water vapour is always isotropic is therefore incorrect. Condensation can be anisotropic, and in the laboratory does cause sustained airflow.

(Table 1) Sulphur concentration and isotopic composition of ODP Leg 126 igneous rocks

Sulfur isotope ratios have been determined in 19 selected igneous rocks from Leg 126. The d34S of the analyzed rocks ranges from -0.1 â to +19.60 â. The overall variation in sulfur isotope composition of the rocks is caused by varying degrees of seawater alteration. Most of the samples are altered by seawater and only five of them are considered to have maintained their magmatic sulfur isotope composition. These samples are all from the backarc sites and have d34S values varying from +0.2 â to +1.6 â, of which the high d34S values suggest that the earliest magmas in the rift are more arc-like in their sulfur isotope composition than the later magmas. The d34S values from the forearc sites are similar to or heavier than the sulfur isotope composition of the present arc.

(Table 1) Geochemistry of realatively fresh rocks of ODP Sites 127-794 and 127-797

The flows and sills drilled at Sites 794 and 797 in the Yamato Basin of the Japan Sea are subalkalic, olivine, and/or plagioclase phyric basalts. Compositionally, the rocks can be divided into a depleted, low-K type and an enriched, relatively high-K type. In addition, two contrasting evolution trends are reflected in the rock compositions, which allow four different magmatic suites to be identified. It is suggested that the depleted or enriched nature of these suites represent primary characteristics, while the different evolution trends are related to fractionation processes in crustal magma chambers. A tholeiitic evolution trend, with increasing FeO and TiO2 and decreasing Al2O3, can be modelled by fractional crystallization of 40%-50% plagioclase, olivine, and augite. A mildly calc-alkalic evolution trend, with decreasing FeO, increasing Al2O3, and nearly constant TiO2, can be modelled by 8%-12% olivine fractionation. Mineralogical evidence suggests that these differences may be related to the effect of small amounts of water during crystallization of the calc-alkalic suites. The tholeiitic suites occur in the lower parts of the drill cores, while the calc-alkalic suites occur in the upper parts. This suggests a complex tectonic and magmatic evolution, perhaps reflecting a transition between calc-alkalic magmatism related to subduction zone activity and tholeiitic magmatism related to back-arc spreading. Furthermore, any magmatic model must be able to account for the range in parental magmas from depleted to enriched throughout the tectonic history of the Yamato Basin.

Physicochemical properties and 3D geostatistical simulations of the Herten and the Descalvado aquifer analogs

This data set presents a comprehensive characterisation of the sedimentary structures from important groundwater hosting formations in Germany (Herten aquifer analog) and Brazil (Descalvado aquifer analog). Multiple 2-D outcrop faces are described in terms of hydraulic, thermal and chemical properties and interpolated in 3D using stochastic techniques. For each aquifer analog, multiple 3D realisations of the facies heterogeneity are provided using different stochastic simulations settings. These are unique analogue data sets that can be used by the wider community to implement approaches for characterising aquifer formations.

Grain size analysis of sediment core CRP-3 from the Ross Sea, Antarctica

Grain-size analyses by sieve and Sedigraph are presented for 115 samples of core from CRP-3, 12 km off the coast of south Victoria Land. The data provide a useful check on visual core descriptions. The geographic setting for the strata sampled, some 790 m of early Oligocene nearshore marine sediments with a persistent glacial influence, is reviewed, and sediment textures interpreted in that context. Sand textures from the CRP-3 samples in the lower part of the core suggest that deposition was initially primarily wave-dominated, but that at times the influence of the waves was over-ridden by episodes of rapid sedimentation. Sedimentary cycles, recognised in the visual description of the core above 485 mbsf, show an increasing proportion of mudstone in the middle of each cycle above 330 mbsf that is interpreted to record periodic sedimentation in deeper water. Sandstone textures in the lower and upper parts of each cycle are interpreted to record departure from and return to shoreface deposition with changes in sea level. Mudstone textures above 176 mbsf indicate sedimentation below wave base. Many of the textures in both sand and mud samples show the coarse 'tail' characteristic of ice-rafted debris, but others do not, indicating ice-free periods. Many sandstones below c. 200 mbsf have virtually no silt, but significant amounts of clay (6 to 17%) that is thought to be of post-depositional origin.

Grain size analysis of sediment coring site CRP-2 from the Ross Sea, Antarctica

The purpose of this note is to present results of grain size analyses from 118 samples of the CRP-2/2A core using sieve and Sedigraph techniques. The samples were selected to represent the range of facies encountered, and tend to become more widely spaced with depth. Fifteen came from the upper 27 m of Quaternary and Pliocene sediments, 62 from the early Miocene-late Oligocene strata (27 to 307 mbsf), and 41 from the early Oligocene strata beneath (307 to 624 mbsf). The results are intended to provide reference data for lithological descriptions in the core logs (Cape Roberts Science Team, 1999), and to help with facies interpretation. The analytical technique used for determining size frequency of the sand fraction in our samples (sieving) is simple, physical and widely practised for over a century. Thus it provides a useful reference point for analyses produced by other faster and more sophisticated techniques, such as the Malvern laser particle size analysis system (Woolfe et al., 2000), and estimates derived from measurements taken with down-hole logging tools (Bücker, pers. com., 1999).

Geochemistry of Cretaceous sills and lava flows of ODP Holes 129-800A and 129-802A

On the basis of their respective eruptive environments and chemical characteristics, alkalic dolerite sills from the northern Pigafetta Basin (Site 800) and tholeiitic pillow lavas from the Mariana Basin (Site 802) sampled during Ocean Drilling Program Leg 129 are considered to represent examples of the widespread mid-Cretaceous volcanic event in the western Pacific. Both groups of basic rocks feature mild, low-grade, anoxic smectite-celadonite-carbonate-pyrite alteration; late-stage oxidation is very limited in extent, with the exception of the uppermost sill unit at Site 800. The aphyric and nonvesicular Site 800 alkalic dolerite sills are all well-evolved mineralogically and chemically, being mainly of hawaiite composition, and are similar to ocean island basalts. They are characterized by high contents of incompatible elements (for example, 300-400 ppm Zr), well-fractionated rare earth element patterns ([La/Yb]N 18-21) and HIMU isotopic characters. They probably represent deep-sea, lateral, intrusive off-shoots from nearby seamounts of similar age. The olivine-plagioclase +/- clinopyroxene phyric tholeiitic pillow lavas and thin flows of Site 802 are nonvesicular and quench-textured throughout. Relative to normal-type mid-ocean ridge basalt, they are enriched in large-ion-lithophile elements, exhibit flat (unfractionated) rare earth element patterns and have distinctive (lower) Zr/Nb, Zr/Ta, La/Ta, and Hf/Th ratios. Overall they are compositionally and isotopically similar to the mid-Cretaceous tholeiites of the Nauru basin and the Ontong-Java and Manihiki plateaus. The Site 802 tholeiites differ from the thickened crustal segments of the oceanic plateaus, however, in apparently representing only a thin veneer over the local basement in an off-axis environment.

Investigation of Pliocene foraminifera on sediment core CRP-1 from the Ross Sea, Antarctica

Mixed assemblages of Pliocene and Quaternary foraminifera occur within the Quaternary succession of the CRP-1 drillhole. Pliocene foraminifera are not present in the lowermost Unit 4.1. are rare in Unit 3.1 and 2.3, are relatively common in Units 2.2 and 2.1, and are absent in Unit 1.1. Fifteen and twelve species were documented in two of the samples from Units 2.2 and 2.1 respectively. A census count of foraminifera in a sample at 26.89 mbsf (Unit 2.2) indicated that 39% of the tests were from a Pliocene source, with the remaining 61% tests assigned to the in situ Quaternary assemblage. There appears to be a close correlation between the stratigraphic distribution of ice-rafted sediments and the test number and diversity of Pliocene taxa. It is concluded that Pliocene assemblages were not derived from submarine outcrops on Roberts Ridge, but are more likely to have been rafted to the site via major trunk valley drainage systems such as operated within the Mackay and Ferrar glacial valleys. The co-occurrence of marine biota (including foraminifera), fossil wood, pollen, and igneous clasts in the Quaternary succession of CRP-l, points to the marine and terrestrial facies of the Pliocene Sirius Group as a likely source. A major episode of erosion and transport of sediment into the offshore marine basins at about ~1 Ma may have been triggered by dynamism in the ice sheet-glacier system, an episode of regional uplift in the Transantarctic Mountains, sea level oscillations and associated changes in the land-to-sea drainage baselines, or some combination of these factors.

Preliminary results and documentation of sediment core CRP-1 from the Ross Sea off Cape Roberts, Antarctica

The first hole of the Cape Roberts Project, CRP-1, was drilled in October, 1997, to a depth of 148 metres below the sea floor (mbsf) before being terminated unexpectedly the loss of fast sea-ice seaward of the rig following a severe storm. The site lies in 150 m of water at 77.008°S and 163.755°E, 16 km off Cape Roberts. This part of the report outlines the geologic setting, a gently tilted sequence near the margin of the Victoria Land Basin, and describes the history of the growth of sea ice, which provided the drilling platform, as well as the history of the drilling itself. Core recovery was around 77% in soft and brittle strata to 100 m and 98% below that. The sequence was found to comprise a Quaternary glacigenic interval down to 43.55 mbsf and below this an early Miocene interval that was also glacigenic. Core properties that were studied include fracture patterns, porosity, sonic velocity and magnetic susceptibility. Velocity in particular was useful in relating the cored sequence to the regional seismic stratigraphy. A preliminary assessment suggests that the bottom of the hole is 15 m short of the boundary between seismic sequences V3 and V4. Analytical facilities new to the Antarctic and used for processing samples for the project are described here and include a bench top palynological processing system and a palaeomagnetic laboratory. The core management and sampling system, which recorded over 2000 samples, is also outlined.