Reduction of non-tidal oceanographic fluctuations in ocean-bottom pressure records of DONET using principal component analysis to enhance transient tectonic detectability

Ocean bottom pressure-gauge (OBP) records play an essential role in seafloor geodesy. Oceanographic fluctuations in OBP data, however, pose as a significant noise source in seafloor transient crustal deformation observations, including slow slip events (SSEs), making it crucial to evaluate them quantitatively. To extract the significant fluctuation phenomena common to multiple observation networks, including oceanographic fluctuations and tectonic signals, we applied principal component analysis (PCA) to the 3-year Dense Oceanfloor Network System for Earthquakes and Tsunamis (DONET) OBP time series for 40 stations during 2016-2019. PCA could separate several oceanographic signals based on the characteristics of their spatial distributions, although evident transient tectonic signals could not be confirmed from the observed pressure records during this observed period. The spatial distribution of the first four principal components (PCs) reflected the common component, inclined component along sea depth, longitude component, and parabola-like pattern, respectively. By subtracting each PC (in particular, PC-2 and PC-4) from the time series, we could significantly reduce the sea depth dependence of OBP records, which has been highlighted in several previous studies and is also evident in this region. We interpreted PCs 2-4 as the reflection of the strength and meandering of ocean geostrophic currents based on a comparison with the PC spatial distribution of the numerical oceanographic models. In addition, to evaluate the ability of PCA to separate transient tectonic signal from OBP time series, including oceanographic fluctuations, we conducted a synthetic ramp assuming an SSE by rectangular fault and then applied PCA. The assumed synthetic tectonic signal could be separated from the oceanographic signals and included in the principal component independently depending on its amplitude, suggesting that the spatial distribution of each PC would change if the amplitude of the synthetic signal were sufficiently large. We propose a transient event-detection method based on the spatial distribution difference of a specific PC with or without a tectonic signal. We used the normalized inner product (NIP) between these PCs as the indicator of their similarities. This method can detect transient tectonic signals more significantly than the moment-magnitude scale of 5.9 from OBP records.