An internally consistent inverse model to calculate ridge-axis hydrothermal fluxes

Fluid and chemical fluxes from high-temperature, on-axis, hydrothermal systems at mid-ocean ridges have been estimated in a number of ways. These generally use simple mass balances based on either vent fluid compositions or the compositions of altered sheeted dikes. Here we combine these approaches in an internally consistent model. Seawater is assumed to enter the crust and react with the sheeted dike complex at high temperatures. Major element fluxes for both the rock and fluid are calculated from balanced stoichiometric reactions. These reactions include end-member components of the minerals plagioclase, pyroxene, amphibole, chlorite and epidote along with pure anhydrite, quartz, pyrite, pyrrhotite, titanite, magnetite, ilmenite and ulvospinel and the fluid species H2O, Mg2+, Ca2+, Fe2+, Na+, Si4+, H2S, H+ and H2. Trace element abundances (Li, B, K, Rb, Cs, Sr, Ba, U, Tl, Mn, Cu, Zn, Co, Ni, Pb and Os) and isotopic ratios (Li, B, O, Sr, Tl, Os) are calculated from simple mass balance of a fluid-rock reaction. A fraction of the Cu, Zn, Pb, Co, Ni, Os and Mn in the fluid after fluid-rock reaction is allowed to precipitate during discharge before the fluid reaches the seafloor. S-isotopes are tied to mineralogical reactions involving S-bearing phases. The free parameters in the model are the amounts of each mineralogical reaction that occurs, the amounts of the metals precipitated during discharge, and the water-to-rock ratio. These model parameters, and their uncertainties, are constrained by: (i) mineral abundances and mineral major element compositions in altered dikes from ODP Hole 504B and the Pito and Hess Deep tectonic windows (EPR crust); (ii) changes in dike bulk-rock trace element and isotopic compositions from these locations relative to fresh MORB glass compositions; and (iii) published vent fluid compositions from basalt-hosted high-temperature ridge axis hydrothermal systems. Using a numerical inversion algorithm, the probability density of different model parameter sets has been computed and thus the probability of different fluid and chemical fluxes. Most data can be fit by the model within their uncertainty. The entire dataset is best-fit with a water-to-rock mass ratio between 1.3 and 2.1 (~1 to 1.5 x10**13 kg yr-1) implying a substantial fraction of the magmatic (latent) heat available to drive the axial hydrothermal system is extracted by these systems. Many element fluxes are better constrained than in previous studies (e.g., Sr: 2 to 7 x10**8 moles yr-1; Ca: 2 to 7 x10**11 moles yr-1). Future developments will use experimental data to further constrain the model.