Earthquake Doublet Research Articles

Persistent asperities at the Kermadec subduction zone controlled by changes in forearc structure: 1976 and 2021 doublet earthquakes

To what degree large earthquakes are ‘characteristic’ is a critical question in understanding the fundamentals of earthquake physics and the seismic cycle. Here we study the ruptures of the 2021 Mw7.4 and Mw8.1 doublet earthquake sequence in the Kermadec subduction zone and compare them with the 1976 doublet that occurred at the same location. We find that although the 2021 mainshock likely re-ruptured the same asperity as the 1976 Mw7.9 event, the detailed slip distribution is different. Other ruptures in the doublets also differ in character and location. Our observations indicate the variability between large earthquakes on the same segment of the plate boundary in each earthquake cycle. All earthquakes occur in an isolated area of the megathrust, which is bounded by changes in the lithospheric structure of the overriding plate as indicated by bathymetric and gravity data. This high-seismicity region is coincident with an isolated forearc sedimentary basin, possibly formed by basal erosion related to seismogenesis, suggesting that seismic slip has persisted here for several million years. Refined up-dip aftershock and background seismicity focal mechanisms have a steeper dip angle than the slab interface, suggesting these events are located within the subducting oceanic slab, possibly forming a rougher plate interface that facilitates basal erosion. We conclude that the stress heterogeneity within this bounded seismogenic zone is long-lived and has produced a rich spectrum of earthquake ruptures.

Lithospheric strength of the Anatolian plateau and implications for strong earthquake activity in Turkey

On February 6, 2023, the doublet earthquake including two main shocks with magnitudes MW7.8 and MW7.5, occurred near the western side of the East Anatolian Fault at the southeast boundary of the Anatolian Plateau in Turkey. Based on the WGM2012 Bouguer gravity anomaly data and the Etopo1 topography data, this study first introduced a joint inversion of admittance and coherence functions and used the Bayesian optimal parameter estimation method to obtain the effective elastic thickness Te and loading ratio F of the lithosphere for various tectonic units in the Anatolian Plateau. Secondly, we discussed the characteristics and influencing factors of the lithospheric mechanical strength and analyzed its relationship with seismic activity. The lithospheric mechanical strength of the Anatolian Plateau showed clear lateral heterogeneity and a "weak-strong-weak" spatial pattern from east to west, reflecting various tectonic processes. At last, the strong seismic activity was found where the lithospheric strength was low in the Anatolian Plate. We also incorporated GPS strain rate and other results to investigate the tectonic background and primary causes of the MW7.8 and MW7.5 doublet earthquakes in Turkey. The results have a good insight into urban safety design in the Turkish region, including post-disaster rehabilitation, earthquake hazard assessment, and loss reduction.

Effect of Continuous Glucose Monitoring Device Assistance on Glycemic Control of 2023 Kahramanmaraş Doublet Earthquake Survivors with Type 1 Diabetes in Adana, Turkey.

The aim of this study was to investigate the effects of the devastating 2023 Kahramanmaraş earthquake on the glycemic control of children with type 1 diabetes (T1DM) in Adana, Turkey. Additionally, the study aimed to assess the impact of continuous glucose monitoring (CGM) device assistance on glycemic control after the earthquake. A retrospective study was conducted involving 134 children with T1DM receiving intensive insulin treatment. Participants were divided into 2 groups: CGM (+) (n = 58), who benefited from CGM device assistance, and CGM (-) (n = 76), who did not utilize CGM device after the earthquake. Glycated hemoglobin (HbA1c) levels were recorded before and after the earthquake. Following the earthquake, the median HbA1c for all participants changed insignificantly from 8.9% to 8.6% (P = .491). However, in the CGM (+) group, HbA1c levels significantly improved post earthquake (P = .001). Conversely, the CGM (-) group experienced a deterioration in glycemic control (P = .027). A 2-way repeated measures ANOVA revealed a significant interaction effect between CGM device usage and the earthquake on HbA1c levels (F = 17.257, P <.001). Subgroup analysis based on age indicated that the effectiveness of CGM was more pronounced in adolescents (≥12 years) than in younger children (<12 years). This study highlights the adverse impact of the earthquake on glycemic control in children with T1DM and underscores the effectiveness of CGM in improving glycemic control, particularly among adolescents. The provision of CGM devices following the earthquake led to enhanced outcomes, mitigating the negative effects of the disaster on glycemic control.

Decoupled co-seismic deformation and stress changes during the 2021 (Mw 6.0, 6.4) North Bandar Abbas doublet earthquakes in Fin region, SE-syntaxis of Zagros, Iran: new insights into the tectonic deformation decoupling process

Decoupled co-seismic deformation and stress changes during the 2021 (Mw 6.0, 6.4) North Bandar Abbas doublet earthquakes in Fin region, SE-syntaxis of Zagros, Iran: new insights into the tectonic deformation decoupling process

Directionality and polarization of response spectral ordinates in the 2023 Kahramanmaras, Türkiye earthquake doublet

Until recently, the orientation of maximum horizontal spectral response was generally believed to not have a predominant orientation at rupture distances greater than 5 km. However, a recent study found that the orientation of maximum spectral response for strike-slip earthquakes in the NGA-West2 database tends to occur close to the epicentral transverse orientation, that is, an orientation perpendicular to a line connecting the epicenter to the station. This article investigates directionality in the 6 February 2023 Türkiye doublet earthquakes ( Mw 7.8 and 7.5) with strike-slip faulting. The orientation of the maximum response of 5%-damped linear elastic oscillators was studied. The spatial distribution of the level of polarization, which in this article refers to the amount of directionality, and intensities at specific orientations were also studied. The maximum spectral response was found to occur systematically close to the epicentral transverse orientation, consistent with previous observations for other strike-slip earthquakes. For the Mw 7.8 event where the location of maximum slip was relatively far from the epicenter, it was found that the orientation of maximum spectral response is, on average, closer to the maximum slip transverse orientation (i.e. perpendicular to a line connecting the station to the surface projection of the point of maximum slip) when compared to the epicentral transverse orientation over most period ranges. This suggests that the maximum slip transverse orientation may be a better estimator for determining the orientation of maximum spectral response in large-magnitude strike-slip earthquakes, although further study using more events is warranted. Polarized motions were observed over large geographical areas, and the orientation of maximum spectral response was found to be close to the epicentral or maximum slip transverse for Joyner–Boore distances up to the farthest studied (400 km). These findings further support the case for the development of orientation-dependent ground motion models for strike-slip earthquakes.

Geometric controls on cascading rupture of the 2023 Kahramanmaraş earthquake doublet

Geometric controls on cascading rupture of the 2023 Kahramanmaraş earthquake doublet

3D Dynamic Rupture Modeling of the 6 February 2023, Kahramanmaraş, Turkey Mw 7.8 and 7.7 Earthquake Doublet Using Early Observations

Abstract The 2023 Turkey earthquake sequence involved unexpected ruptures across numerous fault segments. We present 3D dynamic rupture simulations to illuminate the complex dynamics of the earthquake doublet. Our models are constrained by observations available within days of the sequence and deliver timely, mechanically consistent explanations of the unforeseen rupture paths, diverse rupture speeds, multiple slip episodes, heterogeneous fault offsets, locally strong shaking, and fault system interactions. Our simulations link both earthquakes, matching geodetic and seismic observations and reconciling regional seismotectonics, rupture dynamics, and ground motions of a fault system represented by 10 curved dipping segments and embedded in a heterogeneous stress field. The Mw 7.8 earthquake features delayed backward branching from a steeply branching splay fault, not requiring supershear speeds. The asymmetrical dynamics of the distinct, bilateral Mw 7.7 earthquake are explained by heterogeneous fault strength, prestress orientation, fracture energy, and static stress changes from the previous earthquake. Our models explain the northward deviation of its eastern rupture and the minimal slip observed on the Sürgü fault. 3D dynamic rupture scenarios can elucidate unexpected observations shortly after major earthquakes, providing timely insights for data-driven analysis and hazard assessment toward a comprehensive, physically consistent understanding of the mechanics of multifault systems.

Complex multi-fault rupture and triggering during the 2023 earthquake doublet in southeastern Türkiye

Two major earthquakes (MW 7.8 and MW 7.7) ruptured left-lateral strike-slip faults of the East Anatolian Fault Zone (EAFZ) on February 6, 2023, causing >59,000 fatalities and ~$119B in damage in southeastern Türkiye and northwestern Syria. Here we derived kinematic rupture models for the two events by inverting extensive seismic and geodetic observations using complex 5-6 segment fault models constrained by satellite observations and relocated aftershocks. The larger event nucleated on a splay fault, and then propagated bilaterally ~350 km along the main EAFZ strand. The rupture speed varied from 2.5-4.5 km/s, and peak slip was ~8.1 m. 9-h later, the second event ruptured ~160 km along the curved northern EAFZ strand, with early bilateral supershear rupture velocity (>4 km/s) followed by a slower rupture speed (~3 km/s). Coulomb Failure stress increase imparted by the first event indicates plausible triggering of the doublet aftershock, along with loading of neighboring faults.

Ground motion characteristics of the 2015 central Nepal (Gorkha-Dolakha) Earthquake and 2023 SE Turkey Earthquake doublets

The Alpine-Himalayan orogenic belt, the well-identified seismically active region in the world, has been evidenced by the occurrence of recent two devastating earthquake doublets in 2015 in central Nepal (MW 7.8 and MW 7.3) and in 2023 in southeastern Turkey(MW 7.8 and MW 7.5). In this study, seismic ground motion characteristics of two damaging earthquakes namely the 2015 central Nepal earthquake doublet that occurred in the central Himalaya and the 2023 southeastern Turkey earthquake doublet that occurred in the eastern Mediterranean has been investigated. The observed strong-motion recordings indicate that the peak ground acceleration (PGA) values are about 152.61–154.96 cm/s2 of the 2015 central Nepal earthquake doublet at KATNP station and about 451.834–842.924 cm/s2 of the 2023 southeastern Turkey earthquake doublet at station 3137. The seismic ground motion of both the 2015 central Nepal and the 2023 SE Turkey earthquakes doublet shows that the 2015 central Nepal doublet produced much lower peak ground acceleration than the 2023 southeastern Turkey earthquake doublet.

Satellite Imagery for Rapid Detection of Liquefaction Surface Manifestations: The Case Study of Türkiye–Syria 2023 Earthquakes

The 6 February 2023 earthquake doublet (Mw 7.7 and Mw 7.6) that occurred on the East Anatolian Fault Zone (EAFZ) triggered a significant amount of soil liquefaction phenomena in SE Türkiye and NW Syria. The great areal extent of the affected area and the necessity of rapid response led to the adoption and improvement of a workflow for mapping liquefaction phenomena based on remote sensing data. Using satellite imagery, we identified 1850 sites with liquefaction manifestation and lateral spreading deformation. We acquired a thorough map of earthquake-triggered liquefaction based on visual mapping with optical satellite imagery (high and very high-resolution) and the aid of radar satellite imagery and interferometry. The majority of sites are found along meandering sections of river valleys, coastal plains, drained lakes, swamps, and lacustrine basins along the East Anatolian Fault, highlighting once again the influence of geomorphology/surficial geology on the distribution of liquefaction phenomena. A total of 95% of the liquefaction occurrences were mapped within 25 km from the surface trace of the fault, confirming the distance from fault rupture as a more effective tool for predicting the distribution of liquefaction than epicentral distance. Thus, taking into consideration the rapid documentation of these phenomena without the limitations in terms of time, cost, and accessibility of the field investigation techniques, this desktop-based approach can result in a rapid and comprehensive map of liquefaction from a strong earthquake, and can also be used as a future guide for subsequent field investigations for liquefaction hazard mapping.

The 2023 Mw 7.8 and 7.6 Earthquake Doublet in Southeast Türkiye: Coseismic and Early Postseismic Deformation, Faulting Model, and Potential Seismic Hazard

Abstract On 6 February 2023, Mw 7.8 and 7.6 earthquakes struck southeast Türkiye and northwest Syria. They are the largest earthquakes in Türkiye in over 80 yr, causing significant damage and fatalities. We used Advanced Land Observation Satellite-2 and Sentinel-1 Synthetic Aperture Radar images to obtain near-field coseismic displacements by Differential Interferometric Synthetic Aperture Radar (DInSAR) and pixel offset tracking (POT). Discontinuities of the surface deformation suggest that the Mw 7.8 event ruptured ∼320 km along the East Anatolian fault, and the Mw 7.6 event ruptured ∼150 km near Elbistan, southern Türkiye. We inverted these earthquakes' fault geometry and slip distribution based on Global Navigation Satellite Systems, DInSAR, and POT displacements. The estimated fault slip model shows that the first Mw 7.8 event ruptured a steeply southeast-dipping fault, and the seismogenic fault of the second Mw 7.6 event is north-dipping with complex geometry. The dip angle of subfaults of the Mw 7.6 earthquake decreases from east to west. Faults responsible for the two earthquakes are dominated by left-lateral strike-slip motion, with the maximum slip of ∼9.1 m. Early postseismic deformation within two months exhibits displacement discontinuities in the Amanos and Pazarcık segments and the Çardak fault, suggesting that the afterslip partially compensated the coseismic slip deficit at the shallow depths. Furthermore, static Coulomb failure stress changes induced by the two earthquakes indicate that the southwestern Pütürge segment of the East Anatolian fault has a high risk of future rupture.

Comment on “Inner Core Rotation Captured by Earthquake Doublets and Twin Stations” by Yang and Song

AbstractYang &amp; Song (2022, https://doi.org/10.1029/2022GL098393) first claimed existence of Earth's inner core differential rotation based on the waveform similarity of two neighboring stations AAK and KZA across an earthquake doublet and then postulated a local velocity gradient at the top of the inner core based on the difference of PKiKP‐PKIKP differential times between the stations and inferred inner core differential rotation rate. In this comment, we collectively analyze the seismic data in the region and add the data of another nearby station HORS into analysis. HORS and KZA, located in an opposite direction away from AAK, consistently exhibit high waveform similarity. Collective analysis of seismic data demonstrates the invalidity of both their logic of claiming existence of inner core differential rotation and their postulation of “a local inner core gradient” to infer differential rotation. Localized and episodic inner core surface change provides a physically consistent explanation to the seismic data.

Effect on RC buildings of 6 February 2023 Turkey earthquake doublets and new doctrines for seismic design

Two major catastrophic earthquakes, which occurred 9 h apart on the Eastern Anatolian Fault Zone and one of its branches on February 6, 2023 in Turkey, directly affected 11 provinces in the Eastern and South-Eastern Anatolia regions, where 14 million people live. These earthquakes are among the most destructive earthquakes in the history of Turkey, caused approximately 3 times more property damage and loss of life than the 1999 Marmara Earthquake. More than 24,000 aftershocks with magnitudes of up to 6.7 Mw occurred after these earthquakes. In this study, the causes of heavy destruction in some of the provinces and districts most affected by earthquakes (especially Hatay-Antakya, Kahramanmaraş, Gaziantep-Nurdagi, Gaziantep-Islahiye, Adıyaman-Gölbaşı, Malatya) were examined under sub-headings. The reasons for the damages were supported by the striking images obtained during the site investigation in the earthquake zone, and the issues that should be considered based on both seismic codes and implementations were also stated. In addition, in order to better understand the effects of earthquakes on reinforced concrete structures, the ratios of the loads acting on the structures in both earthquakes to the design load predicted by the code valid at the time the structures were built are presented in graphics. The results obtained show that the effects of design and application mistakes are quite high in the heavy destruction caused by earthquakes. But it is understood that the seismic code design criteria and requirements are insufficient in some regions. In addition, the establishment of residential areas in risky areas without taking any precautions has produced dramatic results.

The complex dynamics of the 2023 Kahramanmaraş, Turkey, Mw 7.8-7.7 earthquake doublet.

The destructive 2023 moment magnitude (Mw) 7.8-7.7 earthquake doublet ruptured multiple segments of the East Anatolian Fault system in Turkey. We integrated multiscale seismic and space-geodetic observations with multifault kinematic inversions and dynamic rupture modeling to unravel the events' complex rupture history and stress-mediated fault interactions. Our analysis reveals three subshear slip episodes during the initial Mw 7.8 earthquake with a delayed rupture initiation to the southwest. The Mw 7.7 event occurred 9 hours later with a larger slip and supershear rupture on its western branch. Mechanically consistent dynamic models accounting for fault interactions can explain the unexpected rupture paths and require a heterogeneous background stress. Our results highlight the importance of combining near- and far-field observations with data-driven and physics-based models for seismic hazard assessment.

Investigating the aftermath of the Türkiye 2023 earthquake: exploring post-disaster uncertainty among Syrian migrants using social network analysis with public health approach.

On February 6th, 2023, a doublet earthquake struck Türkiye, impacting more than 15 million people including migrants, and resulting in over 50,000 deaths. The Syrian migrants experience multiple uncertainties in their daily lives which are further compounded by multifaceted challenges of the post-disaster environment. Social media was used intensively and with impunity in this environment and thereby provides a window into the explicit and implicit dynamics of daily life after a disaster. We aimed to explore how a post-disaster environment potentially generates new uncertainties or exacerbating pre-existing ones for migrants through social media analysis with an indirect perspective, in the context of 2023-Earthquake in Türkiye and Syrian migrants. Social network analysis was used to analyze Twitter-data with the hashtags 'Syrian' and 'earthquake' during a 10-day period beginning on March 22nd, 2023. We calculated network metrics, including degree-values and betweenness-centrality and clustered the network to understand groups. We analyzed a combination of 27 tweets with summative content analysis using a text analysis tool, to identify the most frequently used words. We identified the main points of each tweet and assessed these as possible contributors to post-disaster uncertainty among migrants by using inductive reasoning. There were 1918 Twitter users, 274 tweets, 124 replies and 1726 mentions. Discussions about Syrian migrants and earthquakes were established across various groups (ngroups(edges > 15) = 16). Certain users had a greater influence on the overall network. The nine most frequently used words were included under uncertainty-related category (nmost_frequently_used_words = 20); 'aid, vote, house, citizen, Afghan, illegal, children, border, and leave'. Nine main points were identified as possible post-disaster uncertainties among migrants. The post-disaster environment has the potential to exacerbate existing uncertainties, such as being an undocumented migrant, concerns about deportation and housing, being or having a child, inequality of rights between being a citizen and non-citizen, being in minority within minority, political climate of the host nation and access to education or to generate new ones such equitable distribution of aid, which can lead to poor health outcomes. Recognizing the possible post-disaster uncertainties among migrants and addressing probable underlying factors might help to build more resilient and healthy communities.

Rapid source inversions of the 2023 SE Türkiye earthquakes with teleseismic and strong-motion data

We conducted rapid inversions of rupture process for the 2023 earthquake doublet occurred in SE Türkiye, the first with a magnitude of MW7.8 and the second with a magnitude of MW7.6, using teleseismic and strong-motion data. The teleseismic rupture models of the both events were obtained approximately 88 and 55 minutes after their occurrences, respectively. The rupture models indicated that the first event was an asymmetric bilateral event with ruptures mainly propagating to the northeast, while the second one was a unilateral event with ruptures propagating to the west. This information could be useful in locating the meizoseismal areas. Compared with teleseismic models, the strong-motion models showed relatively higher resolution. A noticeable difference was found for the MW7.6 earthquake, for which the strong-motion models shows a bilateral event, rather than a unilateral event, but the dominant rupture direction is still westward. Nevertheless, all strong-motion models are consistent with the teleseismic models in terms of magnitudes, durations, and dominant rupture directions. This suggests that both teleseismic and strong-motion data can be used for fast determination of major source characteristics. In contrast, the strong-motion data would be preferable in future emergency responses since they are recorded earlier and have a better resolution ability on the source ruptures.

Coseismic Faulting Model and Post-Seismic Surface Motion of the 2023 Turkey–Syria Earthquake Doublet Revealed by InSAR and GPS Measurements

On 6 February 2023 (UTC), an earthquake doublet, consisting of the Mw 7.8 Pazarcik earthquake and the Mw 7.5 Elbistan earthquake, struck south-central Turkey and northwestern Syria, which was the largest earthquake that occurred in Turkey since the 1939 Erzincan earthquake. The faulting model of this earthquake was estimated based on the coseismic InSAR and GPS displacements. In addition, the best-fitting coseismic faulting model indicates that both the Pazarcik earthquake and the Elbistan earthquake were controlled by predominated left-lateral strike-slip motion, with slip peaks of 9.7 m and 10.8 m, respectively. The Coulomb failure stress (CFS) change suggests that the Pazarcik earthquake has a positive effect in triggering the rupture of the seismogenic fault of the Elbistan earthquake. Furthermore, these two main shocks promoted the occurrence of the Mw 6.3 strong aftershock. Additionally, it is found that the 2023 Turkey-Syria earthquake doublet increased the rupture risk of the Puturge segment of the EAF fault and the northern segment of the Dead Sea Fault (DSF). It is crucial to note that the northern segment of the DSF has not experienced a large earthquake in several centuries, highlighting the need for heightened attention to the potential seismic hazard of this segment. Finally, a deformation zone adjacent to the DSF was identified, potentially attributed to the motion of a blind submarine fault.

Multi‐Scale Rupture Growth With Alternating Directions in a Complex Fault Network During the 2023 South‐Eastern Türkiye and Syria Earthquake Doublet

AbstractA devastating doublet of earthquakes with moment magnitude MW 7.9 and MW 7.6 earthquakes contiguously occurred in SE Türkiye near the NW border of Syria. Here we perform a potency‐density tensor inversion to simultaneously estimate rupture evolution and fault geometry for the doublet. We find the initial MW 7.9 earthquake involved discrete episodes of supershear rupture and back‐rupture propagation, and was triggered by initial rupture along a bifurcated splay of the East Anatolian Fault. The second MW 7.6 event was triggered by the earlier MW 7.9 event, and it involved more extensive supershear rupture along a favorably curved fault, and was likely stopped by geometric barriers at the fault ends. Our results highlight the multi‐scale cascading rupture growth across the complex fault network that affects the diverse rupture geometries of the 2023 Türkiye earthquake doublet, contributing to the strong ground shaking and associated devastation.

Source mechanisms and triggering process for the April 12th and 13th, 2014 earthquake doublet in the Solomon Islands

Mw 7.6 and Mw 7.4 earthquakes occurred 16 h apart on April 12 and April 13, 2014, at depths of 15 km and 35 km, respectively, southwest of Makira Island (San Cristobal) in the Solomon Islands. We study the source mechanisms to investigate the interrelation between the two events, and to investigate why the Solomon Islands have a historically high rate of doublet earthquakes. Teleseismic P waveform data for both events were used to obtain the slip distributions for the two earthquakes by a finite fault slip inversion method. The Mw 7.6 event is an east–west left-lateral strike slip, where the maximum slip is observed 8–24 km east of the hypocenter. The Mw 7.4 earthquake is a thrust event with north and south-dipping nodal planes. The fault plane is not clear from the aftershock distribution, but the north-dipping plane considered to be the likely fault plane because of the simpler slip distribution compared to the south-dipping plane. The static Coulomb failure stress changes caused by the first earthquake were calculated in the region of the second earthquake. The results show that there was a + 48.59 kPa change at the hypocenter for the assumed north-dipping plane and + 18.24 kPa for the assumed south-dipping plane of the Mw 7.4 event. The temporal pattern of aftershocks shows a possible rate increase of stress prior to the occurrence of the Mw 7.4 events, which may also contribute to the triggering of the second event by static Coulomb Stress changes. We propose a model that supports the geological complexity of the region that may encourage such doublet events.Graphical

Contemporaneous Thick- and Thin-Skinned Seismotectonics in the External Zagros: The Case of the 2021 Fin Doublet, Iran

In this work, we propose a geodetic model for the seismic sequence, with doublet earthquakes, that occurred in Bandar Abbas, Iran, in November 2021. A dataset of Sentinel-1 images, processed using the InSAR (Interferometric Synthetic Aperture Radar) technique, was employed to identify the surface deformation caused by the major events of the sequence and to constrain their geometry and kinematics using seismological constraints. A Coulomb stress transfer analysis was also applied to investigate the sequence’s structural evolution in space and time. A linear inversion of the InSAR data provided a non-uniform distribution of slip over the fault planes. We also performed an accurate relocation of foreshocks and aftershocks recorded by locally established seismographs, thereby allowing us to determine the compressional tectonic stress regime affecting the crustal volume. Despite the very short time span of the sequence, our results clearly suggest that distinct blind structures that were previously unknown or only suspected were the causative faults. The first Mw 6.0 earthquake occurred on an NNE-dipping, intermediate-angle, reverse-oblique plane, while the Mw 6.4 earthquake occurred on almost horizontal or very low-angle (SSE-dipping) reverse segments with top-to-the-south kinematics. The former, which cut through and displaced the Pan-African pre-Palaeozoic basement, indicates a thick-skinned tectonic style, while the latter rupture(s), which occurred within the Palaeozoic–Cenozoic sedimentary succession and likely exploited the stratigraphic mechanical discontinuities, clearly depicts a thin-skinned style.