the Creative Commons Attribution 4.0 License.

the Creative Commons Attribution 4.0 License.

# The impact of rheological uncertainty on dynamic topography predictions

### Patrice F. Rey

Much effort is being made to extract the dynamic components of the Earth's topography driven by density heterogeneities in the mantle. Seismically mapped density anomalies have been used as an input into mantle convection models to predict the present-day mantle flow and stresses applied on the Earth's surface, resulting in dynamic topography. However, mantle convection models give dynamic topography amplitudes generally larger by a factor of ∼2, depending on the flow wavelength, compared to dynamic topography amplitudes obtained by removing the isostatically compensated topography from the Earth's topography. In this paper, we use 3-D numerical experiments to evaluate the extent to which the dynamic topography depends on mantle rheology. We calculate the amplitude of instantaneous dynamic topography induced by the motion of a small spherical density anomaly (∼100 km radius) embedded into the mantle. Our experiments show that, at relatively short wavelengths (<1000 km), the amplitude of dynamic topography, in the case of non-Newtonian mantle rheology, is reduced by a factor of ∼2 compared to isoviscous rheology. This is explained by the formation of a low-viscosity channel beneath the lithosphere and a decrease in thickness of the mechanical lithosphere due to induced local reduction in viscosity. The latter is often neglected in global mantle convection models. Although our results are strictly valid for flow wavelengths less than 1000 km, we note that in non-Newtonian rheology all wavelengths are coupled, and the dynamic topography at long wavelengths will be influenced.

The Earth's mantle is continuously stirred by hot upwellings from the core–mantle boundary and by subduction of colder plates from the surface into the deep mantle (Pekeris, 1935; Isacks et al., 1968; Molnar and Tapponnier, 1975; Stern, 2002). This introduces temperature and density anomalies that stimulate mantle flow and forces dynamic uplift or subsidence at the plates' surfaces (Gurnis et al., 2000; Braun, 2010; Moucha and Forte, 2011; Flament et al., 2013). Dynamic topography can affect the entire planet's surface with varying magnitudes. Because it is typically a low-amplitude and long-wavelength transient signal, it is often dwarfed by isostatic topography associated with variations in the thickness and density of sediments, crust and mantle lithosphere.

For the present day, the observational constraints on dynamic topography come from residual topography measurements (Hoggard et al., 2016). Residual topography is calculated by removing the isostatically compensated topography from the Earth's topography (Crough, 1983; Cazenave et al., 1989; Davies and Pribac, 1993; Steinberger, 2007, 2016). The comprehensive work from Hoggard et al. (2016) revealed that residual topography varies between ±500 m at very-long wavelengths (i.e. ∼ 10 000 km) and can increase up to ±1000 m at shorter wavelengths (i.e. ∼1000 km). However, these residuals depend on our knowledge of the thermal and mechanical structure of the lithosphere and therefore may not be an accurate estimation of the deeper mantle contribution to the Earth's topography. Another approach to constrain present-day Earth's dynamic topography involves numerical modelling of present-day mantle flow using seismically mapped density anomalies as an input (Steinberger, 2007; Moucha et al., 2008; Conrad and Husson, 2009). However, this method requires a detailed knowledge of the viscosity structure in the Earth's interior (Parsons and Daly, 1983; Hager, 1984; Hager et al., 1985; Hager and Clayton, 1989), and translating seismic velocities to physical properties (e.g. temperature) of the mantle introduces further uncertainties (Cammarano et al., 2003). The problem is that dynamic topography predictions derived from mantle convection models are generally larger by a factor of 2 (more significant at the very-large scales) than estimates from residual topography (Hoggard et al., 2016; Cowie and Kusznir, 2018; Davies et al., 2019; Steinberger et al., 2019). We hypothesise that this could be related to an oversimplification of the mantle rheology. In this paper, we explore how, at wavelengths <1000 km, the magnitude of dynamic topography changes when we use a rheological model in which the viscosity depends on strain rate, temperature, pressure and fluid content. We first summarise the well-established analytical solution for calculating dynamic topography induced by a spherical density anomaly embedded into an isoviscous fluid (Morgan, 1965a; Molnar et al., 2015). Then, assuming isoviscous rheology, we illustrate that the amplitude of dynamic topography depends on the viscosity structure of the Earth's interior as shown by Morgan (1965a) and Molnar et al. (2015). Finally, we use 3-D coupled thermo–mechanical numerical experiments of the Stokes flow to assess the dependence of dynamic topography on nonlinear rheology using viscosity which depends on temperature, pressure, strain rate and fluid content. We show that plausible nonlinear rheologies can induce local variations in viscosity and result in dynamic topography of lower amplitude compared to those derived from models using isoviscous rheology.

## 2.1 Analytical solution for one layer isoviscous fluid

We assume here a simple 2-D model representing a very-viscous spherical
density anomaly embedded into a semi-infinite isoviscous fluid bounded by an
upper free surface. Earliest analytical investigations revealed that, albeit
counter-intuitive, the magnitude of the induced surface deflection due to
the rising sphere is independent of the viscosity of the fluid. The dynamic
topography is a function of the vertical total stress (*σ*_{zz})
applied to the surface which is proportional to the size and depth of the
density anomaly according to Eq. (1) (Morgan, 1965a,
b):

where *g* is the gravitational acceleration, *δ**ρ* is density
difference between the anomaly and the ambient material, *r* is radius of the
sphere, and *D* is distance from the surface to the centre of the anomaly
(modified from Morgan, 1965a, see Fig. 1a). The dynamic topography *e* is
given by the following:

where Δ*ρ* is the density difference between the mantle and air
(or water assuming a sea-load when *e*<0)
(Morgan, 1965a; Houseman and Hegarty, 1987). In
Fig. 1a, we plot the dynamic topography induced by a sphere of 1 %
density anomaly, whose centre is at 372 km depth (*D*=372 km) below the free
surface. We calculate the vertical total stress and convert it to dynamic
topography by using Eq. (2) for different values of the radius of the
sphere. The amplitude of dynamic topography shows an accelerating increase
by cubic dependence on the radius of the spherical density anomaly (Fig. 1a, solid black line). For the same problem, Molnar et al. (2015) provided a
solution by considering a higher-order term, resulting in a slight difference
from the solution of Morgan (1965a) (see Appendix A3 in Molnar et al., 2015),
which allows the consideration of density anomalies of finite viscosity (*η*_{sphere}) (Eq. 3):

where $C=\sqrt{{D}^{\mathrm{2}}+{x}^{\mathrm{2}}}$, and $f=({\mathit{\eta}}_{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}+\frac{\mathrm{3}{\mathit{\eta}}_{\mathrm{sphere}}}{\mathrm{2}})/({\mathit{\eta}}_{\mathrm{1}}+{\mathit{\eta}}_{\mathrm{sphere}})$. One can find that *f*=1.5 if
the sphere is very viscous (*η*_{sphere}≫*η*_{1}), and *f*<1.5 for any other case. In Fig. 1a, we present two more plots of dynamic
topography where *f*=1.5 for hard sphere and *f*=1.25 for *η*_{sphere}=*η*_{1} by using Eqs. (2) and (3). Figure 1a shows that a
rising deformable sphere creates higher dynamic topography compared to a
very-viscous sphere. These show that the viscosity contrast between the
spherical anomaly and the surrounding material can affect the dynamic
topography. In the section that follows, we explore how dynamic topography
varies when there is layering in viscosity such as the presence of a strong
lithosphere above the convective mantle.

## 2.2 The impact of layered viscosity structure on dynamic topography

A more generalised solution has been put forward to accommodate the presence
of a stronger upper layer representing a lithosphere with viscosity *η*_{2} above a weaker layer with viscosity *η*_{1}, and with *η*_{1}<*η*_{2} representing the convective mantle (Fig. 1b). In
this case, Morgan (1965a) showed (Eq. 4) that the total normal stress
induced by the density anomaly is dependent on the mass anomaly per unit
length (*M*_{u}, for point sources integrated along a continuous line), the
depth of the centre of the sphere (*D*) and marginally on the ratio of the
viscosity of the convective mantle to the viscosity of the lithosphere
($R={\mathit{\eta}}_{\mathrm{1}}/{\mathit{\eta}}_{\mathrm{2}}$). The 2-layer problem is treated in Fourier
domain with the resulting total normal stress as below:

where

*C*_{h}=cosh (*n**d*), *S*_{h}=sinh (*n**d*) and *d* is the upper-layer
thickness (modified from Morgan, 1965a). Following Morgan (1965a), Fig. 1b
illustrates the relative importance of *R* as well as the ratio of the thickness
of the upper layer to the depth of the anomaly (*d*∕*D*). As long as the
lithosphere is more viscous than the asthenosphere, the vertical total
stress at the surface has a minor dependence on the viscosity of the
lithosphere (see solid lines with *R*=1 and *R*=0.01 in Fig. 1b). Figure 1b
also shows that the magnitude of dynamic topography increases as the density
anomaly is brought closer to the surface (compare *R*=1, the solid black
line and the dashed black line). Moreover, its sensitivity on the relative
viscosity of the lithosphere also increases. Although an unrealistic
proposition for the Earth, when the lithosphere is less viscous than the
asthenosphere, the normal stress is much reduced and is strongly dependent
on the viscosity of the lithosphere (Fig. 1b). These demonstrate that
layering in viscosity can have a strong impact on the amplitude of dynamic
topography (Sembroni et al., 2017). In the next
section, we use the analytical solutions above to benchmark a numerical
model, which we will then extend to nonlinear viscosity.

## 2.3 Numerical solutions

For comparison with analytical solutions (Morgan, 1965a;
Molnar et al., 2015), we consider 3-D numerical models involving 1, 2 and 3
isoviscous layers. These benchmark experiments will be used as references
for non-isoviscous models discussed in Sect. 3. We use the open-source
code *Underworld* which solves the Stokes equation at insignificant Reynolds number
(Moresi et al., 2003, 2007). The 3-D
computational grid represents a domain 3840 km × 3840 km × 576 km with a
resolution of 6 km along the vertical *z* axis and 10 km along the *x* and *y* axes
(Fig. 2). In all experiments, we include a 42 km thick continental crust
above the upper mantle. The density structure is sensitive to the geotherm
via a coefficient of thermal expansion and compressibility (see Table 1 for
all parameters). The geotherm is defined using a radiogenic heat production
in the crust, a constant temperature of 20 ^{∘}C at the
surface and a constant temperature of 1350 ^{∘}C at 150 km.
We disregard the adiabatic heating, and the asthenosphere is kept at
1350 ^{∘}C. We embed a positive spherical temperature anomaly
of +324 ^{∘}C at a depth of 372 km below the surface, which
delivers a 1 % volumetric density difference. The radius of the sphere is
96 km. In all experiments, we impose free slip velocity boundary conditions
at all walls, such as *V*_{x} and *V*_{y} are set to be free, but *V*_{z}=0 cm yr^{−1} at the top wall. Taking advantage of the symmetry of the
experimental setup, we extract viscosity and velocity fields along a 2-D
cross section passing through the centre of the thermal anomaly, from which
we derive the streamlines and vertical velocity profiles along the vertical
axis at the centre of the models. We calculate the instantaneous dynamic
topography from the normal stress computed at the surface.

n/a – not applicable.
We use the rheological
parameters from ^{a} quartzite (Ranalli, 1995), ^{b} dry or wet
olivine (Hirth and Kohlstedt, 2003).

### 2.3.1 Dynamic topography due to a rising sphere in an isoviscous fluid

In the first experiment (Fig. 3a Experiment 1), we assign the same
depth-independent viscosity of 10^{21} Pa s to the crust, mantle and the
density anomaly. The streamlines for Experiment 1 (Fig. 3a) show formation
of two convective cells at the sides of the sphere covering the entire crust
and mantle. The vertical velocity profile indicates that the thermal anomaly
rises with a peak velocity of ∼2.4 cm yr^{−1}, which is
faster than the 2.0 cm yr^{−1} predicted by the analytical solution
(Fig. 4a). Experiment 1 predicts dynamic topography of 114 m (Fig. 4b) which
is lower than 132 m predicted by the analytical solution of Molnar et al. (2015).
We have verified that increasing the depth of our model from 576 to 864 km increases the dynamic topography from 114 to 122 m. Therefore, we
attribute the misfit in amplitude of dynamic topography to the finite space
in our numerical experiments. Our numerical experiment using isoviscous
material delivers a result globally consistent with the analytical solutions
of Morgan (1965a) and Molnar et al. (2015).

### 2.3.2 Dynamic topography on a strong lithosphere above an isoviscous asthenosphere

In Experiment 2, we assign to the lithosphere a constant viscosity 100 times
larger (10^{23} Pa s) than that of the asthenosphere (10^{21} Pa s, Fig. 3b) between *z*=150 km and base of the model. The convective cells become
narrower by the induced viscosity contrast (Fig. 3b). The streamlines are
deflected across the lithosphere–asthenosphere boundary due to the large
viscosity contrast (Fig. 3b), and there is a sharp variation in vertical
velocity at the base of the lithosphere (Fig. 4a, solid red line). The
maximum vertical velocity ∼2.1 cm yr^{−1} is attained near
the centre of the anomaly. When compared to Experiment 1, the dynamic
topography (Fig. 4b, solid red line) shows a significant increase from
∼114 to ∼174 m. This increase is consistent
with analytical estimations showing an increase in dynamic topography when
viscosity increases toward the surface (Fig. 1b, *R*<1). In
Experiment 2a (not shown here), we tested a different ratio of thickness of
the lithosphere to the depth of the anomaly (see *d*∕*D* in Eq. 4) by
increasing the lithospheric thickness from 150 to 200 km, while keeping
all parameters identical to those of Experiment 2. As predicted by Eq. (4),
Exp. 2 predicted dynamic topography of ∼191 m, being the
largest among all experiments (Fig. 4b, dashed red line). Overall, and
perhaps counter-intuitively, the presence of a thick viscous lithosphere
enhances the dynamic topography. Interestingly, in analogue experiments
where density anomaly is allowed to rise and interact with the lithosphere,
the amplitude of the dynamic topography is inversely correlated with the
thickness of the lithosphere (e.g. Griffiths et al.,
1989; Sembroni et al., 2017).

### 2.3.3 The impact of low-viscosity channel on the dynamic topography

In Experiment 3 (Fig. 3c), we introduce a third 60 km thick low-viscosity
layer (i.e. 10^{19} Pa s) beneath the base of the lithosphere. The
existence of a low-viscosity layer has been discussed in several studies
(e.g. Craig and McKenzie, 1986; Phipps Morgan et al., 1995; Stixrude and
Lithgow-Bertelloni, 2005; Becker, 2017). In this experiment, in order to
prevent large viscosity contrast that can impede the numerical convergence,
the viscosities of the lithosphere and that of the asthenosphere are set to
10^{22} and 10^{21} Pa s, respectively. When compared to Experiment
1, streamlines indicate a further decrease in size of the convective cells
and more importantly, a strong horizontal divergence of the streamlines
within the low-viscosity layer (Fig. 3c). The vertical velocities are also
enhanced in the asthenosphere reaching up to ∼ 2.8 cm yr^{−1} slightly above the centre of the anomaly (Fig. 4a, solid orange
line). When compared to Experiment 1, we observe a strong reduction in
dynamic topography (Fig. 4b, solid orange line) from 114 to 88 m. This is
due to the damping effect of the low-viscosity channel that acts as a
decoupling layer, which reduces the deviatoric stress through its ability to
flow.

Until now, the viscosities were assumed to be constant. However, results from experimental deformation on mantle rocks strongly suggest that the viscosity is highly nonlinear (Hirth and Kohlstedt, 2003). In what follows, we explore the influence of more realistic viscosities on dynamic topography.

## 3.1 Viscosity structure of the Earth's interior

The Earth's mantle is not isoviscous. Geological records of relative sea level changes related to postglacial rebound, geophysical observations of density anomalies inferred from seismic velocity variations in the mantle and satellite measurements of the longest wavelength components of the Earth's geoid have been used to infer the radial viscosity profile of the Earth's interior (Hager et al., 1985; Forte and Mitrovica, 1996; Mitrovica and Forte, 1997; Kaufmann and Lambeck, 2000). Henceforward, beneath the lithosphere, a variation in viscosity up to 2 orders of magnitude has been proposed (e.g. Kaufmann and Lambeck, 2000). Investigations of the rheological properties of crustal and mantle rocks via rock deformation experiments revealed a nonlinear dependence of viscosity on applied deviatoric stress, pressure, temperature, grain size and the presence of fluids (Post and Griggs, 1973; Chopra and Paterson, 1984; Karato, 1992; Karato and Wu, 1993; Gleason and Tullis, 1995; Ranalli, 1995; Hirth and Kohlstedt, 2003; Korenaga and Karato, 2008). These experiments lead to the following relationship:

where $\dot{\mathit{\epsilon}}$ and *A* stands for strain rate and pre-exponential
factor; *r* and *n* are exponents for water fugacity (${f}_{{\mathrm{H}}_{\mathrm{2}}\mathrm{O}}$) and deviatoric
stress, respectively; *V* and *Q* are the volume and energy of activation.

In the case where mantle flow is driven by the temperature difference at the
boundary of the convective layer or by internal heating, the dominant strain
mechanism is diffusion creep because low deviatoric stresses are expected in
the weak convective mantle (Karato and Wu, 1993;
Turcotte and Schubert, 2014). However, mantle flow in the vicinity of a
moving density anomaly is likely driven by deviatoric stresses that exceed
the threshold for dislocation creep. In this case, nonlinear viscosities
lead to strong local variation in viscosity. Are those local variations in
viscosity important for dynamic topography? To answer this question, we need
reasonable constraints on the rheological parameters controlling the
viscosity of mantle rocks. However, the extrapolation from laboratory strain
rates typically in the range of 10^{−6} s^{−1} to 10^{−4} s^{−1} to
mantle conditions where strain rates are typically on the order of
10^{−13} s^{−1} results in significant uncertainties on the activation
volume, activation energy and stress exponent (Hirth and
Kohlstedt, 2003; Korenaga and Karato, 2008). In what follows, we explore how
nonlinear viscosity impacts the dynamic topography and address how the
uncertainties on the activation volume can affect dynamic topography.

In Experiments 4 and 5 (Fig. 5), the viscosity depends on temperature, pressure and strain rate as indicated by Eq. (5), using published visco–plastic rheological parameters for the crust and mantle. Specifically, we use quartzite rheology for the crust (Ranalli, 1995), and we test both dry and wet olivine rheologies for the mantle (Hirth and Kohlstedt, 2003). Other parameters are identical to those in Experiments 1–3. We give all the rheological and thermal parameters in Table 1. For a given olivine rheology (i.e. dry or wet) we vary the activation volume by using the minimum and maximum reported values (Hirth and Kohlstedt, 2003).

In the numerical models, the plastic (i.e. brittle) deformation is described via

where *τ* is the 2nd invariant of the deviatoric stress tensor,
which varies with the coefficient of friction (*μ*), and depth via
lithostatic pressure (*σ*_{n}), as well as the cohesion (*C*_{0}).
Due to strain weakening, the cohesion and coefficient of friction decrease
from *C*_{0}=10 MPa and *μ*_{0}=0.577 to *C*_{0}=2 MPa and
*μ*_{1}=0.017 at which the maximum plastic strain (*ϵ*_{max}) is reached (i.e. 0.2, Table 1). The effective density (*ρ*) of
rocks is determined by the pressure and temperature using the following
equation:

where *ρ*_{0}, *T*_{0}, *α* and *P*_{0} signify the reference
density, reference temperature, thermal expansion coefficient and the
compressibility, respectively.

## 3.2 Numerical results: the case of dry olivine

In Experiments 4a and 4b, we consider dry dislocation creep for olivine (*n*>1, *p*=0, *r*=0). The reported activation volume for this rheology varies
between $\mathrm{6}\times {\mathrm{10}}^{-\mathrm{6}}$ and $\mathrm{27}\times {\mathrm{10}}^{-\mathrm{6}}$ m^{3} mol^{−1}
(Hirth and Kohlstedt, 2003). In Experiment 4a (Fig. 4b), we test
the lower value. The streamlines show a similar a pattern to Experiment 2.
Interestingly, the maximum vertical velocity peaks at 75 cm yr^{−1}, near
the upper boundary of the sphere (Fig. 6a, dashed black line). This is due
to the formation of a low-viscosity region above the rising sphere (Fig. 5a,
Experiment 4a). This experiment gives a dynamic topography of
∼149 m (Fig. 6b, dashed black line). It confirms that a
strong contrast in viscosity between the lithosphere and asthenosphere
enhances the dynamic topography signal. We note that the viscosity contrast
is attained by smoother transition between the lithosphere and asthenosphere
(Fig. 7a, dashed black line). We infer the mechanical thickness of the
lithosphere from the viscosity profiles plotted in Fig. 7a, along which
the lithosphere–asthenosphere transition zone shows a rapid decrease in
viscosity (Conrad and Molnar, 1997). We observe that the effective
mechanical thickness of the lithosphere is reduced to 140 km, compared to
the thickness of the thermal lithosphere (Fig. 7c).

When we increase the activation volume to $\mathrm{27}\times {\mathrm{10}}^{-\mathrm{6}}$ m^{3} mol^{−1},
the convection cells grow much larger and show continuity through the
lithosphere (Fig. 5a, Experiment 4b). The sphere has a very-low rising speed
of ∼ 0.25 cm yr^{−1} (Fig. 6a, solid black line). Compared
to Experiment 4a, the dynamic topography shows a strong decrease from
∼149 to ∼105 m (Fig. 6b, solid black line).
This is an example where the system behaves nearly as a single layer with
homogenous viscosity. The near absence of viscosity contrast between the
lithosphere and asthenosphere explains the smaller magnitude of the dynamic
topography. Moreover, the formation of moderately low-viscosity channel
(Fig. 7a, solid black line) also contributes to the decrease of the dynamic
topography.

## 3.3 Numerical results: the case of wet olivine

In Experiments 5a and 5b, we consider dislocation creep of wet olivine. The
reported activation volume varies between $\mathrm{11}\times {\mathrm{10}}^{-\mathrm{6}}$ and $\mathrm{33}\times {\mathrm{10}}^{-\mathrm{6}}$ m^{3} mol^{−1} (Hirth and Kohlstedt, 2003).
In Experiment 5a, we test the lower value. The streamlines show a pattern
similar to Experiment 4a but with slightly larger convective cells (Fig. 5b, Experiment 5a). The rising velocity of the anomaly exceeds 140 cm yr^{−1} (Fig. 6a, dashed orange line), promoted by the low-viscosity
region sitting above the rising anomaly. The dynamic topography is
∼110 m (Fig. 6b, dashed orange line). This is a bit
surprising given the strong contrast in viscosity (3 orders of magnitude)
between the lithosphere and asthenosphere. However, Fig. 7a shows that the
thickness of the mechanical lithosphere is reduced by about 30 km in
comparison to Experiment 4a (e.g. 10 km reduction from thermal thickness)
which resulted in lower dynamic topography with similar viscosity contrast
(Fig. 7b, c).

In Experiment 5b, we increase the activation volume from $\mathrm{11}\times {\mathrm{10}}^{-\mathrm{6}}$ to $\mathrm{33}\times {\mathrm{10}}^{-\mathrm{6}}$ m^{3} mol^{−1}. The vertical
velocities show significant decrease from 140 to 0.34 cm yr^{−1} (Fig. 6a, solid orange line). This is due to an increase in
viscosity above the rising sphere. Compared to Experiment 5a, the dynamic
topography decreases from ∼110 to ∼90 m
(Fig. 6b, solid orange line). Compared to Experiment 4b, we expect the
dynamic topography to be higher due to slight increase in viscosity contrast
(Fig. 7a, b). However, the increase in thickness of the low-viscosity channel
(Fig. 7a, d) is more effective and thereby causes a greater reduction in
magnitude of the dynamic topography.

In summary, experiments using nonlinear rheology generally give lower amplitudes of dynamic topography compared to experiments using isoviscous rheology (Fig. 8). When we use dry olivine rheology for the upper mantle, the dynamic topography varies between ∼105 and ∼149 m, whereas under wet conditions, the dynamic topography varies between ∼90 and ∼110 m (Fig. 8). These variations are due to uncertainties in the activation volume as well as fluid content in olivine rheologies.

Using coupled 3-D thermo-mechanical numerical experiments, we have modelled
the dynamic topography driven by a rising sphere of 1 % density anomaly,
having 96 km radius and emplaced at 372 km depth. In line with analytical
studies (Morgan, 1965a; Molnar et al., 2015), the
experiments show that dynamic topography is sensitive to viscosity contrast
between the lithosphere and asthenospheric mantle, and the thickness of the
lithosphere (Fig. 7). Higher viscosity contrasts amplify the dynamic
topography (Fig. 7a, b), whereas formation of a low-viscosity channel just
below the lithosphere has the opposite effect (Fig. 7a, d). The experiments
using nonlinear rheologies show local variations in viscosity, which
contribute to the dynamic thinning of the mechanical lithosphere and causes
reduction in dynamic topography. In addition, models using high-activation
volume creates a low-viscosity channel above the density anomaly, which
contributes to decreasing the dynamic topography. Using a larger viscosity
range in the models (10^{18} Pa s ≤ $\mathit{\eta}\left(PT\dot{\mathit{\epsilon}}\right)\le {\mathrm{10}}^{\mathrm{23}}$ Pa s) resulted in ∼5 % variation in the
amplitude of dynamic topography, indicating that the effects of nonlinear
rheology are reasonably captured in our models with smaller viscosity range
(10^{19} Pa s ≤ $\mathit{\eta}\left(PT\dot{\mathit{\epsilon}}\right)\le {\mathrm{10}}^{\mathrm{22}}$ Pa s).

Predictions of dynamic topography derived from mantle convection models are
compared against residual topography which is the component of Earth's
topography that is not compensated by isostasy (Flament et al.,
2013; Hoggard et al., 2016). In a recent work (Cowie and
Kusznir, 2018), it has been argued that dynamic topography predictions
require scaling of amplitudes by ∼0.75 to match the residual
topography, and when density anomalies shallower than 220 km are included,
the misfit requires a scaling factor of ∼0.35. It is also
important to consider that this misfit depends on the flow wavelength and is
suggested to be highest at lowest spherical harmonic degrees (*l*=2) or very-long wavelengths (Steinberger, 2016). Our numerical experiments
show that amplitude of dynamic topography can be nearly halved (e.g. from
∼174 m in Exp. 2 to ∼90 m in Exp. 5b) when we
consider nonlinear mantle rheology. Therefore, we propose that, at shorter
wavelengths (i.e. less than 1000 km), part of the misfit between the
dynamic topography extracted from mantle convection models and dynamic
topography estimated from residual topography can be attributed to the
Newtonian mantle viscosity used in convection models. If the density sources
are shallower, the dynamic topography becomes more sensitive to the
viscosity and density structure (Morgan, 1965a; Hager
and Clayton, 1989; Osei Tutu et al., 2018), and Newtonian viscosity may lead
to higher misfits.

Our models suggest that for shallow density anomalies in the mantle, nonlinear rheologies not only produce lateral variations in viscosity (Richards and Hager, 1989; Moucha et al., 2007) but also additional vertical variations in viscosity that impacts a relatively large area compared to the size of the anomaly in the mantle. We show that this impacts on the thickness of the mechanical lithosphere and predictions of the amplitude of dynamic topography.

As shown in Fig. 8, uncertainties on the activation volume result in variations in dynamic topography which are higher in experiments using dry olivine rheology (i.e. 17 %) compared to experiments using wet olivine rheology (10 %). The comparison between numerical experiments using dry olivine (Exp. 4a) and wet olivine (Exp. 5b) indicates that the variation in dynamic topography can be as much as 25 %. These variations can be lessened if we have better constraints on the mantle rheology, which will advance the dynamic topography models as well as our understanding of the interaction between deep mantle and the Earth's surface.

In our experiments we used Underworld, a free open-source code developed under the Australian AuScope initiative.

The version of Underworld code we used in our study can be found at: https://github.com/OlympusMonds/EarthByte_Underworld (last access: 25 October 2019).

To follow an open-source philosophy and promote reproducible science, we provide our input scripts (a suite of XML files) on the GitHub and EarthByte's freely accessible server: https://github.com/ofbodur/Bodur_and_Rey_EGU_SE_2019_Files (last access: 25 November 2019), https://www.earthbyte.org/ (last access: 19 December 2019).

ÖFB designed the experiments and wrote the paper. PFR contributed to the analysis of numerical modelling results and improved the paper.

The authors declare that they have no conflict of interest.

This article is part of the special issue “Understanding the unknowns: the impact of uncertainty in the geosciences”. It is not associated with a conference.

This research was undertaken with the assistance of resources from the National Computational Infrastructure (NCI) through the National Computational Merit Allocation Scheme supported by the Australian Government, the Pawsey Supercomputing Centre with funding from the Australian Government and the government of Western Australia, and support from the Australian Research Council through the Industrial Transformation Research Hub grant ARC-IH130200012. We benefited from discussions with Gregory Houseman and Nicolas Flament about the Earth's dynamic topography.

This research has been supported by the Australian Research Council (grant no. IH130200012).

This paper was edited by Lucia Perez-Diaz and reviewed by Bernhard Steinberger and Mark Hoggard.

Becker, T. W.: Superweak asthenosphere in light of upper mantle seismic anisotropy, Geochem. Geophy. Geosy., 18, 1986–2003, https://doi.org/10.1002/2017GC006886, 2017.

Braun, J.: The many surface expressions of mantle dynamics, Nat. Geosci., 3, 825–833, https://doi.org/10.1038/ngeo1020, 2010.

Cammarano, F., Goes, S., Vacher, P., and Giardini, D.: Inferring upper-mantle temperatures from seismic velocities, Phys. Earth Planet. In., 138, 197–222, https://doi.org/10.1016/S0031-9201(03)00156-0, 2003.

Cazenave, A., Souriau, A., and Dominh, K.: Global coupling of Earth surface topography with hotspots, geoid and mantle heterogeneities, Nature, 340, 54–57, https://doi.org/10.1038/340054a0, 1989.

Chopra, P. N. and Paterson, M. S.: The role of water in the deformation of dunite, J. Geophys. Res., 89, 7861–7876, 1984.

Conrad, C. P. and Husson, L.: Influence of dynamic topography on sea level and its rate of changeInfluence of dynamic topography on sea level, Lithosphere, 1, 110–120, https://doi.org/10.1130/L32.1, 2009.

Conrad, C. P. and Molnar, P.: The growth of Rayleigh-Taylor-type instabilities in the lithosphere for various rheological and density structures, Geophys. J. Int., 129, 95–112, 1997.

Cowie, L. and Kusznir, N.: Renormalisation of global mantle dynamic topography predictions using residual topography measurements for “normal” oceanic crust, Earth Planet. Sc. Lett., 499, 145–156, https://doi.org/10.1016/j.epsl.2018.07.018, 2018.

Craig, C. H. and McKenzie, D.: The existence of a thin low-viscosity layer beneath the lithosphere, Earth Planet. Sc. Lett., 78, 420–426, 1986.

Crough, S. T.: Hotspot swells, Annu. Rev. Earth Pl. Sc., 11, 165–193, 1983.

Davies, D. R., Valentine, A. P., Kramer, S. C., Rawlinson, N., Hoggard, M. J., Eakin, C. M., and Wilson, C. R.: Earth's multi-scale topographic response to global mantle flow, Nat. Geosci., 12, 845–850, 2019.

Davies, G. F. and Pribac, F.: Mesozoic seafloor subsidence and the Darwin Rise, past and present, Washingt. DC, Am. Geophys. Union Geophys. Monogr. Ser., 77, 39–52, 1993.

Flament, N., Gurnis, M., and Müller, R. D.: A review of observations and models of dynamic topography, Lithosphere, 5, 189–210, https://doi.org/10.1130/L245.1, 2013.

Forte, A. M. and Mitrovica, J. X.: New inferences of mantle viscosity from joint inversion of long-wavelength mantle convection and post-glacial rebound data, Geophys. Res. Lett., 23, 1147–1150, https://doi.org/10.1029/96GL00964, 1996.

Gleason, G. C. and Tullis, J.: A flow law for dislocation creep of quartz aggregates determined with the molten salt cell, Tectonophysics, 247, 1–23, https://doi.org/10.1016/0040-1951(95)00011-B, 1995.

Griffiths, R. W., Gurnis, M., and Eitelberg, G.: Holographic measurements of surface topography in laboratory models of mantle hotspots, Geophys. J. Int., 96, 477–495, 1989.

Gurnis, M., Mitrovica, J. X., Ritsema, J., and Van Heijst, H. J.: Constraining mantle density structure using geological evidence of surface uplift rates: The case of the African Superplume, Geochem. Geophy. Geosy., 1, 1020, https://doi.org/10.1029/1999GC000035, 2000.

Hager, B. H.: Constraints on Mantle Rheology and Flow, J. Geophys. Res., 89, 6003–6015, 1984.

Hager, B. H. and Clayton, R. W.: Constraints on the Structure of Mantle Convection Using Seismic Observations, Flow Models, and the Geoid, in: Mantle Convection: Plate Tectonics and Global Dynamics. Fluid Mechanics of Astrophysics and Geophysics, No. 4, Gordon and Breach Science Publishers, New York, 657–763, ISBN 9780677221205, 1989.

Hager, B., Clayton, R., Richards, Hager, B. H., Clayton, R. W., Richards, M. A., Comer, R. P., and Dziewonski, A. M.: Lower mantle heterogeneity, dynamic topography and the geoid, Nature, 313, 541–545, https://doi.org/10.1038/313541a0, 1985.

Hirth, G. and Kohlstedt, D.: Rheology of the upper mantle and the mantle wedge: A view from the experimentalists, Insid. Subduction Fact., 138, 83–105, 2003.

Hoggard, M. J., White, N., and Al-Attar, D.: Global dynamic topography observations reveal limited influence of large-scale mantle flow, Nat. Geosci., 9, 456–463, https://doi.org/10.1038/ngeo2709, 2016.

Houseman, G. A. and Hegarty, K. A.: Did rifting on Australia's Southern Margin result from tectonic uplift?, Tectonics, 6, 515–527, https://doi.org/10.1029/TC006i004p00515, 1987.

Isacks, B., Oliver, J., and Sykes, L. R.: Seismology and the new global tectonics, J. Geophys. Res., 73, 5855–5899, 1968.

Karato, S.: On the Lehmann discontinuity, Geophys. Res. Lett., 19, 2255–2258, 1992.

Karato, S. and Wu, P.: Rheology of the upper mantle: a synthesis, Science, 260, 771–778, https://doi.org/10.1126/science.260.5109.771, 1993.

Kaufmann, G. and Lambeck, K.: Mantle dynamics, postglacial rebound and the radial viscosity profile, Phys. Earth Planet. In., 121, 301–324, https://doi.org/10.1016/S0031-9201(00)00174-6, 2000.

Korenaga, J. and Karato, S.: A new analysis of experimental data on olivine rheology, J. Geophys. Res.-Sol. Ea., 113, B02403, https://doi.org/10.1029/2007JB005100, 2008.

Mitrovica, J. X. and Forte, A. M.: Radial profile of mantle viscosity: Results from the joint inversion of convection and postglacial rebound observables, J. Geophys. Res., 102, 2751, https://doi.org/10.1029/96JB03175, 1997.

Molnar, P. and Tapponnier, P.: Cenozoic tectonics of Asia: effects of a continental collision, Science, 189, 419–426, 1975.

Molnar, P., England, P. C., and Jones, C. H.: Mantle dynamics, isostasy, and the support of high terrain, J. Geophys. Res.-Sol. Ea., 120, 1932–1957, https://doi.org/10.1002/2014JB011724, 2015.

Moresi, L., Dufour, F., and Mühlhaus, H.-B.: A Lagrangian integration point finite element method for large deformation modeling of viscoelastic geomaterials, J. Comput. Phys., 184, 476–497, https://doi.org/10.1016/S0021-9991(02)00031-1, 2003.

Moresi, L., Quenette, S., Lemiale, V., Mériaux, C., Appelbe, B., and Mühlhaus, H.-B.: Computational approaches to studying non-linear dynamics of the crust and mantle, Phys. Earth Planet. In., 163, 69–82, https://doi.org/10.1016/j.pepi.2007.06.009, 2007.

Morgan, W. J.: Gravity anomalies and convection currents: 1. A sphere and cylinder sinking beneath the surface of a viscous fluid, J. Geophys. Res., 70, 6175–6187, 1965a.

Morgan, W. J.: Gravity anomalies and convection currents: 2. The Puerto Rico Trench and the Mid-Atlantic Rise, J. Geophys. Res., 70, 6189–6204, 1965b.

Moucha, R. and Forte, A. M.: Changes in African topography driven by mantle convection, Nat. Geosci., 4, 707–712, https://doi.org/10.1038/ngeo1235, 2011.

Moucha, R., Forte, A. M., Mitrovica, J. X., and Daradich, A.: Lateral variations in mantle rheology: Implications for convection related surface observables and inferred viscosity models, Geophys. J. Int., 169, 113–135, https://doi.org/10.1111/j.1365-246X.2006.03225.x, 2007.

Moucha, R., Forte, A. M., Mitrovica, J. X., Rowley, D. B., Quéré, S., Simmons, N. A., and Grand, S. P.: Dynamic topography and long-term sea-level variations: There is no such thing as a stable continental platform, Earth Planet. Sc. Lett., 271, 101–108, 2008.

Osei Tutu, A., Steinberger, B., Sobolev, S. V., Rogozhina, I., and Popov, A. A.: Effects of upper mantle heterogeneities on the lithospheric stress field and dynamic topography, Solid Earth, 9, 649–668, https://doi.org/10.5194/se-9-649-2018, 2018.

Parsons, B. and Daly, S.: The relationship between surface topography, gravity anomalies and temperature structure of convection, J. Geophys. Res., 88, 1129–1144, https://doi.org/10.1029/JB088iB02p01129, 1983.

Pekeris, C. L.: Thermal convection in the interior of the Earth, Geophys. Suppl. Mon. Not. R. Astron. Soc., 3, 343–367, 1935.

Phipps Morgan, J., Morgan, W. J., Zhang, Y.-S., and Smith, W. H. F.: Observational hints for a plume-fed, suboceanic asthenosphere and its role in mantle convection, J. Geophys. Res.-Sol. Ea., 100, 12753–12767, https://doi.org/10.1029/95JB00041, 1995.

Post, R. L. and Griggs, D. T.: The earth's mantle: evidence of non-Newtonian flow, Science, 181, 1242–1244, 1973.

Ranalli, G.: Rheology of the Earth, 2nd Edn., Chapman & Hall, London, 1995.

Richards, M. A. and Hager, B. H.: Effects of lateral viscosity variations on long-wavelength geoid anomalies and topography, J. Geophys. Res.-Sol. Ea., 94, 10299–10313, 1989.

Sembroni, A., Kiraly, A., Faccenna, C., Funiciello, F., Becker, T. W., Globig, J., and Fernandez, M.: Impact of the lithosphere on dynamic topography: Insights from analogue modeling, Geophys. Res. Lett., 44, 2693–2702, https://doi.org/10.1002/2017GL072668, 2017.

Steinberger, B.: Effects of latent heat release at phase boundaries on flow in the Earth's mantle, phase boundary topography and dynamic topography at the Earth's surface, Phys. Earth Planet. In., 164, 2–20, 2007.

Steinberger, B.: Topography caused by mantle density variations: observation-based estimates and models derived from tomography and lithosphere thickness, Geophys. Suppl. to Mon. Not. R. Astron. Soc., 205, 604–621, 2016.

Steinberger, B., Conrad, C. P., Osei Tutu, A., and Hoggard, M. J.: On the amplitude of dynamic topography at spherical harmonic degree two, Tectonophysics, 760, 221–228, https://doi.org/10.1016/j.tecto.2017.11.032, 2019.

Stern, R. J.: Subduction zones, Rev. Geophys., 40, 1031–1039, https://doi.org/10.1029/2001RG000108, 2002.

Stixrude, L. and Lithgow-Bertelloni, C.: Mineralogy and elasticity of the oceanic upper mantle: Origin of the low-velocity zone, J. Geophys. Res.-Sol. Ea., 110, 1–16, https://doi.org/10.1029/2004JB002965, 2005.

Turcotte, D. and Schubert, G.: Geodynamics, Cambridge university press, 340–341, 2014.

- Abstract
- Introduction
- Dynamic topography driven by a rising sphere: analytical and numerical solutions
- The impact of nonlinear viscosity on dynamic topography
- Discussion and conclusion
- Code and data availability
- Author contributions
- Competing interests
- Special issue statement
- Acknowledgements
- Financial support
- Review statement
- References

- Abstract
- Introduction
- Dynamic topography driven by a rising sphere: analytical and numerical solutions
- The impact of nonlinear viscosity on dynamic topography
- Discussion and conclusion
- Code and data availability
- Author contributions
- Competing interests
- Special issue statement
- Acknowledgements
- Financial support
- Review statement
- References