Mantle convection and evolution with growing continents

Abstract

[1] We present a three-dimensional spherical shell numerical model for chemical differentiation and redistribution of incompatible elements in a convective Earth's mantle heated mostly from within by U, Th, and K and slightly from below. The evolution-model equations guarantee conservation of mass, momentum, energy, angular momentum, and four sums of the number of atoms of the pairs 238 U-206 Pb, 235 U-207 Pb, 232 Th-208 Pb, and 40 K-40 Ar. The pressure-and temperature-dependent viscosity is supplemented by a viscoplastic yield stress, s y . The lithospheric viscosity is partly imposed to mimic its increase by dehydration of oceanic lithosphere and other effects. Also, the asthenosphere is generated not only by the distribution of temperature and melting temperature, but essentially by the profiles of water solubility and water abundance. Therefore we introduced a radial viscosity profile factor describing that behavior. However, the focus of this paper is the episodic growth of continents and oceanic plateaus. As a complement, the differentiation generates the depleted MORB mantle (DMM) which predominates immediately beneath the lithosphere. Our continents are not artificially imposed on the surface of the spherical shell, but instead they evolve by the interplay between chemical differentiation and convection/mixing. No restrictions are imposed regarding number, size, form, and distribution of continents. However, oceanic plateaus that impinge upon a continent have to be united with it. This mimics the accretion of terranes. The numerical results show an episodic growth of the total mass of the continents and display a plausible time history for the laterally averaged surface heat flow, qob, and the Rayleigh number, Ra. We use our model to explore a moderate region of Ra-s y parameter space. We find Earth-like continent distributions in a central part of the Ra-s y space we explored. We identified a Ra-s y region where the calculated total continental volume is very close to the observed value; another Ra-s y region where the Urey number, Ur, is close to the accepted value; a third Ra-s y area where surface heat flow is very close to the present-day observed mean global heat flow using typical abundances of the heat-producing elements. It is remarkable that these different acceptable Ra-s y regions share a common overlap area, where Earth-like behavior is simultaneously fulfilled.Although the convective flow patterns and the chemical differentiation of oceanic plateaus are coupled, the evolution of time-dependent Rayleigh number, Ra t , is relatively well predictable and the stochastic parts of the Ra t (t) curves are small. Regarding the time distribution of juvenile growth rates of the total mass of the continents, predictions are possible only in the first epoch of the evolution, presumed that the initial conditions are given. Later on, the distribution of the continental growth episodes is increasingly stochastic. Independent of the varying individual runs, our model shows that the total mass of the presen...

Real episodic growth of continental crust or artifact of preservation? A 3‐D geodynamic model

“…If the computed model temperature, T , approaches the melting temperature, T m , in a sufficiently large volume, i.e., if the condition of equation () is fulfilled, then chemical differentiation occurs. The new CC material will be welded to an existent continent if an oceanic plate has carried it to a continental margin [ Walzer and Hendel , ]. Our present model also incorporates those buoyancy forces that are generated by thermally induced deflections of the phase boundaries of the olivine‐wadsleyite transition and of the ringwoodite‐perovskite transition.…”

“…In accordance with Komiya ], we assert that continental growth is determined by both efflux from the mantle and recycling of CC material back into the mantle. The corresponding tracer method is described in Walzer and Hendel ]. The results of our numerical model also show quiescent intervals of magmatic evolution (Figure ) in agreement with observation.…”

Palaeoproterozoic supercontinents and global evolution: correlations from core to atmosphere

“…1i) (Condie 1998;Campbell & Allen 2008), possibly linked to large-scale mantle overturn events (Condie 2000;Rino et al 2004), and the cessation of, or decrease in, subduction flux (Silver & Behn 2008a). Recent models integrating chemical differentiation with mantle convection also predict the episodicity of juvenile crust formation (Walzer & Hendel 2008), as do large-scale mantle overturn events that are thought to take place on a timescale of several hundred million years (Davies 1995).…”

The history of supercontinents and oceans from the standpoint of thermochemical mantle convection

“…In order to overcome the contradictions associated with the geochemical data, we decided (1) to expand the thermal model of mantle convection by including the most significant chemical factors in it (Christensen and Yuen, 1984;Kobozev and Myasnikov, 1987;Christensen, 1997;McNamara and Zhong, 2004;Walzer and Hendel, 2008). We also planned (2) to simulate thermochemical convection mode, which can assemble the supercontinents.…”