Seismic Source Processes of 25 Earthquakes (Mw>5) in the Gulf of California

ABSTRACT The Gulf of California (GoC) is a complex tectonic boundary that has been instrumented in the past several decades to record broadband seismograms. This volume of data has allowed us to study several source parameters systematically. Before, only a few source parameters of earthquakes greater than magnitude five had been studied in the GoC area. We re-examined the focal mechanisms of several earthquakes in the southern GoC that occurred over the last 20 yr using local–regional distance broadband seismograms. These focal mechanisms were then used as input data to retrieve the time–space history of the rupture for each earthquake. This work contributes to the study of 25 rupture-process models computed with the method proposed by Yagi et al. (1999). To investigate more about the nature of the seismicity in the GoC, we also calculated the non-double-couple component of moment tensors for 45 earthquakes. Previous studies (e.g., Ortega et al., 2013, 2016) have shown that non-double-couple components from moment tensors in this region are associated with complex faulting, suggesting that oblique faults or several parallel faults are interacting simultaneously. Our results show that, at least for moderate earthquakes (5 < M < 6), rupture processes in the GoC show a complex interaction between fault systems. It is revealed on the important contribution of non-double-couple component obtained in the full moment tensor analysis.

SEISMICITY of KAMCHATKA and COMMANDER ISLANDS in 2015

The seismicity review of Kamchatka and surrounding territories for 2015 is given. In the Kamchatka earthquake catalogue, the minimum local magnitude of completeness is MLmin=3.5, and for earthquakes with h≥350 km under the Okhotsk sea MLmin=3.6. The Kamchatka earthquake catalogue for 2015 with ML3.5, published in the Appendix to this issue, includes 1213 events. 92 earthquakes of the catalogue with ML=3.0–6.5 were felt in Kamchatka and surrounding areas with seismic intensity I=2–6 according to the MSK-64 scale. For all events with ML5.0 that occurred in 2015 in the KB GS RAS area of responsibility, an attempt to calculate the seismic moment tensor (SMT) was made. There are 32 such events in the regional catalogue. For 28 earthquakes the SMT and depth h of the equivalent point source were calculated successfully. The calculations were performed for the SMT double-couple model using a nonlinear algorithm. In 2015, a typical location of the earthquake epicenters was observed in the Kamchatka zone. In 2015, the seismicity level in all selected zones and in the region as a whole correspond to the background one according to the “SESL’09” scale. The number of recorded events with ML3.5 and strong earthquakes with ML5.0 is close to the average annual value. Anomalous and significant events were not recorded.

2000 Uglegorsk earthquake

The paper presents the results of waveform inversion of the Mw 6.8 August 4 (5), 2000 Uglegorsk earthquake (Sakhalin Island, Russia). The detailed rupture process of the 2000 Uglegorsk earthquake was simulated using the waveform inversion method. The average parameters were calculated for both nodal planes. Waveform inversion was carried out on the basis of Global Seismographic Network (GSN) data. Only P-waves from BHZ channels of all stations from the GSN were used. The simulated source parameters included a double-couple source, the scalar seismic moment, the source time function, and the slip directions. The performed studies made it possible to investigate the features of the rupture development and the amplitude of displacements along the east and west-dipping nodal planes of the August 4 (5), 2000 Uglegorsk earthquake. The obtained P-slip model for the 2000 Uglegorsk earthquake source area is in good agreement with the surface manifestations of the rupture according to the field geology data and the results of geodetic inversion.

Moment tensor catalogue of microearthquakes in West Bohemia from 2008 to 2018

We present a unique catalogue of full moment tensors (MTs) of microearthquakes that occurred in West Bohemia, Czech Republic, in the period from 2008 to 2018. The catalogue is exceptional in several aspects: (1) it represents an extraordinary extensive dataset of more than 5.000 MTs, (2) it covers a long period of seismicity in the studied area, during which several prominent earthquake swarms took place, (3) the locations and retrieved MTs of microearthquakes are of a high accuracy. Additionally, we provide three-component records at the West Bohemia (WEBNET) seismic stations, the velocity model in the region, and the technical specification of the stations. The dataset is ideal for being utilized by a large community of researchers for various seismological purposes, e.g., for studies of (1) the migration of foci and the spatiotemporal evolution of seismicity, (2) redistribution of stress during periods of intense seismicity, (3) the interaction of faults, (4) the Coulomb stress along the faults and local stress anomalies connected to fault irregularities, (5) diffusivity of fluids along the activated faults, or (6) the time-dependent seismic risk due to the migration of seismicity in the region. In addition, the dataset is optimum for developing and testing new inversions for MTs and for tectonic stress. Since most of the earthquakes are non-shear, the dataset can contribute to studies of non-double-couple components of MTs and their relation to shear-tensile fracturing and/or seismic anisotropy in the focal zone.

Uncertainties in Seismic Moment Tensors Inferred from Differences between Global Catalogs

Catalogs of moment tensors form the foundation for a wide variety of seismological studies. However, assessing uncertainties in the moment tensors and the quantities derived from them is difficult. To gain insight, we compare 5000 moment tensors in the U.S. Geological Survey (USGS) and the Global Centroid Moment Tensor (Global CMT) Project catalogs for November 2015–December 2020 and use the differences to illustrate the uncertainties. The differences are typically an order of magnitude larger than the reported errors, suggesting that the errors substantially underestimate the uncertainty. The catalogs are generally consistent, with intriguing differences. Global CMT generally reports larger scalar moments than USGS, with the difference decreasing with magnitude. This difference is larger than and of the opposite sign from what is expected due to the different definitions of the scalar moment. Instead, the differences are intrinsic to the tensors, presumably in part due to different phases used in the inversions. The differences in double-couple components of source mechanisms and the fault angles derived from them decrease with magnitude. Non-double-couple (NDC) components decrease somewhat with magnitude. These components are moderately correlated between catalogs, with correlations stronger for larger earthquakes. Hence, small earthquakes often show large NDC components, but many have large uncertainties and are likely to be artifacts of the inversion. Conversely, larger earthquakes are less likely to have large NDC components, but these components are typically robust between catalogs. If so, these can indicate either true deviation from a double couple or source complexity. The differences between catalogs in scalar moment, source geometry, or NDC fraction of individual earthquakes are essentially uncorrelated, suggesting that the differences reflect the inversion rather than the source process. Despite the differences in moment tensors, the location and depth of the centroids are consistent between catalogs. Our results apply to earthquakes after 2012, before which many moment tensors were common to both catalogs.

Reliability tests of moment tensor inversion of anthropogenic seismicity

<p>Seismic moment tensor becomes part of basic seismic data processing. For anthropogenic events mostly common and available method to determine mechanism is amplitude inversion. However essential for correct amplitude inversion are good quality data. Factors commonly occurring in anthropogenic seismicity like high noise to signal ratio, low magnitude and shortage of seismic stations with unfavorable focal coverage can introduce undetected errors to inversion solution. In this work, synthetic tests for two seismic networks are presented to examine the reliability of P-wave first peak amplitude inversion for these areas. The synthetic tests of the noise influence on the results of full MT solutions were carried out for two surface networks monitoring anthropogenic seismicity: VERIS network in Vietnam and LUMINEOS network in Poland. Various mechanisms with double couple component variability from 10% to 100% were considered to take into account mechanisms caused by different types of human activity. High variability of solutions in tests shows that some spurious components cannot be avoided in full moment tensor solutions obtained for presented networks in certain cases.</p><p>This work was partially supported by research project no. 2017/27/B/ST10/01267, funded by the National Science Centre, Poland, under agreement no. UMO-2017/27/B/ST10/01267.</p>

Analysis of non-double-couple source mechanisms in an area of induced seismicity, West Texas (USA).

<p>In the scope to investigate the possible interactions between injected fluids, subsurface geology, stress field and triggering earthquakes, we investigate seismic source parameters related to the seismicity in West Texas (USA). The analysis of seismic moment tensor is an excellent tool to understand earthquake source process kinematics; moreover, changes in the fluid volume during faulting leads to existence of non-double-couple (NDC) components (Frohlich, 1994; Julian et al., 1998; Miller et al., 1998). The NDC percentage in the source constitutes the sum of absolute ISO and CLVD components so that %NDC= % ISO + %CLVD and %ISO+%CLVD+%DC=100%. It is currently known that the presence of NDC implies more complex sources (mixed shear-tensile earthquakes) correlated to fluid injections, geothermal systems and volcano-seismology where induced and triggered seismicity is observed.</p><p>With this hypothesis, we analyze the micro-earthquakes (M <2 .7) recorded by the Texas Seismological Network (TexNet) and a temporary network constituted by 40 seismic stations (equipped by either broadband or 3 component geophones). Our study area is characterized by Northwest-Southeast faults that follow the local stress/field (SH<sub>max</sub>) and the geological characteristic of the shallow basin structure of the study area. After a selection based on signal-to-noise ratio, we filter (1-50 Hz) the seismograms and estimate P-wave pulse polarities and the first P-wave ground displacement pulse in time domain. Then, we perform the full moment tensor analysis by using hybridMT technique (Andersen, 2001; Kwiatek et al., 2016) with a detailed 1D velocity model. The key parameter is the polarity/area of the first P-wave ground displacement pulse in time domain. Uncertainties of estimated moment tensors are expressed by normalized root-mean-square (RMS errors) between theoretical and estimated amplitudes (Vavricuk et al., 2014). We also evaluate the quality of the seismic moment tensors by bootstrap and resampling. In our preliminary results we obtain NDC percentage (in terms of %ISO and %CLVD components), Mw, seismic moment, P, T and B axes orientation for each source inverted.</p>

Complex magmatic-tectonic interactions during the 2020 Makushin Volcano, Alaska, earthquake swarm

<p>We investigate the processes driving<strong> </strong>a significant earthquake swarm that occurred between June and December 2020 on Unalaska Island, Alaska, ~12 km southeast of the summit of Makushin Volcano. The swarm was energetic, with two M>4 events that were widely felt by the population in Dutch Harbor, ~ 15 km west of the epicenters. This is the strongest seismic activity ever recorded at Makushin since instrumental monitoring began in 1996. To date, no eruptive activity or other surface changes have been observed at the volcano in satellite views, webcam images, GPS or InSAR. Seismic swarms close to volcanoes are often associated with the onset of unrest that can lead to eruption. However, determining whether they reflect magmatic rather than tectonic stresses is challenging. Here, we integrate information from space-time patterns of the hypocenters of the swarm earthquakes with their double-couple fault-plane solutions (FPS). We relocate swarm events using double-difference relocation techniques and a 3D velocity model. We find that most of the events cluster into two perpendicular lineaments with NW-SE and SW-NE orientations, but no apparent migration in time towards a preferred fault. On the one hand, the lack of temporal migration (with both faults slipping concurrently), and FPS for M3+ events consistent with regional stresses, seem to indicate a tectonic driving process. On the other hand, FPS for the lower-magnitude earthquakes have 90°-rotated P-axes perpendicular to the regional principal stress orientation, providing strong evidence for dike inflation/magma intrusion. Coulomb stress modeling indicates that the rotated FPS are best explained by an inflating dike to the SW of the swarm epicenters, in a zone of long-term elevated seismicity. This complex overlapping of regional and magmatic stresses is also evident in the statistical analysis of the sequence, which started as a main-shock/aftershock sequence with the first event having the largest magnitude, and evolved into a swarm sequence indicative of a more pronounced role of magmatic processes.</p>