@article {iwasaki_continuous_2022,
	title = {Continuous Tremor Activity With Stable Polarization Direction Following the 2014 Large Slow Slip Event in the Hikurangi Subduction Margin Offshore New Zealand},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {127},
	number = {2},
	year = {2022},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2021JB022161},
	month = {11/2021},
	pages = {e2021JB022161},
	abstract = {Many types of slow earthquakes have been discovered at subduction zones around the world. However, the physical process of these slow earthquakes is not well understood. To monitor offshore slow earthquakes, a marine seismic and geodetic experiment was conducted at the Hikurangi subduction margin from May 2014 to June 2015. During this experiment, a large slow slip event (Mw 6.8) occurred directly beneath the ocean bottom seismometer (OBS) network. In this study, S-wave splitting and polarization analysis methods, which have been previously used on onshore data to investigate tremor and anisotropy, are applied to continuous OBS waveform data to identify tremors that are too small to detect by the envelope cross correlation method. Continuous tremor activity with stable polarization directions is detected at the end of the 2014 slow slip event and continued for about 2 weeks. The tremors are generated around a southwest bend in the slow slip contours and at the landward edge of a subducted seamount. Our findings corroborate a previous interpretation, based on burst-type repeating earthquakes and intermittent tremor, that localized slow slip and tremor around the seamount was triggered by fluid migration following the large plate boundary slow slip event and indicate tremor occurred continuously rather than as isolated and sporadic individual events.},
	keywords = {New Zealand, polarization, S-wave splitting, Seamount, slow slip, tremor},
	issn = {2169-9356},
	doi = {10.1029/2021JB022161},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2021JB022161},
	author = {Iwasaki, Yuriko and Mochizuki, Kimihiro and Ishise, Motoko and Todd, Erin K. and Schwartz, Susan Y. and Zal, Hubert and Savage, Martha K. and Henrys, Stuart and Sheehan, Anne F. and Ito, Yoshihiro and Wallace, Laura M. and Webb, Spahr C. and Yamada, Tomoaki and Shinohara, Masanao}
}
@article {todd_earthquakes_2018,
	title = {Earthquakes and Tremor Linked to Seamount Subduction During Shallow Slow Slip at the Hikurangi Margin, New Zealand},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {123},
	number = {8},
	year = {2018},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2018JB016136},
	pages = {6769{\textendash}6783},
	abstract = {Shallow slow slip events have been well documented offshore Gisborne at the northern Hikurangi subduction margin, New Zealand, and are associated with tectonic tremor downdip of the slow slip patch and increases in local microseismicity. Tremor and seismicity on the shallow subduction interface are often poorly resolved due to their distance from land-based seismic and geodetic networks. To address this shortcoming, the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip experiment deployed 24 absolute pressure gauges and 15 ocean bottom seismometers on the seafloor above the Gisborne slow slip patch to investigate the spatial and temporal extent of slow slip and associated tremor and earthquake activity. We present a detailed spatiotemporal analysis of the seismic signatures of various interplate slip processes associated with the September/October 2014 Gisborne slow slip event. Tectonic tremor begins toward the end and continues after the geodetically constrained slow slip event and is localized in the vicinity of two subducted seamounts within and updip of the slow slip patch. The subsequent, rather than synchronous occurrence of tremor suggests that tremor may be triggered by stress changes induced by slow slip. However, Coulomb failure stress change models based on the slow slip distribution fail to predict the location of tremor, suggesting that seamount subduction plays a dominant role in the stress state of the shallow megathrust. This and the observed interplay of seismic and aseismic interplate slip processes imply that stress changes from slow slip play a secondary role in the distribution of associated microseismicity.},
	keywords = {earthquake, New Zealand, Seamount, slow slip, subduction, tremor},
	issn = {2169-9356},
	doi = {10.1029/2018JB016136},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018JB016136},
	author = {Todd, Erin K. and Schwartz, Susan Y. and Mochizuki, Kimihiro and Wallace, Laura M. and Sheehan, Anne F. and Webb, Spahr C. and Williams, Charles A. and Nakai, Jenny and Yarce, Jefferson and Fry, Bill and Henrys, Stuart and Ito, Yoshihiro}
}
@article {ball_joint_2014,
	title = {A joint Monte Carlo analysis of seafloor compliance, Rayleigh wave dispersion and receiver functions at ocean bottom seismic stations offshore New Zealand},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {15},
	number = {12},
	year = {2014},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/2014GC005412},
	pages = {5051{\textendash}5068},
	abstract = {Teleseismic body-wave imaging techniques such as receiver function analysis can be notoriously difficult to employ on ocean-bottom seismic data due largely to multiple reverberations within the water and low-velocity sediments. In lieu of suppressing this coherently scattered noise in ocean-bottom receiver functions, these site effects can be modeled in conjunction with shear velocity information from seafloor compliance and surface wave dispersion measurements to discern crustal structure. A novel technique to estimate 1-D crustal shear-velocity profiles from these data using Monte Carlo sampling is presented here. We find that seafloor compliance inversions and P-S conversions observed in the receiver functions provide complimentary constraints on sediment velocity and thickness. Incoherent noise in receiver functions from the MOANA ocean bottom seismic experiment limit the accuracy of the practical analysis at crustal scales, but synthetic recovery tests and comparison with independent unconstrained nonlinear optimization results affirm the utility of this technique in principle.},
	keywords = {Monte Carlo, New Zealand, OBS receiver function, seafloor compliance, sediment shear modes, surface waves},
	issn = {1525-2027},
	doi = {10.1002/2014GC005412},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2014GC005412},
	author = {Ball, Justin S. and Sheehan, Anne F. and Stachnik, Joshua C. and Lin, Fan-Chi and Collins, John A.}
}
@article {zietlow_upper_2014,
	title = {Upper mantle seismic anisotropy at a strike-slip boundary: South Island, New Zealand},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {119},
	number = {2},
	year = {2014},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/2013JB010676},
	pages = {1020{\textendash}1040},
	abstract = {New shear wave splitting measurements made from stations onshore and offshore the South Island of New Zealand show a zone of anisotropy 100{\textendash}200 km wide. Measurements in central South Island and up to approximately 100 km offshore from the west coast yield orientations of the fast quasi-shear wave nearly parallel to relative plate motion, with increased obliquity to this orientation observed farther from shore. On the eastern side of the island, fast orientations rotate counterclockwise to become nearly perpendicular to the orientation of relative plate motion approximately 200 km off the east coast. Uniform delay times between the fast and slow quasi-shear waves of nearly 2.0 s onshore continue to stations approximately 100 km off the west coast, after which they decrease to 1 s at 200 km. Stations more than 300 km from the west coast show little to no splitting. East coast stations have delay times around 1 s. Simple strain fields calculated from a thin viscous sheet model (representing distributed lithospheric deformation) with strain rates decreasing exponentially to both the northwest and southeast with e-folding dimensions of 25{\textendash}35 km (approximately 75\% of the deformation within a zone 100{\textendash}140 km wide) match orientations and amounts of observed splitting. A model of deformation localized in the lithosphere and then spreading out in the asthenosphere also yields predictions consistent with observed splitting if, at depths of 100{\textendash}130 km below the lithosphere, typical grain sizes are 6{\textendash}7 mm.},
	keywords = {mantle lithosphere, MOANA, New Zealand, ocean bottom seismometers, seismic anisotropy, South Island},
	issn = {2169-9356},
	doi = {10.1002/2013JB010676},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2013JB010676},
	author = {Zietlow, Daniel W. and Sheehan, Anne F. and Molnar, Peter H. and Savage, Martha K. and Hirth, Greg and Collins, John A. and Hager, Bradford H.}
}
@article {yang_character_2012,
	title = {The character of seafloor ambient noise recorded offshore New Zealand: Results from the MOANA ocean bottom seismic experiment},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {13},
	number = {10},
	year = {2012},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2012GC004201},
	abstract = {{We analyze the characteristics of ambient noise recorded on ocean-bottom seismographs using data from the 2009{\textendash}2010 MOANA (Marine Observations of Anisotropy Near Aotearoa) seismic experiment deployed west and east of South Island, New Zealand. Microseism and infragravity noise peaks are clear on data recorded on the vertical channel of the seismometer and on the pressure sensor. The noise levels in the infragravity band (<0.03 Hz) on the horizontal seismometer channels are too high to show the infragravity peak. There is a small difference (\~{}0.25 Hz versus \~{}0.2 Hz) in microseism peak frequencies between the two sides of the South Island on all three seismic channels. Our results show clear depth dependence between the peak frequency of infragravity waves and the water depth. We find that the product of water depth and wave number at the peak frequency is a constant},
	keywords = {continental shelf, infragravity wave, New Zealand, noise, OBS},
	issn = {1525-2027},
	doi = {10.1029/2012GC004201},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2012GC004201},
	author = {Yang, Zhaohui and Sheehan, Anne F. and Collins, John A. and Laske, Gabi}
}
