(Table 3) Details of field-based ice and snow profile measurements from the SIMBA 2007 experiment

Although sea-ice extent in the Bellingshausen-Amundsen (BA) seas sector of the Antarctic has shown significant decline over several decades, there is not enough data to draw any conclusion on sea-ice thickness and its change for the BA sector, or for the entire Southern Ocean. This paper presents our results of snow and ice thickness distributions from the SIMBA 2007 experiment in the Bellingshausen Sea, using four different methods (ASPeCt ship observations, downward-looking camera imaging, ship-based electromagnetic induction (EM) sounding, and in situ measurements using ice drills). A snow freeboard and ice thickness model generated from in situ measurements was then applied to contemporaneous ICESat (satellite laser altimetry) measured freeboard to derive ice thickness at the ICESat footprint scale. Errors from in situ measurements and from ICESat freeboard estimations were incorporated into the model, so a thorough evaluation of the model and uncertainty of the ice thickness estimation from ICESat are possible. Our results indicate that ICESat derived snow freeboard and ice thickness distributions (asymmetrical unimodal tailing to right) for first-year ice (0.29 ± 0.14 m for mean snow freeboard and 1.06 ± 0.40 m for mean ice thickness), multi-year ice (0.48 ± 0.26 and 1.59 ± 0.75 m, respectively), and all ice together (0.42 ± 0.24 and 1.38 ± 0.70 m, respectively) for the study area seem reasonable compared with those values from the in situ measurements, ASPeCt observations, and EM measurements. The EM measurements can act as an appropriate supplement for ASPeCt observations taken hourly from the ship's bridge and provide reasonable ice and snow distributions under homogeneous ice conditions. Our proposed approaches: (1) of using empirical equations relating snow freeboard to ice thickness based on in situ measurements and (2) of using isostatic equations that replace snow depth with snow freeboard (or empirical equations that convert freeboard to snow depth), are efficient and important ways to derive ice thickness from ICESat altimetry at the footprint scale for Antarctic sea ice. Spatial and temporal snow and ice thickness from satellite altimetry for the BA sector and for the entire Southern Ocean is therefore possible.

An individual based model spatial distribution, production and biomass for Calanus finmarchicus in the Norwegian Sea

The copepod Calanus finmarchicus is the dominant species of the meso-zooplankton in the Norwegian Sea, and constitutes an important link between the phytoplankton and the higher trophic levels in the Norwegian Sea food chain. An individualbased model for C. finmarchicus, based on super-individuals and evolving traits for behaviour, stages, etc., is two-way coupled to the NORWegian ECOlogical Model system (NORWECOM). One year of modelled C. finmarchicus spatial distribution, production and biomass are found to represent observations reasonably well. High C. finmarchicus abundance is found along the Norwegian shelf-break in the early summer, while the overwintering population is found along the slope and in the deeper Norwegian Sea basins. The timing of the spring bloom is generally later than in the observations. Annual Norwegian Sea production is found to be 29 million tonnes of carbon and a production to biomass (P/B) ratio of 4.3 emerges. Sensitivity tests show that the modelling system is robust to initial values of behavioural traits and with regards to the number of super-individuals simulated given that this is above about 50,000 individuals. Experiments with the model system indicate that it provides a valuable tool for studies of ecosystem responses to causative forces such as prey density or overwintering population size. For example, introducing C. finmarchicus food limitations reduces the stock dramatically, but on the other hand, a reduced stock may rebuild in one year under normal conditions. The NetCDF file contains model grid coordinates and bottom topography.

Treeline dynamics in Siberia under changing climates as inferred from an individual-based model for Larix, supplementary material

Siberian boreal forests are expected to expand northwards in the course of global warming. However, processes of the treeline ecotone transition, as well astiming and related climate feedbacks are still not understood. Here, we present 'Larix Vegetation Simulator' LAVESI, an individual-based spatially-explicit model that can simulate Larix gmelinii (RUPR.) RUPR. stand dynamics in an attempt to improve our understanding about past and future treeline movements under changing climates. The relevant processes (growth, seed production and dispersal, establishment and mortality) are incorporated and adjusted to observation data mainly gained from the literature. Results of a local sensitivity analysis support the robustness of the model's parameterization by giving relatively small sensitivity values. We tested the model by simulating tree stands under modern climate across the whole Taymyr Peninsula, north-central Siberia (c. 64-80° N; 92-119° E). We find tree densities similar to observed forests in the northern to mid-treeline areas, but densities are overestimated in the southern parts of the simulated region. Finally, from a temperature-forcing experiment, we detect that the responses of tree stands lag the hypothetical warming by several decades, until the end of 21st century. With our simulation experiments we demonstrate that the newly-developed model captures the dynamics of the Siberian latitudinal treeline.

(Table T1) Physical properties of minicores from ODP Leg 193 sites, PACMANUS field

Permeability of the ocean crust is one of the most crucial parameters for constraining submarine fluid flow systems. Active hydrothermal fields are dynamic areas where fluid flow strongly affects the geochemistry and biology of the surrounding environment. There have been few permeability measurements in these regions, especially in felsic-hosted hydrothermal systems. We present a data set of 38 permeability and porosity measurements from the PACMANUS hydrothermal field, an actively venting, felsic hydrothermal field in the eastern Manus Basin. Permeability was measured using a complex transient method on 2.54-cm minicores. Permeability varies greatly between the samples, spanning over five orders of magnitude. Permeability decreases with both depth and decreasing porosity. When the alteration intensity of individual samples is considered, relationships between depth and porosity and permeability become more clearly defined. For incompletely altered samples (defined as >5% fresh rock), permeability and porosity are constant with depth. For completely altered samples (defined as <5% fresh rock), permeability and porosity decrease with depth. On average, the permeability values from the PACMANUS hydrothermal field are greater than those in other submarine environments using similar core-scale laboratory measurements; the average permeability, 4.5 x 10-16 m**2, is two to four orders of magnitude greater than in other areas. Although the core-scale permeability is higher than in other seafloor environments, it is still too low to obtain the fluid velocities observed in the PACMANUS hydrothermal field based on simplified analytical calculations. It is likely that core-scale permeability measurements are not representative of bulk rock permeability of the hydrothermal system overall, and that the latter is predominantly fracture controlled.

Geochemistry of microtektites from Victoria Land Transantarctic Mountains

We extended the petrographic and geochemical dataset for the recently discovered Transantarctic Mountain microtektites in order to check our previous claim that they are related to the Australasian strewn field. Based on color and composition, the 465 microtektites so far identified include two groups of transparent glass spheres less than ca. 800 µm in diameter: the most abundant pale-yellow, or normal, microtektites, and the rare pale-green, or high-Mg, microtektites. The major element composition of normal microtektites determined through electron microprobe analysis is characterized by high contents of silica (SiO2 = 71.5 ± 3.6 (1 sigma) wt%) and alumina (Al2O3 = 15.5 ± 2.2 (1 sigma) wt%), low total alkali element contents (0.50-1.85 wt%), and MgO abundances <6 wt%. The high-Mg microtektites have a distinctly higher MgO content >10 wt%. Transantarctic Mountain microtektites contain rare silica-rich (up to 93 wt% SiO2) glassy inclusions similar to those found in two Australasian microtektites analyzed here for comparison. These inclusions are interpreted as partially digested, lechatelierite-like inclusions typically found in tektites and microtektites. The major and trace element (by laser ablation - inductively coupled plasma - mass spectrometry) abundance pattern of the Transantarctic Mountain microtektites matches the average upper continental crust composition for most elements. Major deviations include a strong to moderate depletion in volatile elements including Pb, Zn, Na, K, Rb, Sr and Cs, as a likely result of severe volatile loss during the high temperature melting and vaporization of crustal target rocks. The normal and high-Mg Transantarctic Mountain microtektites have compositions similar to the most volatile-poor normal and high-Mg Australasian microtektites reported in the literature. Their very low H2O and B contents (by secondary ion mass spectrometry) of 85 ± 58 (1 sigma) ?g/g and 0.53 ± 0.21 ?g/g, respectively, evidence the extreme volatile loss characteristically observed in tektites. The Sr and Nd isotopic compositions of multigrain samples of Transantarctic Mountain microtektites are 87Sr/86Sr ~ 0.71629 and 143Nd/144Nd ~ 0.51209, and fall into the Australasian tektite compositional field. The Nd model age calculated with respect to the chondritic uniform reservoir (CHUR) is TNdCHUR ~ 1.1 Ga, indicating a Meso-Proterozoic crustal source rock, as was derived for Australasian tektites as well. Coupled with the Quaternary age from the literature, the extended dataset presented in this work strengthens our previous conclusion that Transantarctic Mountain microtektites represent a major southward extension of the Australasian tektite/microtektite strewn field. Furthermore, the significant depletion in volatile elements (i.e., Pb, B, Na, K, Zn, Rb, Sr and Cs) of both normal and high-Mg Transantarctic Mountain microtektites relative to the Australasian ones provide us with further confirmation of a possible relationship between high temperature-time regimes in the microtektite-forming process and ejection distance.

Stable isotope ratios of foraminifera and sea surface temperature estimation of sediment core W2804A-14

Carbon isotopically based estimates of CO2 levels have been generated from a record of the photosynthetic fractionation of 13C (epsilon p) in a central equatorial Pacific sediment core that spans the last ~255 ka. Contents of 13C in phytoplanktonic biomass were determined by analysis of C37 alkadienones. These compounds are exclusive products of Prymnesiophyte algae which at present grow most abundantly at depths of 70-90 m in the central equatorial Pacific. A record of the isotopic compostion of dissolved CO2 was constructed from isotopic analyses of the planktonic foraminifera Neogloboquadrina dutertrei, which calcifies at 70-90 m in the same region. Values of epsilon p, derived by comparison of the organic and inorganic delta values, were transformed to yield concentrations of dissolved CO2 (c e) based on a new, site-specific calibration of the relationship between epsilon p and c e. The calibration was based on reassessment of existing epsilon p versus c e data, which support a physiologically based model in which epsilon p is inversely related to c e. Values of PCO2, the partial pressure of CO2 that would be in equilibrium with the estimated concentrations of dissolved CO2, were calculated using Henry's law and the temperature determined from the alkenone-unsaturation index UK 37. Uncertainties in these values arise mainly from uncertainties about the appropriateness (particularly over time) of the site-specific relationship between epsilon p and 1/c e. These are discussed in detail and it is concluded that the observed record of epsilon p most probably reflects significant variations in Delta pCO2, the ocean-atmosphere disequilibrium, which appears to have ranged from ~110 µatm during glacial intervals (ocean > atmosphere) to ~60 µatm during interglacials. Fluxes of CO2 to the atmosphere would thus have been significantly larger during glacial intervals. If this were characteristic of large areas of the equatorial Pacific, then greater glacial sinks for the equatorially evaded CO2 must have existed elsewhere. Statistical analysis of air-sea pCO2 differences and other parameters revealed significant (p < 0.01) inverse correlations of Delta pCO2 with sea surface temperature and with the mass accumulation rate of opal. The former suggests response to the strength of upwelling, the latter may indicate either drawdown of CO2 by siliceous phytoplankton or variation of [CO2]/[Si(OH)4] ratios in upwelling waters.

Biostratigraphy, magnetostratigraphy, granulometry, isotopes and clay mineral composition of ODP Site 188-1165

During Ocean Drilling Program Leg 188 to Prydz Bay, East Antarctica, several of the shipboard scientists formed the High-Resolution Integrated Stratigraphy Committee (HiRISC). The committee was established in order to furnish an integrated data set from the Pliocene portion of Site 1165 as a contribution to the ongoing debate about Pliocene climate and climate evolution in Antarctica. The proxies determined in our various laboratories were the following: magnetostratigraphy and magnetic properties, grain-size distributions (granulometry), near-ultraviolet, visible, and near-infrared spectrophotometry, calcium carbonate content, characteristics of foraminifer, diatom, and radiolarian content, clay mineral composition, and stable isotopes. In addition to the HiRISC samples, other data sets contained in this report are subsets of much larger data sets. We included these subsets in order to provide the reader with a convenient integrated data set of Pliocene-Pleistocene strata from the East Antarctic continental margin. The data are presented in the form of 14 graphs (in addition to the site map). Text and figure captions guide the reader to the original data sets. Some preliminary interpretations are given at the end of the manuscript.

Selective mortality associated with variation in CO2 tolerance in a marine fish

Predicted future CO2 levels can affect reproduction, growth, and behaviour of many marine organisms. However, the capacity of species to adapt to predicted changes in ocean chemistry is largely unknown. We used a unique field-based experiment to test for differential survival associated with variation in CO2 tolerance in a wild population of coral-reef fishes. Juvenile damselfish exhibited variation in their response to elevated (700 µatm) CO2 when tested in the laboratory and this influenced their behaviour and risk of mortality in the wild. Individuals that were sensitive to elevated CO2 were more active and move further from shelter in natural coral reef habitat and, as a result, mortality from predation was significantly higher compared with individuals from the same treatment that were tolerant of elevated CO2. If individual variation in CO2 tolerance is heritable, this selection of phenotypes tolerant to elevated CO2 could potentially help mitigate the effects of ocean acidification.

Stable isotope data and calcification temperature from planktic foraminifera in sediment core BOFS31/1K

Determining the past record of temperature and salinity of ocean surface waters is essential for understanding past changes in climate, such as those which occur across glacial-interglacial transitions. As a useful proxy, the oxygen isotope composition (delta18O) of calcite from planktonic foraminifera has been shown to reflect both surface temperature and seawater delta18O, itself an indicator of global ice volume and salinity (Shackleton, 1974; Rostek et al., 1993, doi:10.1038/364319a0). In addition, magnesium/calcium (Mg/Ca) ratios in foraminiferal calcite show a temperature dependence (Nürnberg, 1995, doi:10.2113/gsjfr.25.4.350; Nürnberg et al., 1996, doi:10.1016/0016-7037(95)00446-7; Lea et al., 1999, doi:10.1016/S0016-7037(99)00197-0) due to the partitioning of Mg during calcification. Here we demonstrate, in a field-based calibration experiment, that the variation of Mg/Ca ratios with temperature is similar for eight species of planktonic foraminifera (when accounting for Mg dissolution effects). Using a multi-species record from the Last Glacial Maximum in the North Atlantic Ocean we found that past temperatures reconstructed from Mg/Ca ratios followed the two other palaeotemperature proxies: faunal abundance (CLIMAP, 1981; Mix et al., 1999, doi:10.1029/1999PA900012) and alkenone saturation (Müller et al., 1998, doi:10.1016/S0016-7037(98)00097-0 ). Moreover, combining Mg/Ca and delta18O data from the same faunal assemblage, we show that reconstructed surface water delta18O from all foraminiferal species record the same glacial-interglacial change-representing changing hydrography and global ice volume. This reinforces the potential of this combined technique in probing past ocean-climate interactions.