Geochemistry of the oldest Atlantic oceanic crust

Controversy has surrounded the issue of whether mantle plume activity was responsible for Pangaean continental rifting and massive flood volcanism (resulting in the Central Atlantic Magmatic Province or CAMP, emplaced around 200 Ma) preceding the opening of the central Atlantic Ocean in the Early Mesozoic. Our new Sr-Nd-Pb isotopic and trace element data for the oldest basalts sampled from central Atlantic oceanic crust by deep-sea drilling show that oceanic crust generated from about 160 to 120 Ma displays clear isotopic and chemical signals of plume contamination (e.g., 87Sr/86Sr(i) = 0.7032-0.7036, epsilonNd(t) =+6.2 to +8.2, incompatible element patterns with positive Nb anomalies), but these signals are muted or absent in crust generated between 120 and 80 Ma, which resembles young Atlantic normal mid-ocean ridge basalt. The plume-affected pre-120 Ma Atlantic crustal basalts are isotopically similar to lavas from the Ontong Java Plateau, and may represent one isotopic end-member for CAMP basalts. The strongest plume signature is displayed near the center of CAMP magmatism but the hotspots presently located nearest this location in the mantle reference frame do not appear to be older than latest Cretaceous and are isotopically distinct from the oldest Atlantic crust. The evidence for widespread plume contamination of the nascent Atlantic upper mantle, combined with a lack of evidence for a long-lived volcanic chain associated with this plume, leads us to propose that the enriched signature of early Atlantic crust and possibly the eruption of the CAMP were caused by a relatively short-lived, but large volume plume feature that was not rooted at a mantle boundary layer. Such a phenomenon has been predicted by recent numerical models of mantle circulation.

Age determination of sediment core SO75-26KL, Portuguese margin

A benthic isotope record has been measured for core SO75-26KL from the upper Portuguese margin (1099 m water depth) to monitor the response of thermohaline overturn in the North Atlantic during Heinrich events. Evaluating benthic delta18O in TS diagrams in conjunction with equilibrium deltac fractionation implies that advection of Mediterranean outflow water (MOW) to the upper Portuguese margin was significantly reduced during the last glacial (< 15% compared to 30% today). The benthic isotope record along core SO75-26KL therefore primarily monitors variability of glacial North Atlantic conveyor circulation. The 14C-accelerator mass spectrometry ages of 13.54±.07 and 20.46±.12 ka for two ice-rafted detritus (IRD) layers in the upper core section and an interpolated age of 36.1 ka for a third IRD layer deeper in the core are in the range of published 14C ages for Heinrich events H1, H2, and H4. Marked depletion of benthic delta13C by 0.7-1.1 per mil during the Heinrich events suggests reduced thermohaline overturn in the North Atlantic during these events. Close similarity between meltwater patterns (inferred from planktonic delta18O) at Site 609 and ventilation patterns (inferred from benthic delta13C) in core SO75-26KL implies coupling between thermohaline overturn and surface forcing, as is also suggested by ocean circulation models. Benthic delta13C starts to decrease 1.5-2.5 kyr before Heinrich events Hl and H4, fully increased values are reached 1.5-3 kyr after the events, indicating a successive slowdown of thermohaline circulation well before the events and resumption of the conveyor's full strength well after the events. Benthic delta13C changes in the course of the Heinrich events show subtle maxima and minima suggesting oscillatory behavior of thermohaline circulation, a distinct feature of thermohaline instability in numerical models. Inferrred gradual spin-up of thermohaline circulation after Hl and H4 is in contrast to abrupt wanning in the North Atlantic region that is indicated by sudden increases in Greenland ice core delta18O and in marine faunal records from the northern North Atlantic. From this we infer that thermohaline circulation can explain only in part the rapid climatic oscillations seen in glacial sections of the Greenland ice core record.

(Table 1) Sea ice and snow characteristics and heat fluxes observed during R/V Aurora Australis cruise to East Antarctica in September/October 2007

The properties of snow on East Antarctic sea ice off Wilkes Land were examined during the Sea Ice Physics and Ecosystem Experiment (SIPEX) in late winter of 2007, focusing on the interaction with sea ice. This observation includes 11 transect lines for the measurement of ice thickness, freeboard, and snow depth, 50 snow pits on 13 ice floes, and diurnal variation of surface heat flux on three ice floes. The detailed profiling of topography along the transects and the d18O, salinity, and density datasets of snow made it possible to examine the snow-sea-ice interaction quantitatively for the first time in this area. In general, the snow displayed significant heterogeneity in types, thickness (mean: 0.14 +- 0.13 m), and density (325 +- 38 kg/m**3), as reported in other East Antarctic regions. High salinity was confined to the lowest 0.1 m. Salinity and d18O data within this layer revealed that saline water originated from the surface brine of sea ice in 20% of the total sites and from seawater in 80%. From the vertical profiles of snow density, bulk thermal conductivity of snow was estimated as 0.15 W/K/m on average, only half of the value used for numerical sea-ice models. Although the upward heat flux within snow estimated with this value was significantly lower than that within ice, it turned out that a higher value of thermal conductivity (0.3 to 0.4 W/K/m) is preferable for estimating ice growth amount in current numerical models. Diurnal measurements showed that upward conductive heat flux within the snow and net long-wave radiation at the surface seem to play important roles in the formation of snow ice from slush. The detailed surface topography allowed us to compare the air-ice drag coefficients of ice and snow surfaces under neutral conditions, and to examine the possibility of the retrieval of ice thickness distribution from satellite remote sensing. It was found that overall snow cover works to enhance the surface roughness of sea ice rather than moderate it, and increases the drag coefficient by about 10%. As for thickness retrieval, mean ice thickness had a higher correlation with ice surface roughness than mean freeboard or surface elevation, which indicates the potential usefulness of satellite L-band SAR in estimating the ice thickness distribution in the seasonal sea-ice zone.

Grain entrainment and transport rates of magnetic minerals of coastal sediment samples from the Bay of Plenty, New Zealand

Heavy (magnetic & non-magnetic) minerals are found concentrated by natural processes in many fluvial, estuarine, coastal and shelf environments with a potential to form economic placer deposits. Understanding the processes of heavy mineral transport and enrichment is prerequisite to interpret sediment magnetic properties in terms of hydro- and sediment dynamics. In this study, we combine rock magnetic and sedimentological laboratory measurements with numerical 3D discrete element models to investigate differential grain entrainment and transport rates of magnetic minerals in a range of coastal environments (riverbed, mouth, estuary, beach and near-shore). We analyzed grain-size distributions of representative bulk samples and their magnetic mineral fractions to relate grain-size modes to respective transport modes (traction, saltation, suspension). Rock magnetic measurements showed that distribution shapes, population sizes and grain-size offsets of bulk and magnetic mineral fractions hold information on the transport conditions and enrichment process in each depositional environment. A downstream decrease in magnetite grain size and an increase in magnetite concentration was observed from riverine source to marine sink environments. Lower flow velocities permit differential settling of light and heavy mineral grains creating heavy mineral enriched zones in estuary settings, while lighter minerals are washed out further into the sea. Numerical model results showed that higher heavy mineral concentrations in the bed increased the erosion rate and enhancing heavy mineral enrichment. In beach environments where sediments contained light and heavy mineral grains of equivalent grain sizes, the bed was found to be more stable with negligible amount of erosion compared to other bed compositions. Heavy mineral transport rates calculated for four different bed compositions showed that increasing heavy mineral content in the bed decreased the transport rate. There is always a lag in transport between light and heavy minerals which increases with higher heavy mineral concentration in all tested bed compositions. The results of laboratory experiments were validated by numerical models and showed good agreement. We demonstrate that the presented approach bears the potential to investigate heavy mineral enrichment processes in a wide range of sedimentary settings.

Calcium concentrations and isotopic composition in carbonates and porefluids from ODP Hole 130-807A

The calcium isotopic compositions (d44Ca) of 30 high-purity nannofossil ooze and chalk and 7 pore fluid samples from ODP Site 807A (Ontong Java Plateau) are used in conjunction with numerical models to determine the equilibrium calcium isotope fractionation factor (a_s-f) between calcite and dissolved Ca2+ and the rates of post-depositional recrystallization in deep sea carbonate ooze. The value of a_s-f at equilibrium in the marine sedimentary section is 1.0000+/-0.0001, which is significantly different from the value (0.9987+/-0.0002) found in laboratory experiments of calcite precipitation and in the formation of biogenic calcite in the surface ocean. We hypothesize that this fractionation factor is relevant to calcite precipitation in any system at equilibrium and that this equilibrium fractionation factor has implications for the mechanisms responsible for Ca isotope fractionation during calcite precipitation. We describe a steady state model that offers a unified framework for explaining Ca isotope fractionation across the observed precipitation rate range of ~14 orders of magnitude. The model attributes Ca isotope fractionation to the relative balance between the attachment and detachment fluxes at the calcite crystal surface. This model represents our hypothesis for the mechanism responsible for isotope fractionation during calcite precipitation. The Ca isotope data provide evidence that the bulk rate of calcite recrystallization in freshly-deposited carbonate ooze is 30-40%/Myr, and decreases with age to about 2%/Myr in 2-3 million year old sediment. The recrystallization rates determined from Ca isotopes for Pleistocene sediments are higher than those previously inferred from pore fluid Sr concentration and are consistent with rates derived for Late Pleistocene siliciclastic sediments using uranium isotopes. Combining our results for the equilibrium fractionation factor and recrystallization rates, we evaluate the effect of diagenesis on the Ca isotopic composition of marine carbonates at Site 807A. Since calcite precipitation rates in the sedimentary column are many orders of magnitude slower than laboratory experiments and the pore fluids are only slightly oversaturated with respect to calcite, the isotopic composition of diagenetic calcite is likely to reflect equilibrium precipitation. Accordingly, diagenesis produces a maximum shift in d44Ca of +0.15? for Site 807A sediments but will have a larger impact where sedimentation rates are low, seawater circulates through the sediment pile, or there are prolonged depositional hiatuses.

Compilation of measurements from two interdisciplinary STAR-shaped surveys in the North Sea: The ZISCH - Dataset

The measurements were obtained during two North Sea wide STAR-shaped cruises during summer 1986 and winter 1987, which were performed to investigate the circulation induced transport and biologically induced pollutant transfer within the interdisciplinary research in the project "ZISCH - Zirkulation und Schadstoffumsatz in der Nordsee / Circulation and Contaminant Fluxes in the North Sea (1984-1989)". The inventory presents parameters measured on hydrodynamics, nutrient dynamics, ecosystem dynamics and pollutant dynamics in the pelagic and benthic realm. The research program had the objective of quantifying fluxes of major budgets, especially contaminants in the North Sea. In spring 1986, following the phytoplankton spring bloom, and in late winter 1987, at minimum primary production activity, the North Sea ecosystem was investigated on a station net covering the whole North Sea. The station net was shaped like a star. Sampling started in the centre, followed by the northwest section and moving counter clockwise around the North Sea following the residual currents. By this strategy, a time series was measured in the central North Sea and more synoptic data sets were obtained in the individual sections. Generally advection processes have to be considered when comparing the data from different stations. The entire sampling period lasted for more than six weeks in each cruise. Thus, a time-lag should be considered especially when comparing the data from the eastern and the western part of the central and northern North Sea, where samples were taken at the beginning and at the end of the campaign. The ZISCH investigations represented a qualitatively and quantitatively new approach to North Sea research in several respects. (1) The first simultaneous blanket coverage of all important biological, chemical and physical parameters in the entire North Sea ecosystem; (2) the first simultaneous measurements of major contaminants (metals and organohaline compounds) in the different ecosystem compartments; (3) simultaneous determinations of atmospheric inputs of momentum, energy and matter as important ecosystem boundary conditions; (4) performance of the complex measurement program during two seasons, namely the spring plankton bloom and the subsequent winter period of minimal biological activity; and (5) support of data analysis and interpretation by oceanographic and meteorological numerical models on the same scales.

Coral Sr/Ca records from Bonaire Holocene corals

Proxy reconstructions of tropical Atlantic sea surface temperature (SST) that extend beyond the period of instrumental observations have primarily focused on centennial to millennial variability rather than on seasonal to multidecadal variability. Here we present monthly-resolved records of Sr/Ca (a proxy of SST) from fossil annually-banded Diploria strigosa corals from Bonaire (southern Caribbean Sea). The individual corals provide time-windows of up to 68 years length, and the total number of 295 years of record allows for assessing the natural range of seasonal to multidecadal SST variability in the western tropical Atlantic during snapshots of the mid- to late Holocene. Comparable to modern climate, the coral Sr/Ca records reveal that mid- to late Holocene SST was characterised by clear seasonal cycles, persistent quasi-biennial and prominent interannual as well as inter- to multidecadal-scale variability. However, the magnitude of SST variations on these timescales has varied over the last 6.2 ka. The coral records show increased seasonality during the mid-Holocene consistent with climate model simulations indicating that southern Caribbean SST seasonality is induced by insolation changes on orbital timescales, whereas internal dynamics of the climate system play an important role on shorter timescales. Interannual SST variability is linked to ocean-atmosphere interactions of Atlantic and Pacific origin. Pronounced interannual variability in the western tropical Atlantic is indicated by a 2.35 ka coral, possibly related to a strengthening of the variability of the El Niño/Southern Oscillation throughout the Holocene. Prominent inter- to multidecadal SST variability is evident in the coral records and slightly more pronounced in the mid-Holocene. We finally argue that our coral data provide a target for studying Holocene climate variability on seasonal and interannual to multidecadal timescales, when using further numerical models and high-resolution proxy data.

Geochemistry of sediments beneath the Tonga-Kermadec arc

The fate of subducted sediment and the extent to which it is dehydrated and/or melted before incorporation into arc lavas has profound implications for the thermo-mechanical nature of the mantle wedge and models for crustal evolution. In order to address these issues, we have undertaken the first measurements of 10Be and light elements in lavas from the Tonga-Kermadec arc and the sediment profile at DSDP site 204 outboard of the trench. The 10Be/9Be ratios in the Tonga lavas are lower than predicted from flux models but can be explained if (a) previously estimated sediment contributions are too high by a factor of 2-10, (b) the top 1-22 m of the incoming sediment is accreted, (c) large amounts of sediment erosion are proposed, or (d) the sediment component takes several Myr longer than the subducting plate to reach the magma source region beneath Tonga. The lavas form negative Th/Be-Li/Be arrays that extend from a depleted mantle source composition to lower Th/Be and Li/Be ratios than that of the bulk sediment. Thus, these arrays are not easily explained by bulk sediment addition and, using partition coefficients derived from experiments on the in-coming sediment, we show that they are also unlikely to result from fluid released during dehydration of the sediment (or altered oceanic crust). However, partial melts of the dehydrated sediment residue formed at ~800 °C during the breakdown of amphibole +/- plagioclase and in the absence of cordierite have significantly lowered Th/Be ratios. The lava arrays can be successfully modelled as 10-15% partial melts of depleted mantle after it has been enriched by the addition of 0.2-2% of these partial melts. Phase relations suggest that this requires that the top of the subducting crust reaches temperatures of ~800 °C by the time it attains ~ 80 km depth which is in excellent agreement with the results of recent numerical models incorporating a temperature-dependent mantle viscosity. Under these conditions the wet basalt solidus is also crossed yet there is no recognisable eclogitic signal in the lavas suggesting that on-going dehydration or strong thermal gradients in the upper part of the subducting plate inhibit partialmelting of the altered oceanic crust.

GPS raw data (control points and ground control points) from the Liguria Region, Borghetto Santo Spirito, Italy

Monitoring the impact of sea storms on coastal areas is fundamental to study beach evolution and the vulnerability of low-lying coasts to erosion and flooding. Modelling wave runup on a beach is possible, but it requires accurate topographic data and model tuning, that can be done comparing observed and modeled runup. In this study we collected aerial photos using an Unmanned Aerial Vehicle after two different swells on the same study area. We merged the point cloud obtained with photogrammetry with multibeam data, in order to obtain a complete beach topography. Then, on each set of rectified and georeferenced UAV orthophotos, we identified the maximum wave runup for both events recognizing the wet area left by the waves. We then used our topography and numerical models to simulate the wave runup and compare the model results to observed values during the two events. Our results highlight the potential of the methodology presented, which integrates UAV platforms, photogrammetry and Geographic Information Systems to provide faster and cheaper information on beach topography and geomorphology compared with traditional techniques without losing in accuracy. We use the results obtained from this technique as a topographic base for a model that calculates runup for the two swells. The observed and modeled runups are consistent, and open new directions for future research.

Predicting bed form roughness: the influence of lee side angle

Flow transverse bedforms (ripples and dunes) are ubiquitous in rivers and coastal seas. Local hydrodynamics and transport conditions depend on the size and geometry of these bedforms, as they constitute roughness elements at the bed. Bedform influence on flow energy must be considered for the understanding of flow dynamics, and in the development and application of numerical models. Common estimations or predictors of form roughness (friction factors) are based mostly on data of steep bedforms (with angle-of-repose lee slopes), and described by highly simplified bedform dimensions (heights and lengths). However, natural bedforms often are not steep, and differ in form and hydraulic effect relative to idealised bedforms. Based on systematic numerical model experiments, this study shows how the hydraulic effect of bedforms depends on the flow structure behind bedforms, which is determined by the bedform lee side angle, aspect ratio and relative height. Simulations reveal that flow separation behind bedform crests and, thus, a hydraulic effect is induced at lee side angles steeper than 11 to 18° depending on relative height, and that a fully developed flow separation zone exists only over bedforms with a lee side angle steeper than 24°. Furthermore, the hydraulic effect of bedforms with varying lee side angle is evaluated and a reduction function to common friction factors is proposed. A function is also developed for the Nikuradse roughness (k s), and a new equation is proposed which directly relates k s to bedform relative height, aspect ratio and lee side angle.

(Table 2) Radiocarbon dates from lakes on Fildes Peninsula, King George Island

Precise relative sea level (RSL) data are important for inferring regional ice sheet histories, as well as helping to validate numerical models of ice sheet evolution and glacial isostatic adjustment. Here we develop a new RSL curve for Fildes Peninsula, South Shetland Islands (SSIs), a sub-Antarctic archipelago peripheral to the northern Antarctic Peninsula ice sheet, by integrating sedimentary evidence from isolation basins with geomorphological evidence from raised beaches. This combined approach yields not only a Holocene RSL curve, but also the spatial pattern of how RSL change varied across the archipelago. The curve shows a mid-Holocene RSL highstand on Fildes Peninsula at 15.5 m above mean sea level between 8000 and 7000 cal a BP. Subsequently RSL gradually fell as a consequence of isostatic uplift in response to regional deglaciation. We propose that isostatic uplift occurred at a non-steady rate, with a temporary pause in ice retreat ca. 7200 cal a BP, leading to a short-lived RSL rise of ~1 m and forming a second peak to the mid-Holocene highstand. Two independent approaches were taken to constrain the long-term tectonic uplift rate of the SSIs at 0.22-0.48 m/ka, placing the tectonic contribution to the reconstructed RSL highstand between 1.4 and 2.9 m. Finally, we make comparisons to predictions from three global sea level models.

Documentation of sediment core PS1392-1

The Late Quaternary sediment sequence of the continental margin in the eastern Weddell Sea is well suited for palaeoenvironmental reconstructions. Two cores from the upper slope, which contain the sedimentary record of the last 300 ky, have been sedimentologically investigated. Age models are based on lithostratigraphy and are correlated with the stable isotope record. As a result of a detailed analysis of the clay mineral composition, grain size distributions and structures, this sedimentary record provides the first marine evidence that the Antarctic ice sheet extended to the shelf edge during the last glacial. The variations in volume and size of the ice sheet were also simulated in numerical models. Changes in accumulation rate and ice temperature are of some importance, but the model revealed that fluctuations are primarily driven by changes in eustatic sea-level and that the ice edge extended to the shelf edge during the last glacial maximum. This causal relationship implies that the maximum ice extension strongly depends on the magnitude and duration of the sea-level depression during a glacial period. The results of the sedimentological investigations and of the numerical models show that the Antarctic ice sheet follows glacial events in the northern hemisphere by teleconnections of sea level.