On high-frequency spheroidal modes and the structure of the upper mantle

Summary. The asymptotic properties of spheroidal mode dispersion at high frequency for fixed phase velocity are related to the intercept times τ_β( p ) for P and S waves. If the mode eigenfrequency and the ratio of horizontal to vertical displacement at the surface for the mode are known τ_α( p ) and τ_β( p ) may be separately estimated. If discontinuities exist in the velocity model then 'solotone' effects occur, in frequency at fixed slowness, and in τ_α( p ), τ_β( p ) estimated from the mode dispersion as a function of slowness. The coupling of P and S waves in the spheroidal modes means that the interaction of P waves with upper-mantle discontinuities affects also the estimates of the S wave τ_β( p ) values for which the corresponding turning points lie in the lower mantle. The asymptotic formalism also shows that sharp pulses formed by superposition of spheroidal modes correspond to multiple PS reflections.
A study of τ_α( p ), τ_β( p ) estimates derived from spheroidal modes with periods from 45–50s, calculated for model 1066B, shows that even in the presence of strong upper-mantle discontinuities the errors in intercept time are only about one-tenth of a period. The asymptotic properties may there-for provide a useful means of estimating intercept times from modes with a few seconds period as a supplement to travel-time methods.     

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.     

A new derivation of the ratio Q(wobble)/Qμ (mantle)

Summary An extension of the Love-Larmor theory to a low-loss unelastic earth model, leads to the surprisingly simple approximation
   
where τ_s= 447.4 sidereal day is the static wobble period, τ_R= 306 sidereal day is the rigid-earth wobble period and τ_w= 433 sidereal day is the observed Chandler period. Q _W, Q _μ are the respective average Q values of the wobble and the Earth's mantle at τ_W. The known numerical factor F is only slightly dependent on the Earth structure.     

Seismic sources with observable glut moments of spatial degree two

Let ζ_Λ and r _Λ. be the hypocentral position and time of an extended indigenous seismic source. Backus showed that the force moment tensors of the source, Γ^{( m +1, n )}(ζ_Λ, r _Λ), determine and are determined by the motion which the source produces. For small m + n , only the long-period motion is relevant. The glut moment tensor Λ^{( m,n )} (ζ_Λ, r _Λ.) can be calculated uniquely from γ^{( m +1, n )}(ζ_Λ r _Λ) only if m = 0 or m = 1. The tensor G =Λ^(2,0)(ζ_Λ) gives the spatial variance tensor W_Λ of the source, and W_Λ. roughly describes the size, shape and orientation of the source region. Therefore the failure of the observed F =Γ^(3,0)(ζ_Λ) to determine G uniquely is of seismological interest. In the present paper we show that F determines G uniquely if we assume the source to be a simple straight line source (SSLS) or an ideal fault in an isotropic medium with isotropic prestress (IFIMIP). We give tests on F which determine whether it can come from a SSLS, from an IFIMIP or from a simple plane surface source (SPSS). If we assume the source to be a SPSS then knowing F and the fault plane determines G to within an unknown scalar multiple of a certain tensor tangent to the fault plane. Moreover F determines the fault plane uniquely unless F can come from a SSLS. If it can, then F determines this virtual source line uniquely, and F permits the fault plane to be any plane containing the virtual source line.     

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.     

Stochastic behaviour and scaling laws in geoelectrical signals measured in a seismic area of southern Italy

We explore the inner dynamics of daily geoelectrical time series measured in a seismic area of the southern Apennine chain (southern Italy). Autoregressive models and the Higuchi fractal method are applied to extract maximum quantitative information about the time dynamics from these geoelectrical signals. First, the predictability of the geoelectrical measurements is investigated using autoregressive models. The procedure is based on two forecasting approaches: the global and the local autoregressive approximations. The first views the data as a realization of a linear stochastic process, whereas the second considers the data points as a realization of a deterministic process, which may be non-linear. Comparison of the predictive skills of the two techniques allows discrimination between low-dimensional chaos and stochastic dynamics. Our findings suggest that the physical systems governing electrical phenomena are characterized by a very large number of degrees of freedom and can be described only with statistical laws. Second, we investigate the stochastic properties of the same geoelectrical signals, searching for scaling laws in the power spectrum. The spectrum fits a power law P (  f )∝  f  ^−α , with the scaling exponent α a typical fingerprint of fractional Brownian processes. In this analysis we apply the Higuchi method, which gives a linear relationship between the fractal dimension D _Σ and the spectral power law scaling index α : D _Σ=(3− α )/2. This analysis highlights the stochastic nature of geoelectrical signals recorded in this seismic area of southern Italy.     

Observation of Coriolis coupled modes below 1 mHz

We present observations of spectral energy at toroidal mode frequencies in vertical seismic recordings of the 1998 Balleny Islands earthquake. Since toroidal modes on a spherically symmetric, nonrotating Earth have horizontally polarized particle motion these observations call for an explanation. We first rule out local and instrumental effects as being responsible for the verticalcomponent signal of the toroidal modes ₀ T ₃ (0.59 mHz) and ₀ T ₄ (0.77 mHz). The global effects that we consider are general heterogeneous mantle structure, ellipticity of figure and rotation. We find that rotation through Coriolis coupling of loworder spheroidal and toroidal oscillations is the dominant mechanism.     

Shock-induced temperatures of MgO

Shock-compressed MgO radiates thermally at temperatures between 2900 and 3700 K in the 170–200 GPa pressure range. A simple energy-transport model of the shocked-MgO-targets distinguishes between different shock-induced radiation sources in these targets and provides estimates of spectral absorption-coefficients, α _ΛMgO, for shocked MgO (e.g. at 203 GPa, α _ΛMgO˜ 630, 7500, 4200 and 3800 m⁻¹, at 450, 600, 750 and 900 nm, respectively). The experimentally inferred temperatures of the shock-compressed states of MgO are consistent with temperatures calculated for MgO assuming that (1) it deforms as an elastic fluid, (2) has a Dulong-Petit value for specific heat at constant volume in its shocked state, and (3) undergoes no phase transformation below 200 GPa.     

Asymptotic results for elastodynamic propagator matrices in plane-stratified and spherically-stratified earth models

Summary. Asymptotic results, applicable at high frequencies, are established for elastodynamic propagator or transfer matrices. We are concerned here with a single layer, by which we shall mean a region of a spherically-symmetric earth model in which physical properties vary smoothly with radius, or a region of a plane-stratified earth model in which physical properties vary smoothly with depth. Results are obtained for the spheroidal and toroidal oscillations of a solid spherical layer, the spheroidal oscillations of a spherical fluid layer, P-SV and SH waves in a plane solid layer and P waves in a plane fluid layer. Each of these results is given to first order in ω⁻¹ where ω is the angular frequency.     

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.     

Frequency dependence of the quality factor in the upper crust: a deep seismic sounding approach

Summary. Ray amplitudes are computed in a one-dimensional velocity structure where the quality factor Q varies continuously with depth. An iterative process is then proposed to derive the quality factor distribution in the upper crust from deep sounding data. Results for compressional waves in the French Massif Central and for a signal frequency close to 20 Hz show that Q increases in a rather linear way from about 40 in surface up to 600 at 7 km depth. Q seems to be higher in the Central Alps: 180 in surface, 1600 at 5 km depth (at 20 Hz). Using proper signal processing, the frequency dependence of Q is finally investigated in the frequency range 10–25 Hz. The results indicate a dependence of the form Q = Q ₀ f ^α, where α= 0.25 ± 0.1.     

Palaeomagnetism of the Fongen-Hyllingen gabbro complex, southern Scandinavian Caledonides; plate rotation or polar shift?

Summary. Palaeomagnetic results are presented from the c . 160 km² Caledonian synorogenic layered Fongen-Hyllingen gabbro complex (of probable late Silurian age) located about 75 km SE of Trondheim, Norway, in the allochthonous Seve-Kdli Nappe Complex. A total of 80 oriented samples from eight sites in the northern part of the gabbro were investigated. After detailed af demagnetization two stable high coercivity components emerge: one with a well defined NW direction with D =325°, I =−21° (α₉₅=8°, N =8), and another, less well defined, probably younger, SW direction with D = 237°, I = 6° (α₉₅= 9°, N = 8). Correction for dip of these two directions gives D = 329°, I =−7° (α₉₅= 10°) and D = 238°, I =−11° (α₉₅= 12°), respectively. The corresponding pole positions are P ₁ : 19° N, 225° E and P ₂: 19° S, 308° E, respectively. The reversed pole -P ₂ of the SW direction lies close to other NW European palaeomagnetic poles of Caledonian, Upper Silurian-Lower Devonian age. However, the dominant pole PI is far away from these, and could be due to a late Caledonian geomagnetic excursion of considerable duration; or it could record a c . 90° rotation around a vertical axis of a crustal block within the Scandinavian Caledonides. Block rotation could have been related to nappe translation, although geological observations do not at present appear to support the occurrence of such an event.     

The attenuation mechanism of seismic waves in northwestern Himalayas

We analysed local earthquake waveforms recorded on a broad-band seismic network in northwestern Himalayas to compute the intrinsic and scattered attenuation parameters from coda waves. Similar to other tectonically active and heterogeneous regions, attenuation-frequency relation for western Himalaya is   Q ⁻¹ _c = (113 ± 7)  f ^(1.01±0.05)  where   Q_c   is the coda Q parameter. Intrinsic  ( Q ⁻¹ _i )  and scattering  ( Q ⁻¹ _s )  attenuations was separated using   Q_c   and direct S -wave Q data  ( Q_d )  . It is observed that estimated   Q ⁻¹ _c   is close to   Q ⁻¹ _i   and both of them are much larger than   Q ⁻¹ _s   suggesting that coda decay is predominantly caused by intrinsic attenuation. At higher frequencies, both the attenuation parameters   Q_c   and,   Q_d   are similar indicating that coda is predominantly composed of back-scattered S waves at these frequencies.     

Regional variations of heat flow differences with depth in Alberta, Canada

Summary. Differences between estimated average heat flow values for the Mesozoic and Cenozoic formations ( Q ₁) and estimated average heat flow values for the Palaeozoic formations below the erosional unconformity ( Q ₂) are calculated for the Alberta part of the western Canadian sedimentary basin. Significant heat flow differences exist for these two intervals and the map of Δ Q = Q ₁– Q ₂ shows that Q ₂ is generally greater than Q ₁ in the western and south-western part of Alberta, while in the northern part of the province Q ₂ is generally less than Q ₁. The regional variations of Δ Q are large, with standard deviation of 26 mW m⁻² and average value –13.5 mW m⁻². A regional trend of Δ Q correlates with topographic relief and the hydraulic head variations in the basin. It is shown that there is a heat flow increase with depth in water recharge areas and a decrease in heat flow with depth in the low topographic elevation water discharge areas when comparing the average heat flow in Mesozoic + Cenozoic and Palaeozoic formations.     

Stacking for the frequencies and Qs of 0S0 and 1S0

Summary. Using nine IDA records for the Indonesian earthquake of 1977 August 19, we have formed an optimal linear combination of the records and have measured the frequency and Q of ₀ S ₀ and ₁ S ₀. The frequency was measured using the moment ratio method. The attenuation was measured by the minimum width method and by the time-lapse method. The frequency and attenuation were measured simultaneously by varying them to obtain a best fit to the data. A 2000-hr stack, the sum of nine individual records, for ₀ S ₀ gave a frequency of 0.814664 mHz±4 ppm. The values for the Q of ₀ S ₀ for the three different methods of measurement were 5600,5833 and 5700, respectively. The error in the estimates of Q ^-1 is about 5 per cent for the minimum power method. For ₁ S ₀ a 300-hr stack yielded a frequency of 1.63151 mHz±30 ppm. The values of Q for this mode were 1960, 1800 and 1850, respectively, with an error in Q ^-1 of about 12 per cent for the minimum power method.     

Three-dimensional VP and VP /VS models of the upper crust in the Friuli area (northeastern Italy)

3-D images of P velocity and P - to S -velocity ratio have been produced for the upper crust of the Friuli area (northeastern Italy) using local earthquake tomography. The data consist of 2565 P and 930 S arrival times of high quality. The best-fitting V _P and V _P / V _S 1-D models were computed before the 3-D inversion. V _P was measured on two rock samples representative of the investigated upper layers of the Friuli crust. The tomographic V _P model was used for modelling the gravity anomalies, by converting the velocity values into densities along three vertical cross-sections. The computed gravity anomalies were optimized with respect to the observed gravity anomalies. The crust investigated is characterized by sharp lateral and deep V _P and V _P / V _S anomalies that are associated with the complex geological structure. High V _P / V _S values are associated with highly fractured zones related to the main faulting pattern. The relocated seismicity is generally associated with sharp variations in the V _P / V _S anomalies. The V _P images show a high-velocity body below 6 km depth in the central part of the Friuli area, marked also by strong V _P / V _S heterogeneities, and this is interpreted as a tectonic wedge. Comparison with the distribution of earthquakes supports the hypothesis that the tectonic wedge controls most of the seismicity and can be considered to be the main seismogenic zone in the Friuli area.     

The Influence of Surface Inhomogeneities On Deep Electromagnetic Soundings of the Earth

Summary. Asymptotic expressions for components of the electromagnetic field of a grounded electric dipole are considered for the model consisting of a thin surface-layer overlapping a stratified medium with a highly resistive screen on the roof. It is shown that the method of spatial derivatives makes it possible to obtain proper estimates of the impedance at distances of r ≥|λ₀| from the nearest edge of the surface anomaly (|λ₀| being the effective depth of the field penetration in the underlying section). the magnetotelluric methods allow one to obtain the true values of impedance, provided r ≥ max {|λ₀|, |/( S ⁻¹+ Z ⁻₀|^1/2} where S is the integrated conductivity of the surface layer, is the transverse resistance of the screen, and Z ⁻₀ is the Tikhonov—Cagniard impedance for the medium underlying the surface layer.     

Multichannel appraisal deconvolution

Summary The single channel scalar deconvolution method presented by Oldenburg has been extended to include N channels of data and vector models of the form     ( t ) = ( m ₁( t ), m ₂( t ), …, m _α( t ))^T. The solution has its foundation in the linear inverse theory of Backus & Gilbert and is effected by computing a set of N filters, which, when convolved with the data, yield unique averages of one of the scalar functions of the model. Those averages are the summation of the scalar model convolved with a primary averaging function plus contamination from secondary averaging functions convolved with other model components. It is shown how a set of suitably selected weights can annihilate these secondary averaging functions and thereby greatly simplify the interpretation. The computations are efficiently carried out in the frequency domain and require the inversion of an N × N Hermitian matrix at each frequency. As a type example, we have shown how the time varying elements of a seismic moment tensor might be computed from a set of seismograms.     

Palaeomagnetism of a sequence of red beds of the Middle and Upper Sections of Pagnazo Group (Argentina) and the correlation of Upper Palaeozoic-Lower Mesozoic rocks

Palaeomagnetic data from 182 hand samples collected in a rock sequence of about 620-m of red beds of Late Palaeozoic to Early Triassic age exposed in north-western Argentina (30.3° S 67.7° W), are given.
After cleaning, the majority of the Upper Palaeozoic samples (Middle Section of Paganzo Group) show reversed polarity and yield a palaeomagnetic pole at 78° S 249° E (α₉₅= 3°). They also record a polarity transition which we have correlated with the Middle Permian Quebrada del Pimiento Normal Event. The position of the palaeomagnetic pole and the K-Ar age of a basalatic sill at the base of the sequence support this correlation.
Stable remanent magnetization has been isolated in the majority of samples from the Upper Section of the Paganzo Group; it is predominantly reversed and reveals three normal events and also three geomagnetic excursions suggesting an Illawarra Zone age (post Kiaman, Late Tatarian-Early Scythian). The palaeomagnetic pole of the reversely magnetized samples is located at 75° S 285° E(α₉₅= 13°).
The red beds involved in this study are correlated with red beds from the Corumbataí Formation (State of Paraná, Brazil) and with igneous rocks from the Quebrada del Pimiento Formation (Province of Mendoza, Argentina).
The South American Middle and Upper Permian, Upper Permian—Lower Triassic, Lower, Middle and Upper Triassic and Middle Jurassic palaeomagnetic poles reflect a quasistatic period with mean pole at 82° S 244° E, (α₉₅= 4°) which followed the South American Late Palaeozoic polar shift.     

Experimental study of the anisotropy of longitudinal and transverse waves from local earthquake records

Summary. The paper gives the results of a study of the anisotropy of seismic wave velocities within the Ashkhabad test field in Central Asia. The anisotropy was studied by analysing variations in the values of apparent velocities of first arrivals for epicentral distances ranging from 30 to 130 km and by analysing the delays (Δ t_s1-_s2 ) between the arrival times of shear waves with different polarizations.
The velocities of P -waves vary with azimuth from 5.3 to 6.27 km s^-1 and the velocities of S -waves vary from 3.15 to 3.5 km s^-1.
The delay times Δ t_S1 - _S2 depend on the direction of the propagation. The character of the variation of the propagation velocity of the longitudinal wave, the presence of two differently polarized shear waves S 1 and S 2 propagating at different velocities, and the character of the distribution of Δ t_S1 - _S2 on the stereogram suggest that the symmetry of the anisotropic medium is close to hexagonal with a nearly horizontal symmetry axis coinciding with the direction of maximal velocity. The azimuth of the symmetry axis of the medium is 140° and coincides with the direction of geological faults.