The Valparaiso 2017 sequence occurred in the Central Chile megathrust, an active zone where the last mega‐earthquake occurred in 1730. Intense seismicity started 2 days before the Mw 6.9 mainshock, a ...slow trenchward movement was observed in the coastal GPS antennas and was accompanied by foreshocks and repeater‐type seismicity. To characterize the rupture process of the mainshock, we perform a dynamic inversion using the strong‐motion records and an elliptical patch approach. We suggest that a slow slip event preceded and triggered the Mw 6.9 earthquake, which ruptured an elliptical asperity (semiaxis of 10 km and 5 km, with a subshear rupture, stress drop of 11.71 MPa, yield stress of 17.21 MPa, slip weakening of 0.65 m, and kappa value of 1.98). This earthquake could be the beginning of a long‐term nucleation phase to a major rupture, within the highly coupled Central Chile zone where a megathrust earthquake like 1730 is expected.
Key Points
Valparaiso 2017 Mw 6.9 was preceded by repeaters and a slow slip event
An intense seismicity was recorded before the mainshock
The mainshock broke a small asperity of few kilometers
Low-frequency earthquakes are a particular class of slow earthquakes that provide a unique source of information on the physical processes along a subduction zone during the preparation of large ...earthquakes. Despite increasing detection of these events in recent years, their source mechanisms are still poorly characterised, and the relation between their magnitude and size remains controversial. Here, we present the source characterisation of more than 10,000 low-frequency earthquakes that occurred during tremor sequences in 2012-2016 along the Nankai subduction zone in western Shikoku, Japan. We show that the scaling of seismic moment versus corner frequency for these events is compatible with an inverse of the cube law, as widely observed for regular earthquakes. Their radiation, however, appears depleted in high-frequency content when compared to regular earthquakes. The displacement spectrum decays beyond the corner frequency with an omega-cube power law. Our result is consistent with shear rupture as the source mechanism for low-frequency earthquakes, and suggests a self-similar rupture process and constant stress drop. When investigating the dependence of the stress drop value on the rupture speed, we found that low-frequency earthquakes might propagate at lower rupture velocity than regular earthquakes, releasing smaller stress drop.
A major aim in the study of crustal fluids is the development of automatic methodologies for monitoring deep‐source, non‐volcanic gas emissions’ spatio‐temporal evolution. Crustal fluids play a ...significant role in the generation of large earthquakes and the characterization of their emissions on the surface can be essential for better understanding crustal processes generating earthquakes. We investigate seismic tremors recorded over 4 days in 2019 at the Mefite d’Ansanto (southern Apennines, Italy) that is located at the northern end of the fault system that generated the Mw 6.9 1980 Irpinia Earthquake. The Mefite d’Ansanto is hypothesized to be the largest natural, non‐volcanic, CO2‐rich gas emission on Earth. The seismic tremor is studied by employing a dense temporary seismic network and an automated detection algorithm based on non‐parametric statistics of the recorded signal amplitudes. We extracted signal characteristics (RMS amplitude and statistical moments of amplitudes both in time and frequency domains) for use in the subsequent supervised machine‐learning classification of the target tremor and accidently detected anthropogenic and background noise. The data set is used for the training and optimization of station‐based KNN (k‐Nearest‐Neighbors) binary classifiers obtaining good classification performances with a median overall accuracy across all stations of 92.8%. The classified tremor displayed common features at all stations: variable duration (16 s to 30–40 min), broad peak frequency (4–20 Hz) with varying amplitudes, and two types of signals: (a) long‐duration, high‐amplitude tremor and (b) pulsating tremor. Higher tremor amplitudes recorded at stations closer to local bubbling and pressurized vents suggest multiple local tremor sources.
Key Points
A dense temporary seismic array enabled the acquisition of emission‐related seismic tremor at the Mefite d’Ansanto site
A detection algorithm and supervised machine‐learning approach for the identification and classification of tremors has been developed
New insights on the emission tremor would help to monitor the CO2 degassing at the Mefite d’Ansanto and its temporal evolution
The accurate characterization of microearthquake sequences allows seismologists to better understand the physical processes involved in earthquake nucleation and rupture propagation and to gain ...insights on fault geometry at depth. Standard procedures for seismic sequences analysis are based on manual detection and phase-picking, requiring a huge amount of work from expert seismologists, particularly in the case of microseismic events. Here we show how the investigation of a low-magnitude foreshock-mainshock-aftershock sequence, occurred in August 2020 close to Castelsaraceno village (southern Italy), greatly benefited from the application of a semi-automated template matching and machine-learning based workflow. The phase-picking was automatically performed through a deep-learning algorithm on 202 microearthquakes detected between July and October 2020, followed by an automatic multi-step absolute and relative earthquake location procedure. The 72 relocated events of the seismic sequence were clustered in time (7-12 August) and in a narrow range of depths (10-12 km). The Ml 2.1 foreshock doublet and the Ml 2.9 mainshock identified a persistent asperity. The joint analysis of aftershocks distribution, the mainshock focal mechanism, and the geology of the study area suggest the occurrence of the sequence along a NNE-SSW left-lateral, transtensional fault in the brittle portion of the crystalline basement.
SUMMARY
The retrieval of earthquake finite-fault kinematic parameters after the occurrence of an earthquake is a crucial task in observational seismology. Routinely used source inversion techniques ...are challenged by limited data coverage and computational effort, and are subject to a variety of assumptions and constraints that restrict the range of possible solutions. Back-projection (BP) imaging techniques do not need prior knowledge of the rupture extent and propagation, and can track the high-frequency (HF) radiation emitted during the rupture process. While classic source inversion methods work at lower frequencies and return an image of the slip over the fault, the BP method highlights fault areas radiating HF seismic energy. Patterns in the HF radiation are attributable to the spatial and temporal complexity of the rupture process (e.g. slip heterogeneities, changes in rupture speed). However, the quantitative link between the BP image of an earthquake and its rupture kinematics remains unclear. Our work aims at reducing the gap between the theoretical studies on the generation of HF radiation due to earthquake complexity and the observation of HF emissions in BP images. To do so, we proceed in two stages, in each case analysing synthetic rupture scenarios where the rupture process is fully known. We first investigate the influence that spatial heterogeneities in slip and rupture velocity have on the rupture process and its radiated wave field using the BP technique. We simulate two different rupture processes using a 1-D line source model: a homogeneous process, where the kinematic parameters are constant along the line, and a heterogeneous process, where we introduce a central segment along the line that has a step change in kinematics. For each rupture model, we calculate synthetic seismograms at three teleseismic arrays and apply the BP technique to reveal how HF emissions are influenced by the three kinematic parameters controlling the synthetic model: the rise time, final slip and rupture velocity. Our results show that the HF peaks retrieved from BP analysis are better associated with space–time heterogeneities of slip acceleration. We then build on these findings by testing whether one can retrieve the kinematic rupture parameters along the fault using information from the BP image alone. We apply a machine learning, convolutional neural network (CNN) approach to the BP images of a large set of simulated 1-D rupture processes to assess the ability of the network to retrieve, from the progression of HF emissions in space and time, the kinematic parameters of the rupture. These rupture simulations include along-strike heterogeneities whose size is variable and within which the parameters of rise-time, final slip and rupture velocity change from the surrounding rupture. We show that the CNN trained on 40 000 pairs of BP images and kinematic parameters returns excellent predictions of the rise time and the rupture velocity along the fault, as well as good predictions of the central location and length of the heterogeneous segment. Our results also show that the network is insensitive towards the final slip value, as expected from theoretical results.
We analyze slip distribution and rupture kinematics of a Mw3.3 induced event that occurred in the St. Gallen geothermal reservoir (NE Switzerland) in 2013. We carry out a two‐step procedure: (1) path ...effects are deconvolved from the seismograms using an empirical Green's function, resulting in relative source time functions at all seismic stations; (2) the relative source time functions are back‐projected to the corresponding isochrones on the fault plane. Results reveal that the mainshock rupture propagates toward NNE from the hypocenter with an average velocity of 2,000 m/s. Spatiotemporal organization of foreshocks and aftershocks shows that the mainshock broke a previously less active portion of the fault and suggests that the aftershock sequence could be mainly driven by stress transfer. Applying this method in an operational environment could enable fast retrieval of seismic slip, allowing assessment of fault asperities and structures involved in the reservoir creation process.
Plain Language Summary
In most models and analyses, small earthquakes (i.e., magnitude less than 4) are considered either point sources or homogeneous “penny‐shaped” surfaces. While these assumptions may be valid, details of earthquake ruptures are more complex. Here we study a magnitude 3 induced earthquake that occurred in the St. Gallen geothermal reservoir (NE Switzerland) in 2013. We image the earthquake rupture by refocusing on the fault the seismic energy recorded at six sensors located within 15 km from the source. Our results show that a detailed description of the rupture process of such a small earthquake can indeed be obtained: The rupture propagates from the hypocenter in NNE direction for 150 m, with an average velocity of 2 km/s, breaking into a less active portion of the fault, where no earthquake was previously recorded. The proposed method could be routinely applied during geothermal reservoir operations to allow rapid assessment of fault structures involved in the reservoir creation process.
Key Points
Our approach reveals rupture complexity for a magnitude 3.3 fluid‐induced earthquake at 100‐m scale
The rupture propagated toward NNE, breaking into a previously less active portion of the fault
Early aftershocks distribute around the main slip area, suggesting a stress transfer mechanism
SUMMARY
The M= 8.1, 1 April 2014 Iquique earthquake, which broke part of the northern Chile seismic gap, was preceded by a strong foreshock sequence starting early January 2014. The reported analysis ...of the continuous records of the nearby GPS stations from the Integrated Plate Observatory Chile, North Chile array lead to contradictory results concerning the existence and location of slow slip events (SSEs) on the interplate contact. Resolving this controversy is an important issue, as although many SSEs are reported in subduction zones, only a few were found to be precursory to large earthquakes. Here we show that the records of a long base tiltmeter installed near Iquique, when corrected for coseismic steps, long-term drift, tidal signals and oceanic and atmospheric loading, show significant residual signals. These can be modelled with a sequence of four SSEs located close to Iquique. Their signature was already reported on some GPS stations, but their source was then characterized with a very low resolution in time and space, leading to contradicting models. With the tilt records, we can rule out the previously proposed models with a single large SSE closer to the main shock. Combining tilt with GPS records greatly improves the resolution of GPS alone, and one could locate their sources 100–180 km south–southeast to the main shock epicentre, with moment magnitudes between 5.8 and 6.2, at the edge of the main aftershock asperities. These moderate SSEs thus did not directly trigger the main shock, but contributed to trigger the main foreshock and the main aftershock. Only the sensitivity and resolution of the tiltmeter, added to the GPS records, allowed us to describe with unprecedented accuracy this precursory process as a cascade of cross-triggered, short-term aseismic slip events and earthquakes on the interplate contact. This three months of precursory activation appears to be the final acceleration burst of a weaker, longer term SSE which started mid-2013, already reported, with a moment release history which we could quantify. From the methodological point of view, our study takes advantage of an interesting complementarity of tilt and GPS measurements, due to their different dependence in distance to the source of strain, which turns out to be very efficient for resolving location and moment of strain sources, even when both instruments are close to each other. It finally demonstrates the efficient removal of sequences of small or even undetected coseismic steps from high resolution tilt record signal in order to retrieve the purely aseismic signal, a presently impossible task for high time resolution GPS records due to low signal to noise.
We investigate the effect of extended faulting processes and heterogeneous wave propagation on the early warning system capability to predict the peak ground velocity (PGV) from moderate to large ...earthquakes occurring in the southern Apennines (Italy). Simulated time histories at the early warning network have been used to retrieve early estimates of source parameters and to predict the PGV, following an evolutionary, probabilistic approach. The system performance is measured through the Effective Lead‐Time (ELT), i.e., the time interval between the arrival of the first S‐wave and the time at which the probability to observe the true PGV value within one standard deviation becomes stationary, and the Probability of Prediction Error (PPE), which provides a measure of PGV prediction error. The regional maps of ELT and PPE show a significant variability around the fault up to large distances, thus indicating that the system's capability to accurately predict the observed peak ground motion strongly depends on distance and azimuth from the fault.
A versatile synthetic route based on magnetic Fe
3
O
4
nanoparticle (MNP) prefunctionalization with a phosphonic acid monolayer has been used to covalently bind the gH625 peptide on the nanoparticle ...surface. gH625 is a membranotropic peptide capable of easily crossing the membranes of various cells including the typical human blood-brain barrier components. A similar synthetic route was used to prepare another class of MNPs having a functional coating based on PEG, rhodamine, and folic acid, a well-known target molecule, to compare the performance of the two cell-penetrating systems (i.e., gH625 and folic acid). Our results demonstrate that the uptake of gH625-decorated MNPs in immortalized human brain microvascular endothelial cells after 24 h is more evident compared to folic acid-functionalized MNPs as evidenced by confocal laser scanning microscopy. On the other hand, both functionalized systems proved capable of being internalized in a brain tumor cell line (i.e., glioblastoma A-172). These findings indicate that the functionalization of MNPs with gH625 improves their endothelial cell internalization, suggesting a viable strategy in designing functional nanostructures capable of first crossing the BBB and, then, of reaching specific tumor brain cells.
The analysis of similar earthquakes, such as events in a seismic sequence, is an effective tool with which to monitor and study source processes and to understand the mechanical and dynamic states of ...active fault systems. We are observing seismicity that is primarily concentrated in very limited regions along the 1980 Irpinia earthquake fault zone in Southern Italy, which is a complex system characterised by extensional stress regime. These zones of weakness produce repeated earthquakes and swarm-like microearthquake sequences, which are concentrated in a few specific zones of the fault system. In this study, we focused on a sequence that occurred along the main fault segment of the 1980 Irpinia earthquake to understand its characteristics and its relation to the loading-unloading mechanisms of the fault system.