High-quality lowest-frequency normal mode strain observations at the Black Forest Observatory (SW-Germany) and comparison with horizontal broad-band seismometer data and synthetics

International audienceWe present spectra concentrating on the lowest-frequency normal modes of the Earth obtained from records of the invar-wire strainmeters and STS-1 broad-band seismometers located in the Black Forest Observatory, Germany after the disastrous earthquakes off the NW coast of Sumatra in 2004 and off the coast near Tohoku, Japan in 2011. We compare the spectra to ones obtained from synthetic seismograms computed using a mode summation technique for an anelastic, elliptical, rotating, spherically symmetric Earth model. The synthetics include strain– strain-coupling effects by using coupling coefficients obtained from comparisons between Earth tide signals recorded by the strainmeters and synthetic tidal records. We show that for the low-frequency toroidal and spheroidal modes up to 1 mHz, the strainmeters produce better signal-to-noise ratios than the broad-band horizontal seismometers. Overall, the comparison with the synthetics is satisfactory but not as good as for vertical accelerations. In particular, we demonstrate the high quality of the strainmeter data by showing the Coriolis splitting of toroidal modes for the first time in individual records, the first clear observation of the singlet 2 S 0 1 and the detection of the fundamental radial mode 0 S 0 with good signal-to-noise ratio and with a strain amplitude of 10 −11. We also identify the latter mode in a record of the Isabella strainmeter after the great Chilean quake in 1960, the detection of which was missed by the original studies

Normal mode simulation of prompt elastogravity signals induced by an earthquake rupture

As soon as an earthquake starts, the rupture and the propagation of seismic waves redistribute masses within the Earth. This mass redistribution generates in turn a long-range perturbation of the Earth gravitational field, which can be recorded before the arrival of the direct seismic waves. The recent first observations of such early signals motivate the use of the normal mode theory to model the elastogravity perturbations recorded by a ground-coupled seismometer or gravimeter. Complete modelling by normal mode summation is challenging due to the very large difference in amplitude between the prompt elastogravity signals and the direct P-wave signal. We overcome this problem by introducing a two-step simulation approach. The normal mode approach enables a fast computation of elastogravity signals in layered self-gravitating Earth models. The fast and accurate computation of gravity perturbations indicates instrument locations where signal detection may be achieved, and may prove useful in the implementation of a gravity-based earthquake early warning system

Normal mode simulation of prompt elastogravity signals induced by an earthquake rupture

As soon as an earthquake starts, the rupture and the propagation of seismic waves redistribute masses within the Earth. This mass redistribution generates in turn a long-range perturbation of the Earth gravitational field, which can be recorded before the arrival of the direct seismic waves. The recent first observations of such early signals motivate the use of the normal mode theory to model the elastogravity perturbations recorded by a ground-coupled seismometer or gravimeter. Complete modelling by normal mode summation is challenging due to the very large difference in amplitude between the prompt elastogravity signals and the direct P-wave signal. We overcome this problem by introducing a two-step simulation approach. The normal mode approach enables a fast computation of elastogravity signals in layered self-gravitating Earth models. The fast and accurate computation of gravity perturbations indicates instrument locations where signal detection may be achieved, and may prove useful in the implementation of a gravity-based earthquake early warning system

Time-reversal method and cross-correlation techniques by normal mode theory: a three-point problem

International audienceSince its beginning in acoustics, the Time-Reversal method (hereafter referred as TR) has been explored by different studies to locate and characterize seismic sources in elastic media. But few authors have proposed an analytical analysis of the method, especially in the case of an elastic medium and for a finite body such as the Earth. In this paper, we use a normal mode approach (for general 3-D case and degenerate modes in 1-D reference model) to investigate the convergence properties of the TR method. We first investigate a three-point problem, with two fixed points which are the source and the receiver and a third one corresponding to a changing observation point. We extend the problem of a single channel TR experiment to a multiple channel and multiple station TR experiment. We show as well how this problem relates to the retrieval of Green's function with a multiple source cross-correlation and also the differences between TR method and cross-correlation techniques. Since most of the noise sources are located close to the surface of the Earth, we show that the time derivative of the cross-correlation of long-period seismograms with multiple sources at the surface is different from the Green's function. Next, we show the importance of a correct surface-area weighting of the signal resent by the stations according to a Voronoi tessellation of the Earth surface. We use arguments based on the stationary phase approximation to argue that phase-information is more important than amplitude information for getting a good focusing in TR experiment. Finally, by using linear relationships between the time-reversed displacement (resp. strain wavefields) and the components of a vector force source (resp. a moment tensor source), we show how to retrieve force (or moment tensor components) of any long period tectonic or environmental sources by time reversal

Uncertainty of moment tensor determination from surface wave analysis for shallow earthquakes

International audienceWe treat the problem of the non-uniqueness of the moment tensor inversion for shallow earthquakes from long period surface wave data. We give an existence condition for double couples radiating the same long period surface waves as the deviatoric moment tensor obtained by such inversion. We describe the family of such double couples and show that they may provide better estimates of double couple mechanisms than the traditional “best double-couple” solution. We demonstrate the results of application to four large shallow earthquakes with different source signatures

Source models of great earthquakes from ultra low-frequency normal mode data

International audienceWe present a new earthquake source inversion technique based on normal mode data for the simultaneous determination of the rupture duration, length and moment tensor of large earthquakes with unilateral rupture. We use ultra low-frequency (f w⩾8.5): (i) Sumatra-Andaman (26th December 2004); (ii) Nias, Sumatra (28th March 2005); (iii) Bengkulu (12th September 2007); (iv) Tohoku, Japan (11th March 2011); (v) Maule, Chile (27th February 2010); and, (vi) the 24 May 2013 Mw 8.3 Okhotsk Sea, Russia, deep (607 km) event. While finite source inversions for rupture length, duration, magnitude and fault mechanism are possible for the Sumatra-Andaman and Tohoku events, for all the other events their lower magnitudes only allow stable point source inversions of mode multiplets. We obtain the first normal mode finite source model for the 2011 Tohoku earthquake, which yields a fault length of 461 km, a rupture duration of 151 s, and hence an average rupture velocity of 3.05 km/s, giving an independent confirmation of the compact nature of this event. For all the other earthquakes studied, our new source models agree well with previous studies. We do not find any unexplained systematic differences between our results and those in the literature, suggesting that for the wave frequencies considered, the moment magnitude and the fault mechanism of the earthquakes studied do not show a strong frequency dependence

Time-reversal method and cross-correlation techniques by normal mode theory: a three-point problem

International audienceSince its beginning in acoustics, the Time-Reversal method (hereafter referred as TR) has been explored by different studies to locate and characterize seismic sources in elastic media. But few authors have proposed an analytical analysis of the method, especially in the case of an elastic medium and for a finite body such as the Earth. In this paper, we use a normal mode approach (for general 3-D case and degenerate modes in 1-D reference model) to investigate the convergence properties of the TR method. We first investigate a three-point problem, with two fixed points which are the source and the receiver and a third one corresponding to a changing observation point. We extend the problem of a single channel TR experiment to a multiple channel and multiple station TR experiment. We show as well how this problem relates to the retrieval of Green's function with a multiple source cross-correlation and also the differences between TR method and cross-correlation techniques. Since most of the noise sources are located close to the surface of the Earth, we show that the time derivative of the cross-correlation of long-period seismograms with multiple sources at the surface is different from the Green's function. Next, we show the importance of a correct surface-area weighting of the signal resent by the stations according to a Voronoi tessellation of the Earth surface. We use arguments based on the stationary phase approximation to argue that phase-information is more important than amplitude information for getting a good focusing in TR experiment. Finally, by using linear relationships between the time-reversed displacement (resp. strain wavefields) and the components of a vector force source (resp. a moment tensor source), we show how to retrieve force (or moment tensor components) of any long period tectonic or environmental sources by time reversal

Imaging and modeling the ionospheric airglow response over Hawaii to the tsunami generated by the Tohoku earthquake of 11 March 2011

International audienceAlthough only centimeters in amplitude over the open ocean, tsunamis can generate appreciable wave amplitudes in the upper atmosphere, including the naturally occurring chemiluminescent airglow layers, due to the exponential decrease in density with altitude. Here, we present the first observation of the airglow tsunami signature, resulting from the 11 March 2011 Tohoku earthquake off the eastern coast of Japan. These images are taken using a wide‐angle camera system located at the top of the Haleakala Volcano on Maui, Hawaii. They are correlated with GPS measurements of the total electron content from Hawaii GPS stations and the Jason‐1 satellite. We find waves propagating in the airglow layer from the direction of the earthquake epicenter with a velocity that matches that of the ocean tsunami. The first ionospheric signature precedes the modeled ocean tsunami generated by the main shock by approximately one hour. These results demonstrate the utility of monitoring the Earth's airglow layers for tsunami detection and early warning

Tarapacá intermediate-depth earthquake (Mw 7.7, 2005, northern Chile): A slab-pull event with horizontal fault plane constrained from seismologic and geodetic observations

International audience[1] A large (Mw 7.7) intermediate-depth earthquake occurred on 13 June 2005 in the Tarapacá region of the northern Chile seismic gap. Source parameters are inferred from teleseismic broadbands, strong motions, GPS and InSAR data. Relocated hypocenter is found at $98 km depth within the subducting slab. The 21-days aftershock distribution, constrained by a postseismic temporary array, indicates a sub-horizontal fault plane lying between the planes of the double seismic zone and an upper bound of the rupture area of 60 km  30 km. Teleseismic inversion shows a slab-pull down dip extension mechanism on a nearly horizontal plane. Total seismic and geodetic moments are$5.5  10 20 N.m, with an averaged slip of 6.5 m from geodesy. The earthquake rupture is peculiar in that the effective velocity is slow, 3.5 Km.s À1 for a high stress-drop, 21 –30 MPa. We propose that rupture was due to the reactivation by hydraulic embrittlement of a inherited major lithospheric fault within the subducting plate. The stress-drop suggests that the region of the slab between planes of the double seismic zone can sustain high stresses. Citation: Peyrat, S., et al. (2006), Tarapacá intermediate-depth earthquake (Mw 7.7, 2005, northern Chile): A slab-pull event with horizontal fault plane constrained from seismologic and geodetic observations, Geophys