Mean chemical characteristics and volcano signatures of ice cores from Belukha, the Everest and three Svalbard sites

Ice cores from outside the Greenland and Antarctic ice sheets are difficult to date because of seasonal melting and multiple sources (terrestrial, marine, biogenic and anthropogenic) of sulfates deposited onto the ice. Here we present a method of volcanic sulfate extraction that relies on fitting sulfate profiles to other ion species measured along the cores in moving windows in log space. We verify the method with a well dated section of the Belukha ice core from central Eurasia. There are excellent matches to volcanoes in the preindustrial, and clear extraction of volcanic peaks in the post-1940 period when a simple method based on calcium as a proxy for terrestrial sulfate fails due to anthropogenic sulfate deposition. We then attempt to use the same statistical scheme to locate volcanic sulfate horizons within three ice cores from Svalbard and a core from Mount Everest. Volcanic sulfate is <5% of the sulfate budget in every core, and differences in eruption signals extracted reflect the large differences in environment between western, northern and central regions of Svalbard. The Lomonosovfonna and Vestfonna cores span about the last 1000 years, with good extraction of volcanic signals, while Holtedahlfonna which extends to about AD1700 appears to lack a clear record. The Mount Everest core allows clean volcanic signal extraction and the core extends back to about AD700, slightly older than a previous flow model has suggested. The method may thus be used to extract historical volcanic records from a more diverse geographical range than hitherto.

Biostratigraphic datums and magnetic properties of ODP Leg 178 sites

During Ocean Drilling Program (ODP) Leg 178, eight holes were drilled at three sites (1095, 1096, and 1101) on the continental rise along the western Antarctic Peninsula. The rise sediments proved to be good paleomagnetic recorders and provided continuous magnetostratigraphic records at all three sites. Biosiliceous microfossils, particularly diatoms and radiolarians, were present in the upper Miocene through lower Pliocene sections. In the upper Pliocene to Pleistocene sections, biosiliceous microfossils were rare but calcareous nannofossils and foraminifers were present. This paper summarizes the biostratigraphy and magnetostratigraphy of Leg 178 continental rise sites and is the first attempt at direct calibration of Antarctic biostratigraphic events to the geomagnetic polarity timescale in the Pacific sector of the Southern Ocean.

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.

Time series of ZOOSCAN images and results of zooplankton samples at Villefranche from 1966 to 2003

ZooScan with ZooProcess and Plankton Identifier (PkID) software is an integrated analysis system for acquisition and classification of digital zooplankton images from preserved zooplankton samples. Zooplankton samples are digitized by the ZooScan and processed by ZooProcess and PkID in order to detect, enumerate, measure and classify the digitized objects. Here we present a semi-automatic approach that entails automated classification of images followed by manual validation, which allows rapid and accurate classification of zooplankton and abiotic objects. We demonstrate this approach with a biweekly zooplankton time series from the Bay of Villefranche-sur-mer, France. The classification approach proposed here provides a practical compromise between a fully automatic method with varying degrees of bias and a manual but accurate classification of zooplankton. We also evaluate the appropriate number of images to include in digital learning sets and compare the accuracy of six classification algorithms. We evaluate the accuracy of the ZooScan for automated measurements of body size and present relationships between machine measures of size and C and N content of selected zooplankton taxa. We demonstrate that the ZooScan system can produce useful measures of zooplankton abundance, biomass and size spectra, for a variety of ecological studies.

Physical properties and sedimentology on core GeoB1510-1

Detailed information about the sediment properties and microstructure can be provided through the analysis of digital ultrasonic P wave seismograms recorded automatically during full waveform core logging. The physical parameter which predominantly affects the elastic wave propagation in water-saturated sediments is the P wave attenuation coefficient. The related sedimentological parameter is the grain size distribution. A set of high-resolution ultrasonic transmission seismograms (ca. 50-500 kHz), which indicate downcore variations in the grain size by their signal shape and frequency content, are presented. Layers of coarse-grained foraminiferal ooze can be identified by highly attenuated P waves, whereas almost unattenuated waves are recorded in fine-grained areas of nannofossil ooze. Color-encoded pixel graphics of the seismograms and instantaneous frequencies present full waveform images of the lithology and attenuation. A modified spectral difference method is introduced to determine the attenuation coefficient and its power law a = kfn. Applied to synthetic seismograms derived using a "constant Q" model, even low attenuation coefficients can be quantified. A downcore analysis gives an attenuation log which ranges from ca. 700 dB/m at 400 kHz and a power of n = 1-2 in coarse-grained sands to few decibels per meter and n ? 0.5 in fine-grained clays. A least squares fit of a second degree polynomial describes the mutual relationship between the mean grain size and the attenuation coefficient. When it is used to predict the mean grain size, an almost perfect coincidence with the values derived from sedimentological measurements is achieved.

Physical properties, velocity and attenuation characteristics, and mineralogy of gabbros from ODP Hole 176-735B

Laboratory compressional wave (Vp) and shear wave (Vs) velocities were measured as a function of confining pressure for the gabbros from Hole 735B and compared to results from Leg 118. The upper 500 m of the hole has a Vp mean value of 6895 m/s measured at 200 MPa, and at 500 meters below seafloor (mbsf), Vp measurements show a mean value of 7036 m/s. Vs mean values in the same intervals are 3840 m/s and 3857 m/s, respectively. The mean Vp and Vs values obtained from log data in the upper 600 m are 6520 and 3518 m/s, respectively. These results show a general increase in velocity with depth and the velocity gradients estimate an upper mantle depth of 3.32 km. This value agrees with previous work based on dredged samples and inversion of rare element concentrations in basalts dredged from the conjugate site to the north of the Atlantis Bank. Laboratory measurements show Vp anisotropy ranging between 0.4% and 8.8%, with the majority of the samples having values less than 3.8%. Measurements of velocity anisotropy seem to be associated with zones of high crystal-plastic deformation with predominant preferred mineral orientations of plagioclase, amphiboles, and pyroxenes. These findings are consistent with results on gabbros from the Hess Deep area and suggest that plastic deformation may play an important role in the seismic properties of the lower oceanic crust. In contrast to ophiolite studies, many of the olivine gabbros show a small degree of anisotropy. Log derived Vs anisotropy shows an average of 5.8% for the upper 600 m of Hole 735B and tends to decrease with depth where the overburden pressure and the age of the crustal section suggests closure of cracks and infilling of fractures by alteration minerals. Overall the results indicate that the average shear wave splitting in Hole 735B might be influenced by preferred structural orientations and the average value of shear wave splitting may not be a maximum because structural dips are <90°. The maximum fast-wave orientation values could be influenced by structural features striking slightly oblique to this orientation or by near-field stress concentrations. However, flexural wave dispersion analyses have not been performed to confirm this hypothesis or to indicate to what extent the near-field stresses may be influencing shear wave propagation. Acoustic impedance contrasts calculated from laboratory and logging data were used to generate synthetic seismograms that aid in the interpretation of reflection profiles. Several prominent reflections produced by these calculations suggest that Fe-Ti oxides and shear zones may contribute to the reflective nature of the lower oceanic crust. Laboratory velocity attenuation (Q) measurements from below 500 m have a mean value of 35.1, which is consistent with previous vertical seismic profile (VSP) and laboratory measurements on the upper 500 m.

Iberian Margin 420 kyr geochemical and color variations record

Here we report 420 kyr long records of sediment geochemical and color variations from the southwestern Iberian Margin. We synchronized the Iberian Margin sediment record to Antarctic ice cores and speleothem records on millennial time scales and investigated the phase responses relative to orbital forcing of multiple proxy records available from these cores. Iberian Margin sediments contain strong precession power. Sediment "redness" (a* and 570-560 nm) and the ratio of long-chain alcohols to n-alkanes (C26OH/(C26OH + C29)) are highly coherent and in-phase with precession. Redder layers and more oxidizing conditions (low alcohol ratio) occur near precession minima (summer insolation maxima). We suggest these proxies respond rapidly to low-latitude insolation forcing by wind-driven processes (e.g., dust transport, upwelling, precipitation). Most Iberian Margin sediment parameters lag obliquity maxima by 7-8 ka, indicating a consistent linear response to insolation forcing at obliquity frequencies driven mainly by high-latitude processes. Although the lengths of the time series are short (420 ka) for detecting 100 kyr eccentricity cycles, the phase relationships support those obtained by Shackleton []. Antarctic temperature and the Iberian Margin alcohol ratios (C26OH/(C26OH + C29)) lead eccentricity maxima by 6 kyr, with lower ratios (increased oxygenation) occurring at eccentricity maxima. CO2, CH4, and Iberian SST are nearly in phase with eccentricity, and minimum ice volume (as inferred from Pacific d18Oseawater) lags eccentricity maxima by 10 kyr. The phase relationships derived in this study continue to support a potential role of the Earth's carbon cycle in contributing to the 100 kyr cycle.

Geochemical composition and crystallization conditions of melts from the Mid Atlantic Ridge near Ascension Island (7-11°S)

Phase equilibria simulations were performed on naturally quenched basaltic glasses to determine crystallization conditions prior to eruption of magmas at the Mid-Atlantic Ridge (MAR) east of Ascension Island (7°11°S).The results indicate that midocean ridge basalt (MORB) magmas beneath different segments of the MAR have crystallized over a wide range of pressures (100-900MPa). However, each segment seems to have a specific crystallization history. Nearly isobaric crystallization conditions (100-300MPa) were obtained for the geochemically enriched MORB magmas of the central segments, whereas normal (N)-MORB magmas of the bounding segments are characterized by polybaric crystallization conditions (200-900MPa). In addition, our results demonstrate close to anhydrous crystallization conditions of N-MORBs, whereas geochemically enriched MORBs were successfully modeled in the presence of 0.4-1wt% H2O in the parental melts.These estimates are in agreement with direct (Fourier transform IR) measurements of H2O abundances in basaltic glasses and melt inclusions for selected samples. Water contents determined in the parental melts are in the range 0.04-0.09 and 0.30-0.55 wt% H2O for depleted and enriched MORBs, respectively. Our results are in general agreement (within ±200MPa) with previous approaches used to evaluate pressure estimates in MORB. However, the determination of pre-eruptive conditions of MORBs, including temperature and water content in addition to pressure, requires the improvement of magma crystallization models to simulate liquid lines of descent in the presence of small amounts of water. KEY WORDS: MORB; Mid-Atlantic Ridge; depth of crystallization; water abundances; phase equilibria calculations; cotectic crystallization; pressure estimates; polybaric fractionation

Stable oxygen isotope, UK'37, productivity, CaCO3 and XRF records from sediment core MD06-3075

Proxy records of hydrologic variability in the West Pacific Warm Pool (WPWP) have revealed wide-scale changes in past convective activity in response to orbital and sub-orbital climate forcings. However, attributing proxy responses to regional changes in WPWP hydrology versus local variations in precipitation requires independent records linking the terrestrial and marine realms. We present high-resolution stable isotope, UK'37 sea-surface temperature, X-ray fluorescence (XRF) core scanning and coccolithophore-derived paleoproductivity records covering the past 120 ka from International Marine Global Change (IMAGES) Program Core MD06-3075 (6°29' N, 125°50' E, water depth 1878 m), situated in the Davao Gulf on the southern side of Mindanao. XRF-derived log(Fe/Ca) records provide a robust proxy for runoff-driven sedimentary discharge from Mindanao, whilst past changes in local productivity are associated with variable freshwater runoff and stratification of the surface layer. Significant precessional-scale variability in sedimentary discharge occurred during Marine Isotope Stage (MIS) 5, with peaks in discharge contemporaneous with Northern Hemisphere summer insolation minima. We attribute these changes to the latitudinal migration of the Intertropical Convergence Zone (ITCZ) over the WPWP together with variability in the strength of the Walker circulation acting on precessional timescales. Between 60 and 15 ka sedimentary discharge at Mindanao was muted, displaying little orbital- or millennial-scale variability, likely in response to weakened precessional insolation forcing and lower sea level driving increased subsidence of air masses over the exposed Sunda Shelf. These results highlight the high degree of local variability in the precipitation response to past climate changes in the WPWP.

Grain size analysis of sediment core GeoB1510-2

Detailed information about the sediment properties and microstructure can be provided through the analysis of digital ultrasonic P wave seismograms recorded automatically during full waveform core logging. The physical parameter which predominantly affects the elastic wave propagation in water-saturated sediments is the P wave attenuation coefficient. The related sedimentological parameter is the grain size distribution. A set of high-resolution ultrasonic transmission seismograms (-50-500 kHz), which indicate downcore variations in the grain size by their signal shape and frequency content, are presented. Layers of coarse-grained foraminiferal ooze can be identified by highly attenuated P waves, whereas almost unattenuated waves are recorded in fine-grained areas of nannofossil ooze. Color -encoded pixel graphics of the seismograms and instantaneous frequencies present full waveform images of the lithology and attenuation. A modified spectral difference method is introduced to determine the attenuation coefficient and its power law a = kF. Applied to synthetic seismograms derived using a "constant Q" model, even low attenuation coefficients can be quantified. A downcore analysis gives an attenuation log which ranges from -700 dB/m at 400 kHz and a power of n=1-2 in coarse-grained sands to few decibels per meter and n :s; 0.5 in fine-grained clays. A least squares fit of a second degree polynomial describes the mutual relationship between the mean grain size and the attenuation coefficient. When it is used to predict the mean grain size, an almost perfect coincidence with the values derived from sedimentological measurements is achieved.

3D models created from ROV videos and their corresponding subtransect dimensions

Underwater video transects have become a common tool for quantitative analysis of the seafloor. However a major difficulty remains in the accurate determination of the area surveyed as underwater navigation can be unreliable and image scaling does not always compensate for distortions due to perspective and topography. Depending on the camera set-up and available instruments, different methods of surface measurement are applied, which make it difficult to compare data obtained by different vehicles. 3-D modelling of the seafloor based on 2-D video data and a reference scale can be used to compute subtransect dimensions. Focussing on the length of the subtransect, the data obtained from 3-D models created with the software PhotoModeler Scanner are compared with those determined from underwater acoustic positioning (ultra short baseline, USBL) and bottom tracking (Doppler velocity log, DVL). 3-D model building and scaling was successfully conducted on all three tested set-ups and the distortion of the reference scales due to substrate roughness was identified as the main source of imprecision. Acoustic positioning was generally inaccurate and bottom tracking unreliable on rough terrain. Subtransect lengths assessed with PhotoModeler were on average 20% longer than those derived from acoustic positioning due to the higher spatial resolution and the inclusion of slope. On a high relief wall bottom tracking and 3-D modelling yielded similar results. At present, 3-D modelling is the most powerful, albeit the most time-consuming, method for accurate determination of video subtransect dimensions.

(Table 1) Crystal sizes of measured samples

Due to experimental difficulties grain size distributions of gas hydrate crystallites are largely unknown in natural samples. For the first time, we were able to determine grain size distributions of six natural gas hydrates for samples retrieved from the Gulf of Mexico and from Hydrate Ridge offshore Oregon from varying depths. High-energy synchrotron radiation provides high photon fluxes as well as high penetration depth and thus allows for investigation of bulk sediment samples. The gas hydrate crystallites appear to be (log-) normally distributed in the natural samples and to be of roughly globular shape. The mean grain sizes are in the range from 300-600 µm with a tendency for bigger grains to occur in greater depth, possibly indicating a difference in the formation age. Laboratory produced methane hydrate, starting from ice and aged for 3 weeks, shows half a log-normal curve with a mean value of ~40 µm. This one order-of-magnitude smaller grain sizes suggests that care must be taken when transposing grain-size sensitive (petro-)physical data from laboratory-made gas hydrates to natural settings.

Spatial distribution of key zooplankton species using continuous plankton recorder (CPR) data from the North Atlantic (2000-2009)

The continuous plankton recorder (CPR) survey is an upper layer plankton monitoring program that has regularly collected samples, at monthly intervals, in the North Atlantic and adjacent seas since 1946. Water from approximately 6 m depth enters the CPR through a small aperture at the front of the sampler and travels down a tunnel where it passes through a silk filtering mesh of 270 µm before exiting at the back of the CPR. The plankton filtered on the silk is analyzed in sections corresponding to 10 nautical miles (approx. 3 m**3 of seawater filtered) and the plankton microscopically identified (Richardson et al., 2006 and reference therein). In the present study we used the CPR data to investigate the current basin scale distribution of C. finmarchicus (C5-C6), C. helgolandicus (C5-C6), C. hyperboreus (C5-C6), Pseudocalanus spp. (C6), Oithona spp. (C1-C6), total Euphausiida, total Thecosomata and the presence/absence of Cnidaria and the Phytoplankton Colour Index (PCI). The PCI, which is a visual assessment of the greenness of the silk, is used as an indicator of the distribution of total phytoplankton biomass across the Atlantic basin (Batten et al., 2003). Monthly data collected between 2000 and 2009 were gridded using the inverse-distance interpolation method, in which the interpolated values were the nodes of a 2 degree by 2 degree grid. The resulting twelve monthly matrices were then averaged within the year and in the case of the zooplankton the data were log-transformed (i.e. log10 (x+1).