Average attenuation of 0.7–5.0 Hz Lg waves and magnitude scale determination for the region bounding the western branch of the East African Rift

The investigation of L g attenuation characteristics in the region bounding the western branch of the East African rift system using digital recordings from a seismic network located along the rift between Lake Rukwa and Lake Malawi is reported. A set of 24 recordings of L g waves from 12 regional earthquakes has been used for the determination of anelastic attenuation, Q _Lg , and regional body-wave magnitude, m _{b Lg} , scale. The events used have body-wave magnitudes, m _b , between 4.6 and 5.5, which have been determined teleseismically and listed in ISC bulletins. The data were time-domain displacement amplitudes measured at 10 different frequencies (0.7–5.0  Hz). Q _Lg and its frequency dependence, η , in the region can be represented in the form Q _Lg = (186.2 ± 6.5)  f  ^(0.78±0.05). This model is in agreement with models established in other active tectonic regions. The L g -wave-based magnitude formula for the region is given by m _{b Lg} = log   A + (3.76 ± 0.38)  log   D − (5.72 ± 1.06), where A is a half-peak-to-peak maximum amplitude of the 1  s L g wave amplitude in microns and D is the epicentral distance in kilometres. Magnitude results for the 12 regional earthquakes tested are in good agreement with the ISC body-wave magnitude scale.     

Recognizing explosion sites without seismogram readings: neural network analysis of envelope-transformed multistation SP recordings 3–6 Hz

Seismic-waveform similarities for closely spaced earthquakes and explosions in particular are well established observationally. In many industrialized countries of low seismicity more than 90 per cent of seismic event recordings stem from chemical explosions and thus contribute significantly to the daily analyst workload. In this study we explore the possibility of using envelope waveforms from a priori known explosion sites (learning) for recognizing subsequent explosions from the same site excluding any analyst interference. To ensure high signal correlation while retaining good SNRs we used envelope-transformed waveforms, including both the P and Lg arrivals. To ensure good spatial resolution we used multistation (network) recordings. The interpolation and approximation neural network (IANN) of Winston (1993) was used for teaching the computer to recognize new explosion recordings from a specific site using detector output event files of waveforms only. The IANN output is a single number between 0 and 1, and on this scale an acceptance threshold of 0.4 proved appropriate. We obtained 100 per cent correct decisions between two sets of 'site explosions' and hundreds of 'non-site' explosions/earthquakes using data files from the Norwegian Seismograph Network.     

Smaller source earthquakes and improved measuring techniques allow the largest earthquakes in Iceland to be stress forecast (with hindsight)

A group of three earthquakes in 2000 June in SW Iceland included the two largest earthquakes in Iceland in the past 30 yr. Previously, temporal changes in shear-wave splitting had not been recognized before these earthquakes as there were too few small earthquakes to provide adequate shear-wave data, and they were not stress forecast, even with hindsight. These large earthquakes were subject to a special investigation by the European Community funded PREPARED Project during which the seismic catalogue was extended to include smaller magnitude earthquakes. This more detailed data set, together with a semi-automatic programme for measuring the parameters of shear-wave splitting greatly increased the number of time-delay measurements.
The new measurements displayed the typical temporal variations before larger earthquakes as seen elsewhere: a long-term increase in time delays, interpreted as stress accumulation before the earthquake; and a decrease, interpreted as crack coalescence, immediately prior to the earthquake. The logarithms of the durations of both the implied accumulation of stress and the crack coalescence have the same self-similar relationships to earthquake magnitude as found elsewhere in Iceland. This means that, in principle, the time and magnitude of the larger earthquakes could have been stress forecast in real time had the smaller source earthquakes of the extended catalogue and the improved measuring procedures been available at the time.     

Waveform inversion of mantle Love waves:the Born seismogram approach

Summary. Normal mode theory, extended to the slightly laterally heterogeneous earth by the first-order Born approximation, is applied to the waveform inversion of mantle Love wave (200–500 s) for the Earth's lateral heterogeneity at l = 2 and a spherically symmétric anelasticity ( Q _μ) structure. The data are from the Global Digital Seismograph Network (GDSN). The l =2 pattern is very similar to the results of other studies that used either different méthods, such as phase velocity measurements and multiplet location measurements, or a different data set, such as mantle Rayleigh waves from different instruments. The results are carefully analysed for variance reduction and are most naturally explained by heterogeneity in the upper 420 km. Because of the poor resolution of the data set for the deep interior, however, a fairly large heterogeneity in the transition zones, of the order of up to 3.5 per cent in shear wave velocity, is allowed. It is noteworthy that Love waves of this period range cannot constrain the structure below 420 km and thus any model presented by similar studies below this depth are likely to be constrained by Rayleigh waves (spheroidal modes) only.
The calculated modal Q values for the obtained Q _μ model fall within the error bars of the observations. The result demonstrates the discrepancy of Rayleigh wave Q and Love wave Q and indicates that care must be taken when both Rayleigh and Love wave data, including amplitude information, are inverted simultaneously.
Anomalous amplitude inversions of G2 and G3, for example, are observed for some source-receiver pairs. This is due to multipathing effects. One example near the epicentral region, which is modelled by the obtained l = 2 heterogeneity, is shown.     

Comment on 'Systematic determination of earthquake rupture directivity and fault planes from analysis of long-period P-wave spectra' by L. M. Warren and P. M. Shearer

Warren and Shearer describe a method of estimating the duration of P pulses radiated by earthquakes, on the assumption that the source is a unilateral fracture. The estimates are made in the frequency domain. The estimates obtained by Warren and Shearer for seven of the earthquakes are compared to durations estimated here in the time domain; the time-domain measurements being made on broad-band seismograms (∼0.1–4.0 Hz) derived by filtering from short-period recordings. Overall, the time-domain method indicates that the pulse duration of the earthquakes studied here range from 2.0 to 7.6 s, whereas the estimates from the results of Warren and Shearer range from 7.1 to 9.8 s. This suggests that the method of Warren and Shearer cannot resolve pulse lengths less than about 7 s. The suggestion is supported by the estimates of the duration of rupture—fault length/speed of rupture—implied by the results of Warren and Shearer. For, although the estimated fault-length ranges from 0.8 km to over 40 km most rupture durations are around 8 s.     

Evidence of spatial variations of attenuation in the western Pyrenean range

Summary. Spectral attenuation of coda waves has been studied in the range 2–40 Hz from local events recorded in the western Pyrenean range from 1980 to 1982. Q _c was obtained using a single scattering model of S -waves for different segments of the coda. An increase of Q _c with lapse time was found and attributed to a rapid increase of Q _β with depth.
Three groups of events were selected from distinct focal areas. Two data sets are mainly composed of aftershocks of moderate earthquakes of magnitude 5.1 and 4.8, respectively. No moderate earthquake occurred in the third area in the few years preceding or following the selected events. Use of stations close to epicentres allowed sampling of the coda at very short lapse times and then study of small, distinct scattering volumes. Noticeable differences were found between the three studied areas and attributed to spatial rather than temporal variations.
The Q _c frequency dependence was studied according to Q _c= qf ^α. α is found to range from 0.7 to 1.1 and q from 30 to 140. These values are in agreement with those found in other tectonic areas. It is shown that scattering is the dominant attenuation process below 10Hz.     

Seismotectonics and analysis of earthquakes from the Hindukush region

Summary. Over 80 earthquakes, exclusively from the Hindukush focal region, which were recorded at the Gauribidanur seismic array (GBA) have been used in this study. These events have similar epicentral distances and a narrow azimuthal range from GBA but varying focal depths from 10 to 240 km. A fault plane dipping steeply (75°) in the north-west direction and striking N 66° E has been investigated on the basis of the spatial distribution of earthquakes in two vertical planes through 68° E and 32° N. Short period P -wave recordings up to 30 s were processed using the adaptive cross-correlation filtering technique. Slowness and azimuthal anomalies were obtained for first arrivals. These anomalies show positive as well as negative bias and are attributed to a steep velocity gradient in the upper mantle between the 400–700 km depth range where the seismic rays have their maximum penetration. Relative time residuals between the stations of GBA owe their origin very near to the surface beneath the array. A search of the signals across the array revealed that most of the events occurring at shallower depths had complex signatures as compared to the deeper events. The structure near the source region, complicated source functions and the scattering confined to the crust—upper mantle near source are mainly responsible for the complexity of the Hindukush earthquakes as the transmission zone of the ray tubes from turning point to the recording station is practically the same.     

The effects of seismic rotations on inertial sensors

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.     

Anelasticity of the Crust and Upper Mantle of South America From the Inversion of Observed Surface Wave Attenuation

Fundamental-mode Rayleigh and Love waves generated by several earthquakes situated along great-circle paths between pairs of seismograph stations have been analysed to obtain coefficients of attenuation, group velocities, phase velocities, and specific quality factors in the period range 18–80s in two regions of the South American continent. One set of paths crosses the shield region which lies on the eastern coast and another set traverses the mountainous region inland. the average attenuation coefficient values are clearly higher in the tectonically active western region throughout the entire period range than in the eastern or shield region.
Inversion of the attenuation data yielded shear wave internal friction ( Q ^-1_β) models as a function of depth in the crust and upper mantle in both regions. A low- Q zone below the lithosphere is prominent in both regions. the results show that substantial variations of Q _β occur in the two regions of South America. the Q_β values were found to be inversely related to the heat flow values or to the temperature.     

Seismic reflection studies using local earthquake sources

Summary. A method is presented for processing three-component digital recordings of micro-earthquakes to obtain near-vertical reflection profiles in regions of shallow seismicity. The processing includes magnitude and focal-depth normalization and event stacking, where stacking is by small localized groups, with ray theoretical time and distance corrections applied to compensate for varying focal depths. In areas with high seismicity, this procedure allows earthquakes to be treated as "controlled" sources to probe layered structures of the deep crust and upper mantle. The validity of our approach is demonstrated using S-waves from aftershocks of the Borah Peak, Idaho, earthquake (Ms = 7.3) of 1983.     

Which is the correct stress strain relation for the anelasticity of the Earth's interior?

Summary. The effect of anelasticity on the periods of the Earth's free oscillation (toroidal modes) is examined in the light of several different studies leading to different results. The possibility of discrimination between the different analytic models for the anelasticity of the Earth's interior is discussed. A splitting of the spectral lines caused by Q is also suggested.     

Frequency–size distribution of global seismicity seen frombroad-band radiated energy

The frequency–energy distribution of global seismicity is studied using broad-band radiated energy of shallow earthquakes from January 1987 to December 1994 estimated by NEIC. Rank-ordering statistics are applied to enhance the resolution in retrieving the power-law distribution with undersampled data, namely a few tens of events. Seen in the perspective of broad-band radiated energy with higher resolution, a single (Gutenberg–Richter-type) power-law distribution can fit the data. For earthquakes with energy larger than 10¹⁴ J, the number N of events with energy E depends on E via N∝E ^−B , with the scaling constant B = 0.64 ± 0.04, corresponding to b = 0.95 ± 0.06. This relation is different from that of scalar seismic moment, which shows a transition of power-law distributions between small and large earthquakes. To demonstrate such a difference we use the same set of earthquakes with both broad-band energy estimation and CMT estimation. It is found that for the same data set, the energy distribution and the moment distribution show different patterns. The moment distribution has a clear kink between small and large earthquakes, while the energy distribution shows a single power law with no convincing kink between small and large earthquakes. To investigate the effect of different focal mechanisms and different seismic regions, events with strike-slip mechanisms and events within the Japan–Kuril region are considered. For these subsets of events, a similar pattern exists, in which the moment distribution shows a kink between small and large earthquakes, while the energy distribution shows a single power law.     

Regional mid-ocean stress and a proposed focal mechanism of 'stress-discordant' earthquakes

Summary. Rock stress measurements in Iceland show maximum horizontal compression perpendicular to the trend of Reykjanes Ridge crest and of its extension, the active volcanic zone of Iceland. Fault-plane solutions of dormant stage earthquakes are consistent with the measured stress orientations, but strike—slip earthquakes associated with volcanic surges and some earthquake swarms in active geothermal areas exhibit apparent reversals of mechanism and are here defined as 'stress-discordant' in the sense that they yield deduced stress orientations 90° from the regional stress field as determined by hydrofracturing and strain relief methods. It is proposed, supported by comparison with the pore-pressure induced Denver earthquakes, that the 'stress-discordant' volcanic earthquakes are triggered by increased pore pressure and probably involve stick-slip motion similar to that reported for some laboratory tests of the pore pressure effect, characterized by gradual onset and sudden stopping of each slip episode. The question is raised as to whether stress-discordant earthquakes are dominated by a stopping phase or terminal shock with consequent reversal of the deduced shear couple. A possible stopping mechanism is suggested: the dilatant stiffening of fault gouge during shear.
It is proposed that direct measurements of stress orientation be made by hydrofracturing tests at other places along the mid-ocean ridge crest and on the margins of the Red Sea and East African rifts. The Icelandic stress data indicate the need for sceptical re-examination of some fundamentals of plate tectonics theory.     

Vertical profiles of aeolian sand mass flux

Vertical profiles of the horizontal mass flux of blown sand are investigated experimentally using a passive vertical array in a wind tunnel. Considering lower sampling efficiency of the sand trap in the near-bed region, this investigation is complemented by the measurements of the longitudinal profiles of mass flux made using a horizontal sand trap. The experiments were conducted with two test sands and five different stream velocities.In the upper part of the vertical profile, the measured data exhibit an exponential decay distribution with a positive deviation occurring in the near-bed region. The measured longitudinal profiles are similar to the measured vertical profiles. Linking both profiles and the modes of sand transport, it is possible that saltating sand grains give rise to the well-known exponential decay distribution of mass flux, and that creeping and reptating grains force a deviation from it. A simple equation applicable for both the vertical and the longitudinal sand mass flux variations is introduced and the parameters are estimated from experimental data.     

Improved resolution of complex eigenfrequencies in analytically continued seismic spectra

Summary. The response of the Earth to an earthquake is a transient that is effectively zero several days after the event. A recording of the event, of finite duration in time, has a Fourier spectrum that is an entire, or integral, analytic function of frequency. We present a very simple procedure for computing the Fourier spectrum as a function of complex frequency; the analytically continued spectrum. By investigating the properties of the analytically continued spectrum we show how to extract high- Q modes, how to estimate Q either from the amplitude or from the width of a resonance function, and how to improve the resolution of splitting to the theoretical maximum. Examples of these procedures, using observed data, are presented.     

Coupled normal-mode sensitivity to inner-core shear velocity and attenuation

We have studied the response of normal modes to perturbations in inner-core shear velocity and attenuation, using fully coupled mode synthetics. Our results indicate that (i) mode pairs   _n S _l – _{n ±1} S _l   are strongly coupled by anelasticity, (ii) this coupling causes shear velocity perturbations to strongly affect the Q values of modes through exchange of inner-core characteristics, (iii) there is no evidence for a weakly attenuating inner core in shear, and (iv) the discrepancy between attenuation models returned from normal modes and body waves is small. These results suggest that inversions for inner-core attenuation and shear velocity should be performed simultaneously and should take account of the strong cross-coupling due to attenuation.     

Focal mechanisms of small earthquakes in the southeastern Brazilian shield: a test of stress models of the South American plate

Focal mechanisms of small earthquakes with magnitudes of about 3 in the SE Brazilian shield are calculated using S / P amplitude ratios. Low attenuation ( Q _p from 400 to 800) in the shield upper-crustal layers allowed sharp S arrivals to be recorded up to distances of 100 km. Besides P -wave polarities, SH -wave first motions were also used to constrain the nodal-plane orientations. Normal and reverse faulting mechanisms with strike-slip components were found. The inversion of four mechanisms to estimate the stress tensor indicated a strike-slip stress regime with roughly E–W-orientated σ ₁ and N–S σ ₃. Both the orientations and the shape factor ( φ =0.7) of the inverted stress are in excellent agreement with theoretical predictions for that part of Brazil from the driving-force model of Coblentz & Richardson (1996) . Good agreement with the nature of the stress, as well as its orientation, was also found for the model of Meijer (1995) . Both of these theoretical models include spreading stresses along the continent/ocean lithospheric transition. Because the earthquakes are more than 300 km from the continental shelf they should not be affected by the local flexural forces caused by sediment load in the marginal basins. The agreement between observed and theoretical stresses then confirms the importance of continental spreading forces in modelling intraplate stresses.     

Evidence for different scaling of earthquake source parameters for large earthquakes depending on faulting mechanism

Scaling relationships between seismic moment, rupture length, and rupture width have been examined. For this purpose, the data from several previous studies have been merged into a database containing more than 550 events. For large earthquakes, a dependence of scaling on faulting mechanism has been found. Whereas small and large dip-slip earthquakes scale in the same way, the self-similarity of earthquakes breaks down for large strike-slip events. Furthermore, no significant differences in scaling could be found between normal and reverse earthquakes and between earthquakes from different regions. Since the thickness of the seismogenic layer limits fault widths, most strike-slip earthquakes are limited to rupture widths of between 15 and 30 km while the rupture length is not limited. The aspect ratio of dip-slip earthquakes is similar for all earthquake sizes. Hence, the limitation in rupture width seems to control the maximum possible rupture length for these events. The different behaviour of strike-slip and dip-slip earthquakes can be explained by rupture dynamics and geological fault growth. If faults are segmented, with the thickness of the seismogenic layer controlling the length of each segment, strike-slip earthquakes might rupture connected segments more easily than dip-slip events, and thus could produce longer ruptures than dip-slip events of the same width     

Local and regional components of western Aegean deformation extracted from 100 years of geodetic displacement measurements

Our objectives are as follows. First, we wish to develop a methodology to recover the long-term component of deformation from any set of distributed, time-averaged geodetic strain measurements that were subject to seismic disturbance, given a catalogue of local seismicity that occurred during the measurement period. Second, using seismic and geodetic data sets that span approximately 100 years, we apply this technique in the western Aegean to assess the role of local seismicity in regional deformation. The methodology is developed using a model for crustal deformation constructed from a long-term, smooth regional strain field combined with instantaneous, local perturbations from upper-crustal earthquakes approximated by static elastic dislocations. By inverting geodetic displacements for the smooth field while simultaneously floating influential but uncertain earthquake source parameters, an estimate of the regional component of deformation that is approximately independent of the seismicity can be made. In the western Aegean we find that the horizontal component of regional deformation can be described with minor inaccuracy by a quadratic relative displacement field. The principal horizontal extensional axes calculated from the regionally smooth displacement field agree in orientation with the T-axes of earthquakes in the region. These observations indicate that the instantaneous elastic strain of the 10 km thick seismogenic layer is driven by a stress field that is smooth on the scale of the geodetic network as a whole, 200-300 km.