Studies of Seismic Velocities in Subduction Zones from Continuous Ocean Bottom Seismometer Data

<p><b>This thesis uses continuous ambient noise data recorded by Ocean Bottom Seismometers (OBSs) to study seismic velocities in the upper crust of the overriding plate. The first and second projects (Chapters 3 and 4) focus on temporal seismic velocity variations in the northern Hikurangi subduction zone offshore the North Island, New Zealand, while the third project (Chapter 5) focuses on shear wave velocities in the southwestern Okinawa Trough offshore northeastern Taiwan. In the first project (Chapters 3), we investigate a region of frequent slow slip events (SSEs) offshore Gisborne, North Island, New Zealand. From September to October 2014, an SSE occurred with a slip over 250 mm and was recorded successfully by the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip deployment II (HOBITSS II). We apply coda wave interferometry on the ambient noise data acquired by nine OBSs deployed by the HOBITSS II to study the seismic velocity variations related to the SSE. The average velocity variations display a decrease on the order of 0.05% during the SSE, followed by an increase of similar magnitude afterwards. Two hypotheses are proposed to explain our observation. The first hypothesis, which has been suggested by previous studies, considers that the velocity decrease during the SSE is caused by more fluids migrating into the upper plate as the SSE breaks a low-permeability seal on the plate boundary. After the SSE, the fluids in the upper plate diffuse gradually and the velocity increases; The second hypothesis is that before the SSE, elastic strain accumulates causing contraction and reduction of porosity and therefore increase of velocity (the velocity increase between SSEs). During the SSE, the velocity decrease is caused by increased porosity as the SSE relieves the accumulated elastic strain on the plate interface, which results in dilation. After the SSE, stress and strain accumulate again, causing a porosity decrease and a velocity increase back to the original value. This study demonstrates that the velocity variations related to SSEs are observable and provides evidence for slow slip mechanism hypotheses.</b></p> <p>The second project (Chapter 4) focuses on the temporal seismic velocity variations associated with an SSE in 2019 offshore Gisborne, North Island, New Zealand. This is a later SSE in the same area as the first project (Chapters 3). Based on the success of the HOBITSS II, more ocean bottom instruments were deployed in the northern Hikurangi subduction zone from 2018 to 2019 (HOBITSS V). An SSE lasting approximately one month from the end of March to the beginning of May 2019 occurred during the deployment and was recorded by the network. The main slip was south of the deployment and the slip beneath the deployment was up to 150 mm. This study applies coda wave interferometry on the ambient noise data acquired by five OBSs and computes seismic velocity variations to investigate their relation to the SSE. A velocity decrease on the order of 0.015% during the SSE and an increase back to the original velocity value are observed at 1–2.5 s. This supports the two hypotheses proposed in Chapters 3: fluid migration and/or stain changes through the SSE cycle. In addition, velocity variations computed from individual stations show velocity increases before the SSE, which are destructively interfered in their average. Such a situation could occur if the SSE migrated across the network. If the velocity increases before the SSE from individual stations are real, they can be only explained by the hypothesis of crustal strain changes (the second hypothesis in project 1). However, fluid migration (the first hypothesis in project 1) may still happen concomitantly.</p> <p>The third project focuses on the tectonics in southwestern Okinawa Trough offshore northeastern Taiwan. The southwestern Okinawa Trough is an active back-arc basin, extending and rifting within the continental lithosphere. The tectonic development of the back-arc basin is still not well-understood. This study uses continuous ambient noise data recorded by 34 OBSs deployed by Academia Sinica at various periods from 2010 to 2018. Cross-correlations on vertical seismic components and pressure gauges are computed to construct Rayleigh/Scholte waves to study the shear wave velocity structure in the southwestern Okinawa Trough. Phase velocities are measured from the Rayleigh/Scholte waves. Shear velocities are inverted from the phase velocities. Results show the velocity in the south of the back-arc rifting axis near the axis is slower than the velocity in the north of the rifting axis, suggesting the velocity structure in the southwestern Okinawa Trough is asymmetric along the rifting axis. Previous studies have shown high heat flows (about 110mW/m 2 on average) in the south of the rifting axis. The low velocity in the south could be caused by the high heat flow that may be related to asymmetric back-arc extension and/or rifting. This study presents the shear wave velocity structure in the southwest Okinawa Trough is asymmetric along the rifting axis, which implies the back-arc extending/rifting is asymmetric in the study region. This study also suggests effective techniques for OBS noise corrections and unwrapping the cycle skipping of phase velocity measurements.</p> <p>In summary, this thesis represents three projects focusing on seismic velocities in two subduction zones using ambient noise data collected by OBSs. The first and second projects study the temporal velocity variations and the relation to SSEs. Both studies observe velocity decreases during the SSEs and increases after the SSEs, supporting two hypotheses of fluid migration and/or stain changes through the SSE cycle. The third project finds the shear velocity structure in the southwestern Okinawa Trough is asymmetric along the rifting center, which may imply the back-arc extension is asymmetric.</p>

Studies of Seismic Velocities in Subduction Zones from Continuous Ocean Bottom Seismometer Data

<p><b>This thesis uses continuous ambient noise data recorded by Ocean Bottom Seismometers (OBSs) to study seismic velocities in the upper crust of the overriding plate. The first and second projects (Chapters 3 and 4) focus on temporal seismic velocity variations in the northern Hikurangi subduction zone offshore the North Island, New Zealand, while the third project (Chapter 5) focuses on shear wave velocities in the southwestern Okinawa Trough offshore northeastern Taiwan. In the first project (Chapters 3), we investigate a region of frequent slow slip events (SSEs) offshore Gisborne, North Island, New Zealand. From September to October 2014, an SSE occurred with a slip over 250 mm and was recorded successfully by the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip deployment II (HOBITSS II). We apply coda wave interferometry on the ambient noise data acquired by nine OBSs deployed by the HOBITSS II to study the seismic velocity variations related to the SSE. The average velocity variations display a decrease on the order of 0.05% during the SSE, followed by an increase of similar magnitude afterwards. Two hypotheses are proposed to explain our observation. The first hypothesis, which has been suggested by previous studies, considers that the velocity decrease during the SSE is caused by more fluids migrating into the upper plate as the SSE breaks a low-permeability seal on the plate boundary. After the SSE, the fluids in the upper plate diffuse gradually and the velocity increases; The second hypothesis is that before the SSE, elastic strain accumulates causing contraction and reduction of porosity and therefore increase of velocity (the velocity increase between SSEs). During the SSE, the velocity decrease is caused by increased porosity as the SSE relieves the accumulated elastic strain on the plate interface, which results in dilation. After the SSE, stress and strain accumulate again, causing a porosity decrease and a velocity increase back to the original value. This study demonstrates that the velocity variations related to SSEs are observable and provides evidence for slow slip mechanism hypotheses.</b></p> <p>The second project (Chapter 4) focuses on the temporal seismic velocity variations associated with an SSE in 2019 offshore Gisborne, North Island, New Zealand. This is a later SSE in the same area as the first project (Chapters 3). Based on the success of the HOBITSS II, more ocean bottom instruments were deployed in the northern Hikurangi subduction zone from 2018 to 2019 (HOBITSS V). An SSE lasting approximately one month from the end of March to the beginning of May 2019 occurred during the deployment and was recorded by the network. The main slip was south of the deployment and the slip beneath the deployment was up to 150 mm. This study applies coda wave interferometry on the ambient noise data acquired by five OBSs and computes seismic velocity variations to investigate their relation to the SSE. A velocity decrease on the order of 0.015% during the SSE and an increase back to the original velocity value are observed at 1–2.5 s. This supports the two hypotheses proposed in Chapters 3: fluid migration and/or stain changes through the SSE cycle. In addition, velocity variations computed from individual stations show velocity increases before the SSE, which are destructively interfered in their average. Such a situation could occur if the SSE migrated across the network. If the velocity increases before the SSE from individual stations are real, they can be only explained by the hypothesis of crustal strain changes (the second hypothesis in project 1). However, fluid migration (the first hypothesis in project 1) may still happen concomitantly.</p> <p>The third project focuses on the tectonics in southwestern Okinawa Trough offshore northeastern Taiwan. The southwestern Okinawa Trough is an active back-arc basin, extending and rifting within the continental lithosphere. The tectonic development of the back-arc basin is still not well-understood. This study uses continuous ambient noise data recorded by 34 OBSs deployed by Academia Sinica at various periods from 2010 to 2018. Cross-correlations on vertical seismic components and pressure gauges are computed to construct Rayleigh/Scholte waves to study the shear wave velocity structure in the southwestern Okinawa Trough. Phase velocities are measured from the Rayleigh/Scholte waves. Shear velocities are inverted from the phase velocities. Results show the velocity in the south of the back-arc rifting axis near the axis is slower than the velocity in the north of the rifting axis, suggesting the velocity structure in the southwestern Okinawa Trough is asymmetric along the rifting axis. Previous studies have shown high heat flows (about 110mW/m 2 on average) in the south of the rifting axis. The low velocity in the south could be caused by the high heat flow that may be related to asymmetric back-arc extension and/or rifting. This study presents the shear wave velocity structure in the southwest Okinawa Trough is asymmetric along the rifting axis, which implies the back-arc extending/rifting is asymmetric in the study region. This study also suggests effective techniques for OBS noise corrections and unwrapping the cycle skipping of phase velocity measurements.</p> <p>In summary, this thesis represents three projects focusing on seismic velocities in two subduction zones using ambient noise data collected by OBSs. The first and second projects study the temporal velocity variations and the relation to SSEs. Both studies observe velocity decreases during the SSEs and increases after the SSEs, supporting two hypotheses of fluid migration and/or stain changes through the SSE cycle. The third project finds the shear velocity structure in the southwestern Okinawa Trough is asymmetric along the rifting center, which may imply the back-arc extension is asymmetric.</p>

Spatiotemporal slip distributions associated with the 2018–2019 Bungo Channel long-term slow slip event inverted from GNSS data

Long-term slow slip events (L-SSEs) have repeatedly occurred beneath the Bungo Channel in southwestern Japan with durations of several months to a couple of years, with a recurrence interval of approximately 6 years. We estimated the spatiotemporal slip distributions of the 2018–2019 Bungo Channel L-SSE by inverting processed GNSS time series data. This event was divided into two subevents, with the first on the southwest side of the Bungo Channel from 2018.3 to 2018.7 and the second beneath the Bungo Channel from 2018.8 to 2019.4. Tectonic tremors became active on the downdip side of the L-SSE occurrence region when large slow slips took place beneath the Bungo Channel. Compared with the previous Bungo Channel L-SSEs, this spatiotemporal slip pattern and amount were similar to those of the 2002–2004 L-SSE. However, the slip expanded in the northeast and southwest directions in the latter half of the second subevent. The maximum amount of slip, the maximum slip velocity, the total released seismic moment, and the moment magnitude of the 2018–2019 L-SSE were estimated to be 28 cm, 54 cm/year, $$4.4 \times 10^{19}$$ 4.4 × 10 19 Nm, and 7.0, respectively, all of which were the largest among the 1996–1998, 2002–2004, 2009–2011, and 2018–2019 L-SSEs.

Studies of Seismic Velocities in Subduction Zones from Continuous Ocean Bottom Seismometer Data

<p><b>This thesis uses continuous ambient noise data recorded by Ocean Bottom Seismometers (OBSs) to study seismic velocities in the upper crust of the overriding plate. The first and second projects (Chapters 3 and 4) focus on temporal seismic velocity variations in the northern Hikurangi subduction zone offshore the North Island, New Zealand, while the third project (Chapter 5) focuses on shear wave velocities in the southwestern Okinawa Trough offshore northeastern Taiwan. In the first project (Chapters 3), we investigate a region of frequent slow slip events (SSEs) offshore Gisborne, North Island, New Zealand. From September to October 2014, an SSE occurred with a slip over 250 mm and was recorded successfully by the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip deployment II (HOBITSS II). We apply coda wave interferometry on the ambient noise data acquired by nine OBSs deployed by the HOBITSS II to study the seismic velocity variations related to the SSE. The average velocity variations display a decrease on the order of 0.05% during the SSE, followed by an increase of similar magnitude afterwards. Two hypotheses are proposed to explain our observation. The first hypothesis, which has been suggested by previous studies, considers that the velocity decrease during the SSE is caused by more fluids migrating into the upper plate as the SSE breaks a low-permeability seal on the plate boundary. After the SSE, the fluids in the upper plate diffuse gradually and the velocity increases; The second hypothesis is that before the SSE, elastic strain accumulates causing contraction and reduction of porosity and therefore increase of velocity (the velocity increase between SSEs). During the SSE, the velocity decrease is caused by increased porosity as the SSE relieves the accumulated elastic strain on the plate interface, which results in dilation. After the SSE, stress and strain accumulate again, causing a porosity decrease and a velocity increase back to the original value. This study demonstrates that the velocity variations related to SSEs are observable and provides evidence for slow slip mechanism hypotheses.</b></p> <p>The second project (Chapter 4) focuses on the temporal seismic velocity variations associated with an SSE in 2019 offshore Gisborne, North Island, New Zealand. This is a later SSE in the same area as the first project (Chapters 3). Based on the success of the HOBITSS II, more ocean bottom instruments were deployed in the northern Hikurangi subduction zone from 2018 to 2019 (HOBITSS V). An SSE lasting approximately one month from the end of March to the beginning of May 2019 occurred during the deployment and was recorded by the network. The main slip was south of the deployment and the slip beneath the deployment was up to 150 mm. This study applies coda wave interferometry on the ambient noise data acquired by five OBSs and computes seismic velocity variations to investigate their relation to the SSE. A velocity decrease on the order of 0.015% during the SSE and an increase back to the original velocity value are observed at 1–2.5 s. This supports the two hypotheses proposed in Chapters 3: fluid migration and/or stain changes through the SSE cycle. In addition, velocity variations computed from individual stations show velocity increases before the SSE, which are destructively interfered in their average. Such a situation could occur if the SSE migrated across the network. If the velocity increases before the SSE from individual stations are real, they can be only explained by the hypothesis of crustal strain changes (the second hypothesis in project 1). However, fluid migration (the first hypothesis in project 1) may still happen concomitantly.</p> <p>The third project focuses on the tectonics in southwestern Okinawa Trough offshore northeastern Taiwan. The southwestern Okinawa Trough is an active back-arc basin, extending and rifting within the continental lithosphere. The tectonic development of the back-arc basin is still not well-understood. This study uses continuous ambient noise data recorded by 34 OBSs deployed by Academia Sinica at various periods from 2010 to 2018. Cross-correlations on vertical seismic components and pressure gauges are computed to construct Rayleigh/Scholte waves to study the shear wave velocity structure in the southwestern Okinawa Trough. Phase velocities are measured from the Rayleigh/Scholte waves. Shear velocities are inverted from the phase velocities. Results show the velocity in the south of the back-arc rifting axis near the axis is slower than the velocity in the north of the rifting axis, suggesting the velocity structure in the southwestern Okinawa Trough is asymmetric along the rifting axis. Previous studies have shown high heat flows (about 110mW/m 2 on average) in the south of the rifting axis. The low velocity in the south could be caused by the high heat flow that may be related to asymmetric back-arc extension and/or rifting. This study presents the shear wave velocity structure in the southwest Okinawa Trough is asymmetric along the rifting axis, which implies the back-arc extending/rifting is asymmetric in the study region. This study also suggests effective techniques for OBS noise corrections and unwrapping the cycle skipping of phase velocity measurements.</p> <p>In summary, this thesis represents three projects focusing on seismic velocities in two subduction zones using ambient noise data collected by OBSs. The first and second projects study the temporal velocity variations and the relation to SSEs. Both studies observe velocity decreases during the SSEs and increases after the SSEs, supporting two hypotheses of fluid migration and/or stain changes through the SSE cycle. The third project finds the shear velocity structure in the southwestern Okinawa Trough is asymmetric along the rifting center, which may imply the back-arc extension is asymmetric.</p>

Automatic Detection of Slow Slip Events Using the PICCA: Application to Chilean GNSS Data

The detection of transient events related to slow earthquakes in GNSS positional time series is key to understanding seismogenic processes in subduction zones. Here, we present a novel Principal and Independent Components Correlation Analysis (PICCA) method that allows for the temporal and spatial detection of transient signals. The PICCA is based on an optimal combination of the principal (PCA) and independent component analysis (ICA) of positional time series of a GNSS network. We assume that the transient signal is mostly contained in one of the principal or independent components. To detect the transient, we applied a method where correlations between sliding windows of each PCA/ICA component and each time series are calculated, obtaining the stations affected by the slow slip event and the onset time from the resulting correlation peaks. We first tested and calibrated the method using synthetic signals from slow earthquakes of different magnitudes and durations and modelled their effect in the network of GNSS stations in Chile. Then, we analyzed three transient events related to slow earthquakes recorded in Chile, in the areas of Iquique, Copiapó, and Valparaíso. For synthetic data, a 150 days event was detected using the PCA-based method, while a 3 days event was detected using the ICA-based method. For the real data, a long-term transient was detected by PCA, while a 16 days transient was detected by ICA. It is concluded that simultaneous use of both signal separation methods (PICCA) is more effective when searching for transient events. The PCA method is more useful for long-term events, while the ICA method is better suited to recognize events of short duration. PICCA is a promising tool to detect transients of different characteristics in GNSS time series, which will be used in a next stage to generate a catalog of SSEs in Chile.

Imaging evolution of Cascadia slow-slip event using high-rate GPS

The slip history of short-term slow slip event (SSE) is typically inferred from daily Global Positioning System (GPS) data, which, however, cannot image the sub-daily processes, leaving the underlying mechanisms of SSEs elusive. To address the temporal resolution issue, we attempted to employ the kinematic subdaily GPS analysis, which has never been applied to SSE studies because its signal-to-noise ratio has been believed too low. By carefully post-processing sub-daily positions to remove non-tectonic position fluctuation, our 30-minute kinematic data clearly exhibits the transient motion of a few mm during one Cascadia SSE. A spatiotemporal slip image by inverting the 30-minute data exhibits a multi-stage evolution; it consists of an isotropic growth of SSE followed by an along-strike migration and termination within the rheologically controlled down-dip width. This transition at the slip growth mode is similar to the rupture growth of regular earthquakes, implying the presence of common mechanical factors behind the two distinct slip phenomena. The comparison with a slip inversion of the daily GPS demonstrates the current performance and limitation of the subdaily data in the SSE detection and imaging.Better understanding of the non-tectonic noise in the kinematic GPS analysis will further improve the temporal resolution of SSE.