Free-surface formulation of mantle convection—I. Basic theory and application to plumes

In mantle convection models, the top surface is traditionally approximated as a free-slip boundary, and the dynamic topography is obtained by assuming that the normal stress on the free-slip boundary is compensated instantaneously through surface deformation. It has already been shown that this approximation is-'valid for long-wavelength topography. Based on both viscous and viscoelastic models with a free surface, we have found that the characteristic time for topographic growth is comparable to the timescales of mantle convection (∼10^6 year) for short and intermediate wavelengths (10^3 km or less) and/or a high effective lithospheric viscosity (> 10^(24) Pa s). This suggests that the topography is history-dependent under these conditions and that a free-surface formulation is required to study the topography at these wavelengths. We have developed a new Eulerian finite-element technique to model a free surface. Since the technique is based upon an undeformable Eulerian grid, this enables us to study long-term, free-surface dynamics in the presence of evolving buoyancy. We have compared numerical with analytic solutions of viscous relaxation for fixed buoyancy problems. As long as the magnitude of topography is much smaller than the wavelength, we find that the finite-element method is very accurate, with relative errors of less than 1 per cent. This numerical technique can be applied to a variety of geophysical problems with free surfaces. In applying this technique to dynamic models of mantle plumes, we find that surface relaxation retards the topography at intermediate and short wavelengths and produces a smoother topography, compared with topography from free-slip calculations. This reduced topography has a significant influence on the geoid at the corresponding wavelengths. Moreover, free-surface models, by allowing vertical motion on the free surface, yield a hotter lithosphere over ascending plumes than models with free-slip boundaries.

[1]  Sean C. Solomon,et al.  The evolution of impact basins - Viscous relaxation of topographic relief. [for lunar surface modeling] , 1982 .

[2]  N. A. Haskell The Motion of a Viscous Fluid Under a Surface Load. Part II , 1936 .

[3]  Wing Kam Liu,et al.  Lagrangian-Eulerian finite element formulation for incompressible viscous flows☆ , 1981 .

[4]  D. L. Anderson Lithosphere, asthenosphere, and perisphere , 1995 .

[5]  G. Davies Role of the lithosphere in mantle convection , 1988 .

[6]  R. McConnell,et al.  Viscosity of the mantle from relaxation time spectra of isostatic adjustment , 1968 .

[7]  B. Parsons,et al.  The relationship between surface topography, gravity anomalies, and temperature structure of convection , 1983 .

[8]  B. Hager,et al.  Geoid Anomalies in a Dynamic Earth , 1984 .

[9]  M. Gurnis,et al.  Holographic measurements of surface topography in laboratory models of mantle hotspots , 1989 .

[10]  G. Schubert,et al.  The viscosity of the earth's mantle , 1976 .

[11]  H. Harder,et al.  NUMERICAL SIMULATION OF THERMAL CONVECTION WITH MAXWELLIAN VISCOELASTICITY , 1991 .

[12]  A. Poliakov,et al.  Diapirism and topography , 1992 .

[13]  Patrick Wu,et al.  Rheology of the Upper Mantle: A Synthesis , 1993, Science.

[14]  Shijie Zhong,et al.  Controls on trench topography from dynamic models of subducted slabs , 1994 .

[15]  H. J. Melosh,et al.  The dynamical origin of subduction zone topography , 1980 .

[16]  W. J. Morgan,et al.  A nonlinear rheology model for mid‐ocean ridge axis topography , 1990 .

[17]  B. Hager Subducted slabs and the geoid: Constraints on mantle rheology and flow , 1983 .

[18]  B. Hager,et al.  — a H 03 Dynamically Supported Geoid Highs over Hotspots : Observation and Theory , 2008 .

[19]  N. Ribe,et al.  Observations of flexure and the geological evolution of the Pacific Ocean basin , 1980, Nature.

[20]  C. Eloy,et al.  Free-surface formulation of mantle convection—II. Implication for subduction-zone observables , 1996 .

[21]  Jean Braun,et al.  A numerical method for solving partial differential equations on highly irregular evolving grids , 1995, Nature.

[22]  N. A. Haskell The Motion of a Viscous Fluid Under a Surface Load , 1935 .

[23]  C. Froidevaux,et al.  Geoid heights and lithospheric stresses for a dynamic Earth. , 1984 .

[24]  Bradford H. Hager,et al.  Conman: vectorizing a finite element code for incompressible two-dimensional convection in the Earth's mantle , 1990 .