  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.
