Kinematics and dynamics of the southeastern margin of the Tibetan Plateau

On the southeastern margin of the Tibetan Plateau lies a large region which seismicity and GPS data show to be actively deforming. This paper describes the active faulting in the region, and how it relates to the velocity field observed with GPS. In places the velocity field is accommodated by rotations about vertical axes, and most or all of the strain at the surface in the region appears to be released seismically. GPS velocities are then compared to velocities calculated using a model for deformation driven by gravitational driving forces. Using rheologies estimated from experimentally derived mineral flow laws, the model provides velocities that are in good agreement with observed GPS velocities. It is not possible to uniquely determine the rheology or flow velocity at depth, and there are two forms of model solution which match the observed horizontal surface velocities. In one of these, vertical planes deform by pure shear, and in the other vertical gradients of horizontal velocity are present within the crust. Two distinct regions of normal-faulting earthquakes are present in the region, and have mechanisms which are most easily explained by gravity-driven deformation.     

Post-seismic reloading and temporal clustering on a single fault

Geological studies show evidence for temporal clustering of large earthquakes on individual fault systems. Since post-seismic deformation due to the inelastic rheology of the lithosphere may result in a variable loading rate on a fault throughout the interseismic period, it is reasonable to expect that the rheology of the non-seismogenic lower crust and mantle lithosphere may play a role in controlling earthquake recurrence times. We study this phenomenon using a 2-D, finite element method continuum model of the lithosphere containing a single strike-slip fault. This model builds on a previous study using a 1-D spring-dashpot-slider analogue of a single fault system to study the role of Maxwell viscoelastic relaxation in producing non-periodic earthquakes. In our 2-D model, the seismogenic portion of the fault slips when a predetermined yield stress is exceeded; stress accumulated on the seismogenic fault is shed to the viscoelastic layers below and recycled back to the seismogenic fault through viscoelastic relaxation. We find that random variation of the fault yield stress from one earthquake to the next can cause the earthquake sequence to be clustered; the amount of clustering depends on a non-dimensional number, W , called the Wallace number defined as the standard deviation of the randomly varied fault yield stress divided by the effective viscosity of the system times the tectonic loading rate. A new clustering metric based on the bimodal distribution of interseismic intervals allows us to investigate clustering behaviour of systems over a wide range of model parameters and those with multiple viscoelastic layers. For models with   W ≥ 1  clustering increases with increasing W , while those with   W ≤ 1  are unclustered, or quasi-periodic.     

Effect of 3-D viscoelastic structure on post-seismic relaxation from the 2004 M= 9.2 Sumatra earthquake

The 2004 M = 9.2 Sumatra–Andaman earthquake profoundly altered the state of stress in a large volume surrounding the ∼1400 km long rupture. Induced mantle flow fields and coupled surface deformation are sensitive to the 3-D rheology structure. To predict the post-seismic motions from this earthquake, relaxation of a 3-D spherical viscoelastic earth model is simulated using the theory of coupled normal modes. The quasi-static deformation basis set and solution on the 3-D model is constructed using: a spherically stratified viscoelastic earth model with a linear stress–strain relation; an aspherical perturbation in viscoelastic structure; a 'static' mode basis set consisting of Earth's spheroidal and toroidal free oscillations; a "viscoelastic" mode basis set; and interaction kernels that describe the coupling among viscoelastic and static modes. Application to the 2004 Sumatra–Andaman earthquake illustrates the profound modification of the post-seismic flow field at depth by a slab structure and similarly large effects on the near-field post-seismic deformation field at Earth's surface. Comparison with post-seismic GPS observations illustrates the extent to which viscoelastic relaxation contributes to the regional post-seismic deformation.     

Comment on 'Consequences of progressive eclogitization on crustal exhumation, a mechanical study' by H. Raimbourg, L. Jolivet and Y. Leroy

Two simple end-member models of a subduction channel have been proposed in the literature: (i) the 'pressure-imposed' model for which the pressure within the channel is assumed to be lithostatic, the channel walls have negligible strength with respect to lateral pressure gradients, and the channel geometry therefore varies with time and (ii) the 'geometry-imposed' model of constant channel geometry, rigid walls and resultant lateral variation in pressure. Neither of these models is realistic, but they provide lower and upper bounds to potential pressure distributions in natural subduction zones. The critical parameter is the relative strength of the confining plates, reflected in the effective viscosity ratio between the channel fill and the walls. The assertion that the 'geometry-imposed' model is internally inconsistent is incorrect—it merely represents one bound to possible behaviour and a bound that may be approached for realistic values of the effective viscosity for weak channel fill (e.g. unconsolidated ocean-floor sediments) and relatively cold and strong subducting and overriding lithospheric plates.     

Statistical properties of the Transantarctic Mountains (TAM) micrometeorite collection

Micrometeorites have been recovered from traps located at the summit of nunataks in the Transantarctic Mountains (TAM), Antarctica. They constitute the TAM micrometeorite collection. Micrometeorites accumulated by direct infall for hundreds of thousands of years. This long collection duration is confirmed by the wide range of weathering by dissolution of olivine in the stony micrometeorites from the TAM collection. A statistical study of the size distribution and frequency by type of this collection, and comparison with other Antarctic micrometeorite collections (the South Pole Water Well collection and the Walcott Névé collection), suggest that the TAM collection is essentially unbiased. Thanks to the very long exposure of the traps, large diameter (>1000 μm) micrometeorites are present in sufficiently large numbers to allow a statistically meaningful estimate of their size distribution in this size range for the first time. We found that the slope of the size distribution remains constant in the 100–1600 μm size range. Therefore, the size distribution of micrometeorites in this size range is controlled by a single process.     

Slab pull effects from a flexural analysis of the Tonga and Kermadec trenches (Pacific Plate)

Thin-plate flexure models have been frequently used to explain the mechanical behaviour of the lithosphere at oceanic trenches, but little attention has been paid to using them as a way to check the relative importance of different plate-driving mechanisms. A 2-D numerical algorithm accounting for the flexural deflection of the lithosphere controlled by multilayered elastic–plastic rheology (brittle–elastic–ductile) has been applied to the seaward side of the Tonga and Kermadec trenches. This approach gives a better fit to the bathymetry on both trenches than assuming classical homogeneous plate models, and allows the interplate coupling forces and the lithospheric strength profile to be constrained. Our results show that, in order to fit the observed deflection of the lithosphere, a regional tensile horizontal force must act in both regions. This tensile force and its flexural effects are discussed in terms of slab pull as a main plate-driving mechanism. The predicted stress and yielding distributions partially match the outer-rise earthquake hypocentres within the subducting plate, and thus do not invalidate the model.     

The influence of a stratified rheology on the flexural response of the lithosphere to (un)loading by extensional faulting

We present a two-layered finite difference model for the flexural response of the lithosphere to extensional faulting. The model allows for three modes of flexure: (1) fully coupled, with the upper crust and mantle welded together by the lower crust; (2) fully decoupled, with the upper crust and mantle behaving as independent layers; and (3) partly decoupled, signifying that the response of the upper crust to small-wavelength loads is superimposed on the response of the entire lithosphere to long-wavelength loads. Which of these modes of flexure is to be expected depends on the rheology and especially the thermal state of the lithosphere. Coupled behaviour is related to a cold and strong lithosphere. The Baikal Rift Zone provides a typical example for this mode of flexure. A fully decoupled lithosphere is an exceptional case, related to anomalous high temperatures in the lower crust, and is observed in the Basin and Range province. The most common case is a partly decoupled lithosphere, with the degree of decoupling depending on the thickness and viscosity of the lower crust. This is inferred, for example, for the Bay of Biscay margin.     

Impact of regional mantle flow on subducting plate geometry and interplate stress: insights from physical modelling

Physical models of subduction investigate the impact of regional mantle flow on the structure of the subducted slab and deformation of the downgoing and overriding plates. The initial mantle flow direction beneath the overriding plate can be horizontal or vertical, depending on its location with respect to the asthenospheric flow field. Imposed mantle flow produces either over or underpressure on the lower surface of the slab depending on the initial mantle flow pattern (horizontal or vertical, respectively). Overpressure promotes shallow dip subduction while underpressure tends to steepen the slab. Horizontal mantle flow with rates of 1–10 cm yr⁻¹ provides sufficient overpressure on a dense subducting lithosphere to obtain a subduction angle of  ∼60°  , while the same lithospheric slab sinks vertically when no flow is imposed. Vertical drag force (due to downward mantle flow) exerted on a slab can result in steep subduction if the slab is neutrally buoyant but fails to produce steep subduction of buoyant oceanic lithosphere. The strain regime in the overriding plate due to the asthenospheric drag force depends largely on slab geometry. When the slab dip is steeper than the interplate zone, the drag force produces negative additional normal stress on the interplate zone and tensile horizontal stress in the overriding plate. When the slab dip is shallower than the interplate zone, an additional positive normal stress is produced on the interplate zone and the overriding plate experiences additional horizontal compressive stress. However, the impact of the mantle drag force on interplate pressure is small compared to the influence of the slab pull force since these stress variations can only be observed when the slab is dense and interplate pressure is low.     

Modelling deformation rates in the western Gulf of Corinth: rheological constraints

The Gulf of Corinth is one of the most active extensional regions in the Mediterranean area characterized by a high rate of seismicity. However, there are still open questions concerning the role and the geometry of the numerous active faults bordering the basin, as well as the mechanisms governing the seismicity. In this paper, we use a 2-D plane strain finite element analysis to constrain the upper crust rheology by modelling the available deformation data (GPS and geomorphology). We consider a SSW–NNE cross-section of the rift cutting the main active normal faults (Aigion, West Eliki and Off-Shore faults). The models run for 650 Kyr assuming an elasto-viscoplastic rheology and 1.3 cm yr⁻¹ horizontal extension as boundary condition (resulting from GPS data). We model the horizontal and vertical deformation rates and the accumulation of plastic strain at depth, and we compare them with GPS data, with long term uplift rates inferred from geomorphology and with the distribution of seismicity, respectively. Our modelling results demonstrate that dislocation on high-angle normal faults in a plastic crustal layer plays a key role in explaining the extremely localized strain within the Gulf of Corinth. Conversely, the contribution of structures such as the antithetic Trizonia fault or the buried hypothetical subhorizontal discontinuity are not necessary to model observed data.     

Post-seismic rebound of a spherical Earth: new insights from the application of the Post–Widder inversion formula

The post-seismic response of a viscoelastic Earth to a seismic dislocation can be computed analytically within the framework of normal-modes, based on the application of propagator methods. This technique, widely documented in the literature, suffers from several shortcomings; the main drawback is related to the numerical solution of the secular equation, whose degree increases linearly with the number of viscoelastic layers so that only coarse-layered models are practically solvable. Recently, a viable alternative to the standard normal-mode approach, based on the Post–Widder Laplace inversion formula, has been proposed in the realm of postglacial rebound models. The main advantage of this method is to bypass the explicit solution of the secular equation, while retaining the analytical structure of the propagator formalism. At the same time, the numerical computation is much simplified so that additional features such as linear non-Maxwell rheologies can be simply implemented. In this work, for the first time, we apply the Post–Widder Laplace inversion formula to a post-seismic rebound model. We test the method against the standard normal-mode solution and we perform various benchmarks aimed to tune the algorithm and to optimize computation performance while ensuring the stability of the solution. As an application, we address the issue of finding the minimum number of layers with distinct elastic properties needed to accurately describe the post-seismic relaxation of a realistic Earth model. Finally, we demonstrate the potentialities of our code by modelling the post-seismic relaxation after the 2004 Sumatra–Andaman earthquake comparing results based upon Maxwell and Burgers rheologies.