The alkali element and boron geochemistry of ODP Site 169-1038 in the Escanaba Trough, East Pacific

A suite of conjugate pore fluid and sediment samples were collected during Leg 169 of the ODP from within the clastic sedimentary sequences which host massive sulphides at Central Hill, Escanaba Trough (ODP Site 1038). We report the alkali element and boron, and Li and B isotope data for these samples. Relative to a reference site (Site 1037) located outside the zone of high heat flow, pore fluids from Site 1038 show a wide variation in Cl (300-800 mM), and have far higher concentrations of Li (up to 6.2 mM), B (up to 9.7 mM), Cs (up to 5.0 mM), and Rb (up to 97 mM). We show that the pore fluids are derived from hydrothermal circulation that has extended into the basement oceanic crust, with input of the alkali elements and B as the rising hydrothermal fluids interact geochemically with the overlying clastic sediments. There is, however, no marked depletion of these elements in the conjugate sediments, suggesting that there has been advective transport of fluids away from the primary hydrothermal reaction site. This is supported by modelling of the Li and B isotope systematics of the pore fluids, which shows that they record extensive formation of secondary minerals during cooling of the fluids from ~350 to ~20ºC. Precipitation of metal-rich sulphides would have occurred prior to the formation of these minerals, thus, the pore fluid Li and B isotope data can place important constraints on the locus of sulphide deposition beneath the seafloor at Escanaba.

Physical oceanography and current meter measurements on the Gulf of Lion continental margin between June 1993 and July 1996

Nine hydrographic cruises were performed on the Gulf of Lion continental margin between June 1993 and July 1996. These observations are analysed to quantify the fluxes of particulate matter and organic carbon transported along the slope by the Northern Current and to characterise their seasonal variability. Concentration of particulate matter and organic carbon are derived from light-transmission data and water sample analyses. The circulation is estimated from the geostrophic current field. The uncertainty on the transport estimate, related to the error on the prediction of particle concentrations from light-transmission data and the error on velocities, is assessed. The particulate matter inflow entering the Gulf of Lion off Marseille is comparable to the Rhône River input and varies seasonally with a maximum transport between autumn and spring. These modifications result from variations of the water flux rather than variations of the particulate matter concentration. Residual transports of particulate matter and organic carbon across the entire Gulf of Lion are calculated for two cruises enclosing the domain that were performed in February 1995 and July 1996. The particulate matter budgets indicate a larger export from the shelf to deep ocean in February 1995 (110 ± 20 kg/s) than in July 1996 (25 ± 18 kg/s). Likewise, the mean particulate organic carbon export is 12.8 ± 0.5 kg/s in February 1995 and 0.8 ± 0.2 kg/s in July 1996. This winter increase is due to larger allochthonous and autochthonous inputs and also to enhanced shelf-slope exchange processes, in particular the cascading of cold water from the shelf. The export of particulate matter by the horizontal currents is moreover two orders of magnitude larger than the vertical particulate fluxes measured at the same time with sediment traps on the continental slope.

Sulphate concentrations and sulphur isotopic composition in pore fluids from ODP Leg 168 sites

A numerical model of sulfate reduction and isotopic fractionation has been applied to pore fluid SO4**2- and d34S data from four sites drilled during Ocean Drilling Program (ODP) Leg 168 in the Cascadia Basin at 48°N, where basement temperatures reach up to 62°C. There is a source of sulfate both at the top and the bottom of the sediment column due to the presence of basement fluid flow, which promotes bacterial sulfate reduction below the sulfate minimum zone at elevated temperatures. Pore fluid d34S data show the highest values (135 per mil) yet found in the marine environment. The bacterial sulfur isotopic fractionation factor, a, is severely underestimated if the pore fluids of anoxic marine sediments are assumed to be closed systems and Rayleigh fractionation plots yield erroneous values for a by as much as 15 per mil in diffusive and advective pore fluid regimes. Model results are consistent with a = 1.077+/-0.007 with no temperature effect over the range 1.8 to 62°C and no effect of sulfate reduction rate over the range 2 to 10 pmol/ccm/day. The reason for this large isotopic fractionation is unknown, but one difference with previous studies is the very low sulfate reduction rates recorded, about two orders of magnitude lower than literature values that are in the range of µmol/ccm/day to tens of nmol/ccm/day. In general, the greatest 34S depletions are associated with the lowest sulfate reduction rates and vice versa, and it is possible that such extreme fractionation is a characteristic of open systems with low sulfate reduction rates.

Permeable coral reef sediment dissolution driven by elevated pCO2 and pore water advection

Ocean acidification (OA) is expected to drive the transition of coral reef ecosystems from net calcium carbonate (CaCO3) precipitating to net dissolving within the next century. Although permeable sediments represent the largest reservoir of CaCO3 in coral reefs, the dissolution of shallow CaCO3 sands under future pCO2 levels has not been measured under natural conditions. In situ, advective chamber incubations under elevated pCO2 (~800 µatm) shifted the sediments from net precipitating to net dissolving. Pore water advection more than doubled dissolution rates (1.10 g CaCO3/m**2/day) when compared to diffusive conditions (0.42 g CaCO3/m**2 /day). Sediment dissolution could reduce net ecosystem calcification rates of the Heron Island lagoon by 8% within the next century, which is equivalent to a 25% reduction in the global average calcification rate of coral lagoons. The dissolution of CaCO3 sediments needs to be taken into account in order to address how OA will impact the net accretion of coral reefs under future predicted increases in CO2.

(Table 1) Isotopic and geochemical parameters of gas phase from underground fluids of the Baikal rift and its framing

Results of studying isotopic composition of helium in underground fluids of the Baikal-Mongolian region during the last quarter of XX century are summarized. Determinations of 3He/4He ratio in 139 samples of gas phase from fluids, collected at 104 points of the Baikal rift zone and adjacent structures are given. 3He/4He values lie within the range from 1x10**-8 (typical for crustal radiogenic helium) to 1.1x10**-5 (close to typical MORB reservoir). Repeated sampling in some points during more than 20 years showed stability of helium isotopic composition in time in each of them at any level of 3He/4He values. There is no systematic differences of 3He/4He in samples from surface water sources and deeper intervals of boreholes in the same areas. Universal relationship between isotopic composition of helium and general composition of gas phase is absent either, but the minimum 3He/4He values occurred in methane gas of hydrocarbon deposits, whereas in nitrogen and carbon dioxide gases of helium composition varied (in the latter maximum 3He/4He values have been measured). According to N2/Ar_atm ratio nitrogen gases are atmospheric. In carbonic gas fN2/fNe ratio indicates presence of excessive (non-atmogenic) nitrogen, but the attitude CO2/3He differs from one in MORB. Comparison of helium isotopic composition with its concentration and composition of the main components of gas phase from fluids shows that it is formed under influence of fractionation of components with different solubility in the gas-water system and generation/consumption of reactive gases in the crust. Structural and tectonic elements of the region differ from the spectrum of 3He/4He values. At the pre-Riphean Siberian Platform the mean 3He/4He = (3.6+/-0.9)x10**- 8 is very close to radiogenic one. In the Paleozoic crust of Khangay 3He/4He = (16.3+/-4.6)x10**-8, and the most probable estimate is (12.3+/-2.9)x10**-8. In structures of the eastern flank of the Baikal rift zone (Khentei, Dauria) affected by the Mz-Kz activization 3He/4He values range from 4.4x10**-8 to 2.14x10**-6 (average 0.94x10**-6). Distribution of 3He/4He values across the strike of the Baikal rift zone indicates advective heat transfer from the mantle not only in the rift zone, but also much further to the east. In fluids of the Baikal rift zone range of 3He/4He values is the widest: from 4x10**-8 to 1.1x10**-5. Their variations along the strike of the rift zone are clearly patterned, namely, decrease of 3He/4He values in both directions from the Tunka depression. Accompanied by decrease in density of conductive heat flow and in size of rift basins, this trend indicates decrease in intensity of advective heat transfer from the mantle to peripheral segments of the rift zone. Comparing this trend with data on other continental rift zones and mid-ocean ridges leads to the conclusion about fundamental differences in mechanisms of interaction between the crust and the mantle in these environments.

Noble gases sampled from three cruises in the eastern tropical Atlantic

Oceanic upwelling velocities are too small to be measured directly. Deviations of the He-3/He-4 ratio in the mixed layer from solubility equilibrium provide an indirect means to infer vertical velocities at the base of the mixed layer. This method is applied to the Mauritanian upwelling region for data from three cruises in summer 2006 and winter 2007 and 2008. Diapycnal mixing coefficients are estimated from microstructure measurements, reaching from 10**-3 m**2/s over the shelf break to 10**-5 m**2/s in the open ocean. The resulting upwelling velocities in the onshore region (upto 50 km from the 50 m isobath) are of the order of 2 x 10**-5 m/s}, in agreement with Ekman theory. Further offshore, in some cases the vertical velocities inferred from the helium isotope disequilibrium exceed the values derived from the wind stress curl by one order of magnitude. The Mauritanian coastal area as part of the Canary Current upwelling system belongs to the most productive ocean regions in the world. Nutrient fluxes into the mixed layer (both advective and diffusive) are equivalent to a net community production of about 1 g C/d, and associated heat fluxes vary between 183 +/- 62 W/m**2 in summer and 97 +/- 25 W/m**2 in winter. Regarding the flux into the mixed layer, the contribution of diffusion and advection are of similar magnitude for both heat and nutrients. The upwelling, however, provides the supply of cold and nutrient rich water from below. The large offshore vertical velocities inferred from the helium method are associated with nutrient fluxes of the same order as for the onshore region, and may be responsible for observed patches of high productivity in that area. The offshore heat fluxes due to upwelling and diapycnal mixing are smaller than 70 W/m**2 for all cruises.

Community Structure in Soil Porous System

Diffusion controls the gaseous transport process in soils when advective transport is almost null. Knowledge of the soil structure and pore connectivity are critical issues to understand and modelling soil aeration, sequestration or emission of greenhouse gasses, volatilization of volatile organic chemicals among other phenomena. In the last decades these issues increased our attention as scientist have realize that soil is one of the most complex materials on the earth, within which many biological, physical and chemical processes that support life and affect climate change take place. A quantitative and explicit characterization of soil structure is difficult because of the complexity of the pore space. This is the main reason why most theoretical approaches to soil porosity are idealizations to simplify this system. In this work, we proposed a more realistic attempt to capture the complexity of the system developing a model that considers the size and location of pores in order to relate them into a network. In the model we interpret porous soils as heterogeneous networks where pores are represented by nodes, characterized by their size and spatial location, and the links representing flows between them. In this work we perform an analysis of the community structure of porous media of soils represented as networks. For different real soils samples, modelled as heterogeneous complex networks, spatial communities of pores have been detected depending on the values of the parameters of the porous soil model used. These types of models are named as Heterogeneous Preferential Attachment (HPA). Developing an exhaustive analysis of the model, analytical solutions are obtained for the degree densities and degree distribution of the pore networks generated by the model in the thermodynamic limit and shown that the networks exhibit similar properties to those observed in other complex networks. With the aim to study in more detail topological properties of these networks, the presence of soil pore community structures is studied. The detection of communities of pores, as groups densely connected with only sparser connections between groups, could contribute to understand the mechanisms of the diffusion phenomena in soils.

The Lagrangian description of aperiodic flows: a case study of the Kuroshio Current

This article reviews several recently developed Lagrangian tools and shows how their com- bined use succeeds in obtaining a detailed description of purely advective transport events in general aperiodic flows. In particular, because of the climate impact of ocean transport processes, we illustrate a 2D application on altimeter data sets over the area of the Kuroshio Current, although the proposed techniques are general and applicable to arbitrary time depen- dent aperiodic flows. The first challenge for describing transport in aperiodical time dependent flows is obtaining a representation of the phase portrait where the most relevant dynamical features may be identified. This representation is accomplished by using global Lagrangian descriptors that when applied for instance to the altimeter data sets retrieve over the ocean surface a phase portrait where the geometry of interconnected dynamical systems is visible. The phase portrait picture is essential because it evinces which transport routes are acting on the whole flow. Once these routes are roughly recognised it is possible to complete a detailed description by the direct computation of the finite time stable and unstable manifolds of special hyperbolic trajectories that act as organising centres of the flow.

THE TRANSPORT OF HEAT BY FLOWING GROUNDWATER IN A DISCRETE FRACTURE

In typical theoretical or experimental studies of heat migration in discrete fractures, conduction and thermal dispersion are commonly neglected from the fracture heat transport equation, assuming heat conduction into the matrix is predominant. In this study analytical and numerical models are used to investigate the significance of conduction and thermal dispersion in the plane of the fracture for a point and line sources geometries. The analytical models account for advective, conductive and dispersive heat transport in both the longitudinal and transverse directions in the fracture. The heat transport in the fracture is coupled with a matrix equation in which heat is conducted in the direction perpendicular to the fracture. In the numerical model, the governing heat transport processes are the same as the analytical models; however, the matrix conduction is considered in both longitudinal and transverse directions. Firstly, we demonstrate that longitudinal conduction and dispersion are critical processes that affect heat transport in fractured rock environments, especially for small apertures (eg. 100 μm or less), high flow rate conditions (eg. velocity greater than 50 m/day) and early time (eg. less than 10 days). Secondly, transverse thermal dispersion in the fracture plane is also observed to be an important transport process leading to retardation of the migrating heat front particularly at late time (eg. after 40 days of hot water injection). Solutions which neglect dispersion in the transverse direction underestimate the locations of heat fronts at late time. Finally, this study also suggests that the geometry of the heat sources has significant effects on the heat transport in the system. For example, the effects of dispersion in the fracture are observed to decrease when the width of the heat source expands.

Organic properties of sediments and amino acids in iterstitial waters of ODP Holes 111-677A and 111-678B

The basement at Ocean Drilling Program (ODP) Sites 677 and 678 originated from the Galapagos spreading center of the Costa Rica Rift and has moved about 200 km over the last 6 m.y. (Fig. 1) (Shipboard Scientific Party, 1987, 1988; Scientific Drilling Party, 1987). Sediments about 300 m thick cover basement so young that basal sediments at Sites 677 and 678 have been reheated up to 60?-70?C at Site 677 and altered to limestone and/or chert (Shipboard Scientific Party, 1988). Sediments from both sites indicate (1) a high sedimentation rate (about 48 m/m.y.) and (2) biogenic silica and carbonate as the main constituents of sediments (Table 1) (Shipboard Scientific Party, 1988). Heatflow observations and measurements of interstitial water chemistry around the sites show that Site 677 is in a lower heatflow zone (166 mW/m**2; 1°12.14'N, 83°44.22'W) whereas Site 678 is located in a zone of higher heat flow (250 mW/m**2; 1°13.01'N, 83°43.39'W) (Langseth et al., 1988; Shipboard Scientific Party, 1988). In the flank hydrothermal systems, circulating solution is moving upward through the sedimentary column in zones of higher heat flow while it is moving downward in zones of lower heat flow (Anderson and Skilbeck, 1981). The chemistry of the interstitial waters is modified by several processes such as (1) diagenetic reactions and (2) advective and (3) diffusive transports of dissolved constituents. Analyses of Ca2+ and Mg2+ in interstitial waters from Sites 677 and 678 show that their profiles are mainly controlled by advective transport (Shipboard Scientific Party, 1988). In contrast, the interstitial-water profiles for NH4+, Si, and PO4[3-] are highly affected by reactions in the sediments. Site 677 offers a good opportunity to investigate amino acids in the interstitial waters because sediments of similar compositions have been deposited at constant rates of sedimentation. There are few previous works on amino acid distributions in interstitial waters (Henrichs and Parrington, 1979; Michaelis et al., 1982; Henrichs et al., 1984; Henrichs and Farrington, 1987; Ishizuka et al., 1988). In this chapter, we report (1) Rock-Eval analysis and (2) the composition of total hydrolyzable and dissolved free amino acids (THAA and DFAA, respectively) in the interstitial waters. Our objectives are to discuss (1) the possible origin of organic materials, (2) the characteristics of THAA and DFAA, and (3) their relationships in interstitial waters.

Data compilation of ESOP

The work in this sub-project of ESOP focuses on the advective and convective transforma-tion of water masses in the Greenland Sea and its neighbouring areas. It includes observational work on the sub-mesoscale and analysis of hydrographic data up to the gyre-scale. Observations of active convective plumes were made with a towed chain equipped with up to 80 CTD sensors, giving a horizontal and vertical resolution of the hydrographic fields of a few metres. The observed scales of the penetrative convective plumes compare well with those given by theory. On the mesoscale the structure of homogeneous eddies formed as a result of deep convection was observed and the associated mixing and renewal of the intermediate layers quantified. The relative importance and efficiency of thermal and haline penetrative convection in relation to the surface boundary conditions (heat and salt fluxes and ice cover) and the ambient stratification are studied using the multi year time series of hydro-graphic data in the central Greenland Sea. The modification of the water column of the Greenland Sea gyre through advection from and mixing with water at its rim is assessed on longer time scales. The relative contributions are quantified using modern water mass analysis methods based on inverse techniques. Likewise the convective renewal and the spreading of the Arctic Intermediate Water from its formation area is quantified. The aim is to budget the heat and salt content of the water column, in particular of the low salinity surface layer, and to relate its seasonal and interannual variability to the lateral fluxes and the fluxes at the air-sea-ice interface. This will allow to estimate residence times for the different layers of the Greenland Sea gyre, a quantity important for the description of the Polar Ocean carbon cycle.