TS3.7 | Across the time scales, from earthquakes to earthquake cycle
EDI
Across the time scales, from earthquakes to earthquake cycle
Co-organized by G3/SM4
Convener: Y. Klinger | Co-conveners: Magali RizzaECSECS, Harsha Bhat, Alice-Agnes Gabriel
Orals
| Mon, 24 Apr, 14:00–17:25 (CEST)
 
Room K1
Posters on site
| Attendance Tue, 25 Apr, 16:15–18:00 (CEST)
 
Hall X2
Posters virtual
| Attendance Tue, 25 Apr, 16:15–18:00 (CEST)
 
vHall TS/EMRP
Orals |
Mon, 14:00
Tue, 16:15
Tue, 16:15
During the last decades, methods have significantly improved in geophysics, geodesy, and in paleoseismology-geomorphology. Hence, on one hand the number of earthquakes with well-documented rupture process and deformation pattern has increased significantly. On the other hand, the number of studies documenting long time series of past earthquakes, including quantification of past deformation has also increased. In parallel, the modeling community working on rupture dynamics, including earthquake cycle is also making significant progresses. Thus, this session is the opportunity to bring together these different contributions to foster further collaboration between the different groups focusing all on the same objective of integrating earthquake processes into the earthquake cycle framework. In this session we welcome contributions documenting earthquake ruptures and processes, both for recent or ancient events, from seismological, geodetic, or paleoseismological perspective. Contributions documenting deformation during pre-, post-, or interseismic periods, which are highly relevant to earthquake cycle understanding, are also very welcomed. Finally, we seek for any contribution looking at the earthquake cycle from the modeling perspective, especially including approaches mixing data and modeling.

Orals: Mon, 24 Apr | Room K1

Chairpersons: Y. Klinger, Magali Rizza, Harsha Bhat
14:00–14:05
14:05–14:15
|
EGU23-11854
|
ECS
|
On-site presentation
Sevil Cansu Yıldırım, Fatih Bulut, and Aslı Garagon

Using GPS measurements, historical earthquake records, and instrumental earthquake data, we investigated GPS slip rates along the rupture zone of the 1668 Great Anatolian Earthquake (M8.1). We found three complete and one incomplete earthquake cycles since 1254 compiling all available historical and paleo-earthquake records in the literature. These records verified that a ~750-kilometer section of the North Anatolian Fault Zone was ruptured in 1668.  To simultaneously estimate segment-based slip rates and locking depths, we combined all available GPS measurements and modeled them using an arctangent approach. Slip rates are used to estimate preliminary inter-seismic slip storages assuming fault segments are fully locked after a mainshock. Large residuals between preliminary slip estimates and co-seismic slips indicate that the fault segments do not store slip for some time after a major earthquake. The creeping and locked stages vary in time and space, as our investigation revealed. Our results show that the slip rates along the NAFZ systematically increase from east to west suggesting that the Aegean extensional regime is the main driving force for the westward movement of the Anatolian Plate. Additionally, the locking depths show an east-to-west decreasing pattern verifying east-to-west thinning of crustal thickness along the Anatolian Plate. The earthquakes over the past three complete cycles and the current incomplete cycle indicate that the failure of the NAFZ begins in the east and moves westward reflecting a decelerating pattern. The failure is typically completed within a time period of 239±3 years.

How to cite: Yıldırım, S. C., Bulut, F., and Garagon, A.: East to West Acceleration of the Slip Rates Along the North Anatolian Fault and Its Implıcations Regarding Plate Tectonics and Earthquake Cycle, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-11854, https://doi.org/10.5194/egusphere-egu23-11854, 2023.

14:15–14:25
|
EGU23-7909
|
On-site presentation
Murat utkucu, Hatice durmuş, Fatih uzunca, Süleyman nalbant, and Serap kIZILBUĞA

The M7.4 1999 İzmit earthquake apparently advanced the occurrence of possible future event or events along the segments of North Anatolian Fault Zone (NAFZ) beneath the Eastern Marmara Sea due to the positive stress load. This part of the NAFZ did not produce any large earthquake since the May 1766 earthquake, constituting a seismic gap close to the city of Istanbul. In the present study we constructed a Coulomb stress evolution model for the seismic gap that includes the effect of coseismic, time-dependent postseismic viscoelastic relaxation of the substrate beneath the elastic crust and secular stress loadings through the multiple earthquake cycles since 715 AD. The snapshots of stress changes before and after the large and destructive earthquakes of 740, 989, 1343, 1509, May 1766 and 1999 İzmit have been carefully examined. It has been estimated that the total stress changes before 989, 1343, 1509 and May 1766 earthquakes were in the range from 26 to131 bars. Present stress values at the eastern, middle and western sampling points on the faults within the gap are computed as 115, 131 and 85 bars respectively. Considering that the global mean of stress drops for continental strike-slip faults is about 35 bars, it is suggested that the earthquake hazard for the seismic gap critically high.

How to cite: utkucu, M., durmuş, H., uzunca, F., nalbant, S., and kIZILBUĞA, S.: The time-dependent stress changes within the seismic gap of the Eastern Marmara Sea (NW Türkiye) through multiple earthquake cycles since 715 AD., EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-7909, https://doi.org/10.5194/egusphere-egu23-7909, 2023.

14:25–14:35
|
EGU23-3915
|
ECS
|
Virtual presentation
Ziming Liu and Teng Wang

High-resolution geodetic measurements of the accumulated strains along active faults are important for faulting dynamics studies and seismic hazard evaluation. InSAR has been widely applied to measure the interseismic strain along active strike-slip faults. However, phase unwrapping errors, tropospheric delays, along with over-smooth effects in calculating the strain from velocity limit its capability of mapping the highly localized strain along faults. Phase-gradient stacking that sums up the wrapped phase differences of adjacent pixels has been successfully applied to reveal localized deformation across coseismic fractures and slow-moving landslides, yet lacks application to reveal interseismic strain along faults. Here, we conduct phase-gradient stacking on Sentinel-1 SAR interferograms, for the first time, to map the interseismic strain along the North Anatolian Fault with unprecedented resolution. We reveal several segments with extremely high strain rates attributed to shallow creep of the fault. By comparing with historical earthquake ruptures, we find that the creeps are either related to afterslip of recent earthquakes, or related to slip deficits of earthquakes occurred decades ago, challenging the opinion that the NAF has a uniform surface strain rate, particularly along its eastern portion. Our results show that the phase-gradient stacking can not only reduce the computation burden from phase unwrapping and tropospheric correction, but also achieve a much higher spatial resolution strain map than the traditional InSAR method. The proposed method can be applied to other large strikes-slip faults for distinguishing segments with surface creep and strong coupling and therefore better quantify the shallow strain budget and its associated hazards.

How to cite: Liu, Z. and Wang, T.: High-resolution Interseismic Strain Mapping from InSAR Phase-Gradient Stacking: Application to the North Anatolian Fault with Implications to the Non-uniform Strain Distribution Related to Historical Earthquakes, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-3915, https://doi.org/10.5194/egusphere-egu23-3915, 2023.

14:35–14:45
|
EGU23-11401
|
ECS
|
On-site presentation
Kaan Alper Uçan and Fatih Bulut

This study aims to forecast the magnitude of future strong (6.0≤M<7.0) and major (7.0≤M<8.0) earthquakes along the East Anatolian Fault Zone (EAFZ), a major fault zone of Turkey and an active plate boundary that lies between Arabian and Anatolian plates. We first investigated the segmentation of the EAFZ in this context after compiling the earlier research on its structural setting and historical earthquakes. In order to determine the distribution of slip deficit rates, we analyzed GPS slip rates to obtain back-slips. The current slip budgets on each fault segment are calculated using the resulting slip deficit estimates. To elaborate on whether b-values might be used to distinguish between locked and creeping fault segments, we also examined the distribution of b-values along the fault. As a result, we found a reverse correlation between slip deficit rates and b-values. According to our findings, the EAFZ has currently a slip deficit of 1.51 m. While there is a segment such as Hacılar with no slip deficit, there is enough slip deficit accumulation to generate three strong and three major earthquakes on the other fault segments. Presently, these fault segments have the potential to re-generate previous earthquakes, within the magnitude range of 6.8-7.4. The latest strong earthquake on January 24, 2020, the Elazığ earthquake (M 6.8) verified our magnitude forecasts for the Sivrice-Pütürge segment.

How to cite: Uçan, K. A. and Bulut, F.: Forecasting Earthquake Magnitudes along the East Anatolian Fault Zone using Fault Zone Segmentation, Historical Earthquakes, and GPS Slip Rates, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-11401, https://doi.org/10.5194/egusphere-egu23-11401, 2023.

14:45–14:55
|
EGU23-4233
|
ECS
|
On-site presentation
Wenqian Yao, Jing Liu-Zeng, and Zijun Wang

The propagation of the 2021 Mw7.4 Madoi earthquake rupture from the central Jiangcuo fault (JCF) onto the eastern portion exhibits the most complex geometry with a series of conjugate faults, bends, and stepovers. At the east ~50 km of the 2021 epicenter, the surface rupture along the Jiangcuo eastern branch (JCEB) deviating ~12° anticlockwise from the general strike provides a valuable chance for understanding the particularly complex surface ruptures propagation and the branching behavior of the poorly known JCEB. Using sub-metric orthophotos collected by UAV with a ground resolution of 6 cm, complemented by multiple field investigations, we implemented the surface rupture mapping and coseismic slip distribution of the JCEB in detail associated with this earthquake sequence. Our mapping illuminated the sporadic breaks of the tectonic region in the dune area immediately near the branching point and eastward propagated linear rupture trace. The measurements of the high-resolution coseismic slip along the JCEB show that the slip distribution reveals an approximate dogtail shape to the eastern termination with the maximum left-lateral strike-slip offset of 2.9 m. These data might support the perspective that the rupture propagated with a supershear velocity toward the east. Combined with the accrued displacements along the JCEB, these results indicate that the poorly known divergent branch could accumulate pre-2021 surface breaks as an immature fault and bifurcated in the Madoi quake due to the matched regional stress field. We found linear surface breaks along the NW-strike geologic faults indicating triggered coseismic slip on conjugate faults. In the meantime, the intersections with conjugate faults mark discontinuities in rupture geometry and surface slip on the main fault, suggesting strong fault interaction in the eastern tip zone of the Madoi rupture.

How to cite: Yao, W., Liu-Zeng, J., and Wang, Z.: Rupture Branching and Propagation at the Eastern End of the 2021 Mw 7.4 Madoi Earthquake, North Tibet Plateau, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-4233, https://doi.org/10.5194/egusphere-egu23-4233, 2023.

14:55–15:05
|
EGU23-7209
|
ECS
|
On-site presentation
Liqing Jiao, Teng Wang, Guangcai Feng, Paul Tapponnier, Andrean V. H. Simanjuntak, and Chung-Han Chan

Rupture complexity results in difficulty in quantifying seismic hazards, such as the probability of an earthquake on multiple fault segments and spatial distribution of fault displacement on the surface. Here we tried to propose a dynamic model to fit the rupture behavior of the 2018 Mw7.5 Palu earthquake, which splays along several sub-fault plans on the surface. The Palu event was initiated on an unknown fault and propagated on a curved plane on the Palu-Koro and Matano faults. According to the Interferometric Synthetic Aperture Radar (InSAR) data, both principal (on-fault) and distributed (off-fault) faulting were identified and spatial displacement on the surface could be evaluated. To model the complex geometry of the coseismic rupture plane and corresponding deformation, we proposed a dynamic model through the Discrete Element Method (DEM). Our model demonstrated rupture along a planar fault at depth and several splay faulting with various deformations on the surface, corresponding to the observations. The simulations represented temporal rupture behavior that covers several earthquake cycles and the probability of superficial fault displacement, which shed light on subsequent seismic hazard assessment and probabilistic fault displacement hazard analysis, respectively.

How to cite: Jiao, L., Wang, T., Feng, G., Tapponnier, P., Simanjuntak, A. V. H., and Chan, C.-H.: Impact of rupture complexity on seismic hazard: Case of the 2018 Mw7.5 Palu earthquake, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-7209, https://doi.org/10.5194/egusphere-egu23-7209, 2023.

15:05–15:15
|
EGU23-6567
|
ECS
|
On-site presentation
Rémi Matrau, Yann Klinger, Jonathan Harrington, Thorvaldur Thórdarson, Ármann Höskuldsson, Esther Gudmundsdöttir, Laura Parisi, Margherita Fittipaldi, Ulas Avsar, and Sigurjón Jónsson

Studies of Oceanic Transform Faults (OTFs) usually rely on geophysical data because of the OTF inaccessibility on the seafloor. The Húsavík-Flatey Fault (HFF) in northern Iceland is an OTF connecting the onshore rift in Iceland to an offshore rift segment of the Mid-Atlantic Ridge, accommodating 30% to 50% of the relative plate motion at this latitude between North America and Eurasia. The HFF is unique because its easternmost 25 km-long segment is exposed on land, allowing to study the long-term deformation of the fault. Two historical earthquakes of estimated magnitudes M6.5 - M7 have been reported on the eastern HFF in the last 270 years. However, almost no information exists from prior to the 18th century.

To study the Holocene deformation of the HFF and to build a catalogue of past earthquakes, we excavated 11 paleoseismology trenches at two locations, six on an alluvial fan and five in a pull-apart basin. We also excavated and tracked buried river channels to estimate long-term slip rates and to assess the coseismic displacement of single events. We used radiocarbon dating of birch wood samples together with major element compositions of volcanic ashes (tephras) to constrain the timing of events on the fault.

Trenches at both locations show clear evidence of deformation and surface rupturing events. From offset measurements of glacial morphologies and buried river channels, we calculate a Holocene slip rate of 4 - 6 mm/yr, slightly lower than the estimated present-day geodetic slip rate, suggesting that some of the deformation may be distributed. Based on upward terminations of cracks and faults, we identified eight events in the last ~6000 years, yielding fewer major earthquakes than expected from the 270-year record. We thus suggest that large earthquakes of magnitude ~M7 on the HFF, producing significant surface ruptures, are rare, with a return time of 500 to 600 years. We also propose that the short recurrence times often observed on OTFs may therefore not be representative of the full seismic cycle.

How to cite: Matrau, R., Klinger, Y., Harrington, J., Thórdarson, T., Höskuldsson, Á., Gudmundsdöttir, E., Parisi, L., Fittipaldi, M., Avsar, U., and Jónsson, S.: Holocene deformation on a transform fault: Insights from paleoseismology on the Húsavík-Flatey Fault in North Iceland, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-6567, https://doi.org/10.5194/egusphere-egu23-6567, 2023.

15:15–15:25
|
EGU23-5124
|
ECS
|
On-site presentation
Magali Riesner, Lucilla Benedetti, Stéphane Baize, Stefano Pucci, Matthieu Ferry, Stéphanie Gautier, Régis Braucher, Jules Fleury, Hervé Jomard, Stéphane Mazzotti, and Fabio Villani

Long-term fault escarpments are built by the accumulation of individual earthquakes producing incremental surface displacements on the fault releasing crustal tectonic loading. Cumulative escarpment studies have revealed a spatial slip variability along active faults as well as a temporal variability with the alternation of phases of intense seismic activity over a short period of time followed by long periods of quiescence. Understanding this spatial and temporal slip variability on individual faults and over a complex fault system provide a better knowledge of co-seismic rupture extents, essential for estimating past earthquakes magnitude and for seismic hazard assessment.

Up to now, most studies have focused on a timeframe over few seismic cycles, making it difficult to apprehend the rupture barriers persistence and cumulative slip distribution.  Here, we aim at quantifying the slip variability over several timescales ranging from a few months to a few million years on the same fault.

Our study focusses on the ~50 km-long Liri fault, SW of the Fucino basin. The fault is located at the contact between Cretaceous limestone and patches of Quaternary deposits locally convering Mio-Pliocene flysch sediments. Detailed mapping of the fault trace on high-resolution Digital Elevation Model (DEM) from UAV-acquired images, Pleiades images and Lidar together with field observations revealed changes in the morphological expression of the fault north and south of an important wind gap located at Capistrello. To the north, the faut trace is ~16 km-long located on the eastern side of ~2km-wide limestone ridge, reaching ~1300m asl elevation. Two bends in the fault trace, made of ~5km long segments, can be observed with the fault strike varying between N115° and N140°. In this northern section, the fault scarp appears subtle and we did not observe Quaternary deposits on the hanging wall. In the 30 km-long section, south of Capistrello, the cumulative scarp composed of numerous splays is evidenced by a sharp trace, offsetting several morphological surfaces and associated Quaternary sediment packages. Three major bends are observed in this section of the fault, separating 10 to 30 km-long segments striking between N110° and N160°. An alluvial surface offset by ~14 m of cumulative displacement was dated at ~35kyr using 36Cl cosmogenic exposure dating suggesting a minimum slip rate of 0.4 mm/yr.  Other morphological markers that have accumulated displacement between ~10 and 70 m-high have also been sampled for 36Cl cosmogenic exposure dating. Moreover, we excavated two small trenches at the base of the fault scarp within the Quaternary deposits affected by the fault revealing 3 rupture-surfacing earthquakes over the last 2500 yr, the last one occurring after 1226 CE. 

We will present those results and will discuss how the displacement varies along the fault both in time and space.

How to cite: Riesner, M., Benedetti, L., Baize, S., Pucci, S., Ferry, M., Gautier, S., Braucher, R., Fleury, J., Jomard, H., Mazzotti, S., and Villani, F.: Quantifying the slip over various time scales on active normal faults in the Apennines (Italy):  the Liri fault from paleoearthquakes to long-term slip rate, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-5124, https://doi.org/10.5194/egusphere-egu23-5124, 2023.

15:25–15:35
|
EGU23-3813
|
ECS
|
On-site presentation
Octavi Gómez-Novell, Bruno Pace, and Francesco Visini

A major challenge in seismic hazard research is to quantify the frequency of large earthquakes along active faults, more so when the observational time windows of seismic catalogs are much shorter than the average fault recurrence intervals. In this respect, paleoseismology continues to prove to be an excellent tool to extend the seismic catalogs of faults into prehistorical times. The combination of the ever more advanced trenching surveys and accurate numerical dating techniques allows constraining the timing of paleoearthquakes and, for some datasets, approximating their recurrence models. Despite this, paleoseismic data carries along large uncertainties frequently related to dating technique limitations, poor stratigraphic preservation, and along-strike slip variability that hinder the identification of a complete paleoearthquake record. Subsequently, these issues, among others, challenge the constraint of reliable earthquake chronologies along faults and of the parameters defining their earthquake cycle.

We present an automatic workflow capable to compute and constrain earthquake chronologies along a fault based on the correlation of its available paleoseismic records, including multi-site and poorly constrained datasets. Our inherent premise is that the correlation of paleoseismic data from multiple along-fault locations can help to improve the time constraints and completeness of its paleoearthquake record. Given that paleoseismic records are, by definition, underpopulated, event correlation is not restricted to single occurrences. Instead, an event in a site might be simultaneously correlated with more than one in another if time compatible. Furthermore, to avoid subjectivity biases in event timing estimates and correlation, we exclusively rely on the trench numerical dates limiting each event horizon as the inputs. All earthquake chronologies are modelled probabilistically with a four-step algorithm as we detail. First, all earthquake times in each site are computed as probability density functions (PDFs) using the input numerical dates. The event PDFs from all sites are then integrated to derive a mean curve representing the overall event probabilities for the studied fault in the time span investigated. The probability peaks in this curve, which are assumed as indicative of the event timing at the fault scale, are automatically detected based on peak prominence analysis. A final PDF is then computed for each peak by multiplying all site event PDFs intersecting the peak position. The set of product PDFs constitutes the earthquake chronology of the fault, provided to the user in simple output files that can be externally used to calculate fault parameters for the seismic hazard assessment, and to visualize the modelling.

Preliminary tests on several paleoseismic datasets of the Central Apennines (Italy), the Eastern Betics (Spain), the Dead Sea Fault and the Wasatch Fault (US), have provided good outcomes. The approach significantly reduces the uncertainties in event timing of paleoearthquakes and provides an objective and reliable interpretation of the datasets, especially when these are complex or have wide uncertainties. By extension, the workflow has the potential to reduce the uncertainties in earthquake recurrence estimates and can give insight on the recurrence models that better describe the earthquake cycle in the studied faults.

How to cite: Gómez-Novell, O., Pace, B., and Visini, F.: Automated workflow to compute earthquake chronologies on faults from paleoseismic datasets, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-3813, https://doi.org/10.5194/egusphere-egu23-3813, 2023.

15:35–15:45
|
EGU23-16199
|
ECS
|
On-site presentation
Lucia Andreuttiova, James Hollingsworth, Pieter Vermeesch, and Tom Mitchell

Earthquakes on normal faults in the continental setting are relatively uncommon. The scarcity of surface-rupturing events underpins an absence of surface displacement measurements. It is a common practice to use surface offset as a proxy to understand the fault structure at depth. Hence, the lack of comprehensive surface data impedes the subsurface reconstruction of seismogenic normal faults and prohibits the thorough assessment of earthquake hazards. To supplement the available surface displacement measurements and to make statistically significant inferences, we apply optical image correlation (OIC) methods to historical images from three large continental normal earthquakes in the western United States (1954 Dixie Valley (Mw 6.8) - Fairview Peak (Mw 7.1) earthquake sequence, the 1959 Mw 7.2 Hebgen Lake earthquake and the 1983 Mw 6.9 Borah Peak earthquake). The results of this study are displacement maps with three components of deformation from which we extract high-resolution 3-d measurements everywhere along the surface rupture. 

 

The high-resolution 3-d data are used to quantify the magnitude and direction of the earthquake-related offset, the percentage of off-fault damage as well as the width of the fault zone. These parameters represent the fault maturity, geometric complexity and subsurface structure of the fault. Our observations confirm behaviours previously observed along strike-slip faults (e.g. magnitude of off-fault deformation is proportional to the rupture complexity). In addition, a comparative assessment of the results from the three study areas demonstrates that features such as excess slip detected close to the fault scarp are not unique and can be found along multiple dip-slip faults. Consequently, this study documents the variation of the quantifiable parameters along the normal faults. It suggests that while some parameters are a universal reflection of the fault characteristics, others vary according to the geology or topography in the area and should not be accepted without further investigation.

How to cite: Andreuttiova, L., Hollingsworth, J., Vermeesch, P., and Mitchell, T.: Assessing distribution and pattern of the earthquake-related deformation caused by large continental normal earthquakes using optical image correlation, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-16199, https://doi.org/10.5194/egusphere-egu23-16199, 2023.

Coffee break
Chairpersons: Y. Klinger, Magali Rizza, Alice-Agnes Gabriel
16:15–16:25
|
EGU23-8011
|
ECS
|
On-site presentation
Xiaohang Wang, Mahdi Motagh, and Caijun Xu

The Tianshan range, one of the most active mountain building belts in central Asia, has complex geological structures and frequent strong earthquakes since the Cenozoic. Due to lack of sufficient high-resolution geodetic survey measurments, little is known about detailed fault slip rates and seismic hazards related to main active faults in Tianshan. However, in recent years, the improvements in space-based geodetic technologies (Global Navigation Satellite System (GNSS) and interferometric synthetic aperture radar (InSAR)) with growing coverage and accuracy provide us an opportunity to image more subtle features in this area. In this study, we assesses inter-seismic deformation for the period 2014-2022 over the southwestern Tianshan based on ascending and descending Sentinel-1 SAR data.  Combined with GNSS data, we then constructed the 3D crustal deformation with high precision and high spatial resolution to study the active structures in southwestern Tianshan. The results indicate that: (1) The Tianshan orogenic belt (TSOB) has intense crustal deformation and the shortening rate is approximately 20 mm/yr. The Keping fold-thrust belt (KFB) is the most intensely deformed areas in TSOB, it’s convergence rate accounts for 1/3 of the entire southwestern Tianshan. (2) The South Tianshan fault (STF) and the Piqiang fault (PQF) have obvious left-lateral strike-slip components and the South Tianshan fault also has thrust characteristic. (3) The folds in both western and eastern KFB play an important role in accommodating regional strain, the shortening rate in KFB is accommodated by the thrust-anticlinal zone at the KFB front.

How to cite: Wang, X., Motagh, M., and Xu, C.: Crustal deformation of southwestern Tianshan orogenic belt based on InSAR and GPS observations, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-8011, https://doi.org/10.5194/egusphere-egu23-8011, 2023.

16:25–16:35
|
EGU23-11616
|
On-site presentation
|
Henriette Sudhaus, John Begg, Vasiliki Mouslopoulou, Julia Knüppel, and Tilman May

As well as slip on a primary fault plane, earthquakes can produce slip on neighbouring faults which are not directly linked to the main source. This slip is called syn-seismic. With modern space-borne observation techniques, we observe syn-seismic slip down to a few centimeters on active faults nearby the source. An excellent example is the mapped slip on secondary faults during the 2019 Ridgecrest earthquake sequence in California. The overall spatial pattern of syn-seismic slip with respect to the main fault suggest that these faults respond to local stress changes caused by the main shock.

Data that enable the detection of surface fault slip on such small scale are provided by optical and radar satellites which allow a very high precision with high spatial resolution. In particular, short revisit times of these satellite observations lead to high coherence between images matched in pixel-offset and radar interferometric techniques.

We present further examples of syn-seismic fault slip during ~M6 earthquakes from different regions, such as those recorded in Greece in 2021 (Tyrnavos and Arkalochori) and 2020 in Tibet (W Xizang and near Xegar). We use Sentinel-1 interferometric wide-swath SAR acquisitions, which we process on the highest spatial resolution and apply weak filtering only. Our examples have in common that their syn-seismic fault activation reveals slip of a few centimeters only, persistently along a section of the fault’s length. The slip directions commonly appear to follow the coseismic surface displacement gradients which, in some cases, results in reverse slip on long-term normal faults. The activated faults were either faults previously mapped or concealed faults which were identified due to InSAR.

It is difficult to estimate the depth of syn-seismic fault slip and therefore how much strain has been released due to localized stress changes. We are also uncertain of the extent to which this small slip release contributes to the long-term displacement and displacement rate on faults and whether its contribution should be included in dislocation fault slip models. Our compilation suggests that syn-seismic slip is rather common, despite the rarity of previous observations, and is now detectable only because of improved resolution provided by InSAR data.

How to cite: Sudhaus, H., Begg, J., Mouslopoulou, V., Knüppel, J., and May, T.: InSAR observations of syn-seismic slip on faults due to M~6 earthquakes, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-11616, https://doi.org/10.5194/egusphere-egu23-11616, 2023.

16:35–16:45
|
EGU23-12158
|
ECS
|
On-site presentation
Hugo Boulze, Luce Fleitout, Emilie Klein, Christophe Vigny, and Jean-Didier Garaud

Thanks to space geodesy we know with a millimetric precision how the lithosphere deforms at each stage of the seismic cycle. In particular, during the post-seismic phase, it can deform over thousands of kilometers and for decades. These deformations are partly due to viscoelastic relaxation of the asthenosphere.

In a previous work, we have shown that at the temporal and spatial scale of the seismic cycle, the viscoelastic relaxation can be modeled by a linear creep law [Boulze et al. 2022]. Therefore, because of the linearity of the creep law, the superposition principle applies and the present day deformation is simply the sum of the post-seismic deformations induced by past earthquakes. Based on this result, the objective of our work is to determine what slip history is needed on the Chilean subduction interface to reproduce the current deformation of South America, which is well measured by GNSS.

To investigate this challenging problem, we first develop a 3D spherical finite-element model of the Chilean subduction zone. This model covers the entire South American continent and incorporates a slab with a geometry described by Slab2.0 model [Hayes et al. 2018]. Then, we compare different ways to model the seismic cycle using the backslip theory [Savage 1983]. Finally, by comparing GPS time-series with our seismic cycle model prediction, we discuss many ingredients of the model: e.g. the viscosity of the asthenosphere (Maxwell, Burgers), the impact of a flat slab and low viscosity zones, the magnitude and extent of historical earthquakes.

How to cite: Boulze, H., Fleitout, L., Klein, E., Vigny, C., and Garaud, J.-D.: Modeling surface deformations during the seismic cycle along the Chilean subduction zone, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-12158, https://doi.org/10.5194/egusphere-egu23-12158, 2023.

16:45–16:55
|
EGU23-16
|
On-site presentation
Interseismic locking on a megathrust controls long term strain partitioning in a subduction zone
(withdrawn)
Simon Lamb
16:55–17:05
|
EGU23-15420
|
ECS
|
Virtual presentation
Heng Luo and Teng Wang

Foreshocks are commonly observed before the happening of earthquakes in seismic catalogs. They provide critical precursors to reveal the process for the nucleation and rupture of earthquakes. Two mechanisms, pre-slip and cascade triggering, are thought to be the main physical process to explain the foreshock sequences and the mainshock. However, different from the regular micro-magnitude foreshock sequences (e.g. M1.0-3.0), some moderate-size (e.g. ~M6) foreshocks are also found before the mainshock (e.g. the M6.4 foreshock before the 2017 M7.1 Ridgecrest earthquake). How these moderate-size foreshocks affect the happen of mainshocks as well as their possible triggering mechanisms are still ambiguous and less studied.

In this study, fortunately, we obtain geodetic observations of moderate-size (M5.8 and M6.5) foreshocks for the 2020 M6.0 Turkey and 2022 M6.9 Taiwan earthquakes using Sentinel-1 Synthetic Aperture Radar (SAR) images. It is very rare for the geodetic observations of such foreshocks as they are very temporally close to the mainshocks within one day (i.e. ~10 hours and ~17 hours). We then invert for the fault geometries and slip distributions for these two earthquakes together with their moderate foreshocks constrained by these geodetic observations. Coulomb stresses on the fault planes of mainshocks produced by the moderate-size foreshocks are also calculated as well as the static stress drops of the mainshocks. Our study provides a unique opportunity to explore the possible triggering mechanism between moderate-size foreshocks and mainshocks as well as the conditions for the happening of earthquakes.

How to cite: Luo, H. and Wang, T.: Geodetic modeling of moderate-size foreshocks with mainshocks and the implication to earthquake trigger mechanisms, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-15420, https://doi.org/10.5194/egusphere-egu23-15420, 2023.

17:05–17:15
|
EGU23-2262
|
On-site presentation
Axelle Amon, Ambroise Mathey, David Marsan, Jerome Weiss, and Jerome Crassous

We study an experimental model of a fault consisting in a stationary shear band in a compressed granular sample. To obtain those bands, we perform a biaxial compression of a granular sample constituted of glass beads during which we observe the spontaneous formation of shear planes along the Mohr-Coulomb directions in the sample. We study the post-failure regime during which all the deformation occurs along the stationary shear bands. Using an interferometric method of measurement of micro-deformations based on multiple scattering, we obtain full-field measurements of the local incremental deformation in the sample. The deformation measured are typically of $10^{-5}$ with a resolution of about 300 microns (3 bead diameters). Our technics gives access to the strain fluctuations inside the shear band and we show that the macroscopic mean deformation in the bands is the result of the accumulation of local, intermittent, shear events. The size distribution of those shear events follows the Gutenberg-Richter law. We observe clustering of those events following Omori's law and we apply a declustering method to reveal the causal structure underlying our sequences of events (Houdoux et al. 2021). 

In my talk, I will focus on recent experimental results regarding the dependence of the series statistics on the driving velocity. We have studied sequences of aftershocks for different compression velocities and we have shown that surprinsingly the aftershock sequences we observe are deformation-dependent and not time-dependent. We discuss such a deformation memory effect in the framework of an Olami-Feder-Christensen model.

Houdoux et al. Commun Earth Environ 2, 90 (2021)

How to cite: Amon, A., Mathey, A., Marsan, D., Weiss, J., and Crassous, J.: Deformation-dependent aftershocks in laboratory earthquakes sequences, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-2262, https://doi.org/10.5194/egusphere-egu23-2262, 2023.

17:15–17:25
|
EGU23-9213
|
On-site presentation
Thorsten Becker, Simone Puel, Umberto Villa, Omar Ghattas, and Dunyu Liu

Analysis of geodetic and seismological data helps constrain earthquake dynamics and the physics of lithospheric deformation. Here, we discuss a new modeling approach based on an open-source finite-element framework to invert surface deformation data for constitutive laws and their parameters, such as the Poisson’s ratio or shear modulus in the crust and mantle wedge.

These inversions can be realized by using adjoint-based optimization methods which efficiently reduce the misfit between the calculated and observed displacements. To quantify the associated model uncertainties, we extend the inverse approach to a Bayesian inference problem. Since the data are usually informative only in a few directions in parameter space, we use a low-rank Laplace approximation of the posterior distribution to make the inverse problem computationally tractable. The mean and the posterior covariance are approximated by the solution of the inverse problem (MAP point) and the inverse of the Hessian of the negative log posterior evaluated at the MAP point, respectively. We show how smoothly varying parameter fields can be reconstructed satisfactorily from noisy data.

To improve the spatial resolution of the inverse solution we solve a Bayesian optimal experimental design problem to find the best station configuration by maximizing the expected information gain, defined as the Kullback-Leibler divergence between posterior and prior distributions. We show how and why the optimal network improves the material property inference more than evenly spaced stations. Based on our previous work on inverting for fault slip without Green’s function computations, we combine the two inversion schemes to jointly infer both model parameters, the coseismic slip, and material properties distribution. Lastly, we test this numerical forward/inverse framework with an application, the 2011 Tohoku-oki M9 earthquake. Both continuous land-based and six offshore acoustic GNSS stations located around the earthquake epicenter are inverted to jointly estimate the shear modulus and the fault slip during the megathrust event.

Our results demonstrates the potential of our computational framework and the general approach for inferring constitutive laws to evaluate sensitivity to parameters, and define strategies to improve our understanding of relevant parameters for earthquake dynamics. 

 

How to cite: Becker, T., Puel, S., Villa, U., Ghattas, O., and Liu, D.: An Adjoint-based Method for Inverting for Heterogeneous Material Properties and Fault Slip From Earthquake Surface Deformation Data, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-9213, https://doi.org/10.5194/egusphere-egu23-9213, 2023.

Posters on site: Tue, 25 Apr, 16:15–18:00 | Hall X2

Chairpersons: Y. Klinger, Magali Rizza
X2.261
|
EGU23-16845
|
ECS
Solene L Antoine, Zhen Liu, Yann Klinger, Arthur Delorme, and Jing Liu-Zeng

The 2021 Mw7.4 Maduo earthquake generated a ~160 km-long fault rupture within the Eastern Tibetan plateau, at about 100-150 km to the south-west of the Eastern Kunlun fault. Fault slip measured on the field represents only 20% of the displacements from satellite Interferometric Synthetic Aperture Radar (InSAR) measurements, highlighting the primarily diffuse nature of the surface deformation for this earthquake. Most surface deformation associated with this event corresponds to diffuse shear, occurring over widths of a few hundreds of meters to a few kilometers, and sometimes associated with shearing and tensional cracks mapped in the field. In this study, we use sub-pixel correlation of Pleiades (0.5 m) and SPOT-6/7 (1.6 m) optical images to characterize the near-fault displacement patterns associated with the 2021 Maduo event. We also use other optical data to assess the impact of sensor resolution on the measurements. Our results cover three kilometers on both sides of the rupture area with a resolution of 0.5 m. These results show that, despite the large rupture gaps observed in the field, the shear deformation zone at the surface is continuous along the entire length of the 2021 rupture. Even though, we observe variations in the surface deformation patterns, with regions that present more localized deformation whereas others are primarily characterized by diffuse shear. Using the high-resolution displacement maps, we characterize the transitions from the localized to the diffuse shear along the rupture strike, and investigate the relations with the bulk rock properties, and coseismic slip distribution. We also determine the limit at which deformation starts to localize on fractures that are large enough to be visible in the field.

How to cite: Antoine, S. L., Liu, Z., Klinger, Y., Delorme, A., and Liu-Zeng, J.: Characterizing the transition from diffuse to localized deformation using optical image correlation: the 2021 Mw7.4 Maduo, Tibet, earthquake, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-16845, https://doi.org/10.5194/egusphere-egu23-16845, 2023.

X2.262
|
EGU23-552
|
ECS
Nicolás Pinzon Matapi, Yann Klinger, Xiwei Xu, Jing Liu, and Paul Tapponnier

Identifying earthquake recurrence times and slip distributions over the span of many seismic cycles is key to understand fault-rupture processes and to better assess the seismic hazard. In this study, we used three paleoseismological excavations along the Aksai segment of the Altyn Tagh Fault (ATF) to document preserved evidence of past earthquakes in the sedimentological record such as vertical offset, fault cracks, and folding. We integrated these findings with previous studies on the Annanba and Xorxoli segments in order to build a larger-scale rupture history of the ATF. We reported nine large paleo-earthquakes and three of these with ground rupture expression along the whole three segments (∼ 400 km). Based on a Bayesian approach we present 95-percentile range ages of 6149 – 5285 BC, 5296 – 4563 BC, 3026 – 2677, 2469 - 2254 BC, 2069 - 1964 BC, 1184 – 709 BC, 270 – 635 AD, 875 – 1325 AD and 1491 - 1741. Furthermore, we used high-resolution satellite imagery to measure horizontal offsets recorded in the morphology, which are associated with potential co-seismic deformation. We find that the mean recurrence time is 1171±425yr with a COV of ∼0.31 suggesting a quasi-periodic behavior with a characteristic slip motion based on the similar distribution of fault offsets. The last event seems to be strongly expressed in Xorxoli segment and also found along the Aksai segment, although we could not identify it along the Annanba bend. Whereas, the penultimate event and the two before this appear to well correlate across the Aksai, Annanba and Xorxoli segments. Thus, being strong candidates for the three largest and successive earthquakes along the ATF (roughly rupture longitude ≥ 350 km). Variations in the COVs along the eastern Altyn Tagh Fault accounts for the important control of local structural complexity and/or slip rate variations on the rupture behavior of major fault systems.

How to cite: Pinzon Matapi, N., Klinger, Y., Xu, X., Liu, J., and Tapponnier, P.: Long-term Earthquake Cycle along the eastern Altyn Tagh Fault, China, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-552, https://doi.org/10.5194/egusphere-egu23-552, 2023.

X2.263
|
EGU23-5953
|
ECS
Longfei Han, Jing Liu-Zeng, Guiming Hu, Yann Klinger, Wenxin Wang, Heng Wang, Jing Xu, Bo Zhang, Yunpeng Gao, Zijun Wang, Xianyang Zeng, and Xiaoli Liu

Paleoseismic records are essential for constraining the earthquake recurrence behavior of active faults and evaluating the rupture history. However, paleoseismic studies on the central Altyn Tagh fault (ATF) are still scarce, and previous studies indicate that this fault section with simple geometry is not periodic. In addition, paleoseismic data from two sites along central ATF reveal different amounts of paleoearthquakes and present discordant in time. Therefore, we conducted paleoseismic studies and documented six reliable paleoseismic events at the LaPeiQuan site along the straight section of the central ATF. The results indicate that the most recent event is a small earthquake with a tiny vertical offset. The data A.D. (1752–1880) yr (event A) is significantly later than the last event along the Xorkoli section. The penultimate event at the LaPeiQuan site is a large earthquake for the ages of this event B is A.D. (667–764) yr (event B), which is consistent with the Xorkoli site and Aksay double bend site, producing at least 140 km rupture. In addition, the large vertical offset measurement from the deformed sediment of event B also supports its large one. The ages of Event D are discordant with the adjacent paleoseismic sites. The ages of Event C, Event E and Event F are still in process. The reason earthquake histories are inconsistent may be that small-scale geometrical complexities can prevent earthquake rupture propagation.

How to cite: Han, L., Liu-Zeng, J., Hu, G., Klinger, Y., Wang, W., Wang, H., Xu, J., Zhang, B., Gao, Y., Wang, Z., Zeng, X., and Liu, X.: Paleoseismic Investigation along the straight section of the central Altyn Tagh fault and its constrain on the rupture history, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-5953, https://doi.org/10.5194/egusphere-egu23-5953, 2023.

X2.264
|
EGU23-14568
Newdeskarl Saint Fleur, Yann Klinger, Joseph Emmanuel Dessable, Germain Saint-Preux, Nathalie Feuillet, Dominique Boisson, Eric Calais, and Jean-Bernard de Chabalier

The 14 August 2021 earthquake occurred along the southern peninsula of Haiti only 11 years after the 12 January 2010 devastating earthquake. According to seismological and geodetic data, the events are both complex involving more than one fault. The 2021 rupture mainly portrayed reverse motion to the east near L’Asile town and left-lateral strike-slip motion to the west near Camp-Perrin town and Macaya mountain. A few days after the 2021 event, we conducted the first post-seismic field reconnaissance along the left-lateral Enriquillo-Plantain Garden Fault (EPGF) zone from L’Asile to Macaya mountain. We found numerous fresh cracks and landslides along that fault zone. The 111 cracks are mainly E-W-striking, some are oriented WNW-ESE, consistent with fault orientation in the area. In addition, the biggest cracks are mostly located to the west of the rupture zone, some of them may be potential fault surface rupture as revealed by seismological data. Furthermore, our observations along the northern coast of the southern peninsula revealed no significant coseismic coastal uplift as also suggested by InSAR data. Besides that field reconnaissance, we revisited the fault map around the epicentral area using high-resolution LiDAR data, Pléiades imagery and aerial photographs. We identified several left-lateral offsets of tens of meters corresponding to successive slips along the EPGF from L’Asile to Macaya mountain. In addition to the strike-slip deformation, we identified numerous geomorphic features related to long-term tectonic uplift to the north of the EPGF surface trace near the eastern part of the 2021 rupture. Those features are strikingly rare to the south. Such a pattern may indicate that the EPGF is north-dipping in the area. The 14 August 2021 rupture offers a new opportunity to constrain the kinematics and geometry of the EPGF system in southern Haiti.

How to cite: Saint Fleur, N., Klinger, Y., Dessable, J. E., Saint-Preux, G., Feuillet, N., Boisson, D., Calais, E., and de Chabalier, J.-B.: Active Tectonics in Southern Haiti and Surface Rupture of the 14 August 2021 Earthquake, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-14568, https://doi.org/10.5194/egusphere-egu23-14568, 2023.

X2.265
|
EGU23-17189
|
ECS
Jeffrey Lang, Joel Baker, Julie Rowland, Adam Hartland, Paul Williams, John Hellstrom, Jamie Howarth, Ingrid Ukstins, Travis Cross, and Christopher Wood

Speleoseismology aims to reconstruct palaeoseismic records by dating pre- and post-damage speleothem calcite. A common approach is to infer palaeo-earthquakes from evidence of coinciding damage features (e.g., rockfall and broken speleothems) at multiple locations, which can be challenging in regions of high tectonic strain where short recurrence intervals of large earthquakes require dating of an impractically large number of damage features. Alternative approaches concerned with dating successive growth changes in individual speleothems (e.g., axis changes and growth hiatuses) are better suited to high-seismicity settings, as closely spaced events are more readily resolved. However, the origins of these growth changes can be ambiguous.

This study tested a novel geochemical proxy for quantifying ground shaking that is amenable to high-resolution speleothem studies, and potentially more diagnostic of earthquake damage. We evaluated the hypothesis that past large earthquakes temporarily elevate Mg/Ca in cave drip waters via incongruent carbonate dissolution following host rock fracturing (ICDC), leading to corresponding Mg enrichments in speleothem calcite. To do this, we examined a well-dated Holocene stalagmite (GT1) from a cave near the Alpine Fault, which is Aotearoa/New Zealand’s longest (>500 km) active onshore fault and a major source of seismic hazard. The locality is 4 km from the Alpine Fault’s northern section, which typically ruptures every 414–470 yr in a major (MW >7) to great (MW >8) earthquake, resulting in shaking intensities of MMI >VIII at the study site (MMI: Modified Mercalli Intensity).

We present a record of Mg/Ca variability in GT1 since ~5 ka, obtained by laser ablation inductively coupled plasma mass spectrometry along the stalagmite growth axis, and constrained temporally by >40 U–Th ages. Preliminary data show high baseline Mg concentrations in GT1 that cannot be explained solely by other mechanisms of drip water Mg/Ca enrichment (i.e., prior calcite precipitation), suggesting an ongoing contribution of Mg to drip waters by ICDC. Anomalous Mg peaks are therefore interpreted as high-intensity shaking events that temporarily elevated drip water Mg/Ca above baseline values. Post-2.5 ka Mg peaks are generally more subtle (30–50% enrichment) than pre-2.5 ka peaks (40–100%). Magnesium peaks are also strongly associated with brown-stained laminae inferred to reflect soil-derived organics. We propose that the high-Mg/high-organics horizons represent large earthquakes that both fractured the host rock and enhanced the mobilisation of organics from overlying soil.

We compared the GT1 record with a proximal and independent 1.4-kyr record of well-dated seismically triggered lacustrine turbidites. Given the subtle nature of Mg peaks in this interval, we consider those associated with physical growth changes (i.e., growth onset/cessation and/or axis change) as more likely to represent earthquakes. Of nine Mg peaks identified, five are associated with major physical growth changes. Of the four largest (MMI >VIII) shaking events in the lake turbidite record, which correspond to northern Alpine Fault surface-rupturing earthquakes, three overlap in age with a GT1 Mg peak and physical growth change. Further, two of the three historic earthquakes that generated MMI ≥VII shaking at the study site also overlap in age with a Mg peak.

How to cite: Lang, J., Baker, J., Rowland, J., Hartland, A., Williams, P., Hellstrom, J., Howarth, J., Ukstins, I., Cross, T., and Wood, C.: Testing a novel cave-based proxy for palaeo-earthquake shaking on the Alpine Fault, Aotearoa/New Zealand., EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-17189, https://doi.org/10.5194/egusphere-egu23-17189, 2023.

X2.266
|
EGU23-13616
|
ECS
Sarah Visage, Pauline Souloumiac, Nadaya Cubas, Bertrand Maillot, and Yann Klinger

Large continental strike-slip earthquakes produce spectacular surface deformations. However, ground displacements are only partially measured in comparison with the amount of slip inferred at depth. Relatively few estimates of the proportion of surface deformation accumulated on faults and deformation distributed regionally around faults are available. However, new technological advances such as state-of-the-art space imaging techniques now greatly improve the quality of surface rupture measurements. Their application has revealed that a significant amount of deformation is accommodated as diffuse deformation in an area of several hundred meters around the fault. This distribution is suggested to depend on the fault complexity. It is therefore essential to understand this distribution and its relation with fault segmentation to study the impact of fault complexities in a seismic context, we recently developed an innovative experimental prototype using some granular materials in layers similar to the earth's crust. They consist of a basal layer of rubber powder that stores elastic energy provided by the displacement of a basal plate sliding parallel to a second, fixed plate. The second layer is made of raw, twice broken rice that brings the stick-slip behaviour required for locking the slip between ruptures, and a third layer of sand with the frictional behaviour of cold shallow sediments at the surface. The surface sand layer allows following the evolution of the fault surface trace from the R-shears stage to the anastomosed fault zone composed of a succession of segments separated by zones of complexities. Using image correlation, we analyse the surface displacements. Since the rice layer causes a stick-slip behaviour, the analysis of the surface displacement is done on several seismic cycles: if the surface displacement is lower than the displacement imposed by the motor, it is an inter-seismic period, if the surface displacement is faster than the displacement imposed by the motor then it is a seismic event.

Once this catalog of events is established, the analysis of the gradient of the displacement Ux parallel to the basal enables us to quantify the deformation: localized (On-Fault Deformation) or distributed (Off-Fault Deformation).

At the R-shear stage, we measure [50~80] % of Off-Fault Deformation. Once the strike-slip fault is formed, the percentage of OFD drops between 20 to 30 %. These results are comparable to measurements made by experiments devoid of a stick-slip behaviours (with only of sand). Moreover, if we compare these values to the proportions of OFD estimated for natural earthquakes, we find the same distribution.

These experiments show that the more mature the fault, the more it will rupture seismically, in time as in space.

How to cite: Visage, S., Souloumiac, P., Cubas, N., Maillot, B., and Klinger, Y.: Evolution of the Off-Fault Deformation during experimental strike-slip earthquakes, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-13616, https://doi.org/10.5194/egusphere-egu23-13616, 2023.

X2.267
|
EGU23-174
|
ECS
Juan Portela, Ian J. Hamling, Alejandra Staller, Marta Béjar-Pizarro, Douglas Hernández, Cecilia Polío López, and Manuel Díaz

The country of El Salvador lies on an active tectonic margin, where the Cocos plate is subducting under the Caribbean plate. A crustal fault system, the El Salvador Fault Zone (ESFZ), crosses the country from East to West through the Central American Volcanic Arc, accommodating more than 1 cm/yr of differential deformation between the Chortís block and the volcanic forearc sliver. 

Here we use GNSS and interferometric synthetic aperture radar (InSAR) data to measure interseismic ground deformation across ESFZ. We have processed and updated GNSS data in more than 110 continuous and episodic stations in the region. GNSS results have been useful for determining the broad pattern of the tectonic signal in the area. However, they are scarce and unable to characterise complex behaviour in the intra-fault basins.

SAR data acquired by the ALOS PALSAR L-band satellite (2006-2011), for both the ascending and descending tracks covering El Salvador, were used to form interferograms with a Small Baseline (SBAS) approach. The time series and average velocity were computed. The average coherence obtained for the area is overall good, and the results are coherent with the regional tectonics. 

How to cite: Portela, J., Hamling, I. J., Staller, A., Béjar-Pizarro, M., Hernández, D., Polío López, C., and Díaz, M.: A new velocity field for El Salvador derived from combined InSAR and GNSS data, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-174, https://doi.org/10.5194/egusphere-egu23-174, 2023.

X2.268
|
EGU23-5555
|
ECS
|
Jill Peikert, Andrea Hampel, and Meike Bagge

Earthquakes on faults in the brittle upper crust cause sudden changes in pore fluid pressure as well as postseismic viscoelastic flow in the lower crust. Such transient processes change the velocity field in the crust and cause significant Coulomb stress changes on receiver faults in the vicinity of the source fault, which may trigger or delay next earthquakes. As previous studies focused on natural earthquakes and/or considered poroelastic and viscoelastic processes separately, the combined influence of poroelastic effects and viscoelastic relaxation on postseismic velocity and stress fields has not been systematically studied so far. In a previous study with 2D finite-element models, we showed that postseismic velocity fields contain signals from overlapping poroelastic and viscoelastic effects (Peikert et al., Tectonophysics, 2022). Here, we use 3D finite-element models with arrays of normal and thrust faults, respectively, to analyze the Coulomb stress changes resulting from the interaction between poroelastic effects and viscoelastic relaxation. In different experiments, we vary the permeability of the crust and the viscosity of the lower crust or lithospheric mantle, while keeping the other parameters constant. We also performed experiments with and without pore fluid flow and viscoelastic relaxation, to isolate the effects of fluid flow and viscoelastic relaxation from each other. Our results show that the coseismic (= static) Coulomb stress changes are immediately altered by the signal from poroelastic effects during the first month after the earthquake. In the first postseismic year, Coulomb stress changes arising from poroelastic effects are one order of magnitude stronger than Coulomb stress changes arising from viscoelastic relaxation. In models considering fluid flow, poroelastic effects dominate the stress field in the first two years. Viscoelastic relaxation already occurs in the early postseismic phase, but is overlapped by the strong signal from poroelastic effects and dominates the Coulomb stress change pattern from about the fifth postseismic year onward for several decades.  The Coulomb stress change patterns show a combined signal from poroelastic and viscoelastic effects already during the first postseismic year, if the viscosity is sufficiently low. For sufficiently low permeabilities, Coulomb stress changes induced by poroelastic effects overlap with the signals from viscoelastic relaxation and interseismic stress accumulation for decades. Finally, poroelastic and viscoelastic effects have a strong impact on the magnitudes and patterns of Coulomb stress changes and should therefore be considered together when analyzing Coulomb stress transfer between faults.

How to cite: Peikert, J., Hampel, A., and Bagge, M.: Relative Importance of Poroelastic Effects and Viscoelastic Relaxation for Co- and Postseismic Coulomb Stress Changes on Normal and Thrust faults: Insights from 3D Finite-Element Modeling, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-5555, https://doi.org/10.5194/egusphere-egu23-5555, 2023.

X2.269
|
EGU23-6332
|
ECS
Taco Broerse, Mario D'Acquisto, Rob Govers, Celine Marsman, and Alireza Amiri-Simkooei

Before geodetically derived strain and rotation rates can be robustly compared to geological or seismological observations, we need reliable strain rate uncertainties. Various methods exist to compute strain rates from GNSS-derived interseismic velocities, but a realistic representation of interpolation uncertainties has remained a challenge. The main problem is that commonly used deterministic interpolation methods do not account for uncertainty resulting from the absence of information in between observation sites. We apply stochastic interpolation by means of ordinary kriging to propagate errors both from discontinuous data coverage as well as from observation uncertainties to our strain rate estimates. However, interseismic horizontal surface velocities in tectonically active regions are spatially highly non-stationary, with high spatial variability around active faults and lower velocity variability in tectonically more stable regions. This requires an extension of traditional ordinary kriging approaches. For interpolation uncertainties that reflect the local variability and spatial correlation of the observed surface velocities, we apply a novel method that incorporates the spatially variable statistics of the underlying data. We estimate realistic uncertainties and covariances of the interpolated velocity field. For regions with a high spatial velocity variability, we find a large increase in uncertainty with increasing distance from observation sites, while in areas with little spatial variability, we estimate a small increase in uncertainty with distance. Subsequently, we propagate interpolated velocity covariance to strain rate uncertainties, such that we can assess the statistical significance of the interpolated strain rate field. Applied to a number of actively deforming regions, including the Pacific coast of North America and Japan, we show to what degree we can robustly determine strain rates based on available GNSS-derived velocities. Realistic uncertainties assist the community to better discriminate continuous or localized deformation on active faults from the available geodetic data. 

 
 

How to cite: Broerse, T., D'Acquisto, M., Govers, R., Marsman, C., and Amiri-Simkooei, A.: Realistic interseismic strain rate uncertainties from inherently sparse GNSS-networks, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-6332, https://doi.org/10.5194/egusphere-egu23-6332, 2023.

X2.270
|
EGU23-7301
Emilie Klein, Hugo Boulze, Christophe Vigny, Luce Fleitout, and Jean-Didier Garaud

Ever since the Maule earthquake (Mw8.8, 2010), a quick vertical uplift is measured thanks to GNSS in the Andes, facing the rupture zone (~250 km to the trench). Models built for the Maule earthquake [Klein et al. 2016] have highlighted that a low-viscosity channel is required to explain the post-seismic uplift. This channel is located along the slab between 50 km and 130 km depth and has a viscosity of a few 1017 Pa.s - lower than in the asthenosphere, 1018 Pa.s. 

After the Illapel earthquake (Mw8.3, 2015), simple observations on GNSS time-series show that no uplift occurred in the Andes at an equivalent distance to the trench than in the case of the Maule earthquake. The subduction in the Illapel region is characterized by a flat-slab (called the Pampean flat-slab) in contrast with the normal-dipping subduction in the region of Maule.

Here, we investigate what is the impact of the Pampean flat-slab on the post-seismic deformations of the Illapel earthquake. In particular, we try to understand  whether the presence of the flat-slab inhibits the effect of the low-viscosity channel. For that purpose we compare GNSS vertical displacements with predictions in both regions of Maule and Illapel from a 3D spherical finite-element model that accounts for the slab geometry of the Chilean subduction zone.

How to cite: Klein, E., Boulze, H., Vigny, C., Fleitout, L., and Garaud, J.-D.: Is the Pampean flat-slab responsible for the differences in post-seismic motions between Maule and Illapel earthquakes?, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-7301, https://doi.org/10.5194/egusphere-egu23-7301, 2023.

X2.271
|
EGU23-8837
|
ECS
|
Tobias Köhne, Rishav Mallick, and Mark Simons

Constraining the effective rheology of major faults is crucial to improve our understanding of the physics of plate boundary deformation. Laboratory studies have used analog experiments to propose rheological models based on viscoelasticity or friction that match laboratory-observed behavior under stress-controlled conditions. Such models have since been used to fit real-world observations of deformation near plate interfaces (both for co- and postseismic displacement timeseries), yielding a variety of estimates of key rheological parameters.
However, confidently differentiating between models using purely observations of a single earthquake (coseismic and postseismic deformation) is difficult — especially in the presence of coarse spatiotemporal sampling, inherent observational noise, and the simplifications of our forward models. In this study, we present a framework built on numerical probabilistic simulations aimed at using displacement timeseries across multiple earthquake cycles in a subduction zone, which successfully distinguishes between endmember constitutive models and recovers key rheological properties. Using synthetic Global Navigation Satellite System network datasets, we furthermore investigate the sensitivity of (hyper-)parameters to the recovery of the true underlying rheological models, and present progress made towards using real 3D observations of a megathrust.

How to cite: Köhne, T., Mallick, R., and Simons, M.: Description Of A Framework And Associated Sensitivity Analysis For Recovery Of Rheological Models And Their Key Parameters Using Multi-Cycle Fault Slip Models, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-8837, https://doi.org/10.5194/egusphere-egu23-8837, 2023.

X2.272
|
EGU23-9752
|
ECS
Alejandra Barreto, Renier Viltres, Rémi Matrau, Benedikt G Ófeigsson, and Sigurjón Jónsson

The Tjörnes Fracture Zone poses significant seismic hazard to the town of Húsavík and other nearby coastal communities in North Iceland as it is capable of generating earthquakes of magnitude ~7. The 120 km long offset connects the offshore Kolbeinsey Ridge to the onshore Northern Volcanic Zone and accommodates approximately 18 mm/yr of transform motion between the North American and Eurasian plates. Most of the deformation is taken up by the two main structures of the fracture zone. The Grímsey Oblique Rift exhibits bookshelf faulting and consists of steeply dipping faults, arranged en-echelon and striking roughly N-S, bounding a series of left-stepping basins. The Húsavík-Flatey Fault is a ~100 km-long right-lateral strike-slip fault. It is mostly offshore, except for its easternmost ~25 km that comes onshore just north of Húsavík. To assess how the deformation is partitioned within the Tjörnes Fracture Zone and to calculate the rate of seismic moment accumulation on the Húsavík-Flatey Fault we use geodetic data from our North Iceland GNSS network. The network covers an area of roughly 200 km by 130 km in size and includes 21 continuous and 92 campaign-style GNSS stations. The continuous data now span up to ~21 years from 2001 to 2022. The first campaign measurements that focused on the HFF were carried out in 1995 and since then we have expanded the campaign-station network to the West towards Tröllaskagi and Skagafjörður and remeasured the network on several occasions. Data from the 2002, 2007, 2009, 2010, 2011, 2013, 2016, 2019, and 2022 campaigns are included in our study. In addition, several stations from the nationwide ISNET reference station network within our study area also included. The GNSS data is used to produce the most up to date velocity field from North Iceland. Relative to the North American plate, our results show a gradual increase of East velocities towards the Northeast across the two main transform structures that reach roughly 18 mm/yr on the Eurasian plate. At the northern tip of the Tjörnes peninsula, between the two transform structures, the velocities are roughly at half the total rate seen at the easternmost stations on the Eurasian plate. Limited deformation is found Southwest of the Húsavík-Flatey Fault in Tröllaskagi, within the so-called Dalvík zone, located on the North American plate.  These results are used to study the present day-kinematics of the Tjörnes Fracture Zone and to further improve the locking depth and slip-rate estimates of the main lineaments.

How to cite: Barreto, A., Viltres, R., Matrau, R., Ófeigsson, B. G., and Jónsson, S.: Interseismic deformation in the Tjörnes Fracture Zone, North Iceland from GNSS measurements, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-9752, https://doi.org/10.5194/egusphere-egu23-9752, 2023.

X2.273
|
EGU23-10710
|
ECS
Joaquín Julve, Marcos Moreno, Sylvain Barbot, Andrés Tassara, Rodolfo Araya, and Nicole Catalán

In the last 20 years, the Chile Subduction Zone (CSZ) has hosted two deep-located subduction events, the 2007 Mw 7.7 Tocopilla earthquake at the Mejillones Peninsula, and the 2016 Mw 7.6 Melinka earthquake at the south of the Chiloé Island. Interseismic seismicity at the Northern and Southern segments of the CSZ, show that in both cases, the ruptures initiated at the down-dip limit of the seismogenic zone. Locking degree models suggest that hypocenter location of this kind of megathrust earthquakes is spatially related with the transition from strongly to weakly locked areas. There are major differences in fault geometry, temperature-pressure regime, petrology at the plate interface and forearc structure between the North and South of the CSZ, raising the question about how such different tectonic settings allow a similar style of rupture. By constructing geologically and geophysically constrained dynamic numerical simulations, here we show that moderate-to-large deep nucleated earthquakes are controlled by petrology and pressure-temperature conditions at the plate interface, along with the structure of the forearc wedge. Our results explain the occurrence, recurrence times and coseismic upper crust deformation of both earthquakes, suggesting that blind ruptures are not only generated at specific conditions, but a suitable combination of the aforementioned parameters is needed. Since the Northern Chile subduction zone has no sediments at the megathrust, the frictional behavior is controlled by altered basalt at the seismogenic depth, and seismicity shows a strong temperature-dependence. Once altered basalt no longer behaves as a velocity weakening material, blueschist rocks allow slow-slip events to develop. The Southern Chile subduction zone is filled with Pliocene-to-present sediments feeding a quartz-dominated subduction channel that defines the seismogenic limit. Within this framework, basal accretion structures are overlapped with a fluid concentration zone at the Moho depth, where the Melinka earthquake initiated. These synoptic views of the CSZ manifest a strong interaction between fluid-rock and forearc structures, which explains the occurrence of blind ruptures at the subduction seismic cycle.

How to cite: Julve, J., Moreno, M., Barbot, S., Tassara, A., Araya, R., and Catalán, N.: The deep subduction earthquake machine: A synoptic view of the Chile Subduction Zone., EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-10710, https://doi.org/10.5194/egusphere-egu23-10710, 2023.

X2.274
|
EGU23-12874
Maureen Llinares, Lucilla Benedetti, Ghislain Gassier, and Sophie Viseur

Measuring 36Cl cosmogenic nuclides on exposed bedrock fault scarps has now been used in several places in the Mediterranean to retrieve ages of the fault seismic history (e.g. Mechernich et al. 2022, Iezzi et al. 2021 and Cowie et al. 2017).

In Central Apennines, around the Fucino basin, at least 15 36Cl sampling sites were analyzed in previous studies to interpret the 36Cl data as seismic history or slip-rates. Several codes (e.g., Beck et al. 2018, Shlagenhauf et al. 2010) were used as a basis for solving 36Cl production equations to calculate the 36Cl concentration resulting from bedrock scarp exhumation history. Some codes included an MCMC routine to retrieve the seismic histories the closest to the dataset. The main differences between the various codes lie in: 1-the fault history prior to exhumation, 2-the parameters previous authors decided to inverse (as an example, mean density of the colluvium is inversed in Beck et al. 2018 but not in Tesson et al. 2019) and 3-the a priori distribution of those parameters (for instance, the time between two earthquakes follows an inverse gaussian distribution for Beck et al. 2018 but a uniform distribution for Tesson et al. 2019). I have compared the various codes and run them on the same dataset (one site at Campo Felice, one site at Roccapreturo and one site at Magnola) and found that retrieved seismic histories are similar, although the estimation of uncertainties differs.

Moreover, all previous cited codes run under Matlab or Fortran. Fortran codes have the advantage of fast computing time but could be cumbersomeI here propose a new code, adapted from Tesson et al. 2019, in the more accessible and widely used Python language. The inferred pre-exposure is also inversed and is a function of the height of the fault cumulative escarpment. The parameters considered are the number of events, ages of event, the associated slips, the long term slip rate, the quiescence and the pre-exposure and their optimal evaluation is done with a MCMC algorithm provided by Numpyro (Du Phan et al. 2019).

Using this new code, we have reanalyzed the 15 36Cl sites around the Fucino and, through a gaussian mixture algorithm, checked the hypothesis of common periods of activity throughout all the Fucino basin.

 

How to cite: Llinares, M., Benedetti, L., Gassier, G., and Viseur, S.: A new python code to invert 36Cl cosmogenic nuclide dataset on normal fault bedrock scarps: comparison with previous published codes and results on the accuracy of the retrieved seismic history of two normal fault systems in Central Apennines, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-12874, https://doi.org/10.5194/egusphere-egu23-12874, 2023.

Posters virtual: Tue, 25 Apr, 16:15–18:00 | vHall TS/EMRP

Chairpersons: Y. Klinger, Magali Rizza
vTE.5
|
EGU23-1195
Mohammad Foroutan, Bertrand Meyer, Michel Sébrier, Andrew Murray, Mohammad-Ali Shokri, Shahryar Solaymani Azad, Hamid Nazari, Faezeh Azhandeh, Ailar Sajedi Far, and Mojtaba Bassiri

While long-averaged recurrence times of large earthquakes are documented on many slow-slipping fault zones in intracontinental settings, the variability of the return periods through multiple seismic cycles remains poorly known. Paleoseismic investigations across fault zones with the well-documented instrumental sequence of surface-breaking earthquakes are a way to tackle the problem. In this context, the Gowk fault, a 160-km-long dextral fault in central Iran, that experienced four surface-rupturing earthquakes with magnitudes ranging from Mw 5.8 to 7.0 during a 1981-1998 earthquake sequence offers a case study. The four earthquakes have ruptured a 90-km-stretch of the fault. The most recent one, the 14 March 1998 Fandoqa earthquake of Mw 6.6, produced a 23-km-long surface rupture along the northern part of the Gowk fault with a maximum right-lateral displacement of 3 m. With a Holocene slip-rate between 3.8-5.7 mm yr-1 and several recent seismic events testifying to its high level of seismicity, the Gowk fault is an appropriate target to conduct paleoseismic investigations and address the earthquake behavior of slow-slipping faults activated by a sequence of well-documented instrumental earthquakes. We excavated two neighboring trenches across the 1998 fault breaks and identified at least four Holocene event horizons that preceded the 1981-1998 earthquake sequence. The age of the faulted stratigraphic sequence is constrained by eighteen optically stimulated luminescence samples and one radiocarbon age on charcoal. The ages of the event horizons suggest an irregular seismic behavior of the Gowk fault characterized by significant variability in the return period of surface rupturing earthquakes.

How to cite: Foroutan, M., Meyer, B., Sébrier, M., Murray, A., Shokri, M.-A., Solaymani Azad, S., Nazari, H., Azhandeh, F., Sajedi Far, A., and Bassiri, M.: Irregular recurrence of surface-faulting paleoearthquakes along the Gowk fault, southeast Iran, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-1195, https://doi.org/10.5194/egusphere-egu23-1195, 2023.