Fluid compositions and mineralogy of precipitates from Mid Atlantic Ridge hydrothermal vents at 4°48'S

The effect of volcanic activity on submarine hydrothermal systems has been well documented along fast- and intermediate-spreading centers but not from slow-spreading ridges. Indeed, volcanic eruptions are expected to be rare on slow-spreading axes. Here we report the presence of hydrothermal venting associated with extremely fresh lava flows at an elevated, apparently magmatically robust segment center on the slow-spreading southern Mid-Atlantic Ridge near 5°S. Three high-temperature vent fields have been recognized so far over a strike length of less than 2 km with two fields venting phase-separated, vapor-type fluids. Exit temperatures at one of the fields reach up to 407°C, at conditions of the critical point of seawater, the highest temperatures ever recorded from the seafloor. Fluid and vent field characteristics show a large variability between the vent fields, a variation that is not expected within such a limited area. We conclude from mineralogical investigations of hydrothermal precipitates that vent-fluid compositions have evolved recently from relatively oxidizing to more reducing conditions, a shift that could also be related to renewed magmatic activity in the area. Current high exit temperatures, reducing conditions, low silica contents, and high hydrogen contents in the fluids of two vent sites are consistent with a shallow magmatic source, probably related to a young volcanic eruption event nearby, in which basaltic magma is actively crystallizing. This is the first reported evidence for direct magmatic-hydrothermal interaction on a slow-spreading mid-ocean ridge.

Chemical, and Sr and Li isotopic composition of pore fluids and sediments from ODP Leg 170 and Leg 205 sites

Pore fluid and sediment chemical and isotopic data were obtained for samples from Ocean Drilling Program (ODP) Leg 205 Sites 1253, 1254, and 1255 in the Costa Rica subduction zone. The chemical and isotopic data reported here were generated in our shore-based laboratories to complement shipboard inorganic geochemical data. Li isotopic analyses were carried out by L.-H. Chan at Louisiana State University (USA). The data reported herein include fluoride, bromide, rubidium, cesium, and barium concentrations; Li and Sr isotopic compositions in pore fluids; and Rb, Cs, and Ba concentrations in representative bulk sediments. The data also include new pore fluid fluoride and bromide concentrations from corresponding ODP Leg 170 Sites 1039, 1040, and 1043. O.M. Saether's Site 1039 and 1040 fluoride concentration data are shown for comparison. Basal sediment fluoride concentrations and Li and Sr isotope ratios at both Sites 1253 and 1039 show reversals that approach modern seawater values. Br/Cl ratios are, however, conservative throughout the sediment section at Sites 1039 and 1253. The observed sharp F and Br concentration maxima, Rb and K concentration minima, the most radiogenic 87Sr/86Sr ratios, and highest 7Li values along the décollement and fracture zone (Sites 1040, 1043, 1254, and 1255) strengthen the evidence obtained during Leg 170 that a deeply sourced fluid, originating from fluid-rock reactions at ~150°C and corresponding to between 10 and 15 km depth, is transporting solutes to the ocean.

Major ion concentrations and chlorine isotope composition of water-soluble chloride and structurally bound chloride in DSDP/ODP serpentinized ultramafic rocks

Eight DSDP/ODP cores were analyzed for major ion concentrations and d37Cl values of water-soluble chloride (d37Clwsc) and structurally bound chloride (d37Clsbc) in serpentinized ultramafic rocks. This diverse set of cores spans a wide range in age, temperature of serpentinization, tectonic setting, and geographic location of drilled serpentinized oceanic crust. Three of the cores were sampled at closely spaced intervals to investigate downhole variation in Cl concentration and chlorine isotope composition. The average total Cl content of all 86 samples is 0.26±0.16 wt.% (0.19±0.10 wt.% as water-soluble Cl (Xwsc) and 0.09±0.09 wt.% as structurally bound Cl (Xsbc)). Structurally bound Cl concentration nearly doubles with depth in all cores; there is no consistent trend in water-soluble Cl content among the cores. Chlorine isotope fractionation between the structurally bound Cl**- site and the water-soluble Cl**- site varies from -1.08? to +1.16?, averaging to +0.21?. Samples with negative fractionations may be related to reequilibration of the water-soluble chloride with seawater post-serpentinite formation. Six of the cores have positive bulk d37Cl values (+0.05? to +0.36?); the other two cores (173-1068A (Leg-Hole) and 84-570) have negative bulk d37Cl values (-1.26? and -0.54?). The cores with negative d37Cl values also have variable Cl**-/SO4**2- ratios, in contrast to all other cores. The isotopically positive cores (153-920D and 147-895E) show no isotopic variation with depth; the isotopically negative core (173-1068A) decreases by ~1? with depth for both the water-soluble and structurally bound Cl fractions. Non-zero bulk d37Cl values indicate Cl in serpentinites was incorporated during original hydration and is not an artifact of seawater infiltration during drilling. Cores with positive d37Cl values are most likely explained by open system fractionation during hydrothermal alteration, with preferential incorporation of 37Cl from seawater into the serpentinite and loss of residual light Cl back to the ocean. Fluid / rock ratios were probably low as evidenced by the presence of water-soluble salts. The two isotopically negative cores are characterized by a thick overlying sedimentary package that was in place prior to serpentinization. We believe the low d37Cl values of these cores are a result of hydration of ultramafic rock by infiltrating aqueous pore fluids from the overlying sediments. The resulting serpentinites inherit the characteristic negative d37Cl values of the pore waters. Chlorine stable isotopes can be used to identify the source of the serpentinizing fluid and ultimately discern chemical and tectonic processes involved in serpentinization.

Interstitial water d13C and sediment geochemistry of ODP Leg 166 sites

A review of interstitial water samples collected from Sites 1003-1007 of the Bahamas Transect along with a shore-based analysis of oxygen and carbon isotopes, minor and trace elements, and sediment chemistry are presented. Results indicate that the pore-fluid profiles in the upper 100 meters below seafloor (mbsf) are marked by shifts between 20 and 40 mbsf that are thought to be caused by changes in sediment reactivity, sedimentation rates, and the influence of strong bottom currents that have been active since the late Pliocene. Pore-fluid profiles in the lower Pliocene-Miocene sequences are dominated by diffusion and do not show significant evidence of subsurface advective flow. Deeper interstitial waters are believed to be the in situ fluids that have evolved through interaction with sediments and diffusion. Pore-fluid chemistry is strongly influenced by carbonate recrystallization processes. Increases in pore-fluid Cl- and Na+ with depth are interpreted to result mainly from carbonate remineralization reactions that are most active near the platform margin. A lateral gradient in detrital clay content observed along the transect, leads to an overall lower carbonate reactivity, and enhances preservation of metastable aragonite further away from the platform margin. Later stage burial diagenesis occurs at slow rates and is limited by the supply of reactive elements through diffusion.

Chemistry and age determination of basement fluids from Juan de Fuca Ridge

The geochemical implications of thermally driven flow of seawater through oceanic crust on the mid-ocean ridge flank have been examined on a well-studied 80 km transect across the eastern flank of the Juan de Fuca Ridge at 48°N, using porewater and basement fluid samples obtained on ODP Leg 168. Fluid flow is recognised by near-basement reversals in porewater concentration gradients from altered values in the sediment section to seawater-like values in basaltic basement. In general, the basement fluids become more geochemically evolved with distance from the ridge and broadly follow basement temperature which ranges from not, vert, similar16° to 63°C. Although thermal effects of advective heat exchange are only seen within 20 km east of where basement is exposed near the ridge crest, chemical reactivity extends to all sites. Seawater passing through oceanic crust has reacted with basement rocks leading to increases in Ca2+ and decreases in alkalinity, Mg2+, Na+, K+, SO42- and delta18O. Sr isotope exchange between seawater and oceanic crust off axis is unequivocally demonstrated with endmember 87Sr/86Sr ~ 0.707. Evidence of more evolved fluids is seen at sites where rapid upwelling of fluids through sediments occurs. Chlorinities of the basement fluids are consistent with post-glacial seawater and thus a short residence time in the crust. Rates of lateral flow have been by estimated by modelling porewater sulphate gradients, using Cl as a glacial chronometer, and from radiocarbon dating of basal fluids. All three methods reveal fluid flow with 14C ages less than 10,000 yr and particle velocities of ~1-5 m/yr, in agreement with thermally constrained volumetric flow rates through a ~600 m thick permeable layer of ~10% porosity. Delta(element)/Delta(heat) extraction ratios are similar to values for ridge-crest hydrothermal systems.

Strontium concentrations and isotopic compositions of anhydrites from the TAG active mound

Sr isotope analyses have been conducted on anhydrite samples from the TAG (Trans-Atlantic Geotraverse) active hydrothermal mound (26°08?N, Mid-Atlantic Ridge) that have previously been shown to exhibit two distinct patterns of REE behavior when normalized to TAG end-member hydrothermal fluid. Despite differences in REE patterns, the Sr isotope data indicate that all the anhydrites precipitated from fluids with a similar range of hydrothermal fluid and seawater components, and all but one were seawater-dominated (52%-75%). Speciation calculations using the EQ3/6 software package for geochemical modeling of aqueous systems suggest that the REE complexation behavior in different fluid mixing scenarios can explain the variations in the REE patterns. Anhydrites that exhibit relatively flat REE patterns [(La_bs)/(Yb_bs) = 0.8-2.0; subscript bs indicates normalization to end-member black smoker hydrothermal fluid] and a small or no Eu anomaly [(Eu_bs)/(Eu*_bs) = 0.8-2.0] are inferred to have precipitated from mixes of end-member hydrothermal fluid and cold seawater. REE complexes with hard ligands (e.g., fluoride and chloride) are less stable at low temperatures and trivalent Eu has an ionic radius similar to that of Ca2+ and the other REE, and so they behave coherently. In contrast, anhydrites that exhibit slight LREE-depletion [(La_bs)/(Yb_bs) = 0.4-1.4] and a distinct negative anomaly [(Eu_bs)/(Eu*_bs) = 0.2-0.8] are inferred to have precipitated from mixes of end-member hydrothermal fluid and conductively heated seawater. The LREE depletion results from the presence of very stable LREE chloro-complexes that effectively limit the availability of the LREE for partitioning into anhydrite. Above 250°C, Eu is present only in divalent form as chloride complexes, and discrimination against Eu2+ is likely due to both the mismatch in ionic radii between Eu2+ and Ca2+, and the strong chloro-complexation of divalent Eu which promotes stability in the fluid and inhibits partitioning of Eu2+ into precipitating anhydrite. These variations in REE behavior attest to rapid fluctuations in thermal regime, fluid flow and mixing in the subsurface of the TAG mound that give rise to heterogeneity in the formation conditions of individual anhydrite crystals.

(Table 1) Bulk and grain density, porosity, and p- and s-wave velocity characteristics of ODP Hole 163-990A sediments

Seismic velocities in rocks are influenced by the properties of the solid, the pore fluid, and the pore space. Cracks dramatically affect seismic velocities in rocks; their influence on the effective elastic moduli of rocks depends on their shape and concentration. Thin cracks (or fractures) substantially lower the moduli of a rock relative to the effect of spherical voids (or vesicles), and lower moduli are reflected by lower P- and S-wave velocities. The objective of this research is to determine the types and concentrations of cracks and their influence on the seismic properties of subaerially erupted basalts drilled from Hole 990A on the Southeast Greenland margin during Ocean Drilling Program Leg 163. Ellipsoidal cracks are used to model the voids in the rocks. The elastic moduli of the solid (grains) are also free parameters in the inverse modeling procedure. The apparent grain moduli reflect a weighted average of the moduli of the constituent minerals (e.g., plagioclase, augite, and clay minerals). The results indicate that (1) there is a strong relationship between P-wave velocity and porosity, suggesting a similarity of pore shape distributions, (2) the distribution of crack types within the massive, central region of aa flows from Hole 990A is independent of total porosity, (3) thin cracks are the first to be effectively sealed by alteration products, and (4) grain densities (an alteration index) and apparent grain moduli of the basalt samples are directly related.

(Table 2) Oxygen, carbon, and strontium isotopic composition of carbonate-brucite material from the Lost City hydrothermal field

The isotopic (dD, d18O, d13C, and 87Sr/86Sr) and geochemical characteristics of hydrothermal solutions from the Mid-Atlantic Ridge and the material of brucite-carbonate chimneys at the Lost City hydrothermal field at 30°N, MAR, were examined to assay the role of the major factors controlling the genesis of the fluid and hydrothermal chimneys of the Lost City field. The values of dD and d18O in fluid samples indicates that solutions at the Lost City field were produced during the serpentinization of basement ultramafic rocks at temperatures higher than 200°C and at relatively low fluid/rock ratios (<1). The active role of serpentinization processes in the genesis of the Lost City fluid also follows from the results of the electron-microscopic studying of the material of hydrothermal chimneys at this field. The isotopic (d18O, d13C, and 87Sr/86Sr) and geochemical (Sr/Ca and REE) signatures indicate that, before its submarine discharging at the Lost City field, the fluid filtered through already cold altered outer zones of the Atlantis Massif and cooled via conductive heat loss. During this stage, the fluid could partly dissolve previously deposited carbonates in veins cutting serpentinite at the upper levels of the Atlantis Massif and the carbonate cement of sedimentary breccias underlying the hydrothermal chimneys. Because of this, the age of modern hydrothermal activity at the Lost City field can be much younger than 25 ka.

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.

(Table 2) Mercury and selenium concentration in the brain stem of wild polar bears (Ursus maritimus), east Greenland

Polar bears (Ursus maritimus) are exposed to high concentrations of mercury because they are apex predators in the Arctic ecosystem. Although mercury is a potent neurotoxic heavy metal, it is not known whether current exposures are of neurotoxicological concern to polar bears. We tested the hypotheses that polar bears accumulate levels of mercury in their brains that exceed the estimated lowest observable adverse effect level (20 µg/g dry wt) for mammalian wildlife and that such exposures are associated with subtle neurological damage, as determined by measuring neurochemical biomarkers previously shown to be disrupted by mercury in other high-trophic wildlife. Brain stem (medulla oblongata) tissues from 82 polar bears subsistence hunted in East Greenland were studied. Despite surprisingly low levels of mercury in the brain stem region (total mercury = 0.36 ± 0.12 µg/g dry wt), a significant negative correlation was measured between N-methyl-D-aspartate (NMDA) receptor levels and both total mercury (r = -0.34, p < 0.01) and methylmercury (r = -0.89, p < 0.05). No relationships were observed among mercury, selenium, and several other neurochemical biomarkers (dopamine-2, gamma-aminobutyric acid type A, muscarinic cholinergic, and nicotinic cholinergic receptors; cholinesterase and monoamine oxidase enzymes). These data show that East Greenland polar bears do not accumulate high levels of mercury in their brain stems. However, decreased levels of NMDA receptors could be one of the most sensitive indicators of mercury's subclinical and early effects.

Composition and age of ores and bottom sediments accumulated within and near the Rainbow hydrothermal field

Results of direct geological and geochemical observations of the modern Rainbow hydrothermal field (Mid-Atlantic Ridge, 36°14'N; 33°54'W) carried out from the deep-sea manned Mir submersibles during Cruises 41 and 42 of the R/V Akademik Mstislav Keldysh in 1998-1999 and data of laboratory studies of collected samples are under consideration in the paper. The field lacks neovolcanic rocks and the axial part of the rift is filled in with a serpentinite protrusion. In this field there occur metalliferous sediments, as well as active and relict sulfide edifices composed of sulfide minerals; pyrrhotite, chalcopyrite, isocubanite, sphalerite, marcasite, pyrite, bornite, chalcosine, digenite, magnetite, anhydrite, rare troilite, wurtzite, millerite, and pentlandite have been determined. Sulfide ores are characterized by concentric-zoned textures. During in situ measurements during 35 minutes temperature of hydrothermal fluids was varying within a range from 250 to 350°C. Calculated chemical and isotopic composition of hydrothermal fluid shows elevated concentrations of Cl, Ni, Co, CH4, and H2. Values of d34S of H2S range from +2.4 to +3.1 per mil, of d13C of CH4 from -15.2 to -11.2 per mil, and d13C of CO2 from +1.0 to -4.0 per mil. Fluid inclusions are homogenized at temperatures from 140 to 360°C, whereas salinity of the fluid varies from 4.2 to 8.5 wt %. d34S values of sulfides range from +1.3 to +12.5 per mil. 3He/4He ratio in mineral-forming fluid contained in the fluid inclusions from sulfides of the Rainbow field varies from 0.00000374 to 0.0000101. It is shown that hydrothermal activity in the area continues approximately during 100 ka. It is assumed that the fluid and sulfide edifices contain components from the upper mantle. A hypothesis of phase separation of a supercritical fluid that results in formation of brines is proposed. Hydrothermal activity is related to the tectonic, not volcanic, phase of the Mid-Atlantic Ridge evolution.

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.

Methane, helium isotope and hydrogen data from the Nibelungen hydrothermal area, Mid Atlantic Ridge

Hydrothermal emission of mantle helium appears to be directly related to magma production rate, but other processes can generate methane and hydrogen on mid-ocean ridges. In an on-going effort to characterize these processes in the South Atlantic, the flux and distribution of these gases were investigated in the vicinity of a powerful black smoker recently discovered at 8°17.9' S, 13°30.4' W. The vent lies on the shoulder of an oblique offset in the Mid-Atlantic Ridge and discharges high concentrations of methane and hydrogen. Measurements during expeditions in 2004 and 2006 show that the ratio of CH4 to 3He in the neutrally buoyant plume is quite high, 4 x 10**8. The CTD stations were accompanied by velocity measurements with lowered acoustic Doppler current profilers (LADCP), and from these data we estimate the methane transport to have been 0.5 mol/sec in a WSW-trending plume that seems to develop during the ebb tidal phase. This transport is an order of magnitude greater than the source of CH4 calculated from its concentration in the vent fluid and the rise height of the plume. From this range of methane fluxes, the source of 3He is estimated to be between 0.14 and 1.2 nmol/sec. In either case, the 3He source is significantly lower than expected from the spreading rate of the Mid-Atlantic Ridge. From the inventory of methane in the rift valley adjacent to the vent, it appears that the average specific rate of oxidation is 2.6 to 23/yr, corresponding to a turnover time between 140 and 16 days. Vertical profiles of methane in the surrounding region often exhibited Gaussian-like distributions, and the variances appear to increase with distance from the vent. Using a Gaussian plume model, we obtained a range of vertical eddy diffusivities between 0.009 and 0.08 m2m2/sec. These high values may be due to tidally driven internal waves across the promontory on which the vent is located.

Euechinoidea and Cidaroidea respond differently to ocean acidification

The impact of the chemical changes in the ocean waters due to the increasing atmospheric CO2 depends on the ability of an organism to control extracellular pH. Among sea urchins, this seems specific to the Euechinoidea, sea urchins except Cidaroidea. However, Cidaroidea survived two ocean acidification periods: the Permian-Trias and the Cretaceous-Tertiary crises. We investigated the response of these two sea urchin groups to reduced seawater pH with the tropical cidaroid Eucidaris tribuloides, the sympatric euechinoid Tripneustes ventricosus and the temperate euechinoid Paracentrotus lividus. Both euechinoid showed a compensation of the coelomic fluid pH due to increased buffer capacity. This was linked to an increased concentration of DIC in the coelomic fluid and thus of bicarbonate ions (most probably originating from the surrounding seawater as isotopic signature of the carbon -delta 13C- was similar). On the other hand, the cidaroid showed no changes within the coelomic fluid. Moreover, the delta 13C of the coelomic fluid did not match that of the seawater and was not significantly different between the urchins from the different treatments. Feeding rate was not affected in any species. While euechinoids are able to regulate their extracellular acid-base balance, many questions are still unanswered on the costs of this capacity. On the contrary, cidaroids do not seem affected by a reduced seawater pH. Further investigations need to be undertaken to cover more species and physiological and metabolic parameters in order to determine if energy trade-offs occur and how this mechanism of compensation is distributed among sea urchins.

(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