Slow Tectonic Deformation Research Articles

Exploring the Connection of XVI‐Century Major Historical Earthquakes in the Eastern Betic Cordillera, Spain: Insights From Viscoelastic Relaxation of the Lithosphere

AbstractUnderstanding the crustal fault interaction and connection between earthquakes in areas with slow tectonic deformation, such as Betic Cordillera (South Spain), is challenging. When seismic rates are low and large destructive earthquakes happen less frequently, it is necessary to resort to historical or paleoseismic records. This study investigates the postseismic viscoelastic relaxation mechanism as a potential explanation for the occurrence of three historical earthquakes (IEMSVIII‐IX) in the Eastern Betic Shear Zone during the XVI‐century, all of which occurred within a span of 13 years: 1518 Vera Mw6.2, 1522 Alhama de Almeria Mw6.5, and 1531 Baza Mw6.2 associated with the Palomares, Carboneras, and Baza faults, respectively. The results strongly suggest a sequential stress‐triggering connection between the three events. The northern NS‐oriented section of the Baza fault is found to have experienced a larger positive ΔCFS and, indeed, more prone to rupture in 1531. The study also examines whether the cumulative ΔCFS had influenced the occurrence of further significant earthquakes (≥Mw6.0) in the region. A triggering connection between the cascade and the 1658 Almeria Mw6.2 earthquake is suggested, whereas no indications of similar linkage to the 1674 Lorca Mw6.0 or the 1804 Dalias Mw6.4 events are found. The stress triggering impact of the cascade over nearby active faults is noteworthy. It is expected that this analysis could have future applications for studying other important historical events, and improving seismic hazard analysis in complex fault settings of the Betic Cordillera.

Comparative analysis of the structural pattern over the Indian Ocean basins using EIGEN6C4 Bouguer gravity data

The Indian Ocean witnessed numerous seismic and tectonic activities, making the region an intricate seismotectonic and seismo-genesis potential zone. The high-resolution global satellite-derived Bouguer Gravity Anomaly (BGA) data was used for comparative analysis of the structural patterns of six Indian Ocean basins. THETA (TA) edge enhancement and upward continuation techniques were used to study the Indian Ocean crustal structure by mapping both short and longwave lineaments. Rayleigh wave group velocity (RWGV) tomography at different time periods indicates pseudo depths of ∼30 km and ∼50 km are suitable for the upward continuation (UC) of BGA that perhaps lead to understanding different contrasting boundaries in their TA map. The UC of the BGA was generated at 60 km and 100 km for highlighting different lineaments/structures/contrasting geological boundaries at the depths of about 30 km and 50 km. The orientation of TA lineaments shows NNE-SSW, NE-SW, and EES-WWN structural trends. The statistical analysis of the lineaments shows that the Mid-Indian Basin constitutes the highest values of lineament number, lineament density (LD), circular variance (CV) and circular standard deviation (CSD), indicating a slow tectonic deformation over a long period. The Somalian Basin shows the minimum lineament number and lowest LD with high values of CV and CSD, indicating a relatively stable tectonic set-up compared to other Indian Ocean basins. The Arabian and Bay of Bengal basins offer relatively low lineament statistics than the Mid-Indian, Wharton and Madagascar basins, indicating frequent tectonic activity.

Multitechnology characterization of an unusual surface rupturing intraplate earthquake: the ML 5.4 2019 Le Teil event in France

SUMMARY Metropolitan France is a region of slow tectonic deformation with sparse seismicity. On 11 November 2019, the ML 5.4 Le Teil earthquake became the largest seismic event recorded in the last 16 yr. This event was recorded by the national seismic networks and also by a wide variety of other geophysical techniques including infrasound and InSAR measurements. These complementary technologies offer the opportunity to investigate in detail the earthquake source characteristics and the associated ground motion attenuation. Both seismic waveform inversions and InSAR interferogram reveal a shallow rupture on a reverse fault with an associated moment magnitude of 4.8–4.9. Infrasound signals also provide fast evidences pointing towards the area of ground surface displacements, which coincides with La Rouvière fault, in the Cévennes fault system, known as a formerly active normal fault during the Oligocene. The very significant amount of seismic records also helps towards validating the GMPE laws available for the region. This multitechnology characterization documents the kinematics of this rare example of shallow intraplate fault reactivation.

Assessing post-pliocene deformation in a context of slow tectonic deformation: insights from paleoseismology, remote sensing and shallow geophysics in Provence, France

The Provence region, located in the southeast of France, has experienced a few destructive earthquakes during the last centuries, such as the 1909 Lambesc earthquake or the 1509 and 1708 Manosque earthquakes. However, faults in the area experience slow slip rates (< 0.1 mm/yr) that preclude quantification of deformation using geodetic measurements. Active faults in Provence have long recurrence intervals, and the region experiences erosion under a Mediterranean climate where surface markers of such deformation are rapidly erased. As a consequence, several faults in the region may go unnoticed despite having the potential to generate earthquakes. This work focuses on the Vinon-sur-Verdon area in Provence, where a relatively narrow fault structure, the Maragrate fault, brings into contact Miocene clays and Pliocene conglomerates. Because of its proximity to a major nuclear research facility, this fault poses a potential earthquake hazard that justifies studies aimed at gaining insight into its mechanics and Quaternary seismic activity. In this context, a multidisciplinary approach has been implemented. A LiDAR campaign was undertaken to generate a 30-cm high-resolution Digital Elevation Model that depicts in detail the current fault morphology. In parallel, an electrical resistivity campaign was conducted and three trenches were opened. These investigations revealed a fault contact expressed by a fault gouge zone surmounted and sealed by numerous recent Quaternary periglacial and high-energy channel units. Radiocarbon and optically stimulated luminescence dating allowed the chronological reconstruction of sediment deposition and erosion cycles and of the most recent deformation event, older than 20 ka. These investigations did not allow us to demonstrate with certainty that these deformations are of a co-seismic nature. This work highlights the difficulty of identifying faults and studying the seismic hazard in intra-continental contexts.

Risk faults in stable crust of the eastern Bohemian Massif identified by integrating GNSS, levelling, geological, geomorphological and geophysical data

It is difficult to identify and assess the recently active faults, representing a risk factor for key human activities, in the areas of stable crust with slow tectonic deformation and reduced seismicity. We face the problem in the eastern Bohemian Massif (Central European Variscides) at a contact with the Eastern Alpine – Western Carpathian orogenic belt through building its integrated geodynamic model based on the long-term GNSS data cross-checked with the results from existing levelling, geological surveys, geophysical surveys, and a new geomorphological analysis. The model shows differently moving crustal blocks, determining the intervening boundary zones as the main risk regions, where weak earthquakes are also usually concentrated (primarily Diendorf–Čebín Tectonic Zone, Nectava–Konice Fault, Haná Faults, Bělá Fault and Bulhary Fault Shear Zone). The maximum horizontal GNSS differential velocities in the exposed Bohemian Massif reach up to 1.5–2.0 mm·yr−1. The up-thrusted segments of the Western Carpathians move individually ~2 mm·yr−1. The integration of a number of full-area datasets, corresponding to various dimensions and depth levels of crustal processes, allowed us to highlight the important fault zones as driving elements of regional geodynamics.

Focal mechanisms and inter-event times of low-frequency earthquakes reveal quasi-continuous deformation and triggered slow slip on the deep Alpine Fault

Characterising the seismicity associated with slow deformation in the vicinity of the Alpine Fault may provide constraints on the stresses acting on a major transpressive margin prior to an anticipated great (≥M8) earthquake. Here, we use recently detected tremor and low-frequency earthquakes (LFEs) to examine how slow tectonic deformation is loading the Alpine Fault late in its typical ∼300-yr seismic cycle. We analyse a continuous seismic dataset recorded between 2009 and 2016 using a network of 10–13 short-period seismometers, the Southern Alps Microearthquake Borehole Array. Fourteen primary LFE templates are used in an iterative matched-filter and stacking routine, allowing the detection of similar signals corresponding to LFE families sharing common locations. This yields an 8-yr catalogue containing 10,000 LFEs that are combined for each of the 14 LFE families using phase-weighted stacking to produce signals with the highest possible signal-to-noise ratios. We show that LFEs occur almost continuously during the 8-yr study period and highlight two types of LFE distributions: (1) discrete behaviour with an inter-event time exceeding 2 min; (2) burst-like behaviour with an inter-event time below 2 min. We interpret the discrete events as small-scale frequent deformation on the deep extent of the Alpine Fault and LFE bursts (corresponding in most cases to known episodes of tremor or large regional earthquakes) as brief periods of increased slip activity indicative of slow slip. We compute improved non-linear earthquake locations using a 3-D velocity model. LFEs occur below the seismogenic zone at depths of 17–42 km, on or near the hypothesised deep extent of the Alpine Fault. The first estimates of LFE focal mechanisms associated with continental faulting, in conjunction with recurrence intervals, are consistent with quasi-continuous shear faulting on the deep extent of the Alpine Fault.

Educating for Earthquake Science and Risk in a Tectonically Slowly Deforming Region

Over the past decade, scientists have been called to participate more actively in public education and outreach (E&O). This is particularly true in fields of significant societal impact, such as earthquake science. Local earthquake risk culture plays a role in the way that the public engages in educational efforts. In this article, we describe an adapted E&O program for earthquake science and risk. The program is tailored for a region of slow tectonic deformation, where large earthquakes are extreme events that occur with long return periods. The adapted program has two main goals: (1) to increase the awareness and preparedness of the population to earthquake and related risks (tsunami, liquefaction, fires, etc.), and (2) to increase the quality of earthquake science education, so as to attract talented students to geosciences. Our integrated program relies on activities tuned for different population groups who have different interests and abilities, namely young children, teenagers, young adults, and professionals.

Contemporary kinematics of the Upper Rhine Graben: A 3D finite element approach

The Upper Rhine Graben (URG), a Cenozoic intra-plate rift situated in the Alpine foreland, is presently characterised by relative slow tectonic deformation and low to medium seismicity. Concurrently, it is a region with a significant amount of ongoing subsidence in two recent depocentres (0.1 to 0.2 mm/a geological, 1 mm/a geodetical rate). In this paper, the recent kinematic behaviour of the URG is simulated using a 3D finite element model, containing three lithospheric layers (upper mantle, lower crust and upper crust) with different rheological properties. First order fault structures (e.g. border faults) are implemented as frictional contact surfaces within the upper crustal layer. The stresses generated by applying lateral displacements over a time period of 10 ka are insufficient to obtain a match between predicted and observed stress magnitudes. Therefore, a technique of “combined pre-stressing” has been developed to avoid unrealistic deformation and unrealistic stress magnitudes within the model. The stress magnitudes and stress directions predicted are calibrated against in-situ stress measurements and stress indicator data. For benchmarking of the modelling results, the vertical surface displacements predicted are compared to surface uplift derived from geological and geomorphological data. Furthermore, predicted fault slip rates are compared to available geological and geodetical data. Parameters derived from the calculated stress tensor, such as fracture potential and the regime stress ratio are also analysed in order to describe the possible kinematic behaviour of the URG. The modelling results suggest that the URG is currently being reactivated as a sinistral strike–slip system with the central segment of the URG forming a restraining bend and the two recent depocentres situated in releasing bend settings. The modelling results suggest that both sinistral shearing and mantle uplift are active mechanisms driving the recent kinematics of the URG and that the recent subsidence within the two depocentres is re-enforced by ongoing mantle uplift additionally.

Pleistocene tectonics inferred from fluvial terraces of the northern Upper Rhine Graben, Germany

This study of fluvial terraces of the River Rhine and tributaries aims to search for indications of Pleistocene tectonic activity. The study area includes the northern Upper Rhine Graben (URG), the Mainz Basin and the adjacent Rhenish Massif with the Middle Rhine Valley. High rates of Quaternary surface processes, large amount of human modifications, relatively slow tectonic deformation and presently low intra-plate seismic activity characterize this area. Therefore, the records of relatively slow tectonic deformation are less well preserved and thus difficult to detect. This study uses the relative position of fluvial terraces to determine the more local effects of fault movements on the terraces and to evaluate their displacement rates and patterns. The research is based on a review of previous terrace studies and new terrace mapping from the eastern Mainz Basin and the bordering URG using topographic map interpretations and field observations. This newly mapped sequence of terrace surfaces can be correlated to other terraces in the vicinity on the basis of relative height levels. Terrace correlation between the western Mainz Basin and Middle Rhine Valley relies on a single chronostratigraphic unit (Mosbach sands) and additional relative height correlations. This is the first study to present a continuous correlation of terraces from the western margin of the URG to the Rhenish Massif and enables the study of the transition from the subsiding graben to the uplifted Rhenish Massif. By means of a longitudinal profile, which ranges from the URG to the Rhenish Massif, the influence of individual fault movements on the terrace levels and the large-scale regional uplift is demonstrated. It is evident from the profile that the uplift of Early to Middle Pleistocene terraces increases northwards, towards the Rhenish Massif. The uplift was diachronic, with a significant pulse occurring first in the northern URG (Lower Pleistocene) and later in the Rhenish Massif (Middle Pleistocene). The largest vertical displacements are recorded for the boundary fault separating the Mainz Basin and the Rhenish Massif (Hunsrück–Taunus Boundary Fault) and for faults bounding the northeastern Mainz Basin. The motions and displacement rates calculated for individual faults indicate deformation rates in the order of 0.01–0.08 mm/year. At this stage, the calculation of displacement rates depends mostly on a single dated stratigraphic unit. Additional dating of terrace deposits is urgently needed to better constrain the temporal development of the terrace sequence and the impact of tectonic movements.

Late Cenozoic geomorphologic signal of Andean forearc deformation and tilting associated with the uplift and climate changes of the Southern Atacama Desert (26°S–28°S)

We analyze remarkable examples of the large (∼ 10,000 km 2) and local-scale (∼ 100 km 2) landscape forms related to Late Cenozoic geomorphologic evolution of the Andean forearc region in the Southern Atacama Desert. We also consider the continental sedimentary deposits, so-called “Atacama Gravels”, which are related to the degradation of the landscape during the Neogene. Our analysis integrates 1:50,000 field cartography, Landsat TM images observations, ∼ 1:1000 sedimentary logging data, and 50 m horizontal resolution topographic data to reconstruct the Late Cenozoic geomorphologic evolution of this region and discuss the factors that control it, i.e., Miocene aridification of the climate and Neogene Central Andean uplift. We determine that the Precordillera was already formed in the Oligocene and most of the present-day altitude of the Precordillera was reached before that time. Afterward, five episodes of geomorphologic evolution can be differentiated: (1) the development of an Oligocene deep incised drainage system cutting the uplifted Precordillera (up to 2000 m of vertical incision) and connecting it to the Ocean; followed by (2) the infilling of deep incised valleys by up to 400 m of Atacama Gravels. This infill started in the Early Miocene with the development of fluvial deposition and finished in the Middle Miocene with playa and playa lake depositions. We propose that playa-related deposition occurs in an endorheic context related to tectonic activity of the Atacama Fault System and Coastal Cordillera uplift. However, the upward sedimentologic variation in the Atacama Gravels evidences a progressive aridification of the climate. Subsequently, we have identified the effects of the Middle–Upper Miocene slow tectonic deformation: the Neogene Andean uplift is accommodated by a tilting or flexuring of the inner-forearc (Central Depression and Precordillera) related to some hundreds of meters of uplift in the Precordillera. This tilting or flexuring results in (3) the Middle Miocene re-opening of the valley network to the Pacific Ocean. Upper Miocene aridification, from arid to hyperarid, induces alluvial fans backfilling in the Central Depression (4) resulting in up to 50 m of Atacama Gravel deposition. Finally, in response to an increase in the rate of tilting, a new phase of vertical incision (up to 800 m in the Precordillera) allows the development of the canyon that crosses the forearc (5).

Experimental reproduction of tectonic deformation lamellae in quartz and comparison to shock‐induced planar deformation features

Abstract— Planar features can develop in quartz during comparatively slow tectonic deformation and during very fast dynamic shock metamorphism. Despite their very different structural nature, tectonically induced deformation lamellae have sometimes been mistaken as shock‐induced planar deformation features (PDFs). To understand the formation of deformation lamellae and to address the substantial differences between them and PDFs, we have conducted deformation experiments on single crystals of quartz in a Griggs‐type apparatus, at a temperature of 800 °C, a confining pressure of 12 kbar, and a strain rate of 0.7–1.1 · 10−6. The deformed samples were analyzed with transmission electron microscopy (TEM) and compared to natural PDFs from the Ries Crater, Germany. TEM revealed that tectonic deformation lamellae are associated with numerous sub‐parallel curved subgrain walls, across which the orientation of the crystal changes slightly. The formation of deformation lamellae is due to glide‐ and climb‐controlled deformation in the exponential creep regime. In contrast, the PDFs in shocked quartz from the Ries are perfectly planar, crystallographically controlled features that originally represented amorphous lamellae. Due to post‐shock annealing and hydrothermal activity they are recrystallized and decorated with fluid inclusions.

Interplay between tectonic, fluvial and erosional processes along the Western Border Fault of the northern Upper Rhine Graben, Germany

The northern Upper Rhine Graben, situated in the central part of the European Cenozoic rift system, is currently characterized by low intra-plate seismicity. Historical earthquakes have not been large enough to produce surface rupturing. Moreover, the records of Quaternary surface processes and human modifications are presumably better preserved than the record of the relatively slow tectonic deformation. In order to gain information on the neotectonic activity and paleoseismicity in this setting, the geological and geomorphological records of fault movements along a segment of the Western Border Fault (WBF) were studied using an integration of techniques in paleoseismology, structural analysis and shallow geophysics. The WBF segment investigated follows a 20 km long linear scarp of unclear origin. A series of geophysical measurements were performed and the results suggested that near-surface deformation structures are present at the segments' southern end. Several trenches opened at this location revealed fault structures with consistent extensional style and a maximum vertical displacement of 0.5 m. In one trench, the deformation structures were dated between 19 and 8 ka. Assuming the deformation has been caused by an earthquake, a M w ∼6.5 earthquake would be implied. Aseismic deformation would point to a fault creep rate ≥ 0.04 mm/yr. A reconstruction of the sequence of events at the trench site, from Middle Pleistocene to Present, demonstrates that the morphology at the base of the scarp is the result of interplay between tectonic activity and fluvial and erosional processes. At the regional scale, a mixed origin for the WBF scarp is proposed, combining the effects of fluvial dynamics, erosion, regional uplift and localized tectonic activity on the WBF.

The influence of water on the stress supported by experimentally faulted Westerly granite

Summary. Fault zones in wet Westerly granite deformed at temperatures of 300 and 400°C require markedly lower shear stresses for sliding than when dry, even when the effective confining pressure is held constant between the wet and dry tests, provided that the strain rate is lower than lO-'s-l. The rate of strength reduction is enhanced by increasing the pore water pressure. The deformation rate is a power function of the applied stress where the stress exponent is approximately 7 for pore water pressure of l00MPa and 21 for pore water pressure of 20 MPa. The experimental results are extrapolated to conditions believed to occur at depths of 10 km along the San Andreas Fault Zone. It is suggested that for slow tectonic deformation at strain rates of lo-'' and 10-14s-' the shear stress for sliding on faults in granite would be approximately 60 and 20MPa, respectively, at pore water pressures equal to the hydrostatic head. Fluid overpressures of c. 0.8 lithostatic pressure are required to lower the shear stress for sliding into the 10MPa range at the slower strain rate.

Strainmeters and tiltmeters in geophysics

Several types of sensitive strainmeters and tiltmeters have been developed, and it is now becoming clear which geophysical applications are most suitable for these instruments. In general, strainmeters and tiltmeters are used for observing ground deformation at periods of minutes to days. Small-scale lateral inhomogeneities at the instrument sites distort signals by a few percent, although the effects of large structures can be calculated. In earth tide work these lateral inhomogeneities and unknown ocean loading signals prevent accurate values of the regional tide from being obtained. This limits tidal investigations to looking for temporal variations, possibly associated with pre-earthquake dilatancy, and spatial variations caused by gross elasticity contrasts in the local geological structure. Strainmeters and tiltmeters are well suited for observing long-period seismic waves, seismic slip events on faults and volcano tumescence, where small site-induced distortions in the measured signals are seldom important. Geodetic techniques are superior for monitoring slow tectonic deformation over large areas, but strainmeters and tiltmeters may be used to complement geodetic networks by providing continuous records. The recording instruments would show whether deformation proceeded episodically, but could not be relied upon to give accurate values for secular strain and tilt rates because of instrumental drift.