TARDI VARISCAN DEFORMATION IN THE IBERIAN PENINSULA; A LATE FEATURE IN THE LAURASIA-GONDWANA DEXTRAL COLLISION

The Late Variscan deformation event in Iberia, is characterized by an intraplate deformation regime induced by the oblique collision between Laurentia and Gondwan. This episode in Iberia is characterized by NNE-SSW strike-slip faults, which are considered by the classic works as sinistral strike-slips. However, the absence of Mesozoic formations constraining the age of this sinistral kinematics, led some authors to consider it as the result of Alpine reworking. Structural studies in Almograve and Ponta Ruiva sectors (SW Portugal), not only shows that NNE-SSW faults presents a clear sinistral kinematics and are occasionally associated with E-W dextral shears, but also that this kinematics is related to the late deformation episodes of Variscan Orogeny. In Almograve sector, the late Variscan structures are characterized by NNE-SSW sinistral kink-bands, spatially associated with E-W dextral faults. These structures are contemporaneous and affect the previously deformed Carboniferous units. The Ponta Ruiva Sector constrains the age of deformation because the E-W dextral shears affect the Late Carboniferous (late Moscovian) units, but not the overlying Triassic series. The new exposed data shows that the NNE-SSW and the E-W faults are dynamically associated and results from the same deformation event. The NNE-SSW sinistral faults could be considered as second order dominoes structures related with first order E-W dextral shears, related with Laurasia-Gondwana collision during Late Carboniferous-Permian Times.

Tectonic evolution of Variscan Iberia: Gondwana–Laurussia

An integrated interpretation of the late Paleozoic structural and geochronological record of the Iberian Massif is presented and discussed under the perspective of a Gondwana-Laurussia collision giving way to the Variscan orogen. Compressional and extensional structures developed during the building of the Variscan orogenic crust of Iberia are linked together into major tectonic events operating at lithosphere scale. A review of the tectonometamorphic and magmatic evolution of the IberianMassif reveals backs and forths in the overall conver- gence between Gondwana and Laurussia during theamalgamation of Pangea in late Paleozoic times. Stages dom- inated by lithosphere compression are characterized by subduction, both oceanic and continental, development of magmatic arcs, (over- and under-) thrusting of continental lithosphere, and folding. Variscan convergence re- sulted in the eventual transference of a large allochthonous set of peri-Gondwanan terranes, the Iberian Allochthon, onto the Gondwana mainland. The Iberian Allochthon bears the imprint of previous interaction be- tween Gondwana and Laurussia, including their juxtaposition after the closure of the Rheic Ocean in Lower De- vonian times. Stages governed by lithosphere extension are featured by the opening of two short-lived oceanic basins that dissected previous Variscan orogenic crust, first in the Lower-Middle Devonian, following the closure of the Rheic Ocean, and then in the early Carboniferous, following the emplacement of the peri-Gondwanan allochthon. An additional, major intra-orogenic extensional event in the early-middle Carboniferous dismem- bered the Iberian Allochthon into individual thrust stacks separated by extensional faults and domes. Lateral tec- tonics played an important role through the Variscan orogenesis, especially during the creation of new tectonic blocks separated by intracontinental strike-slip shear zones in the late stages of continental convergence.

Superposition relations of microfabrics in the northern hanging-wall block

This study reports alternation of D2 extension-related and D3 contraction-related microfabrics in the northern hanging wall block of a gneiss dome-like structure recognized in the Évora Massif (Ossa-Morena Zone). In the Arraiolos – Santo Antonio de Alcorrego traverse high- to low-grade mylonites are dominant. Microfabrics related to D2 ductile deformation and M2 high-amphibolite to greenschist facies characterize an extensional shear zone with telescoping metamorphic isograds. D2 microstructures indicate shear sense with top-to-SE. Superposition of D3 contraction developed under greenschist facies (M3) producing folding of D2 microfabrics, mylonitization of granites along strike-slip shear zones and retrogression of M2 mineral assemblages.

Slip complexity in earthquake fault models.

We summarize studies of earthquake fault models that give rise to slip complexities like those in natural earthquakes. For models of smooth faults between elastically deformable continua, it is critical that the friction laws involve a characteristic distance for slip weakening or evolution of surface state. That results in a finite nucleation size, or coherent slip patch size, h*. Models of smooth faults, using numerical cell size properly small compared to h*, show periodic response or complex and apparently chaotic histories of large events but have not been found to show small event complexity like the self-similar (power law) Gutenberg-Richter frequency-size statistics. This conclusion is supported in the present paper by fully inertial elastodynamic modeling of earthquake sequences. In contrast, some models of locally heterogeneous faults with quasi-independent fault segments, represented approximately by simulations with cell size larger than h* so that the model becomes "inherently discrete," do show small event complexity of the Gutenberg-Richter type. Models based on classical friction laws without a weakening length scale or for which the numerical procedure imposes an abrupt strength drop at the onset of slip have h* = 0 and hence always fall into the inherently discrete class. We suggest that the small-event complexity that some such models show will not survive regularization of the constitutive description, by inclusion of an appropriate length scale leading to a finite h*, and a corresponding reduction of numerical grid size.

3D geometry of growth strata in a fault-propagation fold: insights into space–time evolution of the Crevillente Fault (Abanilla-Alicante sector), Betic Cordillera, Spain

This work presents a 3D geometric model of growth strata cropping out in a fault-propagation fold associated with the Crevillente Fault (Abanilla-Alicante sector) from the Bajo Segura Basin (eastern Betic Cordillera, southern Spain). The analysis of this 3D model enables us to unravel the along-strike and along-section variations of the growth strata, providing constraints to assess the fold development, and hence, the fault kinematic evolution in space and time. We postulate that the observed along-strike dip variations are related to lateral variation in fault displacement. Along-section variations of the progressive unconformity opening angles indicate greater fault slip in the upper Tortonian–Messinian time span; from the Messinian on, quantitative analysis of the unconformity indicate a constant or lower tectonic activity of the Crevillente Fault (Abanilla-Alicante sector); the minor abundance of striated pebbles in the Pliocene-Quaternary units could be interpreted as a decrease in the stress magnitude and consequently in the tectonic activity of the fault. At a regional scale, comparison of the growth successions cropping out in the northern and southern limits of the Bajo Segura Basin points to a southward migration of deformation in the basin. This means that the Bajo Segura Fault became active after the Crevillente Fault (Abanilla-Alicante sector), for which activity on the latter was probably decreasing according to our data. Consequently, we propose that the seismic hazard at the northern limit of the Bajo Segura Basin should be lower than at the southern limit.

Finite element analysis of fault bend influence on stick-slip instability along an intra-plate fault

Earthquakes have been recognized as resulting from stick-slip frictional instabilities along the faults between deformable rocks. A three-dimensional finite-element code for modeling the nonlinear frictional contact behaviors between deformable bodies with the node-to-point contact element strategy has been developed and applied here to investigate the fault geometry influence on the nucleation and development process of the stick-slip instability along an intra-plate fault through a typical fault bend model, which has a pre-cut fault that is artificially bent by an angle of 5.6degrees at the fault center. The numerical results demonstrate that the geometry of the fault significantly affects nucleation, termination and restart of the stick-slip instability along the intra-plate fault, and all these instability phenomena can be well simulated using the current finite-element algorithm.

On the stacking fault energies of some close-packed hexagonal metals

It is now realised (1,2,3) that a knowledge of stacking fault energy is fundamental for an understanding of the mechanical behaviour of metals. There are several processes in which the imperfect dislocations have to recombine locally to form an unextended dislocation . For intersection of two dislocations it is, for example, necessary to form 'constrictions'. Cross slip of extended dislocations also involves constriction. The onset of stage llI work hardening in a crystal with close-packed structure is attributed to cross slip and hence is controlled by the stacking fault energy (SPE). Methods of estimation of SFE are based on either the direct observation of stacking faults in an electron microscope or their effects on the deformation processes.

Deformation of nanograined Ni-60Co alloy with low stacking fault energy

The evolution of deformation texture in a Ni-60Co alloy with low stacking fault energy and a grain size in the nanometre range has been investigated. The analyses of texture and microstructure suggest different mechanisms of deformation in nanocrystalline as compared to microcrystalline Ni-60Co alloy. In nanocrystalline material, the mechanism responsible for texture formation has been identified as partial slip, whereas in microcrystalline material, a characteristic texture forms due to twinning and shear banding.

Role of stacking fault energy on texture evolution revisited

Three materials, pure aluminium, Al-4 wt.% Mg, alpha-brass have been chosen to understand the evolution of texture and microstructure during rolling. Pure Al develops a strong copper-type rolling texture and the deformation is entirely slip dominated. In Al-4Mg alloy, texture is copper-type throughout the deformation. The advent of Cu-type shear bands in the later stages of deformation has a negligible effect on the final texture. alpha-brass shows a characteristic brass-type texture from the early stages of rolling. Extensive twinning in the intermediate stages of deformation (epsilon(t) similar to 0.5) causes significant texture reorientation towards alpha-fiber. Beyond 40% reduction, deformation is dominated by Bs-type shear bands, and the banding coincides with the evolution of <111>parallel to ND components. The crystallites within the bands preferentially show <110>parallel to ND components. The absence of the Cu component throughout the deformation process indicates that, for the evolution of brass-type texture, the presence of Cu component is not a necessary condition. The final rolling texture is a synergistic effect of deformation twinning and shear banding.

Mechanism on Progressive Failure of a Faulted Rock Slope Due to Slip-Weakening

Slip-weakening is one of the characteristics of geological materials under certain loadings. Non-uniform rock structure may exist in the vicinity of the slip surface for a rock slope. Some portion of the slip surface may be penetrated but the other not. For the latter case, the crack or the fault surface will undergo shear deformation before it becomes a successive surface under a certain loading. As the slipped portion advances,slip-weakening occurs over a distance behind the crack tip. In the weakening zone, the shear strength will decrease from its peak value to residual friction level. The stress will redistribute along the surface of crack and in the weakening zone. Thus the changed local stress concentration leads the crack to extend and the ratio of penetration of the slip surface to increase. From the view of large-scale for the whole slip surface, the shear strength will decrease due to the damage of interior rock structure, and the faulted rock behaves as a softening material. Such a kind of mechanism performs in a large number of practical landslides in the zones experienced strong earthquakes. It should be noted that the mechanism mentioned above is different from that of the breakage of structural clay,in which the geological material is regarded as a medium containing structural lumps and structural bands. In this paper, the softening behavior of a faulted rock should be regarded as a comprehensive result of the whole complicated process including slip-weakening, redistribution of stress, extension of crack tip, and the penetration of the slip surface. This process is accompanied by progressive failure and abrupt structural damage. The size of slip-weakening zone is related to the undergoing strain. Once the relative slide is initiated (local or integrated), the effect of slip-weakening will behave in a certain length behind the crack tip until the formation of the whole slip surface.

Earthquakes of the Nepal Himalaya : towards a physical model of the seismic cycle

Home to hundreds of millions of souls and land of excessiveness, the Himalaya is also the locus of a unique seismicity whose scope and peculiarities still remain to this day somewhat mysterious. Having claimed the lives of kings, or turned ancient timeworn cities into heaps of rubbles and ruins, earthquakes eerily inhabit Nepalese folk tales with the fatalistic message that nothing lasts forever. From a scientific point of view as much as from a human perspective, solving the mysteries of Himalayan seismicity thus represents a challenge of prime importance. Documenting geodetic strain across the Nepal Himalaya with various GPS and leveling data, we show that unlike other subduction zones that exhibit a heterogeneous and patchy coupling pattern along strike, the last hundred kilometers of the Main Himalayan Thrust fault, or MHT, appear to be uniformly locked, devoid of any of the “creeping barriers” that traditionally ward off the propagation of large events. The approximately 20 mm/yr of reckoned convergence across the Himalaya matching previously established estimates of the secular deformation at the front of the arc, the slip accumulated at depth has to somehow elastically propagate all the way to the surface at some point. And yet, neither large events from the past nor currently recorded microseismicity nearly compensate for the massive moment deficit that quietly builds up under the giant mountains. Along with this large unbalanced moment deficit, the uncommonly homogeneous coupling pattern on the MHT raises the question of whether or not the locked portion of the MHT can rupture all at once in a giant earthquake. Univocally answering this question appears contingent on the still elusive estimate of the magnitude of the largest possible earthquake in the Himalaya, and requires tight constraints on local fault properties. What makes the Himalaya enigmatic also makes it the potential source of an incredible wealth of information, and we exploit some of the oddities of Himalayan seismicity in an effort to improve the understanding of earthquake physics and cipher out the properties of the MHT. Thanks to the Himalaya, the Indo-Gangetic plain is deluged each year under a tremendous amount of water during the annual summer monsoon that collects and bears down on the Indian plate enough to pull it away from the Eurasian plate slightly, temporarily relieving a small portion of the stress mounting on the MHT. As the rainwater evaporates in the dry winter season, the plate rebounds and tension is increased back on the fault. Interestingly, the mild waggle of stress induced by the monsoon rains is about the same size as that from solid-Earth tides which gently tug at the planets solid layers, but whereas changes in earthquake frequency correspond with the annually occurring monsoon, there is no such correlation with Earth tides, which oscillate back-and-forth twice a day. We therefore investigate the general response of the creeping and seismogenic parts of MHT to periodic stresses in order to link these observations to physical parameters. First, the response of the creeping part of the MHT is analyzed with a simple spring-and-slider system bearing rate-strengthening rheology, and we show that at the transition with the locked zone, where the friction becomes near velocity neutral, the response of the slip rate may be amplified at some periods, which values are analytically related to the physical parameters of the problem. Such predictions therefore hold the potential of constraining fault properties on the MHT, but still await observational counterparts to be applied, as nothing indicates that the variations of seismicity rate on the locked part of the MHT are the direct expressions of variations of the slip rate on its creeping part, and no variations of the slip rate have been singled out from the GPS measurements to this day. When shifting to the locked seismogenic part of the MHT, spring-and-slider models with rate-weakening rheology are insufficient to explain the contrasted responses of the seismicity to the periodic loads that tides and monsoon both place on the MHT. Instead, we resort to numerical simulations using the Boundary Integral CYCLes of Earthquakes algorithm and examine the response of a 2D finite fault embedded with a rate-weakening patch to harmonic stress perturbations of various periods. We show that such simulations are able to reproduce results consistent with a gradual amplification of sensitivity as the perturbing period get larger, up to a critical period corresponding to the characteristic time of evolution of the seismicity in response to a step-like perturbation of stress. This increase of sensitivity was not reproduced by simple 1D-spring-slider systems, probably because of the complexity of the nucleation process, reproduced only by 2D-fault models. When the nucleation zone is close to its critical unstable size, its growth becomes highly sensitive to any external perturbations and the timings of produced events may therefore find themselves highly affected. A fully analytical framework has yet to be developed and further work is needed to fully describe the behavior of the fault in terms of physical parameters, which will likely provide the keys to deduce constitutive properties of the MHT from seismological observations.

Characterization of diverse megathrust fault behavior related to seismic supercycles, Mentawai Islands, Sumatra

Long paleoseismic histories are necessary for understanding the full range of behavior of faults, as the most destructive events often have recurrence intervals longer than local recorded history. The Sunda megathrust, the interface along which the Australian plate subducts beneath Southeast Asia, provides an ideal natural laboratory for determining a detailed paleoseismic history over many seismic cycles. The outer-arc islands above the seismogenic portion of the megathrust cyclically rise and subside in response to processes on the underlying megathrust, providing uncommonly good illumination of megathrust behavior. Furthermore, the growth histories of coral microatolls, which record tectonic uplift and subsidence via relative sea level, can be used to investigate the detailed coseismic and interseismic deformation patterns. One particularly interesting area is the Mentawai segment of the megathrust, which has been shown to characteristically fail in a series of ruptures over decades, rather than a single end-to-end rupture. This behavior has been termed a seismic “supercycle.” Prior to the current rupture sequence, which began in 2007, the segment previously ruptured during the 14th century, the late 16th to late 17th century, and most recently during historical earthquakes in 1797 and 1833. In this study, we examine each of these previous supercycles in turn.

First, we expand upon previous analysis of the 1797–1833 rupture sequence with a comprehensive review of previously published coral microatoll data and the addition of a significant amount of new data. We present detailed maps of coseismic uplift during the two great earthquakes and of interseismic deformation during the periods 1755–1833 and 1950–1997 and models of the corresponding slip and coupling on the underlying megathrust. We derive magnitudes of Mw 8.7–9.0 for the two historical earthquakes, and determine that the 1797 earthquake fundamentally changed the state of coupling on the fault for decades afterward. We conclude that while major earthquakes generally do not involve rupture of the entire Mentawai segment, they undoubtedly influence the progression of subsequent ruptures, even beyond their own rupture area. This concept is of vital importance for monitoring and forecasting the progression of the modern rupture sequence.

Turning our attention to the 14th century, we present evidence of a shallow slip event in approximately A.D. 1314, which preceded the “conventional” megathrust rupture sequence. We calculate a suite of slip models, slightly deeper and/or larger than the 2010 Pagai Islands earthquake, that are consistent with the large amount of subsidence recorded at our study site. Sea-level records from older coral microatolls suggest that these events occur at least once every millennium, but likely far less frequently than their great downdip neighbors. The revelation that shallow slip events are important contributors to the seismic cycle of the Mentawai segment further complicates our understanding of this subduction megathrust and our assessment of the region’s exposure to seismic and tsunami hazards.

Finally, we present an outline of the complex intervening rupture sequence that took place in the 16th and 17th centuries, which involved at least five distinct uplift events. We conclude that each of the supercycles had unique features, and all of the types of fault behavior we observe are consistent with highly heterogeneous frictional properties of the megathrust beneath the south-central Mentawai Islands. We conclude that the heterogeneous distribution of asperities produces terminations and overlap zones between fault ruptures, resulting in the seismic “supercycle” phenomenon.

Frictional properties of faults : from observation on the Longitudinal Valley Fault, Taiwan, to dynamic simulations

Faults can slip either aseismically or through episodic seismic ruptures, but we still do not understand the factors which determine the partitioning between these two modes of slip. This challenge can now be addressed thanks to the dense set of geodetic and seismological networks that have been deployed in various areas with active tectonics. The data from such networks, as well as modern remote sensing techniques, indeed allow documenting of the spatial and temporal variability of slip mode and give some insight. This is the approach taken in this study, which is focused on the Longitudinal Valley Fault (LVF) in Eastern Taiwan. This fault is particularly appropriate since the very fast slip rate (about 5 cm/yr) is accommodated by both seismic and aseismic slip. Deformation of anthropogenic features shows that aseismic creep accounts for a significant fraction of fault slip near the surface, but this fault also released energy seismically, since it has produced five M_w>6.8 earthquakes in 1951 and 2003. Moreover, owing to the thrust component of slip, the fault zone is exhumed which allows investigation of deformation mechanisms. In order to put constraint on the factors that control the mode of slip, we apply a multidisciplinary approach that combines modeling of geodetic observations, structural analysis and numerical simulation of the "seismic cycle". Analyzing a dense set of geodetic and seismological data across the Longitudinal Valley, including campaign-mode GPS, continuous GPS (cGPS), leveling, accelerometric, and InSAR data, we document the partitioning between seismic and aseismic slip on the fault. For the time period 1992 to 2011, we found that about 80-90% of slip on the LVF in the 0-26 km seismogenic depth range is actually aseismic. The clay-rich Lichi M\'elange is identified as the key factor promoting creep at shallow depth. Microstructural investigations show that deformation within the fault zone must have resulted from a combination of frictional sliding at grain boundaries, cataclasis and pressure solution creep. Numerical modeling of earthquake sequences have been performed to investigate the possibility of reproducing the results from the kinematic inversion of geodetic and seismological data on the LVF. We first investigate the different modeling strategy that was developed to explore the role and relative importance of different factors on the manner in which slip accumulates on faults. We compare the results of quasi dynamic simulations and fully dynamic ones, and we conclude that ignoring the transient wave-mediated stress transfers would be inappropriate. We therefore carry on fully dynamic simulations and succeed in qualitatively reproducing the wide range of observations for the southern segment of the LVF. We conclude that the spatio-temporal evolution of fault slip on the Longitudinal Valley Fault over 1997-2011 is consistent to first order with prediction from a simple model in which a velocity-weakening patch is embedded in a velocity-strengthening area.

Experimental Investigation of Thrust Fault Rupture Mechanics

Thrust fault earthquakes are investigated in the laboratory by generating dynamic shear ruptures along pre-existing frictional faults in rectangular plates. A considerable body of evidence suggests that dip-slip earthquakes exhibit enhanced ground motions in the acute hanging wall wedge as an outcome of broken symmetry between hanging and foot wall plates with respect to the earth surface. To understand the physical behavior of thrust fault earthquakes, particularly ground motions near the earth surface, ruptures are nucleated in analog laboratory experiments and guided up-dip towards the simulated earth surface. The transient slip event and emitted radiation mimic a natural thrust earthquake. High-speed photography and laser velocimeters capture the rupture evolution, outputting a full-field view of photo-elastic fringe contours proportional to maximum shearing stresses as well as continuous ground motion velocity records at discrete points on the specimen. Earth surface-normal measurements validate selective enhancement of hanging wall ground motions for both sub-Rayleigh and super-shear rupture speeds. The earth surface breaks upon rupture tip arrival to the fault trace, generating prominent Rayleigh surface waves. A rupture wave is sensed in the hanging wall but is, however, absent from the foot wall plate: a direct consequence of proximity from fault to seismometer. Signatures in earth surface-normal records attenuate with distance from the fault trace. Super-shear earthquakes feature greater amplitudes of ground shaking profiles, as expected from the increased tectonic pressures required to induce super-shear transition. Paired stations measure fault parallel and fault normal ground motions at various depths, which yield slip and opening rates through direct subtraction of like components. Peak fault slip and opening rates associated with the rupture tip increase with proximity to the fault trace, a result of selective ground motion amplification in the hanging wall. Fault opening rates indicate that the hanging and foot walls detach near the earth surface, a phenomenon promoted by a decrease in magnitude of far-field tectonic loads. Subsequent shutting of the fault sends an opening pulse back down-dip. In case of a sub-Rayleigh earthquake, feedback from the reflected S wave re-ruptures the locked fault at super-shear speeds, providing another mechanism of super-shear transition.