On February 6, 2023, two destructive earthquakes struck southern and central Turkey and northern and western Syria. The epicenter of the first event(MW7.8)was 37km west-northwest of Gaziantep. The earthquake had a maximum Mercalli intensity of Ⅻ around the epicenter and in Antakya. It was followed by a MW7.7 earthquake nine hours later. This earthquake was centered 95km north-northeast from the first one. There was widespread damage and tens of thousands of fatalities. In response to these catastrophic events, in March 2023, a seismic scientific expedition led by China Earthquake Administration(CEA)was promptly organized to investigate the surface ruptures caused by these earthquakes. Here, we focus on the surface ruptures of the second earthquake, known as the Elbistan earthquake. The post-earthquake field survey revealed that the Elbistan earthquake occurred on the East Anatolian fault zone's northern branch(the Cardak Fault). This event resulted in forming a main surface rupture zone approximately 140km long and a secondary fault rupture zone approximately 20km long, which is nearly perpendicular to the main rupture.
We combined the interpretation of high-resolution satellite imagery and geomorphic investigations along the fault to determine the fault geometry and kinematics of the second earthquake event. The Elbistan earthquake formed a main surface rupture zone approximately 140km long, which strikes in an east-west direction along the Cardak Fault. The main rupture zone starts from Göksun in the west and extends predominantly eastward until the western end of the Sürgü Fault. It then propagates northeast along the southern segment of the Malatya fault zone. The entire Cardak Fault and the Malatya fault zone's southern segment are considered seismic structures for this earthquake. The overall surface rupture zone exhibits a linear and continuous distribution. Secondary ruptures show a combination of left-lateral strike-slip or left-lateral oblique-thrust deformation. Along the rupture zone, a series of en echelon fractures, moletracks, horizontal fault striations, and numerous displaced piercing markers, such as mountain ridges, wheat fields, terraces, fences, roads, and wheel ruts, indicate the predominance of pure left-lateral strike-slip motion for most sections. The maximum measured horizontal displacement is(7.6±0.3)m. According to the empirical relationship between the seismic moment magnitude of strike-slip faulting earthquakes and the length of surface rupture(SRL), a main rupture zone of 140km in length corresponds to a moment magnitude of approximately 7.6. Based on the relationship between the seismic moment magnitude and the maximum coseismic displacement, a maximum coseismic displacement of(7.6±0.3)m corresponds to a moment magnitude of about 7.5. The magnitudes derived from the two empirical relationships are essentially consistent, and they also agree with the moment magnitude provided by the USGS. Besides the main surface rupture zone, a secondary fault rupture zone extends nearly north-south direction for approximately 20km long. Unfortunately, due to the limited time and traffic problem, we did not visit this north-south-trending secondary fault rupture zone.
According to the summary of the history of earthquakes, it is evident that the main surface rupture zone has only recorded one earthquake in history, the 1544 MS6.8 earthquake, which indicates significantly less seismic activity compared to the main East Anatolian Fault. Moreover, the “earthquake doublet” will inevitably significantly impact the stress state and seismic hazard of other faults in the region. Seismic activity in this area remain at a relatively high level for years or even decades to come. The east-west striking fault, which has not been identified on the published active fault maps at the western end of the surface rupture zone, and the north-east striking Savrun Fault, which did not rupture this time, will experience destructive earthquakes in the future. It remains unknown why the east-west striking rupture did not propagate to the Sürgü Fault this time. More detailed paleoearthquake studies are needed to identify whether it is due to insufficient energy accumulation or because this section acts as a barrier. If the Sürgü Fault, about 40km long, was to rupture entirely in the future, the magnitude could reach 7 based on the empirical relationship.
Considering the distribution of historical earthquakes along the East Anatolian fault zone, as well as the geometric distribution of the surface ruptures from the recent “earthquake doublet” and the surrounding active faults, it is believed that the future earthquake hazards in the northeastern segment of the East Anatolian fault zone, the northern segment of the Dead Sea Fault, and the Malatya Fault deserve special attention.
The ML5.7 Xingwen earthquake on December 16, 2018 is very likely induced by shale gas hydraulic fracturing, which caused not only massive landslides and rock collapse, but also some casualties in the surrounding area, with the direct economic loss of about 50 million CNY. It is of great significance to determine the source rupturing process of such an induced earthquake with large magnitude.
Finite fault inversion is one of the commonly adopted methods to determine coseismic slip displacement distribution. For finite fault inversion, various data have different sensitivities to various aspects of the rupture process. The seismic data can provide the historical information about the earthquake rupture process because it contains the Doppler effect of the space-time rupture behavior on the fault. In comparison, the near-field geodetic data(such as InSAR and GPS)can constrain the fault parameters and the static slip distribution well because they contain the surface motion information. Therefore, the reliability of the inversion for the complex seismic rupture process can greatly be improved by combined use of seismicdata and InSAR data.
In this study, strong-motion seismic data recorded at 8 near-field stations are chosen and filtered by a band-pass of 0.15-0.60Hz. The same InSAR data used in Wang et al.(2022)is adopted in this joint study. For inversion, a sufficiently large potential fault plane of 15km long and 10km wide is chosen and divided into 15×10 subfaults. Finally, the rupture process is obtained by joint inversion of strong-motion seismic data and InSAR data. The results show that the earthquake is characterzied by a typical unilateral rupture with the rupturing direction nearly towards the north. The duration of the rupture process was 6s, and the energy release was mainly concentrated in the first 5s. The rupture process is segmented and can be divided into two stages. The first stage is distributed from 1-3s and is located in the range of 0~5km from the source; and the 2nd stage is distributed from 3-5s and is located between 6 and 8km from the source. The coseismic slip is mainly concentrated in areas shallower than 5km, with a peak slip of approximately 0.27m. This can be used to explain why the Xingwen earthquake with a magnitude of ML5.7 caused relatively serious damages.
Combined with the distribution of foreshocks and aftershocks, it can be seen that the foreshocks were mainly concentrated to the eastern edge of the major coseismicslip zone, which are close to some hydraulic fracturing wells. This suggests that these foreshocks occuring at the edge of the main rupture zone has a certain correlation with fluids, and the presence of fluids further leads to the fault weakening of the mainshock due to the increase of pore pressure and the decrease of effective compressive stress, which plays a triggering role in the occurrence of the Xingwen earthquake. The aftershocks are mainly distributed around the main slip zone, which are caused by after slips after the mainshock. The results from seismic inversion, InSAR inversion and joint inversion of the two data types reveal that the Xingwen earthquake is a northward unilateral rupture. The rupture propagation direction and coseismic slip distribution may be related to the physical property changes along the fault plane.
Compared with the two single inversion results, the joint inversion overcomes the influence of uneven distribution of seismic stations, improves the resolution of slip distribution, and produces results that are more consistent with the real physical process. The slip model obtained by joint inversion in this study can be helpful for further understanding the mechanisms of induced earthquake, the correlation between induced earthquake and geological structure, earthquake disaster assessment and post-earthquake disaster prevention and hazard mitigation.
The Tan-Lu fault zone is a huge seismic-tectonic belt in the eastern China. It can be generally divided into three segments: the north, the middle, and the south segment. Among them, recent activity of the middle segment has been most thoroughly studied. The junction section between Jiangsu and Anhui Province is located in the transition zone between the middle and the south segment of the fault zone. Due to the complex tectonic structure, unevenly distributed Quaternary deposits and severely transformed surface landscape, it is difficult to study the recent Quaternary activity of the fault. Research in recent years have shown that the faults in the Fushan and Ziyang areas to the south of the Huaihe River were active during late Pleistocene-early Holocene, and their activities were characterized by thrusting, normal faulting, tension and twisting. How is the fault activity extending southwards to Nüshan Lake and whether the late Quaternary activity occurred at Nüshan Lake are issues worthy of attention.
Geomorphology of the study area is characterized by slope plains and uplands. The uplands mostly extend in near north-south direction and are obviously controlled by the faults. In the remote sensing satellite images, linearity features of the fault from Huaihe River to Nüshan Lake are distinct. Field investigations confirmed that in the farmland to the east of Liugudui Village, north of Nüshan Lake, there are scarps extending in NNE direction and distributing intermittently due to faulting. In this study, we chose relatively clear scarps and excavated trenches across the fault. The trench revealed abundant faulting phenomena. The trench wall revealed a fault deformation zone as wide as 2~4 meters, consisting of 3 fault branches. Among them, faults f1 and f3 are the boundary faults while fault f2 is developed within the deformation zone. The latest activity of fault f3 on the west side has ruptured the overlying horizon of late Pleistocene strata, and the rupture extended upwards to the surface. OSL dating samples were collected in the uppermost layer of the faulted horizons. Dating results show that the fault has been active at least in late Pleistocene. The scratches and steps developed on the fault plane indicate that the fault has experienced thrusting and dextral faulting. The deformation zone appears dark brown, which is conspicuously different from the horizons on both sides. Materials in the fault zone are compacted, crumpled and deformed, and the alignment direction is consistent with the fault. The deformation zone contains gravels and calcium tuberculosis of different sizes. Two brownish-yellow clay masses in irregular shape are deposited near the upper part of the fault plane. Among them, the clay mass tk1 on the south wall of the trench is quite clear, with the upper part connected with f1 and the middle part obliquely cut by f2. OSL dating samples were collected from clay masses from two trench walls. The dating results are consistent with the late Pleistocene horizons, indicating that the brownish-yellow clay masses were involved in the fault zone when faulting occurred in the middle-late Pleistocene, and the faulting event occurred roughly between(50.92±4.65)kaBP and(27.12±2.26)kaBP. Our research shows that late Quaternary activity of the most active fault of the eastern branch of the Tanlu fault zone extended southwards to Nüshan Lake in Mingguang, but intensity of the fault activity has weakened.
The segment from Sihong in Jiangsu Province to Mingguang in Anhui Province is the structural node between the middle segment and the southern segment of the Tanlu fault zone. Trench exposures in Wangqian, Sunpaifang, Dahongshan in Sihong and Santang, Ziyang, Zhuliu in Mingguang and other places revealed a variety of faulting phenomena such as wedges, wedge-shaped mass, normal faulting, negative flower-shaped structure, clay mass, etc. These show that faults that were dominantly thrusting led to the local and abundant phenomena near surface in this region. The reasons for these different phenomena may be related to the influence of regional complex stresses and their changes on large-scale fault systems at different time and spaces scales.
The September 5, 2022, M6.8 Luding earthquake occurred along the southeastern segment of the Xianshuihe fault zone. Tectonics around the epicenter area is complicated and several faults had been recognized. Focal mechanisms of the main shock and inversions from earthquake data suggest that the earthquake occurred on a northwest-trending, steeply dipping strike-slip fault, which is consistent with the strike and slip of the Xianshuihe fault zone. We conducted a field investigation along the fault sections on both sides of the epicenter immediately after the earthquake. NW-trending fractures that were recognized as surface ruptures during the earthquake, and heavy landslides along the fault section between Ertaizi-Aiguocun village were observed during the field investigations. There are no surface ruptures developed along the fault sections north of the epicenter and south of Aiguocun village. Thus it can be concluded that there is a 15.5km-long surface rupture zone developed along the Moxi Fault(the section between Ertaizi and Aiguo village). The surface rupture zone trends northwest and shows a left-lateral strike slip, which is consistent with the strike and motion constrained by the focal mechanism. The coseismic displacements were measured to 20~30cm. Field observations, focal fault plane, distribution of the aftershocks, GNSS, and InSAR observation data suggest that the seismogenic structure associated with the M6.8 Luding earthquake is the Moxi Fault that belongs to the southeastern segment of the Xianshuihe fault zone. Slip along the segment south of the epicenter generated this earthquake, and also triggered slip along a northeast-trending fault and the northwestern section of the Moxi Fault in the epicenter. So, the M6.8 Luding earthquake is an event that is nucleated on the section south of the epicenter and then triggered an activity of the whole fault segment.
This paper studies the linear concentrated distribution of geomagnetic diurnal induced current and the spatial distribution characteristics of short-term in-situ recurrence anomalies in the 1~3 years before the 2016 Zadoi M6.2 earthquake, the 2017 Jiuzhaigou M7.0 earthquake and the 2017 Milin M6.9 earthquake. The main conclusions are as follows:(1)The overlapping segment anomalies occurring within 1~3 years before strong earthquakes usually have the phenomenon of seismic gaps. The overlapping segment gap is a large high-resistance body and also a hard body(hard inclusion)in the seismogenic model of hard body. Overlapping segment gap and seismic gap are the detection results of hard body with different depth distribution range by different physical detection methods. The distribution range of hard body is delineated by seismic gap in the upper and middle crust and overlapping segment gap in the middle-lower crust and upper mantle. The overlapping segment gaps occurred within 1~3 years before strong earthquakes, which are the anomalies in the third stage of the seismogenic model. The seismic gap before strong earthquakes has different stages. For example, background gaps are formed decades of years before strong earthquakes, and the gaps formed again about 1 year before the earthquakes.(2)The overlapping segment anomalies occurring within 1~3 years before strong earthquakes reflect the formation of short-time high conductive current channels in the high conductive layers among high resistance bodies. These short-time high conductive current channels are caused by the mutually detached slip events with up-arching property among the high-resistivity bodies located in the middle-lower crust and upper mantle, resulting from the upwelling of deep thermal fluid. They are the events in which the energy in the middle-lower crust and upper mantle migrates to the hard body in the seismogenic model of hard body, while the seismic gap events are the ones in which the energy in the middle-upper crust migrates to the hard body before the earthquake.(3)Based on the results of seismic high-pressure fluid experiments in recent years, and combined with the mechanism of overlapping segment seismic anomalies, it is considered that each sealed high-pressure fluid in the seismogenic fault of the source body will not rupture at the same time in the impending earthquake stage. The original free water in the fault, the sealed high-pressure fluid broken in the earlier stages, and the high-pressure thermal fluid upwelling into the fault in the deep may make the seismogenic fault of the source body full of free water, and may form a high-conductivity current channel in the fault with impending earthquake significance. The high-conductivity current channel may be a real impending earthquake anomaly. Obviously, it is found that the high-conductivity current channel in the fault in the impending earthquake stage has practical significance for the short-term and impending earthquake prediction.(4)The detachment slip events detected from the overlapping segment anomalies are located below the strong earthquake source, which is similar to the phenomenon that slow earthquake zone is located below the earthquake zone. Although the relationship between slow earthquakes and earthquake above them is unclear, some scholars believe that slip events produce stress accumulation on the surface of locked plates. A slip event may trigger a destructive earthquake, that is, a high-incidence period of intermittent tremors and slips can produce a peak period of seismicity. The above views on slow earthquakes are similar to the relationship between the linear overlapping segment anomalies of induced current associated with geomagnetic diurnal variation and earthquakes. The detachment slip events detected from the overlapping segment anomalies may be similar to the inter-plate slow earthquake or slow slip involving the upwelling and migration of thermal fluid under the continent, but this speculation needs to be demonstrated based on the research results of seismology.
The Bolokenu-Aqikekuduke Fault(Bo-A Fault)is a large-scale right-lateral strike-slip fault zone, which starts in Kazakhstan in the west, enters China along the NW direction, passes eastward through Alashankou, Lake Aibi and the southwestern margin of Turpan Basin, and terminates in the Jueluotage Mountain, with a total length of about 1 000km. At present, researches on the fault mainly focus on the area from Lake Alakol to Jinghe.Through satellite images, it can be found that the Bo-A Fault enters the southwestern margin of the Turpan Basin in the SE direction, and offset various landforms such as river terraces and alluvial fans, forming clear linear features on the surface, which indicates that there have been obvious activities since late Quaternary in this fault section. However, no detailed research has been carried out on the tectonic deformation characteristics of the Bo-A Fault in this area. The active characteristics of the faults in the southwestern margin of the Turpan Basin are studied, and the results are helpful to understand the role of the Bo-A Fault in the Cenozoic tectonic deformation of the Tianshan Mountains.The study area is located in the southwestern margin of the Turpan Basin, where three stages of alluvial-proluvial fans are developed. The first-stage alluvial-proluvial fan is called Fan3, which was formed earlier and its distribution is relatively limited, formed roughly in the early late Pleistocene; The second-stage alluvial-proluvial fan is called Fan2, which is the most widely distributed geomorphological surface in the study area. The geomorphic surface in this period was roughly formed from the late Pleistocene to the early Holocene. The third-stage alluvial-proluvial fan is called Fan1, which belongs to the Holocene accumulation, most of which are located at the outlet of gullies near the mountain passes, forming irregular fan-shaped inclined surfaces.To the west of Zulumutaigou, the fault offset the Fan3 alluvial-proluvial fan, forming dextral dislocation and fault scarp of the gully on the surface. The measurement shows that the amount of the dextral dislocation produced by the fault is between 22m and 40m. The height of the scarp is 3.9~4.2m. The section exposed by the fault shows that the Paleozoic bedrock thrust northward onto the Quaternary gravel layer, and the fault fracture width is about 1m, which reflects that the Bo-A Fault also has a certain thrust component. On the east bank of Zulu Mutaigou, the fault offset the Fan3 alluvial-proluvial fan, and the measurement results show that the offset of the gully is between 46.3m and 70.2m. To sum up, the movement mode of the Bo-A Fault in the study area is dominated by dextral strike-slip.On the Fan2 alluvial-proluvial fan at the northwest of Zulu Mutaigou, there are two secondary faults arranged in a right-step en-echelon pattern, forming high scarps with a height of 1.6~3.9m on the surface. Trench profiles reveal that both faults are SW-dipping thrust faults, thrusting from south to north, and they are preliminarily judged to be formed by the expansion of the Bo-A Fault into the basin.There are mainly three stages of alluvial-proluvial fans developed in the study area. Although no specific dating results have been obtained in this work, we believe that the age of the Quaternary landforms in the study area is the same as that in the Chaiwopu Basin, which is only separated by a mountain. Quaternary geomorphological ages are basically the same. Through geomorphological comparison, we believe that the age of Fan2 alluvial-proluvial fan is 12~15ka, and the age of Fan3 alluvial-proluvial fan is 74ka. It is estimated that the dextral slip rate of the Bo-A Fault is about 1mm/a since the formation of Fan3, and the vertical movement rate of the fault is about 0.13~0.32mm/a since the formation of Fan2.According to GPS observations and geological data, the NS-direction shortening rate in the East Tianshan area can reach 2~5mm/a. Through this study, it can be found that the Bo-A Fault also plays a role in regulating the near-NS-trending compressive stress in the East Tianshan area by accommodating the compression strain inside the Tianshan Mountains mainly through the NWW-directed right-lateral strike-slip motion. In addition, in the study area, the youngest fault scarp is located on the Fan2 alluvial-proluvial fan at the north of the main fault. It is preliminarily judged that the latest activity of the Bo-A Fault has a tendency to migrate from the mountain front to the basin.
The attenuation characteristics and site response are calculated respectively for each individual tectonic unit in Sichuan (Sichuan Basin,west Sichuan plateau and Panzhihua-Xichang area),using digital waveform data recorded by regional seismic networks and relevant seismic phase data collected from China Seismograph Network.The frequency dependent Q(f) is obtained by the iterative grid-search technique described by Atkinson and Mereu based on trilinear geometrical spreading model.The source spectra are determined by the model of Brune and the site responses of seismic stations are derived by Moya's method using genetic algorithms.Comparison to conventional ML estimates shows that the network local magnitude bias is quite significant at low and intermediate magnitudes.The bias at the jth station for the ith event is defined as ΔMij=Mij-Mi, where ΔMij is the station magnitude and Mi the network-average value.For comparison,we mapped the spatial distribution of biases by digital seismograms recorded from 10535 earthquakes of magnitude 2.5≤ML≤4.9 that occurred in Sichuan from January 1,2009 to June 30,2015.Based on the above data,the attenuation characteristics,site response and their effects on magnitude determination in Sichuan are analyzed.Our results demonstrate that the associated model for regional quality factor for frequencies can be expressed as Q1(f)=450.6f0.513 4 for Sichuan Basin,Q2(f)=136.6f0.581 3 for west Sichuan Plateau and Q3(f)=101.9f0.666 3 for Panzhihua-Xichang area.Site response results indicate that different stations show different amplifications.Maps of biases appear to be different,but with similar dominant spatial distribution.For stations in Sichuan Basin,their greater magnitudes are functions of low attenuation in structure and amplification effects of both seismic stations and basin effects.For stations in west Sichuan Plateau,the possible causes of these lower magnitudes are severe dependence upon source region due to extreme lateral variations in either structure or path effect attenuation.For stations in Panzhihua-Xichang area,broken medium caused by strong tectonic activity or large earthquakes and heat flow up-welling along active faults may be the main reasons of low magnitude values when earthquakes occur in western Sichuan and eastern Tibetan region.And the greater magnitudes for earthquakes along the Longmen Mountains appear to be well correlated with edge effect of sedimentary basin on strong ground motion.In our study,stations magnitude biases appear to be extremely correlated with tectonic structures and different regions when seismic rays passing through,magnitudes are affected significantly by lateral variations in attenuation characters rather than site responses.
The Wangjiagou Fault set,a set of Holocene active faults,is located at western suburbs of Urumqi City.The faults dislocated the gravel platform of the mid Pleistocene and the third level terrace of the Wangjiagou east bank,generating apparent fault scarps of opposite-slope direction on the surface with clear geomorphic traces.There are a series of deformation indications on landform,such as seismic fault,scarp and upheaval.In the field,thirty-nine groups of data were measured by using line tape along the fault.Among them,six were measured on the third level terrace of the Wangjiagou,and the others on the mid Pleistocene platform.Based on the data measured across the fault,we obtain that the height of the scarps is 0.4~1.6m and the width of the fault deformation is about 50m on the third level terrace.And on the mid Pleistocene platform,the height of scarps is 1.5~5.0m and the width of the deformation is about 90m.After comparing the profile of strong topographic deformation zone with the trench section,we primarily recognize that the ratio of hanging wall to foot wall deformation width is 2: 1approximately.The widths of strong surface deformation belt on the mid Pleistocene platform and the third level terrace on the two walls are 60m,30m and 33m,17m,respectively.