Multiscale Modeling of Glacial Loading by a 3D Thermo-Hydro-Mechanical Approach Including Erosion and Isostasy

We present a computational framework that allows investigating the Thermo-Hydro- Mechanical response of a representative part of a sedimentary basin during a glaciation cycle. We tackle the complexity of the problem, arising by the mutual interaction among several phenomena, by means of a multi-physics, multi-scale model with respect to both space and time. Our contribution addresses both the generation of the computational grid and the algorithm for the numerical solution of the problem. In particular we present a multi-scale approach accounting for the global deformation field of the lithosphere coupled with the Thermo-Hydro-Mechanical feedback of the ice load on a representative part of the domain at a finer scale. In the fine scale model we also include the erosion possibly caused by the ice melting. This methodology allows investigating the evolution of the sedimentary basin as a response to glaciation cycle at a fine scale, taking also into account the large spatial scale movement of the lithosphere due to isostasy. The numerical experiments are based on the analysis of simple scenario, and show the emergence of effects due to the multi-physics nature of the problem that are barely captured by simpler approaches.

Large fault slip peaking at trench in the 2011 Tohoku-oki earthquake

During the 2011 magnitude 9 Tohoku-oki earthquake, very large slip occurred on the shallowest part of the subduction megathrust. Quantitative information on the shallow slip is of critical importance to distinguishing between different rupture mechanics and understanding the generation of the ensuing devastating tsunami. However, the magnitude and distribution of the shallow slip are essentially unknown due primarily to the lack of near-trench constraints, as demonstrated by a compilation of 45 rupture models derived from a large range of data sets. To quantify the shallow slip, here we model high-resolution bathymetry differences before and after the earthquake across the trench axis. The slip is determined to be about 62 m over the most near-trench 40 km of the fault with a gentle increase towards the trench. This slip distribution indicates that dramatic net weakening or strengthening of the shallow fault did not occur during the Tohoku-oki earthquake.

The 2011 Tohoku-oki earthquake slip occurred on the shallowest part of the megathrust, but the nature of the shallow slip has been poorly constrained. Here, the authors model bathymetry differences before and after the earthquake to determine that the slip exceeded 60 m increasing towards the trench.

Rheology of the upper mantle: a synthesis

Rheological properties of the upper mantle of the Earth play an important role in the dynamics of the lithosphere and asthenosphere. However, such fundamental issues as the dominant mechanisms of flow have not been well resolved. A synthesis of laboratory studies and geophysical and geological observations shows that transitions between diffusion and dislocation creep likely occur in the Earth's upper mantle. The hot and shallow upper mantle flows by dislocation creep, whereas cold and shallow or deep upper mantle may flow by diffusion creep. When the stress increases, grain size is reduced and the upper mantle near the transition between these two regimes is weakened. Consequently, deformation is localized and the upper mantle is decoupled mechanically near these depths.