Methane discharge in sediments of the Dvurechenskii mud volcano

During the MARGASCH cruise M52/1 in 2001 with RV Meteor we sampled surface sediments from three stations in the crater of the Dvurechenskii mud volcano (DMV, located in the Sorokin Trough of the Black Sea) and one reference station situated 15 km to the northeast of the DMV. We analysed the pore water for sulphide, methane, alkalinity, sulphate, and chloride concentrations and determined the concentrations of particulate organic carbon, carbonate and sulphur in surface sediments. Rates of anaerobic oxidation of methane (AOM) were determined using a radiotracer (14CH4) incubation method. Numerical transport-reaction models were applied to derive the velocity of upward fluid flow through the quiescently dewatering DMV, to calculate rates of AOM in surface sediments, and to determine methane fluxes into the overlying water column. According to the model, AOM consumes 79% of the average methane flux from depth (8.9 x 10**+ 6 mol a**-1), such that the resulting dissolved methane emission from the volcano into the overlying bottom water can be determined as 1.9 x 10**+ 6 mol a**-1. If it is assumed that all submarine mud volcanoes (SMVs) in the Black Sea are at an activity level like the DMV, the resulting seepage represents less than 0.1% of the total methane flux into this anoxic marginal sea. The new data from the DMV and previously published studies indicate that an average SMV emits about 2.0 x 10**+ 6 mol a**-1 into the ocean via quiescent dewatering. The global flux of dissolved methane from SMVs into the ocean is estimated to fall into the order of 10**+10 mol a**-1. Additional methane fluxes arise during periods of active mud expulsion and gas bubbling occurring episodically at the DMV and other SMVs.

Coral-algae metabolism and diurnal changes in the CO2-carbonate system of bulk sea water

Precise measurements were conducted in continuous flow seawater mesocosms located in full sunlight that compared metabolic response of coral, coral-macroalgae and macroalgae systems over a diurnal cycle. Irradiance controlled net photosynthesis (Pnet), which in turn drove net calcification (Gnet), and altered pH. Pnet exerted the dominant control on [CO3]2- and aragonite saturation state (Omega arag) over the diel cycle. Dark calcification rate decreased after sunset, reaching zero near midnight followed by an increasing rate that peaked at 03:00 h. Changes in Omega arag and pH lagged behind Gnet throughout the daily cycle by two or more hours. The flux rate Pnet was the primary driver of calcification. Daytime coral metabolism rapidly removes dissolved inorganic carbon (DIC) from the bulk seawater and photosynthesis provides the energy that drives Gnet while increasing the bulk water pH. These relationships result in a correlation between Gnet and Omega arag, with Omega arag as the dependent variable. High rates of H+ efflux continued for several hours following mid-day peak Gnet suggesting that corals have difficulty in shedding waste protons as described by the Proton Flux Hypothesis. DIC flux (uptake) followed Pnet and Gnet and dropped off rapidly following peak Pnet and peak Gnet indicating that corals can cope more effectively with the problem of limited DIC supply compared to the problem of eliminating H+. Over a 24 h period the plot of total alkalinity (AT) versus DIC as well as the plot of Gnet versus Omega arag revealed a circular hysteresis pattern over the diel cycle in the coral and coral-algae mesocosms, but not the macroalgae mesocosm. Presence of macroalgae did not change Gnet of the corals, but altered the relationship between Omega arag and Gnet. Predictive models of how future global changes will effect coral growth that are based on oceanic Omega arag must include the influence of future localized Pnet on Gnet and changes in rate of reef carbonate dissolution. The correlation between Omega arag and Gnet over the diel cycle is simply the response of the CO2-carbonate system to increased pH as photosynthesis shifts the equilibria and increases the [CO3]2- relative to the other DIC components of [HCO3]- and [CO2]. Therefore Omega arag closely tracked pH as an effect of changes in Pnet, which also drove changes in Gnet. Measurements of DIC flux and H+ flux are far more useful than concentrations in describing coral metabolism dynamics. Coral reefs are systems that exist in constant disequilibrium with the water column.

(Table 2) Composition of the deapest pore fluids from ODP Holes 168-1030B and 168-1031A

On the eastern flank of the Juan de Fuca Ridge, reaction between upwelling basement fluid and sediment alters hydrothermal fluxes of Ca, SiO2(aq), SO4, PO4, NH4, and alkalinity. We used the Global Implicit Multicomponent Reactive Transport (GIMRT) code to model the processes occurring in the sediment column (diagenesis, sediment burial, fluid advection, and multicomponent diffusion) and to estimate net seafloor fluxes of solutes. Within the sediment section, the reactions controlling the concentrations of the solutes listed above are organic matter degradation via SO4 reduction, dissolution of amorphous silica, reductive dissolution of amorphous Fe(III)-(hydr)oxide, and precipitation of calcite, carbonate fluorapatite, and amorphous Fe(II)-sulfide. Rates of specific discharge estimated from pore-water Mg profiles are 2 to 3 mm/yr. At this site the basement hydrothermal system is a source of NH4, SiO2(aq), and Ca, and a sink of SO4, PO4, and alkalinity. Reaction within the sediment column increases the hydrothermal sources of NH4 and SiO2(aq), increases the hydrothermal sinks of SO4 and PO4, and decreases the hydrothermal source of Ca. Reaction within the sediment column has a spatially variable effect on the hydrothermal flux of alkalinity. Because the model we used was capable of simulating the observed pore-water chemistry by using mechanistic descriptions of the biogeochemical processes occurring in the sediment column, it could be used to examine the physical controls on hydrothermal fluxes of solutes in this setting. Two series of simulations in which we varied fluid flow rate (1 to 100 mm/yr) and sediment thickness (10 to 100 m) predict that given the reactions modeled in this study, the sediment section will contribute most significantly to fluxes of SO4 and NH4 at slow flow rates and intermediate sediment thickness and to fluxes of SiO2(aq) at slow flow rates and large sediment thickness. Reaction within the sediment section could approximately double the hydrothermal sink of PO4 over a range of flow rates and sediment thickness, and could slightly decrease (by

Lithologic description and vertical permeability measurements of ODP Site 180-1109

Vertical permeability testing was conducted on four samples collected from Site 1109, a borehole advanced during Ocean Drilling Program Leg 180. Closed conditions were applied during each test, and the samples were measured using a constant flow approach and permeant solutions that matched the geochemistry of nearby interstitial waters. Vertical permeabilities measured at 34.5 kPa effective stress generally decreased with depth and ranged from 10**-14 m**2 at 212.53 meters below seafloor (mbsf) to 10**-18 m**2 at 698.10 mbsf. The three deepest samples differed in permeability by less than one order of magnitude. Reconsolidation testing on the shallowest sample yielded a minimum permeability of 1.56 x 10**-16 m**2 at 276 kPa effective stress. Subsequent rebound testing yielded a hysteresis-type curve, with the final permeability measuring lower than the initial permeability by nearly 1.5 orders of magnitude. Dilution experiments indicated that use of a permeant solution matching the geochemistry of the interstitial waters may be necessary for accuracy in measurements and mitigation of clay swellage and collapse during testing, but further research is mandated.

Corrected ICESat altimetry data, surface mass balance, and firn elevation change on Antarctic ice shelves

Accurate prediction of global sea-level rise requires that we understand the cause of recent, widespread and intensifying glacier acceleration along Antarctic ice-sheet coastal margins. Floating ice shelves buttress the flow of grounded tributary glaciers and their thickness and extent are particularly susceptible to changes in both climate and ocean forcing. Recent ice-shelf collapse led to retreat and acceleration of several glaciers on the Antarctic Peninsula. However, the extent and magnitude of ice-shelf thickness change, its causes and its link to glacier flow rate are so poorly understood that its influence on the future of the ice sheets cannot yet be predicted. Here we use satellite laser altimetry and modelling of the surface firn layer to reveal for the first time the circum-Antarctic pattern of ice-shelf thinning through increased basal melt. We deduce that this increased melt is the primary driver of Antarctic ice-sheet loss, through a reduction in buttressing of the adjacent ice sheet that has led to accelerated glacier flow. The highest thinning rates (~7 m/a) occur where warm water at depth can access thick ice shelves via submarine troughs crossing the continental shelf. Wind forcing could explain the dominant patterns of both basal melting and the surface melting and collapse of Antarctic ice shelves, through ocean upwelling in the Amundsen and Bellingshausen Seas and atmospheric warming on the Antarctic Peninsula. This implies that climate forcing through changing winds influences Antarctic Ice Sheet mass balance, and hence global sea-level, on annual to decadal timescales.

Shipborne ocean color remote sensing

In this paper, modernized shipborne procedures are presented to collect and process above-water radiometry for remote sensing applications. A setup of five radiometers and a bidirectional camera system, which provides panoramic sea surface and sky images, is proposed for the collection of high-resolution radiometric quantities. Images from the camera system can be used to determine sky state and potential glint, whitecaps, or foam contamination. A peak in the observed remote sensing reflectance RRS spectra between 750-780 nm was typically found in spectra with relatively high surface reflected glint (SRG), which suggests this waveband could be a useful SRG indicator. Simplified steps for computing uncertainties in SRG corrected RRS are proposed and discussed. The potential of utilizing "unweighted multimodel averaging," which is the average of four or more common SRG correction models, is examined to determine the best approximation RRS. This best approximation RRS provides an estimate of RRS based on various SRG correction models established using radiative transfer simulations and field investigations. Applying the average RRS provides a measure of the inherent uncertainties or biases that result from a user subjectively choosing any one SRG correction model. Comparisons between inherent and apparent optical property derived observations were used to assess the robustness of the SRG multimodel averaging ap- proach. Correlations among the standard SRG models were completed to determine the degree of association or similarities between the SRG models. Results suggest that the choice of glint models strongly affects derived RRS values and can also influence the blue to green band ratios used for modeling biogeochemical parameters such as for chlorophyll a. The objective here is to present a uniform and traceable methodology for determining ship- borne RRS measurements and its associated errors due to glint correction and to ensure the direct comparability of these measurements in future investigations. We encourage the ocean color community to publish radiometric field measurements with matching and complete metadata in open access repositories.

Consolidation properties and permeability values of sediments from three ODP Leg 205 sites

Vertical permeability and sediment consolidation measurements were taken on seven whole-round drill cores from Sites 1253 (three samples), 1254 (one sample), and 1255 (three samples) drilled during Ocean Drilling Program Leg 205 in the Middle America Trench off of Costa Rica's Pacific Coast. Consolidation behavior including slopes of elastic rebound and virgin compression curves (Cc) was measured by constant rate of strain tests. Permeabilities were determined from flow-through experiments during stepped-load tests and by using coefficient of consolidation (Cv) values continuously while loading. Consolidation curves and the Casagrande method were used to determine maximum preconsolidation stress. Elastic slopes of consolidation curves ranged from 0.097 to 0.158 in pelagic sediments and 0.0075 to 0.018 in hemipelagic sediments. Cc values ranged from 1.225 to 1.427 for pelagic carbonates and 0.504 to 0.826 for hemipelagic clay-rich sediments. In samples consolidated to an axial stress of ~20 MPa, permeabilities determined by flow-through experiments ranged from a low value of 7.66 x 10**-20 m**2 in hemipelagic sediments to a maximum value of 1.03 x 10**-16 m**2 in pelagic sediments. Permeabilities calculated from Cv values in the hemipelagic sediments ranged from 4.81 x 10**-16 to 7.66 x 10**-20 m**2 for porosities 49.9%-26.1%.

Effective stress, porosity, p-wave velocity and mineral composition of ODP Hole 174A-1073A sediments

Porosity, permeability, and compressional (P-wave) velocity were measured as a function of stress on sediments from Ocean Drilling Program Site 1073, U.S. Mid-Atlantic continental slope. Thin sections, scanning electron microscopy, and X-ray diffraction analyses provided mineralogical characteristics of the samples. Uniaxial strain boundary conditions were imposed on the samples during consolidation tests with the maximum effective axial stress reaching 13 MPa. The maximum effective radial stress necessary to maintain uniaxial strain was 7.6 MPa. Over an effective axial stress interval of 0 to 5.2 MPa, Sample 174A-1073A-26X-2, 82-89 cm (226.65 meters below seafloor [mbsf]), exhibited the largest decrease in porosity (51% to 41%), whereas Sample 71X-1, 2-8 cm (644.70 mbsf), exhibited the smallest decrease in porosity (48% to 45%). All samples showed negligible porosity increases during unloading. The permeability (on the order of 1 x 10-17 m**2) of Sample 174A-1073A-71X-1, 2-8 cm, was twice that measured on Sample 8H-1, 23-26 cm (63.75 mbsf), even though the former was considerably deeper and older. The differences in porosity-stress behavior and permeability between shallow and deep samples is related to lithologic, mineralogic, and diagenetic differences between the sediments above and below the Pliocene-Pleistocene to Miocene unconformity. P-wave velocity for Samples 174A-1073A-41X-5, 97-103 cm (372.35 mbsf), and 71X-1, 2-8 cm, increased with decreasing porosity, but did not change significantly during unloading.