Pub Date : 2024-07-15DOI: 10.26443/seismica.v3i2.1125
Oliver Lamb, S. Bannister, J. Ristau, Craig Miller, S. Sherburn, Katie Jacobs, Jonathan Hanson, E. D'Anastasio, S. Hreinsdóttir, E. Snee, Mike Ross, Eleanor R. H. Mestel, F. Illsley‐Kemp
Taupō is a large caldera volcano located beneath a lake in the centre of the North Island of New Zealand and most recently erupted ~1800 years ago. The volcano has experienced at least 16 periods of unrest since 1872, each of which were characterised by increased seismic activity. Here we detail seismic activity during the most recent period of unrest from May 2022 to May 2023. The unrest was notable for the highest number of earthquakes detected during instrumented unrest episodes, and for one of the largest magnitude earthquakes detected beneath the lake for at least 50 years (ML 5.7). Relocated earthquakes indicate seismic activity was focused around an area hosting overlapping caldera structures and a hydrothermal system. Moment tensor inversion for the largest earthquake includes a non-negligible inflationary isotropic component. We suggest the seismic unrest was caused by the reactivation of faults due to an intrusion of magma at depth.
{"title":"Seismic characteristics of the 2022-2023 unrest episode at Taupō volcano, Aotearoa New Zealand","authors":"Oliver Lamb, S. Bannister, J. Ristau, Craig Miller, S. Sherburn, Katie Jacobs, Jonathan Hanson, E. D'Anastasio, S. Hreinsdóttir, E. Snee, Mike Ross, Eleanor R. H. Mestel, F. Illsley‐Kemp","doi":"10.26443/seismica.v3i2.1125","DOIUrl":"https://doi.org/10.26443/seismica.v3i2.1125","url":null,"abstract":"Taupō is a large caldera volcano located beneath a lake in the centre of the North Island of New Zealand and most recently erupted ~1800 years ago. The volcano has experienced at least 16 periods of unrest since 1872, each of which were characterised by increased seismic activity. Here we detail seismic activity during the most recent period of unrest from May 2022 to May 2023. The unrest was notable for the highest number of earthquakes detected during instrumented unrest episodes, and for one of the largest magnitude earthquakes detected beneath the lake for at least 50 years (ML 5.7). Relocated earthquakes indicate seismic activity was focused around an area hosting overlapping caldera structures and a hydrothermal system. Moment tensor inversion for the largest earthquake includes a non-negligible inflationary isotropic component. We suggest the seismic unrest was caused by the reactivation of faults due to an intrusion of magma at depth.","PeriodicalId":509514,"journal":{"name":"Seismica","volume":"34 35","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-07-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141645411","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-07-05DOI: 10.26443/seismica.v3i2.1202
Dannielle Fougere, James Dolan, Edward Rhodes, Sally McGill
We use lidar- and field-based mapping coupled with single-grain infrared-stimulated luminescence dating to constrain three new slip rate estimates from the western and central segments of the Garlock fault in southern California, revealing a more complete picture of incremental slip rate in time and space for this major plate-boundary fault. These new rates reinforce and refine previous evidence showing that the Garlock fault experiences significant temporal variations in slip rates that span multiple earthquake cycles, with multi-millennial periods of very fast (13-14 mm/yr) early and late Holocene slip separated by a mid-Holocene period of slow slip (3 mm/yr). Similar ca. 8 ka slip rates for the central Garlock fault of 8.8 ± 1.0 mm/yr and 8.2 +1.0/-0.8 mm/yr for the western Garlock fault demonstrate that the fault has slipped at a faster long-term average rate than suggested by previous studies. These fast rates are consistent with kinematic models in which the western and central Garlock fault segments are driven primarily by lateral extrusion associated with N-S contractional shortening, with additional slip driven by WNW-ENE Basin and Range extension north of the fault and minor rotation of the Garlock within the N-S zone of dextral ECSZ shear.
{"title":"Refined Holocene Slip Rate for the Western and Central Segments of the Garlock Fault: Record of Alternating Millennial-Scale Periods of Fast and Slow Fault Slip","authors":"Dannielle Fougere, James Dolan, Edward Rhodes, Sally McGill","doi":"10.26443/seismica.v3i2.1202","DOIUrl":"https://doi.org/10.26443/seismica.v3i2.1202","url":null,"abstract":"We use lidar- and field-based mapping coupled with single-grain infrared-stimulated luminescence dating to constrain three new slip rate estimates from the western and central segments of the Garlock fault in southern California, revealing a more complete picture of incremental slip rate in time and space for this major plate-boundary fault. These new rates reinforce and refine previous evidence showing that the Garlock fault experiences significant temporal variations in slip rates that span multiple earthquake cycles, with multi-millennial periods of very fast (13-14 mm/yr) early and late Holocene slip separated by a mid-Holocene period of slow slip (3 mm/yr). Similar ca. 8 ka slip rates for the central Garlock fault of 8.8 ± 1.0 mm/yr and 8.2 +1.0/-0.8 mm/yr for the western Garlock fault demonstrate that the fault has slipped at a faster long-term average rate than suggested by previous studies. These fast rates are consistent with kinematic models in which the western and central Garlock fault segments are driven primarily by lateral extrusion associated with N-S contractional shortening, with additional slip driven by WNW-ENE Basin and Range extension north of the fault and minor rotation of the Garlock within the N-S zone of dextral ECSZ shear.","PeriodicalId":509514,"journal":{"name":"Seismica","volume":" 7","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-07-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141676976","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-07-02DOI: 10.26443/seismica.v2i4.1177
N. Harrichhausen, K. Morell, C. Regalla
We employ numerical models to explore the connection between subduction zone coupling or megathrust rupture and upper plate faults in the northern Cascadia forearc. Active forearc faults north of the Olympic Peninsula exhibit similar characteristics: west-northwest strike, oblique right-lateral slip senses, and low slip rates (<1 mm/yr), but a potential to generate large (M ~ 7) earthquakes. Previous hypotheses suggest that stress in the upper plate due to interseismic coupling or coseismic rupture along the subduction zone interface could drive permanent forearc strain. To test these hypotheses, we used a 3D boundary element method model to predict slip on the LRDM if interseismic coupling or coeseismic rupture cause deformation. Our model predicts reverse left-lateral slip if the strain results solely from subduction zone coupling, or normal right-lateral slip if these faults accommodate strain during a megathrust rupture. These results contradict the observed fault kinematics. Additionally, if we use our model to mimic strain partitioning, where only the strain from the strike-slip component of subduction zone coupling is accommodated in the forearc, our results are also inconsistent observed fault kinematics. These models challenge the hypothesis that subduction zone coupling or coseismic rupture are the primary driver of permanent forearc deformation in northern Cascadia.
{"title":"Forearc faults in northern Cascadia do not accommodate elastic strain driven by the megathrust seismic cycle","authors":"N. Harrichhausen, K. Morell, C. Regalla","doi":"10.26443/seismica.v2i4.1177","DOIUrl":"https://doi.org/10.26443/seismica.v2i4.1177","url":null,"abstract":"We employ numerical models to explore the connection between subduction zone coupling or megathrust rupture and upper plate faults in the northern Cascadia forearc. Active forearc faults north of the Olympic Peninsula exhibit similar characteristics: west-northwest strike, oblique right-lateral slip senses, and low slip rates (<1 mm/yr), but a potential to generate large (M ~ 7) earthquakes. Previous hypotheses suggest that stress in the upper plate due to interseismic coupling or coseismic rupture along the subduction zone interface could drive permanent forearc strain. To test these hypotheses, we used a 3D boundary element method model to predict slip on the LRDM if interseismic coupling or coeseismic rupture cause deformation. Our model predicts reverse left-lateral slip if the strain results solely from subduction zone coupling, or normal right-lateral slip if these faults accommodate strain during a megathrust rupture. These results contradict the observed fault kinematics. Additionally, if we use our model to mimic strain partitioning, where only the strain from the strike-slip component of subduction zone coupling is accommodated in the forearc, our results are also inconsistent observed fault kinematics. These models challenge the hypothesis that subduction zone coupling or coseismic rupture are the primary driver of permanent forearc deformation in northern Cascadia.","PeriodicalId":509514,"journal":{"name":"Seismica","volume":"63 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141687089","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-06-10DOI: 10.26443/seismica.v3i1.980
J. Jara, R. Jolivet, A. Socquet, D. Comte, E. Norabuena
Detections of slow slip events (SSEs) are now common along most plate boundary fault systems at the global scale. However, no such event has been described in the south Peru - north Chile subduction zone so far, except for the early preparatory phase of the 2014 Iquique earthquake. We use geodetic template matching on GNSS-derived time series of surface motion in Northern Chile to extract SSEs hidden within the geodetic noise. We detect 33 events with durations ranging from 9 to 40 days and magnitudes from Mw 5.6 to 6.2. The moment released by these aseismic events seems to scale with the cube of their duration, suggesting a dynamic comparable to that of earthquakes. We compare the distribution of SSEs with the distribution of coupling along the megathrust derived using Bayesian inference on GNSS- and InSAR-derived interseismic velocities. From this comparison, we obtain that most SSEs occur in regions of intermediate coupling where the megathrust transitions from locked to creeping or where geometrical complexities of the interplate region have been proposed. We finally discuss the potential role of fluids as a triggering mechanism for SSEs in the area.
{"title":"Detection of slow slip events along the southern Peru - northern Chile subduction zone","authors":"J. Jara, R. Jolivet, A. Socquet, D. Comte, E. Norabuena","doi":"10.26443/seismica.v3i1.980","DOIUrl":"https://doi.org/10.26443/seismica.v3i1.980","url":null,"abstract":"Detections of slow slip events (SSEs) are now common along most plate boundary fault systems at the global scale. However, no such event has been described in the south Peru - north Chile subduction zone so far, except for the early preparatory phase of the 2014 Iquique earthquake. We use geodetic template matching on GNSS-derived time series of surface motion in Northern Chile to extract SSEs hidden within the geodetic noise. We detect 33 events with durations ranging from 9 to 40 days and magnitudes from Mw 5.6 to 6.2. The moment released by these aseismic events seems to scale with the cube of their duration, suggesting a dynamic comparable to that of earthquakes. We compare the distribution of SSEs with the distribution of coupling along the megathrust derived using Bayesian inference on GNSS- and InSAR-derived interseismic velocities. From this comparison, we obtain that most SSEs occur in regions of intermediate coupling where the megathrust transitions from locked to creeping or where geometrical complexities of the interplate region have been proposed. We finally discuss the potential role of fluids as a triggering mechanism for SSEs in the area.","PeriodicalId":509514,"journal":{"name":"Seismica","volume":" 48","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141365090","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-06-07DOI: 10.26443/seismica.v2i4.1097
C. Nuyen, David Schmidt
Episodic tremor and slip (ETS) is well-documented along the entire length of the Cascadia subduction zone. We explore how the occurrence of ETS varies at the southernmost edge of the subduction zone, where geometric complexity and a slab window likely alter conditions along the plate interface. This work uses tremor and GNSS time series data to identify nineteen of the largest ETS events in southern Cascadia between 2016.5-2022 and document source properties for events approaching the slab edge. Distributed slip models for these events show that cumulative fault slip along the megathrust reaches a maximum near 40.5° N latitude and that large ETS events accommodate up to 85% of plate convergence at this location. However, ETS fault slip and tremor terminate near 40° N latitude, some 50 km before the southern lateral edge of the subducting plate. After considering a range of explanations, we propose that the complex geometry and progressive heating of the subducting plate modifies ETS behavior and does not allow seismic slip to occur along the plate interface in southernmost Cascadia below 35 km depth.
{"title":"Along-strike changes in ETS behavior near the slab edge of Southern Cascadia","authors":"C. Nuyen, David Schmidt","doi":"10.26443/seismica.v2i4.1097","DOIUrl":"https://doi.org/10.26443/seismica.v2i4.1097","url":null,"abstract":"Episodic tremor and slip (ETS) is well-documented along the entire length of the Cascadia subduction zone. We explore how the occurrence of ETS varies at the southernmost edge of the subduction zone, where geometric complexity and a slab window likely alter conditions along the plate interface. This work uses tremor and GNSS time series data to identify nineteen of the largest ETS events in southern Cascadia between 2016.5-2022 and document source properties for events approaching the slab edge. Distributed slip models for these events show that cumulative fault slip along the megathrust reaches a maximum near 40.5° N latitude and that large ETS events accommodate up to 85% of plate convergence at this location. However, ETS fault slip and tremor terminate near 40° N latitude, some 50 km before the southern lateral edge of the subducting plate. After considering a range of explanations, we propose that the complex geometry and progressive heating of the subducting plate modifies ETS behavior and does not allow seismic slip to occur along the plate interface in southernmost Cascadia below 35 km depth.","PeriodicalId":509514,"journal":{"name":"Seismica","volume":" 16","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-06-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141373515","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-06-05DOI: 10.26443/seismica.v3i1.1186
Guy Salomon, Theron Finley, Edwin Nissen, Roger Stephen, Brian Menounos
The advent of sub-meter resolution topographic surveying has revolutionized active fault mapping. Light detection and ranging (lidar) collected using crewed airborne laser scanning (ALS) can provide ground coverage of entire fault systems but is expensive, while Structure-from-Motion (SfM) photogrammetry from uncrewed aerial vehicles (UAVs) is popular for mapping smaller sites but cannot image beneath vegetation. Here, we present a new UAV laser scanning (ULS) system which overcomes these limitations to survey fault-related topography cost-effectively, at desirable spatial resolutions, and even beneath dense vegetation. In describing our system, data acquisition and processing workflows, we provide a practical guide for other researchers interested in developing their own ULS capabilities. We showcase ULS data collected over faults from a variety of terrain and vegetation types across the Canadian Cordillera and compare them to conventional ALS and SfM data. Due to the lower, slower UAV flights, ULS offers improved ground return density (~260 points/m2 for the capture of a paleoseismic trenching site and ~10–72 points/m2 for larger, multi-kilometer fault surveys) over conventional ALS (~3–9 points/m2) as well as better vegetation penetration than both ALS and SfM. The resulting ~20–50 cm-resolution ULS terrain models reveal fine-scale tectonic landforms that would otherwise be challenging to image.
{"title":"Mapping fault geomorphology with drone-based lidar","authors":"Guy Salomon, Theron Finley, Edwin Nissen, Roger Stephen, Brian Menounos","doi":"10.26443/seismica.v3i1.1186","DOIUrl":"https://doi.org/10.26443/seismica.v3i1.1186","url":null,"abstract":"The advent of sub-meter resolution topographic surveying has revolutionized active fault mapping. Light detection and ranging (lidar) collected using crewed airborne laser scanning (ALS) can provide ground coverage of entire fault systems but is expensive, while Structure-from-Motion (SfM) photogrammetry from uncrewed aerial vehicles (UAVs) is popular for mapping smaller sites but cannot image beneath vegetation. Here, we present a new UAV laser scanning (ULS) system which overcomes these limitations to survey fault-related topography cost-effectively, at desirable spatial resolutions, and even beneath dense vegetation. In describing our system, data acquisition and processing workflows, we provide a practical guide for other researchers interested in developing their own ULS capabilities. We showcase ULS data collected over faults from a variety of terrain and vegetation types across the Canadian Cordillera and compare them to conventional ALS and SfM data. Due to the lower, slower UAV flights, ULS offers improved ground return density (~260 points/m2 for the capture of a paleoseismic trenching site and ~10–72 points/m2 for larger, multi-kilometer fault surveys) over conventional ALS (~3–9 points/m2) as well as better vegetation penetration than both ALS and SfM. The resulting ~20–50 cm-resolution ULS terrain models reveal fine-scale tectonic landforms that would otherwise be challenging to image.","PeriodicalId":509514,"journal":{"name":"Seismica","volume":"61 6","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-06-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141381573","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-05-15DOI: 10.26443/seismica.v3i1.1098
J. Asplet, J. Wookey, Micheal Kendall, Mark Chapman, Ritima Das
The behaviour of fluids in preferentially aligned fractures plays an important role in a range of dynamic processes within the Earth. In the near-surface, understanding systems of fluid-filled fractures is crucial for applications such as geothermal energy production, monitoring CO2 storage sites, and exploration for metalliferous sub-volcanic brines. Mantle melting is a key geodynamic process, exerting control over its composition and dynamic processes. Upper mantle melting weakens the lithosphere, facilitating rifting and other surface expressions of tectonic processes.Aligned fluid-filled fractures are an efficient mechanism for seismic velocity anisotropy, requiring very low volume fractions, but such rock physics models also predict significant shear-wave attenuation anisotropy. In comparison, the attenuation anisotropy expected for crystal preferred orietation mechanisms is negligible or would only operate outside of the seismic frequency band.Here we demonstrate a new method for measuring shear-wave attenuation anisotropy, apply it to synthetic examples, and make the first measurements of SKS attenuation anisotropy using data recorded at the station FURI, in Ethiopia. At FURI we measure attenuation anisotropy where the fast shear-wave has been more attenuated than the slow shear-wave. This can be explained by the presence of aligned fluids, most probably melts, in the upper mantle using a poroelastic squirt flow model. Modelling of this result suggests that a 1% melt fraction, hosted in aligned fractures dipping ca. 40° that strike perpendicular to the Main Ethiopian Rift, is required to explain the observed attenuation anisotropy. This agrees with previous SKS shear-wave splitting analysis which suggested a 1% melt fraction beneath FURI. The interpreted fracture strike and dip, however, disagrees with previous work in the region which interprets sub-vertical melt inclusions aligned parallel to the Main Ethiopian Rift which only produce attenuation anisotropy where the slow shear-wave is more attenuated. These results show that attenuation anisotropy could be a useful tool for detecting mantle melt, and may offer strong constraints on the extent and orientation of melt inclusions which cannot be achieved from seismic velocity anisotropy alone.
流体在优先排列的裂缝中的行为在地球内部的一系列动态过程中发挥着重要作用。在近地表,了解充满流体的裂缝系统对于地热能源生产、监测二氧化碳封存地点以及勘探火山下卤水等应用至关重要。地幔熔化是一个关键的地球动力学过程,对地幔的组成和动态过程具有控制作用。上地幔熔化削弱了岩石圈,促进了断裂和其他构造过程的地表表现形式。排列整齐的充满流体的裂缝是地震速度各向异性的有效机制,所需的体积分数非常低,但此类岩石物理模型也预测了显著的剪切波衰减各向异性。在这里,我们展示了一种测量剪切波衰减各向异性的新方法,并将其应用于合成实例,利用在埃塞俄比亚 FURI 站记录的数据首次测量了 SKS 衰减各向异性。在 FURI 站,我们测量到快速剪切波比慢速剪切波衰减更多的衰减各向异性。这可以用上地幔中存在排列整齐的流体(很可能是熔体)来解释,使用的模型是孔弹性喷流模型。对这一结果的建模表明,要解释所观测到的衰减各向异性,需要有1%的熔体成分,寄存在垂直于埃塞俄比亚主裂谷的倾角约为40°的排列断裂中。这与之前的 SKS 剪切波分裂分析一致,该分析表明 FURI 地下有 1%的熔融部分。然而,对断裂走向和倾角的解释与该地区以前的研究不一致,以前的研究解释了与埃塞俄比亚主裂谷平行的次垂直熔融包裹体,而这种包裹体只在慢剪切波衰减较大的地方产生衰减各向异性。这些结果表明,衰减各向异性可作为探测地幔熔体的有用工具,并可对熔体包裹体的范围和方位提供有力的约束,而这一点仅靠地震速度各向异性是无法实现的。
{"title":"Shear-wave attenuation anisotropy: a new constraint on mantle melt near the Main Ethiopian Rift","authors":"J. Asplet, J. Wookey, Micheal Kendall, Mark Chapman, Ritima Das","doi":"10.26443/seismica.v3i1.1098","DOIUrl":"https://doi.org/10.26443/seismica.v3i1.1098","url":null,"abstract":"The behaviour of fluids in preferentially aligned fractures plays an important role in a range of dynamic processes within the Earth. In the near-surface, understanding systems of fluid-filled fractures is crucial for applications such as geothermal energy production, monitoring CO2 storage sites, and exploration for metalliferous sub-volcanic brines. Mantle melting is a key geodynamic process, exerting control over its composition and dynamic processes. Upper mantle melting weakens the lithosphere, facilitating rifting and other surface expressions of tectonic processes.Aligned fluid-filled fractures are an efficient mechanism for seismic velocity anisotropy, requiring very low volume fractions, but such rock physics models also predict significant shear-wave attenuation anisotropy. In comparison, the attenuation anisotropy expected for crystal preferred orietation mechanisms is negligible or would only operate outside of the seismic frequency band.Here we demonstrate a new method for measuring shear-wave attenuation anisotropy, apply it to synthetic examples, and make the first measurements of SKS attenuation anisotropy using data recorded at the station FURI, in Ethiopia. At FURI we measure attenuation anisotropy where the fast shear-wave has been more attenuated than the slow shear-wave. This can be explained by the presence of aligned fluids, most probably melts, in the upper mantle using a poroelastic squirt flow model. Modelling of this result suggests that a 1% melt fraction, hosted in aligned fractures dipping ca. 40° that strike perpendicular to the Main Ethiopian Rift, is required to explain the observed attenuation anisotropy. This agrees with previous SKS shear-wave splitting analysis which suggested a 1% melt fraction beneath FURI. The interpreted fracture strike and dip, however, disagrees with previous work in the region which interprets sub-vertical melt inclusions aligned parallel to the Main Ethiopian Rift which only produce attenuation anisotropy where the slow shear-wave is more attenuated. These results show that attenuation anisotropy could be a useful tool for detecting mantle melt, and may offer strong constraints on the extent and orientation of melt inclusions which cannot be achieved from seismic velocity anisotropy alone.","PeriodicalId":509514,"journal":{"name":"Seismica","volume":"69 4","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140971658","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-05-09DOI: 10.26443/seismica.v3i1.1142
Dario Jozinović, John Clinton, F. Massin, Maren Böse, C. Cauzzi
It is increasingly common for seismic networks to operate multiple independent automatic algorithms to characterise earthquakes in real-time, such as in earthquake early warning (EEW) or even standard network practice. Commonly used methods to select the best solution at a given time are simple and use ad hoc rules. An absolute measure of how well a solution (event origin and magnitude) matches the observations by the goodness-of-fit between the observed and predicted envelopes is a robust and independent metric to select optimal solutions. We propose such a measure that is calculated as a combination of amplitude and cross-correlation fit. This metric can be used to determine when a preferred solution reaches an appropriate confidence level for alerting, or indeed to compare two (or more) different event characterisations directly. We demonstrate that our approach can also be used to suppress false alarms commonly seen at seismic networks. Tests using the 10 largest earthquakes in Switzerland between 2013 and 2020, and events that caused false alarms demonstrate that our approach can effectively prefer solutions with small errors in location and magnitude, and can clearly identify and discard false origins or incorrect magnitudes, at all time scales, starting with the first event characterisation.
{"title":"Realtime Selection of Optimal Source Parameters Using Ground Motion Envelopes","authors":"Dario Jozinović, John Clinton, F. Massin, Maren Böse, C. Cauzzi","doi":"10.26443/seismica.v3i1.1142","DOIUrl":"https://doi.org/10.26443/seismica.v3i1.1142","url":null,"abstract":"It is increasingly common for seismic networks to operate multiple independent automatic algorithms to characterise earthquakes in real-time, such as in earthquake early warning (EEW) or even standard network practice. Commonly used methods to select the best solution at a given time are simple and use ad hoc rules. An absolute measure of how well a solution (event origin and magnitude) matches the observations by the goodness-of-fit between the observed and predicted envelopes is a robust and independent metric to select optimal solutions. We propose such a measure that is calculated as a combination of amplitude and cross-correlation fit. This metric can be used to determine when a preferred solution reaches an appropriate confidence level for alerting, or indeed to compare two (or more) different event characterisations directly. We demonstrate that our approach can also be used to suppress false alarms commonly seen at seismic networks. Tests using the 10 largest earthquakes in Switzerland between 2013 and 2020, and events that caused false alarms demonstrate that our approach can effectively prefer solutions with small errors in location and magnitude, and can clearly identify and discard false origins or incorrect magnitudes, at all time scales, starting with the first event characterisation.","PeriodicalId":509514,"journal":{"name":"Seismica","volume":" 8","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140997371","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-05-08DOI: 10.26443/seismica.v3i1.1145
P. Bonatis, Vasileios Karakostas, C. Kourouklas, A. Kostoglou, E. Papadimitriou
The two moderate earthquakes that occurred close and to the north of the North Aegean Trough (NAT) on 26 September 2020 (Mw5.3) and 16 January 2022 (Mw5.4), both followed by aftershock activity, are examined. Seismic activity along the NAT and its parallel branches is continuous and remarkable, with numerous strong instrumental (M≥6.0) earthquakes. Yet, the frequency of moderate (5.0≤M<6.0) earthquakes outside these major fault branches is rather rare and therefore their investigation provides the optimal means to decipher the seismotectonic properties of the broader area. The temporal and spatial proximity of the two seismic excitations from late September of 2020 through early 2022, intrigues for exhaustive investigation of seismic activity with the employment of earthquake relocation techniques, moment tensor solutions and statistical analysis. Our research revealed that this seismic activity purely falls inside the Mainshock – Aftershock type, with fast aftershock decay rates and moderate productivity. According to our findings, the two seismic sequences, despite their close proximity, exhibit distinctive features as a result of the intricate stress field generated at the western termination of the NAF system in an extensional domain.
{"title":"Spatiotemporal characteristics and earthquake statistics of the 2020 and 2022 adjacent earthquake sequences in North Aegean Sea (Greece)","authors":"P. Bonatis, Vasileios Karakostas, C. Kourouklas, A. Kostoglou, E. Papadimitriou","doi":"10.26443/seismica.v3i1.1145","DOIUrl":"https://doi.org/10.26443/seismica.v3i1.1145","url":null,"abstract":"The two moderate earthquakes that occurred close and to the north of the North Aegean Trough (NAT) on 26 September 2020 (Mw5.3) and 16 January 2022 (Mw5.4), both followed by aftershock activity, are examined. Seismic activity along the NAT and its parallel branches is continuous and remarkable, with numerous strong instrumental (M≥6.0) earthquakes. Yet, the frequency of moderate (5.0≤M<6.0) earthquakes outside these major fault branches is rather rare and therefore their investigation provides the optimal means to decipher the seismotectonic properties of the broader area. The temporal and spatial proximity of the two seismic excitations from late September of 2020 through early 2022, intrigues for exhaustive investigation of seismic activity with the employment of earthquake relocation techniques, moment tensor solutions and statistical analysis. Our research revealed that this seismic activity purely falls inside the Mainshock – Aftershock type, with fast aftershock decay rates and moderate productivity. According to our findings, the two seismic sequences, despite their close proximity, exhibit distinctive features as a result of the intricate stress field generated at the western termination of the NAF system in an extensional domain.","PeriodicalId":509514,"journal":{"name":"Seismica","volume":" 6","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140997870","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-05-03DOI: 10.26443/seismica.v3i1.887
O. Fadugba, V. Sahakian, D. Melgar, A. Rodgers, Roey Shimony
Accurately modeling time-dependent coseismic crustal deformation as observed on high-rate Global Navigation Satellite System (HR-GNSS) lends insight into earthquake source processes and improves local earthquake and tsunami early warning algorithms. Currently, time-dependent crustal deformation modeling relies most frequently on simplified 1D radially symmetric Earth models. However, for shallow subduction zone earthquakes, even low-frequency shaking is likely affected by the many strongly heterogeneous structures such as the subducting slab, mantle wedge, and the overlying crustal structure. We demonstrate that including 3D structure improves the estimation of key features of coseismic HR-GNSS time series, such as the peak ground displacement (PGD), the time to PGD (tPGD), static displacements (SD), and waveform cross-correlation values. We computed synthetic 1D and 3D, 0.25 Hz and 0.5 Hz waveforms at HR-GNSS stations for four M7.3+ earthquakes in Japan using MudPy and SW4, respectively. From these synthetics, we computed intensity-measure residuals between the synthetic and observed GNSS waveforms. Comparing 1D and 3D residuals, we observed that the 3D simulations show better fits to the PGD and SD in the observed waveforms than the 1D simulations for both 0.25 Hz and 0.5 Hz simulations. We find that the reduction in PGD residuals in the 3D simulations is a combined effect of both shallow and deep 3D structures; hence incorporating only the upper 30 km of 3D structure will still improve the fit to the observed PGD values. Our results demonstrate that 3D simulations significantly improve models of GNSS waveform characteristics and will not only help understand the underlying processes, but also improve local tsunami warning.
{"title":"The Impact of the Three-Dimensional Structure of a Subduction Zone on Time-dependent Crustal Deformation Measured by HR-GNSS","authors":"O. Fadugba, V. Sahakian, D. Melgar, A. Rodgers, Roey Shimony","doi":"10.26443/seismica.v3i1.887","DOIUrl":"https://doi.org/10.26443/seismica.v3i1.887","url":null,"abstract":"Accurately modeling time-dependent coseismic crustal deformation as observed on high-rate Global Navigation Satellite System (HR-GNSS) lends insight into earthquake source processes and improves local earthquake and tsunami early warning algorithms. Currently, time-dependent crustal deformation modeling relies most frequently on simplified 1D radially symmetric Earth models. However, for shallow subduction zone earthquakes, even low-frequency shaking is likely affected by the many strongly heterogeneous structures such as the subducting slab, mantle wedge, and the overlying crustal structure. We demonstrate that including 3D structure improves the estimation of key features of coseismic HR-GNSS time series, such as the peak ground displacement (PGD), the time to PGD (tPGD), static displacements (SD), and waveform cross-correlation values. We computed synthetic 1D and 3D, 0.25 Hz and 0.5 Hz waveforms at HR-GNSS stations for four M7.3+ earthquakes in Japan using MudPy and SW4, respectively. From these synthetics, we computed intensity-measure residuals between the synthetic and observed GNSS waveforms. Comparing 1D and 3D residuals, we observed that the 3D simulations show better fits to the PGD and SD in the observed waveforms than the 1D simulations for both 0.25 Hz and 0.5 Hz simulations. We find that the reduction in PGD residuals in the 3D simulations is a combined effect of both shallow and deep 3D structures; hence incorporating only the upper 30 km of 3D structure will still improve the fit to the observed PGD values. Our results demonstrate that 3D simulations significantly improve models of GNSS waveform characteristics and will not only help understand the underlying processes, but also improve local tsunami warning.","PeriodicalId":509514,"journal":{"name":"Seismica","volume":"77 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141015961","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}