Geochemistry, isotopes, trace elements and fluid inclusion geothermometry at DSDP Hole 83-504B

Stockwork-like metal sulfide mineralizations were found at 910-928 m below seafloor (BSF) in the pillow/dike transition zone of Hole 504B. This is the same interval where most physical properties of the 5.9-m.y.-old crust of the Costa Rica Rift change from those characteristic of Layer 2B to those of Layer 2C. The pillow lavas, breccias, and veins of the stockwork-like zone were studied by transmitted and reflected light microscopy, X-ray diffraction, and electron microprobe analysis. Bulk rock oxygen isotopic analyses as well as isolated mineral oxygen and sulfur isotopic analyses and fluid inclusion measurements were carried out. A complex alteration history was reconstructed that includes three generations of fissures, each followed by precipitation of characteristic hydrothermal mineral parageneses: (1) Minor and local deposition of quartz occurred on fissure walls; adjacent wall rocks were silicified, followed by formation of chlorite and minor pyrite I in the veins, whereas albite, sphene, chlorite and chlorite-expandable clay mixtures, actinolite, and pyrite replaced igneous phases in the host rocks. The hydrothermal fluids responsible for this first stage were probably partially reacted seawater, and their temperatures were at least 200-250° C. (2) Fissures filled during the first stage were reopened and new cracks formed. They were filled with quartz, minor chlorite and chlorite-expandable clay mixtures, traces of epidote, common pyrite, sphalerite, chalcopyrite, and minor galena. During the second stage, hydrothermal fluids were relatively evolved metal- and Si-rich solutions whose temperatures ranged from 230 to 340° C. The fluctuating chemical composition and temperature of the solutions produced a complex depositional sequence of sulfides in the veins: chalcopyrite I, ± Fe-rich sphalerite, chalcopyrite II ("disease"), Fe-poor sphalerite, chalcopyrite III, galena, and pyrite II. (3) During the last stage, zeolites and Mg-poor calcite filled up the remaining spaces and newly formed cracks and replaced the host rock plagioclase. Analcite and stilbite were first to form in veins, possibly at temperatures below 200°C; analcite and earlier quartz were replaced by laumontite at 250°C, whereas calcite formation temperature ranged from 135 to 220°C. The last stage hydrothermal fluids were depleted in Mg and enriched in Ca and 18O compared to seawater and contained a mantle carbon component. This complex alteration history paralleling a complex mineral paragenesis can be interpreted as the result of a relatively long-term evolution of a hydrothermal system with superimposed shorter term fluctuations in solution temperature and composition. Hydrothermal activity probably began close to the axis of the Costa Rica Rift with the overall cooling of the system and multiple fracturing stages due to movement of the crust away from the axis and/or cooling of a magmatic heat source.

Chemical and Sr isotopic compositions of pore fluids and marine carbonates from ODP Hole 130-807A

The 87Sr/86Sr ratios and Sr concentrations in sediment and pore fluids are used to evaluate the rates of calcite recrystallization at ODP Site 807A on the Ontong Java Plateau, an 800-meter thick section of carbonate ooze and chalk. A numerical model is used to evaluate the pore fluid chemistry and Sr isotopes in an accumulating section. The deduced calcite recrystallization rate is 2% per million years (%/Myr) near the top of the section and decreases systematically in older parts of the section such that the rate is close to 0.1/age (in years). The deduced recrystallization rates have important implications for the interpretation of Ca and Mg concentration profiles in the pore fluids. The effect of calcite recrystallization on pore fluid chemistry is described by the reaction length, L, which varies by element, and depends on the concentration in pore fluid and solid. When L is small compared to the thickness of the sedimentary section, the pore fluid concentration is controlled by equilibrium or steady-state exchange with the solid phase, except within a distance L of the sediment-water interface. When L is large relative to the thickness of sediment, the pore fluid concentration is mostly controlled by the boundary conditions and diffusion. The values of L for Ca, Sr, and Mg are of order 15, 150, and 1500 meters, respectively. L_Sr is derived from isotopic data and modeling, and allows us to infer the values of L_Ca and L_Mg. The small value for L_Ca indicates that pore fluid Ca concentrations, which gradually increase down section, must be equilibrium values that are maintained by solution-precipitation exchange with calcite and do not reflect Ca sources within or below the sediment column. The pore fluid Ca measurements and measured alkalinity allow us to calculate the in situ pH in the pore fluids, which decreases from 7.6 near the sediment-water interface to 7.1+/-0.1 at 400-800 mbsf. While the calculated pH values are in agreement with some of the values measured during ODP Leg 130, most of the measurements are artifacts. The large value for L_Mg indicates that the pore fluid Mg concentrations at 807A are not controlled by calcite-fluid equilibrium but instead are determined by the changing Mg concentration of seawater during deposition, modified by aqueous diffusion in the pore fluids. We use the pore fluid Mg concentration profile at Site 807A to retrieve a global record for seawater Mg over the past 35 Myr, which shows that seawater Mg has increased rapidly over the past 10 Myr, rather than gradually over the past 60 Myr. This observation suggests that the Cenozoic rise in seawater Mg is controlled by continental weathering inputs rather than by exchange with oceanic crust. The relationship determined between reaction rate and age in silicates and carbonates is strikingly similar, which suggests that reaction affinity is not the primary determinant of silicate dissolution rates in nature.

Properties, sedimentation and accumulation rates, and chemical composition of bottom sediments from the Deryugin Basin, Sea of Okhotsk

Results of a study of contents and accumulation rates of Fe, Mn, and some trace elements in Upper Quaternary sediments of the Deryugin Basin are presented. Maps of average contents and accumulation rates of excessive Fe, Mn, Zn, Ba, Ni, Pb, Cu, and Mo in sediments of the first oxygen isotope stage (OIS) have been plotted. Anomalous contents and accumulation rates are confined to peripheral zones of the Deryugin sedimentary basin and large fracture zones. Different mechanisms of influence of fluid-dynamic processes on rate of hydrogenic and biogenic accumulation of ore elements are assumed.

Strontium, oxygen and hydrogen isotope ratios of basalts from DSDP Hole 504B

DSDP Hole 504B was drilled into 6 Ma crust, about 200 km south of the Costa Rica Rift, Galapagos Spreading Center, penetrating 1.35 km into a section that can be divided into four zones-Zone I: oxic submarine weathering; Zone II: anoxic alteration; Zones III and IV: hydrothermal alteration to greenschist facies. In Zone III there is intense veining of pillow basalts. Zone IV consists of altered sheeted dikes. Isotopic geochemical signatures in relation to the alteration zones are recorded in Hole 504B, as follows: Zone Depth(m) Average87Sr/86Sr Average delta18O (?) Average deltaD (?) I 275-550 0.7032 7.3 -63 II 550-890 0.7029 6.5 -45 III 890-1050 0.7035 5.6 -31 IV 1050-1350 0.7032 5.5 -36 Alteration temperatures are as low as 10°C in Zones I and II based on oxygen isotope fractionation. Strontium isotopic data indicate that a circulation of seawater is much more restricted in Zone II than in Zone I. Fluid inclusion measurements of vein quartz indicate the alteration temperature was mainly 300 +/- 20°C in Zones III and IV, which is consistent with secondary mineral assemblages. The strontium, oxygen, and hydrogen isotopic compositions of hydrothermal fluids which were responsible for the greenschist facies alteration in Zones III and IV are estimated to be 0.7037, 2?, and 3?, respectively. Strontium and oxygen isotope data indicate that completely altered portions of greenstones and vein minerals were in equilibrium with modified seawater under low water/rock ratios (in weight) of about 1.6. This value is close to that of the end-member hydrothermal fluids issuing at 21°N EPR. Basement rocks are not completely hydrothermally altered. About 32% of the greenstones in Zones III and IV have escaped alteration. Thus 1 g of fresh basalt including the 32% unaltered portion are required in order to make 1 g of end-member solution from fresh seawater in water-rock reactions.

Chemical composition of plagioclase from DSDP Legs 45 and 46, Atlantic Ocean

Phyric basalts recovered from DSDP Legs 45 and 46 contain abundant plagioclase phenocrysts which occur as either discrete single grains (megacrysts) or aggregates (glomerocrysts) and which are too abundant and too anorthitic to have crystallized from a liquid with the observed bulk rock composition. Almost all the plagioclase crystals are complexly zoned. In most cases two abrupt and relatively large compositional changes associated with continuous internal morphologic boundaries divide the plagioclase crystals into three parts: core, mantle and rim. The cores exhibit two major types of morphology: tabular, with a euhedral to slightly rounded outline; or a skeletal inner core wrapped by a slightly rounded homogeneous outer core. The mantle region is characterized by a zoning pattern composed of one to several spikes/plateaus superimposed on a gently zoned base line, with one large plateau always at the outside of the mantle, and by, in most cases, a rounded internal morphology. The inner rim is typically oscillatory zoned. The width of the outer rim can be correlated with the position of the individual crystal in the basalt pillow. The presence of a skeletal inner core and the concentration of glass inclusions in low-An zones in the mantle region suggest that the liquid in which these parts of the crystals were growing was undercooled some amount. The resorption features at the outer margins of low-An zones indicate superheating of the liquid with respect to the crystal. It is proposed that the plagioclase cores formed during injection of primitive magma into a previously existing magma chamber, that the mantle formed during mixing of a partially mixed magma and the remaining magma already in the chamber, and that the inner rim formed when the mixed magma was in a sheeted dike system. The large plateau at the outside of the mantle may have formed during the injection of the next batch of primitive magma into the main chamber, which may trigger an eruption. This model is consistent with fluid dynamic calculations and geochemically based magma mixing models, and is suggested to be the major mechanism for generating the disequilibrium conditions in the magma.

In-situ geotechnical measurements with the dynamic penetrometer Nimrod on the continental shelf off the Coromandel Peninsula, New Zealand

We tested the ability of a small dynamic penetrometer, Nimrod, to infer geotechnical properties of sediment mixtures in the inner shelf. The penetrometer is light and easy to operate, and its operation by scuba divers ensures a greater degree of precision than ship-based penetrometer deployments. We have studied selected positions along a sorted bedform (~ 100 m wide) on the continental shelf off the Coromandel Peninsula close to Tairua, North Island of New Zealand, and additionally took sediment samples at the exact positions of penetrometer impact, also by scuba divers. The derived dynamic penetrometer signatures (i) measured deceleration of the probe and estimated quasi-static bearing capacity as a measure of sediment strength, (ii) reflected changes in grain-size distribution ranging from very fine to very coarse sands, and (iii) revealed the uppermost seafloor stratification (top layer 2-6 cm) potentially being an indicator for sediment dynamics. In this manner, the device proved to be suitable for spatially fine-scaled surveys using divers' support and might deliver complementary information about sediment dynamics, in this case sorted-bedform maintenance.

Nimrod (dynamic penetrometer) measurements of a sandbar at Raglan, New Zealand in 2009

Dynamic penetrometer data obtained with the Nimrod penetrometer (MARUM). Data is presented as (i) penetration depth (including for different layers if present), (ii) measured deceleration and (iv) estimated quasi-static bearing capacity including range of uncertainty due to the processing method. Lat/Long coordinates are given.

Permeability measurements on sheared and unsheared samples from different IODP Holes of Leg 316 and from IODP Hole 315-C0001E and ODP Hole 190-1174B

At subduction zones, the permeability of major fault zones influences pore pressure generation, controls fluid flow pathways and rates, and affects fault slip behavior and mechanical strength by mediating effective normal stress. Therefore, there is a need for detailed and systematic permeability measurements of natural materials from fault systems, particularly measurements that allow direct comparison between the permeability of sheared and unsheared samples from the same host rock or sediment. We conducted laboratory experiments to compare the permeability of sheared and uniaxially consolidated (unsheared) marine sediments sampled during IODP Expedition 316 and ODP Leg 190 to the Nankai Trough offshore Japan. These samples were retrieved from: (1) The décollement zone and incoming trench fill offshore Shikoku Island (the Muroto transect); (2) Slope sediments sampled offshore SW Honshu (the Kumano transect) ~ 25 km landward of the trench, including material overriden by a major out-of-sequence thrust fault, termed the "megasplay"; and (3) A region of diffuse thrust faulting near the toe of the accretionary prism along the Kumano transect. Our results show that shearing reduces fault-normal permeability by up to 1 order of magnitude, and this reduction is largest for shallow (< 500 mbsf) samples. Shearing-induced permeability reduction is smaller in samples from greater depth, where pre-existing fabric from compaction and lithification may be better developed. Our results indicate that localized shearing in fault zones should result in heterogeneous permeability in the uppermost few kilometers in accretionary prisms, which favors both the trapping of fluids beneath and within major faults, and the channeling of flow parallel to fault structure. These low permeabilities promote the development of elevated pore fluid pressures during accretion and underthrusting, and will also facilitate dynamic hydrologic processes within shear zones including dilatancy hardening and thermal pressurization.

Temperature measurements, gas chemistry, and gas hydrates of the Dvurechenskii mud volcano

The sediment temperature distribution at mud volcanoes provides insights into their activity and into the occurrence of gas hydrates. If ambient pressure and temperature conditions are close to the limits of the gas hydrate stability field, the sediment temperature distribution not only limits the occurrence of gas hydrates, but is itself influenced by heat production and consumption related to the formation and dissociation of gas hydrates. Located in the Sorokin Trough in the northern Black Sea, the Dvurechenskii mud volcano (DMV) was in the focus of detailed investigations during the M72/2 and M73/3a cruises of the German R/V Meteor and the ROV Quest 4000 m in February and March 2007. A large number of in-situ sediment temperature measurements were conducted from the ROV and with a sensor-equipped gravity corer. Gas hydrates were sampled in pressurized cores using a dynamic autoclave piston corer (DAPC). The thermal structure of the DMV suggests a regime of fluid flow at rates decreasing from the summit towards the edges of the mud volcano, accompanied by intermittent mud expulsion at the summit. Modeled gas hydrate dissociation temperatures reveal that the gas hydrates at the DMV are very close to the stability limits. Changes in heat flow due to variable seepage rates probably do not result in changes in sediment temperature but are compensated by gas hydrate dissociation and formation.

Mean dynamic topography and oceanographic parameters estimated from an inverse model and satellite geodesy, with link to model result in one single NetCDF file (392 MB), including the inverse data error covariance

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography-the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99-110, 2012, doi:10.1016/j.jog.2011.07.0069; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie, 2012, http://nbn-resolving.de/nbn:de:hbz:5n-29199). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.

In-situ Iron isotope measurements and mineral compositions from ODP Hole 206-1256D and IODP Hole 335-U1256D

In-situ Fe isotope measurements have been carried out to estimate the impact of the hydrothermal metamorphic overprint on the Fe isotopic composition of Fe-Ti-oxides and Fe-sulfides of the different lithologies of the drilled rocks from IODP Hole 1256D (eastern equatorial Pacific; 15 Ma crust formed at the East Pacific Rise). Most igneous rocks normally have a very restricted range in their 56Fe/54Fe ratio. In contrast, Fe isotope compositions of hot fluids (> 300 °C) from mid-ocean-ridge spreading centers define a narrow range that is shifted to lower delta 56Fe values by 0.2 per mil - 0.5 per mil as compared to igneous rocks. Therefore, it is expected that mineral phases that contain large amounts of Fe are especially affected by the interaction with a fluid that fractionates Fe isotopes during exsolution/precipitation of those minerals. We have used a femtosecond UV-Laser ablation system to determine mineral 56Fe/54Fe ratios of selected samples with a precision of < 0.1 per mil (2 sigma level) at micrometer-scale. We have found significant variations of the delta 56Fe (IRMM-014) values in the minerals between different samples as well as within samples and mineral grains. The overall observed scale of delta 56Fe (magnetite) in 1256D rocks ranges from - 0.12 to + 0.64 per mil, and of delta 56Fe (ilmenite) from - 0.77 to + 0.01 per mil. Pyrite in the lowermost sheeted dike section is clearly distinguishable from the other investigated lithological units, having positive delta 56Fe values between + 0.29 and + 0.56 per mil, whereas pyrite in the other samples has generally negative delta 56Fe values from - 1.10 to - 0.59 permil. One key observation is that the temperature dependent inter-mineral fractionations of Fe isotopes between magnetite and ilmenite are systematically shifted towards higher values when compared to theoretically expected values, while synthesized, well equilibrated magnetite-ilmenite pairs are compatible with the theoretical predictions. Theoretical considerations including beta-factors of different aqueous Fe-chlorides and Rayleigh-type fractionations in the presence of a hydrous, chlorine-bearing fluid can explain this observation. The disagreement between observed and theoretical equilibrium fractionation, the fact that magnetite, in contrast to ilmenite shows a slight downhole trend in the delta 56Fe values, and the observation of small scale heterogeneities within single mineral grains imply that a general re-equilibration of the magnetite-ilmenite pairs is overprinted by kinetic fractionation effects, caused by the interaction of magnetite/ilmenite with hydrothermal fluids penetrating the upper oceanic crust during cooling, or incomplete re-equilibration at low temperatures. Furthermore, the observation of significant small-scale variations in the 56Fe/54Fe ratios of single minerals in this study highlights the importance of high spatial-resolution-analyses of stable isotope ratios for further investigations.

Pore fluid pressure at DSDP Leg 84 Holes

Evidence of considerable overpressuring of pore fluids in the sediment drilled during Leg 84 was obtained from direct measurement of pressure by two methods. The first involved measurement of back pressure when the annulus of the drill hole became constricted with unremoved drill cuttings or constriction was caused by plastic inflow of the drill hole walls. The second involved measurement of pressure ahead of the bit in conjunction with in situ water samples and heat flow. All measurements indicated abnormally high pore pressure even in slope deposits of the Middle America Trench off Guatemala.