Search

Microseismic monitoring has been used in geo-energy related activities, such as shale-gas exploitation, mining, deep geothermal exploitation, geotechnical and structural engineering, for detecting and locating fractures, rock failures and micro-earthquakes. The success of microseismic monitoring depends on reliable detection and location of the recorded microseismicity. Multichannel coherence migration (MCM) is a detection and location waveform migration-based approach which does not require phase picking, identification and association and performs well on noisy data. Its caveat is a high computational cost, which impedes its application of MCM on large data sets or for real-time monitoring. To address this issue, we propose an improved approach, the multichannel coherence migration grid search (MCMgs), by introducing an adaptive grid optimization technique. Based on results from synthetic and real data, we show that MCMgs reduces the computation time up to 64 times. In addition, MCMgs generates multiple maximum coherence values with various grid sizes instead of a single (maximum) coherence value that links to a single gridpoint and size, thus resulting in more accurate locations. Our simulation results on different deployment geometries demonstrate that MCMgs is effective even with a small number of recordings available—a minimum of seven. We conduct a sensitivity analysis to assess how the detectability of events is affected by the spatial arrangement of the deployed monitoring array. If a limited number of seismometers are available for deployment, our analysis favours a patch array deployment geometry. We show that 12 seismometers deployed at a patch array geometry can have similar detection and localization capability as a large rectangular array of more than 100 seismometers but at a much lower computational and deployment cost. [ABSTRACT FROM AUTHOR]

In this work, we present both 1-D and 3-D shear wave velocity (Vs) models of the oceanic crust and uppermost mantle below the Deep OCean Test ARray area, located ∼ 70 km north of the central section of the Gloria Fault, in the eastern North-Atlantic Ocean. The velocity models are inferred from the dispersion of surface waves recorded on ocean bottom seismometers. Dispersion measurements are obtained from the analysis of ambient seismic noise at short periods (< 14 s) and teleseismic surface waves at long periods (> 14 s) using the two-station method. The 1-D Vs model is inferred from the joint inversion of Rayleigh wave phase velocities and Love wave group and phase velocities. The 3-D tomographic model is obtained by inversion of 2-D Love wave group velocity maps as a function of depth, further constrained by the average of Love wave phase velocities obtained from ambient noise (4–9 s) and the average Rayleigh and Love wave phase velocities calculated from teleseismic data (14–44 s). The 1-D Vs model shows a sediment layer with a low velocity of 1.05 km s−1, similar to previous studies in the region. Below the sediments, we find an oceanic crust with velocities ranging from 3.3 to 4.5 km s−1. The model reaches an unusually high velocity of 4.9 km s−1 in a 20 km thick layer at depths between 16 and 36 km. We interpret this fast velocity layer as indicative of the presence of harzburgite, a residue of enhanced melting that might have been formed by the proximity between the Mid-Atlantic Ridge and the Azores mantle plume. At greater depths the velocity decreases, forming a low-velocity zone that reaches a minimum at ∼ 70 km depth, which we interpret as the maximum depth for the lithosphere–asthenosphere boundary. The 3-D model shows a structure that is mostly horizontally layered, with Vs isocontours at 3.5–4.5 km s−1 highlighting oscillations of the crustal structure with wavelengths of ∼25–30 km. These oscillations may be due to changes in the rate of mantle upwelling and magma supply rate. [ABSTRACT FROM AUTHOR]

This study presents the shear wave velocity (VS) structures of sedimentary sequences and a section of the upper crustal layer in the Fukushima forearc region of the Japan Trench subduction zone, which were obtained by analysing the horizontal-to-vertical (H/V) spectral ratios of ambient vibration records. The H/V curves were derived using 31 d of continuous seismic data from 3 broad-band and 16 short-period ocean bottom seismometer (OBS) stations. Using the broad-band data, H/V ratios from 0.01 to 10 Hz were derived, but the ratios below 0.1 Hz frequencies were unusually large and temporally unstable. Characterization of seismic noise energy from ∼1 yr of seismic data of three broad-band OBSs revealed variable and elevated energy conditions below 0.1 Hz due to typical long-period oceanic noise; we link these observations with the unstable H/V ratios below this frequency. Therefore, H/V analysis was performed in the frequency range of 0.1–10 Hz for both broad-band and short-period OBSs to obtain subsurface VS profiles. For the forward calculation of the H/V ratios in the inversion process, we used the recently developed 'hvgeneralized' method, which is based on the diffuse field assumption, and accounts for the water layer on top of stratified media. Moreover, available prior geological and geophysical information was utilized during the inversion of the H/V curves. We found that subsurface VS ranged from approximately 30 m s−1 at the seabed to approximately 4900 m s−1 at 7000 m below the sea floor (mbsf). Starting with the best model candidate at each OBS location, the effect of the water layer on the H/V curve in the deep ocean was investigated by comparing synthetic H/V curves with and without the water layer. The synthetic H/V analysis revealed that the water layer had a significant effect on H/V amplitudes at higher frequencies (>1 Hz), whereas comparatively little effect was observed at lower frequencies (<1 Hz). This study provides an empirical basis for H/V analysis using OBS data to determine VS down to several kilometres of sedimentary sequences to the upper crust with high-resolution. [ABSTRACT FROM AUTHOR]

We present a new compilation and analysis of broad-band ocean bottom seismometer noise properties from 15 yr of seismic deployments. We compile a comprehensive data set of representative four-component (seismometer and pressure gauge) noise spectra and cross-spectral properties (coherence, phase and admittance) for 551 unique stations spanning 18 U.S.-led experiments. This is matched with a comprehensive compilation of metadata parameters related to instrumentation and environmental properties for each station. We systematically investigate the similarity of noise spectra by grouping them according to these metadata parameters to determine which factors are the most important in determining noise characteristics. We find evidence for improvements in similarity of noise properties when grouped across parameters, with groupings by seismometer type and deployment water depth yielding the most significant and interpretable results. Instrument design, that is the entire deployed package, also plays an important role, although it strongly covaries with seismometer and water depth. We assess the presence of traditional sources of tilt, compliance, and microseismic noise to characterize their relative role across a variety of commonly used seismic frequency bands. We find that the presence of tilt noise is primarily dependent on the type of seismometer used (covariant with a particular subset of instrument design), that compliance noise follows anticipated relationships with water depth, and that shallow, oceanic shelf environments have systematically different microseism noise properties (which are, in turn, different from instruments deployed in shallow lake environments). These observations have important implications for the viability of commonly used seismic analysis techniques. Finally, we compare spectra and coherences before and after vertical channel tilt and compliance noise removal to evaluate the efficacy and limitations of these now standard processing techniques. These findings may assist in future experiment planning and instrument development, and our newly compiled noise data set serves as a building block for more targeted future investigations by the marine seismology community. [ABSTRACT FROM AUTHOR]

Secondary microseisms are the most energetic noise in continuous seismometer recordings. They are generated by interactions between ocean waves, including between gravity waves incident on and reflected from the coast. Coastal reflections of ocean waves leading to coastal microseismic sources are hard to estimate in various global numerical wave models, and independent quantification of these coastal sources through direct measurements can greatly improve these models. Here, we exploit a 41-km-long submarine optical fibre (OF) cable located offshore Toulon, France, using distributed acoustic sensing (DAS). We record both the amplitude and frequency of seafloor strains induced by ocean surface gravity waves, as well as secondary microseisms caused by the interaction of gravity waves incident and reflected from the coast. By leveraging the spatially distributed nature of DAS measurements, additional fundamental information is recovered such as the velocity and azimuth of the waves. We find that on average 30 per cent of the gravity waves are reflected at the coast generating local sources of secondary microseisms that manifest as Scholte waves. These local sources represent the most energetic contribution to the seismic noise recorded along the OF and by an onshore broad-band station located near the DAS interrogator. Furthermore, we estimate a coastal reflection coefficient of ocean surface gravity waves R 2 of about 0.07, which provides improved constraints for seismic noise generation models. In addition, we show that new local sources of microseisms can be generated when gravity waves characteristics (azimuth and frequency content) change and lead to some delays between the OF cable and buoy recordings. These analyses pave the way for a wide use of DAS data to monitor ocean–solid earth interactions as they provide a wealth of information on the reflection of gravity waves, coastal microseismic sources, and new constraints for numerical models of microseismic noise. [ABSTRACT FROM AUTHOR]

Air-to-ground coupled seismic arrivals, resulting from infrasound interacting with the ground surface, can be identified on seismometers co-located with microbarometers. Dense networks of co-located sensors provide detailed spatial information about the infrasonic wavefield, allowing investigation of infrasonic propagation and an assessment of the variability in air-to-ground coupling mechanisms. Signals generated by an explosion at Camp Minden, Louisiana, United States, on 2012 October 16 were recorded across the eastern United States on the USArray Transportable Array. 233 infrasound signals were identified at distances of up to 1566 km from the source, with 105 of these sites, at distances of up to 1147 km from the source, recording above-noise air-to-ground coupled seismic arrivals in the 2–4 Hz passband (chosen to maximise seismic arrival signal-to-noise ratios). The spatially dense infrasound recordings allowed a unique investigation of the azimuthal variation in infrasound amplitudes, showing that stratospheric arrival amplitudes increased by a factor greater than 4 between the edge and centre of the stratospheric propagation duct. Air-to-ground coupling coefficients, calculated as the ratios of temporally coincident peak-to-trough seismic and infrasound time-domain amplitudes, span over two orders of magnitude ([9.0 × 10−8, 5.0 × 10−5] ms−1 Pa−1). Sites exhibiting high coupling coefficients are predominantly located on alluvial sediments that support the generation of air-to-ground coupled Rayleigh waves, resulting in poor coherence between the seismic and infrasonic waveforms. In contrast, sites exhibiting low coupling coefficients are predominantly located on chemically weathered bedrock and only support direct coupling, leading to high coherence between air-to-ground coupled seismic and infrasonic recordings. The Camp Minden explosion observations, and a consideration of expected seismic noise and air-to-ground coupled arrival amplitudes, suggest that co-located pairs of seismometers and microbarometers are likely to record signals from small (tens of tons) explosions at distances of over 1000 km under favourable propagation conditions. Co-located deployments may therefore provide a cost-effective method for augmenting current state-of-the-art infrasound array networks, by providing a higher spatial density of recordings; this may assist in correctly associating recorded signals with, and estimating the location of, source events. [ABSTRACT FROM AUTHOR]

Tilting of the ground due to loading by the variable atmosphere is known to corrupt very long period horizontal seismic records (below 10 mHz) even at the quietest stations. At BFO (Black Forest Observatory, SW-Germany), the opportunity arose to study these disturbances on a variety of simultaneously operated state-of-the-art broad-band sensors. A series of time windows with clear atmospherically caused effects was selected and attempts were made to model these 'signals' in a deterministic way. This was done by simultaneously least-squares fitting the locally recorded barometric pressure and its Hilbert transform to the ground accelerations in a bandpass between 100 and 3600 s periods. Variance reductions of up to 97 per cent were obtained. We show our results by combining the 'specific pressure induced accelerations' for the two horizontal components of the same sensor as vectors on a horizontal plane, one for direct pressure and one for its Hilbert transform. It turned out that at BFO the direct pressure effects are large, strongly position dependent and largely independent of atmospheric events for instruments installed on piers, while three post-hole sensors are only slightly affected. The infamous 'cavity effects' are invoked to be responsible for these large effects on the pier sensors. On the other hand, in the majority of cases all sensors showed very similar magnitudes and directions for the vectors obtained for the regression with the Hilbert transform, but highly variable from event to event especially in direction. Therefore, this direction most certainly has to do with the gradient of the pressure field moving over the station which causes a larger scale deformation of the crust. The observations are very consistent with these two fundamental mechanisms of how fluctuations of atmospheric surface pressure causes tilt noise. The results provide a sound basis for further improvements of the models for these mechanisms. The methods used here can already help to reduce atmospherically induced noise in long-period horizontal seismic records. [ABSTRACT FROM AUTHOR]

In late June 2016, the Harmalière clayey landslide (located 30 km south of the city of Grenoble, French Alps) was dramatically reactivated at the headscarp after a 35-yr-long period of continuous but limited activity. The total involved volume, which moved as sliding blocks of various sizes, was estimated to be about 2 × 10 6 m3. Two seismometers were installed at the rear of the main headscarp in August 2016, on both sides of a developing fracture delineating a block with a volume of a few hundred cubic metres. For 4 months, they continuously recorded seismic ambient vibrations and microearthquakes until the block broke. Five seismic parameters were derived from the monitoring: the cumulative number of microearthquakes (CNe), the seismic energy (SE), the block resonance frequency (fB), the relative variation in Rayleigh wave velocity (d V/V) deduced from noise cross-correlations between the two sensors and the associated correlation coefficient (CC). All parameters showed a significant precursory signal before the rupture, but at very different times, which indicates the complexity of the rupture mechanism in this clay material. [ABSTRACT FROM AUTHOR]

Station noise levels play a fundamental limitation in our ability to detect seismic signals. These noise levels are frequency-dependent and arise from a number of physically different drivers. At periods greater than 100 s, station noise levels are often limited by the self-noise of the instrument as well as the sensitivity of the instrument to non-seismic noise sources. Recently, station operators in the Global Seismographic Network (GSN) have deployed several Streckeisen STS-6A very broad-band borehole seismometers. These sensors provide a potential replacement for the no-longer-produced Streckeisen STS-1 seismometer and the GeoTech KS-54 000 borehole seismometer. Along with showing some of the initial observational improvements from installing modern very broad-band seismometers at depth, we look at current limitations in the seismic resolution from earth tide periods 100 000 s (0.01 mHz) to Nyquist at most GSN sites (0.02 s or 50 Hz). Finally, we show the potential for improved observations of continuously excited horizontal Earth hum as well as the splitting of very long-period torsional modes. Both of these observations make use of the low horizontal noise levels which are obtained by installing very broad-band borehole seismometers at depth. [ABSTRACT FROM AUTHOR]

We use seismic noise cross-correlations to obtain a 3-D tomography model of SV -wave velocities beneath the western Indian Ocean, in the depth range of the oceanic crust and uppermost mantle. The study area covers 2000 × 2000 km2 between Madagascar and the three spreading ridges of the Indian Ocean, centred on the volcanic hotspot of La Réunion. We use seismograms from 38 ocean bottom seismometers (OBSs) deployed by the RHUM-RUM project and 10 island stations on La Réunion, Madagascar, Mauritius, Rodrigues, and Tromelin. Phase cross-correlations are calculated for 1119 OBS-to-OBS, land-to-OBS, and land-to-land station pairs, and a phase-weighted stacking algorithm yields robust group velocity measurements in the period range of 3–50 s. We demonstrate that OBS correlations across large interstation distances of >2000 km are of sufficiently high quality for large-scale tomography of ocean basins. Many OBSs yielded similarly good group velocity measurements as land stations. Besides Rayleigh waves, the noise correlations contain a low-velocity wave type propagating at 0.8–1.5 km s−1 over distances exceeding 1000 km, presumably Scholte waves travelling through seafloor sediments. The 100 highest-quality group velocity curves are selected for tomographic inversion at crustal and lithospheric depths. The inversion is executed jointly with a data set of longer-period, Rayleigh-wave phase and group velocity measurements from earthquakes, which had previously yielded a 3-D model of Indian Ocean lithosphere and asthenosphere. Robust resolution tests and plausible structural findings in the upper 30 km validate the use of noise-derived OBS correlations for adding crustal structure to earthquake-derived tomography of the oceanic mantle. Relative to crustal reference model CRUST1.0, our new shear-velocity model tends to enhance both slow and fast anomalies. It reveals slow anomalies at 20 km depth beneath La Réunion, Mauritius, Rodrigues Ridge, Madagascar Rise, and beneath the Central Indian spreading ridge. These structures can clearly be associated with increased crustal thickness and/or volcanic activity. Locally thickened crust beneath La Réunion and Mauritius is probably related to magmatic underplating by the hotspot. In addition, these islands are characterized by a thickened lithosphere that may reflect the depleted, dehydrated mantle regions from which the crustal melts where sourced. Our tomography model is available as electronic supplement. [ABSTRACT FROM AUTHOR]

We propose a new Bayesian method to reveal the Vs structure of the near surface of the earth using spatial autocorrelation (SPAC) functions and apply this new method to synthetic, broad-band and geophone data sets. The principle of SPAC is introduced, and an implementation of the Bayesian Monte Carlo inversion (BMCI) for modelling SPAC coherency functions is described. To demonstrate its effectiveness, BMCI is applied to synthetic tests, data from 14 SPAC array sites in the Salt Lake Valley (SLV), Utah, and two arrays (one broad-band and one geophone) located in south central Utah. The Vs models derived from previous SPAC analysis of the 14 SLV sites differ by 10 per cent at most from those determined by BMCI and lie within uncertainties determined for the BMCI models. These agreements demonstrate the effectiveness of the BMCI method. The synthetic tests and applications to the SLV SPAC data show BMCI has great potential to resolve Vs structure down to at least 400 m. To achieve resolution for deeper Vs structure, longer duration deployments, wider array apertures, and additional seismometers or geophones can be used. Additionally, when the target frequencies are greater than 0.1 Hz, there is no apparent disadvantage in using geophone data for BMCI compared to broad-band data. Most significantly, BMCI places a quantifiable constraint on the uncertainties of the Vs models as well as Vs30. [ABSTRACT FROM AUTHOR]

We formulate a full wavefield inversion (FWFI) for ambient seismic noise recorded by large and dense seismograph arrays. FWFI exploits the constraints on the gradients of the wavefield that array data inherently possess. We pose FWFI as a partial differential equation (PDE) constrained inverse problem resulting in a joint estimation of a reconstructed wavefield and the medium parameters. The inverse problem is solved by variable projection. We explicitly allow for non-unique solutions to the PDE system that is imposed as a constraint. The boundary conditions of the wavefield do not need to be specified, and can remain unknown. This makes the algorithm suitable for inverting observations of ambient seismic noise by dense arrays. The result is that the inverse problem for subsurface properties becomes insensitive to the character and distribution of the noise sources that excited the seismic wavefield. In principle the formulation holds equally for ambient noise wavefields and for wavefields excited by controlled sources. The theory is supported with examples in one dimension in the time domain, and in two dimensions in the frequency domain. The latter are of interest in the inversion of surface wave ambient noise for phase velocity maps. [ABSTRACT FROM AUTHOR]

Many applications in seismology rely on the accurate absolute timing of seismograms. However, both seismological land stations and ocean bottom seismometers (OBSs) can be affected by clock errors, which cause the absolute timing of seismograms to deviate from a highly accurate reference time signal, usually provided by GPS satellites. Timing problems can occur in land stations when synchronization with a GPS signal is temporarily or permanently lost. This can give rise to complicated, time-dependent clock drifts relative to GPS time, due to varying environmental conditions. Seismometers at the ocean bottom cannot receive GPS satellite signals, but operate in more stable ambient conditions than land stations. The standard protocol is to synchronize an OBS with a GPS signal immediately before deployment and after recovery. The measured timing deviation, called ‘skew’, is assumed to have accumulated linearly over the deployment interval, an assumption that is plausible but usually not verifiable. In recent years, cross-correlations of ambient microseismic noise have been put to use for correcting timing errors, but have been limited to interstation distances of at most a few tens of kilometres without reducing the temporal resolution. We apply noise cross-correlations to the evaluation of clock errors in four broad-band land stations and 53 wideband and broad-band OBSs, which were installed on and around the island of La Réunion in the western Indian Ocean during the RHUM-RUM (Réunion Hotspot and Upper Mantle-Réunions Unterer Mantel) experiment. We correlate all three seismic components, plus a hydrophone channel in OBS stations. Daily cross-correlation functions are derived for intermediate distances (∼20 km) for land-to-land station pairs; stable, 10 d stacks are obtained for very large interstation distances up to >300 km for land-to-OBS and OBS-to-OBS configurations. Averaging over multiple station pairs, and up to 16 component pairs per station, improves the accuracy of the method by a factor of four compared to the single-channel approaches of prior studies. The timing accuracy of our method is estimated to be ∼20 ms standard deviation or one sample at a sampling rate of 50 Hz. In land stations, nonlinear clock drifts and clock jumps of up to 6 min are detected and successfully corrected. For 52 out of 53 OBSs, we successfully obtain drift functions over time, which validate the common assumption of linear clock drift. Skew values that were available for 29 of these OBSs are consistent with our independent estimates within their observational error bars. For 23 OBSs that lacked skew measurements, linear OBS clock drifts range between 0.2 and 8.8 ms d−1. In addition to linear drift, three OBSs are affected by clock jumps of ∼1 s, probably indicating a missing sample problem that would otherwise have gone undetected. Thus we demonstrate the routine feasibility of high-accuracy clock corrections in land and OBSs over a wide range of interstation distances. [ABSTRACT FROM AUTHOR]

In this paper, an underground experiment at the Äspö Hard Rock Laboratory (HRL) is described. Main goal is optimizing geothermal heat exchange in crystalline rock mass at depth by multistage hydraulic fracturing with minimal impact on the environment, that is, seismic events. For this, three arrays with acoustic emission, microseismicity and electromagnetic sensors are installed mapping hydraulic fracture initiation and growth. Fractures are driven by three different water injection schemes (continuous, progressive and pulse pressurization). After a brief review of hydraulic fracture operations in crystalline rock mass at mine scale, the site geology and the stress conditions at Äspö HRL are described. Then, the continuous, single-flow rate and alternative, multiple-flow rate fracture breakdown tests in a horizontal borehole at depth level 410 m are described together with the monitoring networks and sensitivity. Monitoring results include the primary catalogue of acoustic emission hypocentres obtained from four hydraulic fractures with the in situ trigger and localizing network. The continuous versus alternative water injection schemes are discussed in terms of the fracture breakdown pressure, the fracture pattern from impression packer result and the monitoring at the arrays. An example of multistage hydraulic fracturing with several phases of opening and closing of fracture walls is evaluated using data from acoustic emissions, seismic broad-band recordings and electromagnetic signal response. Based on our limited amount of in situ tests (six) and evaluation of three tests in Ävrö granodiorite, in the multiple-flow rate test with progressively increasing target pressure, the acoustic emission activity starts at a later stage in the fracturing process compared to the conventional fracturing case with continuous water injection. In tendency, also the total number and magnitude of acoustic events are found to be smaller in the progressive treatment with frequent phases of depressurization. [ABSTRACT FROM AUTHOR]

Obtaining new seismic responses from existing recordings is generally referred to as seismic interferometry (SI). Conventionally, the SI responses are retrieved by simple crosscorrelation of recordings made by separate receivers: one of the receivers acts as a 'virtual source' whose response is retrieved at the other receivers. When SI is applied to recordings of ambient seismic noise, mostly surface waves are retrieved. The newly retrieved surface wave responses can be used to extract receiver-receiver phase velocities. These phase velocities often serve as input parameters for tomographic inverse problems. Another application of SI exploits the temporal stability of the multiply scattered arrivals of the newly retrieved surface wave responses. Temporal variations in the stability and/or arrival time of these multiply scattered arrivals can often be linked to temporally varying parameters such as hydrocarbon production and precipitation. For all applications, however, the accuracy of the retrieved responses is paramount. Correct response retrieval relies on a uniform illumination of the receivers: irregularities in the illumination pattern degrade the accuracy of the newly retrieved responses. In practice, the illumination pattern is often far from uniform. In that case, simple crosscorrelation of separate receiver recordings only yields an estimate of the actual, correct virtual-source response. Reformulating the theory underlying SI by crosscorrelation as a multidimensional deconvolution (MDD) process, allows this estimate to be improved. SI by MDD corrects for the non-uniform illumination pattern by means of a so-called point-spread function (PSF), which captures the irregularities in the illumination pattern. Deconvolution by this PSF removes the imprint of the irregularities on the responses obtained through simple crosscorrelation. We apply SI by MDD to surface wave data recorded by the Malargue seismic array in western Argentina. The aperture of the array is approximately 60 km and it is located on a plateau just east of the Andean mountain range. The array has a T-shape: the receivers along one of the two lines act as virtual sources whose responses are recorded by the receivers along the other (perpendicular) line. We select time windows dominated by surface wave noise travelling in a favourable direction, that is, traversing the line of virtual sources before arriving at the receivers at which we aim to retrieve the virtual-source responses. These time windows are selected through a frequency-dependent slowness analysis along the two receiver lines. From the selected time windows, estimates of virtual-source responses are retrieved by means of crosscorrelations. Similarly, crosscorrelations between the positions of the virtual sources are computed to build the PSF. We use the PSF to deconvolve the effect of illumination irregularities and the source function from the virtual-source responses retrieved by crosscorrelation. The combined effect of time-window selection and MDD results in more accurate and temporally stable surface wave responses. [ABSTRACT FROM AUTHOR]

Ocean wave activity excites seismic waves that propagate through the solid earth, known as microseisms, which, once recorded on oceanic islands, can be used to analyse the swell. Here, we analyse the microseismic noise recorded in different period ranges by the permanent seismic station RER on La Réunion Island and by a temporary network of 10 broad-band seismic stations deployed on the island to analyse extreme swell events. We perform a comparative analysis of cyclonic and austral swell events by analysing not only the primary (PM, ~10-20 s period) and secondary (SM, ~3-10 s) microseisms but also the long-period secondary microseisms (LPSMs, ~ 7-10 s), which may result from the interaction between incident ocean waves and the reflected waves off the coast. We compare the microseismic observations with buoy data when available and with hindcasts from numerical ocean wave models. We show that each cyclone is characterized by its own individual signature in the SM, which depends not only on its distance and intensity but also on its dynamics and trajectory. Thus, the SMcontains relevant information for cyclone detection and monitoring. Analysing the PM and the LPSM, and comparing it to direct buoy observations and/or wave numerical models allows characterizing the local impact of the swell with the island in terms of amplitude, period, and sometimes, direction of propagation, making possible to use a seismic station as an ocean wave gauge. The microseisms, which link the atmosphere, the ocean and the solid Earth, can thus provide valuable observations on extreme swell events, in addition to oceanic and meteorological data. [ABSTRACT FROM AUTHOR]

Microseisms in the period of 2-10 s are generated in deep oceans and near coastal regions. It is common for microseisms from multiple sources to arrive at the same time at a given seismometer. It is therefore desirable to be able to measure multiple slowness vectors accurately. Popular ways to estimate the direction of arrival of ocean induced microseisms are the conventional (fk) or adaptive (Capon) beamformer. These techniques give robust estimates, but are limited in their resolution capabilities and hence do not always detect all arrivals. One of the limiting factors in determining direction of arrival with seismic arrays is the array response, which can strongly influence the estimation of weaker sources. In this work, we aim to improve the resolution for weaker sources and evaluate the performance of two deconvolution algorithms, Richardson-Lucy deconvolution and a new implementation of CLEAN-PSF. The algorithms are tested with three arrays of different aperture (ASAR, WRA and NORSAR) using 1 month of real data each and compared with the conventional approaches. We find an improvement over conventional methods from both algorithms and the best performance with CLEAN-PSF. We then extend the CLEAN-PSF framework to three components (3C) and evaluate 1 yr of data from the Pilbara Seismic Array in northwest Australia. The 3C CLEAN-PSF analysis is capable in resolving a previously undetected Sn phase. [ABSTRACT FROM AUTHOR]

Inhomogeneous noise sources surrounding stations produce asymmetric amplitudes in crosscorrelation functions that yield preferential source directions. Here we show that preprocessing biases the dominant source direction estimate towards the source producing long-duration signals by down-weighting high-amplitude signals. Tests with both synthetic data and observations show that conventional preprocessing, where only earthquakes and local transients (e.g. trawling, fish impacts) are removed, is more sensitive to coherent energy, while one-bit preprocessing and running-absolute-mean preprocessing are more influenced by signal duration. Comparisons between different preprocessing methods are made on data from the Cascadia Initiative ocean bottom seismometer array, where we find that the total energy arriving from pelagic and coastal areas is similar. Moreover, pelagic-generated signals tend to be weaker but have longer duration, in contrast to coastal-generated signals that tend to be stronger but have shorter duration. [ABSTRACT FROM AUTHOR]

Ocean waves activity is a major source of microvibrations that travel through the solid Earth, known as microseismic noise and recorded worldwide by broadband seismometers. Analysis of microseismic noise in continuous seismic records can be used to investigate noise sources in the oceans such as storms, and their variations in space and time, making possible the regional and global-scale monitoring of the wave climate. In order to complete the knowledge of the Atlantic and Pacific oceans microseismic noise sources, we analyse 1 yr of continuous data recorded by permanent seismic stations located in the Indian Ocean basin.We primarily focus on secondary microseisms (SM) that are dominated by Rayleigh waves between 6 and 11 s of period. Continuous polarization analyses in this frequency band at 15 individual seismic stations allow us to quantify the number of polarized signal corresponding to Rayleigh waves, and to retrieve their backazimuths (BAZ) in the time-frequency domain. We observe clear seasonal variations in the number of polarized signals and in their frequencies, but not in their BAZ that consistently point towards the Southern part of the basin throughout the year. This property is very peculiar to the Indian Ocean that is closed on its Northern side, and therefore not affected by large ocean storms during Northern Hemisphere winters. We show that the noise amplitude seasonal variations and the backazimuth directions are consistent with the source areas computed from ocean wave models. [ABSTRACT FROM AUTHOR]

Based on cross-correlations of ambient seismic noise computed using 61 ocean bottom seismometers (OBSs) within the Juan de Fuca (JdF) plate from the Cascadia Initiative experiment and 42 continental stations near the coast of the western United States, we investigate the locations of generation of the primary (11-20 s period) and secondary (5-10 s period) microseisms in the northern Pacific Ocean by analysing the directionality and seasonality of the microseism (Rayleigh wave) signals received in this region. We conclude that (1) the ambient noise observed across the array is much different in the primary and secondary microseism bands, both in its azimuthal content and seasonal variation. (2) The principal secondary microseism signals propagate towards the east, consistent with their generation in deep waters of the North Pacific, perhaps coincident both with the region of observed body wave excitation and the predicted wave--wave interaction region from recent studies. (3) The primary microseism, as indicated by observations of the azimuthal dependence of the fundamental mode Rayleigh wave as well as observations of precursory arrivals, derives significantly from the shallow waters of the eastern Pacific near to the JdF plate but also has a component generated at greater distance of unknown origin. (4) These observations suggest different physical mechanisms for generating the two microseisms: the secondary microseisms are likely to be generated by non-linear wave--wave interaction over the deep Pacific Ocean, while the primary microseism may couple directly into the solid earth locally in shallow waters from ocean gravity waves. (5) Above 5 s period, high quality empirical Green's functions are observed from cross-correlations between deep water OBSs and continental stations, which illustrates that microseisms propagate efficiently from either deep or shallow water source regions onto the continent and are well recorded by continental seismic stations. [ABSTRACT FROM AUTHOR]

We analysed background surface waves in seismic ambient noise by cross-correlating continuous records of eight ocean bottom seismometers and nine differential pressure gauges deployed in the northwestern Pacific Ocean by the PLATE project. After estimating the clock delay and instrumental phase responses of differential pressure gauges by using cross-correlation functions, we measured average phase velocities in the area of the array for the fundamental-, first higher- and second higher-mode Rayleigh waves, and the fundamental-mode Love waves at a period range of 3–40 s by waveform fitting. We then measured phase-velocity anomalies of fundamental-mode and first higher-mode Rayleigh waves for each pair of stations at a period range of 5–25 s, and corrected the effect of variation in water-depths. The seismic anomalies imply the presence of strong azimuthal anisotropy beneath the eastern part of array. The direction of maximum velocity is approximately N35°E in the fossil seafloor spreading direction perpendicular to magnetic lineations from the ancient triple junction at this location. The peak-to-peak intensity of shear-wave velocity anisotropy in the mantle is ∼7 per cent. [ABSTRACT FROM PUBLISHER]

Seismic noise, generated by oceanic microseisms and other sources, illuminates the crust in a manner different from tectonic sources, and therefore provides independent information. The primary measurable is the two-point cross-correlation, evaluated using traces recorded at a pair of seismometers over a finite-time interval. However, raw seismic traces contain intermittent large-amplitude perturbations arising from tectonic activity and instrumental errors, which may corrupt the estimated cross-correlations of microseismic fluctuations. To diminish the impact of these perturbations, the recorded traces are filtered using the non-linear one-bit digitizer, which replaces the measurement by its sign. Previous theory shows that for stationary Gaussian-distributed seismic noise fluctuations one-bit and raw correlation functions are related by a simple invertible transformation. Here we extend this to show that the simple correspondence between these two correlation techniques remains valid for non-stationary Gaussian and a very broad range of non-Gaussian processes as well. For a limited range of stationary and non-stationary Gaussian fluctuations, we find that one-bit filtering performs at least as well as spectral whitening. We therefore recommend using one-bit filtering when processing terrestrial seismic noise, with the substantial benefit that the measurements are fully compatible with current theoretical interpretation (e.g. adjoint theory). Given that seismic records are non-stationary and comprise small-amplitude fluctuations and intermittent, large-amplitude tectonic/other perturbations, we outline an algorithm to accurately retrieve the correlation function of the small-amplitude signals. [ABSTRACT FROM AUTHOR]

In this study, the crustal and uppermost mantle shear wave velocities beneath the Japanese islands have been determined by inversion from seismic ambient noise tomography using data recorded at 75 Full Range Seismograph Network of Japan broad-band seismic stations, which are uniformly distributed across the Japanese islands. By cross-correlating 2 yr of vertical component seismic ambient noise recordings, we are able to extract Rayleigh wave empirical Green's functions, which are subsequently used to measure phase velocity dispersion in the period band of 6–50 s. The dispersion data are then inverted to yield 2-D tomographic phase velocity maps and 3-D shear wave velocity models. Our results show that the velocity variations at short periods (∼10 s), or in the uppermost crust, correlate well with the major known surface geological and tectonic features. In particular, the distribution of low-velocity anomalies shows good spatial correlation with active faults, volcanoes and terrains of sediment exposure, whereas the high-velocity anomalies are mainly associated with the mountain ranges. We also observe that large upper crustal earthquakes (5.0 ≤ M ≤ 8.0, depth ≤ 25 km) mainly occurred in low-velocity anomalies or along the boundary between low- and high-velocity anomalies, suggesting that large upper crustal earthquakes do not strike randomly or uniformly; rather they are inclined to nucleate within or adjacent to low-velocity areas. [ABSTRACT FROM AUTHOR]

We present an approach for seismic noise which improves Rayleigh-wave ellipticity estimates by reducing the influences from body waves and Love waves. The method requires three-component seismographs and uses the phase-shift information between the vertical and horizontal components. We select data that show the phase-shift of about 90° in order to separate Rayleigh waves from body waves and Love waves. In comparison to the usual H/V approach, the estimates for Rayleigh-wave ellipticity are systematically different by about 20 per cent because the existence of S waves and Love waves systematically increases horizontal amplitudes. The differences in the inverted S-wave velocity structure reach up to about 20 per cent in the upper 20 km. The method is limited to a relatively low-frequency range, below 0.3 Hz, as the 90° phase-shift peaks become difficult to identify above 0.3 Hz. Despite this limitation in the frequency range, this approach can improve the S-wave velocity structure in the upper 5–10 km, which should lead to a better prediction of long-period ground motions for periods of about 5–10 s. An improved prediction of such long-period ground motion is important for hazard mitigation in the world's metropolitan areas with high-rise buildings. This approach is applied to the structure in the broader Los Angeles basin in Southern California and improvement of the reference velocity structure in the upper 5–10 km is demonstrated. [ABSTRACT FROM AUTHOR]

SUMMARY The accuracy with which geophysical observations are made is inherently determined by the geometry of the observation network, and typically depends on a highly non-linear relationship between data and earth parameters. Statistical experimental design provides a means of optimizing the network geometry to provide maximum information about parameters of interest. Here, we re-derive the nonlinear experimental design D N optimization method, without the need for the usual assumption of a multivariate normal model of data uncertainties. We demonstrate the criterion's utility by applying it to the problem of seismic network expansion in the active Kawerau geothermal field, Taupo Volcanic Zone, New Zealand. The design calculations maximize the ratio of the hypocentre data generalized variance (attributable to resolvable spatial separation of earthquakes) to the measurement error generalized variance (attributable to observational uncertainties), and incorporate realistic 3-D velocity and attenuation models, surface noise sources, and both P- and S-wave data. In geologically complex areas, statistical experimental design provides a means of objectively deploying finite observational resources to target areas of particular interest while taking into account environmental and logistical factors. [ABSTRACT FROM AUTHOR]

SUMMARY A limitation of most forms of passive seismic tomography using distant earthquakes lies in the fact that crustal structure is poorly resolved. An attempt is made here to address this issue by modelling teleseismic receiver functions (RFs) and dispersion curves derived from ambient noise through a multistep approach. The SEAL3 experiment in central and southern New South Wales (NSW) used here, represents one of 13 array deployments that so far comprise the large WOMBAT project, which aims to cover a significant portion of the Australian continent with a rolling array of seismometers. An interactive, forward-modelling software package (IRFFM2) is introduced and applied to the observed RFs and surface wave dispersion curves to define a prior, physically acceptable range of elastic parameters in the lithosphere, which is combined with a grid-search and a linearized inversion. Our results emphasize the importance of a joint treatment of RFs and dispersion data as the predictions from 1-D velocity models at individual stations derived from only RFs display large departures from the observed ambient noise dispersion curves. In total, 27 jointly constrained 1-D shear wave models are produced, which provide sufficient sampling of the crust beneath SEAL3 to permit detailed inferences about lateral variations in structure to be made. Of particular note is the observation that the Moho deepens towards the mountainous southeast, where it exceeds 50 km in depth beneath the Southern Highlands of NSW, thus marking out some of the thickest crust in Australia. The complex lateral variations in mid-lower crustal velocity that we observe probably reflect the manifold interactions of a thinning lithosphere, associated igneous underplating, recent hot-spot-related volcanism and uplift. Our results image an important part of the lithosphere that is poorly constrained by regional and teleseismic tomography, and contribute to the understanding of the formation of the southern highlands and the Palaeozoic Lachlan Orogen. [ABSTRACT FROM AUTHOR]

In this work, we analyse continuous measurements of microseisms to assess the reliability of the fundamental resonance frequency estimated by means of the horizontal-to-vertical ( H/ V) spectral ratio within the 0.1–1 Hz frequency range, using short-period sensors (natural period of 1 s). We apply the H/ V technique to recordings of stations installed in two alluvial basins with different sedimentary cover thicknesses—the Lower Rhine Embayment (Germany) and the Gubbio Plain (Central Italy). The spectral ratios are estimated over the time–frequency domain, and we discuss the reliability of the results considering both the variability of the microseism activity and the amplitude of the instrumental noise. We show that microseisms measured by short period sensors allow the retrieval of fundamental resonance frequencies greater than about 0.1–0.2 Hz, with this lower frequency bound depending on the relative amplitude of the microseism signal and the self-noise of the instruments. In particular, we show an example where the considered short-period sensor is connected to instruments characterized by an instrumental noise level which allows detecting only fundamental frequencies greater than about 0.4 Hz. Since the frequency at which the peak of the H/ V spectral ratio is biased depends upon the seismic signal-to-instrument noise ratio, the power spectral amplitude of instrumental self-noise should be always considered when interpreting the frequency of the peak as the fundamental resonance frequency of the investigated site. [ABSTRACT FROM AUTHOR]

With the important increase in the number of instruments, especially in the near field of quite significant earthquakes, unsaturated traces have led seismologists to question what they are really measuring. We have performed a review of previous studies related to the effects of rotations on both horizontal and vertical components of various sensors. Illustrations of near-field records show that the recovering static displacement requires an accurate rotation estimation that present-day seismometers cannot achieve. Estimations of coseismic tilts (rotation around an horizontal axis) using seismometers cannot be achieved independently of translational motion. Therefore, for specific configurations such as near-field or long-period far-field, reconstruction of the ground motion requires specific rotation measurements. Effects of the Chi-Chi earthquake (1999) on accelerograms make the static displacement estimation unreliable. For long-period background noise, far-field horizontal seismic signals present unexpected N45° polarization, which is explained by the similar influence of the rotation around the vertical axis on the two horizontal components. For the GEOSCOPE network, this feature has been seen for stations with a Streckeisen STS-1 sensor, while those with a Streckeisen STS-2 sensor do not show this polarization feature. This study suggests that sensor installation should follow a protocol that will better guarantee the verticality of the sensors. Moreover, rotational recordings with ad hoc sensors are necessary, and by adequately correcting the traces, they will enable us to reconstruct the translational motion from recorded seismic signals. [ABSTRACT FROM AUTHOR]

Acceleration power spectral densities of vertical seismic noise at the best seismic stations show a minimum near 3 mHz. We suggest that this minimum is caused by a cancellation near this frequency of Newtonian attraction vs. free air and inertial effects exerted by atmospheric phenomena on the sensor mass. Simplistic models of atmospheric phenomena are used to quantify this effect and examples are shown for special atmospheric events. [ABSTRACT FROM AUTHOR]

A method for using PS converted waves to detect and characterize azimuthal seismic anisotropy is developed, and applied to four-component ocean-bottom seismometer data acquired near the base of the continental slope west of Svalbard. The method determines the distribution of coherent signal vibrating perpendicular to the sagittal plane (the transverse direction), and employs an optimized rotation technique to separate the effects of anisotropy from other contributions to the seismic amplitude. The magnitude of the PS-wave transverse amplitude, when normalized by its radial counterpart, is shown to correspond to the degree of ellipticity of particle motions and is used to constrain the anisotropic symmetry class in the near-surface sediments. Results indicate the presence of orthorhombic or horizontal, transversely isotropic symmetry that is consistent with the presence of near-vertical, aligned cracks. From its orientation, the observed anisotropy is associated with: (1) gravitational, downslope stresses in the sediment sequence, which dips by approximately two degrees to the SW and (2) tectonic stresses associated with the segmented spreading centres of the northern Knipovich Ridge, which spread in the NW–SE direction. Aligned cracks are likely to influence the local distribution of gas hydrate and free gas in these sediments. An investigation of the effect of seismic noise on optimized rotation shows that the vectorial combination of independent geophone records has variable sensitivity to source-related noise. This sensitivity peaks where the sagittal azimuth bisects the deployment orientation of the geophones on the seabed. A method to reduce the sensitivity to directional noise is developed and shown to improve the transverse amplitude plots. This approach to noise reduction may be used to improve the quality of seismic images derived from multicomponent data. [ABSTRACT FROM AUTHOR]

The horizontal to vertical spectral ratio (HVSR) technique of microtremor measurement is widely applied in microzonation and site response studies. Given the small amplitude of the measured signal, unwanted disturbances may adversely affect the results. We analysed the influence of local wind on microtremor measurements as follows: first, we studied variation in HVSR at a permanent, three-component seismological station under various wind patterns; second, we conducted an experiment under controlled condition in a wind room equipped with a laser particle image velocimeter, analysing the effect of increasing wind speed on various sensor/digitizer configurations; third, we undertook numerical modelling to explain the observed data. Our main conclusions are that wind increases the amplitude of all components of the microtremor wavefield, but does not affect HVSR; however, certain sensor/digitizer configurations can be adversely affected by instrumental noise due to turbulence at the air/instrument/ground interface. [ABSTRACT FROM AUTHOR]

Paul Denton describes an experiment to determine whether it's the band or the dancing that generates seismic noise at music festivals.Paul Denton describes an experiment to determine whether it's the band or the dancing that generates seismic noise at music festivals. [ABSTRACT FROM AUTHOR]