A three-dimensional radially anisotropic model of shear velocity in the whole mantle

We present a 3-D radially anisotropic S velocity model of the whole mantle (SAW642AN), obtained using a large three component surface and body waveform data set and an iterative inversion for structure and source parameters based on Non-linear Asymptotic Coupling Theory (NACT). The model is parametrized in level 4 spherical splines, which have a spacing of ∼ 8°. The model shows a link between mantle flow and anisotropy in a variety of depth ranges. In the uppermost mantle, we confirm observations of regions with   V_SH > V_SV   starting at ∼80 km under oceanic regions and ∼200 km under stable continental lithosphere, suggesting horizontal flow beneath the lithosphere. We also observe a   V_SV > V_SH   signature at ∼150–300 km depth beneath major ridge systems with amplitude correlated with spreading rate for fast-spreading segments. In the transition zone (400–700 km depth), regions of subducted slab material are associated with   V_SV > V_SH   , while the ridge signal decreases. While the mid-mantle has lower amplitude anisotropy (<1 per cent), we also confirm the observation of radially symmetric   V_SH > V_SV   in the lowermost 300 km, which appears to be a robust conclusion, despite an error in our previous paper which has been corrected here. The 3-D deviations from this signature are associated with the large-scale low-velocity superplumes under the central Pacific and Africa, suggesting that   V_SH > V_SV   is generated in the predominant horizontal flow of a mechanical boundary layer, with a change in signature related to transition to upwelling at the superplumes.     

Seismic anisotropy within the uppermost mantle of southern Germany

This paper presents an updated interpretation of seismic anisotropy within the uppermost mantle of southern Germany. The dense network of reversed and crossing refraction profiles in this area made it possible to observe almost 900 traveltimes of the P_n phase that could be effectively used in a time-term analysis to determine horizontal velocity distribution immediately below the Moho. For 12 crossing profiles, amplitude ratios of the P_n phase compared to the dominant crustal phase were utilized to resolve azimuthally dependent velocity gradients with depth. A P -wave anisotropy of 3–4 per cent in a horizontal plane immediately below the Moho at a depth of 30 km, increasing to 11 per cent at a depth of 40 km, was determined. For the axis of the highest velocity of about 8.03 km s⁻¹ at a depth of 30 km a direction of N31°F was obtained. The azimuthal dependence of the observed P_n amplitude is explained by an azimuth-dependent sub-Moho velocity gradient decreasing from 0.06 s⁻¹ in the fast direction to 0 s⁻¹ in the slow direction of horizontal P -wave velocity. From the seismic results in this study a petrological model suggesting a change of modal composition and percentage of oriented olivine with depth was derived.     

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.     

Localization of gravity and topography: constraints on the tectonics and mantle dynamics of Venus

We develop a method for spatio-spectral localization of harmonic data on a sphere and use it to interpret recent high-resolution global estimates of the gravity and topography of Venus in the context of geodynamical models. Our approach applies equally to the simple spatial windowing of harmonic data and to variable-length-scale analyses, which are analogous to a wavelet transform in the Cartesian domain. Using the variable-length-scale approach, we calculate the localized RMS amplitudes of gravity and topography, as well as the spectral admittance between the two fields, as functions of position and wavelength. The observed admittances over 10 per cent of the surface of Venus (highland plateaus and tessera regions) are consistent with isostatic compensation of topography by variations in crustal thickness, while admittances over the remaining 90 per cent of the surface (rises, plains and lowlands) indicate that long-wavelength topography is dominantly the result of vertical convective tractions at the base of the lithosphere. The global average crustal thickness is less than 30 km, but can reach values as large as 40 km beneath tesserae and highland plateaus. We also note that an Earth-like radial viscosity structure cannot be rejected by the gravity and topography data and that, without a mechanical model of the lithosphere, admittance values cannot constrain the thickness of the thermal boundary layer of Venus. Modelling the lithosphere as a thin elastic plate indicates that at the time of formation of relief in highland plateaus and tesserae, the effective elastic plate thickness, T_e , was less than 20 km. Estimates of T_e at highland rises are consistently less than 30 km. Our inability to find regions with T_e > 30 km is inconsistent with predictions made by a class of catastrophic resurfacing models.     

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.     

Coercive force of single crystals of magnetite at low temperatures

The temperature dependence of coercive force H _c was studied on well-characterized and stoichiometric millimetre-sized single crystals of magnetite at a series of 16 temperatures from 300 to 10 K using a SQUID magnetometer. H _c decreases gradually with cooling to the isotropic temperature, T _i = 130 K, where the first magnetocrystalline anisotropy constant K ₁ becomes zero. H _c exhibits a sharp increase at the Verwey transition, T _v = 120 K, where the structure changes from cubic to monoclinic. In crossing the Verwey transition, H _c increases by more than two orders of magnitude, from 20 μT to 2.4 mT, and the shape of the hysteresis loops becomes wasp-waisted.
Observed coercivity between 300 K and 170 K varies with temperature as λ _s / M _s , where λ _s is the magnetostriction constant and M _s is the saturation magnetization, indicating that the coercivity in MD magnetite is controlled mainly by internal stress associated with dislocations or other crystal defects. It seems likely that the stable single-domain-like magnetic memory observed in large MD magnetite crystals is due to magnetoelastically pinned domain walls. The discontinuous change in H _c at the Verwey transition is controlled by abrupt changes in magnetocrystalline and magnetostriction constants due to crystal deformation from cubic to monoclinic structure.     

Thrust faulting and crust—upper mantle structure in East Australia

Summary. The mid-crustal earthquake of 1973 March 9 (m_b= 5.5, h ≤ 20 km) located 60 km south-west of Sydney, Australia, provides unambiguous evidence of contemporary thrust faulting in South-eastern Australia — a region of high heat flow and Cenozoic basaltic volcanism. Aftershock locations suggest a steeply dipping fault in the depth range from 8 to 24 km with a lateral extent of about 8 km. The mechanism solution is consistent with a tectonic stress field that is dominated by east—west horizontal compression. A seismic moment of 5.7 ± 10²³± 20 per cent dyne-cm was computed from surface-wave amplitudes. Minimum values of slip and stress drop, 2 cm and 1 bar respectively, were estimated from the moment and a fault size taken' from aftershock locations.
Refinement modelling by a controlled Monte Carlo technique was used to provide unbiased models directly from multimode group velocities. The dispersion of fundamental and higher mode surface waves recorded at the digital high-gain station at Charters Towers, Queensland, and the WWSSN station at Adelaide, South Australia, is satisfied by crust- and upper-mantle models which have neither pronounced S-wave low-velocity zones nor thick high-velocity lids within 140 km of the Earth's surface. These models have subcrustal shear velocities of 4.20–4.32 km/s which are 0.4–0.5 km/s slower than Canadian shield shear velocities (CANSD).     

Magnetic fabric and its significance in the Florianópolis dyke swarm, southern Brazil

Magnetic fabric was determined by applying the anisotropy from the low-field magnetic susceptibility (AMS) technique in 62 mafic dykes from the Mesozoic Florianópolis (Santa Catarina Island) dyke swarm, southern Brazil. These dykes cut the crystalline basement rocks, which are mainly Proterozoic. They are vertical or subvertical in dip and trend mainly NE, although NW-trending dykes are also found. Dykes are tholeiitic in composition and are geochemically similar to those from the Ponta Grossa swarm. Thicknesses vary from 0.3 to 60 m. Polished sections show that titanomagnetites carry the AMS in these dykes. Hysteresis parameters show that the magnetic minerals fall in the PSD range. Two types of magnetic fabric are recognized. Type I is characterized by K ₁- K ₂ parallel to the dyke wall, representing magma flow within the dykes; type II, with K ₁- K ₃ parallel to the dyke wall, was found in four dykes. Type I is found in 94 per cent of the dykes, and approximately 20 per cent of these have K ₁ inclinations of less than 30°, suggesting horizontal or subhorizontal flow. About 80 per cent have K ₁ inclinations of greater than 30°, due to inclined to vertical flow. The comparison of AMS studies from both the Florianópolis and the Ponta Grossa dykes suggests a source position closer to Santa Catarina Island than the Ponta Grossa arch.     

Seismic image of a CO2 reservoir beneath a seismically active volcano

Mammoth Mountain is a seismically active volcano 200 000 to 50 000 years old, situated on the southwestern rim of Long Valley caldera, California. Since 1989 it has shown evidence of unrest in the form of earthquake swarms (Hill et al. 1990), volcanic 'long-period' earthquakes (Pitt & Hill 1994), increased output of magmatic ³He (Sorey et al. 1993) and the emission of about 500 tonnes day ⁻¹ of CO₂ (Farrar et al. 1995; Hill 1996; M. Sorey, personal communication, 1997), which has killed trees and poses a threat to human safety. Local-earthquake tomography shows that in mid-1989 areas of subsequent tree-kill were underlain by extensive regions where the ratio of the compressional and shear elastic-wave speeds V_P/V_S was about 9 per cent lower than in the surrounding rocks. Theory (Mavko & Mukerji 1995), experiment (Ito, DeVilbiss & Nur 1979), and experience at other geothermal/volcanic areas (Julian et al. 1996) and at petroleum reservoirs (Harris et al. 1996) indicate that V_P/V_S is sensitive to pore-fluid compressibility, through its effect on V_P . The observed V_P/V_S anomaly is probably caused directly by CO₂, and seismic V_P/V_S tomography is thus a promising tool for monitoring gas concentration and movement in volcanoes, which may in turn be related to volcanic activity.     

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.     

Approaches to Modelling the Surface Albedo of a High Arctic Glacier

Broadband surface albedo measurements, made during the summer melt season at three weather stations on John Evans Glacier (79°40 ' N, 74°00 ' W), varied strongly with the solar zenith angle ( θ _z ). Tests were carried out to assess the impact of diurnal variations in surface albedo on seasonal net shortwave radiation ( K ^* ) totals. Removing the diurnal signal from albedo measurements by daily averaging of hourly measurements, or by applying midday measurements to all hours of the day, changed K ^* by up to 16%. Ignoring measurements made at θ _z & 75°, to account for measurement (cosine) error at high θ _z , decreased K ^* by between 5 and 18%. Given the sensitivity of K ^* to diurnal patterns in surface albedo, experiments were carried out with two albedo models. One model accounted for albedo variations with θ _z and one did not. The model driven by θ _z , when implemented within a surface energy balance model for John Evans Glacier, produced better melt estimates. This suggests that diurnal variations in surface albedo should be accounted for in energy balance models of glacier melt.     

Tidal gravity and ocean tide loading in Europe

b
The results are presented from tidal gravity measurements at five sites in Europe using LaCoste and Romberg ET gravimeters. Improvements that we have made to the accuracies of these gravimeters are discussed. It is shown that the 'standard' calibration of the International Center for Earth Tides, used for worldwide tidal gravity profiles, is 1.2 per cent too high. The M₂ and O₁ observations are compared with model calculations of the Earth's body tide and ocean tide loading and it is shown that there is a very significant improvement in the agreement between observations and models compared to that obtained with previous tidal gravity measurements. For O₁, where the ocean tide loading and attraction in central Europe is only 0.4 per cent of the body tide, our measurements verify that the Dehant-Wahr anelastic body tide model gravimetric factor is accurate to 0.2 per cent. It is also shown that the effects of lateral heterogeneities in Earth structure on tidal gravity are too small to explain the large anomalies in previously published tidal gravity amplitudes. The observations clearly show the importance of conserving tidal mass in the Schwiderski ocean tide model. For sites in central Europe, the M₂ and O₁ observations and the models are in agreement at the 0.1 μgal (10⁻⁹ m s⁻²) level and tidal corrections to this accuracy can now be made to absolute gravity measurements.     

Observation of Tidal Tilt On Kilauea Volcano, Hawaii

Summary. We have analysed the east-west tilt components, O₁, K₁, N₂, M₂ and S₂ from a continuously recording tiltmeter located in Uwekahuna Vault on Kilauea Volcano, Hawaii, for the period 1971—79. Detailed analysis of the M₂ component gives values of 30.9 ± 2.0 (95 per cent) nrad and 116.0 ± 2.0° for the amplitude and phase, respectively, compared to values of 48.5 nrad and 139.4° for the equilibrium tide. the total theoretical tide, found by summing the equilibrium and load tides, amounts to 37.2 nrad at a phase of 121.7°. the 20 per cent discrepancy with that observed may be due to an inaccurate cotical chart, cavity effects in the vault, strain—tilt coupling or an inappropriate solid earth model. In the vicinity of Hawaii (≤ 3°) two independent cotidal charts give almost identical results for the near field ocean load. At greater distances, we use the Schwiderski (1978) cotidal chart. We estimate that local cavity and strain—tilt coupling effects are less than 12 per cent owing to the agreement between geodetically determined and instrumental tilt but we can not rule out regional effects. Assuming these are small and the cotical charts correct, we find that the M₂ results are brought into satisfactory agreement if, instead of using an average oceanic earth model in the (< 75 km) vicinity of Hawaii, we use an earth model with nearly one-half the oceanic rigidity. Such a low upper mantle and crustal rigidity is consistent with Kilauea's position above the thermal upwelling associated with the Hawaiian hotspot.     

Magnetotelluric response of a uniformly stratified earth containing a magnetized layer

Summary The magnetotelluric (MT) response is studied of a uniformly stratified earth which contains a magnetized layer. The impedance as a function of the layer parameters (resistivity, ρ permeability, μ and thickness h ) is discussed. The MT response from a layer (μ, ρ, h ) is equivalent to that from a layer (μ/μ_r=μ₀, μ_rρ, μ_r h ) where μ_r is the relative permeability of the layer. Thus the effect of a magnetized layer is to make it apppear μ_r times more resistive and μ_r times thicker than an unmagnetized equivalent layer. Master curves of apparent resistivity and phase are computed for three-layer models with varying permeability associated with varying resistivity in each layer. An example of MT field data is presented in which the most reasonable interpretation is that a magnetized layer exists beneath the observatory site.     

3-D seismic velocities calculated from lattice-preferred orientation and reflectivity of a lower crustal section: examples of the Val Sesia section (Ivrea zone, northern Italy)

To quantify the seismic properties of lower crustal rocks and to better constrain the origin of the lower crustal seismic reflectivity, we determined the complete 3-D seismic properties of a lower crustal section. Eight representative samples of the main lithologic and structural units outcropping in the Val Sesia (Ivrea zone) were studied in detail. The seismic velocities were calculated using the single crystal stiffness coefficients and the lattice preferred orientation (LPO) of each mineral in all samples. The 21 stiffness coefficients characterizing the elastic behaviour of each rock are determined. Mafic and ultramafic rocks such as pyroxenite and pyroxene-bearing gabbros display complex shear wave properties. These rocks are weakly birefringent (maximum 0.1 kms⁻¹) and it is difficult to find consistent relationships between the seismic properties and the rock structure. On the other hand, seismic properties of deformed felsic rocks are essentially controlled by mica. They display strong S -wave birefringence (0.3 km s⁻¹) and relatively high V _p anisotropy (7.6 per cent). Amphibole also strongly influences the rock birefringence patterns. For both kind of rocks, the foliation is highly birefringent and the fast polarized shear wave is systematically oriented parallel to the foliation. We show that the number of mineral phases in the rock strongly controls the anisotropy. The seismic anisotropy has a complex role in the P -wave reflectivity. Compared to the isotropic case, anisotropy enhances the reflection coefficient for about 60 per cent of the possible lithological interfaces. For 40 per cent of the interfaces, the reflection coefficient is much lower when one considers the medium as anisotropic.     

Error estimates in the measurement of anisotropic magnetic susceptibility

Hext's (1963 ) first-order analysis for error estimation in measurements of second-rank tensors is applied to the torque meter and the orthogonal-coil induction instrument, both representatives of the class of instruments that measures differences in magnetic susceptibility. The conventional measurement design for the torque meter is rotatable; an efficient and nearly rotatable design is proposed for the orthogonal-coil instrument. Error estimates are derived for a set of anisotropy parameters, K _mean , H and μ . The applicability of a first-order analysis is discussed and illustrated by simulation studies.     

Seismic and electromagnetic evidence of dehydration as a free water source in the reactivated crust

According to recent estimates, the continental mid-crust contains 35–40 per cent amphibolites. Heating of the crust by an underlying mantle plume, for example beneath continental rifts, high plateaus, and areas of intraplate volcanic activity, releases water. Dehydration of amphibole-bearing rocks at depths of 20–40  km occurs mainly in the temperature range 650–700 °C, and this releases about 0.4  wt per cent of water.
  Seismic tomography studies of the crust in the Kirgyz Tien Shan Range, where the age of the tectonic activity is less than 30  Ma, revealed a low-velocity zone in the mid-crust. The velocity of P waves was 0.4  km  s⁻¹ lower than in normal crust. MT sounding data in the region show the existence of a low-resistivity layer with an average resistivity of about 25  Ω  m at the depth of the low-velocity layer. The spatial correlation of the observed anomalous layers and calculated effect of fluid phase on seismic and electric parameters of rocks suggests the presence of aqueous fluids released by the heating of the mid-crust.     

Upper mantle shear structure of North America

Summary. The waveforms and travel times of S and SS phases in the range 10°–60° have been used to derive upper mantle shear velocity structures for two distinct tectonic provinces in North America. Data from earthquakes on the East Pacific Rise recorded at stations in western North America were used to derive a tectonic upper mantle model. Events on the north-west coast of North America and earthquakes off the coast of Greenland provided the data to investigate the upper mantle under the Canadian shield. All branches from the triplications due to velocity jumps near 400 and 660 km were observed in both areas. Using synthetic seismograms to model these observations placed tight constraints on heterogeneity in the upper mantle and on the details of its structure. SS–S travel-time differences of 30 s along with consistent differences in waveforms between the two data sets require substantial heterogeneity to at least 350 km depth. Velocities in the upper 170 km of the shield are about 10 per cent higher than in the tectonic area. At 250 km depth the shield velocities are still greater by about 4.5 per cent and they gradually merge near 400 km. Below 400 km no evidence for heterogeneity was found. The two models both have first-order discontinuities of 4.5 per cent at 405 km and 7.5 per cent at 695 km. Both models also have lids with lower velocities beneath. In the western model the lid is very thin and of relatively low velocity. In the shield the lid is 170 km thick with very high elocity (4.78 km s^-1); below it the velocity decreases to about 4.65 km s^-1. Aside from these features the models are relatively smooth, the major difference between them being a larger gradient in the tectonic region from 200 to 400 km.     

Upper crustal shear velocity models from higher mode Rayleigh wave dispersion in Scotland

Summary. Group velocities for first and second higher mode Rayleigh waves, in the frequency range 0.8–4.8 Hz, generated from a local earthquake of magnitude 3.7 M _L in western Scotland, are measured at stations along the 1974 LISPB line. These provide detailed information about the crustal structure west of the line. The data divide the region into seven apparently homogeneous provinces. Averaged higher mode velocity dispersion curves for each province are analysed simultaneously using a linearized inversion technique, yielding regionalized shear velocity profiles down to a depth of 17 km into the upper crust. Shear wave velocity is between 3.0 and 3.4 km s⁻¹ in the upper 2 km, with a slow increase to around 3.8 km s⁻¹. P -wave models computed using these results agree with profiles from the LISPB and LUST refraction experiments.