Gravity constraints on the dynamics of the crust–mantle system during Calabrian subduction

The thermomechanic evolution of the lithosphere–upper mantle system during Calabrian subduction is analysed using a 2-D finite element approach, in which the lithosphere is compositionally stratified into crust and mantle. Gravity and topography predictions are cross-checked with observed gravity and topography patterns of the Calabrian region. Modelling results indicate that the gravity pattern in the arc-trench region is shaped by the sinking of light material, belonging to both the overriding and subduction plates. The sinking of light crustal material, up to depths of the order of 100–150 km is the ultimate responsible for the peculiar gravity signature of subduction, characterized by a minimum of gravity anomaly located at the trench, bounded by two highs located on the overriding and subducting plates, with a variation in magnitude of the order of 200 mGal along a wavelength of 200 km, in agreement with the isostatically compensated component of gravity anomaly observed along a transect crossing the Calabrian Arc, from the Tyrrhenian to the Ionian Seas. The striking agreement between the geodetic retrieved profiles and the modelled ones in the trench region confirms the crucial role of compositional stratification of the lithosphere in the subduction process and the correctness of the kinematic hypotheses considered in our modelling, that the present-day configuration of crust–mantle system below the Calabrian arc results from trench's retreat at a rate of about 3 cm yr⁻¹, followed by gravitational sinking of the subducted slab in the last 5 Myr.     

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.     

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.     

Numerical modelling study on the flexural uplift of the Transantarctic Mountains

In this study, based on a 2-D thermomechanical finite element model, the uplift of the Transantarctic Mountains (TAM) is discussed in relation to the flexural uplift of a rheologically layered lithosphere, which is described by Vening-Meinesz's cantilever kinematics. The general model behaviour shows that the thickness of the crust and the geothermal gradient in the lithosphere are the principal factors in controlling the effective elastic thickness ( T _e). Although T _e is also significantly dependent on the magnitude of the uplift and the wet or dry rheological condition of rocks, these two factors do not have a dominant influence on the half-wavelength of the TAM. The model with a plausible crustal structure beneath Antarctica shows that the thermal structure beneath East Antarctica is the critical factor, controlling the half-wavelength of the TAM. If there is a significant radiogenic heat source in the Antarctic lithosphere, T _e beneath East Antarctica is estimated to be 50 km, at most, and the lithosphere has no potential to explain an exceptionally large-scale half-wavelength of the TAM. Even for the model without a heat source, if East Antarctica is significantly thermally influenced by West Antarctica, T _e is estimated to be about 60 km, and it is difficult to reproduce the half-wavelength of the TAM. Contrarily, when a radiogenic heat source is absent and the thermal structure beneath East Antarctica is not significantly affected by that beneath West Antarctica, the rheological structure beneath East Antarctica has the potential to reproduce the half-wavelength of the TAM ( T _e∼ 100 km). Thus, the presence of a radiogenic heat source in the crust and mantle and the thermal influence of West Antarctica on East Antarctica are crucial factors in the reproduction of the half-wavelength found in the TAM.     

A multidisciplinary study of a slow-slipping fault for seismic hazard assessment: the example of the Middle Durance Fault (SE France)

Assessing seismic hazard in continental interiors is difficult because these regions are characterized by low strain rates and may be struck by infrequent destructive earthquakes. In this paper, we provide an example showing that interpretations of seismic cross sections combined with other kinds of studies such as analysis of microseismicity allow the whole seismogenic source area to be imaged in this type of region. The Middle Durance Fault (MDF) is an 80-km-long fault system located southeastern France that has a moderate but regular seismicity and some palaeoseismic evidence for larger events. It behaves as an oblique ramp with a left-lateral-reverse fault slip and has a low strain rate. MDF is one of the rare slow active fault system monitored by a dedicated dense velocimetric short period network. This study showed a fault system segmented in map and cross section views which consists of staircase basement faults topped by listric faults ramping off Triassic evaporitic beds. Seismic sections allowed the construction of a 3-D structural model used for accurate location of microseismicity. Southern part of MDF is mainly active in the sedimentary cover. In its northern part and in Alpine foreland, seismicity deeper than 8 km was also recorded meaning active faults within the crust cannot be excluded. Seismogenic potential of MDF was roughly assessed. Resulting source sizes and estimated slip rates imply that the magnitude upper limit ranges from 6.0 to 6.5 with a return period of a few thousand years. The present study shows that the coupling between 3-D fault geometry imaging and accurate location of microseismicity provides a robust approach to analyse active fault sources and consequently a more refined seismic hazard assessment.     

Active tectonics of an apparently aseismic region: distributed active strike-slip faulting in the Hangay Mountains of central Mongolia

We identify and describe a series of east–west left-lateral strike-slip faults (named the Songino-Margats, the Hag Nuur, the Uliastay and the South Hangay fault systems) in the Hangay mountains of central Mongolia: an area that has little in the way of recorded seismicity and which is often considered as a rigid block within the India–Eurasia collision zone. The strike-slip faults of central Mongolia constitute a previously unrecognized hazard in this part of Mongolia. Each of the strike-slip faults show indications of late Quaternary activity in the form of aligned sequences of sag-ponds and pressure-ridges developed in alluvial deposits. Total bed-rock displacements of ∼3 km are measured on both the Songino-Margats and South Hangay fault systems. Bed-rock displacements of 11 km are observed across the Hag Nuur fault. Cumulative offset across the Uliastay fault systems are unknown but are unlikely to be large. We have no quantitative constraint on the age of faulting in the Hangay. The ≤20 km of cumulative slip on the Hangay faults might, at least in part, be inherited from earlier tectonic movements. Our observations show that, despite the almost complete absence of instrumentally recorded seismicity in the Hangay, this part of Mongolia is cut through by numerous distributed strike-slip faults that accommodate regional left-lateral shear between Siberia and China. Central Mongolia is thus an important component of the India–Eurasia collision that would be overlooked in models of the active tectonics based on the distribution of seismicity. We suggest that active faults such as those identified in the Hangay of Mongolia might exist in other, apparently aseismic, regions within continental collision zones.     

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.     

Plate generation and two-phase damage theory in a model of mantle convection

The formation of narrow, rapidly deforming plate boundaries separating strong plate interiors are integral components of the generation of plate tectonics from mantle convection. The development of narrow plate boundaries requires the interaction of a non-linear rheology and convection. One such non-linear rheology is two-phase damage theory which employs a non-equilibrium relation between interfacial surface energy, pressure and viscous deformation, thereby forming a theoretical model for void generation. Two-phase damage theory was recently extended to allow for deformational work to increase the fineness (reduce the grain size) of the matrix phase. We present results testing two-phase damage theory in a 2-D convectively driven system where we allow for (1) pure void-generating damage, (2) pure fineness-generating damage and (3) combined void- and fineness-generating damage. Pure void-generating damage is found to be unsuccessful at producing plate-like features. Fineness-generating damage is successful at inducing plate-like behaviour in certain circumstances, including increasing viscosity sensitivity to fineness and certain regimes of damage input and healing rate. Cases with combined void- and fineness-generating damage produce significantly more localization than the end-members due to the apparent increase of deformational work input into fineness generation. The interaction of microcracks and grain size reduction in two-phase damage theory suggests a rheological model for shear localization necessary for the formation of plate tectonic boundaries.     

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.     

Palaeomagnetism and the age of the Cracow volcanic rocks (S Poland)

Permian rhyodacites, melaphyres and tuffs from the Cracow area (South Poland) were sampled for the palaeomagnetic and isotope studies. Single-grain U-Pb dating of most zircon grains separated from the rhyodacites gave mean age of magma emplacement of 294.2 ± 2.1 Ma. Some zircons, however, displayed younger ages (268.7 ± 3.4 Ma), probably related to the metasomatic alterations of these rocks. Two Permian components of magnetizations related to these processes were isolated and together with previously defined Late Carboniferous–Permian palaeomagnetic poles from South Poland were used for construction of the regional apparent polar wander path (APWP). The Early Permian segment of this APWP shows a certain departure from the coeval part of the Fennoscandian APWP due to anticlockwise rotations of studied rocks most probably caused by mid-Permian sinistral tectonic movements along reactivated prominent Variscan faults of Central Europe. This sense of tectonic mobility does not support the hypothesis about transformation from Pangea 'B' to Pangea 'A' along an intra-Pangea dextral megashear during the Permian. Older than previously assumed ages of the post-Variscan igneous rocks of Central Europe reduce overlap of Gondwana's and Laurussia's parts of the Early Permian Pangea 'A'.     

Insights into the lithospheric structure and tectonic setting of the Barents Sea region from isostatic considerations

We study the tectonic setting and lithospheric structure of the greater Barents Sea region by investigating its isostatic state and its gravity field. 3-D forward density modelling utilizing available information from seismic data and boreholes shows an apparent shift between the level of observed and modelled gravity anomalies. This difference cannot be solely explained by changes in crustal density. Furthermore, isostatic calculations show that the present crustal thickness of 35–37 km in the Eastern Barents Sea is greater than required to isostatically balance the deep basins of the area (>19 km). To isostatically compensate the missing masses from the thick crust and deep basins and to adequately explain the gravity field, high-density material (3300–3350 kg m⁻³) in the lithospheric mantle below the Eastern Barents Sea is needed. The distribution of mantle densities shows a regional division between the Western and Eastern Barents and Kara Seas. In addition, a band of high-densities is observed in the lower crust along the transition zone from the Eastern to Western Barents Sea. The distribution of high-density material in the crust and mantle suggests a connection to the Neoproterozoic Timanide orogen and argues against the presence of a Caledonian suture in the Eastern Barents Sea. Furthermore, the results indicate that the basins of the Western Barents Sea are mainly affected by rifting, while the Eastern Barents Sea basins are located on a stable continental platform.