Geodynamic modeling of thermal structure of subduction zones

Abstract

During subduction processes, slabs continuously have heat exchange with the ambient mantle, including both conduction and advection effects. The evolution of slab thermal structure is one of the dominant factors in controlling physical and chemical property changes in subduction zones. It also affects our understanding of many key geological processes, such as mineral dehydration, rock partial melting, arc volcanism, and seismic activities in subduction zones. There are mainly two ways for studying thermal structure of subduction zones with geodynamic models: analytical model and numerical model. Analytical model provides insights into the most dominant controlling physical parameters on the thermal structure, such as slab age, velocity and dip angle, shear stress and thermal conductivity, etc. Numerical model can further deal with more complicated environments, such as viscosity change in the mantle wedge, coupling process between slabs and the ambient mantle, and incorporation of petrology and mineralogy. When applying geodynamic modeling results to specific subduction zones on the Earth, there are many factors which may complicate the process, therefore it is difficult to precisely constrain the thermal structure of subduction zones. With the development of new quantitative methods in geophysics and geochemistry, we may obtain more observational constraints for thermal structure of subduction zones, thus providing more reasonable explanations for geological processes related to subduction zones.