An environmental magnetic study, sedimentological and numerical methods around Tauranga Harbour, New Zealand

Magnetic iron minerals are widespread and indicative sediment constituents in estuarine, coastal and shelf systems. We combine environmental magnetic, sedimentological and numerical methods to identify magnetite-enriched placer-like zones in a complex coastal system and delineate their formation mechanisms. Magnetic susceptibility and remanence measurements on 245 surficial sediment samples collected in and around Tauranga Harbour, the largest barrier-enclosed tidal estuary of New Zealand, reveal several discrete enrichment zones controlled by local hydrodynamic conditions. Active magnetite enrichment takes place in tidal channels, which feed into two coast-parallel nearshore magnetite-enriched belts centered at water depths of 6-10 m and 10-20 m. A close correlation between magnetite content and magnetic grain size was found, where higher susceptibility values are associated within coarser magnetic crystal sizes. Two key mechanisms for magnetite enrichment are identified. First, tide-induced residual currents primarily enable magnetite enrichment within the estuarine channel network. A coast-parallel, fine sand magnetite enrichment belt in water depths of less than 10 m along the barrier island has a strong decrease in magnetite content away from the southern tidal inlet and is apparently related to active coast-parallel transport combined with mobilizing surf zone processes. A second, less pronounced, but more uniform magnetite enrichment belt at 10-20 m water depth is composed of non-mobile, medium-coarse-grained relict sands, which have been reworked during post-glacial sea level transgression. We demonstrate the potential of magnetic methods to reveal and differentiate coastal magnetite enrichment patterns and investigate their formative mechanisms.

Major, trace element and XRF data for sediments and fish tooth samples, from a shallow marine setting with varying redox conditions (Site U1356, Wilkes Land Margin, E.Antarctica)

Fossil fish teeth from pelagic open ocean settings are considered a robust archive for preserving the neodymium (Nd) isotopic composition of ancient seawater. However, using fossil fish teeth as an archive to reconstruct seawater Nd isotopic compositions in different sedimentary redox environments and in terrigenous-dominated, shallow marine settings is less proven. To address these uncertainties, fish tooth and sediment samples from a middle Eocene section deposited proximal to the East Antarctic margin at Integrated Ocean Drilling Program Site U1356 were analyzed for major and trace element geochemistry, and Nd isotopes. Major and trace element analyses of the sediments reveal changing redox conditions throughout deposition in a shallow marine environment. However, variations in the Nd isotopic composition and rare earth element (REE) patterns of the associated fish teeth do not correspond to redox changes in the sediments. REE patterns in fish teeth at Site U1356 carry a typical mid-REE-enriched signature. However, a consistently positive Ce anomaly marks a deviation from a pure authigenic origin of REEs to the fish tooth. Neodymium isotopic compositions of cleaned and uncleaned fish teeth fall between modern seawater and local sediments and hence could be authigenic in nature, but could also be influenced by sedimentary fluxes. We conclude that the fossil fish tooth Nd isotope proxy is not sensitive to moderate changes in pore water oxygenation. However, combined studies on sediments, pore waters, fish teeth and seawater are needed to fully understand processes driving the reconstructed signature from shallow marine sections in proximity to continental sources. This article is protected by copyright. All rights reserved.

(Appendix) Petrographic information, trace element and Sr isotope ratios of Pito Deep samples

Fluid flow through the axial hydrothermal system at fast spreading ridges is investigated using the Sr-isotopic composition of upper crustal samples recovered from a tectonic window at Pito Deep (NE Easter microplate). Samples from the sheeted dike complex collected away from macroscopic evidence of channelized fluid flow, such as faults and centimeter-scale hydrothermal veins, show a range of 87Sr/86Sr from 0.7025 to 0.7030 averaging 0.70276 relative to a protolith with 87Sr/86Sr of ~0.7024. There is no systematic variation in 87Sr/86Sr with depth in the sheeted dike complex. Comparison of these new data with the two other localities that similar data sets exist for (ODP Hole 504B and the Hess Deep tectonic window) reveals that the extent of Sr-isotope exchange is similar in all of these locations. Models that assume that fluid-rock reaction occurs during one-dimensional (recharge) flow lead to significant decreases in the predicted extent of isotopic modification of the rock with depth in the crust. These model results show systematic misfits when compared with the data that can only be avoided if the fluid flow is assumed to be focused in isolated channels with very slow fluid-rock exchange. In this scenario the fluid at the base of the crust is little modified in 87Sr/86Sr from seawater and thus unlike vent fluids. Additionally, this model predicts that some rocks should show no change from the fresh-rock 87Sr/86Sr, but this is not observed. Alternatively, models in which fluid-rock reaction occurs during upflow (discharge) as well as downflow, or in which fluids are recirculated within the hydrothermal system, can reproduce the observed lack of variation in 87Sr/86Sr with depth in the crust. Minimum time-integrated fluid fluxes, calculated from mass balance, are between 1.5 and 2.6 * 10**6 kg/m**2 for all areas studied to date. However, new evidence from both the rocks and a compilation of vent fluid compositions demonstrates that some Sr is leached from the crust. Because this leaching lowers the fluid 87Sr/86Sr without changing the rock 87Sr/86Sr, these mass balance models must underestimate the time-integrated fluid flux. Additionally, these values do not account for fluid flow that is channelized within the crust.

Nd and Sr isotope compositions, REE concentrations and Al/Ca ratios of different phases of surface sediments in the South Pacific

The radiogenic isotope composition of neodymium (Nd) and strontium (Sr) are useful tools to investigate present and past oceanic circulation or input of terrigenous material. We present Nd and Sr isotope compositions extracted from different sedimentary phases, including early diagenetic Fe-Mn coatings, "unclean" foraminiferal shells, fossil fish teeth, and detritus of marine surface sediments (core-tops) covering the entire midlatitude South Pacific. Comparison of detrital Nd isotope compositions to deep water values from the same locations suggests that "boundary exchange" has little influence on the Nd isotope composition of western South Pacific seawater. Concentrations of Rare Earth Elements (REE) and Al/Ca ratios of "unclean" planktonic foraminifera suggest that this phase is a reliable recorder of seawater Nd isotope composition. The signatures obtained from fish teeth and "nondecarbonated" leachates of bulk sediment Fe-Mn oxyhydroxide coatings also agree with "unclean" foraminifera. Direct comparison of Nd isotope compositions extracted using these methods with seawater Nd isotope compositions is complicated by the low accumulation rates yielding radiocarbon ages of up to 24 kyr, thus mixing the signal of different ocean circulation modes. This suggests that different past seawater Nd isotope compositions have been integrated in authigenic sediments from regions with low sedimentation rates. Combined detrital Nd and Sr isotope signatures indicate a dominant role of the Westerly winds transporting lithogenic material from South New Zealand and Southeastern Australia to the open South Pacific. The proportion of this material decreases toward the east, where supply from the Andes increases and contributions from Antarctica cannot be ruled out.

Barium intensities and tie points between holes of Leg 208

The Paleocene - Eocene thermal maximum (PETM) is one of the best known examples of a transient climate perturbation, associated with a brief, but intense, interval of global warming and a massive perturbation of the global carbon cycle from injection of isotopically light carbon into the ocean-atmosphere system. One key to quantifying the mass of carbon released, identifying the source(s), and understanding the ultimate fate of this carbon is to develop high-resolution age models. Two independent strategies have been employed, cycle stratigraphy and analysis of extraterrestrial Helium (HeET), both of which were first tested on Ocean Drilling Program (ODP) Site 690. Both methods are in agreement for the onset of the PETM and initial recovery, or the clay layer ("main body"), but seem to differ in the final recovery phase of the event above the clay layer, where the carbonate contents rise and carbon isotope values return toward background values. Here we present a state-of-the-art age model for the PETM derived from a new orbital chronology developed with cycle stratigraphic records from sites drilled during ODP Leg 208 (Walvis Ridge, Southeastern Atlantic) integrated with published records from Site 690 (Weddell Sea, Southern Ocean, ODP Leg 113). During Leg 208, five Paleocene - Eocene (P-E) boundary sections (Sites 1262 to 1267) were recovered in multiple holes over a depth transect of more than 2200 m at the Walvis Ridge yielding the first stratigraphically complete P-E deep-sea sequence with moderate to relatively high sedimentation rates (1 to 3 cm/kyr). A detailed chronology was developed with non-destructive X-ray fluorescence (XRF) core scanning records on the scale of precession cycles, with a total duration of the PETM now estimated to be ~ 170 kyr. The revised cycle stratigraphic record confirms original estimates for the duration of the onset and initial recovery, but suggests a new duration for the final recovery that is intermediate to the previous estimates by cycle stratigraphy and HeET.

K-Ar dating of altered deep-sea igneous rocks from DSDP Legs 2, 34, and 35

In an attempt to establish criteria for obtaining reliable K-Ar dates, conventional K-Ar studies of several Deep Sea Drilling Project sites were undertaken. K-Ar dates of these rocks may be subject to inaccuracies as the result of sea-water alteration. Inaccuracies may also result from the presence of excess radiogenic 40Ar trapped in rapidly cooled rocks at the time of their formation. The results obtained for DSDP Leg 34 basalts indicate that lowering of K-Ar dates, which is related to potassium addition by weathering, is a major cause of uncertainty in obtaining reliable K-Ar dates for deep-sea rocks. It could not be determined if the potassium addition to the basalts occurred at the time of formation, t_o, or continuously from t_o to the present. Calculations show that sediment cover is not a significant barrier to the diffusion of potassium into the basalt. 40Ar loss contributes, at least in part, to the lowering of the K-Ar date in rocks that have added potassium. The meaning of the K-Ar results obtained for DSDP Legs 35 and 2 basalts could not be unambiguously established. Because of the problems involved, caution must be used in interpreting the meaning of conventional K-Ar dates for deep-sea rocks.

Automated Tools for Layer Recognition (BMPix) and Counting (PEAK)

We present tools for rapid and quantitative detection of sediment lamination. The BMPix tool extracts color and gray-scale curves from images at pixel resolution. The PEAK tool uses the gray-scale curve and performs, for the first time, fully automated counting of laminae based on three methods. The maximum count algorithm counts every bright peak of a couplet of two laminae (annual resolution) in a smoothed curve. The zero-crossing algorithm counts every positive and negative halfway-passage of the curve through a wide moving average, separating the record into bright and dark intervals (seasonal resolution). The same is true for the frequency truncation method, which uses Fourier transformation to decompose the curve into its frequency components before counting positive and negative passages. We applied the new methods successfully to tree rings, to well-dated and already manually counted marine varves from Saanich Inlet, and to marine laminae from the Antarctic continental margin. In combination with AMS14C dating, we found convincing evidence that laminations in Weddell Sea sites represent varves, deposited continuously over several millennia during the last glacial maximum. The new tools offer several advantages over previous methods. The counting procedures are based on a moving average generated from gray-scale curves instead of manual counting. Hence, results are highly objective and rely on reproducible mathematical criteria. Also, the PEAK tool measures the thickness of each year or season. Since all information required is displayed graphically, interactive optimization of the counting algorithms can be achieved quickly and conveniently.

(Figures S2a-c) Phase diagrams of the upper crust, the lower crust and the upper mantle based on averaged compositions of the oceanic crust and mantle

The age of the subducting Nazca Plate off Chile increases northwards from 0 Ma at the Chile Triple Junction (46°S) to 37 Ma at the latitude of Valparaíso (32°S). Age-related variations in the thermal state of the subducting plate impact on (a) the water influx to the subduction zone, as well as on (b) the volumes of water that are released under the continental forearc or, alternatively, carried beyond the arc. Southern Central Chile is an ideal setting to study this effect, because other factors for the subduction zone water budget appear constant. We determine the water influx by calculating the crustal water uptake and by modeling the upper mantle serpentinization at the outer rise of the Chile Trench. The water release under forearc and arc is determined by coupling FEM thermal models of the subducting plate with stability fields of water-releasing mineral reactions for upper and lower crust and hydrated mantle. Results show that both the influx of water stored in, and the outflux of water released from upper crust, lower crust and mantle vary drastically over segment boundaries. In particular, the oldest and coldest segments carry roughly twice as much water into the subduction zone as the youngest and hottest segments, but their release flux to the forearc is only about one fourth of the latter. This high variability over a subduction zone of < 1500 km length shows that it is insufficient to consider subduction zones as uniform entities in global estimates of subduction zone fluxes. This article is protected by copyright. All rights reserved.