Abstract. We resolve a previously unrecognized shallow subducting seamount from a re-processed multichannel seismic profile crossing the 1994 Mw 7.8 Java tsunami earthquake rupture area. Seamount subduction occurs where the overriding plate experiences uplift by lateral shortening and vertical thickening. Pronounced back-thrusting at the landward slope of the forearc high and the formation of splay faults branching off the landward flank of the subducting seamount are observed. The location of the seamount in relation to the 1994 earthquake hypocentre and its co-seismic slip model suggests that the seamount acted as a seismic barrier to the up-dip co-seismic rupture propagation of this moderate-size earthquake.

Abstract. We resolve a previously unrecognized shallow subducting seamount from a re-processed multichannel seismic depth image crossing the 1994 M7.8 Java tsunami earthquake slip area. Seamount subduction is related to the uplift of the overriding plate by lateral shortening and vertical thickening, causing pronounced back-thrusting at the landward slope of the forearc high and the formation of splay faults branching off the landward flank of the subducting seamount. The location of the seamount in relation to the 1994 earthquake hypocentre and its co-seismic slip model suggests that the seamount acted as a seismic barrier to the up-dip co-seismic rupture propagation of this moderate size earthquake. The wrapping of the co-seismic slip contours around the seamount indicates that it diverted rupture propagation, documenting the control of forearc structures on seismic rupture.

No abstract

Abstract. Accurate subsurface velocity models are crucial for geological interpretations based on seismic depth images. Seismic reflection tomography is an effective iterative method to update and refine a preliminary velocity model for depth imaging. Based on residual move-out analysis of reflectors in common image point gathers, an update of the velocity is estimated by a ray-based tomography. To stabilize the tomography, several preconditioning strategies exist. Most critical is the estimation of the depth error to account for the residual move-out of the reflector in the common image point gathers. Because the depth errors for many closely spaced image gathers must be picked, manual picking is extremely time-consuming, human biased, and not reproducible. Data-driven picking algorithms based on coherence or semblance analysis are widely used for hyperbolic or linear events. However, for complex-shaped depth events, purely data-driven picking is difficult. To overcome this, the warping method named non-rigid matching is used to estimate a depth error displacement field. Warping is used, for example, to merge photographic images or to match two seismic images from time-lapse data. By matching a common image point gather against its duplicate that has been shifted by one offset position, a locally smooth-shaped displacement field is calculated for each data sample by gather matching. Depending on the complexity of the subsurface, sample tracking through the displacement field along predefined horizons or on a simple regular grid yields discrete depth error values for the tomography. The application to a multi-channel seismic line across the Sunda subduction zone offshore Lombok island, Indonesia, illustrates the approach and documents the advantages of the method to estimate a detailed velocity structure in a complex tectonic regime. By incorporating the warping scheme into the reflection tomography, we demonstrate an increase in the velocity resolution and precision by improving the data-driven accuracy of depth error picks with arbitrary shapes. This approach offers the possibility to use the full capacities of tomography and further leads to more accurate interpretations of complex geological structures.

The updip limit of seismic rupture during a megathrust earthquake exerts a major control on the size of the resulting tsunami. Offshore Northern Chile, the 2014 Mw 8.1 Iquique earthquake ruptured the plate boundary between 19.5° and 21°S. Rupture terminated under the mid-continental slope and did not propagate updip to the trench. Here, we use state-of-the-art seismic reflection data to investigate the tectonic setting associated with the apparent updip arrest of rupture propagation at 15 km depth during the Iquique earthquake. We document a spatial correspondence between the rupture area and the seismic reflectivity of the plate boundary. North and updip of the rupture area, a coherent, highly reflective plate boundary indicates excess fluid pressure, which may prevent the accumulation of elastic strain. In contrast, the rupture area is characterized by the absence of plate boundary reflectivity, which suggests low fluid pressure that results in stress accumulation and thus controls the extent of earthquake rupture. Generalizing these results, seismic reflection data can provide insights into the physical state of the shallow plate boundary and help to assess the potential for future shallow rupture in the absence of direct measurements of interplate deformation from most outermost forearc slopes.