Open AccessA thermomechanical numerical model for crustal accretion of medium to fast spreading mid‐ocean ridges
Author(s) -
Machetel P.,
Garrido C. J.
Publication year - 2009
Publication title -
geochemistry, geophysics, geosystems
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.928
H-Index - 136
ISSN - 1525-2027
DOI - 10.1029/2008gc002270
Subject(s) - geology , sill , geophysics , thermal , vorticity , mid ocean ridge , hydrothermal circulation , thermal diffusivity , stream function , mechanics , petrology , mantle (geology) , thermodynamics , vortex , seismology , physics
We propose a new thermomechanical numerical model of mid‐ocean ridge accretion aimed at investigating asymmetric spreading rates, diverse configurations of lens and sill magma injections, crystallization and depth, and on‐ and off‐axis patterns of hydrothermal cooling. The numerical algorithm iteratively resolves temperature and motion equations until it reaches a stationary solution. The motion equation was written in a vorticity‐stream function formalism, with boundary and internal conditions applied to the stream function to impose the style of magma injection. Unlike in previous models, our model does not assume an a priori shape for the temperature field, which is initiated by an initial half‐space cooling according to the left and right spreading rates. Complex patterns of hydrothermal cooling are simulated by enhanced thermal diffusivity. The model succeeds in describing the dynamic and thermal effects of spreading rates, the style of magma intrusion, and the hydrothermal cooling. Accurate descriptions of these are essential to study the cooling histories of crustal rocks and geophysical observables.