Detection of slow slip events using wavelet analysis of GNSS recordings
In many places, tectonic tremor is observed in relation to slow slip and can be used as a proxy to study slow slip events of moderate magnitude where surface deformation is hidden in Global Navigation Satellite System (GNSS) noise. However, when no clear relationship between tremor and slow slip occurrence is observed, these methods cannot be applied, and we need other methods to be able to better detect and quantify slow slip. Wavelets methods such as the Discrete Wavelet Transform (DWT) and the Maximal Overlap Discrete Wavelet Transform (MODWT) are mathematical tools for analyzing time series simultaneously in the time and the frequency domain by observing how weighted differences of a time series vary from one period to the next. We use wavelet methods to analyze GNSS time series of slow slip events in Cascadia. We use detrended GNSS data, apply the MODWT transform and stack the wavelet details over several nearby GNSS stations. As an independent check on the timing of slow slip events, we also compute the cumulative number of tremor in the vicinity of the GNSS stations, detrend this signal, and apply the MODWT transform. In both time series, we can then see simultaneous waveforms whose timing corresponds to the timing of slow slip events. We assume that there is a slow slip event whenever there is a positive peak followed by a negative peak in the wavelet signal. We verify that there is a good agreement between slow slip events detected with only GNSS data, and slow slip events detected with only tremor data. The wavelet-based detection method effectively detects events of magnitude higher than 6 as determined by independent event catalogs. As a demonstration of using the wavelet analysis in a region without significant tremor, we also analyze GNSS data from New Zealand and detect slow slip events that are spatially and temporally close to those detected previously by other studies.