Modelling Earth's surface topography: decomposition of the static and dynamic components

Mattia Guerri; Fabio Cammarano; Paul J. Tackley

doi:10.1016/j.pepi.2016.10.009

Modelling Earth's surface topography: decomposition of the static and dynamic components

Mattia Guerri, Fabio Cammarano, Paul J. Tackley

Geology

11 Citations (Scopus)

Abstract

Contrasting results on the magnitude of the dynamic component of topography motivate us to analyse the sources of uncertainties affecting long wavelength topography modelling. We obtain a range of mantle density structures from thermo-chemical interpretation of available seismic tomography models. We account for pressure, temperature and compositional effects as inferred by mineral physics to relate seismic velocity with density. Mantle density models are coupled to crustal density distributions obtained with a similar methodology. We compute isostatic topography and associated residual topography maps and perform instantaneous mantle flow modelling to calculate the dynamic topography. We explore the effects of proposed mantle 1-D viscosities and also test a 3D pressure- and temperature-dependent viscosity model. We find that the patterns of residual and dynamic topography are robust, with an average correlation coefficient (r) of respectively ∼0.74 and ∼0.71, upper-lower quartile ranges of 0.86–0.65 for residual topography and 0.83–0.62 for dynamic topography maps. The amplitudes are, on the contrary, poorly constrained. For the static component, the inferred density models of lithospheric mantle give an interquartile range of isostatic topography that is always higher than 100 m, reaching 1.7 km in some locations, and averaging ∼720 m. Crustal density models satisfying the same compressional velocity structure lead to variations in isostatic topography averaging 350 m, with peaks of 1 km in thick crustal regions, and always higher than 100 m. The uncertainties on isostatic topography are strong enough to mask, if present, the contribution of mantle convection to surface topography. For the dynamic component, we obtain a peak-to-peak dynamic topography amplitude exceeding 3 km for all our mantle density and viscosity models. These extremely high values would be associated with a magnitude of geoid undulations that is not in agreement with observations. Considering chemical heterogeneities in correspondence with the lower mantle Large Low Shear wave Velocity Provinces (LLSVPs) helps to decrease the peak-to-peak amplitudes of dynamic topography and geoid, but significantly reduces the correlation between synthetic and observed geoid. The correlation coefficients between all our residual and dynamic topography maps (a total of 220 and 198, respectively) is <0.55 (average = ∼0.19). The correlation slightly improves when considering only the very long-wavelength components of the maps (average = ∼0.23). We therefore conclude that a robust determination of dynamic topography is not feasible since current uncertainties affecting crustal density, mantle density and mantle viscosity are still too large. A truly interdisciplinary approach, combining constraints from the geological record with a multi-methodological interpretation of geophysical observations, is required to tackle the challenging task of linking the surface topography to deep processes.

Original language	English
Journal	Physics of the Earth and Planetary Interiors
Volume	261
Pages (from-to)	172–186
ISSN	0031-9201
DOIs	https://doi.org/10.1016/j.pepi.2016.10.009
Publication status	Published - 1 Dec 2016

Access to Document

10.1016/j.pepi.2016.10.009

Cite this

@article{a6ef6d1a46284db79eb39c97c2e5ee59,

title = "Modelling Earth's surface topography: decomposition of the static and dynamic components",

abstract = "Contrasting results on the magnitude of the dynamic component of topography motivate us to analyse the sources of uncertainties affecting long wavelength topography modelling. We obtain a range of mantle density structures from thermo-chemical interpretation of available seismic tomography models. We account for pressure, temperature and compositional effects as inferred by mineral physics to relate seismic velocity with density. Mantle density models are coupled to crustal density distributions obtained with a similar methodology. We compute isostatic topography and associated residual topography maps and perform instantaneous mantle flow modelling to calculate the dynamic topography. We explore the effects of proposed mantle 1-D viscosities and also test a 3D pressure- and temperature-dependent viscosity model. We find that the patterns of residual and dynamic topography are robust, with an average correlation coefficient (r) of respectively ∼0.74 and ∼0.71, upper-lower quartile ranges of 0.86–0.65 for residual topography and 0.83–0.62 for dynamic topography maps. The amplitudes are, on the contrary, poorly constrained. For the static component, the inferred density models of lithospheric mantle give an interquartile range of isostatic topography that is always higher than 100 m, reaching 1.7 km in some locations, and averaging ∼720 m. Crustal density models satisfying the same compressional velocity structure lead to variations in isostatic topography averaging 350 m, with peaks of 1 km in thick crustal regions, and always higher than 100 m. The uncertainties on isostatic topography are strong enough to mask, if present, the contribution of mantle convection to surface topography. For the dynamic component, we obtain a peak-to-peak dynamic topography amplitude exceeding 3 km for all our mantle density and viscosity models. These extremely high values would be associated with a magnitude of geoid undulations that is not in agreement with observations. Considering chemical heterogeneities in correspondence with the lower mantle Large Low Shear wave Velocity Provinces (LLSVPs) helps to decrease the peak-to-peak amplitudes of dynamic topography and geoid, but significantly reduces the correlation between synthetic and observed geoid. The correlation coefficients between all our residual and dynamic topography maps (a total of 220 and 198, respectively) is <0.55 (average = ∼0.19). The correlation slightly improves when considering only the very long-wavelength components of the maps (average = ∼0.23). We therefore conclude that a robust determination of dynamic topography is not feasible since current uncertainties affecting crustal density, mantle density and mantle viscosity are still too large. A truly interdisciplinary approach, combining constraints from the geological record with a multi-methodological interpretation of geophysical observations, is required to tackle the challenging task of linking the surface topography to deep processes.",

author = "Mattia Guerri and Fabio Cammarano and Tackley, {Paul J.}",

year = "2016",

month = dec,

day = "1",

doi = "10.1016/j.pepi.2016.10.009",

language = "English",

volume = "261",

pages = "172–186",

journal = "Physics of the Earth and Planetary Interiors",

issn = "0031-9201",

publisher = "Elsevier",

}

TY - JOUR

T1 - Modelling Earth's surface topography: decomposition of the static and dynamic components

AU - Guerri, Mattia

AU - Cammarano, Fabio

AU - Tackley, Paul J.

PY - 2016/12/1

Y1 - 2016/12/1

N2 - Contrasting results on the magnitude of the dynamic component of topography motivate us to analyse the sources of uncertainties affecting long wavelength topography modelling. We obtain a range of mantle density structures from thermo-chemical interpretation of available seismic tomography models. We account for pressure, temperature and compositional effects as inferred by mineral physics to relate seismic velocity with density. Mantle density models are coupled to crustal density distributions obtained with a similar methodology. We compute isostatic topography and associated residual topography maps and perform instantaneous mantle flow modelling to calculate the dynamic topography. We explore the effects of proposed mantle 1-D viscosities and also test a 3D pressure- and temperature-dependent viscosity model. We find that the patterns of residual and dynamic topography are robust, with an average correlation coefficient (r) of respectively ∼0.74 and ∼0.71, upper-lower quartile ranges of 0.86–0.65 for residual topography and 0.83–0.62 for dynamic topography maps. The amplitudes are, on the contrary, poorly constrained. For the static component, the inferred density models of lithospheric mantle give an interquartile range of isostatic topography that is always higher than 100 m, reaching 1.7 km in some locations, and averaging ∼720 m. Crustal density models satisfying the same compressional velocity structure lead to variations in isostatic topography averaging 350 m, with peaks of 1 km in thick crustal regions, and always higher than 100 m. The uncertainties on isostatic topography are strong enough to mask, if present, the contribution of mantle convection to surface topography. For the dynamic component, we obtain a peak-to-peak dynamic topography amplitude exceeding 3 km for all our mantle density and viscosity models. These extremely high values would be associated with a magnitude of geoid undulations that is not in agreement with observations. Considering chemical heterogeneities in correspondence with the lower mantle Large Low Shear wave Velocity Provinces (LLSVPs) helps to decrease the peak-to-peak amplitudes of dynamic topography and geoid, but significantly reduces the correlation between synthetic and observed geoid. The correlation coefficients between all our residual and dynamic topography maps (a total of 220 and 198, respectively) is <0.55 (average = ∼0.19). The correlation slightly improves when considering only the very long-wavelength components of the maps (average = ∼0.23). We therefore conclude that a robust determination of dynamic topography is not feasible since current uncertainties affecting crustal density, mantle density and mantle viscosity are still too large. A truly interdisciplinary approach, combining constraints from the geological record with a multi-methodological interpretation of geophysical observations, is required to tackle the challenging task of linking the surface topography to deep processes.

AB - Contrasting results on the magnitude of the dynamic component of topography motivate us to analyse the sources of uncertainties affecting long wavelength topography modelling. We obtain a range of mantle density structures from thermo-chemical interpretation of available seismic tomography models. We account for pressure, temperature and compositional effects as inferred by mineral physics to relate seismic velocity with density. Mantle density models are coupled to crustal density distributions obtained with a similar methodology. We compute isostatic topography and associated residual topography maps and perform instantaneous mantle flow modelling to calculate the dynamic topography. We explore the effects of proposed mantle 1-D viscosities and also test a 3D pressure- and temperature-dependent viscosity model. We find that the patterns of residual and dynamic topography are robust, with an average correlation coefficient (r) of respectively ∼0.74 and ∼0.71, upper-lower quartile ranges of 0.86–0.65 for residual topography and 0.83–0.62 for dynamic topography maps. The amplitudes are, on the contrary, poorly constrained. For the static component, the inferred density models of lithospheric mantle give an interquartile range of isostatic topography that is always higher than 100 m, reaching 1.7 km in some locations, and averaging ∼720 m. Crustal density models satisfying the same compressional velocity structure lead to variations in isostatic topography averaging 350 m, with peaks of 1 km in thick crustal regions, and always higher than 100 m. The uncertainties on isostatic topography are strong enough to mask, if present, the contribution of mantle convection to surface topography. For the dynamic component, we obtain a peak-to-peak dynamic topography amplitude exceeding 3 km for all our mantle density and viscosity models. These extremely high values would be associated with a magnitude of geoid undulations that is not in agreement with observations. Considering chemical heterogeneities in correspondence with the lower mantle Large Low Shear wave Velocity Provinces (LLSVPs) helps to decrease the peak-to-peak amplitudes of dynamic topography and geoid, but significantly reduces the correlation between synthetic and observed geoid. The correlation coefficients between all our residual and dynamic topography maps (a total of 220 and 198, respectively) is <0.55 (average = ∼0.19). The correlation slightly improves when considering only the very long-wavelength components of the maps (average = ∼0.23). We therefore conclude that a robust determination of dynamic topography is not feasible since current uncertainties affecting crustal density, mantle density and mantle viscosity are still too large. A truly interdisciplinary approach, combining constraints from the geological record with a multi-methodological interpretation of geophysical observations, is required to tackle the challenging task of linking the surface topography to deep processes.

U2 - 10.1016/j.pepi.2016.10.009

DO - 10.1016/j.pepi.2016.10.009

M3 - Journal article

SN - 0031-9201

VL - 261

SP - 172

EP - 186

JO - Physics of the Earth and Planetary Interiors

JF - Physics of the Earth and Planetary Interiors

ER -

Modelling Earth's surface topography: decomposition of the static and dynamic components

Abstract

Access to Document

Fingerprint

Cite this