Shallow Seismic Structure Research Articles

Array analysis of P-wave polarization data for constraining near-surface shear wave speeds: Aplication to theValley of Mexico seismic network

SUMMARY We present an array method for constraining near-surface shear wave speed values from P-wave polarization data. As the direction of P-wave particle motion is determined both by the incident ray parameter and the local surficial shear wave speed, the latter can be inferred from three-component seismic records if the ray parameter is estimated first. Previous studies have analysed data recorded at a single station from various earthquakes in order to estimate near-surface shear wave speed. In contrast, this paper analyses data recorded at an array of stations. As such, this method can estimate both ray parameter data for the recorded earthquakes and near-surface shear wave speed values at the studied sites. And as an advantage, the method can analyse records of local and regional earthquakes where ray parameter information cannot be reliably estimated using 1D reference velocity models. Additionally, by processing polarization data spanning multiple events and stations, our method provides robust estimates for near-surface shear wave speed values. We apply this array-based method to a seismic array covering the Valley of Mexico. Analysed records correspond to regional earthquakes with epicentral distances of roughly 300 km. The corresponding ray parameter values for all earthquakes were estimated from the delay of the seismic waves arrival throughout the network as all stations have a reliable time reference. Near-surface shear wave speed values estimated from the array-based polarization analysis are correlated with major geological features of the area such as the Trans-Mexican Volcanic Belt and the local basin. We also apply the presented method to a minimal data scenario. The flexibility of this method allows it to employ a single ray parameter value as input, producing robust estimates for near-surface shear wave speed values within our case study. This feature becomes useful in a variety of scenarios, such as utilizing seismic data lacking timing information. The presented method constitutes an accessible and inexpensive alternative to study shallow seismic structure, opening opportunities to expand and complement current tools and techniques to assess seismic hazards.

CSRM‐1.0: A China Seismological Reference Model

AbstractHigh‐resolution seismic model is crucial for advancing our understandings on geological processes and enhancing seismic hazard mitigation programs. We construct a high‐resolution China Seismological Reference Model (CSRM‐1.0) in the top 100 km of the crust and uppermost mantle in continental China following a top‐down construction process. The employed seismic constraints include P‐wave polarization angle from tele‐seismic event, short‐period Rayleigh wave ellipticity from ambient noise, long‐period Rayleigh wave ellipticity from earthquake data, receiver function, empirical Green's function from ambient noise, Rayleigh wave phase/group velocity dispersion curves from regional earthquakes, and Pn‐wave travel time extracted from seismic data of 4,435 stations. CSRM‐1.0 has a spatial crustal resolution of ∼60 km beneath the north‐south seismic belt and trans‐North China orogen regions and ∼120 km beneath the rest of continental China, and a spatial mantle resolution of ∼300 km. CSRM‐1.0 exhibits prominent velocity heterogeneities in the crust and uppermost mantle and an eastward thinning of the crust, geographically correlating with geological settings. CSRM‐1.0 improvements include accurate estimation of shallow seismic structure, increased spatial resolution and improved model accuracy. Crustal composition inferred from CSRM‐1.0 exhibits a general transition from a felsic upper crust to a mafic lower crust. Mafic rocks in the lower crust are found predominantly along inter‐block boundaries and sporadically within the interiors of blocks, likely resulted from preferential inter‐block intrusions of magmas related to various oceanic plate subductions and the Emeishan mantle plume. This study contributes seismic constraints and CSRM‐1.0 to the CSRM product center (http://chinageorefmodel.org) as a backbone open‐access geophysical cyberinfrastructure.

Differences in Scattering Properties of the Shallow Crusts of Earth, Mars, and the Moon Revealed by P‐Wave Receiver Functions

AbstractThe scattering properties of terrestrial planetary bodies can provide valuable insights into their shallow seismic structure, meteoritic impact history, and geological activity. Scattering properties of the shallow crusts of Earth, Mars, and the Moon are investigated by constructing P‐wave receiver functions (PRFs) from teleseismic waveforms with high signal‐to‐noise ratios. The authors’ analysis reveals that strong coda waves lead to significant variations in the PRF waveforms calculated using different time windows, and the stability of the PRF is primarily influenced by the fractional velocity fluctuation. Synthetic PRFs for various scattering media confirm these observations. Comparing the observed and synthetic PRFs, it is found that the fractional velocity fluctuation in the shallow crust is greater than ∼0.2 for the Moon but less than ∼0.2 for Earth and Mars. The authors further discuss possible mechanisms that could have affected the fractional velocity fluctuation and suggest that the distinct fractional velocity fluctuation between the Moon and Earth/Mars is mainly due to differences in the water content of the crustal rocks of the three planetary bodies.

Lateral Variations Across the Southern San Andreas Fault Zone Revealed From Analysis of Traffic Signals at a Dense Seismic Array

AbstractWe image the shallow seismic structure across the Southern San Andreas Fault (SSAF) using signals from freight trains and trucks recorded by a dense nodal array, with a linear component perpendicular to SSAF and two 2D subarrays centered on the Banning Fault and Mission Creek Fault (MCF). Particle motion analysis in the frequency band 2–5 Hz shows that the examined traffic sources can be approximated as moving single‐ or multi‐point sources that primarily induce Rayleigh waves. Using several techniques, we resolve strong lateral variations of Rayleigh wave velocities and Q‐values across the SSAF, including 35% velocity reduction across MCF toward the northeast and strong attenuation around the two fault strands. We further resolve 10% mass density reduction and 45% shear modulus decrease across the MCF. These findings suggest that the MCF is currently the main strand of the SSAF in the area with important implications for seismic hazard assessments.

Shallow seismic structure around the Vrancea Seismic Zone from joint inversion of ambient noise H/V ratios and surface wave dispersion

The South-East Carpathians in Romania host the Vrancea Seismic Zone (VSZ), an intraplate seismic nest of frequent high magnitude earthquakes that cause unusually high damage over extended areas. Estimating the amplification of seismic signals strongly depends on the depth to the geophysical bedrock and the complex near surface geological structures. To improve seismic hazard estimates and ultimately mitigate seismic risk, we determine the frequency of resonance at 81 broadband seismic stations in southeast Romania using the horizontal-to-vertical spectral ratio (HVSR) of ambient seismic noise. We jointly invert HVSR functions with Rayleigh wave phase velocity dispersion curves using a Markov-chain Monte Carlo algorithm under the assumption of a diffuse field. Our results show excellent correlation between the fundamental frequency of resonance and local surface geology: stations located on foreland sedimentary basins show HVSR peaks in the range 0–1 Hz, while those located in orogenic nappes reach values as high as 10 Hz. Seismic bedrock topography mirrors the crystalline bedrock depth inferred from gravity modeling, but is consistently ∼2000 m shallower, corresponding to the interface between Paleozoic-Triassic to Jurassic-Cretaceous sedimentary levels on the Moesian Platform constrained by borehole data. The near surface seismic model shows anomalously low velocities in the Focsani forearc basin, adjacent to the VSZ, consistent with extremely low resonance frequencies and deep sedimentary layer successions (>5 km). In the back-arc of the Carpathian Orogen, at the contact with the Transylvanian Basin, average resonance frequencies are <2 Hz and inferred seismic bedrock depths vary laterally, likely affected by buried volcanic structures and diapiric salt domes. Our study presents the largest endeavor to study anomalous intensity patterns and local site effects in and around the VSZ and has fundamental implications for seismic risk evaluation and future civil engineering design planning.

Subsurface Imaging With Ocean‐Bottom Distributed Acoustic Sensing and Water Phases Reverberations

AbstractSeismic waves from earthquakes recorded on the seafloor are composed of complex multiple arrivals. Here, distributed acoustic sensing (DAS) observations along a cable located offshore the Sanriku Coast, Japan, show that the local earthquake wavefield is particularly rich in Scholte waves. We introduce a processing pipeline to extract these surface waves from DAS records. We then invert hundreds of dispersion curves along a section of the cable to form a shallow high‐resolution shear‐wave velocity model. Moreover, we focus on the possible generation mechanisms of Scholte waves through a series of 2D and 3D full‐wavefield numerical simulations. We show that water phase reverberations greatly contribute to the generation of Scholte waves on the ocean floor. This study demonstrates the potential of DAS to observe and better understand a poorly known marine wave phenomenon and image the offshore shallow seismic structure with an unprecedented spatial resolution.

Shallow seismic structure beneath the continental China revealed byP-wave polarization, Rayleigh wave ellipticity and receiver function

SUMMARYWe construct a high-resolution shallow 3-D seismic model in the top 10 km of the upper crust in the continental China, with constraints of P polarization, Rayleigh wave ellipticity and receiver function obtained from records of 3848 seismic stations. Our 3-D seismic model has a spatial resolution of 0.6–1.2° in the north–south seismic belt and the trans-north China orogen, and 1–2° in the rest of the continental China (except the Tarim basin and the southwest Tibet). The seismic model exhibits low velocity anomalies of deposits in major sedimentary basins and high velocity anomalies of crustal bedrocks in young orogenic belts and old tectonic blocks. The inferred sediment thickness maps display thick deposits in major sedimentary basins, some compacted sediments in the intermontane basins in young orogenic belts and little sediments in old tectonic blocks. We also discuss compaction effects of the sediments and implications of tectonic history and geological evolution of the major basins in the continental China based on the inferred seismic models. This study provides an effective mean of seismic imaging through joint inversion of various seismic constraints and establishes a framework of seismic data sharing for future studies in the seismological community in a first step of developing a China Seismological Reference Model.

Surface-wave dispersion spectrum inversion method applied to Love and Rayleigh waves recorded by distributed acoustic sensing

Recently, distributed acoustic sensing (DAS) has been applied to shallow seismic structure imaging providing dense spatial sampling at a relatively low cost. DAS on a standard straight fiber-optic cable mostly records axial dynamic strain, which makes it difficult to separate the Rayleigh and Love wavefields. As a result, the mixed Rayleigh and Love wave signals cannot be used in the conventional surface-wave dispersion inversion method. Therefore, it is often ensured that the source and the cable are in the same line and only Rayleigh wave dispersion is used, which limits the constraints on structure and model resolution. We have inverted surface-wave dispersion spectra instead of dispersion curves. This inversion method can use mixed Rayleigh and Love waves recorded when the source and receiver array are not aligned. The multiple-channel records are transformed to the frequency domain, and a slant stack method is used to construct the dispersion spectra. The genetic algorithm method is used to obtain an optimal S-wave velocity model that minimizes the difference between theoretical and observed dispersion spectra. A series of synthetic tests are conducted to validate our method. The results suggest that our method not only improves the flexibility of the acquisition system design, but the Love wave data also provide additional constraints on the structure. Our method is applied to the active source and ambient noise data sets acquired at a geothermal site and provides consistent results for different data sets and acquisition geometries. The sensitivity of the dispersion spectra to layer thickness, density, and P-wave velocity is also discussed. With our method, the amount of usable data can be increased, helping deliver better subsurface images.

Seismogenic Source Model of the 2019, Mw 5.9, East-Azerbaijan Earthquake (NW Iran) through the Inversion of Sentinel-1 DInSAR Measurements

In this work, we investigate the Mw 5.9 earthquake occurred on 7 November 2019 in the East-Azerbaijan region, in northwestern Iran, which is inserted in the tectonic framework of the East-Azerbaijan Plateau, a complex mountain belt that contains internal major fold-and-thrust belts. We first analyze the Differential Synthetic Aperture Radar Interferometry (DInSAR) measurements obtained by processing the data collected by the Sentinel-1 constellation along ascending and descending orbits; then, we invert the achieved results through analytical modelling, in order to better constrain the geometry and characteristics of the seismogenic source. The retrieved fault model shows a rather shallow seismic structure, with a center depth at about 3 km, approximately NE–SW-striking and southeast-dipping, characterized by a left-lateral strike-slip fault mechanism (strike = 29.17°, dip = 79.29°, rake = −4.94°) and by a maximum slip of 0.80 m. By comparing the inferred fault with the already published geological structures, the retrieved solution reveals a minor fault not reported in the geological maps available in the open literature, whose kinematics is compatible with that of the surrounding structures, with the local and regional stress states and with the performed field observations. Moreover, by taking into account the surrounding geological structures reported in literature, we also use the retrieved fault model to calculate the Coulomb Failure Function at the nearby receiver faults. We show that this event may have encouraged, with a positive loading, the activation of the considered receiver faults. This is also confirmed by the distribution of the aftershocks that occurred near the considered surrounding structures. The analysis of the seismic events nucleated along the left-lateral strike-slip minor faults of the East-Azerbaijan Plateau, such as the one analyzed in this work, is essential to improve our knowledge on the seismic hazard estimation in northwestern Iran.

A single Rayleigh mode may exist with multiple values of phase-velocity at one frequency

SUMMARYOther than commonly assumed in seismology, the phase velocity of Rayleigh waves is not necessarily a single-valued function of frequency. In fact, a single Rayleigh mode can exist with three different values of phase velocity at one frequency. We demonstrate this for the first higher mode on a realistic shallow seismic structure of a homogeneous layer of unconsolidated sediments on top of a half-space of solid rock (LOH). In the case of LOH a significant contrast to the half-space is required to produce the phenomenon. In a simpler structure of a homogeneous layer with fixed (rigid) bottom (LFB) the phenomenon exists for values of Poisson’s ratio between 0.19 and 0.5 and is most pronounced for P-wave velocity being three times S-wave velocity (Poisson’s ratio of 0.4375). A pavement-like structure (PAV) of two layers on top of a half-space produces the multivaluedness for the fundamental mode. Programs for the computation of synthetic dispersion curves are prone to trouble in such cases. Many of them use mode-follower algorithms which loose track of the dispersion curve and miss the multivalued section. We show results for well established programs. Their inability to properly handle these cases might be one reason why the phenomenon of multivaluedness went unnoticed in seismological Rayleigh wave research for so long. For the very same reason methods of dispersion analysis must fail if they imply wave number kl(ω) for the lth Rayleigh mode to be a single-valued function of frequency ω. This applies in particular to deconvolution methods like phase-matched filters. We demonstrate that a slant-stack analysis fails in the multivalued section, while a Fourier–Bessel transformation captures the complete Rayleigh-wave signal. Waves of finite bandwidth in the multivalued section propagate with positive group-velocity and negative phase-velocity. Their eigenfunctions appear conventional and contain no conspicuous feature.

Frequency‐Dependent P Wave Polarization and Its Subwavelength Near‐Surface Depth Sensitivity

AbstractNear‐surface structure is crucial to assessing seismic hazards and understanding earthquakes and surface processes yet is a major challenge to robustly image. Recently, an approach based on body‐wave polarization was introduced for constraining shallow seismic structure, but the depth sensitivity of the polarization measurement has remained unclear. Using waveform simulations based on a layer over a half space, we find that the depth sensitivity of P wave polarization peaks at the surface and decreases abruptly over a depth range shorter than its wavelength. A strong frequency dependence provides constraints on local 1‐D structure, with frequencies between 0.1 and 10 Hz illuminating structure at depths of 10 m to several kilometers. Applying these results to teleseismic recordings in Japan provides constraints on structure at about 120 to 750 m, including a distinctive weak zone along the Median Tectonic Line in the Kii peninsula and Awaji Island.

Shallow Shear-Wave Velocity Beneath Jakarta, Indonesia Revealed by Body-Wave Polarization Analysis

Noting the importance of evaluating near-surface geology in earthquake risk assessment, we explored the application to the Jakarta Basin of a relatively new and simple technique to map shallow seismic structure using body-wave polarization. The polarization directions of P-waves are sensitive to shear-wave velocities (Vs), while those of S-waves are sensitive to both body-wave velocities. Two dense, temporary broadband seismic networks covering Jakarta city and its vicinity were operated for several months, firstly, from October 2013 to February 2014 consisting of 96 stations, and secondly, between April and October 2018 consisting of 143 stations. By applying the polarization technique to earthquake signals recorded during these deployments, the apparent half-space shear-wave velocity (Vsahs) beneath each station is obtained, providing spatially dense coverage of the sedimentary deposits and the edge of the basin. The results showed that spatial variations in Vsahs obtained from polarization analysis are compatible with previous studies, and appear to reflect the average Vs of the top 150 m. The low Vs that characterizes sedimentary deposits dominates most of the area of Jakarta, and mainly reaches the outer part of its administrative margin to the southwest, more than 10 km away. Further study is required to obtain a complete geometry of the Jakarta Basin. In agreement with previous studies, we found that the polarization technique was indeed a simple and effective method for estimating near-surface Vs that can be implemented at very low-cost wherever three-component seismometers are operated, and it provides an alternative to the use of borehole and active source surveys for such measurements. However, we also found that for deep basins such as Jakarta, care must be taken in choosing window lengths to avoid contamination of basement converted phases.

Application of seismic interferometric migration for shallow seismic high precision data processing: A case study in the Shenhu area

The stability of submarine geological structures has a crucial influence on the construction of offshore engineering projects and the exploitation of seabed resources. Marine geologists should possess a detailed understanding of common submarine geological hazards. Current marine seismic exploration methods are based on the most effective detection technologies. Therefore, current research focuses on improving the resolution and precision of shallow stratum structure detection methods. In this article, the feasibility of shallow seismic structure imaging is assessed by building a complex model, and differences between the seismic interferometry imaging method and the traditional imaging method are discussed. The imaging effect of the model is better for shallow layers than for deep layers because coherent noise produced by this method can result in an unsatisfactory imaging effect for deep layers. The seismic interference method has certain advantages for geological structural imaging of shallow submarine strata, which indicates continuous horizontal events, a high resolution, a clear fault, and an obvious structure boundary. The effects of the actual data applied to the Shenhu area can fully illustrate the advantages of the method. Thus, this method has the potential to provide new insights for shallow submarine strata imaging in the area.

Verification of the shallow seismic crustal structure of the western Krušné Hory crystalline unit, Czech Republic

We analyze refraction measurements along a short profile in western Krusne hory crystalline unit. The profile passed close to the seismically active zone of Nový Kostel. The measurements were carried out to distances of about 15 km using quarry blasts near the village of Horni Rozmysl, located at the eastern margin of the crystalline unit. Smoothed P-wave travel times were interpreted using the Wiechert-Herglotz method, which yielded a 1-D velocity model of the shallow crustal structure of the crystalline unit down to a depth of 1.7 km. The P-wave velocity of the model increases from about 4.0 km/s at the surface to 5.9 km/s at a depth of 1.7 km. The superficial velocities of our model are somewhat higher than the superficial velocities of the model that is routinely used for earthquake location in the region.

The Shallow Seismic Reflection Structure in the Northeastern Part of the Sakurajima Volcano with the Pseudo-Reflection Method

The Shallow Seismic Reflection Structure in the Northeastern Part of the Sakurajima Volcano with the Pseudo-Reflection Method

Heterogeneous seismic velocity structure of the upper lithosphere at Kane oceanic core complex, Mid‐Atlantic Ridge

The Kane oceanic core complex (OCC) is a large, corrugated megamullion that was formed by a long‐lived detachment fault at the axis of the Mid‐Atlantic Ridge adjacent to Kane Fracture Zone between 2.1 and 3.3 Ma. We use refracted arrivals recorded along a 6‐km‐long hydrophone streamer during a multichannel seismic survey to constrain the shallow seismic velocity structure of the OCC. Results are presented in high‐resolution traveltime seismic tomographic models along six lines that cover all of the main morphological features of the megamullion. The models show large lateral variability in P wave velocity within the upper ∼0.5–2.0 km of the lithosphere, and these variations correlate to first order with observed variations in lithology, documented by in situ basement samples and seafloor morphology. Lithological interpretation of the velocity models indicates that there is marked lateral variability in distribution of gabbroic intrusions, serpentinized peridotites, and basalts at scales of a few kilometers to ∼10 km. Serpentinized peridotites appear to dominate the central and older parts of the OCC. High‐velocity gabbros are consistently (but not exclusively) present closer to the termination of the Kane detachment fault and toward the ends of the OCC. The structure of the lithosphere exhumed by the Kane detachment fault is far from the standard ophiolite‐based Penrose model, and it does not show segment‐centered magmatism that is commonly interpreted at slow spreading ridges. If the gabbros exhumed toward the termination of the OCC were emplaced deep (∼10 km) beneath the spreading axis, they may have constituted a weak zone that focused initiation of the Kane detachment fault. Alternately, as the OCC footwall was being exhumed the gabbros may have been emplaced because of dynamic changes in melt supply, changes in mantle fertility, or decompression melting. Late stage volcanism is clearly associated with a major high‐angle normal fault that cuts the detachment surface; this volcanism may have been stimulated or enhanced by bending stresses in the bending footwall. The shape of the large‐scale corrugated morphology of the OCC is nearly invariant in the dip direction across major changes in basement lithology, indicating that once established, the form of the Kane detachment fault was highly resistant to modification.

Shallow seismic structure of Kunlun fault zone in northern Tibetan Plateau, China: implications for the 2001Ms8.1 Kunlun earthquake

SUMMARY The shallow seismic velocity structure of the Kunlun fault zone (KLFZ) was jointly deduced from seismic refraction profiling and the records of trapped waves that were excited by five explosions. The data were collected after the 2001 Kunlun Ms8.1 earthquake in the northern Tibetan Plateau. Seismic phases for the in-line record sections (26 records up to a distance of 15 km) along the fault zone were analysed, and 1-D P- and S-wave velocity models of shallow crust within the fault zone were determined by using the seismic refraction method. Sixteen seismic stations were deployed along the off-line profile perpendicular to the fault zone. Fault-zone trapped waves appear clearly on the record sections, which were simulated with a 3-D finite difference algorithm. Quantitative analysis of the correlation coefficients of the synthetic and observed trapped waveforms indicates that the Kunlun fault-zone width is 300 m, and S-wave quality factor Q within the fault zone is 15. Significantly, S-wave velocities within the fault zone are reduced by 30–45 per cent from surrounding rocks to a depth of at least 1–2 km, while P-wave velocities are reduced by 7–20 per cent. A fault-zone with such Pand S-low velocities is an indication of high fluid pressure because Vs is affected more than Vp. The low-velocity and low-Q zone in the KLFZ model is the effect of multiple ruptures along the fault trace of the 2001 Ms8.1 Kunlun earthquake.

Inversion of Scholte wave dispersion and waveform modeling for shallow structure of the Ninetyeast Ridge

The construction of S-wave velocity models of marine sediments down to hundreds of meters below the seafloor is important in a number of disciplines. One of the most significant trends in marine geophysics is to use interface waves to estimate shallow shear velocities which play an important role in determining the shallow crustal structure. In marine settings, the waves trapped near the fluid–solid interface are called Scholte waves, and this is the subject of the study. In 1998, there were experiments on the Ninetyeast Ridge (Central Indian Ocean) to study the shallow seismic structure at the drilled site. The data were acquired by both ocean bottom seismometer and ocean bottom hydrophone. A new type of seafloor implosion sources has been used in this experiment, which successfully excited fast and high frequency (>500 Hz) body waves and slow, intermediate frequency (<20 Hz) Scholte waves. The fundamental and first higher mode Scholte waves have both been excited by the implosion source. Here, the Scholte waves are investigated with a full waveform modeling and a group velocity inversion approach. Shear wave velocities for the uppermost layers of the region are inferred and results from the different methods are compared. We find that the full waveform modeling is important to understand the intrinsic attenuation of the Scholte waves between 1 and 20 Hz. The modeling shows that the S-wave velocity varies from 195 to 350 m/s in the first 16 m of the uppermost layer. Depths levels of high S-wave impedance contrasts compare well to the layer depth derived from a P-wave analysis as well as from drilling data. As expected, the P- to S-wave velocity ratio is very high in the uppermost 16 m of the seafloor and the Poisson ratio is nearly 0.5. Depth levels of high S-wave impedance contrasts are comparable to the layer depth derived from drilling data.

Morphogenesis of the SW Balearic continental slope and adjacent abyssal plain, Western Mediterranean Sea

We present the seafloor morphology and shallow seismic structure of the continental slope south-east of the Balearic promontory and of the adjacent Algero-Balearic abyssal plain from multibeam and chirp sonar data. The main purpose of this research was to identify the sediment pathways from the Balearic promontory to the Algero-Balearic deep basin from the Early Pliocene to the Present. The morphology of the southern Balearic margin is controlled by a SW–NE structural trend, whose main expressions are the Emile Baudot Escarpment transform fault, and a newly discovered WSW–ENE trend that affects the SW end of the escarpment and the abyssal plain. We relate the two structural trends to right-lateral simple shear as a consequence of the Miocene westward migration of the Gibraltar Arc. Newly discovered steep and narrow volcanic ridges were probably enabled to grow by local transtension along the transform margin. Abyssal plain knolls and seahills relate to the subsurface deformation of early stage halokinetic structures such as salt rollers, salt anticlines, and salt pillows. The limited thickness of the overburden and the limited amount of deformation in the deep basin prevent the formation of more mature halokinetic structures such as diapirs, salt walls, bulbs, and salt extrusions. The uppermost sediment cover is affected by a dense pattern of sub-vertical small throw normal faults resulting from extensional stress induced in the overburden by subsurface salt deformation structures. Shallow gas seismic character and the possible presence of an active polygonal fault system suggest upward fluid migration and fluid and sediment expulsion at the seafloor through a probable mud volcano and other piercement structures. One large debris flow deposit, named Formentera Debris Flow, has been identified on the lower slope and rise of the south Formentera margin. Based on current observations, we hypothesize that the landslide originating the Formentera Debris Flow occurred in the Holocene, perhaps in historical times.

Using Ambient Seismic Noise to Determine Short-Period Phase Velocities and Shallow Shear Velocities in Young Oceanic Lithosphere

Using 10 broadband ocean bottom seismometers from the 11-month-long deployment of the Gravity Lineations Intraplate Melting Petrologic and Seismologic Expedition (GLIMPSE) passive seismic experiment located in the south central Pacific, we have estimated the seismic impulse responses from ambient seismic noise for 45 station-to-station paths. The raw impulse responses show moveout with station-to-station distance, and there is a trend of decreasing signal-to-noise ratios as the station-to-station distance increases. The decrease in signal-to-noise ratio is expected as a smaller range of azimuths of propagating surface waves will contribute constructively to the cross-correlated signal with increasing distance, although scattering may also play a role in the coherence of seismic noise at periods less than 16 sec. From these station-to-station paths, we determined group velocities for the fundamental mode Rayleigh waves of a 2–16-sec period and the second mode Rayleigh wave of a 3.5–7-sec period. We calculate phase velocities for the fundamental mode and second mode Rayleigh wave over the same period range as the group velocities by applying a time variable filter to the noise correlation function and carefully unwrapping the phase spectrum of the resulting filtered impulse responses. Within this period range, there is a transition from waves at short periods, whose energy is mostly in the water column, to waves sensitive to crustal and upper mantle structure. The phase velocities for the second mode, which have peak sensitivity in the lower crust and shallow mantle, show evidence for azimuthal anisotropy. The average phase velocities of the station-to-station paths in the east–west direction are 2% faster than the north–south paths at the 4–7-sec period, consistent with the fast directions determined from SKS wave splitting measurements of N100°E. By incorporating the short-period fundamental and higher mode phase velocities from ambient seismic noise with longer period (16–100 sec) fundamental mode Rayleigh-wave phase velocities determined from teleseismic events, we inverted for the average crustal and upper mantle shear velocity structure and water column depth and velocity. The predicted phase velocities are extremely sensitive to the water column compressional velocity. We determined the average water column velocity to be 1466±3 m/sec, in contrast to the average of 1500 m/sec estimated from shipboard measurements weighted according to the Rayleigh-wave sensitivity kernel. The difference may be due to the dispersive effects of scattering by bathymetry or by the thin variable thickness sediment layer. The inversion also produces a Vp/Vs ratio of 1.88 for the crust. This method can provide useful information about the shallow seismic structure of the oceanic crust and uppermost mantle and is an important complement to longer period studies.