Seismic source analysis of the destructive earthquake November 21, 2022, Mw 5.6 Cianjur (Indonesia) from relocated aftershock

A destructive shallow earthquake with a magnitude of 5.6 struck Cianjur, West Java, Indonesia on November 21, 2022. This earthquake resulted in 602 casualties and the collapse of over 67,504 residences. The day after the mainshock, we deployed 19 temporary seismic stations to monitor aftershocks for a period of 30 days. We manually picked arrival times for 4499 P-waves and 3419 S-waves and determined locations for 514 events. Following the velocity model update, phase refinement through waveform cross correlation, and relocation using double-difference methods, we were able to determine 442 well-defined hypocenters of the aftershocks. We identified two clusters of aftershocks: one in the NNW-SSE direction, with a length of about 8 km, and another in the WSW-ENE direction, with a length of around 6 km. The seismogenic zone of these clusters ranges from a depth of 3 to 13 km. Our interpretation suggests that these clusters may indicate a conjugate fault. It is possible that the mainshock (Mw5.6) Cianjur earthquake on November 21, 2022 occurred on the WSW-ENE direction with sinistral movement.

Cluster Analysis of Seismicity in the Eastern Gulf of Corinth Based on a Waveform Template Matching Catalog

The Corinth Rift, in Central Greece, is one of the most seismically active areas in Europe. In the eastern part of the Gulf of Corinth, which has been the site of numerous large and destructive earthquakes in both historic and modern times, a pronounced earthquake swarm occurred in 2020-2021 at the Perachora peninsula. Herein, we present an in-depth analysis of this sequence, employing a high-resolution relocated earthquake catalog, further enhanced by the application of a multi-channel template matching technique, producing additional detections of over 7600 events between January 2020 and June 2021. Single-station template matching enriches the original catalog thirty-fold, providing origin times and magnitudes for over 24,000 events. We explore the variable levels of spatial and temporal resolution in the catalogs of different completeness magnitudes and also of variable location uncertainties. We characterize the frequency-magnitude distributions using the Gutenberg-Richter scaling relation and discuss possible b-value temporal variations that appear during the swarm and their implications for the stress levels in the area. The evolution of the swarm is further analyzed through spatiotemporal clustering methods, while the temporal properties of multiplet families indicate that short-lived seismic bursts, associated with the swarm, dominate the catalogs. Multiplet families present clustering effects at all time scales, suggesting triggering by aseismic factors, such as fluid diffusion, rather than constant stress loading, in accordance with the spatiotemporal migration patterns of seismicity.

Tectonic Implication of the 2022 MS 6.9 Earthquake in Menyuan, Qinghai, China: Analysis of Precise Earthquake Locations and InSAR

Precise earthquake locations and InSAR (Interferometric Synthetic Aperture Radar) deformation observation are the major methods to understand the earthquake occurrence and disaster-causing process. This paper proposes a processing framework for analyzing strong earthquake mechanisms from one-dimensional velocity inversion to precise earthquake locations combined with InSAR deformation observation, and discusses earthquake-generating fault and dynamic mechanisms of tectonic deformation. We analyzed the Menyuan Ms 6.9 earthquake in 2022 and discuss the historical seismic activities and corresponding stress adjustment processes in the research region. To analyze and study the seismogenic structure and mechanism of the earthquake, we investigated the spatial and temporal distribution characteristics of the Menyuan earthquake sequence and analyzed the InSAR coseismic deformation field. We obtained the precise locations of the main shock and aftershocks and the coseismic InSAR deformation field of the main shock. It was confirmed that the Ms 6.9 earthquake was a shallow sinistral strike-slip earthquake, which led to the sequential activation of the Tuolaishan and Lenglongling faults. The main seismogenic fault of the mainshock was the northwestern end of the Lenglongling fault, and the earthquake rupture was segmented. It can be inferred that the earthquake was a stress-adjusted event triggered in the Qilian-Haiyuan tectonic belt caused by the northeasterly push of the Qinghai-Tibet Plateau. The risk of moderate to high earthquakes in the region remains high in the future, requiring enhanced seismic observations.

Microseismicity and lithosphere thickness at a nearly-amagmatic oceanic detachment fault system

Oceanic detachment faults play a central role in accommodating the plate divergence at slow-ultraslow spreading mid-ocean ridges. Successive flip-flop detachment faults in a nearly-amagmatic region of the ultraslow spreading Southwest Indian Ridge (SWIR) at 64°30'E accommodate ~100% of plate divergence, with mostly ultramafic smooth seafloor. Here we present microseismicity data, recorded by ocean bottom seismometers, showing that the axial brittle lithosphere is on the order of 15 km thick under the nearly-amagmatic smooth seafloor, which is no thicker than under nearby volcanic seafloor or at more magmatic SWIR detachment systems. Our data reveal that microearthquakes with normal focal mechanisms are colocated with seismically-imaged damage zones of the active detachment fault and of antithetic hanging-wall faults. The level of the hanging-wall seismicity is significantly higher than that documented at more magmatic detachments of slow-ultraslow ridges, which may be a unique feature of nearly-amagmatic flip-flop detachment systems.

A comprehensive suite of earthquake catalogues for the 2016-2017 Central Italy seismic sequence

The protracted nature of the 2016-2017 central Italy seismic sequence, with multiple damaging earthquakes spaced over months, presented serious challenges for the duty seismologists and emergency managers as they assimilated the growing sequence to advise the local population. Uncertainty concerning where and when it was safe to occupy vulnerable structures highlighted the need for timely delivery of scientifically based understanding of the evolving hazard and risk. Seismic hazard assessment during complex sequences depends critically on up-to-date earthquake catalogues—i.e., data on locations, magnitudes, and activity of earthquakes—to characterize the ongoing seismicity and fuel earthquake forecasting models. Here we document six earthquake catalogues of this sequence that were developed using a variety of methods. The catalogues possess different levels of resolution and completeness resulting from progressive enhancements in the data availability, detection sensitivity, and hypocentral location accuracy. The catalogues range from real-time to advanced machine-learning procedures and highlight both the promises as well as the challenges of implementing advanced workflows in an operational environment.

Monitoring microseismicity of the Hengill Geothermal Field in Iceland

Induced seismicity is one of the main factors that reduces societal acceptance of deep geothermal energy exploitation activities, and felt earthquakes are the main reason for closure of geothermal projects. Implementing innovative tools for real-time monitoring and forecasting of induced seismicity was one of the aims of the recently completed COSEISMIQ project. Within this project, a temporary seismic network was deployed in the Hengill geothermal region in Iceland, the location of the nation’s two largest geothermal power plants. In this paper, we release raw continuous seismic waveforms and seismicity catalogues collected and prepared during this project. This dataset is particularly valuable since a very dense network was deployed in a seismically active region where thousand of earthquakes occur every year. For this reason, the collected dataset can be used across a broad range of research topics in seismology ranging from the development and testing of new data analysis methods to induced seismicity and seismotectonics studies.

Measurement(s)	Seismic waveforms (seismograms) • Seismicity (Origin time, location and magnitude of earthquakes)
Technology Type(s)	Seismic stations (velocity sensors) • SeisComP data acquisition and processing system

Microseismicity Appears to Outline Highly Coupled Regions on the Central Chile Megathrust

We compiled a novel microseismicity catalog for the Central Chile megathrust (29°-35°S), comprising 8,750 earthquakes between April 2014 and December 2018. These events describe a pattern of three trenchward open half-ellipses, consisting of a continuous, coast-parallel seismicity band at 30-45 km depth, and narrow elongated seismicity clusters that protrude to the shallow megathrust and separate largely aseismic regions along strike. To test whether these shapes could outline highly coupled regions ("asperities") on the megathrust, we invert GPS displacement data for interplate locking. The best-fit locking model does not show good correspondence to seismicity, possibly due to lacking resolution. When we prescribe high locking inside the half-ellipses, however, we obtain models with similar data fits that are preferred according to the Bayesian Information Criterion (BIC). We thus propose that seismicity on the Central Chile megathrust may outline three adjacent highly coupled regions, two of them located between the rupture areas of the 2010 Maule and the 2015 Illapel earthquakes, a segment of the Chilean margin that may be in a late interseismic stage of the seismic cycle.

New seismological data from the Calabrian arc reveal arc-orthogonal extension across the subduction zone

The Calabrian Arc subduction-rollback system along the convergent Africa/Eurasia plate boundary is among the most active geological structures in the Mediterranean Sea. However, its seismogenic behaviour is largely unknown, mostly due to the lack of seismological observations. We studied low-to-moderate magnitude earthquakes recorded by the seismic network onshore, integrated by data from a seafloor observatory (NEMO-SN1), to compute a lithospheric velocity model for the western Ionian Sea, and relocate seismic events along major tectonic structures. Spatial changes in the depth distribution of earthquakes highlight a major lithospheric boundary constituted by the Ionian Fault, which separates two sectors where thickness of the seismogenic layer varies over 40 km. This regional tectonic boundary represents the eastern limit of a domain characterized by thinner lithosphere, arc-orthogonal extension, and transtensional tectonic deformation. Occurrence of a few thrust-type earthquakes in the accretionary wedge may suggest a locked subduction interface in a complex tectonic setting, which involves the interplay between arc-orthogonal extension and plate convergence. We finally note that distribution of earthquakes and associated extensional deformation in the Messina Straits region could be explained by right-lateral displacement along the Ionian Fault. This observation could shed new light on proposed mechanisms for the 1908 Messina earthquake.

A Method to Improve the Distribution of Observations in GNSS Water Vapor Tomography

Water vapor is an important driving factor in the related weather processes in the troposphere, and its temporal-spatial distribution and change are crucial to the formation of cloud and rainfall. Global Navigation Satellite System (GNSS) water vapor tomography, which can reconstruct the water vapor distribution using GNSS observation data, plays an increasingly important role in GNSS meteorology. In this paper, a method to improve the distribution of observations in GNSS water vapor tomography is proposed to overcome the problem of the relatively concentrated distribution of observations, enable satellite signal rays to penetrate more tomographic voxels, and improve the issue of overabundance of zero elements in a tomographic matrix. Numerical results indicate that the accuracy of the water vapor tomography is improved by the proposed method when the slant water vapor calculated by GAMIT is used as a reference. Comparative results of precipitable water vapor (PWV) and water vapor density (WVD) profiles from radiosonde station data indicate that the proposed method is superior to the conventional method in terms of the mean absolute error (MAE), standard deviations (STD), and root-mean-square error (RMS). Further discussion shows that the ill-condition of tomographic equation and the richness of data in the tomographic model need to be discussed separately.

Dehydration of subducting slow-spread oceanic lithosphere in the Lesser Antilles

Subducting slabs carry water into the mantle and are a major gateway in the global geochemical water cycle. Fluid transport and release can be constrained with seismological data. Here we use joint active-source/local-earthquake seismic tomography to derive unprecedented constraints on multi-stage fluid release from subducting slow-spread oceanic lithosphere. We image the low P-wave velocity crustal layer on the slab top and show that it disappears beneath 60–100 km depth, marking the depth of dehydration metamorphism and eclogitization. Clustering of seismicity at 120–160 km depth suggests that the slab’s mantle dehydrates beneath the volcanic arc, and may be the main source of fluids triggering arc magma generation. Lateral variations in seismic properties on the slab surface suggest that serpentinized peridotite exhumed in tectonized slow-spread crust near fracture zones may increase water transport to sub-arc depths. This results in heterogeneous water release and directly impacts earthquakes generation and mantle wedge dynamics.

During subduction water is transported into the mantle, but constraining its release remains challenging. Here, using seismic tomography of the Lesser Antilles arc, the authors track the multistage dehydration of the slab and its lateral variations associated with heterogeneous slab composition.

Induced earthquakes. Sharp increase in central Oklahoma seismicity since 2008 induced by massive wastewater injection

Unconventional oil and gas production provides a rapidly growing energy source; however, high-production states in the United States, such as Oklahoma, face sharply rising numbers of earthquakes. Subsurface pressure data required to unequivocally link earthquakes to wastewater injection are rarely accessible. Here we use seismicity and hydrogeological models to show that fluid migration from high-rate disposal wells in Oklahoma is potentially responsible for the largest swarm. Earthquake hypocenters occur within disposal formations and upper basement, between 2- and 5-kilometer depth. The modeled fluid pressure perturbation propagates throughout the same depth range and tracks earthquakes to distances of 35 kilometers, with a triggering threshold of ~0.07 megapascals. Although thousands of disposal wells operate aseismically, four of the highest-rate wells are capable of inducing 20% of 2008 to 2013 central U.S. seismicity.