Seismic Observations Research Articles

SEISMICITY of the CARPATHIANS in 2018–2019

The article describes seismic observations in the Carpathian region in 2018–2019, which were carried out, as before, by two organizations from two states: in Ukraine – the Seismicity department of the Carpathian region of the Institute of Geophysics of the National Academy of Sciences of Ukraine, in Moldova – the seismology laboratory of the Institute of Geology and Seismology of the Academy of Sciences of Moldova. The monitoring of seismicity was carried out 21 stationary digital and 3 temporary stations of Ukraine with a processing centre in Lviv and six stations of Moldova with a processing centre in Chisinau. The consolidated catalog of earthquakes was created in Lviv. The total number of Carpathian earthquakes in 2018–2019 was NΣ=99 in the range KR=4.5–13.6 within the interval of hypocenter depths h=0.4–157.5 km, with total seismic energy ΣE=4.2·1013 J. The maximum earthquake with KR=13.6, MSH=5.3 and hрР=150.6 km was registered on October 28, 2018 in the Vrancea zone. The epicentral zone was subjected to shaking with an intensity up to Ii=V by MM (Mercally) scale. The earthquake was felt by the population of Romania, Bulgaria, Ukraine, Moldova and other European countries. In the Precarpathian and Transcarpathian regions, all earthquakes were weaker. For Transcarpathia, their energy class was in the range of KR=5.5–7.9, and the depths was in the range of h=0.6–8.0 km. The strongest event in Precarpathian was a tangible earthquake that occurred in the Drohobych district of the Lviv region on September 27, 2018 with KR=7.3 at a depth of 1.7±0.9 km. The strongest event in Bukovina was an earthquake of March 24, 2018 with KR=9.4, МSH=2.4 and h=1 km near Novodnestrovsk city, Ukraina, which was felt by the population of the village of Lomachintsy with an intensity of 4 points. In Novodnestrovsk, the earthquake was felt only on the upper floors. In general, in 2018–2019 there was a slight decrease in the level of seismicity in seismically active areas: № 1 (North-Western), № 2 (Vrancea), № 5 (Bukovina) and № 8 (Bacau), compared with that in the previous period 2016–2017. In areas № 3 (South Carpathians), № 4 (Banat) and № 7 (Transylvania) there was no earthquakes, and in area № 6 (Krishana) only there was one earthquake with КР=10.8, MSH=3.2.

Plumbing System Architecture and Differentiation Processes of the Nyiragongo Volcano, DR Congo

ABSTRACT The Nyiragongo volcano is one of the most alkali-rich volcanic centres on the planet (Na2O + K2O generally >10 wt.%, agpaitic index up to 1.34), characterized by a semi-permanently active lava lake which hosts silica-undersaturated (SiO2 < 40 wt.%), low viscosity lavas. To improve our understanding of this unique magmatic system, we present a set of 291 samples, acquired during new field excursions between 2017 and 2021. The major and trace element composition of all samples was measured, revealing a lithological range extending from primitive picrites (Mg# 82) erupted from parasitic cones to a variety of highly evolved nephelinites, leucitites, and melilitites erupted from the main edifice as recently as 2002, 2016, and 2021. We measured major and trace element compositions from the full spectrum of minerals present in all sampled lithologies. From these we calculated that the main magma reservoirs feeding Nyiragongo are at approximately 9–15 and 21–33 km depth, in agreement with recent seismic observations. Fractional crystallization modelling using observed mineral compositions and proportions was performed to quantitatively link the lithologies to specific residual liquid fractions assuming evolution from an olivine-melilite parental melt. Our modelling indicates that fractionation and cumulate formation in deep chambers reduces the remaining melt fraction to ~60%, after which melts are injected into upper, liquid dominated magma chambers where fractionation and accumulation of clinopyroxene, melilite, and feldspathoids dominates. Characterisation of mineral textures and geochemistry reveals high crystal mobility in a repeatedly recharging plumbing system split between liquid-dominated, evolved magma chambers and more solid-dominated, primitive mushes, decreasing in liquid fraction with depth.

GEOSCOPE Network: 40 Yr of Global Broadband Seismic Data

Abstract The GEOSCOPE observatory (Institut de physique du globe de Paris [IPGP] and École et Observatoire des Sciences de la Terre de Strasbourg, 1982) provides more than four decades of high-quality continuous broadband data to the scientific community. Started in 1982 with only two stations, the network has grown over the years thanks to numerous international partnerships. At present, 34 stations operate in 18 countries across all continents and on islands throughout the oceans, filling important gaps in global Earth coverage. Most of the first installed stations are still running today, allowing for long-term observations, and new sites are being prospected to further improve global coverage. Over the years, GEOSCOPE has contributed to defining today’s global seismology standards (data format, data quality level, instrumentation requirements), being the French contribution to the international effort for global seismic observations. The stations are instrumented with the best quality seismometers (from the very first STS-1 in the early 80s to the last STS-6A and Trillium T360 today) and digitizers (Q330HR and Centaur) to record with high fidelity the ground motions generated by all types of seismic sources. Real-time data are sent to the tsunami warning centers and both validated and real-time data are available at the IPGP, Epos-France and Earthscope data centers. The quality of GEOSCOPE data and metadata is ensured by daily and yearly validation that enables issue detection and mitigation. GEOSCOPE, in collaboration with the other global networks, has played and continues to play a crucial role in the study of Earth’s structure and global dynamics and the characterization of all types of seismic sources.

Present-day movement characteristics of the Qinghai Nanshan fault and its surrounding area from GPS observation

The Qinghai Nanshan fault is a larger fault in the Northeastern Xizang Plateau. In previous studies, its movement characteristics are mainly investigated with geological and seismic observations, and the tectonic transformation role of the fault on its east is not yet clear. This study uses data fusion to obtain denser GPS observations near the Qinghai Nanshan fault. Based on tectonic characteristics, we establish a block model to investigate the fault slip rate, locking degree, and slip deficit. The results show that the Qinghai Nanshan fault slip rate is characterized by sinistral and convergent movement. Both the sinistral and convergent rates display a decreasing trend from west to east. The locking degree and slip deficit are higher in the western segment (with an average of about 0.74 and 1.1 mm/a) and lower in the eastern segment. Then, we construct a strain rate field using GPS observations to analyze the regional strain characteristics. The results indicate that along the fault, the western segment shows a larger shear strain rate and negative dilation rate. Regional earthquake records show that the frequency of earthquakes is lower near the fault. The joint results suggest that the western segment may have a higher earthquake risk. In addition, the insignificant fault slip rate in the eastern segment may indicate that it does not participate in the tectonic transformation among the Riyueshan, Lajishan, and West Qinling faults.

Gravitational Constraints on the Earth's Inner Core Differential Rotation

AbstractThe differential axial rotation of the solid inner core (IC) is suggested by seismic observations and expected from core dynamics models. A rotation of the IC by an angle α takes its degree 2, order 2 topography (peak‐to‐peak amplitude δh) out of its gravitational alignment with the mantle. This creates a gravity variation of degree 2, order 2 proportional to δh and to α. Here, we use gravity observations from Satellite Laser Ranging, the Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow‐On to reconstruct the time‐variable S2,2 Stokes coefficient. We show that for δh = 90 m, S2,2 provides upper bounds on α of 0.09°, 0.3°, and 0.4° at periods of ∼4, ∼6, and ∼12 years, respectively. These are overestimates, as our reconstructed S2,2 signal likely remains polluted by hydrology, although viscous relaxation of the IC can permit larger amplitudes.

Assimilating Seismic Data And Soil Basic Information to Predict Near-surface Permeability And Water Content

Soil permeability plays a key role in the study of earth sciences for agriculture and forestry. Because the traditional laboratory measurement is time-consuming and laborious, and only discrete profile data can be obtained. In this paper, a simple and efficient method using soil transfer function (PTFs), is proposed to estimate permeability and other parameters indirectly by seismic wave observation. Based on Biot’s theory establishing the relationship between shear-wave velocity and permeability, and the basic soil information in the study area, including land survey, logging and drilling data, using the PTFs constraint to estimate saturated hydraulic conductivity, soil components, saturated water content and residual water content may favor the physical corrections of the PTFs. The results show that there are three distinct layers of permeability in the study area, consistent with the velocity structure. The measured basic information further verifies the accuracy of the saturated permeability coefficient estimation, which shows the reliability of the improved PTFs. The established near surface permeability estimation method by integrating seismic wave observation and soil basic information may provide further assimilation for parameter estimation combining geophysical observation and hydrogeographic statistics.

Low‐Velocity Structure of Subducted Oceanic Crust in the Upper Mantle: Insights From High Pressure and Temperature Elasticity Measurements of Aragonite

AbstractThis study examines the effects of water and carbon on the velocity profiles of subducted oceanic crust in the upper mantle. High pressure and temperature Brillouin measurements were conducted to determine the single‐crystal elasticity of aragonite (CaCO3) up to 20 GPa and 600 K, respectively. Using the finite‐strain method, we determined the elastic parameters: KS0 = 70.7 GPa, G0 = 36.1(3) GPa with KS′ = 5.0(1), G′ = 1.3(1), ∂KS/∂T = −0.020(2) GPa/K, and ∂G/∂T = −0.015(1) GPa/K. When combined with literature results, our findings reveal that neither hydrous minerals nor carbonate alone can explain the observed 3%–4% and 3%–7% low velocity anomalies in the compressional (VP) and shear‐wave (VS) velocities of subducted oceanic crust in the circum‐Pacific region at 150–250 km depth. Considering the combined effect of water and carbonate, the addition of 5.9–6.9 vol.% aragonite together with 10 vol.% lawsonite in the oceanic crust can produce a 3%–4% and 6.1%–8.1% low velocity anomalies in the VP and VS at 150–250 km depth, respectively, consistent with the seismic observations in the region. Complete dehydration of lawsonite accompanied by the gradual decarbonization of the subducted oceanic crust explains the absence of low‐velocity anomalies below 300‐km depth. Our findings help for a better understanding of the possible influence of carbon and water on seismic velocities of the mantle. As a result, the circulation of carbon and water may be better understood while taking the complex velocity structure of subduction zones and mineral physics findings into account.

3D near-surface P-wave velocity structure imaging with Distributed Acoustic Sensing and electric hammer source

Distributed Acoustic Sensing (DAS) is an emerging technique for ultra-dense seismic observation, which provides a new method for high-resolution sub-surface seismic imaging. Recently a large number of linear DAS arrays have been used for two-dimensional S-wave near-surface imaging in urban areas. In order to explore the feasibility of three-dimensional (3D) structure imaging using a DAS array, we carried out an active source experiment at the Beijing National Earth Observatory. We deployed a 1 ​km optical cable in a rectangular shape, and the optical cable was recast into 250 sensors with a channel spacing of 4 ​m. The DAS array clearly recorded the P, S and surface waves generated by a hammer source. The first-arrival P wave travel times were first picked with a Short-Term Average/Long-Term Average (STA/LTA) method and further manually checked. The P-wave signals recorded by the DAS are consistent with those recorded by the horizontal components of short-period seismometers. At shorter source-receiver distances, the picked P-wave arrivals from the DAS recording are consistent with vertical component recordings of seismometers, but they clearly lag behind the latter at greater distances. This is likely due to a combination of the signal-to-noise ratio and the polarization of the incoming wave. Then, we used the TomoDD software to invert the 3D P-wave velocity structure for the uppermost 50 ​m with a resolution of 10 ​m. The inverted P-wave velocity structures agree well with the S-wave velocity structure previously obtained through ambient noise tomography. Our study indicates the feasibility of 3D near-surface imaging with the active source and DAS array. However, the inverted absolute velocity values at large depths may be biased due to potential time shifts between the DAS recording and seismometer at large source-receiver distances.

Thermal pattern of Tatun volcanic system by satellite-observed temperatures and its correlation with earthquake magnitudes

The land surface temperature (LST) of volcanoes detected from satellite sensors reflects the thermal status of heat sources in the subsurface. Volcanic earthquakes occur as magma and volcanic fluids transport to the surface from depth. Thus, both LST and earthquake magnitude are key parameters to the study of active volcanoes. Here we investigate the volcanic status of Tatun Volcanic Group (TVG) based on LST and seismic observations. The Earth-observing satellites onboard thermal sensor derived land surface temperature, and the seismic records retrieved volcanic earthquake magnitude are used to delineate the past and current pattern of volcanic activity plus the future trend of the TVG. The spatiotemporal distribution of LST and volcanic earthquake magnitude in TVG are analyzed. The high-similarity trends of the 4-decade LST time series and 3-decade earthquake magnitude time series are inspected. The retrieved surface thermal pattern shows the non-steady-state nature of the subsurface thermal sources at this volcanic complex. The LST trend exhibits a rather positive correlation with the energy released from volcanic earthquakes and consequently, the presumption on the connection between LSTs and earthquakes is validated.

Mantle Deformation Processes During the Rift‐To‐Drift Transition at Magma‐Poor Margins

AbstractThe rift‐to‐drift transition at rifted margins is an area of active investigation due to an incomplete understanding of the spatial and temporal evolution of physical and chemical processes at the ocean‐continent transition (OCT). Deep structures that characterize modern OCTs are often difficult to identify by seismic observations, while terrestrial exposures are preserved in fragments separated by tectonic discontinuities. Numerical modeling is a powerful method for contextualizing physical processes and observations relevant to rifted margin evolution. We synthesize results from geological observations of fossil OCTs preserved in ophiolites, a recent seismic experiment on the Ivorian margin, and numerical modeling to characterize mantle deformation and melt production for magma‐poor margins. Across varied surface heat fluxes, mantle potential temperatures, and extension rates, our model results show homologies with geological observations. We propose that the development of large shear zones in the mantle, melt infiltration, grain size reduction, and anastomosing detachment faults control the structure of OCTs. We also infer that a hot, upwelling, melt‐rich asthenosphere is an important control on the local stress environment. During the exhumation phase, continentward‐dipping shear zones couple with seaward‐dipping detachment faults to exhume the subcontinental and formerly asthenospheric mantle. The mantle forms core‐complex‐like domes of peridotite at or near the surface. The faults that exhume these peridotite bodies are largely anastomosing and exhibit magmatic accretion in their footwalls. A combination of magmatic accretion and volcanic activity derived from the shallow melt region constructs the oceanic lithosphere in the footwalls of the out‐of‐sequence continentward‐dipping detachment faults in the oceanic crust and subcontinental mantle.

Mapping Glacier Basal Sliding Applying Machine Learning

AbstractDuring the RESOLVE project (“High‐resolution imaging in subsurface geophysics: development of a multi‐instrument platform for interdisciplinary research”), continuous surface displacement and seismic array observations were obtained on Glacier d’Argentière in the French Alps for 35 days in May 2018. The data set is used to perform a detailed study of targeted processes within the highly dynamic cryospheric environment. In particular, the physical processes controlling glacial basal motion are poorly understood and remain challenging to observe directly. Especially in the Alpine region for temperate based glaciers where the ice rapidly responds to changing climatic conditions and thus, processes are strongly intermittent in time and heterogeneous in space. Spatially dense seismic and Global Positioning System (GPS) measurements are analyzed applying machine learning to gain insight into the processes controlling glacial motions of Glacier d’Argentière. Using multiple bandpass‐filtered copies of the continuous seismic waveforms, we compute energy‐based features, develop a matched field beamforming catalog and include meteorological observations. Features describing the data are analyzed with a gradient boosting decision tree model to directly estimate the GPS displacements from the seismic noise. We posit that features of the seismic noise provide direct access to the dominant parameters that drive displacement on the highly variable and unsteady surface of the glacier. The machine learning model infers daily fluctuations and longer term trends. The results show on‐ice displacement rates are strongly modulated by activity at the base of the glacier. The techniques presented provide a new approach to study glacial basal sliding and discover its full complexity.

Anomalous Sound Velocities of Talc at High Pressure and Implications for Estimating Water Content in Mantle Wedge

AbstractTalc is widely distributed in metamorphic and hydrothermally altered rocks, and its velocity is useful for constraining the velocity of mantle wedge. In this study, we used ultrasonic interferometry to study the elastic properties of talc polycrystalline samples at high pressure (up to 5.6 GPa). VP increases monotonically with increasing pressure, and VS increases with pressure up to ∼3.2 GPa, but negative pressure dependence occurs above ∼3.2 GPa. Talc has a VP/VS ratio of 1.63–2.14 and Poisson's ratio of 0.20–0.36 at the pressure of 1.2–5.6 GPa, which can be used to explain some anomalous seismic observations in subduction zone. Combining the elastic properties of talc with those of olivine, orthopyroxene, and antigorite, the wave velocity of mantle wedge was calculated as a function of water content. The results suggest that lower water content is needed to match the observed values of seismic tomography when talc is present, compared with the antigorite scenario.

Active Dipping Interface of the Southern San Andreas Fault Revealed by Space Geodetic and Seismic Imaging

AbstractThe Southern San Andreas Fault (SSAF) in California is one of the most thoroughly studied faults in the world, but its configuration at seismogenic depths remains enigmatic in the Coachella Valley. We use a combination of space geodetic and seismic observations to demonstrate that the relatively straight southernmost section of the SSAF, between Thousand Palms and Bombay Beach, is dipping to the northeast at 60–80° throughout the upper crust (<10 km), including the shallow aseismic layer. We constrain the fault attitude in the top 2–3 km using inversions of surface displacements associated with shallow creep, and seismic data from a dense nodal array crossing the fault trace near Thousand Palms. The data inversions show that the shallow dipping structure connects with clusters of seismicity at depth, indicating a continuous throughgoing fault surface. The dipping fault geometry has important implications for the long‐term fault slip rate, the intensity of ground shaking during future large earthquakes, and the effective strength of the southern SAF.

A Single‐Station Method for Seismic Detection of Slow Earthquakes: Applications to Japan and the Mexican Subduction Zone

AbstractSlow‐earthquake signals are generally smaller than or comparable to noise levels at almost all seismological frequencies. Comprehensive detection of these events requires continuous waveforms from many stations, but such data are not always available, even in regions with high slow‐earthquake activity. We therefore need a simple and stable detection method that is also applicable to regions with sparse seismic observation networks to truly advance our understanding of slow earthquakes. Here, we utilize the proportionality between the seismic energy rate and seismic moment rate of slow earthquakes to develop a slow‐earthquake detection method using broadband waveforms from a single station. We introduce the method and estimate its performance using continuous waveform data in Japan. The new method only detects events when the tectonic tremors occur near the station, for example, 89.1% of the detections have corresponding tectonic tremor activities within 60 km, which suggests that the false‐positive rate is low. We then apply it to waveform data from the Guerrero, Oaxaca, and Jalisco regions along the Mexican subduction zone. The results for the Guerrero and Oaxaca regions are largely consistent with the geodetic and seismological results from a previous study. We also detect many events in the Jalisco region using a permanent station and provide the first seismological report of long‐term slow‐earthquake activity in this region. Some of the large‐scale slow‐earthquake activity that is detected in this study is consistent with the timings of known slow‐slip events, which implies that other unknown slow‐slip events may have occurred when we detected many events.

Tracking changes in the co-eruptive seismic tremor associated with magma degassing at Piton de la Fournaise volcano

Volcanic eruptions in basaltic systems represent the most frequent expression of volcanic activity in the World, whose surface manifestations range from effusive lava flows to more explosive events. Their eruptive dynamics and style are largely controlled by the degassing of magma. In this study, we analyze seismic tremors recorded during 23 eruptions that occurred at Piton de la Fournaise volcano (La Réunion, France) during 2014–2021 and show that their properties are strongly linked to the magma degassing. We apply the network covariance matrix method to build a catalog of tremor sources associated with all 23 eruptions, to measure their time-frequency properties, and to locate their sources. We then conduct a multi-disciplinary analysis by comparing the seismic tremor observations with images of the eruptive sites, lava flow rate measurements, impulsive earthquakes, as well as deformation and magma composition data. The tremor depth distribution is found to be well correlated with magma gas content, indicating that the tremor generating mechanism is largely controlled by the shallow magma degassing under the eruptive site. The resulting seismic tremor signals analysis enable to track changes in degassing regimes and associated eruptive styles at the surface, thus improving monitoring and knowledge of the eruptive dynamics of the studied volcano.

Potential Seismic Hazard in Seoul, South Korea: A Comprehensive Analysis of Geology, Seismic, and Geophysical Field Observations, Historical Earthquakes, and Strong Ground Motions

ABSTRACT A series of moderate-size (Mw 4.0–6.0) earthquakes occurred in South Korea after the 2011 Mw 9.0 Tohoku–Oki megathrust earthquake, incurring public concern about possible occurrence of devastating earthquakes in Seoul—the capital city of South Korea, where historical seismic damage was reported. The seismicity is distributed in Seoul, being dominated by strike-slip earthquakes. The fault planes are oriented in north-northeast–south-southwest, which is a favorable direction to respond to the ambient stress field. Higher rates of seismicity are observed in the northwestern Seoul at depths of <10 km. Micro-to-small earthquakes occur episodically in the central Seoul along the Chugaryeong fault system that traverses Seoul in north–south. Seismic, geophysical, and geological properties illuminate the fault structures. Stochastic modeling of ground motions reproduces the seismic damages of historical earthquakes reasonably, supporting the occurrence of devastating historical earthquakes in Seoul. The seismicity distribution, focal mechanism solutions, geological features, and seismic and geophysical properties suggest the possible presence of earthquake-spawning blind faults in Seoul. The peak ground motions are assessed for moderate-size scenario earthquakes (Mw 5.4 with focal depth of 7 km) at six representative subregions in Seoul. The upper bounds of peak ground accelerations reach ∼11 m/s2. The seismic damage potentials for moderate-size earthquakes are high in most areas of Seoul, particularly around river sides covered by alluvium.

Active subduction and strain partitioning in western Myanmar revealed by a dense survey GNSS network

The active tectonics of the Myanmar (formerly Burma) microplate in central Myanmar are caused by the oblique collision between the Sunda and Eurasia plates to the east and the Indian plate to the west. The kinematics of this plate boundary zone remain a subject of controversy, due in part to limited geodetic and seismic observations in the region. Here, we report 41 new and 67 updated survey-mode and continuous Global Navigation Satellite System (GNSS) velocities spanning central and western Myanmar, between 15 – 23°N and 93 – 98°E. For the first time, our network allows us to directly observe the occurrence of active convergence along the central Rakhine coast in a manner consistent with active subduction and kinematic coupling of the megathrust, at a rate of 15 – 20 mm/yr. The data also establish a nearly constant strike-slip rate of 20 – 23 mm/yr along the entire length of the Sagaing fault. Using an elastic block-modeling approach, we examine several proposed models for the partitioning of oblique convergence into strike-slip and thrust motion across the Indo-Myanmar Ranges and the Myanmar Central Basin. We find evidence for shortening perpendicular to local topography throughout the Indo-Myanmar Ranges, indicating active convergence and nearly complete strain partitioning between thrust and transform faults north of 19°N, and estimate a rate of moment-deficit accumulation along the megathrust of 0.8−1.5×1019 N m/yr, equivalent to one MW 9.0 earthquake every 250 to 450 years, assuming all elastic strain is released seismically. We also find that a relatively small set of major active crustal faults can explain much of the data, though the data cannot resolve smaller faults that may be active at rates lower than a few mm/yr. Our results have important implications for earthquake and tsunami preparedness across the region, and indicate areas for future study where the geodetic observations remain ambiguous about ongoing fault activity and seismic hazard, particularly within the densely populated Ayeyarwady Delta region and along the southern Rakhine coast.

Collective motion in hcp-Fe at Earth’s inner core conditions

Earth's inner core is predominantly composed of solid iron (Fe) and displays intriguing properties such as strong shear softening and an ultrahigh Poisson's ratio. Insofar, physical mechanisms to explain these features coherently remain highly debated. Here, we have studied longitudinal and shear wave velocities of hcp-Fe (hexagonal close-packed iron) at relevant pressure-temperature conditions of the inner core using in situ shock experiments and machine learning molecular dynamics (MLMD) simulations. Our results demonstrate that the shear wave velocity of hcp-Fe along the Hugoniot in the premelting condition, defined as T/Tm (Tm: melting temperature of iron) above 0.96, is significantly reduced by ~30%, while Poisson's ratio jumps to approximately 0.44. MLMD simulations at 230 to 330 GPa indicate that collective motion with fast diffusive atomic migration occurs in premelting hcp-Fe primarily along [100] or [010] crystallographic direction, contributing to its elastic softening and enhanced Poisson's ratio. Our study reveals that hcp-Fe atoms can diffusively migrate to neighboring positions, forming open-loop and close-loop clusters in the inner core conditions. Hcp-Fe with collective motion at the inner core conditions is thus not an ideal solid previously believed. The premelting hcp-Fe with collective motion behaves like an extremely soft solid with an ultralow shear modulus and an ultrahigh Poisson's ratio that are consistent with seismic observations of the region. Our findings indicate that premelting hcp-Fe with fast diffusive motion represents the underlying physical mechanism to help explain the unique seismic and geodynamic features of the inner core.

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

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

Rapid Source Characterization of the 2023 Mw 6.8 Al Haouz, Morocco, Earthquake

Abstract The U.S. Geological Survey (USGS) National Earthquake Information Center (NEIC) estimates source characteristics of significant damaging earthquakes, aiming to place events within their seismotectonic framework. Contextualizing the 8 September 2023, Mw 6.8 Al Haouz, Morocco, earthquake is challenging, because it occurred in an enigmatic region of active surface faulting, and low seismicity yet produced significant damage and loss of life. Here, we present the rapid earthquake source products produced by the USGS NEIC, describing how the source model was derived using both seismic and geodetic observations. Our analysis indicates that the earthquake was the result of oblique-reverse faulting in the lower crust on either a steeply north-dipping fault or a moderately south-dipping fault. Finite-slip models using seismic and geodetic data reveal a compact source, with slip occurring at depths of 15–35 km. The causative fault is not apparent, because the rupture did not break the surface, and it is not possible to definitively attribute the earthquake to a known structure. The earthquake centroid depth of 25 km is noteworthy, because it shows slip extending beyond common estimates of seismogenic depth. This earthquake highlights that the seismogenic processes associated with mountain building in this wide plate boundary region are poorly understood.