Ocean Bottom Seismometers Research Articles

Seismic anisotropy indicates organized melt beneath the Mid-Atlantic Ridge aids seafloor spreading

Abstract Lithospheric plates diverge at mid-ocean ridges and asthenospheric mantle material rises in response. The rising material decompresses, which can result in partial melting, potentially impacting the driving forces of the system. Yet the geometry and spatial distribution of the melt as it migrates to the ridge axis are debated. Organized melt fabrics can cause strong seismic anisotropy, which can be diagnostic of melt, although this is typically not found at ridges. We present anisotropic constraints from an array of 39 ocean-bottom seismometers deployed on 0–80 Ma lithosphere from March 2016 to March 2017 near the equatorial Mid-Atlantic Ridge (MAR). Local and SKS measurements show anisotropic fast directions away from the ridge axis, which are consistent with strain and associated fabric caused by plate motions with short delay times, δt (<1.1 s). Near the ridge axis, we find several ridge-parallel fast splitting directions, φ, with SKS δt that are much longer (1.7–3.8 s). This is best explained by ridge-parallel sub-vertical orientations of sheet-like melt pockets. This observation is much different than anisotropic patterns observed at other ridges, which typically reflect fabric related to plate motions. One possibility is that thicker sub-ridge lithosphere with steep sub-ridge topography beneath slower spreading centers focuses melt into vertical, ridge-parallel melt bands, which effectively weakens the plate. Associated buoyancy forces elevate the sub-ridge plate, providing greater potential energy and enhancing the driving forces of the plates.

Optimizing Ocean Bottom Seismometers for Shallow Water Using Shielding and Horizontal Pressure Gradient Data

AbstractWe show large “shields” can significantly reduce horizontal component noise levels for ocean bottom seismometers by shielding the seismic sensors from shaking due to ocean currents. However, at short period (<1 s) some noise from vibration or rocking of the shield under current can partly couple through the ground causing peaks in the noise spectra. Deformation of the seafloor under ocean wave loading (“compliance”) raises noise levels at long period (>5 s) for both vertical and horizontal components at shallow water sites. Pressure gauge data is now routinely used to predict and remove this source of noise on vertical component data. We have developed horizontal pressure gradient (HPG) gauges to measure and remove this source of noise from horizontal component data. The combination of HPG data and shielding enables low noise long period seismic observations at otherwise very noisy shallow seafloor sites.

Discontinuous Igneous Addition Along the Eastern North American Margin Beneath the East Coast Magnetic Anomaly

AbstractDetailed models of crustal structure at volcanic passive margins offer insight into the role of magmatism and the distribution of igneous addition during continental rifting. The Eastern North American Margin (ENAM) is a volcanic passive margin that formed during the breakup of Pangea ∼200 Myr ago. The offshore, margin‐parallel East Coast Magnetic Anomaly (ECMA) is thought to mark the locus of syn‐rift magmatism. Previous widely spaced margin‐perpendicular studies seismically imaged igneous addition as seaward dipping reflectors (SDRs) and high velocity lower crust (HVLC; >7.2 km/s) beneath the ECMA. Along‐strike imaging is necessary to more accurately determine the distribution and volume of igneous addition during continental breakup. We use wide‐angle, marine active‐source seismic data from the 2014–2015 ENAM Community Seismic Experiment to determine crustal structure beneath a ∼370‐km‐long section of the ECMA. P‐wave velocity models based on data from short‐period ocean bottom seismometers reveal a ∼21‐km‐thick crust with laterally variable lower crust velocities ranging from 6.9 to 7.5 km/s. Sections with HVLC (>7.2 km/s) alternate with two ∼30‐km‐wide areas where the average velocities are less than 7.0 km/s. This variable structure indicates that HVLC is discontinuous along the margin, reflecting variable amounts of intrusion along‐strike. Our results suggest that magmatism during rifting was segmented. The HVLC discontinuities roughly align with locations of Mid‐Atlantic Ridge fracture zones, which may suggest that rift segmentation influenced later segmentation of the Mid‐Atlantic Ridge.

Crustal structure beneath the Alcock Seamount in the Andaman Sea

Alcock and Sewell seamounts are important geomorphic features in the Andaman Sea due to their association with the tectonic history of the Sunda subduction zone. Despite the availability of global geophysical data, the nature of crust beneath these seamounts remains debated due to the paucity of local geophysical data. Therefore, a passive ocean bottom seismometer (OBS) experiment was conducted by deploying broadband OBSs in the Andaman region from 25th December 2013 to 10th May 2014. Based on the local earthquakes recorded by OBS, we attempted to analyse the waveforms of the events using the state-of-the-art transfer function, and estimated the crustal velocities beneath the Alcock seamount. The inferred P wave velocity increases from 5.7 km/s to 7.3 km/s at a depth of 9 km and from 7.3 km/s to 8.1 km/s at a depth of 15 km. We validated the upper velocity model using an active OBS record obtained from a concurrent multichannel seismic survey acquired onboard RV Samudra Ratnakar. Further, we carried out traveltime streamer tomography by downward continuing the shot gathers of MCS data and obtaining high-resolution velocity-depth models over two shallow ridges on the Alcock Seamount. Tomography-derived velocity-depth model matches with that estimated from the active and passive OBS records. The presence of high subsurface velocities beneath the Alcock Seamount suggests the occurrence of extensive volcanoclastic deposits, possibly formed as a result of the magmatic intrusion. Higher velocities for deep crustal layers suggest an oceanic or island arc crust beneath the Alcock Seamount, which has probably evolved due to backarc volcanism in the Andaman Sea.

A new method of alternating shooting of two different seismic sources for deep geophysical surveys

In the latest geophysical survey crossing the Ninety East Ridge of the Indian Ocean, a new method was employed to perform proportional double seismic source excitation and synchronously receive signals from the sea surface and the seabed. The two seismic sources used for excitation were two sets of gun arrays with different energies and dominant frequencies, a G gun array and a Bolt gun array. The G gun array consisted of 3 G.II guns with a total capacity of 450 in3 and a dominant frequency of 20–100 Hz. The Bolt gun array consisted of 4 Bolt 1500LL air guns with a total capacity of 6000 in3 and a dominant frequency of 10–40 Hz. The seismic receiving system comprised a 480-channel seismic streamer towed from the sea surface and 21 ocean bottom seismometers (OBS). During offshore operations, the integrated navigation system produced equidistant trigger signals at an interval of 50 m. The trigger signals were distributed to the G gun array and Bolt gun array at a ratio of 3:1 after passing through a pulse signal proportional distributor. The two sets of gun arrays fired alternatingly at a given ratio. The receiving equipment on the sea surface and seabed simultaneously received the seismic signals excited by the two sets of gun arrays. After targeted data processing, in addition to the seismic profile generated by the conventional G gun seismic source, the deep seismic profile generated by the Bolt gun seismic source and the survey profile of the active-source OBS were obtained simultaneously. The penetration depths of the three sets of profiles reach 2 km, 6 km, and 30 km, respectively, greatly improving the efficiency of offshore deep-sea seismic surveys.

Mantle flow underneath the South China Sea revealed by seismic anisotropy.

It has long been established that plastic flow in the asthenosphere interacts constantly with the overlying lithosphere and plays a pivotal role in controlling the occurrence of geohazards such as earthquakes and volcanic eruptions. Unfortunately, accurately characterizing the direction and lateral extents of the mantle flow field is notoriously difficult, especially in oceanic areas where deployment of ocean bottom seismometers (OBSs) is expensive and thus rare. In this study, by applying shear wave splitting analyses to a dataset recorded by an OBS array that we deployed between mid-2019 and mid-2020 in the South China Sea (SCS), we show that the dominant mantle flow field has a NNW-SSE orientation, which can be attributed to mantle flow extruded from the Tibetan Plateau by the ongoing Indian-Eurasian collision. In addition, the results suggest that E-W oriented flow fields observed in South China and the Indochina Peninsula do not extend to the central SCS.

Analysis of clock drift based on a three‐dimensional controlled‐source ocean bottom seismometer experiment

AbstractClock drift is common for the autonomous ocean bottom seismometer, and a technique of linear drift correction is usually used to fix this kind of timing error assuming that the internal clock has constant drift rate during deployment. However, non‐linear component of clock drift has also been found recently. In order to test the behaviours of the clock drift for 19 Chinese ocean bottom seismometers from a three‐dimensional controlled‐source seismic survey, we have analysed two types of differences in direct wave travel‐times, including one from all shots to the ocean bottom seismometer onboard and another one from the nearest two shots on different shooting lines to the ocean bottom seismometer on the seafloor. Comparison of the differences in direct wave travel‐times from seismic data with and without clock drift correction shows that linear drift correction can commendably eliminate them for most of the ocean bottom seismometers, but some large residuals remain for part of the ocean bottom seismometers. These residuals are inferred to be possibly related to the non‐linear component of clock drift caused by the temperature shock during deployment. The non‐linear clock drift usually varies from different ocean bottom seismometers and thus difficult to correct completely and resultantly will influence the timing accuracy somewhat. Hence, higher uncertainties are suggested to assign to the picked travel‐times for the ocean bottom seismometer with possibly larger non‐linear clock drift in the further tomographic modelling and the internal clock for this kind of ocean bottom seismometer should be recalibrated prior to the next experiment if which requires higher timing accuracy.

Olivine Fabrics in the Oceanic Lithosphere Constrained by Pn Azimuthal Anisotropy

AbstractWe constrain olivine fabrics in the oceanic lithosphere using active and passive seismic observations of Pn azimuthal anisotropy. We first analyze active‐source data derived from a broadband ocean bottom seismometer array that was deployed in the Northwestern Pacific. We infer the azimuthal dependence of the Pn velocity, including the 2ϕ and 4ϕ terms of the sinusoidal functions, where ϕ is the back azimuth. We observe a skewed azimuthal dependence of the Pn velocity, with a large peak‐to‐peak amplitude of about 12%. Similar features are observed using an independent data set of teleseismic Pn waves. We constrain the direction of the crystallographic axes of olivine to explain the observed azimuthal dependence and identify A‐type olivine with a slightly dipping a‐axis and slightly tilting b‐axis being compatible with our observations. In contrast, we find that D‐type olivine with any direction of crystallographic axes cannot explain our observations. Secondary deformation and recrystallization in the older Pacific may be responsible for this strong and complex seismic anisotropy in the study region.

Seismic observation of an active detachment faulting system beneath the Longqi hydrothermal field at the ultraslow spreading southwest Indian ridge

Hydrothermal processes in detachment settings at slow and ultraslow spreading ridges differ greatly from those at melt-rich faster spreading ridges. Active detachment faulting provides the possibility for off-axis high-temperature hydrothermal vents located far away from the heat source beneath the axial volcanic ridge. Seismic data from the slow spreading Mid-Atlantic ridge revealed that hydrothermal fluids may exploit detachment faults to extract heat from a melt zone near the crust-mantle interface. However, knowledge of the subsurface structure and the kinematic processes of detachment faults, and their interaction with hydrothermal fields at the slowest spreading ridges is still insufficient. Here, we report new insights on the seismicity beneath the Longqi hydrothermal field at the ultraslow spreading southwest Indian ridge from three ocean bottom seismometer monitoring experiments. The seismicity outlines the subsurface geometry of a detachment faulting system (DF1 and DF2). The strongly flexed DF1 (45°) is a mature detachment fault with a domed-shaped OCC where ultramafic rocks are exposed on the seafloor. An active young DF2 with intense earthquake activity along the steep subsurface (67°) suggests the initiation phase of rotation. The diversity of hydrothermal activities in the Longqi detachment-hydrothermal systems is closely related to the evolution of detachment faults. Additionally, both the non-transform discontinuity on the western margin of the Longqi-1 field and local faults facilitate hydrothermal circulation. Our study provides baseline observations for a Longqi-type of hydrothermal circulation at the inside corner associated with detachment faulting and non-transform offsets.

Event recognition in marine seismological data using Random Forest machine learning classifier

SUMMARYAutomatic detection of seismic events in ocean bottom seismometer (OBS) data is difficult due to elevated levels of noise compared to the recordings from land. Popular deep-learning approaches that work well with earthquakes recorded on land perform poorly in a marine setting. Their adaptation to OBS data requires catalogues containing hundreds of thousands of labelled event examples that currently do not exist, especially for signals different than earthquakes. Therefore, the usual routine involves standard amplitude-based detection methods and manual processing to obtain events of interest. We present here the first attempt to utilize a Random Forest supervised machine learning classifier on marine seismological data to automate catalogue screening and event recognition among different signals [i.e. earthquakes, short duration events (SDE) and marine noise sources]. The detection approach uses the short-term average/long-term average method, enhanced by a kurtosis-based picker for a more precise recognition of the onset of events. The subsequent machine learning method uses a previously published set of signal features (waveform-, frequency- and spectrum-based), applied successfully in recognition of different classes of events in land seismological data. Our workflow uses a small subset of manually selected signals for the initial training procedure and we then iteratively evaluate and refine the model using subsequent OBS stations within one single deployment in the eastern Fram Strait, between Greenland and Svalbard. We find that the used set of features is well suited for the discrimination of different classes of events during the training step. During the manual verification of the automatic detection results, we find that the produced catalogue of earthquakes contains a large number of noise examples, but almost all events of interest are properly captured. By providing increasingly larger sets of noise examples we see an improvement in the quality of the obtained catalogues. Our final model reaches an average accuracy of 87 per cent in recognition between the classes, comparable to classification results for data from land. We find that, from the used set of features, the most important in separating the different classes of events are related to the kurtosis of the envelope of the signal in different frequencies, the frequency with the highest energy and overall signal duration. We illustrate the implementation of the approach by using the temporal and spatial distribution of SDEs as a case study. We used recordings from six OBSs deployed between 2019 and 2020 off the west-Svalbard coast to investigate the potential link of SDEs to fluid dynamics and discuss the robustness of the approach by analysing SDE intensity, periodicity and distance to seepage sites in relation to other published studies on SDEs.

Characteristics of the Oceanic Ambient Seismic Noise Around Tristan da Cunha in the South Atlantic From OBS Data

AbstractIt has long been recognized that specific ambient seismic noise is generated within the ocean. Therefore, ocean bottom seismometers have advantages to reveal the noise characteristics directly in the deep sea. We use the seismic recordings around Tristan da Cunha collected by 20 OBSs and two island stations to investigate the ambient noise sources in the middle of the South Atlantic. We analyzed the power spectral density, ambient noise cross‐correlations, and correlations with wave heights model. While primary microseisms (∼10–20 s) in this region are deficient, the secondary microseisms are strong and dominate the ambient noise at periods of ∼1–10 s, with the peak broadening to shorter periods (1–3 s). The secondary microseisms are strongly correlated with the local wave heights, especially at shorter periods (∼1–3 s), suggesting that the microseisms are affected mostly by local wave interactions. We observed clear Empirical Green's functions at periods of ∼3–10 s, indicating the existence of distant sources of secondary microseisms. Results of asymmetry analysis of cross‐correlations suggest the prominent direction of the distant source should be located to the Southwest and show barely seasonal variations. The correlations with wave heights indicate a possible source region in the southeast Pacific Ocean. Therefore, besides the dominating local wave‐wave interactions, the secondary microseisms around Tristan region are affected simultaneously by persistent distant sources from the deep basin in the Southeast Pacific Ocean, which might be related to the intense wave‐wave interactions caused by the blocking effect of the Drake Passage on ocean currents and waves.

On the seismic response and earthquake-triggered failures of subaqueous slopes in Swiss lakes

SUMMARY Seismically triggered subaqueous mass movements in lakes may generate tsunamis that can cause significant damage on the shore. In this study, we assess the seismic response and stability of subaqueous slopes in Swiss lakes based on recorded seismological data, historical and geological information and geotechnical surveys. We performed seismic investigations at multiple locations in Lake Lucerne using Ocean Bottom Seismometers (OBS). For these locations, we derived ground-motion amplification functions from local and regional earthquakes and horizontal-to-vertical spectral ratios (H/V) from the earthquake and ambient vibration recordings. The results show (1) very high amplification levels, often exceeding values of 50–100 in the frequency range between 1 and 10 Hz, (2) the fundamental frequency of resonance in the range of 0.5–3.5 Hz and (3) laterally variable site response even for closely located stations. We sought also the signatures of non-linear site response in the H/V curves or ground-motion amplification functions but found only weak indicative effects and no clear evidence. This is most likely due to the low levels of ground motion recorded during the OBS campaigns. We conducted back analyses of historical earthquakes in Switzerland with available documental and/or geological evidence of induced (tsunamigenic) subaqueous slope failures in Swiss lakes. The data set of historical events was complemented with a selection of instrumentally recorded earthquakes in Switzerland. For the analyses, we selected multiple sites in Swiss lakes which failed in the past or are prone to failure in the future. We modelled the ground motion at these locations assuming Swiss standard reference rock conditions (vs30 = 1105 m s−1). The modelled ground motion intensity measures (IM) included peak ground acceleration (PGA), peak ground velocity (PGV) and pseudospectral acceleration (PSA) at 0.3, 1 and 2 s. We estimated the minimum ground motion and macroseismic intensity at reference rock conditions required to trigger the failures of subaqueous slopes. In addition, we defined a threshold for the seismic triggering of such failures in terms of moment magnitude (Mw) and epicentral distance (Re) as: $$\begin{eqnarray} M_{\rm w}=2.891+1.904\log_{10}(R_e+5.166)\: {\rm for}\: R_e\ge 3.7\: {\rm km}. \end{eqnarray}$$ Our results are consistent with previous studies based on worldwide observations. Furthermore, we related the modelled ground motions to the Swiss seismic hazard products and estimated the return period of critical ground shaking responsible for triggering subaqueous slope failures (with potential for tsunami generation) to be in the range of 36–224 yr. Finally, based on previously collected geotechnical data (in situ Cone Penetration Testing and laboratory sediment analysis), we determined the most likely values of the seismic coefficient k to be used with the ground motion IMs modelled at reference rock conditions in infinite slope stability analyses to estimate the factor of safety (FS). For PGA, we found a k = 1; for PGV, k = 2; for PSA0.3s, k = 0.6; for PSA1s, k = 2 and for PSA2s, k = 5.5. These estimates are conservative and affected by the trade-off between the thickness of unconsolidated sediments and the slope angle. Thus, we recommend applying them to slopes with a low-to-moderate gradient (<15°) and sediment thickness of more than 2 m.

CASIE21-OBS: An Open-Access, OBS Controlled-Source Seismic Data Set for Investigating the Structure and Properties of the Cascadia Accretionary Wedge and the Downgoing Explorer-Juan de Fuca-Gorda Plate System

Abstract Geological processes at subduction zones and their associated geohazards (e.g., megathrust earthquakes, submarine landslides, tsunamis, and arc volcanism) are, to a large extent, controlled by the structure, physical properties and fluid content of the subducting plate, the accreted sediments, and the overriding plate. In these settings, modern seismic modeling and imaging techniques based on controlled-source, multicomponent ocean-bottom seismometer (OBS) data are some of the best tools available for determining the subseafloor elastic properties, which can be linked to the aforementioned properties. Here, we present CASIE21-OBS, a controlled-source marine wide-angle OBS data set recently collected across the Cascadia convergent margin as part of the larger CAscadia Seismic Imaging Experiment 2021 (CASIE21). The main component of CASIE21 is a long-offset multichannel seismic (MCS) survey of the Cascadia margin conducted in June–July 2021 onboard R/V M.G. Langseth (cruise MGL2104) aiming to characterize the incoming plate, the plate interface geometry and properties, and the overlying sediment stratigraphy and physical properties. CASIE21-OBS was conducted during R/V M.G. Langseth cruise MGL2103 (May 2021) and R/V Oceanus cruise OC2106A (June–July 2021). It consisted of 63 short-period four-component OBSs deployed at a total 120 stations along 10 across-trench profiles extending from ∼50 km seaward of the deformation front to the continental shelf, and from offshore northern Vancouver Island to offshore southern Oregon. The OBSs recorded the airgun signals of the CASIE21-MCS survey as well as natural seismicity occurring during the deployment period (24 May 2021 19:00 UTC–9 July 2021 09:00 UTC). The OBS data are archived and available at the Incorporated Research Institutions for Seismology Data Management Center under network code YR_2021 for continuous time series (miniSEED) and identifier 21-008 for assembled data set (SEG-Y).

Deep crustal structures with reverse time migration applied to offshore wide-angle seismic data: Equatorial and North-West Brazilian margins

We image the Moho discontinuity and deep crustal layers by applying the Reverse Time Migration (RTM) method to wide-angle seismic (WAS) data that were acquired along two different profiles in NW and SE offshore Brazil, where the spacing between OBS is 12.5 km and 150 m between shots. The application of this method is quite uncommon to ocean bottom seismometers (OBS) data due to the OBS wide spacing deployment and low folds. We analyze the effectiveness of the RTM method when applied to the reflectivity of the WAS data. We imaged the structures by cross-correlating the forward and backward wavefields given by the acoustic seismic equation. This allowed us to use each interface crossed by a ray as both a source and a receiver. The velocity models used to perform RTM were previously obtained by applying a procedure of two-dimensional forward ray-tracing followed by a damped least-squares travel time inversion. The results obtained have an unexpected large contribution from the wavefield traveling as refractions within the Earth and, because of that, we can talk about recovering refractors from the RTM results. We obtained strong and continuous refractors for depths of 7–15 km that correspond to the basement and the Moho discontinuity. As we move landwards, the refractors that correspond to the Moho discontinuity disappear due to the deepening of the surface. Finally, we investigate how the slight change on the velocity model influences the result of the RTM method. The obtained results are promising for a wide range of applications at a crustal scale seismic exploration, with wide-angle seismic data.

Mantle heterogeneity caused by trapped water in the Southwest Basin of the South China Sea

Water is the most common volatile component inside the Earth. A substantial amount of water can be carried down to the interior of the Earth by subducting plates. However, how the subducted water evolves after the subducting slab breaks off remains poorly understood. Here we use the data from a passive seismic experiment using ocean bottom seismometers (OBSs) together with the land stations to determine the high-resolution, three-dimensional seismic structure of the Southwest Sub-basin (SWSB) of the South China Sea (SCS). At depths below 40 km, the mantle shear velocity (Vsv) beneath the northern side of the SWSB is similar to that of the conventional oceanic pyrolite mantle, but roughly 3% shear-velocity reduction is found beneath the southern side of the SWSB. Results of thermal dynamic modeling reveal that the observed shear-velocity reduction could be explained by the presence of 150–300 ppm of water and 5–10% of lower continental crust. The inferred high-water content at the southern side of the SWSB is consistent with a model in which the Proto-SCS plate subducted southward prior to and during the formation of the SCS basin, releasing water into the upper mantle of the SWSB.

Wide-angle seismic imaging of divergent and transform segments of the Pará-Maranhão-Barreirinhas-Ceará margin, NW Brazil

The structure of the North-East equatorial Brazilian margin was investigated during the MAGIC (Margins of brAzil, Ghana and Ivory Coast) seismic wide-angle experiment. This study focuses on the MC5 profile, that spans NW-SE 720 km in length, from the São Paulo Double Fracture Zone to the Barreirinhas margin and continental Borborema province. Its main objective is to understand the fundamental processes which lead to the thinning and finally to the breakup of the continental crust in a specific context of a divergent pull-apart system with two strike-slip borders. The experiment was devised to obtain the 2D structure along this profiles from joint pre-stack depth migration of streamer data and travel-time inversion by forward modeling of 43 Ocean Bottom Seismometers and 21 Land Seismic Stations records. Along the MC5 wide-angle transect, 4 major sectors are identified: 1) the São Paulo Double Fracture Zone presenting a 4.5 km thick volcano-sedimentary Basin on top of a 5.5 km thick basement; 2) a volcanic alignment and intermediate domain SE-ward, formed by the 4.5 km thick Basin III; 3) the 7.5 km thick Basin II, and the 5.5 km thick Basin I composing the continental slope and shelf. While all the offshore basement remains about 6 km thick in the deep-sea domains, acoustic velocity evolves from two-layer 4.8–6 km/s and 6.1–6.8 km/s beneath Basin III to two-layer high velocity 6.1–6.8 km/s and 7.2–7.4 km/s beneath Basin II and I, The necking zone, forming the Parnaiba Platform and associated Piaui-Camocim and Ceará Basins, is 50 km wide; 4) the Médio Coreaú and Ceará Central thrust belt, where the un-thinned continental crust thickness reaches 32 km. Finally, a schematic kinematic reconstruction that satisfies these observation is argumented.

A Novel Method to Estimate Orientations of an Ocean-Bottom Seismometer Array for Accurate Measurement of Waveform Phases and Amplitudes

Abstract The misorientation of three-component seismometers restricts the application of relevant seismic experiments such as ocean-bottom seismometer (OBS) arrays. Previous orientation determination relied on an assumption that the propagation azimuth of seismic waves follows the great-circle path (GCP) azimuth. This assumption may yield systematic errors in the estimated orientation, particularly when the ray paths are bent laterally due to velocity heterogeneity in the Earth. Here, we develop a new method for unbiasedly estimating the horizontal orientations of seismic sensors and apply this method to the Blanco transform fault OBS experiment. We first retrieve the orientations relative to the propagation azimuths from the recorded Rayleigh and P waveforms, and then determine the geographic north orientations by calculating the propagation azimuths via an Eikonal-equation-based phase-tracking method that theoretically accounts for the effect of ray bending. Synthetics test validates that the phase-tracking method can retrieve unbiased propagation azimuths of seismic waves. The final results derived from Rayleigh- and P-wave polarization analyses with the respective phase-tracked propagation azimuths are more consistent and the orientation errors are smaller, indicating the robustness and accuracy of this method. Comparing the orientations from our phase-tracking method to those from the GCP assumption, the deviation can reach up to 8° between these two techniques in the study region. Subsequently, when orientations of the synthetics modeled from three-dimensional elastic waveform simulation are deviated according to the GCP-predicted orientations, we find nonnegligible bias in the phase and amplitude measurements that could reduce the accuracy and resolution of following inversion, which indicates the significance of our phase-tracking method in accurate orientation of OBS arrays as well as inland seismic experiments.

Characteristics of Background Noise in the Oldest-1 Array Deployed on the Oldest Part of the Pacific Plate

ABSTRACT We characterize the sources of background noise recorded by an ocean-bottom seismometer array deployed on the oldest oceanic lithosphere of the western Pacific plate. This igneous oceanic crust is blanketed up to 1000 m cover of sediments, making it a unique place to investigate correlations between sediment thickness and background noise levels in the marine environment. Although the data are limited, we found that the thicker sediments attenuated short-period noise levels and amplified long-period noise levels, in agreement with observations made for seismic arrivals in onshore sedimentary basins. Comparison of the noise levels from our Oldest-1 deployment with those (1) in the western Pacific basin, (2) near the Mariana trench, and (3) offshore New Zealand indicates that the Oldest-1 has the highest noise levels in the double-frequency (secondary) microseism band by up to 7 dB. This result was unexpected given the greater water depths (averaging 5800 m) for the Oldest-1 array. However, even at these great depths, we find a strong correlation between significant wave heights and background noise levels in the short-period microseism band (2–5 s). In addition, the Oldest-1 array was subjected to a great number of intense close typhoons than the other arrays and was proximal to a great number of taller seamounts than the other arrays. Finally, tracking noise source azimuths using station pairs led to the identification of a possible new source for the single-frequency microseism band in the coastal areas of Papua New Guinea and Australia.

Gravity inversion constrained by OBS receiver function reveals crustal structure in Ryukyu Trench

The Ryukyu Trench is located in the northern Philippine Sea. The subduction of the Philippine Sea Plate to the Eurasia Plate along the Ryukyu Trench in the NNW direction leads to extremely complex tectonic processes such as subduction, collision, accretion, and back-arc expansion. It is of great significance to obtain the crustal structure of the Ryukyu Trench for understanding the crustal deformation characteristics, subduction direction and the scale of the Philippine Sea Plate. The ocean bottom seismometer (OBS) could provide important information on the deep structure of the oceanic area, and the wide coverage of gravity anomaly can provide more constraints for the regional crustal structure. In this study, we put the time-delay term of the sedimentary layer into the H-κ stacking formula, which eliminates the influence of the sedimentary layer in traditional H-κ stacking. Based on the accurate crustal thickness at OBS locations, a gravity inversion method is proposed to solve the problem of non-uniqueness in gravity data and limited distribution of OBSs. By using this method, we obtain a crustal structure with higher resolution than Crust1.0 model in Ryukyu Trench. The results show that the crustal thickness in the Ryukyu Trench region thickens along the NNW subduction direction, revealing that the Philippine Plate’s subduction direction towards the Eurasian Plate is in the NNW direction.

The effect of deep ocean currents on ocean- bottom seismometers records

Abstract. Ocean-bottom seismometers (OBSs) are usually deployed for seismological investigations, but these objectives are impaired by noise resulting from the ocean environment. We split the OBS-recorded seismic noise into three bands: short periods, microseisms and long periods, also known as tilt noise. We show that bottom currents control the first and third bands, but these are not always a function of the tidal forcing. Instead, we suggest that the ocean bottom has a flow regime resulting from two possible contributions: the permanent low-frequency bottom current and the tidal current. The recorded noise displays the balance between these currents along the entire tidal cycle, between neap and spring tides. In the short-period noise band, the ocean current generates harmonic tremors corrupting seismic dataset records. We show that, in the investigated cases, the harmonic tremors result from the interaction between the ocean current and mechanical elements of the OBS that are not essential during the sea bottom recording and thus have no geological origin. The data from a new broadband OBS type, designed and built at Instituto Dom Luiz (IDL – University of Lisbon)/Centre of Engineering and Product Development (CEIIA), hiding non-essential components from the current flow, show how utmost harmonic noise can be eliminated.