[1] Attempting to reconstruct the thermal history of the Earth from a geophysical point of view has for a long time been in disagreement with geochemical data. The geophysical approach uses parameterized models of mantle cooling. The rate of cooling of the Earth at the beginning of its history obtained in these models is generally too rapid to allow a sufficient present-day secular cooling rate. Geochemical estimates of radioactive element concentrations in the mantle then appear too low to explain the observed present mantle heat loss. Cooling models use scaling laws for the mean heat flux out of the mantle as a function of its Rayleigh number of the form Q / Ra b . Recent studies have introduced very low values of the exponent b, which can help reduce the cooling rate of the mantle. The present study instead focuses on the coefficient C in the relation Q = C Ra b and, in particular, on its variation with the wavelength of convection. The heat transfer strongly depends on the wavelength of convection. The length scale of convection in Earth's mantle is that of plate tectonics, implying convective cells of wide aspect ratio. Taking into account the long wavelength of convection in Earth's mantle can significantly reduce the efficiency of heat transfer. The likely variations of this wavelength with the Wilson cycle thus imply important variations of the heat flow out of the Earth on a intermediate timescale of 100 Ma, which renders parameterized models of thermal evolution inaccurate for quantitative predictions.

Abstractfor present time, continental growth must begin at least 3 Gy ago and steady-state for continental area must be reached for at least 1.5 Gy in our cooling model.

[1] A scaling law for the heat flux out of a convective fluid covered totally or partially by a finitely conducting lid is proposed. This scaling is constructed in order to quantify the heat transfer out of the Earth's mantle, taking into account the effect of the dichotomy between oceans and continents, which imposes heterogeneous thermal boundary conditions at the surface of the mantle. The effect of these heterogeneous boundary conditions is studied here using simple two-dimensional models, with the mantle represented by an isoviscous fluid heated from below and continents represented by nondeformable lids of finite thermal conductivity set above the surface of the model. We use free-slip boundary conditions under the oceanic and continental zones in order to study in an isolated way the possible thermal effect of continents, independently of all mechanical effect. A systematic study of the heat transfer as a function of the Rayleigh number of the fluid, of the width of the lid, and of its thermal properties is carried out. We show that estimates of continental lithosphere thickness imply a strong insulating effect from continents on mantle heat loss, at least locally. The heat flux below continents was low in the past and of the order of the present one if the continental thickness has remained broadly constant over the Earth's history.

No abstract

No abstract