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.     

Seismic anisotropy as measured under high-pressure, high-temperature conditions

Summary. Four types of crustal and upper-mantle rocks have been used for the investigation of seismic P -wave velocities in three mutually perpendicular directions. Hydrostatic pressure, up to 6 kbar and temperatures up to 500°C were applied to the samples. Measurements of the ultrasonic P -wave travel times and velocities were carried out along two geotherms. All rock types show an anisotropic behaviour which is caused by the orientation of certain minerals. The anisotropy is not dependent on temperature and pressure. Gneiss and peridotite have 5–6 per cent anisotropy whereas granite and a metagabbro show values of only 2–3 per cent. The smallest velocity is always in the z direction, perpendicular to a schistocity or foliation. It is shown that the data agree with those of field observation. We conclude that anisotropy caused by preferred orientation of minerals must be expected in the whole lithosphere. Additional effects of layering, of cracks, and of nonhydrostatic stresses are estimated.     

Waveform modelling of teleseismic S, Sp, SsPmP, and shear-coupled PL waves for crust- and upper-mantle velocity structure beneath Africa

We describe a waveform modelling technique and demonstrate its application to determine the crust- and upper-mantle velocity structure beneath Africa. Our technique uses a parallelized reflectivity method to compute synthetic seismograms and fits the observed waveforms by a global optimization technique based on a Very Fast Simulated Annealing (VFSA). We match the S , Sp, SsPmP and shear-coupled PL phases in seismograms of deep (200–800 km), moderate-to-large magnitude (5.5–7.0) earthquakes recorded teleseismically at permanent broad-band seismic stations in Africa. Using our technique we produce P - and S -wave velocity models of crust and upper mantle beneath Africa. Additionally, our use of the shear-coupled PL phase, wherever observed, improves the constraints for lower crust- and upper-mantle velocity structure beneath the corresponding seismic stations. Our technique retains the advantages of receiver function methods, uses a different part of the seismogram, is sensitive to both P - and S -wave velocities directly, and obtains helpful constraints in model parameters in the vicinity of the Moho. The resulting range of crustal thicknesses beneath Africa (21–46 km) indicates that the crust is thicker in south Africa, thinner in east Africa and intermediate in north and west Africa. Crustal P - (4.7–8 km s⁻¹) and S -wave velocities (2.5–4.7  km s⁻¹) obtained in this study show that in some parts of the models, these are slower in east Africa and faster in north, west and south Africa. Anomalous crustal low-velocity zones are also observed in the models for seismic stations in the cratonic regions of north, west and south Africa. Overall, the results of our study are consistent with earlier models and regional tectonics of Africa.     

The Agulhas Plateau: structure and evolution of a Large Igneous Province

Large Igneous Provinces (LIP) are of great interest due to their role in crustal generation, magmatic processes and environmental impact. The Agulhas Plateau in the southwest Indian Ocean off South Africa has played a controversial role in this discussion due to unclear evidence for its continental or oceanic crustal affinity. With new geophysical data from seismic refraction and reflection profiling, we are able to present improved evidence for its crustal structure and composition. The velocity–depth model reveals a mean crustal thickness of 20 km with a maximum of 24 km, where three major units can be identified in the crust. In our seismic reflection records, evidence for volcanic flows on the Agulhas Plateau can be observed. The middle crust is thickened by magmatic intrusions. The up to 10 km thick lower crustal body is characterized by high seismic velocities of 7.0–7.6 km s⁻¹. The velocity–depth distribution suggests that the plateau consists of overthickened oceanic crust similar to other oceanic LIPs such as the Ontong-Java Plateau or the northern Kerguelen Plateau. The total volume of the Agulhas Plateau was estimated to be 4 × 10⁶ km³ of which about 10 per cent consists of extruded igneous material. We use this information to obtain a first estimate on carbon dioxide and sulphur dioxide emission caused by degassing from this material. The Agulhas Plateau was formed as part of a larger LIP consisting of the Agulhas Plateau itself, Northeast Georgia Rise and Maud Rise. The formation time of this LIP can be estimated between 100 and 94 (± 5) Ma.     

Ivrea seismic array: a study of continental crust and upper mantle

A seismic-array study of the continental crust and upper mantle in the Ivrea-Yerbano and Strona-Ceneri zones (northwestern Italy) is presented. A short-period network is used to define crustal P - and S -wave velocity models from earthquakes. The analysis of the seismic-refraction profile LOND of the CROP-ECORS project provided independent information and control on the array-data interpretation.
Apparent-velocity measurements from both local and regional earthquakes, and time-term analysis are used to estimate the velocity in the lower crust and in the upper mantle. The geometry of the upper-lower crust and Moho boundaries is determined from the station delay times.
We have obtained a three-layer crustal seismic model. The P -wave velocity in the upper crust, lower crust and upper mantle is 6.1±0.2 km s⁻¹, 6.5±0.3 km s⁻¹ and 7.8±0.3 km s⁻¹ respectively. Pronounced low-velocity zones in the upper and lower crust are not observed. A clear change in the velocity structure between the upper and lower crust is documented, constraining the petrological interpretation of the Ivrea-type reflective lower continental crust derived from small-scale petrophysical data. Moreover, we found a V _P/ V _S ratio of 1.69±0.04 for the upper crust and 1.82±0.08 for the lower crust and upper mantle. This is consistent with the structural and petrophysical differences between a compositionally uniform and seismically transparent upper crust and a layered and reflective lower crust. The thickness of the lower crust ranges from about 8 km in front of the Ivrea body (ARVO, Arvonio station) in the northern part of the array to a maximum of about 15 km in the southern part of the array. The lower crust reaches a minimum depth of 5 km below the PROV (Provola) station.     

First-order P-wave ray synthetic seismograms in inhomogeneous weakly anisotropic media

We propose approximate equations for P -wave ray theory Green's function for smooth inhomogeneous weakly anisotropic media. Equations are based on perturbation theory, in which deviations of anisotropy from isotropy are considered to be the first-order quantities. For evaluation of the approximate Green's function, earlier derived first-order ray tracing equations and in this paper derived first-order dynamic ray tracing equations are used.
The first-order ray theory P -wave Green's function for inhomogeneous, weakly anisotropic media of arbitrary symmetry depends, at most, on 15 weak-anisotropy parameters. For anisotropic media of higher-symmetry than monoclinic, all equations involved differ only slightly from the corresponding equations for isotropic media. For vanishing anisotropy, the equations reduce to equations for computation of standard ray theory Green's function for isotropic media. These properties make the proposed approximate Green's function an easy and natural substitute of traditional Green's function for isotropic media.
Numerical tests for configuration and models used in seismic prospecting indicate negligible dependence of accuracy of the approximate Green's function on inhomogeneity of the medium. Accuracy depends more strongly on strength of anisotropy in general and on angular variation of phase velocity due to anisotropy in particular. For example, for anisotropy of about 8 per cent, considered in the examples presented, the relative errors of the geometrical spreading are usually under 1 per cent; for anisotropy of about 20 per cent, however, they may locally reach as much as 20 per cent.     

Wave Propagation In A Vertical Transversely Isotropic Medium: Field Experiment and Model Study

Summary. The continent-ocean transition adjacent to Hatton Bank was studied using a dense grid of single-ship and two-ship multichannel seismic profiles. The interpretation of the explosive expanding spread profiles (ESPs) which were shot as part of this survey are discussed here in detail. Extensive seaward dipping reflectors are developed in the upper crust across the entire margin. These seaward dipping reflectors continue northwards on the Faeroes and Vøring margins, where they have been shown to be caused by basaltic lavas, as well as on the conjugate margin of East Greenland. The dipping reflectors are an important feature of the rifting history of the margin and show that extensive volcanism was associated with the extension. The ESPs show clear seismic arrivals out to ranges of 100 km. Wide-angle Moho reflections can be seen on all the lines as well as good mid and lower crustal arrivals. The determination of seismic velocity structure was constrained by ray tracing and by amplitude modelling using reflectivity synthetic seismograms. The results from the ESPs show that there is a thick region of lower crustal material beneath the margin with an unusually high crustal velocity of 7.3–7.4 km s⁻¹. This lower crustal material reaches a maximum thickness of 14 km beneath the central part of the margin and is terminated at depth by the Moho. The lower crustal lens of high-velocity material is interpreted as underplated or intruded igneous rocks associated with the large volumes of extrusive basaltic lavas, now seen as dipping reflectors on the margin.     

Interpretation of subhorizontal crustal reflections by metamorphic and rheologic effects in the eastern part of the Pannonian Basin

The geologic origin of subhorizontal reflections, often observed in crustal seismic sections, was investigated by establishing metamorphic facies and strength of rocks in depth, and correlating these properties to seismic reflection sections from eastern Hungary. Estimation of the depths of metamorphic mineral stability zones utilized the principles developed by Fyfe et al. and known geothermal data of the area. The strength versus depth profile was derived by relating local seismic P -wave interval velocities to Meissner et al. 's activation energy. The results show that the series of subhorizontal reflections, observed in the Pannonian Basin, are a consequence of combined metamorphic and rheologic changes in depths. The synthesis of the integrated data set suggests that the retrograde alteration of the pre-Tertiary basement above the percolation threshold was made possible by the softening effect of shear zones and their water-conducting capacity. The subhorizontal reflections of highest energy, of the consolidated crust below the percolation threshold, originate in the depths of greenschist, amphibolite and granulite metamorphic mineral facies, which were formed in geothermal and pressure conditions similar to those existing today. These results imply the overprint of earlier (Variscan) metamorphic sequences of the crust by more recent retrograde metamorphic processes.     

P-wave velocity anisotropy in crystalline rocks

Summary. Compressional wave velocities and anisotropy coefficients determined at high hydrostatic pressures are compiled from the data published for the main types of crystalline rocks. The crack-free elastic anisotropy of igneous crustal rocks is generally very low, between 1 and 3 per cent on average. The anisotropy of metamorphic rocks is higher (up to 22 per cent), but very variable. The average anisotropy coefficients in schists and amphibolites are about 10 per cent, in gneisses between 3 and 7 per cent, and in granulites less than 3 per cent. The average anisotropy of olivine ultramafites is between 7 and 12 per cent, whereas in pyroxenites and eclogites it is usually less than 4 per cent. A comparison of ranges of average velocities and average anisotropies for the individual rock groups suggests that, whereas in the crust the lateral velocity variations are mainly due to compositional changes, in the olivine of the uppermost mantle the velocity variations due to anisotropic structures could be of the same magnitude as the variations due to inhomogeneities.     

Detailed S-wave structure in the Dublin Basin and its northern margin

Summary. The seismic structure has been measured to a depth of about 3 km along a 30 km seismic profile in east central Ireland. This profile is unusual in that it is the S -wave velocity—depth structure that has been measured to a degree of precision more normally associated with P -wave results. One reason for this is that the sources used were quarry blasts which generated strong S -waves and short-period surface waves but rather weak P -waves.
The results show a layer of Carboniferous limestone with shear velocity 2.65 km⁻¹ s overlying a layer with a velocity of 3.06 km s⁻¹. This second layer was interpreted as Lower Palaeozoic strata (Silurian/Ordovician) since this velocity was evident in an inlier seen at the surface at the northern end of the line. A third refraction horizon, shear velocity 3.45 km s⁻¹ and displaying a basinal structure, was also recognized. This may be Cambrian or Precambrian basement.     

Phase-array decomposition of a seismic section

A seismic re fraction/wide-angle reflection profile is analysed for the presence of correlated events ('phases'). The correlation problem is formulated in terms of temporally, spatially and frequency-local complex covariances. For robustness, the method concentrates on phase rather than amplitude information. This allows a computationally efficient algorithm that can make allowance for signal correlation length and can model curved wavefronts. A statistical test based on residual phase misfit across the analysed subarray is used to assess the probability that a detected event represents a real correlated signal.
With our chosen analysis parameters and confidence level (over 99.9 per cent). 1222 events were detected in the data. Using simple techniques based on 1-D earth models, detected events are associated with a small number of particular wave types. In this way, we have succeeded in classifying almost 95 per cent of the detected events. Those that remain describe those components of the data that are inconsistent with our simple ray paths in the 1-D assumption and with our prescribed tolerance. These include reverberations, near-surface guided waves and reflected waves from strongly laterally inhomogeneous structures. According to our modelling, about 25 per cent of the detected events are consistent with simple P -wave reflected energy, and these are to a very large extent (over 85 per cent) distinct from all the other wave-type models we have used. A direct mapping of the detected events into the offset-depth domain reveals dear internal and external consistencies among the detections for the various wave types. Estimated earth structure is consistent with models from previous analyses based on much larger data sets.
We have thus succeeded in extracting correlated events from the data and decomposing these, approximately but meaningfully, into distinct classes (ray paths)     

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.     

Shallow velocity structure along the Hirapur–Mandla profile using traveltime inversion of wide-angle seismic data, and its tectonic implications

In order to investigate the velocity structure, and hence shed light on the related tectonics, across the Narmada–Son lineament, traveltimes of wide-angle seismic data along the 240 km long Hirapur–Mandla profile in central India have been inverted. A blocky, laterally heterogeneous, three-layer velocity model down to a depth of 10 km has been derived. The first layer shows a maximum thickness of the upper Vindhyans (4.5 km s⁻¹ ) of about 1.35 km and rests on top of normal crystalline basement, represented by the 5.9 km s⁻¹ velocity layer. The anomalous feature of the study is the absence of normal granitic basement in the great Vindhyan Graben, where lower Vindhyan sediments (5.3 km s⁻¹ ) were deposited during the Precambrian on high-velocity (6.3 km s⁻¹ ) metamorphic rock. The block beneath the Narmada–Son lineament represents a horst feature in which high-velocity (6.5 km s⁻¹ ) lower crustal material has risen to a depth of less than 2 km. South of the lineament, the Deccan Traps were deposited on normal basement during the upper Cretaceous period and attained a maximum thickness of about 800 m.     

Finite-frequency traveltime tomography of San Francisco Bay region crustal velocity structure

Seismic velocity structure of the San Francisco Bay region crust is derived using measurements of finite-frequency traveltimes. A total of 57 801 relative traveltimes are measured by cross-correlation over the frequency range 0.5–1.5 Hz. From these are derived 4862 'summary' traveltimes, which are used to derive 3-D P -wave velocity structure over a 341 × 140 km² area from the surface to 25 km depth. The seismic tomography is based on sensitivity kernels calculated on a spherically symmetric reference model. Robust elements of the derived P -wave velocity structure are: a pronounced velocity contrast across the San Andreas fault in the south Bay region (west side faster); a moderate velocity contrast across the Hayward fault (west side faster); moderately low velocity crust around the Quien Sabe volcanic field and the Sacramento River delta; very low velocity crust around Lake Berryessa. These features are generally explicable with surface rock types being extrapolated to depth ∼10 km in the upper crust. Generally high mid-lower crust velocity and high inferred Poisson's ratio suggest a mafic lower crust.     

Transitional continental–oceanic structure beneath the Norwegian Sea from inversion of surface wave group velocity data

We have analysed the fundamental mode of Love and Rayleigh waves generated by 12 earthquakes located in the mid-Atlantic ridge and Jan Mayen fracture zone. Using the multiple filter analysis technique, we isolated the Rayleigh and Love wave group velocities for periods between 10 and 50  s. The surface wave propagation paths were divided into five groups, and average group velocities calculated for each group. The average group velocities were inverted and produced shear wave velocity models that correspond to a quasi-continental oceanic structure in the Greenland–Norwegian Sea region. Although resolution is poor at shallow depth, we obtained crustal thickness values of about 18  km in the Norwegian Sea area and 9  km in the region between Svalbard and Iceland. The abnormally thick crust in the Norwegian Sea area is ascribed to magmatic underplating and the thermal blanketing effect of sedimentary layers. Maximum crustal shear velocities vary between 3.5 and 3.9  km  s⁻¹ for most paths. An average lithospheric thickness of 60  km was observed, which is lower than expected for oceanic-type structure of similar age. We also observed low shear wave velocities in the lower crust and upper mantle. We suggest that high heat flow extending to depths of about 30  km beneath the surface can account for the thin lithosphere and observed low velocities. Anisotropy coefficients of 1–5 per cent in the shallow layers and >7 per cent in the upper mantle point to the existence of polarization anisotropy in the region.     

Precise measurement of seismic travel times from an air gun source

Summary. Seismic travel times for extrema, zero-crossings, or entire body waves need to be determined precisely to one part in 10³ or better in several varieties of seismic studies employing an impulsive artificial source. Examples are crosshole surveys which delineate rock crack distribution separating the holes and monitoring of crustal seismic travel times in earthquake precursor studies. A timing resolution of one part in 10³ has been achieved previously using digitally recorded seismic data. These methods, however, do not use interpolation between digitized data points as a method to increase the timing resolution. We report travel-time determinations based on interpolation between digitized points which achieve a precision of two parts in 10⁴, a five-fold improvement over the existing methods. In addition, the effects of seismic noise on travel-time measurement have been compared for the extremum location, the unnormalized correlation, and the normalized correlation method. The following conclusions are drawn from this comparison: (1) the normalized correlation method provides an 18–55 per cent improvement in the standard deviation of the mean over the extremum location method, and (2) results as accurate as those by the normalized correlation can be obtained by the unnormalized correlation if a complete up-and-down swing of the waveform is used as the master trace and if the master trace is close to being sinusoidal. The advantage of the unnormalized correlation over the normalized correlation is speed; the unnormalized correlation is faster by a factor of 28 in computing time.     

Crustal structure from an OBS survey of the Nootka fault zone off western Canada

Summary. The Nootka fault zone is the boundary between the small Explorer and Juan de Fuca plates which are situated between the America and Pacific plates off western Canada. To investigate the crustal structure in the region, three explosive/large airgun refraction lines were shot into three ocean bottom seismometers (OBSs) with three-component geophone assemblies. In this phase of the study, P -wave velocity—depth models are interpreted by comparison of the travel time and amplitude characteristics of the observed data with theoretical seismograms computed using a WKBJ algorithm. The interpretation gives relatively consistent results for the upper crust. However, the structure of the lower crust is significantly different among the various profiles. Upper mantle velocities range from 7.5 to 8.3 kms⁻¹ and the sub-bottom crustal thickness vanes from 6.4 to 11 km. Nevertheless, these seismic models are consistent in general terms with oceanic crustal models represented by ophiolite complexes. Some aspects of the differences among profiles can be explained by consideration of a recent tectonic model for the development of the fault zone. This requires, within a 1 Myr time interval, variations in the process of crustal formation at the ridge, crustal 'maturing', or both. The abnormally thick crust near a spreading centre may result in part from the complex interaction of the Juan de Fuca and Explorer plates with the larger and older America and Pacific plates. Upper mantle velocity variations are consistent with the concept of velocity anisotropy. The different record sections show that seismic energy is attenuated for ray paths traversing the Nootka fault zone.     

Anisotropy in multi-offset deep-crustal seismic experiments

Modelling of deep-seismic wide-angle data commonly assumes that the Earth is heterogeneous and isotropic. It is important to know the magnitudes of errors that may be introduced by isotropic-based wide-angle models when the Earth is anisotropic. It is equally important to find ways of detecting anisotropy and determining its properties.
  This paper explores the errors introduced by interpreting anisotropic seismic data with isotropic models. Errors in P -wave reflector depths are dependent on the magnitude of the velocity anisotropy and the direction of the fast axis. The interpreted, isotropic, model velocity function is found to correspond closely to the horizontal velocity of the anisotropic medium. An additional observed parameter is the time mismatch , which we define to be the difference between the vertical two-way traveltime to a reflector and the time-converted wide-angle position of the reflector. The magnitude of the time mismatch is typically <1.0  s (when the whole crust is anisotropic) and is found to be closely related to the magnitude and sign of the anisotropic anellipticity. The relationships are extendible to more complicated models, including those with vertical velocity gradients, crustal zonation, and lower symmetry orders.
  A time mismatch may be symptomatic of the presence of anisotropy. We illustrate the observation of a time mismatch for a real multi-offset seismic data set collected north of Scotland and discuss the implications for crustal anisotropy in that region.     

The Hatton Bank continental margin—I. Shallow structure from two-ship expanding spread seismic profiles

Summary. The Hatton Bank passive continental margin exhibits thick seaward dipping reflector sequences which consist of basalts extruded during rifting between Greenland and Rockall Plateau. Multichannel seismic reflection profiling across the margin reveals three reflector wedges with a maximum thickness near 7 km, extending from beneath the upper continental slope to the deep ocean basin. We present results of the velocity structure within the dipping reflector sequences at eight locations across the margin, interpreted by synthetic seismogram modelling a set of multichannel expanding spread profiles parallel to the margin. At the top of some reflector sequences, we observe a series of 100 m thick high- and low-velocity zones, which are interpreted as basalt flows alternating with sediments or weathered and rubble layers. At the profile locations, the base of the dipping reflectors correlates with P -wave velocities near 6.5 km s⁻¹. However, elsewhere the reflectors appear to extend significantly deeper than the inferred 6.5 km s⁻¹ velocity contour, indicating that the velocity structure may not be controlled solely by lithological boundaries but also by metamorphic effects. Shear-waves were observed on two lines, permitting the calculation of Poisson's ratio. The decrease in Poisson's ratio from 0.28 to near 0.25 in the upper 5 km of crust may also indicate the effect of metamorphism on seismic properties, or alternatively may be explained by crack closure under load.     

Anisotropy in the oceanic lithosphere of the North-western Pacific Basin

Summary. A series of long-range explosion seismological experiments has been conducted by the use of specially designed ocean bottom seismographs (OBSs) in the Western Pacific. OBS studies of apparent velocity measurements by the use of natural earthquakes have also been made. The experiments have made clear that large-scale P -wave anisotropy exists in the entire thickness of the oceanic lithosphere. The existence of the large-scale anisotropy in the oceanic lithosphere has been demonstrated for the first time by seismic body-wave studies. Previously, anisotropy had been found only in the uppermost oceanic mantle in the Eastern Pacific.
The azimuth of the maximum velocity, 8.6 km s^-1, is about 155° clock-wise from north. The direction is perpendicular to the magnetic lineation of the region, however, the direction differs from the direction of the present plate motion by about 30°. So it appears that the anisotropy has been 'frozen' at least since the change of the plate motion that occurred 40 Myr ago. The frozen anisotropy should set important constraints on the mechanical properties of the lithosphere such as the viscosity and temperature of the lower lithosphere.