Observations of Earth’s Normal Modes on Broadband Ocean Bottom Seismometers

It is generally thought that high noise levels in the oceans inhibit the observation of long-period earthquake signals such as Earth’s normal modes on ocean bottom seismometers (OBSs). Here, we document the observation of Earth’s gravest modes at periods longer than 500 s (or frequencies below 2 mHz). We start with our own 2005–2007 Plume-Lithosphere-Undersea-Mantle Experiment (PLUME) near Hawaii that deployed a large number of broadband OBSs for the first time. We collected high-quality normal mode spectra for the great November 15, 2006 Kuril Islands earthquake on multiple OBSs. The random deployment of instruments from different OBS groups allows a direct comparison between different broadband seismometers. For this event, mode S06 (1.038 mHz) consistently rises above the background noise at all OBSs that had a Nanometrics Trillium T-240 broadband seismometer. We also report observations of other deployments in the Pacific ocean that involved instruments of the U.S. OBS Instrument Pool (OBSIP) where we observe even mode S04 (0.647 mHz). Earth’s normal modes were never the initial target of any OBS deployment, nor was any other ultra-low-frequency signal. However, given the high costs of an OBS campaign, the fact that data are openly available to future investigators not involved in the campaign, and the fact that seismology is evolving to investigate ever-new signals, this paper makes the case that the investment in a high-quality seismic sensor may be a wise one, even for a free-fall OBS

Finite-frequency effects in global surface-wave tomography

We compare traditional ray-theoretical surface-wave tomography with finite-frequency tomography, using 3-D Born sensitivity kernels for long-period, fundamental-mode dispersion measurements. The 3-D kernels preserve sidelobes beyond the first Fresnel zone, and fully account for the directional dependence of surface-wave scattering, and the effects of time-domain tapering and seismic source radiation. Tomographic inversions of Love and Rayleigh phase-delay measurements and synthetic checkerboard tests show that (1) small-scale S-wave velocity anomalies are better resolved using finite-frequency sensitivity kernels, especially in the lowermost upper mantle; (2) the resulting upper-mantle heterogeneities are generally stronger in amplitude than those recovered using ray theory and (3) finite-frequency tomographic models fit long-period dispersion data better than ray-theoretical models of comparable roughness. We also examine the reliability of 2-D, phase-velocity sensitivity kernels in global surface-wave tomography, and show that phase-velocity kernels based upon a forward-scattering approximation or previously adopted geometrical simplifications do not reliably account for finite-frequency wave-propagation effects. 3-D sensitivity kernels with full consideration of directional-dependent seismic scattering are the preferred method of inverting long-period dispersion data. Finally, we derive 2-D boundary sensitivity kernels for lateral variations in crustal thickness, and show that finite-frequency crustal effects are not negligible in long-period surface-wave dispersion studies, especially for paths along continent-ocean boundaries. Unfortunately, we also show that, in global studies, linear perturbation theory is not sufficiently accurate to make reliable crustal corrections, due to the large difference in thickness between oceanic and continental crust

Asymmetric shallow mantle structure beneath the Hawaiian Swell—evidence from Rayleigh waves recorded by the PLUME network

Author Posting. © The Author(s), 2011. This article is posted here by permission of John Wiley & Sons for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 187 (2011): 1725–1742, doi:10.1111/j.1365-246X.2011.05238.x.We present models of the 3-D shear velocity structure of the lithosphere and asthenosphere beneath the Hawaiian hotspot and surrounding region. The models are derived from long-period Rayleigh-wave phase velocities that were obtained from the analysis of seismic recordings collected during two year-long deployments for the Hawaiian Plume-Lithosphere Undersea Mantle Experiment. For this experiment, broad-band seismic sensors were deployed at nearly 70 seafloor sites as well as 10 sites on the Hawaiian Islands. Our seismic images result from a two-step inversion of path-averaged dispersion curves using the two-station method. The images reveal an asymmetry in shear velocity structure with respect to the island chain, most notably in the lower lithosphere at depths of 60 km and greater, and in the asthenosphere. An elongated, 100-km-wide and 300-km-long low-velocity anomaly reaches to depths of at least 140 km. At depths of 60 km and shallower, the lowest velocities are found near the northern end of the island of Hawaii. No major velocity anomalies are found to the south or southeast of Hawaii, at any depth. The low-velocity anomaly in the asthenosphere is consistent with an excess temperature of 200–250 °C and partial melt at the level of a few percent by volume, if we assume that compositional variations as a result of melt extraction play a minor role. We also image small-scale low-velocity anomalies within the lithosphere that may be associated with the volcanic fields surrounding the Hawaiian Islands.This research was financed by the National Science Foundation under grants OCE-00-02470 and OCE-00-02819. Markee was partly sponsored by a SIO graduate student fellowship

Shear wave splitting at the Hawaiian hot spot from the PLUME land and ocean bottom seismometer deployments

Author Posting. © American Geophysical Union, 2012. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 13 (2012): Q02007, doi:10.1029/2011GC003881.We examine upper mantle anisotropy across the Hawaiian Swell by analyzing shear wave splitting of teleseismic SKS waves recorded by the PLUME broadband land and ocean bottom seismometer deployments. Mantle anisotropy beneath the oceans is often attributed to flow-induced lattice-preferred orientation of olivine. Splitting observations may reflect a combination of both fossil lithospheric anisotropy and anisotropy due to present-day asthenospheric flow, and here we address the question whether splitting provides diagnostic information on possible asthenospheric plume flow at Hawaii. We find that the splitting fast directions are coherent and predominantly parallel to the fossil spreading direction, suggesting that shear wave splitting dominantly reflects fossil lithospheric anisotropy. The signature of anisotropy from asthenospheric flow is more subtle, although it could add some perturbation to lithospheric splitting. The measured delay times are typically 1 s or less, although a few stations display larger splitting delays of 1–2 s. The variability in the delay times across the different stations indicates differences in the degree of anisotropy or in the thickness of the anisotropic layer or in the effect of multilayer anisotropy. Regions with smaller splitting times may have experienced processes that modified the lithosphere and partially erased the fossil anisotropy; alternatively, asthenospheric splitting may either constructively add to or destructively subtract from lithospheric splitting to produce the observed variability in delay times.The PLUME project was supported by NSF.2012-08-1

The character of seafloor ambient noise recorded offshore New Zealand : results from the MOANA ocean bottom seismic experiment

Author Posting. © American Geophysical Union, 2012. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 13 (2012): Q10011, doi:10.1029/2012GC004201.We analyze the characteristics of ambient noise recorded on ocean-bottom seismographs using data from the 2009–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, koH = 1.5. This relationship can be used to determine the variation of phase and group velocity of infragravity waves with water depth, and the location of the infragravity peak and corresponding noise notch at any water depth. These estimates of spectral characteristics, particularly low noise bands, are useful for future OBS deployments.The collection of OBS data was supported by the National Science Foundation Continental Dynamics program under grants EAR-0409564, EAR-0409609, and EAR-0409835.2013-04-1

Underplating of the Hawaiian Swell : evidence from teleseismic receiver functions

Author Posting. © The Authors, 2010. This article is posted here by permission of John Wiley & Sons for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 183 (2010): 313-329, doi:10.1111/j.1365-246X.2010.04720.x.The Hawaiian Islands are the canonical example of an age-progressive island chain, formed by volcanism long thought to be fed from a hotspot source that is more or less fixed in the mantle. Geophysical data, however, have so far yielded contradictory evidence on subsurface structure. The substantial bathymetric swell is supportive of an anomalously hot upper mantle, yet seafloor heat flow in the region does not appear to be enhanced. The accumulation of magma beneath pre-existing crust (magmatic underplating) has been suggested to add chemical buoyancy to the swell, but to date the presence of underplating has been constrained only by local active-source experiments. In this study, teleseismic receiver functions derived from seismic events recorded during the PLUME project were analysed to obtain a regional map of crustal structure for the Hawaiian Swell. This method yields results that compare favourably with those from previous studies, but permits a much broader view than possible with active-source seismic experiments. Our results indicate that the crustal structure of the Hawaiian Islands is quite complicated and does not conform to the standard model of sills fed from a central source. We find that a shallow P-to-s conversion, previously hypothesized to result from the volcano-sediment interface, corresponds more closely to the boundary between subaerial and subaqueous extrusive material. Correlation between uplifted bathymetry at ocean-bottom-seismometer locations and presence of underplating suggests that much of the Hawaiian Swell is underplated, whereas a lack of underplating beneath the moat surrounding the island of Hawaii suggests that underplated crust outward of the moat has been fed from below by dykes through the lithosphere rather than by sills spreading from the island centre. Local differences in underplating may reflect focusing of magma-filled dykes in response to stress from volcanic loading. Finally, widespread underplating adds chemical buoyancy to the swell, reducing the amplitude of a mantle thermal anomaly needed to match bathymetry and supporting observations of normal heat flow.We are grateful to the Ocean Sciences Division of the U.S. National Science Foundation for their support of this project under grants OCE-0002470, OCE-0002552 and OCE-0002819

Glaciohydraulic seismic tremors on an Alpine glacier

Hydraulic processes impact viscous and brittle ice deformation. Water-driven fracturing as well as turbulent waterflow within and beneath glaciers radiate seismic waves which provide insights into otherwise hard-to-access englacial and subglacial environments. In this study, we analyze glaciohydraulic tremors recorded by four seismic arrays installed in different parts of Glacier de la Plaine Morte, Switzerland. Data were recorded during the 2016 melt season including the sudden subglacial drainage of an ice-marginal lake. Together with our seismic data, discharge, lake level, and ice flow measurements provide constraints on glacier hydraulics. We find that the tremors are generated by subglacial water flow, in moulins, and byicequake bursts. The dominating process can vary on sub-kilometer and sub-daily scales. Consistent with field observations,continuous source tracking via matched-field processing suggests a gradual upglacier progression of an efficient drainage sys-tem as the melt season progresses. The ice-marginal lake likely connects to this drainage system via hydrofracturing, whichis indicated by sustained icequake signals emitted from the proximity of the lake basin and starting roughly 24 hours prior to the lake drainage. To estimate the hydraulics associated with the drainage, we use tremor-discharge scaling relationships. Our analysis suggests a pressurization of the subglacial environment at the drainage onset, followed by an increase in the hydraulicradii of the conduits and a subsequent decrease in the subglacial water pressure as the capacity of the drainage system increases.The pressurization is in phase with the drop in the lake level and its retrieved maximum coincides with ice-uplift measured via GPS. Our results highlight the use of cryo-seismology for monitoring glacier hydraulicsISSN:1994-0432ISSN:1994-044

Tilt corrections for normal mode observations on ocean bottom seismic data, an example from the PI-LAB experiment

Earth's normal modes are fundamental observations used in global seismic tomography to understand Earth structure. Land seismic station coverage is sufficient to constrain the broadest scale Earth structures. However, 70% of Earth's surface is covered by the oceans, hampering our ability to observe variations in local mode frequencies that contribute to imaging small-scale structures. Broadband ocean bottom seismometers can record spheroidal modes to fill in gaps in global data coverage. Ocean bottom recordings are contaminated by signals from complex interactions between ocean and solid Earth dynamics at normal mode frequencies. We present a method for correcting tilt on broadband ocean bottom seismometers by rotation. The correction improves the ability of some instruments to observe spheroidal modes down to 0S4. We demonstrate this method using 15 broadband ocean bottom seismometers from the PI-LAB array. We measure normal mode peak frequency shifts and compare with 1-D reference mode frequencies and predictions from 3-D global models. Our measurements agree with the 3-D models for modes between 0S14 - 0S37 with small but significant differences. These differences likely reflect real Earth structure. This suggests incorporating ocean bottom normal mode measurements into global inversions will improve models of global seismic velocity structure