OS43A-2028
Numerical Modeling of the Hydrothermal System at East Pacific Rise 9°50’N Including Anhydrite Precipitation

Thursday, 17 December 2015
Poster Hall (Moscone South)
Kannikha Parameswari Kolandaivelu, Virginia Polytechnic Institute and State University, Blacksburg, VA, United States and Robert P Lowell, Virginia Polytechnic Institute and State University, Geosciences, Blacksburg, VA, United States
Abstract:
To better understand the effects of anhydrite precipitation on mid-ocean ridge hydrothermal systems, we conducted 2-D numerical simulations of two-phase hydrothermal circulation in a NaCl-H2O fluid at the East Pacific Rise 9°50’N. The simulations were constrained by key observational thermal data and seismicity that suggests the fluid flow is primarily along axis with recharge focused into a small zone near a 4th order discontinuity. The simulations considered an open-top square box with a fixed seafloor pressure of 25 MPa, and nominal seafloor temperature of 10 °C. The sides of the box were assumed to be impermeable and insulated. We considered two models: a homogeneous model with a permeability of 10-13 m2 and a heterogeneous model in which layer 2A extrusives were given a higher permeability. Both models had a fixed bottom temperature distribution and initial porosity of 0.1. Assuming that anhydrite precipitation resulted from the decrease in solubility with increasing temperature as downwelling fluid gets heated, we calculated the rate of porosity decrease and sealing times in each cell at certain time snapshots in the simulations. The results showed that sealing would occur most rapidly in limited regions near the base of the high-temperature plumes, where complete sealing could occur on decadal time scales. Though more detailed analysis is needed, it appeared that the areas of rapid sealing would likely have negligible impact on the overall circulation pattern and hydrothermal vent temperatures. The simulations also indicated that sealing due to anhydrite precipitation would occur more slowly at the margins of the ascending plumes. The sealing times in the deep recharge zone determined in these simulations were considerably greater than estimated from 1D analytical calculations, suggesting that with a 2D model, focused recharge at the EPR 9°50’N site may occur, at least on a decadal time scale.