No abstract

[1] Mid-ocean ridges host vigorous hydrothermal systems that remove large quantities of heat from the oceanic crust. Inferred Nusselt numbers (Nu), which are the ratios of the total heat flux to the heat flux that would be transported by conduction alone, range from 8 to several hundred. Such vigorous convection is not fully described by most numerical models of hydrothermal circulation. A major difficulty arises at high Nu from the numerical solution of the temperature equation. To avoid classical numerical artifacts such as nonphysical oscillatory behavior and artificial diffusion, we implement the Multidimensional Positive Definite Advection Transport Algorithm (MPDATA) technique, which solves the temperature equation using an iterated upwind corrected scheme. We first validate the method by comparing results for models with uniform fluid properties in closed-and open-top systems to existing solutions with Nu $20. We then incorporate realistic fluid properties and run models for Nu up to 50-60. Solutions are characterized by an unstable bottom thermal boundary layer where thermal instabilities arise locally. The pattern of heat extraction is periodic to chaotic. At any Nu > $13 the venting temperatures in a given plume are chaotic and oscillate from $350°to 450°C. Individual plumes can temporarily stop short of the surface for intervals ranging from tens to hundreds of years at times when other plumes vent with an increased flow rate. The solutions also display significant recirculation, and as a result large areas of downflow are relatively warm with temperatures commonly exceeding 150°C at middepths. Our results have important implications for mid-ocean ridge hydrothermal systems and suggest the following: (1) The reaction zones of mid-ocean ridge hydrothermal systems are enlarged by thermal instabilities that migrate laterally toward upflow zones. This will substantially increase the volume of rock involved in chemical reactions compared to steady state configurations. (2) Hydrothermal discharge can stop temporarily as zones of venting are dynamically replaced by zones of seawater recharge. (3) Anhydrite precipitation occurring at temperatures exceeding $150°C will likely occur throughout a large portion of recharge zone and will not necessarily clog downflow pathways as efficiently as has been recently inferred.

[1] Mid-ocean ridge hydrothermal systems are known to vent fluids with salinities substantially different from seawater as a result of phase separation and segregation of the resulting vapor and brine phases. Time series of vent temperature and salinity (chlorinity) show that some black-smoker vent fields such as the Main Endeavour Field on the Juan de Fuca Ridge have vented fluids with salinities well below seawater for over a decade, which raises important questions concerning the fate of brines in these systems. One widely accepted model is that high-density brines formed by supercritical phase separation sink to the base of hydrothermal systems, leading to the development of a two-layer system in which a recirculating brine layer underlies a single-pass seawater cell. We first present theoretical arguments to constrain the dynamics of such a deep brine layer in a system still undergoing phase separation, and we conclude that if brines are stored in a basal layer, they are unlikely to convect because they will be stably stratified. One consequence of this result is that the brine layer beneath black smoker systems has to be thin (<10 m) to match the high heat fluxes. However, estimates of the rate at which brines are accumulating in the crust below the main field on the Endeavour segment of the Juan de Fuca Ridge suggest that the brine layer is likely at least 100 m thick. To resolve this apparent paradox, we propose an alternative model. We argue that interfacial tensions between fluid and solid phases will likely favor the segregation of vapor into the main fractures and brine into the smaller fissures and backwaters. This allows the vapor to flow efficiently through the system and transport large heat fluxes while most of the porosity in the lower part of the system fills with brines. It is generally believed that the pressure gradients in mid-ocean ridge hydrothermal systems are close to cold hydrostatic. At the high temperatures and pressures characteristic of the deeper parts of these systems, brines with salinities as high as 20 wt % NaCl have densities around 800-900 kg m À3 and will be buoyant in a cold-hydrostatic system. Rather than sinking to the base of the system, it is possible that brines produced by supercritical phase separation rise slowly until they reach a level of neutral buoyancy as they cool or enter high-permeability regions in which the pressure gradients decrease.

[1] The chemical composition of mid-ocean ridge hydrothermal vent fluids is thought to reflect conditions within a deep-seated reaction zone. Although temperature and pressure conditions within this region are key parameters that characterize the subseafloor hydrothermal regime and the cooling of mid-ocean ridges, they are poorly constrained. In this paper, we developed a model in which high-temperature, vapor-type (low-salinity) vent fluid silica (Si) and chlorine (Cl) concentrations can be used to define lines in pressuretemperature space whose intersection is used to estimate conditions at the top of the reaction zone, under the simplifying assumption that Si and Cl reflect a common point of equilibration. We apply this model to various basaltic-hosted mid-ocean ridge sites. Results suggest a minimal variation in inferred temperatures, ranging from 415 to 445°C. This lends support to the fluxibility model in which upwelling hydrothermal plumes rise at temperatures that maximize the energy flux. Quartz precipitation due to reequilibration during upflow tends to lower temperature and pressure estimates and can artificially indicate shallower transition from reaction to upflow zone. However, maximum equilibration pressures are site-dependent and compare well with depth to magma chamber imaged by seismic studies. This suggests that vapors circulate close to magma chambers and is difficult to reconcile with models in which mid-ocean ridge hydrothermal circulation occurs in two layers with a substantial layer of convecting brine. Accordingly, equilibration pressure predicted by our model can also be used to infer the depth of the magma chamber at sites where seismic data are not available but where vapor-like fluids have been collected and analyzed.Components: 5327 words, 5 figures, 1 table.Keywords: hydrothermal vents; mid-ocean ridges; reaction zone; silica; chlorine; AMC depth.