Pore water methane characteristics and sulphate reduction rates of ODP Leg 164 sediments

Anaerobic methane oxidation (AMO) was characterized in sediment cores from the Blake Ridge collected during Ocean Drilling Program (ODP) Leg 164. Three independent lines of evidence support the occurrence and scale of AMO at Sites 994 and 995. First, concentration depth profiles of methane from Hole 995B exhibit a region of upward concavity suggestive of methane consumption. Diagenetic modeling of the concentration profile indicates a 1.85-m-thick zone of AMO centered at 21.22 mbsf, with a peak rate of 12.4 nM/d. Second, subsurface maxima in tracer-based sulfate reduction rates from Holes 994B and 995B were observed at depths that coincide with the model-predicted AMO zone. The subsurface zone of sulfate reduction was 2 m thick and had a depth integrated rate that compared favorably to that of AMO (1.3 vs. 1.1 nmol/cm**2/d, respectively). These features suggest close coupling of AMO and sulfate reduction in the Blake Ridge sediments. Third, measured d13CH4 values are lightest at the point of peak model-predicted methane oxidation and become increasingly 13C-enriched with decreasing sediment depth, consistent with kinetic isotope fractionation during bacterially mediated methane oxidation. The isotopic data predict a somewhat (60 cm) shallower maximum depth of methane oxidation than do the model and sulfate reduction data.

Chemical and isotopic compositions of pore fluids from ODP Sites 169-1034 and 169-1037

We have analysed the concentrations of Li, K, Rb, Cs, and B, and the isotopic ratios of Li and B of a suite of pore fluids recovered from ODP Sites 1037 (Leg 169; Escanaba Trough) and 1034 (Leg 169S; Saanich Inlet). In addition, we have analysed dissolved K, Rb, and Cs concentrations for estuarine mixing of the Ganges-Brahmaputra river system. Together, these data sets have been used to assess the role of sediments in the marine geochemical cycles of the alkali elements and boron. Uptake onto clay minerals during estuarine mixing removes 20-30% of the riverine input of dissolved Cs and Rb to the oceans. Prior to this study, the only other recognised sink of Rb and Cs was uptake during low-temperature alteration of the oceanic crust. Even with this additional sink there is an excess of inputs over outputs in their modern oceanic mass balance. Pore fluid data show that Li and Rb are transferred into marine sediments during early diagenesis. However, modeling of the Li isotope systematics of the pore fluids from Site 1037 shows that seawater Li taken up during marine sedimentation can be readily returned to solution in the presence of less hydrated cations, such as NH4+. This process also appears to result in high concentrations of pore fluid Cs (relative to local seawater) due to expulsion of adsorbed Cs from cation exchange sites. Flux calculations based on pore fluid data for a series of ODP sites indicate that early diagenesis of clay sediments removes around 8% of the modern riverine input of dissolved Li. Although NH4+-rich fluids do result in a flux of Cs to the oceans, on the global scale this input only augments the modern riverine Cs flux by ~3%. Nevertheless, this may have implications for the fate of radioactive Cs in the natural environment and waste repositories.

A suite of global, cross-scale topographic variables for environmental and biodiversity modeling, links to files in GeoTIFF format

Topographic variation, the spatial variation in elevation and terrain features, underpins a myriad of patterns and processes in geography and ecology and is key to understanding the variation of life on the planet. The characterization of this variation is scale-dependent, i.e. it varies with the distance over which features are assessed and with the spatial grain (grid cell resolution) of analysis. A fully standardized and global multivariate product of different terrain features has the potential to support many large-scale basic research and analytical applications, however to date, such technique is unavailable. Here we used the digital elevation model products of global 250 m GMTED and near-global 90 m SRTM to derive a suite of topographic variables: elevation, slope, aspect, eastness, northness, roughness, terrain roughness index, topographic position index, vector ruggedness measure, profile and tangential curvature, and 10 geomorphological landform classes. We aggregated each variable to 1, 5, 10, 50 and 100 km spatial grains using several aggregation approaches (median, average, minimum, maximum, standard deviation, percent cover, count, majority, Shannon Index, entropy, uniformity). While a global cross-correlation underlines the high similarity of many variables, a more detailed view in four mountain regions reveals local differences, as well as scale variations in the aggregated variables at different spatial grains. All newly-developed variables are available for download at http://www.earthenv.org and can serve as a basis for standardized hydrological, environmental and biodiversity modeling at a global extent.

Strontium isotopes in pore waters of east equatorial Pacific sediments

Pore-water samples from the equatorial sedimentary bulge area show reversals in depth profiles of 87Sr/86Sr ratios at the sediment/basement interface. Results of this work support inferences made from previous pore-water data (from DSDP drilling in the area) that large-scale horizontal advection of seawater has occurred through the basement underlying the thick sedimentary sequence in this region. The area of apparent advection includes the eastern part of the equatorial high-productivity zone and part of the Guatemala Basin. We attempted to find links between the observed near-basement reversals in pore-water chemistry and sedimentary thickness, age, and topography of the area. Most of the sites that show horizontal advection have disturbed basement topography or outcrops within 10 to 20 km, suggesting that the cooling effects of outcrops may extend for at least 20 km horizontally. Heat-flow data from the area were compared to determine whether sites showing near-bottom chemistry reversals were consistent with areas of low conductive heat flow. This was generally true for the area of the sedimentary bulge and Guatemala Basin. Not enough pore-water data from the Nazca Plate were available to establish any reliable systematics. Because the high-productivity area is well-sealed from hydrothermal circulation, the missing heat must be lost by horizontal advective heat transport. From profiles of strontium isotopes and other elements that show departure from seawater values with increasing depth in the sediments, but return to seawater values near the basement, it appears that water flows relatively freely through much of the oceanic crust, even when sealed by considerable sedimentary cover.

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.

Pore water isotopic ratios of CO2 and SO4, and sediment CH4 content of ODP Leg 164 sites

A unique set of geochemical pore-water data, characterizing the sulfate reduction and uppermost methanogenic zones, has been collected at the Blake Ridge (offshore southeastern North America) from Ocean Drilling Program (ODP) Leg 164 cores and piston cores. The d13C values of dissolved CO2 (sum CO2) are as 13C-depleted as -37.7 per mil PDB (Site 995) at the sulfate-methane interface, reflecting a substantial contribution of isotopically light carbon from methane. Although the geochemical system is complex and difficult to fully quantify, we use two methods to constrain and illustrate the intensity of anaerobic methane oxidation in Blake Ridge sediments. An estimate using a two-component mixing model suggests that ~24% of the carbon residing in the sum CO2 pool is derived from biogenic methane. Independent diagenetic modeling of a methane concentration profile (Site 995) indicates that peak methane oxidation rates approach 0.005 µmol/cm**3/yr, and that anaerobic methane oxidation is responsible for consuming ~35% of the total sulfate flux into the sediments. Thus, anaerobic methane oxidation is a significant biogeochemical sink for sulfate, and must affect interstitial sulfate concentrations and sulfate gradients. Such high proportions of sulfate depletion because of anaerobic methane oxidation are largely undocumented in continental rise sediments with overlying oxic bottom waters. We infer that the additional amount of sulfate depleted through anaerobic methane oxidation, fueled by methane flux from below, causes steeper sulfate gradients above methane-rich sediments. Similar pore water chemistries should occur at other methane-rich, continental-rise settings associated with gas hydrates.

Depth profiles of pore-water and solid-phase constituents, respiration rates and grain size distribution in sediments of the Clarion-Clipperton Fracture Zone

Manganese nodules of the Clarion-Clipperton Fracture Zone (CCFZ) in the NE Pacific Ocean are highly enriched in Ni, Cu, Co, Mo and rare-earth elements, and thus may be the subject of future mining operations. Elucidating the depositional and biogeochemical processes that contribute to nodule formation, as well as the respective redox environment in both, water column and sediment, supports our ability to locate future nodule deposits and evaluates the potential ecological and environmental effects of future deep-sea mining. For these purposes we evaluated the local hydrodynamics and pore-water geochemistry with respect to the nodule coverage at four sites in the eastern CCFZ. Furthermore, we carried out selective leaching experiments at these sites in order to assess the potential mobility of Mn in the solid phase, and compared them with the spatial variations in sedimentation rates. We found that the oxygen penetration depth is 180 - 300 cm at all four sites, while reduction of Mn and NO3- is only significant below the oxygen penetration depth at sites with small or no nodules on the sediment surface. At the site without nodules, potential microbial respiration rates, determined by incubation experiments using 14C-labelled acetate, are slightly higher than at sites with nodules. Leaching experiments showed that surface sediments covered with big or medium-sized nodules are enriched in mobilizable Mn. Our deep oxygen measurements and pore-water data suggest that hydrogenetic and oxic-diagenetic processes control the present-day nodule growth at these sites, since free manganese from deeper sediments is unable to reach the sediment surface. We propose that the observed strong lateral contrasts in nodule size and abundance are sensitive to sedimentation rates, which in turn, are controlled by small-scale variations in seafloor topography and bottom-water current intensity.

Pore water chemistry of sediment core GeoB15101-7 and water chemistry of CTD casts GeoB15101-1 and GeoB15101-6

Submarine brine lakes feature sharp and persistent concentration gradients between seawater and brine, though these should be smoothed out by free diffusion in open ocean settings. The anoxic Urania basin of the Eastern Mediterranean contains an ultra sulfidic, hypersaline brine of Messinian origin above a thick layer of suspended sediments. With a dual modeling approach we reconstruct its contemporary stratification by geochemical solute transport fundamentals, and show that thermal convection is required to maintain mixing in the brine and mud layer. The origin of the Urania basin stratification was dated to 1650 years before present, which may be linked to a major earthquake in the region. The persistence of the chemoclines may be key to the development of diverse and specialized microbial communities. Ongoing thermal convection in the fluid mud layer may have important, yet unresolved consequences for sedimentological and geochemical processes, also in similar environments.

Geochemical data from pore water and solids from cores GS08-153-91GC12 and GS08-153-53-GC6, Arctic Mid-Ocean Ridge

Microbial communities and their associated metabolic activity in marine sediments have a profound impact on global biogeochemical cycles. Their composition and structure are attributed to geochemical and physical factors, but finding direct correlations has remained a challenge. Here we show a significant statistical relationship between variation in geochemical composition and prokaryotic community structure within deep-sea sediments. We obtained comprehensive geochemical data from two gravity cores near the hydrothermal vent field Loki's Castle at the Arctic Mid-Ocean Ridge, in the Norwegian-Greenland Sea. Geochemical properties in the rift valley sediments exhibited strong centimeter-scale stratigraphic variability. Microbial populations were profiled by pyrosequencing from 15 sediment horizons (59,364 16S rRNA gene tags), quantitatively assessed by qPCR, and phylogenetically analyzed. Although the same taxa were generally present in all samples, their relative abundances varied substantially among horizons and fluctuated between Bacteria- and Archaea-dominated communities. By independently summarizing covariance structures of the relative abundance data and geochemical data, using principal components analysis, we found a significant correlation between changes in geochemical composition and changes in community structure. Differences in organic carbon and mineralogy shaped the relative abundance of microbial taxa. We used correlations to build hypotheses about energy metabolisms, particularly of the Deep Sea Archaeal Group, specific Deltaproteobacteria, and sediment lineages of potentially anaerobic Marine Group I Archaea. We demonstrate that total prokaryotic community structure can be directly correlated to geochemistry within these sediments, thus enhancing our understanding of biogeochemical cycling and our ability to predict metabolisms of uncultured microbes in deep-sea sediments.

Chemical composition of sediment core GeoB12309-5 from the northern Arabian Sea, its pore water profile and age depth relationship as well as CTD data obtained during cruise M74/3

The Arabian Sea off the Pakistan continental margin is characterized by one of the world's largest oxygen minimum zones (OMZ). The lithology and geochemistry of a 5.3 m long gravity core retrieved from the lower boundary of the modern OMZ (956 m water depth) were used to identify late Holocene changes in oceanographic conditions and the vertical extent of the OMZ. While the lower part of the core (535 - 465 cm, 5.04 - 4.45 cal kyr BP, Unit 3) is strongly bioturbated indicating oxic bottom water conditions, the upper part of the core (284 - 0 cm, 2.87 cal kyr BP to present, Unit 1) shows distinct and well-preserved lamination, suggesting anoxic bottom waters. The transitional interval from 465 to 284 cm (4.45 - 2.87 cal kyr BP, Unit 2) contains relicts of lamination which are in part intensely bioturbated. These fluctuations in bioturbation intensity suggest repetitive changes between anoxic and oxic/suboxic bottom-water conditions between 4.45 - 2.87 cal kyr BP. Barium excess (Baex) and total organic carbon (TOC) contents do not explain whether the increased TOC contents found in Unit 1 are the result of better preservation due to low BWO concentrations or if the decreased BWO concentration is a result of increased productivity. Changes in salinity and temperature of the outflowing water from the Red Sea during the Holocene influenced the water column stratification and probably affected the depth of the lower boundary of the OMZ in the northern Arabian Sea. Even if we cannot prove certain scenarios, we propose that the observed downward shift of the lower boundary of the OMZ was also impacted by a weakened Somali Current and a reduced transport of oxygen-rich Indian Central Water into the Arabian Sea, both as a response to decreased summer insolation and the continuous southward shift of the Intertropical Convergence Zone during the late Holocene.

Uranium contents and isotope ratios of uranium and strontium in sediments, leachates and pore fluids from ODP Hole 162-984A

Bulk dissolution rates for sediment from ODP Site 984A in the North Atlantic are determined using the 234U/238U activity ratios of pore water, bulk sediment, and leachates. Site 984A is one of only several sites where closely spaced pore water samples were obtained from the upper 60 meters of the core; the sedimentation rate is high (11-15 cm/ka), hence the sediments in the upper 60 meters are less than 500 ka old. The sediment is clayey silt and composed mostly of detritus derived from Iceland with a significant component of biogenic carbonate (up to 30%). The pore water 234U/238U activity ratios are higher than seawater values, in the range of 1.2 to 1.6, while the bulk sediment 234U/238U activity ratios are close to 1.0. The 234U/238U of the pore water reflects a balance between the mineral dissolution rate and the supply rate of excess 234U to the pore fluid by a-recoil injection of 234Th. The fraction of 238U decays that result in a-recoil injection of 234U to pore fluid is estimated to be 0.10 to 0.20 based on the 234U/238U of insoluble residue fractions. The calculated bulk dissolution rates, in units of g/g/yr are in the range of 0.0000004 to 0.000002 1/yr. There is significant down-hole variability in pore water 234U/238U activity ratios (and hence dissolution rates) on a scale of ca. 10 m. The inferred bulk dissolution rate constants are 100 to 1000 times slower than laboratory-determined rates, 100 times faster than rates inferred for older sediments based on Sr isotopes, and similar to weathering rates determined for terrestrial soils of similar age. The results of this study suggest that U isotopes can be used to measure in situ dissolution rates in fine-grained clastic materials. The rate estimates for sediments from ODP Site 984 confirm the strong dependence of reactivity on the age of the solid material: the bulk dissolution rate (R_d) of soils and deep-sea sediments can be approximately described by the expression R_d ~ 0.1 1/age for ages spanning 1000 to 500,000,000 yr. The age of the material, which encompasses the grain size, surface area, and other chemical factors that contribute to the rate of dissolution, appears to be a much stronger determinant of dissolution rate than any single physical or chemical property of the system.