Characterizing the large M4.7+ seismic events during the 2018 Kīlauea eruption is important to understand the complex subsurface deformation at the Kīlauea summit. The first 12 events (May 17 - May 26) are associated with long-duration seismic signals and the remaining 50 events (May 29 - August 02) are accompanied by large-scale caldera collapses. Resolving the source location and mechanism is challenging because of the shallow source depth, significant non double-couple components, and complex velocity structure. We demonstrate that combining multiple geophysical data from broadband seismometers, accelerometers and infrasound is essential to resolve different aspects of the seismic source. Seismic moment tensor solutions using near-field summit stations show the early events are highly volumetric. Infrasound data and particle motion analysis identify the inflation source as the Halema’uma’u reservoir. For the later collapse events, two independent moment tensor inversions using local and global stations consistently show that asymmetric slips occur on inward-dipping normal faults along the northwest corner of the caldera. While the source mechanism from May 29 onwards is not fully resolvable seismically using far-field stations, infrasound records and simulations suggest there may be inflation during the collapse. The summit events are characterized by both inflation and asymmetric slip, which are consistent with geodetic data. Based on the location of the slip and microseismicity, the caldera may have failed in a ‘see-saw’ manner: small continuous slips in the form of microseismicity on the southeast corner of the caldera, compensated by large slips on the northwest during the large collapse events.

We present a new 3-D seismic structural model of the eastern Indonesian region and its surroundings from full-waveform inversion (FWI) that exploits seismic data filtered at periods between 15 - 150 s. SASSY21 - a recent 3-D FWI tomographic model of Southeast Asia - is used as a starting model, and our study region is characterized by particularly good data coverage, which facilitates a more refined image. We use the spectral-element solver Salvus to determine the full 3-D wavefield, accounting for the fluid ocean explicitly by solving a coupled system of acoustic and elastic wave equations. This is computationally more expensive but allows seismic waves within the water layer to be simulated, which becomes important for periods ≤ 20 s. We investigate path-dependent effects of surface elevation (topography and bathymetry) and the fluid ocean on synthetic waveforms, and compare our final model to the tomographic result obtained with the frequently used ocean loading approximation. Furthermore, we highlight some of the key features of our final model - SASSIER22 - after 34 L-BFGS iterations, which reveals detailed anomalies down to the mantle transition zone, including a convergent double-subduction zone along the southern segment of the Philippine Trench, which was not evident in the starting model. A more detailed illumination of the slab beneath the North Sulawesi Trench reveals a pronounced positive wavespeed anomaly down to 200 km depth, consistent with the maximum depth of seismicity, and a more diffuse but aseismic positive wavespeed anomaly that continues to the 410 km discontinuity.

The style of convective force transmission to plates and strain-localization within and underneath plate boundaries remain debated. To address some of the related issues, we analyze a range of deformation indicators in southern California from the surface to the asthenosphere. Present-day surface strain rates can be inferred from geodesy. At seismogenic crustal depths, stress can be inferred from focal mechanisms and splitting of shear waves from local earthquakes via crack-dependent seismic velocities. At larger depths, constraints on rock fabrics are obtained from receiver function anisotropy, Pn and P tomography, surface wave tomography, and splitting of SKS and other teleseismic core phases. We construct a synthesis of deformation-related observations focusing on quantitative comparisons of deformation style. We find consistency with roughly N-S compression and E-W extension near the surface and in the asthenospheric mantle. However, all lithospheric anisotropy indicators show deviations from this pattern. Pn fast axes and dipping foliations from receiver functions are fault-parallel with no localization to fault traces and match post-Farallon block rotations in the Western Transverse Ranges. Local shear wave splitting inferences deviate from the stress orientations inferred from focal mechanisms in significant portions of the area. We interpret these observations as an indication that lithospheric fabric, developed during Farallon subduction and subsequent extension, has not been completely reset by present-day transform motion and may influence the current deformation behavior. This provides a new perspective on the timescales of deformation memory and lithosphere-asthenosphere interactions.

We present the first continental-scale seismic model of the lithosphere and underlying mantle beneath Southeast Asia obtained from adjoint waveform tomography (often referred to as full-waveform inversion or FWI), using seismic data filtered at periods from 20 - 150s. Based on >3,000h of analyzed waveform data gathered from ~13,000 unique source-receiver pairs, we image isotropic P-wave velocity, radially anisotropic S-wave velocity and density via an iterative non-linear inversion that begins from a 1-D reference model. At each iteration, the full 3-D wavefield is determined through an anelastic Earth, accommodating effects of topography, bathymetry and ocean load. Our data selection aims to maximize sensitivity to deep structure by accounting for body-wave arrivals separately. SASSY21, our final model after 87 iterations, is able to explain true-amplitude data from events and receivers not included in the inversion. The trade-off between inversion parameters is estimated through an analysis of the Hessian-vector product. SASSY21 reveals detailed anomalies down to the mantle transition zone, including multiple subduction zones. The most prominent feature is the (Indo-)Australian plate descending beneath Indonesia, which is imaged as one continuous slab along the 180-degree curvature of the Banda Arc. The tomography confirms the existence of a hole in the slab beneath Mount Tambora and locates a high S-wave velocity zone beneath northern Borneo that may be associated with subduction termination in the mid-late Miocene. A previously undiscovered feature beneath the east coast of Borneo is also revealed, which may be a signature of post-subduction processes, delamination or underthrusting from the formation of Sulawesi.

Volcanic arcs consist of many distinct vents that are ultimately fueled by the common process of melting in the subduction zone mantle wedge. Seismic imaging of crustal scale magmatic systems can provide insight into how melt is organized in the deep crust and eventually focused beneath distinct vents as it ascends and evolves. Here we investigate the crustal-scale structure beneath a section of the Cascades arc spanning four major stratovolcanoes: Mt. Hood, Mt. St. Helens, Mt. Adams, and Mt. Rainier, based on ambient noise interferometry measurements from 234 seismographs. Simultaneous inversion of Rayleigh and Love wave dispersion better constrain the isotropic shear velocity (Vs) and identify the unusual occurrence of radially anisotropic structures. Isotropic Vs shows two sub-parallel low-Vs zones at ~15-30 km depth with one connecting Mt. Rainier to Mt. Adams, and another connecting Mt. St. Helens to Mt. Hood, which are interpreted as deep crustal magma reservoirs containing up to ~2.5-6% melt, assuming near-equilibrium melt geometry. Negative radial anisotropy is prevalent in this part of the Cascadia margin, but is interrupted by positive radial anisotropy extending vertically beneath Mt. Adams and Mt. Rainier at ~10-30 km depth and weaker positive anisotropy beneath Mt. St. Helens with a west dipping. The positive anisotropy regions are adjacent to rather than co-located with the isotropic low-Vs anomalies. Ascending melt that stalled and mostly crystallized in sills with possible compositional difference from the country rock may explain the near-average Vs and positive radial anisotropy adjacent to the active deep crustal magma reservoirs.