Outline
Comptes Rendus

Geodynamics / Géodynamique
Global thermal and dynamical perturbations due to Cretaceous mantle avalanche
Comptes Rendus. Géoscience, Volume 335 (2003) no. 1, pp. 91-97.

Abstracts

The numerical models of mantle convection agree to depict avalanches behaviour according to the level of endothermicity of the spinel → perovskite phase change. Their potential effects on the global thermal and dynamical states of the mantle have been computed thanks to a numerical code, which takes into account both the 400-km exothermic and the 660-km endothermic phase changes. The cycle followed by the avalanches is: local layering, destabilization of the 660-km thermal layer, travelling and spreading on the core, and reappearing of the local layering. Therefore, mantle convection is characterized by quiet periods of partial layering embedded in catastrophic events. During the avalanche, the amplitude of the surface velocity is multiplied by two, which would imply an enhanced plate tectonic and ridge activities. The global thermal effects of the avalanche are compatible with a high mantle temperature and an acceleration of Earth's rotation during the Cretaceous. They also offer a coherent explanation to locate the origin of mantle plumes both within the CMB and just below the transition zone.

Les modèles numériques de convection dans le manteau de la Terre sont en accord pour décrire des phénomènes d'avalanches qui dépendent du caractère endothermique de la transition de phase spinelle → pérovskite. Les effets potentiels de ces avalanches sur l'état thermique et la dynamique du manteau ont été calculés en prenant en compte à la fois la transition de phase exothermique à 400 km et la transition endothermique à 660 km. Les avalanches suivent un cycle de mise en couche partielle de la convection, puis de déstabilisation de la couche thermique à 660 km, de voyage, puis d'écrasement de l'avalanche sur le noyau, puis de restauration de la mise en couche partielle. Ainsi, la convection dans le manteau se traduit par une succession de périodes calmes délimitées par des événements catastrophiques. Pendant les avalanches, l'amplitude des vitesses observées en surface est multipliée par deux, ce qui devrait se traduire par une activité tectonique accrue et des dorsales plus actives. Les effets globaux des avalanches sont compatibles avec une période de plus forte température dans le manteau au moment du Crétacé, ainsi qu'avec une accélération de la rotation de la Terre. L'évolution du champ de température pendant ces événements permet également de proposer que certains panaches du manteau puissent naı̂tre au niveau de la CMB et que d'autres proviennent de dessous la zone de transition.

Metadata
Received:
Accepted:
Published online:
DOI: 10.1016/S1631-0713(03)00007-5
Keywords: Mantle convection, Phase changes, Temperature, Length of day, Heat flux, Tectonics, Plumes
Mot clés : Convection du manteau, Changement de phase, Température, Longueur du jour, Flux de chaleur, Tectonique, Panache

Philippe Machetel 1

1 Laboratoire de tectonophysique, UMR 5568, CNRS/UM2 CC056, Institut des sciences de la Terre, de l'Eau et de l'Espace de Montpellier (Isteem), place Eugène-Bataillon, 34095 Montpellier cedex, France
@article{CRGEOS_2003__335_1_91_0,
     author = {Philippe Machetel},
     title = {Global thermal and dynamical perturbations due to {Cretaceous} mantle avalanche},
     journal = {Comptes Rendus. G\'eoscience},
     pages = {91--97},
     publisher = {Elsevier},
     volume = {335},
     number = {1},
     year = {2003},
     doi = {10.1016/S1631-0713(03)00007-5},
     language = {en},
}
TY  - JOUR
AU  - Philippe Machetel
TI  - Global thermal and dynamical perturbations due to Cretaceous mantle avalanche
JO  - Comptes Rendus. Géoscience
PY  - 2003
SP  - 91
EP  - 97
VL  - 335
IS  - 1
PB  - Elsevier
DO  - 10.1016/S1631-0713(03)00007-5
LA  - en
ID  - CRGEOS_2003__335_1_91_0
ER  - 
%0 Journal Article
%A Philippe Machetel
%T Global thermal and dynamical perturbations due to Cretaceous mantle avalanche
%J Comptes Rendus. Géoscience
%D 2003
%P 91-97
%V 335
%N 1
%I Elsevier
%R 10.1016/S1631-0713(03)00007-5
%G en
%F CRGEOS_2003__335_1_91_0
Philippe Machetel. Global thermal and dynamical perturbations due to Cretaceous mantle avalanche. Comptes Rendus. Géoscience, Volume 335 (2003) no. 1, pp. 91-97. doi : 10.1016/S1631-0713(03)00007-5. https://comptes-rendus.academie-sciences.fr/geoscience/articles/10.1016/S1631-0713(03)00007-5/

Version originale du texte intégral

1 Introduction

For approximately ten years, the numerical models of convection have shown the importance of the main mantle phase transitions for global convection studies. The olivine-spinel phase change at 400 km depth is exothermic, which facilitates convection through upper mantle. On the contrary, the spinel–perovskite phase change at 660 km is endothermic, which can slow down or even stop the rising and the sinking convection currents at the discontinuity level [15].

The numerical models, which have taken these phase changes into account, agree to depict intermittent pattern of mantle convection according to the level of endothermicity of the spinel → perovskite phase change. Such a regime, which has been found with independent 2D, spherical [15,16], 3D, spherical [14,22], and 2D and 3D Cartesian geometry [8,25], is very robust, since it has been encountered with internal heating rates varying from 0 to chondritic, viscosity stratification and pressure and temperature dependence of viscosity [2]. Within a 2D or 3D spherical, inelastic framework, whole mantle convection is prevailing for a Clapeyron slope lower than -2×106 Pa K-1. Conversely, it is two-layered for more endothermic reaction. Thus, for a sinking (resp. rising) convection current, the density jump at 660 km resists to flow penetration. Indeed, as the Clapeyron's pressure, necessary to complete the phase change, corresponds, at this colder (resp. higher) temperature, to a depth deeper (resp. shallower) than the actual surrounding density jump, the density contrast may block convection that may spread horizontally.

The agreement between these numerical models does not mean that the answer about the mantle structure is definitively known, but that all these models give a coherent representation of the phase-change effects in terms of physic of fluid. The geoid retrieving by circulation models [23] and the tomographic imaging of deep slabs plunging into the lower mantle argue in favour of a present one- or partially-layered mantle structure [18,24]. This present state results from the long process that has been acting since the early Earth. At this epoch, the huge quantity of heat accumulated during the accretion and the higher radioactive release of energy were heating the mantle, which was necessarily hotter and less viscous than today. Then, the Rayleigh number of this primitive earth was therefore higher than the present one, which has been shown to induce layering of convection [20,21,27,28]. A continuous evolution, from early layered mantle convection to whole mantle, or partly layered mantle structure, implies that we are presently in an avalanche-like regime or, at least, that it has been encountered during the past.

2 Mechanism of the mantle avalanches

The resolution of the full set of convection equations, including the Navier–Stokes equation, but also the non-linear energy equation, is very expensive in terms of computer time and generates huge quantities of data, which cannot be recorded in details. In order to deepen the preliminary results obtained during previous studies [2,3], a new series of axisymmetrical, spherical convection models have been computed to quantify the potential effects of mantle avalanches on the global thermal and dynamical states of the mantle. The numerical code takes into account an exothermic phase change at 400 km, an endothermic phase change at 660 km, a temperature- and pressure-dependent viscosity and a viscosity jump of 30 between the upper and lower mantle. A more detailed description of the physical hypothesis and of the numerical method can be found in our previous work [2].

Fig. 1 displays the various types of avalanches that have been encountered during the computation. An animated sequence can be found at the web address: http://www.isteem.univ-montp2.fr/PERSO/machetel/anim.html. The occurrence of the avalanches is a temporally chaotic process, which illustrates the compromise between a global thermal state of the mantle and its local recent history. Fig. 1 (and animation) shows that an avalanche may occur at a fixed location (Fig. 1a), cease and reappear at almost the same place (Fig. 1b). However, it also can start as an unexpected event over a weak zone created by plumes rising from lower mantle (Fig. 1c).

Fig. 1

Computed temperature field at 210 Myr (computation time) (panel a), 624 Myr (panel b) and 810 Myr (panel c). Computation has been restarted from (time arbitrary set to 0) after the occurrence of a first very powerful avalanche. The theoretical model includes compressibility; internal and basal heating; a viscosity jump of 30 between the upper and lower mantle; pressure and temperature dependent viscosity; an endothermic phase change at 660 km and an exothermic phase change at 400 km. Panel a shows that an avalanche may occur at a fix location, cease and reappear at an almost contiguous place (panel b). It also can also start, as an unexpected event, over a weak zone destabilized by a hot lower mantle rising plume (panel c).

Champ de température obtenu après 210 Ma (a), 624 Ma (b) et 810 Ma (c). Le calcul a été réinitialisé (origine arbitraire du temps) après une première avalanche très puissante. Le modèle numérique prend en compte la compressibilité du manteau, un mode de chauffage mixte radioactif et basal, un contraste de viscosité de 30 entre le manteau supérieur et le manteau inférieur, un changement de phase endothermique à 660 km de profondeur et un changement de phase exothermique à 400 km de profondeur. La partie (a) montre que les avalanches peuvent survenir de manière répétitive au même emplacement, cesser, puis réapparaı̂tre presque au même endroit (b). Elles peuvent également démarrer de manière inattendue à partir d'une zone de faiblesse, rendue instable par la présence d'un panache chaud du manteau inférieur (c).

In all the cases, the mantle avalanches follow cycles based on four main phases.

1. A local layering of convection occurs, which is illustrated by the stopping of the upper mantle sinking currents. Simultaneously, the insulation due to the layering prevents the heat from escaping, which induces a warming and a softening of the lower mantle. Simultaneously, convection cools and hardens the upper mantle. Both effects result in the growing of a thermal layer at the 660-km level. A new avalanche will shortly be generated by this increasingly unstable evolution (Fig. 2A) [4].

Fig. 2

The global sequence that characterizes the avalanche phenomena is: (A) local layering of convection; (B) destabilization of the 670 boundary and brutal advection of cold upper mantle convective structure (the avalanche itself); (C) spreading of the avalanche on the core and surface arrival of return flows; (D) restoration of local layering which marks the end of the avalanche.

Les événements qui caractérisent le phénomène d'avalanche se déroulent ainsi : (A) mise en couches partielle de la convection entre manteau supérieur et manteau inférieur ; (B) déstabilisation de la couche limite thermique à 670 km et advection brutale dans le manteau inférieur de matière froide en provenance du manteau supérieur ; (C) arrivée de l'avalanche au niveau du noyau et des courants de retour au niveau de la surface ; (D) le retour d'une mise en couche partielle marque la fin de l'avalanche.

2. The destabilization of the 660-km thermal layer, which has occurred, is now accelerated by the brutal advection of cold upper mantle material (the avalanche itself). The surrounding upper mantle is drawn toward the avalanche epicentre, while hot rising currents, originating from lower mantle, cross upward the 660-km discontinuity to compensate the effect of the avalanche (Fig. 2B). This mechanism generates violent, rising widely distributed plumes, originating just below the 660-km phase-change level.

3. The sinking avalanche travels the lower mantle and, after a few tens of million years, spreads on the core by sweeping the (softer) CMB boundary layer. Conversely, the return flows provoke the re-heating of the upper mantle and local decreases of the viscosity, which amplify the viscosity contrasts through the 660-km phase change. The approximate synchronization of the avalanche spreading onto the core and of the return flow emergence at the surface is due to a focalization of the sinking avalanche, while the return flows are diffuse over all the 660-km phase-change spherical surface (Fig. 2C).

4. The local layering reappears, while no more upper mantle material sinks into the lower mantle (or at low rate), which also stops (or reduces) the return flows. As the avalanche is no longer fed by cold upper-mantle cold matter, the internal heating and the core heat reheat the lower mantle and reduce the CMB thermal gradient (Fig. 2D). A new cycle is ready to start.

3 Thermal consequences of the avalanche

In a recent work, we have monitored in details the evolution of the mantle temperature during an avalanche. The return flows rise the upper mantle temperature by about 50 to 70 °C, while the lower mantle is simultaneously cooled by the avalanche [12]. It is interesting to note that this process has probably acted as a regulator for the global mantle temperature through the long-term history of the Earth. Indeed, it is probable that as early as the avalanche regime has been established in the past, the thermal layer at the 660-km phase-change level has been recursively destroyed as soon as it was becoming unstable. Within this context, the mantle convection is characterized by quiet periods of partial layering embedded in catastrophic events. Fig. 3 (top) displays the flux of matter that has crossed the 660-km phase change (solid line), and the bottom heat flux versus time (dashed line). During the avalanche, the bottom heat flux remains strong, as a consequence of steep thermal gradient at the CMB level [2]. However, it decreases to reach a low level during the quiet period. Bottom panel of Fig. 3 allows comparing the heat fluxes that simultaneously cross the upper surface and the bottom. Both heat fluxes display the same decreasing trend during the quiet period, but their global shapes are very different, since the signal seems smoother for the bottom heat flux than for the surface heat flux. The bottom heat flux varies in response to the arrival of the cold material of the avalanche, which lowers mantle temperature, and consequently increases the heat flux through the bottom surface. The upper surface heat flux displays a more complex behaviour, because two different kinds of lower mantle return flows act to reheat the upper mantle. Those of the first type start, from the 660-km phase-change level, simultaneously with the beginning of the avalanche, for mass conservation reasons. This can be deciphered on the bottom part of Fig. 3, around time 600 Myr, when the surface heat flux increases in response to the starting of a new avalanche. The delay between the surface heat flow and the flux through the 660-km phase change corresponds to the transfer time of return flow from 660 km to the surface. The explanation of the increase of heat flux around 350 Myr is less obvious. It is due to the arrival at the surface of plumes created at the CMB by the sweeping effect of the sinking avalanche. Indeed, when the cold material of the avalanche arrives on the CMB surface, it creates a strong temperature gradient that results in violent plume rising from the CMB. We are going to verify this hypothesis on a longer computation time, in order to check this behaviour of the surface heat flux for several avalanches. The geomagnetic field reversal stopping during the superchron periods may probably be explained by the temperature perturbations of the CMB, occurring with the spreading of the avalanche onto the core [7]. This event may be embedded by superplume activity, which is due to return flows starting from transition zone. Then, the plumes reach the surface just before the beginning of the superchron. Conversely, the superplume activity, due to the spreading of the avalanche on the CMB thermal layer, reaches the surface after the end of the superchron.

Fig. 3

Top: flux of matter through the 660-km phase change (solid line) and heat flux through the bottom surface (dashed line); bottom: comparison between the bottom and surface heat fluxes.

Haut : flux de matière au travers du changement de phase des 660 km (ligne continue) et flux de chaleur à travers la surface inférieure. Bas : comparaison des flux de chaleur au travers de la surface supérieure et de la surface inférieure.

4 Effects of the avalanches on surface velocity

The physical principles that govern convection rely on energy conservation, but also on mass conservation. This implies that avalanches have not only thermal, but also dynamical consequences, which can be seen in Fig. 4. In order to monitor this effect, the maximum of the horizontal velocity at the surface and the mean value of the amplitude of the velocity on the sphere have been computed and drawn versus time. The surface horizontal velocity is very sensitive to the avalanche occurrence. Indeed, its amplitude is twice during the paroxysmal phase. Furthermore, Fig. 4 shows that it is not a local effect, but a global one, since, even after an averaging over the spherical surface, the mean value of the module of the horizontal velocity displays the same trend. Finally, there is an almost perfect synchronization between the evolution of the velocity surface and the mass transfer through the 660-km phase change. This is an important result, since it implies that an avalanche, when it occurs, must go with a simultaneously enhanced plate-tectonic and ridge activity.

Fig. 4

Time evolution of the maximum of surface velocity Vtmax (thin solid line), flux of matter crossing the 660-km phase change (dotted line) and evolution of the mean value of the surface velocity (heavy solid line).

Évolution temporelle de la vitesse horizontale de surface (ligne fine continue), flux de matière au travers du changement de phase des 660 km (ligne pointillée) et évolution de la vitesse moyenne de surface (ligne continue épaisse).

5 Summary and discussion

During our previous studies [12], we have concluded that a high mantle temperature during the Cretaceous was compatible with the global thermal effect of an avalanche. This results was deduced from isostasy considerations and from ODP-sample Na8/Fe8 temperature geochemistry, which both support independently evidence for a mantle thermal high around 125 Myr with amplitude (50 °C) and duration (30 Myr) comparable to those obtained from theoretical convection models.

Furthermore, the potential impact of global thermal anomalies due to an avalanche has also been discussed in term of ability to accelerate the Earth rotation, since the global thermal fluctuations change the depth of the phase changes through the Clapeyron pressure–temperature relationship [13].

These two results have been widely discussed within the context of the large-scale events that have occurred during the Cretaceous: Strong True Polar Wander event [1,17], superchron occurrence and magnetic field reversal perturbations [6], possible earth rotation acceleration [26], enhanced large scale tectonic and ridge activity [5], lower mantle contaminations of Pacific and Atlantic MORB [10,11], evidence for anomalous temperature in hot spots and MORB [9,19].

Within this context, the mantle avalanches offer a coherent explanation for both CMB and transition zone origins for plumes. Such a phenomenon would also be coherent with a global increase of the plate tectonic activity (ridge accretion and continental tearing) as a direct consequence of the surface velocity during the avalanche.

A more philosophical result of this work could be to emphasize the richness of the mantle behaviours, which would consist in alternating quiet periods of partial layering and catastrophic avalanche events. The coupling of the surface and bottom heat fluxes, the thermal regime of the mantle, the origin of plumes, and the surface activity would obviously depend on the instantaneous style of regime with long-term cycles through the Earth's history. We have to keep in mind that this richness has been obtained by including the phase change physics within the convection but without ad-hoc assumptions. On the other hand, the real Earth displays such (and may be worse) complexities that are surprisingly matched by the consequences of mantle avalanches.

Acknowledgements

The author thanks René Guiraud, Vincent Courtillot, and Jean Besse for useful comments and discussions. This work was supported by the CNRS/INSU programme ‘Intérieur de la Terre’.


References

[1] J. Besse, V. Courtillot, Apparent and true polar wander and the geometry of the geomagnetic field in the last 200 million years, J. Geophys. Res. (in press)

[2] D. Brunet; P. Machetel Large-scale tectonic features induced by mantle avalanches with phase, temperature, and pressure lateral variations of viscosity, J. Geophys. Res., Volume 103 (1998), pp. 4929-4945

[3] D. Brunet; P. Machetel Slab weakening by the exothermic olivine–spinel phase change, Geophys. Res. Lett., Volume 25 (1998), pp. 3231-3234

[4] S.L. Butler; W.R. Peltier On scaling relations in time-dependent mantle convection and the heat transfer constraint on layering, J. Geophys. Res., Volume 105 (2000), pp. 3175-3208

[5] V. Courtillot; C. Jaupart; I. Manighetti; P. Tapponnier; J. Besse On causal links between flood basalts and continental breakup, Earth Planet. Sci. Lett., Volume 166 (1999), pp. 177-195

[6] Y. Gallet; G. Hulot Stationary and non-stationary behaviour within the geomagnetic polarity time scale, Geophys. Res. Lett., Volume 24 (1997), pp. 1875-1878

[7] G.A. Glatzmaier; R.S. Coe; L. Hongre; P.H. Roberts The role of the Earth's mantle in controlling the frequency of geomagnetic reversals, Nature, Volume 401 (1999), pp. 885-890

[8] S. Honda; S. Balachandar; D.A. Yuen; D. Reuteler Three-dimensional mantle dynamics with an endothermic phase transition, Geophys. Res. Lett., Volume 20 (1993), pp. 221-224

[9] E. Humler; C. Langmuir; V. Daux Depth versus age: new perspectives from the chemical compositions of ancient crust, Earth Planet. Sci. Lett., Volume 173 (1999), pp. 7-23

[10] P.E. Janney; P.R. Castillo Geochemistry of Mesozoic Pacific mid-ocean ridge basalt: constraints on melt generation and the evolution of the Pacific upper mantle, J. Geophys. Res., Volume 102 (1997), pp. 5207-5229

[11] P.E. Janney; P.R. Castillo Geochemistry of the oldest Atlantic oceanic crust suggests mantle plume involvement in the early history of the central Atlantic Ocean, Earth Planet. Sci. Lett., Volume 192 (2001), pp. 291-302

[12] P. Machetel, E. Humler, High mantle temperature during Cretaceous, Earth Planet. Sci. Lett. (in press)

[13] P. Machetel; E. Thomassot Cretaceous length of day perturbation by mantle avalanche, Earth Planet. Sci. Lett., Volume 202 (2002), pp. 379-386

[14] P. Machetel; C. Thoraval; D. Brunet Spectral and geophysical consequences of 3D spherical mantle convection with an endothermic phase change at the 670-km discontinuity, Phys. Earth Planet. Inter., Volume 88 (1995), pp. 43-51

[15] P. Machetel; P. Weber Intermittent layered convection in a model mantle with an endothermic phase change at 670 km, Nature, Volume 350 (1991), pp. 55-57

[16] W.R. Peltier; L.P. Solheim Mantle phase transitions and layered chaotic convection, Geophys. Res. Lett., Volume 19 (1992), pp. 321-324

[17] M. Prévot; E. Mattern; P. Camps; M. Daignières Evidence for a 20° tilting of the Earth's rotation axis 110 million years ago, Earth Planet. Sci. Lett., Volume 179 (2000), pp. 517-528

[18] Y. Ricard; H.-C. Nataf; J.-P. Montagner The three-dimensional seismological model a-priori constrained: confrontation with seismic data, J. Geophys. Res., Volume 101 (1996), pp. 8457-8472

[19] K. Ricciardi; D. Abbott Increased mantle convection during the mid-Cretaceous: a comparative study of mantle potential temperature, J. Geophys. Res., Volume 101 (1996), pp. 8673-8684

[20] L.P. Solheim; W.R. Peltier Mantle phase transitions and layered convection, Can. J. Earth Sci., Volume 30 (1993), pp. 881-892

[21] V. Steinbach; D.A. Yuen Effects of depth-dependent properties on the thermal anomalies produced in flush instabilities from phase transitions, Phys. Earth Planet. Inter., Volume 86 (1994), pp. 165-183

[22] P.J. Tackley; D.J. Stevenson; G.A. Glatzmaier; G. Schubert Effects of an endothermic phase transition at 670-km depth in a spherical model of convection in the Earth's mantle, Nature, Volume 361 (1993)

[23] C. Thoraval; P. Machetel; A. Cazenave Locally layered convection inferred from dynamic models of the Earth's mantle, Nature, Volume 375 (1995), pp. 777-780

[24] R. van der Hilst; T. Seno Effects of relative plate motion on the deep structure and penetration depth of slabs below the Izu–Bonin ans Mariana Island arcs, Earth Planet. Sci. Lett., Volume 120 (1993), pp. 395-407

[25] S.A. Weinstein Catastrophic overturn of the Earth's mantle driven by multiple phase changes and internal heat generation, Geophys. Res. Lett., Volume 20 (1993), pp. 101-104

[26] G.E. Williams Geological constraints on the Precambrian history of Earth's rotation, Rev. Geophys., Volume 38 (2000), pp. 37-59

[27] D.A. Yuen; S. Balachandar; V.C. Steinbach; S. Honda; D.M. Reuteler; J.J. Smedsmo; G.S. Lauer Non-equilibrium effects of core-cooling and time dependent internal heating on mantle flush events, Nonlinear Processes in Geophysics, Volume 2 (1995), pp. 206-221

[28] D.A. Yuen; U. Hansen; W. Zhao; A.P. Vincent; A.V. Malevsky Hard turbulent thermal convection and thermal evolution of the mantle, J. Geophys. Res., Volume 98 (1993), pp. 5355-5373


Comments - Politique