@article{TutuSteinbergerSobolevetal.2018,
  author    = {Tutu, Anthony Osei and Steinberger, Bernhard and Sobolev, Stephan Vladimir and Rogozhina, Irina and Popov, Anton A.},
  title     = {Effects of upper mantle heterogeneities on the lithospheric stress field and dynamic topography},
  series = {Solid earth},
  volume    = {9},
  journal   = {Solid earth},
  number    = {3},
  publisher = {Copernicus},
  address   = {G{\"o}ttingen},
  issn      = {1869-9510},
  doi       = {10.5194/se-9-649-2018},
  pages     = {649 -- 668},
  year      = {2018},
  abstract  = {The orientation and tectonic regime of the observed crustal/lithospheric stress field contribute to our knowledge of different deformation processes occurring within the Earth's crust and lithosphere. In this study, we analyze the influence of the thermal and density structure of the upper mantle on the lithospheric stress field and topography. We use a 3-D lithosphere-asthenosphere numerical model with power-law rheology, coupled to a spectral mantle flow code at 300 km depth. Our results are validated against the World Stress Map 2016 (WSM2016) and the observation-based residual topography. We derive the upper mantle thermal structure from either a heat flow model combined with a seafloor age model (TM1) or a global S-wave velocity model (TM2). We show that lateral density heterogeneities in the upper 300 km have a limited influence on the modeled horizontal stress field as opposed to the resulting dynamic topography that appears more sensitive to such heterogeneities. The modeled stress field directions, using only the mantle heterogeneities below 300 km, are not perturbed much when the effects of lithosphere and crust above 300 km are added. In contrast, modeled stress magnitudes and dynamic topography are to a greater extent controlled by the upper mantle density structure. After correction for the chemical depletion of continents, the TM2 model leads to a much better fit with the observed residual topography giving a good correlation of 0.51 in continents, but this correction leads to no significant improvement of the fit between the WSM2016 and the resulting lithosphere stresses. In continental regions with abundant heat flow data, TM1 results in relatively small angular misfits. For example, in western Europe the misfit between the modeled and observation-based stress is 18.3°. Our findings emphasize that the relative contributions coming from shallow and deep mantle dynamic forces are quite different for the lithospheric stress field and dynamic topography.},
  language  = {en}
}
@article{TutuSobolevSteinbergeretal.2018,
  author    = {Tutu, Anthony Osei and Sobolev, Stephan Vladimir and Steinberger, Bernhard and Popov, A. A. and Rogozhina, Irina},
  title     = {Evaluating the Influence of Plate Boundary Friction and Mantle Viscosity on Plate Velocities},
  series = {Geochemistry, geophysics, geosystems},
  volume    = {19},
  journal   = {Geochemistry, geophysics, geosystems},
  number    = {3},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {1525-2027},
  doi       = {10.1002/2017GC007112},
  pages     = {642 -- 666},
  year      = {2018},
  abstract  = {Lithospheric plates move over the low-viscosity asthenosphere balancing several forces, which generate plate motions. We use a global 3-D lithosphere-asthenosphere model (SLIM3D) with visco-elasto-plastic rheology coupled to a spectral model of mantle flow at 300 km depth to quantify the influence of intra-plate friction and asthenospheric viscosity on plate velocities. We account for the brittle-ductile deformation at plate boundaries (yield stress) using a plate boundary friction coefficient to predict the present-day plate motion and net rotation of the lithospheric plates. Previous modeling studies have suggested that small friction coefficients (mu < 0.1, yield stress similar to 100 MPa) can lead to plate tectonics in models of mantle convection. Here we show that in order to match the observed present-day plate motion and net rotation, the frictional parameter must be less than 0.05. We obtain a good fit with the magnitude and orientation of the observed plate velocities (NUVEL-1A) in a no-net-rotation (NNR) reference frame with mu < 0.05 and a minimum asthenosphere viscosity of similar to 5 . 10(19) Pas to 10(20) Pas. Our estimates of net rotation (NR) of the lith-osphere suggest that amplitudes similar to 0.1-0.2 (degrees/Ma), similar to most observation-based estimates, can be obtained with asthenosphere viscosity cutoff values of similar to 10(19) Pas to 5 . 10(19) Pas and friction coefficients mu < 0.05.},
  language  = {en}
}