Full Moment Tensor Inversion Research Articles

Crustal seismicity and updated stress mapping in the Bolivian Orocline, Central Andes

The Central Andes region, with significant seismic activity, has distinct faulting mechanisms across its varied topography. Notably, the sub-Andean province predominantly features reverse faulting, whereas the high Andean plateau shows a predominance of normal and strike-slip faults. Despite the importance of this extensive orogenic plateau, Bolivia remains under-documented in seismic studies. We used data from recently installed seismic stations in Bolivia and regional stations in Brazil to determine 13 new focal mechanisms of shallow crustal earthquakes in Bolivia, some of them felt in major cities. Employing probabilistic full waveform moment tensor inversion and P-wave polarities, we mapped the stress distribution across the Central Andes. The main patterns of faulting mechanisms were: reverse faulting along the NW-SE trending sub-Andean belt north of the Orocline (north of Cochabamba) with NE-SW compression; reverse faulting along the N-S trending sub-Andean belt south of the Orocline (south of Santa Cruz) with ∼ E-W oriented compression; strike-slip faulting within the Eastern Cordillera with NE-SW P axes. The results indicate that the sub-Andean belt experiences compressional forces perpendicular to the front of the Andean plateau, arising from both the spreading stresses of the plateau and the broader plate-wide compression. On the other hand, the high plateau (Altiplano) is characterized by normal and strike-slip mechanisms, suggesting a dynamic equilibrium between local extensional gravitational stresses and regional compressional forces due to the Nazca plate convergence.

Injection-induced fault slip and associated seismicity in the lab: Insights from source mechanisms, local stress states and fault geometry

Probing source mechanisms of natural and induced earthquakes is a powerful tool to unveil associated rupture kinematics. The source processes of failure and slip instability driven by stress loading are affected by fault geometry, but the source ruptures of injection-induced seismicity in relation to fault structures and local stress states remain poorly understood. We have conducted a series of fault reactivation and slip experiments on sandstone samples containing faults with different surface roughness (smooth saw-cut fault and fractured rough fault). We impose progressive fluid injection to induce fault slip, and simultaneously monitor the associated acoustic emission (AE) activity. Using high-resolution AE recordings, we perform full moment tensor inversion of all located AE sources, and investigate the changes of AE source characteristics associated with induced fault slip and their relation to fault roughness. For the complex and rough fault, we observe significant non-double-couple components of AE sources and a high degree of focal mechanism heterogeneity. The temporal changes of AE mechanisms associated with injection-induced fault slip on the smooth fault reveal increasing proportions of double-couple components and decreasing variability of AE focal mechanisms when approaching the onset of slip events. The observed inconsistency between the nodal planes of AE sources and the macroscopic fault plane orientation is attributed to the development of secondary fracture networks surrounding the principal slip surface. We analyze changes in the magnitude-frequency characteristics and source mechanisms of AEs with fault-normal distance, showing that for the smooth (mature) fault, Gutenberg–Richter b-value of on-fault seismicity is lower and focal mechanisms are less heterogeneous, compared to off-fault seismicity. Our results emphasize the important role of roughness-related changes in local fault geometry and associated stress heterogeneity for source mechanisms and rupture kinematics of injection-induced seismicity.

Active Deformation Constraints on the Nubia‐Somalia Plate Boundary Through Heterogenous Lithosphere of the Turkana Depression

AbstractThe role of lithospheric heterogeneities, presence or absence of melt, local and regional stresses, and gravitational potential energy in strain localization in continental rifts remains debated. We use new seismic and geodetic data to identify the location and orientation of the modern Nubia‐Somalia plate boundary in the 300‐km‐wide zone between the southern Main Ethiopian Rift (MER) and Eastern Rift (ER) across the Mesozoic Anza rift in the Turkana Depression. This region exhibits lithospheric heterogeneity, 45 Ma‐Recent magmatism, and more than 1,500 m of base‐level elevation change, enabling the assessment of strain localization mechanisms. We relocate 1716 earthquakes using a new 1‐D velocity model. Using a new local magnitude scaling with station corrections, we find 1 ≤ ML ≤ 4.5, and a b‐value of 1.22 ± 0.06. We present 59 first motion and 3 full moment tensor inversions, and invert for opening directions. We use complementary geodetic displacement vectors and strain rates to describe the geodetic strain field. Our seismic and geodetic strain zones demonstrate that only a small part of the 300 km‐wide region is currently active; low elevation and high‐elevation regions are active, as are areas with and without Holocene magmatism. Variations in the active plate boundary's location, orientation and strain rate appear to correspond to lithospheric heterogeneities. In the MER‐ER linkage zone, a belt of seismically fast mantle lithosphere generally lacking Recent magmatism is coincident with diffuse crustal deformation, whereas seismically slow mantle lithosphere and Recent magmatism are characterized by localized crustal strain; lithospheric heterogeneity drives strain localization.

NA-GCAP: A New Full Moment Tensor Inversion Method Based on the Renewed GRTM and NA Method

Abstract Moment tensor inversion plays a crucial role in determining earthquake types, magnitude, and source geometry. Compared to polarity and amplitude methods, full-waveform approaches provide more comprehensive constraints on the complex full moment tensor (FMT). In this study, we propose a novel FMT inversion method named neighborhood algorithm-generalized cut-and-paste (NA-GCAP) method. Similar to the “cut-and-paste” method, our approach divides seismograms into Pnl and surface-wave segments. To enhance inversion efficiency, we employ the NA in conjunction with a newly developed strategy for efficient Green’s functions computation based on the renewed fast generalized reflection and transmission method (GRTM). Our method requires only a single forward computation and stores ten Green’s functions, significantly improving computational efficiency. We validate the robustness of our approach through synthetic tests using a velocity model perturbed by 5% relative to the input model. Furthermore, we apply the NA-GCAP method to the 2019 Changning earthquake sequence comprising 16 earthquakes, where twelve events exhibit double-couple (DC) components larger than 0.95, indicating a simple dislocation source, and four events display significant non-DC components. Our results align well with the previous studies and demonstrate the potential for widespread application to other earthquake sequences in the future.

Improved full waveform moment tensor inversion of Cratonic intraplate earthquakes in southwest Australia

SUMMARYIn contrast to global observations in stable continental crust, the present-day orientation of the maximum horizontal stress in Western Australia is at a high angle to plate motion, suggesting that in addition to large-scale plate driving forces, local factors also play an important role in stress repartitioning. As a reliable stress indicator, full waveform moment tensor solutions are calculated for earthquakes that occurred between 2010 and 2018 in the southern Yilgarn Craton and the adjacent Albany-Fraser Orogen in southwestern Australia. Due to regional velocity heterogeneities in the crust, we produced two geographically distinct shear wave velocity models by combining published crustal velocity models with new ambient noise tomography results. We applied a full waveform inversion technique to 15 local earthquakes and obtained 10 robust results. Three of these events occurred near Lake Muir in the extreme south of the study area within the Albany-Fraser Orogen. The focal mechanism of the 16th September 2018 Lake Muir event is thrust; two ML≥ 4.0 aftershocks are normal and strike-slip. Our results are consistent with field observations, fault orientations inferred from aeromagnetic data and surface displacements mapped by Interferometric Synthetic Aperture Radar which are all consistent with reactivation of existing faults. The other seven solutions are in the southeastern Yilgarn Craton. These solutions show that the faulting mechanisms are predominantly thrust and strike-slip. This kinematic framework is consistent with previous studies that linked active seismicity in the Yilgarn Craton to the reactivation of the NNW–SSE oriented Neoarchean structures by an approximately E–W oriented regional stress field. Our results suggest that the kind of faulting that occurs in southwest Australia is critically dependent on the local geological structure. Thrust faulting is the dominant rupture mechanism, with some strike-slip faulting occurring on favourably oriented structures.

Orientations of Broadband Stations of the KOERI Seismic Network (Turkey) from Two Independent Methods: P- and Rayleigh-Wave Polarization

Abstract The correct orientation of seismic sensors is critical for studies such as full moment tensor inversion, receiver function analysis, and shear-wave splitting. Therefore, the orientation of horizontal components needs to be checked and verified systematically. This study relies on two different waveform-based approaches, to assess the sensor orientations of the broadband network of the Kandilli Observatory and Earthquake Research Institute (KOERI). The network is an important backbone for seismological research in the Eastern Mediterranean Region and provides a comprehensive seismic data set for the North Anatolian fault. In recent years, this region became a worldwide field laboratory for continental transform faults. A systematic survey of the sensor orientations of the entire network, as presented here, facilitates related seismic studies. We apply two independent orientation tests, based on the polarization of P waves and Rayleigh waves to 123 broadband seismic stations, covering a period of 15 yr (2004–2018). For 114 stations, we obtain stable results with both methods. Approximately, 80% of the results agree with each other within 10°. Both methods indicate that about 40% of the stations are misoriented by more than 10°. Among these, 20 stations are misoriented by more than 20°. We observe temporal changes of sensor orientation that coincide with maintenance work or instrument replacement. We provide time-dependent sensor misorientation correction values for the KOERI network in the supplemental material.

Seismic source characterization of the Arabian Peninsula and Zagros Mountains from regional moment tensor and coda envelopes

Reliable estimates of moment magnitude and source mechanism for seismic events in the Middle East can be challenging due to a small number of openly available stations, the complex tectonic setting, and regions of high attenuation. Access to high-quality waveform data from well-calibrated regional seismic stations is fundamental in producing robust and stable estimates of earthquake source parameters, particularly when measurements of absolute ground motion amplitudes are required. Earthquake source mechanism and moment magnitude are invaluable information in the assessment of seismic hazard, plate motions, and the characterization of faults and regional stress field. The expansion of the Kingdom of Saudi Arabia’s national seismic network in the past several decades provides an opportunity to develop the capabilities of routine focal mechanism and moment magnitude estimations. In this study, we use time-domain full waveform moment tensor inversion and coda envelope–derived amplitude measurements to solve for earthquake source mechanism, moment magnitudes, and their source type. We compared the moment magnitudes calculated from the two methods and publicly available earthquake catalogs and discuss the implications of the obtained source parameters.

Towards a regional, automated full moment tensor inversion for medium to large magnitude events in the Iranian plateau

Considering the high seismicity rate of the Iranian plateau, the development and implementation of an automated routine for near-real-time determination of earthquake point source parameters would be extremely beneficial, e.g., for dedicated seismological studies, constraining the regional stress field, and shakemap assessment. We discuss the implementation of a fast, automated regional moment tensor inversion routine, based on SeisComP3 and KIWI tool packages. To assess the routine efficiency, we used a dataset including 90 earthquakes spatially distributed within all seismotectonic provinces of Iran, recorded at local to regional distances, with magnitudes in the range MN 4.5–7.3. Resolved moment tensor solutions and source parameters are in most cases in good agreement with reference ones from global and regional catalogs. However, we identify consistent differences in magnitudes and depths, with respect to global catalogs, which can be attributed to a higher resolution of our approach, using local and regional data and considering regional velocity models. The depth estimations are shallower than in reference global catalogs but compatible with the tectonic settings of the Iranian region. The analysis of the selected dataset is accompanied by a test run of automated processing for 8 months of online operation, from July 2019 to February 2020. The processing time in online routine, with an average of ~ 23 min, mostly depends on the magnitude, number of stations, azimuthal coverage, data quality, and complexity of the crustal model, so that computation times may vary significantly for earthquakes located in different Iranian provinces.

GCAPjoint, A Software Package for Full Moment Tensor Inversion of Moderately Strong Earthquakes with Local and Teleseismic Waveforms

AbstractEarthquake moment tensors and focal depths are crucial to assessing seismic hazards and studying active tectonic and volcanic processes. Although less powerful than strong earthquakes (M 7+), moderately strong earthquakes (M 5–6.5) occur more frequently and extensively, which can cause severe damages in populated areas. The inversion of moment tensors is usually affected by insufficient local waveform data (epicentral distance <5°) in sparse seismic networks. It would be necessary to combine local and teleseismic data (epicentral distance 30°–90°) for a joint inversion. In this study, we present the generalized cut-and-paste joint (gCAPjoint) algorithm to jointly invert full moment tensor and centroid depth with local and teleseismic broadband waveforms. To demonstrate the effectiveness and explore the limitations of this algorithm, we perform case studies on three earthquakes with different tectonic settings and source properties. Comparison of our results with global centroid moment tensor and other catalog solutions illustrates that both non-double-couple compositions of the focal mechanisms and centroid depths can be reliably recovered for very shallow (<10 km) earthquakes and intermediate-depth events with this software package.

Full moment tensor inversion constrained by double-couple focal mechanism for induced seismicity*

In this study, we propose a new method to determine full moment tensor solution for induced seismicity. This method generalizes the full waveform matching algorithm we have developed to determine the double-couple (DC) focal mechanism based on the neighbourhood algorithm. One major difference between the new method and the former one is that we adopt a new misfit function to constrain the candidate moment tensor solutions with respect to a reference DC solution in addition to other misfit terms characterizing the waveform matching. Through synthetic tests using a real passive seismic survey geometry, the results show the new constraint can help better recover the DC components of inverted moment tensors. We further investigate how errors in the velocity model and source location affect the moment tensor solution. The synthetic test results indicate that the constrained inversion is robust in recovering both the DC and non-DC components. We also test the proposed method on several real induced events in an oil/gas field in Oman using the same observation system as synthetic tests. While it is found that the full moment tensor solutions without using the DC constraints have much larger non-DC components than solutions with the DC constraints, both solutions are able to fit the observed waveforms at similar levels. The synthetic and real test results suggest the proposed DC constrained inversion method can reliably retrieve full moment tensor solutions for the induced seismicity.

Discriminating a Nuclear Blast from A Natural Seismic Event Using A Full Moment-Tensor Inversion Method

Identification of seismic signals if they are generated by tectonic or volcanic origin or underground nuclear tests is of primary importance for impacts on environmental issues. The recent blast on September 3, 2017 has been in question due to its potential destruction on the global environment. Here, we applied a method of full moment-tensor inversions to analyse broadband waveforms generated by such an explosion. The signals were recorded by a network of 4 regional seismic stations (MAJO, INCN, MDJ, BJT), situated between 360 and 1,100 km away from the explosion epicenter. The results showed that a variance reduction of 72.6%, providing high confidence for all values obtained from the method used: a seismic moment of Mo = 6.25 · 1024 dyne cm, corresponding to a moment magnitude of Mw = 5.8, with a dominant volumetric component (ISO) of 50% relatively compared to a 6.8% double-couple (DC) part and a 43.2% compensated linear vector dipole (CLVD) part, and a centroid-depth of 1.0 km. The relative dominance of the ISO seismic component over the other two components, along with the shallow source extracted from the cross-correlation and source mechanism analysis, suggests that the September 3, 2017 explosion was induced by a suspected nuclear test.

1-D velocity model for the North Korean Peninsula from Rayleigh wave dispersion of ambient noise cross-correlations

Monitoring seismic activity in the north Korean Peninsula (NKP) is important not only for understanding the characteristics of tectonic earthquakes but also for monitoring anthropogenic seismic events. To more effectively investigate seismic properties, reliable seismic velocity models are essential. However, the seismic velocity structures of the region have not been well constrained due to a lack of available seismic data. This study presents 1-D velocity models for both the inland and offshore (western East Sea) of the NKP. We constrained the models based on the results of a Bayesian inversion process using Rayleigh wave dispersion data, which were measured from ambient noise cross-correlations between stations in the southern Korean Peninsula and northeast China. The proposed models were evaluated by performing full moment tensor inversion for the 2013 Democratic People’s Republic of Korea (DPRK) nuclear test. Using the composite model consisting of both inland and offshore models resulted in consistently higher goodness of fit to observed waveforms than previous models. This indicates that seismic monitoring can be improved by using the proposed models, which resolve propagation effects along different paths in the NKP region.

Investigating Factors Influencing Moment Tensor Inversion of Induced Seismicity in Virtual IoT

The Internet of Things (IoT) for seismicity in underground mine consists of a combination of sensor network hardware and software, which allows the collection of seismic data and data processing for basic event parameters. A seismic moment tensor is one of the essential parameters in evaluating the seismic hazard. The reliability of the resolved seismic moment tensor depends on numerous factors, and it is important to understand the influence of these disturbances. We focused on influencing factors, including the azimuthal coverage, the accuracy of source coordinate, the mismodeling of velocity structure, and the noise contamination. A set of moment tensor inversions using synthetic data simulated for a double couple source and a superimposed source from virtual IoT was carried out to test the effects of different factors. We concluded that the sensitivity of the moment tensor inversion to the azimuthal coverage is closely related to the type of the source. For the shear event, the inversion does not require a good azimuthal coverage, but for the non-pure shear event, at least 90° azimuthal coverage is required to get the reliable solution. The full moment tensor is not sensitive to the location accuracy when the source error is less than 60 m, regardless of the type of source. But the double couple deviation increases with the growing error of the location, especially for the fault-slip-related or the shear rupture-dominated event. The erroneous velocity models show nearly no influence on the full moment tensor inversion results. The fit of amplitude spectra does not require a precise alignment of the observed and the synthetic waveforms and is less dependent on the chosen velocity model. The superimposed source is more sensitive to noise than the pure shear source. Although the noise does not necessarily affect the fault plane solution of the event, it is capable to lead an inaccurate seismic moment tensor and mislead the interpretation of the source type.

Source characteristics of the upper mantle 21 May, 2014 Bay of Bengal earthquake of $${{\varvec{M}}}_{\!{{\varvec{w}}}}$$ M w 5.9

We measure source parameters for the 21 May, 2014 Bay of Bengal earthquake through inversion modeling of S-wave displacement spectra from radial–transverse–vertical (RTZ) components recorded at ten broadband stations in the eastern Indian shield. The average source parameters are estimated using estimates from seven near stations (within epicentral distances \({\le }500\,\hbox {km}\)). The average seismic moment and source radius are determined to be \(1.0{\times }10^{18}\,\hbox {N-m}\) and 829 m, respectively, while average stress drop is found to be 76.5 MPa. The mean corner frequency and moment magnitude are calculated to be \(1.6\pm 0.1\) and \(5.9\pm 0.2\) Hz, respectively. We also estimated mean radiated energy and apparent stress, which are found to be \(6.1{\times }10^{13}\) joules and 1.8 MPa, respectively. We observe that mean \(E_{s}/M_{o}\) estimate of \(5.5{\times }10^{-5}\) is found to be larger than the global average for oceanic strike-slip events. This observation along with large stress drop and apparent stress estimates explains the observed remarkably felt intensity data of the 2014 event. The full waveform moment tensor inversion of the band-passed (0.03–0.12 Hz) broadband displacement data suggests the best fit for the multiple point sources on a plane located at 65 km depth, with a moment magnitude 6.4, and a focal mechanism with strike \(318^{\mathrm{o}}\), dip \(87^{\mathrm{o}}\), and rake \(34^{\mathrm{o}}\).

Magmatic or Not Magmatic? The 2015–2016 Seismic Swarm at the Long-Dormant Jailolo Volcano, West Halmahera, Indonesia

Seismic swarms close to volcanoes often signal the onset of unrest. Establishing whether magma is the culprit and the unrest can be flagged as magmatic may be challenging. Here we analyze the spatio-temporal pattern of a seismic swarm that occurred in November 2015 – February 2016 around Jailolo volcano, a long-dormant and poorly studied volcano located on Halmahera island, North Moluccas, Indonesia. The swarm included four Mw>5 earthquakes and hundreds of events were felt by the population. We relocate the earthquakes using both the Indonesian Seismic Network and single-station location techniques. We find that the earthquakes cluster in a narrow strip, stretching 5 km E–W and 20 km N–S, migrating southward away from Jailolo volcano at ~10 km/d. We investigate the source mechanisms of the largest earthquakes via full moment tensor inversion. The non-double-couple component is around 50%, such that the earthquakes, besides normal faulting, included a relatively large opening component. After a thorough examination of the possible causes of the Jailolo swarm we conclude that a laterally propagating dike of tens of millions of cubic meters is the triggering mechanism for the seismicity. The swarm marks the first documented magmatic unrest at Jailolo. We find that there is a probability >0.1 that the unrest will last for more than two years. This magmatic unrest calls for the classification of Jailolo volcano as active, and for an urgent assessment of the associated volcanic hazard.

Waveform-based Bayesian full moment tensor inversion and uncertainty determination for the induced seismicity in an oil/gas field

Waveform-based Bayesian full moment tensor inversion and uncertainty determination for the induced seismicity in an oil/gas field

Aftershocks of the 2012Mw 7.6 Nicoya, Costa Rica, Earthquake and Mechanics of the Plate Interface

Abstract Subduction of the Cocos plate beneath the Nicoya Peninsula, Costa Rica, generates large underthrusting earthquakes with a recurrence interval of about 50 yrs. The most recent of these events occurred on 5 September 2012 ( M w 7.6). A vigorous sequence of more than 6400 aftershocks was recorded by a local seismic network within the first four months after the mainshock. We determine locations and focal mechanisms for as many aftershocks as possible with M ≥1.5 occurring within the first nine days of the mainshock, all aftershocks with M ≥3 through the end of 2012, and all events with M ≥4 through the end of 2015. We determine faulting geometries using regional full waveform moment tensor (MT) inversion for the largest events ( M ≥4) and P ‐wave first‐motion polarities for smaller events, producing a mechanism catalog with 347 earthquakes. Sixty percent of these events are identified as underthrusting, and their locations are compared with spatial distributions of mainshock slip, afterslip, prior interplate seismicity, and slow‐slip phenomena to better understand the mechanical behavior of the plate interface. Most of the aftershocks on the megathrust occur up‐dip of the coseismic slip, where afterslip is large, and between coseismic slip and shallow slow‐slip patches. The pattern of interplate seismicity during the interseismic period is similar to that for the aftershocks but does not extend to as great a depth. The coseismic slip extends even deeper than the interplate aftershocks, suggesting that the mainshock ruptured a strongly locked patch driving down‐dip slip into the conditionally stable part of the deep plate interface that also hosts slow slip. About 80% of the aftershocks have one nodal plane oriented favorably to promote failure from static stress changes following the mainshock and early afterslip, whereas most others occur in regions of large afterslip. Electronic Supplement: Tables of hypocenter location and focal mechanisms and selection criteria for underthrusting events and figures showing the Gutenberg–Richter distribution, regional moment tensor inversion, spatial distribution of earthquake activity, and temporal distribution of underthrusting events.

Modeling of Source Parameters of the 15 December 2015 Deogarh Earthquake of Mw 4.0

ABSTRACT We herein present source parameters and focal mechanism of a rare cratonic upper crustal earthquake of Mw4.0, which occurred at 8 km depth (centroid depth) below a region near Deogarh, Jharkhand. For our study, we used broadband waveform data from a seismic network of 15 three-component seismographs in the eastern Indian craton. The average seismic moment, moment magnitude and source radius are estimated to be 1.1 x 1015 N-m, 4.0 and 180.6 m, respectively. The high average stress drop of 14.27 MPa could be attributed to its lower-crustal origin. The mean corner frequency is calculated to be 4.1 Hz. To study the source mechanism, we perform a deviatoric constrained full waveform moment tensor inversion of multiple point sources on the band-passed (0.06 - 0.14 Hz) broadband displacement data of the Deogarh event, using ISOLA software. The best fit is obtained for the source at 8 km centroid depth, with a moment magnitude 3.7, and a right-lateral strike-slip mechanism with strike 162°, dip 72° and rake 169°. The P-axis orients N24°E, which is parallel to the direction of the absolute plate motion direction of Indian plate, while T-axis orients E-W, which is parallel to the strike of the pre-existing Damodar Graben (DG) of Gondwana age. The occurrence of this earthquake is attributed to the neotectonic reactivation of a fault associated with the E-W trending DG shear zone.

A one-dimensional seismic model for Uturuncu volcano, Bolivia, and its impact on full moment tensor inversions

Using receiver functions, Rayleigh wave phase velocity dispersion determined from ambient noise and teleseismic earthquakes, and Rayleigh wave horizontal to vertical ground motion amplitude ratios from earthquakes observed across the PLUTONS seismic array, we construct a one-dimensional (1‑D) S-wave velocity (Vs) seismic model with uncertainties for Uturuncu volcano, Bolivia, located in the central Andes and overlying the eastward-subducting Nazca plate. We find a fast upper crustal lid placed upon a low-velocity zone (LVZ) in the mid-crust. By incorporating all three types of measurements with complimentary sensitivity, we also explore the average density and Vp/Vs (ratio of P-wave to S-wave velocity) structures beneath the young silicic volcanic field. We observe slightly higher Vp/Vs and a decrease in density near the LVZ, which implies a dacitic source of the partially molten magma body. We exploit the impact of the 1-D model on full moment tensor inversion for the two largest local earthquakes recorded (both magnitude ∼3), demonstrating that the 1-D model influences the waveform fits and the estimated source type for the full moment tensor. Our 1-D model can serve as a robust starting point for future efforts to determine a three-dimensional velocity model for Uturuncu volcano.

Moment Tensor Analysis of Very Shallow Sources

An issue for moment tensor (MT) inversion of shallow seismic sources is that some components of the Green’s functions have vanishing amplitudes at the free surface, which can result in bias in the MT solution. The effects of the free surface on the stability of the MT method become important as we continue to investigate and improve the capabilities of regional full MT inversion for source‐type identification and discrimination. It is important to understand free‐surface effects on discriminating shallow explosive sources for nuclear monitoring purposes. It may also be important in natural systems that have very shallow seismicity, such as volcanic and geothermal systems. We examine the effects of the free surface on the MT via synthetic testing and apply the MT‐based discrimination method to three quarry blasts from the HUMMING ALBATROSS experiment. These shallow chemical explosions at ∼10 m depth and recorded up to several kilometers distance represent rather severe source–station geometry in terms of free‐surface effects. We show that the method is capable of recovering a predominantly explosive source mechanism, and the combined waveform and first‐motion method enables the unique discrimination of these events. Recovering the design yield using seismic moment estimates from MT inversion remains challenging, but we can begin to put error bounds on our moment estimates using the network sensitivity solution technique (Ford et al. , 2010). Online Material: Figures showing synthetic tests for a pure explosion and a composite source at local distances and table of moment tensor components.