Continental emergence and growth on a cooling earth

Isostasy considerations are connected to a 1-D model of mantle differentiation due to pressure release partial melting to obtain a model for the evolution of the relative sea level with respect to the continent during the earth secular cooling. In this context, a new mechanism is derived for the selective exhumation of exposed ancient cratons. The model results in a quantitative scenario for sea-level fall due to the changing thicknesses of the oceanic basaltic crust and its harzburgite residual layer as a function of falling mantle temperature. It is also shown that the buoyancy of the harzburgite root of a stabilized continental craton has an important effect on sea-level and on the isostatic readjustment and exhumation of exposed continental surface during the earths secular cooling. The model does not depend on the usual assumption of constant continental freeboard and crustal thickness and its application is not restricted to the post-Archaean. It predicts large-scale continental emergence near the end of the Archaean and the early Proterozoic. This provides an explanation for reported late Archaean emergence and the subsequent formation of late Archaean cratonic platforms and early Proterozoic sedimentary basins. For a period of secular cooling of 3.8 Ga, corresponding to the length of the geological record, the model predicts a fall of the ocean floor of some 4 km or more. For a constant ocean depth, this implies a sea-level fall of the same magnitude. A formula is derived that allows for an increasing ocean depth due to either the changing ratio of continental with respect to oceanic area, or to a possible increase of the oceanic volume during the geological history. Increasing ocean depth results in a later emergence of submarine ancient geological formations compared to the case when ocean depth is constant. Selective exhumation is studied for the case of constant ocean depth. It is shown that for this case, early exposed continental crust can be exhumed to a lower crustal depth, which explains the relative vertical displacement of low-grade- with respect to high-grade terrain. Increasing ocean depth is not expected to result in diminished exhumation. © 2000 Elsevier Science B.V. All rights reserved

Early formation and long-term stability of continents resulting from decompression melting in a convecting mantle

The origin of stable old continental cratonic roots is still debated. We present numerical modelling results which show rapid initial formation during the Archaean of continental roots of ca. 200 km thick. These results have been obtained from an upper mantle thermal convection model including differentiation by pressure release partial melting of mantle peridotite. The upper mantle model includes time-dependent radiogenic heat production and thermal coupling with a heat reservoir representing the Earths lower mantle and core. This allows for model experiments including secular cooling on a time-scale comparable to the age of the Earth. The model results show an initial phase of rapid continental root growth of ca. 0.1 billion year, followed by a more gradual increase of continental volume by addition of depleted material produced through hot diapiric, convective upwellings which penetrate the continental root from below. Within ca. 0.6 Ga after the start of the experiment, secular cooling of the mantle brings the average geotherm below the peridotite solidus thereby switching off further continental growth. At this time the thickness of the continental root has grown to ca. 200 km. After 1 Ga of secular cooling small scale thermal instabilities develop at the bottom of the continental root causing continental delamination without breaking up the large scale layering. This delaminated material remixes with the deeper layers. Two more periods, each with a duration of ca. 0.5 Ga and separated by quiescent periods were observed when melting and continental growth was reactivated. Melting ends at 3 Ga. Thereafter secular cooling proceeds and the compositionally buoyant continental root is stabilized further through the increase in mechanical strength induced by the increase of the temperature dependent mantle viscosity. Fluctuating convective velocity amplitudes decrease to below 10 mma Õ 1 and the volume average temperature of the sub-continental convecting mantle has decreased ca. 340 K after 4 Ga. Surface heatflow values decrease from 120 to 40 mW m Õ 2 during the 4 Ga model evolution. The surface heatflow contribution from an almost constant secular cooling rate was estimated to be 6 mW m Õ 2, in line with recent observational evidence. The modelling results show that the combined effects of compositional buoyancy and strong temperature dependent rheology result in continents which overall remain stable for a duration longer than the age of the Earth. Tracer particles have been used for studying the patterns of mantle differentiation in greater detail. The observed ( p, T, F, t)-paths are consistent with proposed stratification and thermo-mechanical history of the depleted continental root, which have been inferred from mantle xenoliths and other upper mantle samples. In addition, the particle tracers have been used to derive the thermal age of the modelled continental root, defined by a hypothetical closing temperature. © 2000 Elsevier Science B.V. All rights reserved

Stability and growth of continental shields in mantle convection models including recurrent melt production

The long-term growth and stability of compositionally layered continental upper mantle has been investigated by numerical modelling. We present the first numerical model of a convecting mantle including differentiation through partial melting resulting in a stable compositionally layered continental upper mantle structure. This structure includes a continental root extending to a depth of about 200 km. The model covers the upper mantle including the crust and incorporates physical features important for the study of the continental upper mantle during secular cooling of the Earth since the Archaean. Among these features are: a partial melt generation mechanism allowing consistent recurrent melting, time-dependent non-uniform radiogenic heat production, and a temperature- and pressure-dependent rheology. The numerical results reveal a long-term growth mechanism of the continental compositional root. This mechanism operates through episodical injection of small diapiric upwellings from the deep layer of undepleted mantle into the continental root which consists of compositionally distinct depleted mantle material. Our modelling results show the layered continental structure to remain stable during at least 1.5 Ga. After this period mantle differentiation through partial melting ceases due to the prolonged secular cooling and small-scale instabilities set in through continental delamination. This stable period of 1.5 Ga is related to a number of limitations in our model. By improving on these limitations in the future this stable period will be extended to more realistic values. Ó 1998 Elsevier Science B.V. All rights reserved

Can Lower Mantle Slab-like Seismic Anomalies be Explained by Thermal Coupling Between the Upper and Lower Mantles?

Below subduction zones, high resolution seismic tomographic models resolve fast anomalies that often extend into the deep lower mantle. These anomalies are generally interpreted as slabs penetrating through the 660-km seismic discontinuity, evidence in support of whole-mantle convection. However, thermal coupling between two ow systems separated by an impermeable interface might provide an al ternative explanation of the tomographic results. We have tested this hypothesis within the context of an axisymmet ric model of mantle convection in which an impermeable boundary is imposed at a depth of 660 km. When an increase in viscosity alone is imposed across the impermeable interface, our results demonstrate the dominant role of mechanical coupling between shells, producing lower mantle upwellings (downwellings) below upper mantle downwellings (upwellings). However, we find that the effect of mechanical coupling can be significantly weakened if a narrow low viscosity zone exists beneath the 660-km discontinuity. In such a case, both thermally induced `slabs' in the lower mantle and thermally activated plumes that rise from the upper/lower mantle boundary are observed even though mass transfer between the shells does not exist

Continental emergence and growth on a cooling earth

Isostasy considerations are connected to a 1-D model of mantle differentiation due to pressure release partial melting to obtain a model for the evolution of the relative sea level with respect to the continent during the earth secular cooling. In this context, a new mechanism is derived for the selective exhumation of exposed ancient cratons. The model results in a quantitative scenario for sea-level fall due to the changing thicknesses of the oceanic basaltic crust and its harzburgite residual layer as a function of falling mantle temperature. It is also shown that the buoyancy of the harzburgite root of a stabilized continental craton has an important effect on sea-level and on the isostatic readjustment and exhumation of exposed continental surface during the earths secular cooling. The model does not depend on the usual assumption of constant continental freeboard and crustal thickness and its application is not restricted to the post-Archaean. It predicts large-scale continental emergence near the end of the Archaean and the early Proterozoic. This provides an explanation for reported late Archaean emergence and the subsequent formation of late Archaean cratonic platforms and early Proterozoic sedimentary basins. For a period of secular cooling of 3.8 Ga, corresponding to the length of the geological record, the model predicts a fall of the ocean floor of some 4 km or more. For a constant ocean depth, this implies a sea-level fall of the same magnitude. A formula is derived that allows for an increasing ocean depth due to either the changing ratio of continental with respect to oceanic area, or to a possible increase of the oceanic volume during the geological history. Increasing ocean depth results in a later emergence of submarine ancient geological formations compared to the case when ocean depth is constant. Selective exhumation is studied for the case of constant ocean depth. It is shown that for this case, early exposed continental crust can be exhumed to a lower crustal depth, which explains the relative vertical displacement of low-grade- with respect to high-grade terrain. Increasing ocean depth is not expected to result in diminished exhumation. © 2000 Elsevier Science B.V. All rights reserved

Passive margin evolution, initiation of subduction and the Wilson cycle

We have constructed finite element models at various stages of passive margin evolution, in which we have incorporated the system of forces acting on the margin, depth-dependent rheological properties and lateral variations across the margin. We have studied the interrelations between age-dependent forces, geometry and rheology, to decipher their net effect on the state of stress at passive margins. Lithospheric flexure induced by sediment loading dominates the state of stress at passive margins. This study has shown that if after a short evolution of the margin (time span a few tens of million years) subduction has not yet started, continued aging of the passive margin alone does not result in conditions more favourable for transformation into an active margin. Although much geological evidence is available in support of the key role small ocean basins play in orogeny and ophiolite emplacement, evolutionary frameworks of the Wilson cycle usually are cast in terms of opening and closing of wide ocean basins. We propose a more limited role for large oceans in the Wilson cycle concept

Various mechanisms to induce present-day shallow flat subduction and implications for the younger Earth: a numerical parameter study

Shallow flat subduction is a relatively common feature at present-day subduction zones. Several mechanisms to explain this feature have been proposed, and can be subdivided into three groups: overthrusting of the subducting plate, subduction of a plume-generated oceanic plateau, and slab suction forces. We developed a numerical model to investigate these mechanisms and tested it through a comparison of the model results with the observations of the Peru flat slab where all three mechanisms seem to be contributing. The ratio of contributions of overthrusting continent to plateau subduction is in the range of 1:1 to 1:2, and the role of plate suction forces is likely to be significant. By applying the overthrusting continent and plateau subduction mechanisms separately, we were able to determine the most important model parameters for each of the mechanisms. Overthrusting easily results in flat subduction, and the flat slab length is primarily a function of slab age, overriding plate motion and mantle viscosity. An oceanic plateau is much less likely to cause flat subduction, and favorable conditions for flat subduction include a young slab age, long-lived plateau buoyancy after subduction, a strong mantle, and addition of slab suction forces that are large enough to further reduce the subduction dip angle, once the plateau initiates this flattening. Furthermore, we found that even though today flat subduction can be explained with the dominant model parameters within a reasonable range, for a slightly hotter, younger Earth, these flat subduction conditions are much less favorable, and so this style of subduction was probably not present in the past. This contradicts earlier predictions that flat subduction was a more wide-spread phenomenon in the early stages of plate tectonics in a younger earth

On the role of subducting oceanic plateaus in the development of shallow flat subduction

Oceanic plateaus, aseismic ridges or seamount chains all have a thickened crust and their subduction has been proposed as a possible mechanism to explain the occurrence of flat subduction and related absence of arc magmatism below Peru, Central Chile and at the Nankai Trough (Japan). Their extra compositional buoyancy could prohibit the slab from sinking into the mantle. With a numerical thermochemical convection model, we simulated the subduction of an oceanic lithosphere that contains an oceanic crustal plateau of 18-km thickness. With a systematic variation, we examined the required physical parameters to obtain shallow flat subduction. Metastability of the basaltic crust in the eclogite stability field is of crucial importance for the slab to remain buoyant throughout the subduction process. In a 44-Ma-old subducting plate, basalt must be able to survive a temperature of 600–700 jC to keep the plate buoyant sufficiently long to cause a flat-slab segment. We found that the maximum yield stress in the slab must be limited to about 600 MPa to allow for the necessary bending to the horizontal. Young slabs show flat subduction for larger parameter ranges than old slabs, since they are less gravitationally unstable and show less resistance against bending. Hydrous weakening of the mantle wedge area and lowermost continent are required to allow for the necessary deformation of a change in subduction style from steep to flat. The maximum flat slab extent is about 300 km, which is sufficient to explain the observed shallow flat subduction near the Nankai Trough (Japan). However, additional mechanisms, such as active overthrusting by an overriding continental plate, need to be invoked to explain the flatslab segments up to 500 km long below Peru and Central Chile