Tectonophysics最新文献

GNSS (Global Navigation Satellite System) data in northern Sulawesi and western Halmahera reveals a pattern of coseismic displacement that was caused by the 7 July 2019 (Mw 6.9) and 14 November 2019 (Mw 7.1) Molucca Sea earthquakes. The coseismic slip of these earthquakes are obtained via inversion on rectangular fault planes of surface GNSS coseismic deformation offsets. The 7 July 2019 earthquake ruptured on an east-dipping fault with a maximum slip of ∼35 cm located at ∼4 km depth and ∼ 100 km north-west of the epicenter. The 14 November 2019 earthquake also ruptured on an east-dipping fault, which has a maximum slip of ∼64 cm located at ∼22 km depth and ∼ 20 km south-west of the epicenter. The coseismic slip distribution of the 14 November earthquake is spatially aligned to that of an earthquake of similar magnitude that took place on 15 November 2014 in the same region. This observation points to the possibility of synchronization, thus raising the prospect of a future earthquake exceeding magnitude 7. If the stress was completely released during the 2014 event, it would have been necessary to reload that portion of the fault at a rate significantly larger than the observed convergence rate of ∼5.8 cm/yr. This suggests that partial ruptures are likely controlling the recurrence time of large earthquakes in the region.

We performed an analysis of present-day stress profiles for four wells penetrating the Lower Palaeozoic shale sequence of the Baltic Basin (northern Poland). Breakouts, hydraulic fracturing and leak-off tests were used to calibrate stress models based on anisotropic elastic shale properties. Initial stress models, balancing the lengths of the modelled and observed breakouts, indicated a degradation of the mechanical properties reconstructed from the density and dipole acoustic tool logs. Taking this adverse effect into account, final stress models were calculated which showed a stratification of the stress regime consistent with the lithostratigraphic shale units. A clear dominance of the horizontal stress-generating gravity factor over the horizontal tectonic strain was demonstrated. The obtained values of tectonic strain in the shale sequence compared to the previously determined strain in the crystalline basement of the same study area suggest a significant role of viscoelastic relaxation of the shale sequence.

The Haiyuan Arcuate Tectonic Belt (HATB) in the northeastern Tibetan Plateau features the interactions of three intersecting blocks: the eastern Qilian Shan, the Alxa Block, and the Ordos Block. While the HATB has displayed active responses to the ongoing collision between the Indian and Eurasian plates, the exact process behind the formation of this arcuate belt remains unclear. In pursuit of further insights into this topic and a deeper comprehension of the tectonic responses in NE Tibet, we conducted receiver function calculations using teleseismic waveforms recorded by two seismic short-period dense arrays spanning the western and eastern HATB, respectively, extending into the Alxa and Ordos Blocks. The CCP results in the HATB show major structural features that are different from those of adjacent blocks, mainly characterized by structural discontinuities in the crust due to severe deformation, including bending and uplifting in the lower crust. Together with previous geological studies, the bending interfaces in the lower crust of the HATB illuminate the existence of a crustal-scale tectonic accretionary wedge within the HATB, which originated in the Early Paleozoic. Furthermore, a decoupled deformation is seen within the HATB, with the lower crust undergoing shortening and the upper crust experiencing sequential stepwise thrusting towards the north. These scenarios, coupled with the resistance from the rigid Alxa and Ordos Blocks, lead to the conclusion that the arcuate shape of this belt is influenced by the weak crust of the HATB, which primarily orients the northeast, where the weak lithosphere of the Helan tectonic belt is situated between the Alxa and Ordos Blocks. Meanwhile, the progression of a series of thrusting faults in the upper crust within the HATB extends outward, involving adjacent blocks in plateau's growth.

The size-frequency distribution of many geological and geophysical variables, in relation to fractures, faulting and seismicity, is well described by a statistical distribution of the power law type which is characterized by its exponent. For earthquake magnitudes, the exponent is the well-known b-parameter of the Gutenberg-Richter scaling law. In this paper we:

• •

  provide a strict statistical derivation of the distribution law of earthquake magnitudes,

• •

  show that the maximum likelihood estimator of the b-parameter is unbiased,

• •

  demonstrate that the maximum likelihood estimator is invariant to the value chosen as the minimum magnitude threshold in so far as it is larger than the magnitude of completeness of the earthquake catalogue, and

• •

  provide a new estimator based on the minimization of the Kolmogorov-Smirnov statistic and provide a strategy for detecting and mapping the spatio-temporal variation of the b-parameter in seismic swarms.

The findings are illustrated with simulated data and a case study with real data.

The Three Gorges Reservoir, one of the largest water conservation system in the world, has been of keen interest to scientists globally since its impoundment. After construction of the dam, there has been a significant increase in seismic activity in the head area of the reservoir. It is generally accepted that earthquakes in this region are predominantly caused by the Jiuwanxi and Xiannvshan faults. This study focused on the stress changes occurring in the research area. A three-dimensional finite element model of the reservoir area was constructed using the geological structure and digital ground elevation data of the reservoir area. The fluid-solid coupling theory was applied to calculate the dynamic spatial changes in pore pressure and Coulomb stress in the faults and surrounding rocks during reservoir impoundment. The findings indicated that the added head pore water pressure at the bottom of the reservoir had a maximum impact range of approximately −2800 m on the surrounding rock, whereas the Xiannvshan and Jiuwanxi faults had a maximum diffusion range of approximately −4300 m. Rock permeability also played a significant role in the water storage process. During the 1 56 m water impoundment stage, owing to rapid water storage activity, stress could not be transmitted to both sides in a timely manner, resulting in the formation of an extreme stress change zone at −4000 m inside the fault. This may have been the reason for the frequent earthquakes during this stage. The 17 5 m cycle water storage stage also exhibited a significant degree of seismicity, potentially attributable to the long-term infiltration of reservoir water and accumulation of stress in the previous stage. The stress in the study area at the four stages are in a process of accumulation-release-accumulation-release.

The co-seismic and post-seismic deformation field associated with the M_w5.9 Arkalochori main shock that occurred in central Crete (Greece) on 27 September 2021 is analyzed using Copernicus Sentinel-1A & 1B images, GNSS measurements and seismological data. The fault geometry is constrained through the joint inversion of multiple datasets and the slip distribution for the co-seismic and post-seismic period is obtained using a homogeneous half-space elastic model and the Steepest Descent Method. The results indicate a blind normal fault striking 215° with a 55° dip to the northwest and the co-seismic slip model suggests a nearly circular main slip patch (8 × 6 km²) with a maximum slip of 0.98 m. Post-seismic displacements started rapidly after the main shock followed by a gradual decay as highlighted by the calculated InSAR time series. The temporal evolution of post-seismic slip is described by a simple logarithmic function, decaying faster at the southwest part of the fault. The cumulative afterslip model suggests that the maximum post-seismic slip of 0.23 m occurred within a similar depth range compared to the co-seismic one, yet with a shift towards the southwest. Post-seismic slip inside the main shock rupture area is sustained, highlighting the slow recovery of locking in the co-seismic slip region. Afterslip (seismic or aseismic) played a dominant role in the early post-seismic period acting complementarily to the main rupture. Indications suggest that the spatiotemporal evolution of the productive aftershock sequence may be driven afterslip, alongside other potential factors.

Although drilled samples of fault rocks have yielded information on frictional features of shallow subduction zones, the relationship of rupture propagation to the levels of friction and pore-fluid pressure remains uncertain. To investigate this topic, we performed dynamic rupture simulations along the megasplay fault that slipped during the 1944 M_w 8.0 Tonankai earthquake in the Nankai Trough. We used actual data from friction experiments on rocks from the fault segment and pre-existing pore pressures deduced from geophysical surveys for the shallow portion of 0–10 km depth along the fault. Simulations of low friction (friction coefficient ca. 0.04) produced large slip (about 30 m), whereas simulations using higher friction (friction coefficient ca. 0.2) suppressed the rupture. In simulations with low friction in which the pore-fluid pressure was nearly equal to the lithostatic stress, the slip decreased to about 25 m. However, when the simulations included slip-strengthening at shallow depth and higher friction, the slip still reached roughly 20 m. Such variability in slip during rupture propagation is caused by differences in the friction features and fluid pressure conditions of fault rocks, in which the friction features might be related to the mineral composition. Spatiotemporal heterogeneity in fault-rock type and their physical and hydraulic properties may fundamentally produce the complexity and variability of earthquake rupture propagation along the Nankai plate-subduction boundary.

Detachment faulting related to oceanic core complexes (OCCs) has been suggested to be a manifestation of slow seafloor spreading. Although numerical models suggest OCCs form under low magma supply, the specific interaction between magmatism and tectonic faulting remains elusive. This paper examines seismic observations detailing the spatiotemporal interactions between magmatism, high-angle faulting, and detachment faulting at a slow-spreading mid-ocean ridge in the West Philippine Basin. We identified a magma-rich spreading phase, indicated by a magmatic top basement and oceanic crust with shallow-penetrating high-angle normal faults. An axial valley reveals an along-strike transition from magmatically-dominated to highly tectonized oceanic crust over a distance of 70 km. Two older OCCs with concave-down fault geometries and a younger OCC with steep-dipping faulting suggest sequential detachments with the same polarity. Our findings suggest: (1) slow seafloor spreading alternates between high-angle faulting with a relatively high magma supply and detachment faulting with a limited magma supply; (2) sequential development of younger detachments in the footwall of its predecessor leads to an asymmetric split in the newly accreted crust; and (3) the life cycle of OCC ends with high-angle faults that overprint the detachment and act as magma pathways, sealing the OCC. Our study captures the dynamic interaction between high-angle and detachment faults and their concurrent and subsequent relationship to magmatic systems. This reveals that strain distribution along strike is critical to OCC formation, thus enriching our understanding beyond conventional considerations such as spreading rates and melt budgets at mid-ocean ridges.

Continental rifts normally initiate within previously deformed lithosphere and thus their evolution and architecture can be largely controlled by inherited weak zones in the pre-rift crust. Here, we quantify the role of the strength and obliquity of pre-existing crustal-scale weak zones in the evolution of continental rift systems. We use a 3D numerical geodynamic model to assess strain localization, associated fault development, and rift segmentation during the early stages of tectonic extension. We find that both the strength and obliquity of the weak zones significantly influence the patterns of strain localization. A pre-existing very weak zone with low obliquity can promote the development of continuous and long-lived border faults parallel to the rift axis. Conversely, a comparatively strong weak zone with high obliquity leads to a staggered en-echelon rift geometry that lacks rectilinear laterally persistent strain localization. Furthermore, we find that rift obliquity and weak zone strength may modulate rift fault length, throw, and azimuth. These results provide new and compelling insights into the structure and evolution of natural active rifts that develop within orogenic basement terranes.

Reactivation processes play a significative role in the localization of deformation but still remain hard to establish at the lithospheric scale. In this work, we built a 3D structural model, which enables to bridge the gap between the main tectonic structures observed at the surface and the geometry of the major interfaces (the Mohorovičić-discontinuity (hereafter Moho) and top of the basement) inferred from geophysical data acquired in the external Western Alps and their foreland. The geometry of these tectonic structures is interpreted in relation to their geodynamic evolution. The main results of this study highlight: (1) a strong contribution of thick-skinned Pyrenean-Provence and Alpine tectonics, (2) a lithospheric rooting of Variscan shear zones and related faults, and (3) the regional-scale influence of these inherited structures on the post-Paleozoic strain localization in the study area. Our 3D model shows that the pattern of Variscan shear zones that were developed at the end of the Paleozoic involved the whole crust, as emphasized by the Moho offsets. These shear zones were reactivated and localized Meso-Cenozoic deformation. The Variscan deformation pattern controlled the geometry of extensional basins, the propagation of Pyrenean-Provence deformation, and finally the Alpine deformation at crustal scale. Our 3D model shows minor crustal thickening (ca. 40 km) located below the Pelvoux External Crystalline Massif, which probably resulted from both Pyrenean and Alpine tectonic phases. In contrast, the southern part of the Alpine front shows a thinned crust (ca. 18 km) resulting from extensional Meso-Cenozoic phases between the Cévennes margin and the Durance basin.