Henri-Claude Nataf
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Research Interests
Henri-Claude Nataf has a total of 94 co-authors. Here are the top co-authors by number of shared publications:
- Stéphane Labrosse
ENS de Lyon, Université Lyon-1, LGL-TPE , 46 allée d’Italie, F-69007 Lyon
15 shared publications - Thomas Frasson
Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, IRD, Univ. Gustave Eiffel, ISTerre, 38000 Grenoble, France
15 shared publications Nicolas Coltice
Université Côte d'Azur, CNRS, Observatoire de la Côte d'Azur, IRD, Géoazur, France
13 shared publications- Nicolas Flament
Environmental Futures, School of Science, University of Wollongong, Wollongong, Australia
11 shared publications - Nathanaël Schaeffer
Université Grenoble Alpes
10 shared publications - Philippe Cardin
ISTerre, Université de Grenoble 1, CNRS, Grenoble, France
10 shared publications
Latest Publications
- Open Access
- 8 May 2026
- Published by Elsevier in Earth and Planetary Science Letters
- Vol. 687
Effect of realistic heat flux patterns on geodynamo simulations
- 8 May 2026
- Published by Elsevier in Earth and Planetary Science Letters
- Vol. 687
Open AccessEffect of realistic heat flux patterns on geodynamo simulations
- Open Access
- 6 November 2025
- Published by Springer Nature in Experiments in Fluids
- Vol. 66 (12) , 1
Correction to: Modal acoustic velocimetry in libration-driven flows
- 6 November 2025
- Published by Springer Nature in Experiments in Fluids
- Vol. 66 (12) , 1
Open AccessCorrection to: Modal acoustic velocimetry in libration-driven flows
- 15 September 2025
- Published by Springer Nature in Experiments in Fluids
- Vol. 66 (10) , 1-25
Modal acoustic velocimetry in libration-driven flows
- 15 September 2025
- Published by Springer Nature in Experiments in Fluids
- Vol. 66 (10) , 1-25
Modal acoustic velocimetry in libration-driven flows
- Open Access
- 28 December 2024
- Published by Oxford University Press (OUP) in Geophysical Journal International
- Vol. 240 (3) , 1481-1504
Geomagnetic dipole stability and zonal flows controlled by mantle heat flux heterogeneities
SUMMARY: Palaeomagnetic evidence shows that the behaviour of the geodynamo has changed during geological times. These changes are visible through variations in the strength and stability of the magnetic dipole. Variations in the heat flux at the core–mantle boundary (CMB) due to mantle convection have been suggested as one possible mechanism capable of driving such a change of behaviour. This work aims at acquiring a more complete understanding of how lateral heterogeneities of the CMB heat flux affect the geodynamo while other relevant parameters are pushed towards realistic values. For this purpose, we ran geodynamo simulations with degree 1 and 2 spherical harmonic patterns of heat flux at the CMB. Several geodynamo models are used, ranging from standard numerical dynamos to more extreme parameters, including strong field cases and turbulent cases. We show that heat flux heterogeneities with amplitudes compatible with our knowledge of mantle convection history can favour multipolar dynamos. The multipolar transition is associated with a disruption of westward flows either through eastward thermal winds or through a loss of equatorial symmetry. Strong field dynamo models are found to have larger westward flows and are less sensitive to heat flux heterogeneities. Furthermore, we find that the dipolar fraction of the magnetic field correlates with $M_{Za}^{*}=\dfrac{\Lambda _{Za}}{Rm_{Za}^2}$ where $\Lambda _{Za}$ is the zonal antisymmetric Elsasser number and $Rm_{Za}$ is the zonal antisymmetric magnetic Reynolds number. Importantly, $M_{Za}^{*}$ estimated for the Earth’s core is consistent with a reversing dipolar magnetic field. Within the range of $M_{Za}^{*}$ susceptible to reversals, breaking the equatorial symmetry or forcing eastward zonal flows through an equatorial cooling of the core consistently triggers reversals or a transition towards multipolar dynamos in our simulations. Our results support that time variations of heat-flux heterogeneities driven by mantle convection through Earth’s history are capable of inducing the significant variations in the reversal frequency observed in the palaeomagnetic record.- 28 December 2024
- Published by Oxford University Press (OUP) in Geophysical Journal International
- Vol. 240 (3) , 1481-1504
Open AccessGeomagnetic dipole stability and zonal flows controlled by mantle heat flux heterogeneities
SUMMARY: Palaeomagnetic evidence shows that the behaviour of the geodynamo has changed during geological times. These changes are visible through variations in the strength and stability of the magnetic dipole. Variations in the heat flux at the core–mantle boundary (CMB) due to mantle convection have been suggested as one possible mechanism capable of driving such a change of behaviour. This work aims at acquiring a more complete understanding of how lateral heterogeneities of the CMB heat flux affect the geodynamo while other relevant parameters are pushed towards realistic values. For this purpose, we ran geodynamo simulations with degree 1 and 2 spherical harmonic patterns of heat flux at the CMB. Several geodynamo models are used, ranging from standard numerical dynamos to more extreme parameters, including strong field cases and turbulent cases. We show that heat flux heterogeneities with amplitudes compatible with our knowledge of mantle convection history can favour multipolar dynamos. The multipolar transition is associated with a disruption of westward flows either through eastward thermal winds or through a loss of equatorial symmetry. Strong field dynamo models are found to have larger westward flows and are less sensitive to heat flux heterogeneities. Furthermore, we find that the dipolar fraction of the magnetic field correlates with $M_{Za}^{*}=\dfrac{\Lambda _{Za}}{Rm_{Za}^2}$ where $\Lambda _{Za}$ is the zonal antisymmetric Elsasser number and $Rm_{Za}$ is the zonal antisymmetric magnetic Reynolds number. Importantly, $M_{Za}^{*}$ estimated for the Earth’s core is consistent with a reversing dipolar magnetic field. Within the range of $M_{Za}^{*}$ susceptible to reversals, breaking the equatorial symmetry or forcing eastward zonal flows through an equatorial cooling of the core consistently triggers reversals or a transition towards multipolar dynamos in our simulations. Our results support that time variations of heat-flux heterogeneities driven by mantle convection through Earth’s history are capable of inducing the significant variations in the reversal frequency observed in the palaeomagnetic record. - Open Access
- 16 September 2024
- Published by Copernicus GmbH
Potential links between internal and external geodynamics at periods around 6 years
The GRACEFUL project aims to better understand dynamic processes in the Earth’s deep interior using a combination of satellite observations of the Earth’s magnetic field, gravity field, and rotation. Significant oscillations at periods ranging from 6 to 8 years have been linked to dynamical processes in the fluid outer core and at the core-mantle boundary, influencing the three aforementioned observables. However, the occurrence of a 6-year cycle was also recently evidenced in the climate system, expressed in the global mean surface temperature, zonal winds, precipitation, terrestrial water storage changes, sea level changes, and ice mass changes. This study provides a comprehensive review of the 6-year cycle observed across the entire Earth system, from its deep interior to its fluid external envelopes. We propose several mechanisms that could explain small amplitude variations in geodetic observations, such as the length of day and gravity field, highlighting potential links between internal and external geodynamics at periods of around 6 years.
- 16 September 2024
- Published by Copernicus GmbH
Open AccessPotential links between internal and external geodynamics at periods around 6 years
The GRACEFUL project aims to better understand dynamic processes in the Earth’s deep interior using a combination of satellite observations of the Earth’s magnetic field, gravity field, and rotation. Significant oscillations at periods ranging from 6 to 8 years have been linked to dynamical processes in the fluid outer core and at the core-mantle boundary, influencing the three aforementioned observables. However, the occurrence of a 6-year cycle was also recently evidenced in the climate system, expressed in the global mean surface temperature, zonal winds, precipitation, terrestrial water storage changes, sea level changes, and ice mass changes. This study provides a comprehensive review of the 6-year cycle observed across the entire Earth system, from its deep interior to its fluid external envelopes. We propose several mechanisms that could explain small amplitude variations in geodetic observations, such as the length of day and gravity field, highlighting potential links between internal and external geodynamics at periods of around 6 years.
- 21 May 2024
- Published by MathDoc/Centre Mersenne in Comptes Rendus. Géoscience
- Vol. 356 (G1) , 1-30
Dynamic regimes in planetary cores: τ – ℓ diagrams
Planetary cores are the seat of rich and complex fluid dynamics, in which the effects of rotation and magnetic field combine. The equilibria governing the strength of the magnetic field produced by the dynamo effect, the organisation and amplitude of the flow, and those of the density field, remain debated despite remarkable progress made in their numerical simulation. This paper describes an approach based on the explicit consideration of the variation of time scales τ with spatial scales ℓ for the different physical phenomena involved. The τ – ℓ diagrams thus constructed constitute a very complete graphic summary of the dynamics of the object under study. We highlight the role of the available convective power in controlling this dynamics, together with the relevant force balance, for which we derive a very telling τ – ℓ translation. Several scenarios are constructed and discussed for the Earth’s core, shedding new light on the width of convective columns, and on the force equilibria to be considered. A QG-MAC scenario adapted from [Aubert, 2019] gives a good account of the observations. A diversion to Venus reveals the subtlety and relativity of the notion of “fast rotator”. We discuss scaling laws and their validity domain, and illustrate “path strategies”. A complete toolbox is provided, allowing everyone to construct a τ – ℓ diagram of a numerical simulation, a laboratory experiment, a theory, or a natural object.
- 21 May 2024
- Published by MathDoc/Centre Mersenne in Comptes Rendus. Géoscience
- Vol. 356 (G1) , 1-30
Dynamic regimes in planetary cores: τ – ℓ diagrams
Planetary cores are the seat of rich and complex fluid dynamics, in which the effects of rotation and magnetic field combine. The equilibria governing the strength of the magnetic field produced by the dynamo effect, the organisation and amplitude of the flow, and those of the density field, remain debated despite remarkable progress made in their numerical simulation. This paper describes an approach based on the explicit consideration of the variation of time scales τ with spatial scales ℓ for the different physical phenomena involved. The τ – ℓ diagrams thus constructed constitute a very complete graphic summary of the dynamics of the object under study. We highlight the role of the available convective power in controlling this dynamics, together with the relevant force balance, for which we derive a very telling τ – ℓ translation. Several scenarios are constructed and discussed for the Earth’s core, shedding new light on the width of convective columns, and on the force equilibria to be considered. A QG-MAC scenario adapted from [Aubert, 2019] gives a good account of the observations. A diversion to Venus reveals the subtlety and relativity of the notion of “fast rotator”. We discuss scaling laws and their validity domain, and illustrate “path strategies”. A complete toolbox is provided, allowing everyone to construct a τ – ℓ diagram of a numerical simulation, a laboratory experiment, a theory, or a natural object.
- Open Access
- 14 May 2024
- Published by Copernicus GmbH in Solid Earth
- Vol. 15 (5) , 617-637
On the impact of true polar wander on heat flux patterns at the core–mantle boundary
The heat flux across the core–mantle boundary (CMB) is a fundamental variable for Earth evolution and internal dynamics. Seismic tomography provides access to seismic heterogeneities in the lower mantle, which can be related to present-day thermal heterogeneities. Alternatively, mantle convection models can be used to either infer past CMB heat flux or to produce statistically realistic CMB heat flux patterns in self-consistent models. Mantle dynamics modifies the inertia tensor of the Earth, which implies a rotation of the Earth with respect to its spin axis, a phenomenon called true polar wander (TPW). This rotation must be taken into account to link the dynamics of the mantle to the dynamics of the core. In this study, we explore the impact of TPW on the CMB heat flux over long timescales (∼1 Gyr) using two recently published mantle convection models: one model driven by a plate reconstruction and a second that self-consistently produces a plate-like behaviour. We compute the geoid in both models to correct for TPW. In the plate-driven model, we compute a total geoid and a geoid in which lateral variations of viscosity and density are suppressed above 350 km depth. An alternative to TPW correction is used for the plate-driven model by simply repositioning the model in the original paleomagnetic reference frame of the plate reconstruction. The average TPW rates range between 0.4 and 1.8° Myr−1, but peaks up to 10° Myr−1 are observed. We find that in the plate-driven mantle convection model used in this study, the maximum inertia axis produced by the model does not show a long-term consistency with the position of the magnetic dipole inferred from paleomagnetism. TPW plays an important role in redistributing the CMB heat flux, notably at short timescales (≤10 Myr). Those rapid variations modify the latitudinal distribution of the CMB heat flux, which is known to affect the stability of the magnetic dipole in geodynamo simulations. A principal component analysis (PCA) is computed to obtain the dominant CMB heat flux pattern in the different cases. These heat flux patterns are representative of the mantle convection cases studied here and can be used as boundary conditions for geodynamo models.
- 14 May 2024
- Published by Copernicus GmbH in Solid Earth
- Vol. 15 (5) , 617-637
Open AccessOn the impact of true polar wander on heat flux patterns at the core–mantle boundary
The heat flux across the core–mantle boundary (CMB) is a fundamental variable for Earth evolution and internal dynamics. Seismic tomography provides access to seismic heterogeneities in the lower mantle, which can be related to present-day thermal heterogeneities. Alternatively, mantle convection models can be used to either infer past CMB heat flux or to produce statistically realistic CMB heat flux patterns in self-consistent models. Mantle dynamics modifies the inertia tensor of the Earth, which implies a rotation of the Earth with respect to its spin axis, a phenomenon called true polar wander (TPW). This rotation must be taken into account to link the dynamics of the mantle to the dynamics of the core. In this study, we explore the impact of TPW on the CMB heat flux over long timescales (∼1 Gyr) using two recently published mantle convection models: one model driven by a plate reconstruction and a second that self-consistently produces a plate-like behaviour. We compute the geoid in both models to correct for TPW. In the plate-driven model, we compute a total geoid and a geoid in which lateral variations of viscosity and density are suppressed above 350 km depth. An alternative to TPW correction is used for the plate-driven model by simply repositioning the model in the original paleomagnetic reference frame of the plate reconstruction. The average TPW rates range between 0.4 and 1.8° Myr−1, but peaks up to 10° Myr−1 are observed. We find that in the plate-driven mantle convection model used in this study, the maximum inertia axis produced by the model does not show a long-term consistency with the position of the magnetic dipole inferred from paleomagnetism. TPW plays an important role in redistributing the CMB heat flux, notably at short timescales (≤10 Myr). Those rapid variations modify the latitudinal distribution of the CMB heat flux, which is known to affect the stability of the magnetic dipole in geodynamo simulations. A principal component analysis (PCA) is computed to obtain the dominant CMB heat flux pattern in the different cases. These heat flux patterns are representative of the mantle convection cases studied here and can be used as boundary conditions for geodynamo models.
- Open Access
- 21 December 2023
- Published by Copernicus GmbH
- Vol. 2023, 1-24
Reply on AC5
Heat flux across the core-mantle boundary (CMB) is an important variable of Earth's thermal evolution and dynamics. Seismic tomography provides access to seismic heterogeneities in the lower mantle, which can be related to present-day thermal heterogeneities. Alternatively, mantle convection models can be used to either infer past CMB heat flux or to produce statistically realistic CMB heat flux patterns in self-consistent models. Mantle dynamics modifies the inertia tensor of the Earth, which implies a rotation of the Earth with respect to its spin axis, a phenomenon called true polar wander (TPW). This rotation must be taken into account to link the dynamics of the mantle to the dynamics of the core. In this study, we use two recently published mantle convection models to explore the impact of TPW on the CMB heat flux over long timescales (~ 1 Gyr). One of the mantle convection models is driven by a plate reconstruction, while the other self-consistently produces a plate-like behavior. We compute the geoid in both models to correct for TPW. In the plate-driven model, we compute a total geoid and a geoid in which lateral variations of viscosity and temperature are suppressed above 350 km depth. We show that TPW plays an important role in redistributing the CMB heat flux, notably at short time scales (≤ 10 Myr). Those rapid variations modify the latitudinal distribution of the CMB heat flux, which is known to affect the stability of the magnetic dipole in geodynamo simulations. A principal component analysis (PCA) is computed to obtain the dominant CMB heat flux pattern in the different cases. These heat flux patterns can be used as boundary conditions for geodynamo models as representative of the mantle convection cases studied here. We note that the geoids produced by the two models are widely different from each other and from the observed present-day geoid. Work thus still needs to be done to improve the computation of the geoid in mantle convection models related to plate-tectonics.
- 21 December 2023
- Published by Copernicus GmbH
- Vol. 2023, 1-24
Open AccessReply on AC5
Heat flux across the core-mantle boundary (CMB) is an important variable of Earth's thermal evolution and dynamics. Seismic tomography provides access to seismic heterogeneities in the lower mantle, which can be related to present-day thermal heterogeneities. Alternatively, mantle convection models can be used to either infer past CMB heat flux or to produce statistically realistic CMB heat flux patterns in self-consistent models. Mantle dynamics modifies the inertia tensor of the Earth, which implies a rotation of the Earth with respect to its spin axis, a phenomenon called true polar wander (TPW). This rotation must be taken into account to link the dynamics of the mantle to the dynamics of the core. In this study, we use two recently published mantle convection models to explore the impact of TPW on the CMB heat flux over long timescales (~ 1 Gyr). One of the mantle convection models is driven by a plate reconstruction, while the other self-consistently produces a plate-like behavior. We compute the geoid in both models to correct for TPW. In the plate-driven model, we compute a total geoid and a geoid in which lateral variations of viscosity and temperature are suppressed above 350 km depth. We show that TPW plays an important role in redistributing the CMB heat flux, notably at short time scales (≤ 10 Myr). Those rapid variations modify the latitudinal distribution of the CMB heat flux, which is known to affect the stability of the magnetic dipole in geodynamo simulations. A principal component analysis (PCA) is computed to obtain the dominant CMB heat flux pattern in the different cases. These heat flux patterns can be used as boundary conditions for geodynamo models as representative of the mantle convection cases studied here. We note that the geoids produced by the two models are widely different from each other and from the observed present-day geoid. Work thus still needs to be done to improve the computation of the geoid in mantle convection models related to plate-tectonics.
- Open Access
- 24 November 2023
- Published by Copernicus GmbH
- Vol. 2023, 1-24
Reply on RC3
Heat flux across the core-mantle boundary (CMB) is an important variable of Earth's thermal evolution and dynamics. Seismic tomography provides access to seismic heterogeneities in the lower mantle, which can be related to present-day thermal heterogeneities. Alternatively, mantle convection models can be used to either infer past CMB heat flux or to produce statistically realistic CMB heat flux patterns in self-consistent models. Mantle dynamics modifies the inertia tensor of the Earth, which implies a rotation of the Earth with respect to its spin axis, a phenomenon called true polar wander (TPW). This rotation must be taken into account to link the dynamics of the mantle to the dynamics of the core. In this study, we use two recently published mantle convection models to explore the impact of TPW on the CMB heat flux over long timescales (~ 1 Gyr). One of the mantle convection models is driven by a plate reconstruction, while the other self-consistently produces a plate-like behavior. We compute the geoid in both models to correct for TPW. In the plate-driven model, we compute a total geoid and a geoid in which lateral variations of viscosity and temperature are suppressed above 350 km depth. We show that TPW plays an important role in redistributing the CMB heat flux, notably at short time scales (≤ 10 Myr). Those rapid variations modify the latitudinal distribution of the CMB heat flux, which is known to affect the stability of the magnetic dipole in geodynamo simulations. A principal component analysis (PCA) is computed to obtain the dominant CMB heat flux pattern in the different cases. These heat flux patterns can be used as boundary conditions for geodynamo models as representative of the mantle convection cases studied here. We note that the geoids produced by the two models are widely different from each other and from the observed present-day geoid. Work thus still needs to be done to improve the computation of the geoid in mantle convection models related to plate-tectonics.
- 24 November 2023
- Published by Copernicus GmbH
- Vol. 2023, 1-24
Open AccessReply on RC3
Heat flux across the core-mantle boundary (CMB) is an important variable of Earth's thermal evolution and dynamics. Seismic tomography provides access to seismic heterogeneities in the lower mantle, which can be related to present-day thermal heterogeneities. Alternatively, mantle convection models can be used to either infer past CMB heat flux or to produce statistically realistic CMB heat flux patterns in self-consistent models. Mantle dynamics modifies the inertia tensor of the Earth, which implies a rotation of the Earth with respect to its spin axis, a phenomenon called true polar wander (TPW). This rotation must be taken into account to link the dynamics of the mantle to the dynamics of the core. In this study, we use two recently published mantle convection models to explore the impact of TPW on the CMB heat flux over long timescales (~ 1 Gyr). One of the mantle convection models is driven by a plate reconstruction, while the other self-consistently produces a plate-like behavior. We compute the geoid in both models to correct for TPW. In the plate-driven model, we compute a total geoid and a geoid in which lateral variations of viscosity and temperature are suppressed above 350 km depth. We show that TPW plays an important role in redistributing the CMB heat flux, notably at short time scales (≤ 10 Myr). Those rapid variations modify the latitudinal distribution of the CMB heat flux, which is known to affect the stability of the magnetic dipole in geodynamo simulations. A principal component analysis (PCA) is computed to obtain the dominant CMB heat flux pattern in the different cases. These heat flux patterns can be used as boundary conditions for geodynamo models as representative of the mantle convection cases studied here. We note that the geoids produced by the two models are widely different from each other and from the observed present-day geoid. Work thus still needs to be done to improve the computation of the geoid in mantle convection models related to plate-tectonics.
- Open Access
- 24 November 2023
- Published by Copernicus GmbH
- Vol. 2023, 1-24
Reply on EC2
Heat flux across the core-mantle boundary (CMB) is an important variable of Earth's thermal evolution and dynamics. Seismic tomography provides access to seismic heterogeneities in the lower mantle, which can be related to present-day thermal heterogeneities. Alternatively, mantle convection models can be used to either infer past CMB heat flux or to produce statistically realistic CMB heat flux patterns in self-consistent models. Mantle dynamics modifies the inertia tensor of the Earth, which implies a rotation of the Earth with respect to its spin axis, a phenomenon called true polar wander (TPW). This rotation must be taken into account to link the dynamics of the mantle to the dynamics of the core. In this study, we use two recently published mantle convection models to explore the impact of TPW on the CMB heat flux over long timescales (~ 1 Gyr). One of the mantle convection models is driven by a plate reconstruction, while the other self-consistently produces a plate-like behavior. We compute the geoid in both models to correct for TPW. In the plate-driven model, we compute a total geoid and a geoid in which lateral variations of viscosity and temperature are suppressed above 350 km depth. We show that TPW plays an important role in redistributing the CMB heat flux, notably at short time scales (≤ 10 Myr). Those rapid variations modify the latitudinal distribution of the CMB heat flux, which is known to affect the stability of the magnetic dipole in geodynamo simulations. A principal component analysis (PCA) is computed to obtain the dominant CMB heat flux pattern in the different cases. These heat flux patterns can be used as boundary conditions for geodynamo models as representative of the mantle convection cases studied here. We note that the geoids produced by the two models are widely different from each other and from the observed present-day geoid. Work thus still needs to be done to improve the computation of the geoid in mantle convection models related to plate-tectonics.
- 24 November 2023
- Published by Copernicus GmbH
- Vol. 2023, 1-24
Open AccessReply on EC2
Heat flux across the core-mantle boundary (CMB) is an important variable of Earth's thermal evolution and dynamics. Seismic tomography provides access to seismic heterogeneities in the lower mantle, which can be related to present-day thermal heterogeneities. Alternatively, mantle convection models can be used to either infer past CMB heat flux or to produce statistically realistic CMB heat flux patterns in self-consistent models. Mantle dynamics modifies the inertia tensor of the Earth, which implies a rotation of the Earth with respect to its spin axis, a phenomenon called true polar wander (TPW). This rotation must be taken into account to link the dynamics of the mantle to the dynamics of the core. In this study, we use two recently published mantle convection models to explore the impact of TPW on the CMB heat flux over long timescales (~ 1 Gyr). One of the mantle convection models is driven by a plate reconstruction, while the other self-consistently produces a plate-like behavior. We compute the geoid in both models to correct for TPW. In the plate-driven model, we compute a total geoid and a geoid in which lateral variations of viscosity and temperature are suppressed above 350 km depth. We show that TPW plays an important role in redistributing the CMB heat flux, notably at short time scales (≤ 10 Myr). Those rapid variations modify the latitudinal distribution of the CMB heat flux, which is known to affect the stability of the magnetic dipole in geodynamo simulations. A principal component analysis (PCA) is computed to obtain the dominant CMB heat flux pattern in the different cases. These heat flux patterns can be used as boundary conditions for geodynamo models as representative of the mantle convection cases studied here. We note that the geoids produced by the two models are widely different from each other and from the observed present-day geoid. Work thus still needs to be done to improve the computation of the geoid in mantle convection models related to plate-tectonics.