Normal fault earthquakes or graviquakes

Earthquakes are dissipation of energy throughout elastic waves. Canonically is the elastic energy accumulated during the interseismic period. However, in crustal extensional settings, gravity is the main energy source for hangingwall fault collapsing. Gravitational potential is about 100 times larger than the observed magnitude, far more than enough to explain the earthquake. Therefore, normal faults have a different mechanism of energy accumulation and dissipation (graviquakes) with respect to other tectonic settings (strike-slip and contractional), where elastic energy allows motion even against gravity. The bigger the involved volume, the larger is their magnitude. The steeper the normal fault, the larger is the vertical displacement and the larger is the seismic energy released. Normal faults activate preferentially at about 60° but they can be shallower in low friction rocks. In low static friction rocks, the fault may partly creep dissipating gravitational energy without releasing great amount of seismic energy. The maximum volume involved by graviquakes is smaller than the other tectonic settings, being the activated fault at most about three times the hypocentre depth, explaining their higher b-value and the lower magnitude of the largest recorded events. Having different phenomenology, graviquakes show peculiar precursor

Ultra-thin clay layers facilitate seismic slip in carbonate faults

Many earthquakes propagate up to the Earth's surface producing surface ruptures. Seismic slip propagation is facilitated by along-fault low dynamic frictional resistance, which is controlled by a number of physico-chemical lubrication mechanisms. In particular, rotary shear experiments conducted at seismic slip rates (1 ms(-1)) show that phyllosilicates can facilitate co-seismic slip along faults during earthquakes. This evidence is crucial for hazard assessment along oceanic subduction zones, where pelagic clays participate in seismic slip propagation. Conversely, the reason why, in continental domains, co-seismic slip along faults can propagate up to the Earth's surface is still poorly understood. We document the occurrence of micrometer-thick phyllosilicate-bearing layers along a carbonate-hosted seismogenic extensional fault in the central Apennines, Italy. Using friction experiments, we demonstrate that, at seismic slip rates (1 ms(-1)), similar calcite gouges with pre-existing phyllosilicate-bearing (clay content ≤3 wt.%) micro-layers weaken faster than calcite gouges or mixed calcite-phyllosilicate gouges. We thus propose that, within calcite gouge, ultra-low clay content (≤3 wt.%) localized along micrometer-thick layers can facilitate seismic slip propagation during earthquakes in continental domains, possibly enhancing surface displacement

Tectonically asymmetric Earth: From net rotation to polarized westward drift of the lithosphere

AbstractThe possibility of a net rotation of the lithosphere with respect to the mantle is generally overlooked since it depends on the adopted mantle reference frames, which are arbitrary. We review the geological and geophysical signatures of plate boundaries, and show that they are markedly asymmetric worldwide. Then we compare available reference frames of plate motions relative to the mantle and discuss which is at best able to fit global tectonic data. Different assumptions about the depths of hotspot sources (below or within the asthenosphere, which decouples the lithosphere from the deep mantle) predict different rates of net rotation of the lithosphere relative to the mantle. The widely used no-net-rotation (NNR) reference frame, and low (1°/Ma) net rotation (shallow hotspots source), all plates, albeit at different velocity, move westerly along a curved trajectory, with a tectonic equator tilted about 30° relative to the geographic equator. This is consistent with the observed global tectonic asymmetries

Complex geometry and kinematics of subsidiary faults within a carbonate-hosted relay ramp

Minor fault geometry and kinematics within relay ramps is strongly related to the stress field perturbations that can be produced when two major fault segments overlap and interact. Here we integrate classical fieldwork and interpretation of a virtual outcrop to investigate the geometry and kinematics of subsidiary faults within a relay ramp along the Tre Monti normal fault in the Central Apennines. Although the Tre Monti fault strikes parallel to the regional extension (NE-SW) it shows predominant dip-slip kinematics, suggesting a NW-SE oriented extension acting at sub-regional scale (1–10 km). Conversely, the slickenlines collected on the front segment of the relay ramp highlight right-lateral kinematics. The subsidiary faults in the relay ramp show a complex geometry (variable attitudes) and slickenlines describe multiple kinematics (left-lateral, dip-slip, right-lateral), independently of their orientation. Our fault slip analysis indicates that a local stress field retrieved from the kinematic inversion of the slickenlines collected on the front segment, and likely promoted by the interaction between the overlapping fault segments that bound the relay zone, can explain most of the geometry and kinematics of the subsidiary faults. Further complexity is added by the temporal interaction with both the regional and sub-regional stress fields

Ground deformation and source geometry of the 30 October 2016 Mw 6.5 Norcia earthquake (Central Italy) investigated through seismological data, DInSAR measurements, and numerical modelling

We investigate the Mw 6.5 Norcia (Central Italy) earthquake by exploiting seismological data, DInSAR measurements, and a numerical modelling approach. In particular, we first retrieve the vertical component (uplift and subsidence) of the displacements affecting the hangingwall and the footwall blocks of the seismogenic faults identified, at depth, through the hypocenters distribution analysis. To do this, we combine the DInSAR measurements obtained from coseismic SAR data pairs collected by the ALOS-2 sensor from ascending and descending orbits. The achieved vertical deformation map displays three main deformation patterns: (i) a major subsidence that reaches the maximum value of about 98 cm near the epicentral zones nearby the town of Norcia; (ii) two smaller uplift lobes that affect both the hangingwall (reaching maximum values of about 14 cm) and the footwall blocks (reaching maximum values of about 10 cm). Starting from this evidence, we compute the rock volumes affected by uplift and subsidence phenomena, highlighting that those involved by the retrieved subsidence are characterized by significantly higher deformation values than those affected by uplift (about 14 times). In order to provide a possible interpretation of this volumetric asymmetry, we extend our analysis by applying a 2D numerical modelling approach based on the finite element method, implemented in a structural-mechanic framework, and exploiting the available geological and seismological data, and the ground deformation measurements retrieved from the multi-orbit ALOS-2 DInSAR analysis. In this case, we consider two different scenarios: the first one based on a single SW-dipping fault, the latter on a main SW-dipping fault and an antithetic zone. In this context, the model characterized by the occurrence of an antithetic zone presents the retrieved best fit coseismic surface deformation pattern. This result allows us to interpret the subsidence and uplift phenomena caused by the Mw 6.5 Norcia earthquake as the result of the gravitational sliding of the hangingwall along the main fault plane and the frictional force acting in the opposite direction, consistently with the double couple fault plane mechanism

Numerical analysis of interseismic, coseismic and postseismic phases for normal and reverse faulting earthquakes in Italy

The preparation, initiation, and occurrence dynamics of earthquakes in Italy are governed by several frequently unknown physical mechanisms and parameters. Understanding these mechanisms is crucial for developing new techniques and approaches for earthquake monitoring and hazard assessments. Here, we develop a first-order numerical model simulating quasi-static crustal interseismic loading, coseismic brittle episodic dislocations, and postseismic relaxation for extensional and compressional earthquakes in Italy based on a common framework of lithostatic and tectonic forces. Our model includes an upper crust, where the fault is locked, and a deep crust, where the fault experiences steady shear. The results indicate that during the interseismic phase, the contrasting behavior between the upper locked fault segment and lower creeping fault segment generates a stretched volume at depth in the hanging wall via extensional tectonics while a contracted volume forms via compressional tectonics. The interseismic stress and strain gradients invert at the coseismic stage, with the interseismic dilated volume contracting during the coseismic stage, and vice versa. Moreover, interseismic stress gradients promote coseismic gravitational subsidence of the hanging wall for normal fault earthquakes and elastic uplift for reverse fault earthquakes. Finally, the postseismic relaxation is characterized by further ground subsidence and uplift for normal and reverse faulting earthquakes, respectively, which is consistent with the faulting style. The fault is the passive feature, with slipping generating the seismic waves, whereas the energy activating the movement is stored mostly in the hanging wall volume. The main source of energy for normal faulting and thrust is provided by the lithostatic load and elastic load, respectively

The tectonic puzzle of the Messina area (Southern Italy): Insights from new seismic reflection data

The Messina Strait, that separates peninsular Italy from Sicily, is one of the most seismically active areas of the Mediterranean. The structure and seismotectonic setting of the region are poorly understood, although the area is highly populated and important infrastructures are planned there. New seismic reflection data have identified a number of faults, as well as a crustal scale NE-trending anticline few km north of the strait. These features are interpreted as due to active right-lateral transpression along the north-eastern Sicilian offshore, coexisting with extensional and right-lateral transtensional tectonics in the southern Messina Strait. This complex tectonic network appears to be controlled by independent and overlapping tectonic settings, due to the presence of a diffuse transfer zone between the SE-ward retreating Calabria subduction zone relative to slab advance in the western Sicilian side

Are normal fault earthquakes due to elastic rebound or gravitational collapse?

n/

On the geodynamics of the northern Adriatic plate

The northern Adriatic plate underwent Permian-Mesozoic rifting and was later shortened by three orogenic belts (i.e., Apennines, Alps and Dinarides) developed along three independent subduction zones. The inherited Mesozoic horst and graben grain determined structural undulations of the three thrust belts. Salients developed in grabens or more shaly basins, whereas recesses formed regularly around horsts. A new interpretation of seismic reflection profiles, subsidence rates from stratigraphic analysis, and GPS data prove that the three orogens surrounding the northern Adriatic plate are still active. The NE-ward migration of the Apennines subduction hinge determines the present-day faster subsidence rate in the western side of the northern Adriatic (> 1 mm/year). This is recorded also by the SW-ward dip of the foreland regional monocline, and the SW-ward increase of the depth of the Tyrrhenian sedimentary layer, as well as the increase in thickness of the Pliocene and Pleistocene sediments. These data indicate the dominant influence of the Apennines subduction, which controls the asymmetric subsidence in the northern Adriatic realm. The Dinarides front has been tilted by the Apennines subduction hinge, as shown by the eroded Dalmatian anticlines subsiding in the eastern Adriatic Sea. GPS data suggest that southward tilting of the western and central Southern Alps, whereas the eastern Southern Alps are uplifting. The obtained strain rates are on average within 20 nstrain/year. The horizontal shortening obtained from GPS velocities at the front of the three belts surrounding the northern Adriatic plate are about 2-3 mm/year (Northern Apennines), 1-2 mm/year (Southern Alps), and < 1 mm/year (Dinarides). The shortening directions tend to be perpendicular to the thrust belt fronts. The areas where the strain rate sharply decreases along a tectonic feature (e.g., the Ferrara salient, the Venetian foothills front) are proposed to be occupied by locked structures where stress is accumulating in the brittle layer and thus seismically prone. Finally, we speculate that, since the effects of three independent subduction zones coexist and overlap in the same area, plate boundaries are passive features