TY - JOUR

T1 - Influence of lattice thermal conductivity on thermal convection with strongly temperature-dependent viscosity

AU - Yanagawa, Tomohiko K.B.

AU - Nakada, Masao

AU - Yuen, David A.

N1 - Funding Information:
Acknowledgments. The authors wish to thank Paul J. Tackley and Masaki Ogawa for constructive and helpful reviews. The numerical calculations were performed on the Fujitsu VPP5000 Supercomputer in the data processing center of the National Astronomical Observatory of Japan. This research is supported by the geophysics and CSEDI programs of the National Science Foundation.

PY - 2005

Y1 - 2005

N2 - To examine the effects of temperature-dependent lattice thermal conductivity (lattice-k) on the thermal convection with a temperature-dependent viscous fluid, particularly on the upper thermal boundary layer, we have studied numerically 2-D thermal convective flows with both constant thermal conductivity (constant-k) and lattice-k models. Numerical experiments with large viscosity contrasts, greater than a million, produce a cooler and thinner upper thermal boundary layer for the lattice-k compared with those for the constant-k, implying that thermal convection with lattice-k produces a much sharper boundary between the lithosphere and asthenosphere. The differences between the constant-k and lattice-k can be reasonably explained by the following two causes: (i) the decreasing lattice-k with depth increases an effective Rayleigh number around the bottom of the thermal boundary layer, and (ii) the distribution of lattice-k and uniform vertical heat flux within the thermal boundary layer determine the temperature distribution. The predicted sharper boundary, i.e. sharper vertical viscosity gradient near the bottom of the lithosphere, may play an important role on controlling the amount of lithospheric deformation associated with the downwelling.

AB - To examine the effects of temperature-dependent lattice thermal conductivity (lattice-k) on the thermal convection with a temperature-dependent viscous fluid, particularly on the upper thermal boundary layer, we have studied numerically 2-D thermal convective flows with both constant thermal conductivity (constant-k) and lattice-k models. Numerical experiments with large viscosity contrasts, greater than a million, produce a cooler and thinner upper thermal boundary layer for the lattice-k compared with those for the constant-k, implying that thermal convection with lattice-k produces a much sharper boundary between the lithosphere and asthenosphere. The differences between the constant-k and lattice-k can be reasonably explained by the following two causes: (i) the decreasing lattice-k with depth increases an effective Rayleigh number around the bottom of the thermal boundary layer, and (ii) the distribution of lattice-k and uniform vertical heat flux within the thermal boundary layer determine the temperature distribution. The predicted sharper boundary, i.e. sharper vertical viscosity gradient near the bottom of the lithosphere, may play an important role on controlling the amount of lithospheric deformation associated with the downwelling.

UR - http://www.scopus.com/inward/record.url?scp=21844466520&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=21844466520&partnerID=8YFLogxK

U2 - 10.1186/BF03351802

DO - 10.1186/BF03351802

M3 - Article

AN - SCOPUS:21844466520

SN - 1343-8832

VL - 57

SP - 15

EP - 28

JO - Earth, Planets and Space

JF - Earth, Planets and Space

IS - 1

ER -