Optimal Multichannel Stretch Factors for Estimating Changes in Seismic Velocity: Application to the 2012 Mw 7.8 Haida Gwaii Earthquake

Abstract

We propose new methods for assessing temporal changes in seismic velocity using the S-wave coda for repeating earthquakes and cross-correlation functions of ambient noise. For a pair of seismic waveforms representing a common source–receiver path, the relative change in path-averaged velocity over the corresponding time interval is directly proportional to the factor by which one waveform needs to be stretched or compressed with respect to the other to achieve maximum coherence. For an arbitrary number of waveforms, initial pair-wise stretch factors determined through standard approaches can be improved through solution of an overdetermined system and further refined through an iterative approach exploiting the singular value decomposition to minimize rank of the stretched waveform section. We apply this combined approach to both repeating earthquakes and ambient noise correlations for Haida Gwaii in western Canada, the site of a Mw 7.8 thrust earthquake in 2012. Optimal stretch factors for repeating earthquake families indicate that path-averaged S velocities dropped by up to 0.16% after the earthquake. Ambient noise correlations indicate that velocities dropped by between 0.26% and 0.39%, which we interpret to be more pronounced in the uppermost levels of the crust. We explore these results in terms of changes in crustal porosity and hydrogeologic conditions by considering the observation that hot spring activity on Haida Gwaii ceased following the 2012 mainshock and recovered over the next several years.