Upper mantle structure beneath the Mid-Atlantic Ridge north of the Azores based on observations of compressional wave

Summary. Four seismic refraction profiles have been interpreted which serve to indicate the structure of the lithosphere near the Mid-Atlantic Ridge close to the Azores. An east–west profile which crosses the ridge axis yields a crustal structure. Although energy is propagated across the ridge axis within the crust the axial region marks a clear barrier to propagation within the mantle. A profile parallel to the axis (4 my isochron) shows, below a 7.6 km/s layer, a low-velocity zone underlain by an 8.3 km/s refractor 9 km below the sea bed. On profies normal to the ridge axis higher velocities, which are observed on lines shot towards the ridge, can be attributed to this refractor if it has a dip of several degrees away from the ridge. On another profile parallel to the axis (9 my isochron) a velocity of about 8.3 km/s is only found to exist much deeper at about 30 km depth. These observations are interpreted in the light of seismic refraction results recently obtained by Lewis & Snydsman and of quantitative petrological models, such as that of Bottinga & Allègre. A velocity model based on Bottinga & Allègre's model allows us to understand our results qualitatively. In particular the two 8.3 km/s refractors at 9 and 30 km depth correspond to two different residual peridotite layers. The upper layer contains 1.5–2 wt per cent water and as the lithosphere moves away from the ridge axis the temperature in this layer becomes low enough to start hydration reactions. These cause the low-velocity zone observed at 4 my and the total disappearance of the shallow level refractor before 9 my.     

Crustal tomographic imaging of a transitional continental rift: the Ethiopian rift

In this study we image crustal structure beneath a magmatic continental rift to understand the interplay between crustal stretching and magmatism during the late stages of continental rifting: the Main Ethiopian Rift (MER). The northern sector of this region marks the transition from continental rifting in the East African Rift to incipient seafloor spreading in the southern Red Sea and western Gulf of Aden. Our local tomographic inversion exploits 172 broad-band instruments covering an area of 250 × 350 km of the rift and adjacent plateaux. The instruments recorded a total of 2139 local earthquakes over a 16-month period. Several synthetic tests show that resolution is good between 12 and 25 km depth (below sea level), but some horizontal velocity smearing is evident along the axis of the Main Ethiopian Rift below 16 km. We present a 3-D P -wave velocity model of the mid-crust and present the first 3-D Vp / Vs model of the region. Our models show high P -wave velocities (6.5 km s⁻¹) beneath the axis of the rift at a depth of 12–25 km. The presence of high Vp / Vs ratios (1.81–1.84) at the same depth range suggest that they are cooled mafic intrusions. The high Vp / Vs values, along with other geophysical evidence, suggest that dyking is pervasive beneath the axis of the rift from the mid-crustal depths to the surface and that some portion of partial melt may exist at lower crustal depths. Although the crustal stretching factor across the Main Ethiopian Rift is ∼1.7, our results indicate that magma intrusion in narrow zones accommodates a large proportion of extensional strain, with similarities to slow-spreading mid-ocean ridge processes.     

The Mid-Atlantic Ridge: structure at 45° N

Summary. A structural model of the Mid-Atlantic Ridge at 45° N is proposed on the basis of travel-time data, amplitudes and synthetic seismograms. The crustal structure seems to be similar to that in the FAMOUS area (Fowler). At the ridge axis there is an absorptive zone in the upper mantle, the depth below the seabed to the top of this zone being about 6 km. Away from the ridge axis there is a positive velocity gradient of about 0.04 to 0.05 km/(skm) in the top 5 to 8 km of the upper mantle. Shear waves propagate across the ridge axis, suggesting that there is no sizeable crustal magma chamber. The shear-wave velocity of the uppermost mantle is 4.35 km/s.     

Crustal structure of the Mid-Atlantic ridge crest at 37° N

Summary. A structural model of the Mid-Atlantic Ridge at 37° N is proposed on the basis of travel-time data and synthetic seismograms. At the ridge axis the crust is only 3 km thick and overlies material with an anomalously low'upper mantle'velocity of 7.2 km s⁻¹. Crustal thickening and the formation of layer 3 and a layer with velocity 7.2–7.3 km s⁻¹ takes place within a few kilometres of the axis, producing a 6–7 km thick crust by less than 10 km from the axis. A normal upper mantle velocity of 8.1 km s⁻¹ exists within 10 km of the axis. Shear waves propagate across the axis, thus precluding the existence of any sizeable magma chamber at shallow depth.     

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.     

Global mapping of upper mantle reflectors from long-period SS precursors

Long-period precursors to SS resulting from underside reflections off upper mantle discontinuities ( SdS where d is the discontinuity depth) can be used to map the global distribution and depth of these reflectors. We analyse 5,884 long-period seismograms from the Global Digital Seismograph Network (1976-1987, shallow sources, transverse component) in order to identify SdS arrivals. Corrections for velocity dispersion, topography and crustal thickness at the SS bounce point, and lateral variation in mantle velocity are critical for obtaining accurate estimates of discontinuity depths. The 410 and 660 km discontinuities are observed at average depths of 413 and 653 km, and exhibit large-scale coherent patterns of topography with depth variations up to 40 km. These patterns are roughly correlated with recent tomographic models, with fast anomalies in the transition zone associated with highs in the 410 km discontinuity and lows in the 660 km discontinuity, a result consistent with laboratory measurements of Clapeyron slopes for the appropriate phase changes. The best resolved feature in these maps is a trough in the 660 km discontinuity in the northwest Pacific, which appears to be associated with the subduction zones in this region. Amplitude variations in SdS arrivals are not correlated with discontinuity depths and probably result from focusing and defocusing effects along the ray paths. The SdS arrivals suggest the presence of regional reflectors in the upper mantle above 400 km. However, only the strongest of these features are above probable noise levels due to sampling inadequacies.     

Slow-ridge/hotspot interactions from global gravity, seismic tomography and 87Sr/86Sr isotope data

Among the mantle hotspots present under oceanic areas, a large number are located on—or close to—active oceanic ridges. This is especially true in the slow-spreading Atlantic and Indian oceans. The recent availability of worldwide gravity grids and the increasing coverage of geochemical data sets along active spreading centres allow a fruitful comparison of these data with global geoid and seismic tomography models, and allow one to study interactions between mantle plumes and active slow-spreading ridges. The observed correlations allow us to draw preliminary conclusions on the general links between surficial processes, which shape the detailed morphology of the ridge axes, and deeper processes, active in the upper mantle below the ridge axial domains as a whole. The interactions are first studied at the scale of the Atlantic (the Mid-Atlantic Ridge from Iceland to Bouvet Island) from the correlation between the zero-age free-air gravity anomaly, which reflects the zero-age depth of the ridge axis, and Sr isotopic ratios of ridge axis basalts. The study is then extended to a more global scale (the slow ridges from Iceland to the Gulf of Aden) by including geoid and upper-mantle tomography models. The interactions appear complex, ranging from the effect of large and very productive plumes, almost totally overprinting the long-wavelength segmentation pattern of the ridge, to that of weaker hotspots, barely marking some of the observables in the ridge axial domain. Intermediate cases are observed, in which hotspots of medium activity (or whose activity has gradually decreased) located at some distance from the ridge axis produce geophysical or geochemical signals whose variation along the axis can be correlated with the geometry of the plume head in the upper mantle. Such observations tend to preclude the use of a single hotspot/ridge interaction model and stress the need for additional observations in various plume/ridge configurations.     

Anisotropy of the mantle inferred from observations of P to S converted waves

Summary. Seismic anisotropy has been previously studied at depths usually not exceeding 100 or 150 km. In this paper we present a method of analysis of seismic records which is very sensitive to azimuthal anisotropy and is applicable at almost any depth range. The idea of the method is to detect and analyse the SH -component of the waves, converted from P to S in the mantle. The procedure of record processing includes frequency filtering, axis rotation, transformation of the record to a standard form, stacking the standardized SH -component records of many seismic events, and the harmonic analysis of amplitude as a function of the direction of wave propagation. When applied to the long-period records of NORSAR the procedure detected a converted wave with the properties implying the possibility of its propagation in a transversely isotropic medium with a horizontal axis of symmetry . Our preferred model postulates anisotropy of ∼ 1 per cent in a layer 50 km thick at the base of the upper mantle.     

Some ocean-bottom seismograph observations on the Reykjanes Ridge at 59° N

Summary. Epicentres of microearthquakes detected by two ocean-bottom seismographs deployed on the Reykjanes Ridge at latitude 59° N can be interpreted as trending between 010° and 020° E of N, parallel to detailed morphological features seen in the area, rather than the overall trend of the ridge (036° E of N). Travel times from an airgun source to the seismographs indicate much slower P wave velocities for the top 4 km at the ridge crest compared with that 5–10 km off axis.     

Constraints on a plume in the mid-mantle beneath the Iceland region from seismic array data

Teleseismic P waves passing through low-wave-speed bodies in the mantle are refracted, causing anomalies in their propagation directions that can be measured by seismometer arrays. Waves from earthquakes in the eastern Pacific and western North America arriving at the NORSAR array in Norway and at seismic stations in Scotland pass beneath the Iceland region at depths of ∼ 1000–2000 km. Waves arriving at NORSAR have anomalous arrival azimuths consistent with a low-wave-speed body at a depth of ∼ 1500 km beneath the Iceland–Faeroe ridge with a maximum diameter of ∼250 km and a maximum wave-speed contrast of ∼ 1.5 per cent. This agrees well with whole-mantle tomography results, which image a low-wave-speed body at this location with a diameter of ∼ 500 km and a wave-speed anomaly of ∼ 0.5 per cent, bearing in mind that whole-mantle tomography, because of its limited resolution, broadens and weakens small anomalies. The observations cannot resolve the location of the body, and the anomaly could be caused in whole or in part by larger bodies farther away, for example by a body imaged beneath Greenland by whole-mantle tomography.     

Crustal model of the Bransfield Rift, West Antarctica, from detailed OBS refraction experiments

The first detailed deep seismic refraction study in the Bransfield Strait, West Antarctica, using sensitive OBSs (ocean bottom seismographs) was carried out successfully during the Antarctic summer of 1990/1991. The experiment focused on the deep crustal structure beneath the axis of the Bransfield Rift. Seismic profile DSS-20 was located exactly in the Bransfield Trough, which is suspected to be a young rift system. Along the profile, five OBSs were deployed at spacings of 50-70 km. 51 shots were fired along the 310 km profile. This paper gives the first presentation of the results. A detailed model of the crustal structure was obtained by modelling the observed traveltimes and amplitudes using a 2-D ray-tracing technique. The uppermost (sedimentary?) cover, with velocities of 2.0-5.5 km s⁻¹, reaches a depth of up to 8 km. Below this, a complex with velocities of 6.4-6.8 km s⁻¹ is observed. The presence of a high-velocity body, with V _p= 7.3-7.7 km s⁻¹, was detected in the 14-32 km depth range in the central part of the profile. These inhomogeneities can be interpreted as a stage of back-arc spreading and stretching of the continental crust, coinciding with the Deception-Bridgeman volcanic line. Velocities of 8.1 km s⁻¹, characteristic of the Moho, are observed along the profile at a depth of 30-32 km.     

Lower-crustal rifting in the Rukwa Graben, East Africa

An M_w 5.9 earthquake occurred in the Lake Rukwa rift, Tanzania, on 1994 August 18, and was well recorded by 20 broad-band seismic stations at distances of 160 to 800 km and 21 broad-band stations at teleseismic distances. The regional and teleseismic waveforms have been used to investigate the source characteristics of the main shock, and also to locate aftershocks that occurred within three weeks of the main shock. Teleseismic body-wave modelling yields the following source parameters for the main shock: source depth of 25 ± 2 km, a normal fault orientation, with a horizontal tension axis striking NE-SW and an almost vertical pressure axis (Nodal Plane I: strike 126°–142°, dip 63°–66°, and rake 280°–290°; Nodal Plane II: strike 273°–289°, dip 28°–31°, and rake 235°–245°), a scalar moment of 4.1 times 10¹⁷ N m, and a 2 s impulsive source time function. Four of the largest aftershocks also nucleated at depths of 25 km, as deduced from regional sPmp–Pmp times. The nodal planes are broadly consistent with the orientation of both the Lupa and Ufipa faults, which bound the Rukwa rift to the northeast and southwest, respectively. The rupture radius of the main shock, assuming a circular fault, is estimated to be 4 km with a corresponding stress drop of 6.5 MPa. Published estimates of crustal thickness beneath the Rukwa rift indicate that the foci of the main shock and aftershocks lie at least 10 km above the Moho. The presence of lower-crustal seismicity beneath the Rukwa rift suggests that the pre-rift thermal structure of the rifted crust has not been strongly modified by the rifting, at least to depths of 25 km.     

A comparison of the isostatic response of bathymetric features in the north Pacific Ocean and Philippine Sea

Summary. This paper concerns the calculation and analysis of admittance functions from large and uniform data sets of gravity and topography in four regions of the northern and western Pacific Ocean. The purpose is to separate and describe possible differences in isostatic compensation between several 'type' regions of oceanic crust: a mid-ocean ridge (Juan de Fuca), a mid-plate seamount chain (Hawaiian Ridge), fracture zone topography on old crust (north of Hawaii) and a marginal basin (Philippine Sea). Results suggest that there are significant differences in the degree to which long wavelength topography has been compensated which can be distinguished between regions. These differences are set in the perspective of three simple compensation mechanisms. Two of these consider local Airy models in which raised topography is compensated at depth either by crustal roots or low density mantle. A third considers the effects of an elastic plate of variable thickness supporting crustal variations. Conclusions are that: (a) a thick plate possibly in excess of 30 km supports the Hawaiian Ridge; (b) a much thinner plate of 5 to 15 km existed when the fracture zone topography was formed; (c) the Juan de Fuca Ridge is compensated either regionally by a plate 5 to 10 km thick or locally by sub-crustal low densities at depths of 15 to 20 km; and (d) the Philippine Sea shows no evidence for regional support: ridges are compensated locally by differences in crustal thickness whereas the basins are underlain by density variations at depths comparable to those of the much younger Juan de Fuca Ridge. The major difference between admittance functions for the Philippine Sea and comparably aged regions of the north Pacific Ocean adds further new evidence of possible evolutionary differences between it and normal ocean basins.     

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.     

Crustal root beneath the highlands of southern Norway resolved by teleseismic receiver functions

Teleseismic data have been collected with temporary seismograph stations on two profiles in southern Norway. Including the permanent arrays NORSAR and Hagfors the profiles are 400 and 500 km long and extend from the Atlantic coast across regions of high topography and the Oslo Rift. A total of 1071 teleseismic waveforms recorded by 24 temporary and 8 permanent stations are analysed. The depth-migrated receiver functions show a well-resolved Moho for both profiles with Moho depths that are generally accurate within ±2 km.
For the northern profile across Jotunheimen we obtain Moho depths between 32 and 43 km (below sea level). On the southern profile across Hardangervidda, the Moho depths range from 29 km at the Atlantic coast to 41 km below the highland plateau. Generally the depth of Moho is close to or above 40 km beneath areas of high mean topography (>1 km), whereas in the Oslo Rift the crust locally thins down to 32 km. At the east end of the profiles we observe a deepening Moho beneath low topography. Beneath the highlands the obtained Moho depths are 4–5 km deeper than previous estimates. Our results are supported by the fact that west of the Oslo Rift a deep Moho correlates very well with low Bouguer gravity which also correlates well with high mean topography.
The presented results reveal a ca . 10–12 km thick Airy-type crustal root beneath the highlands of southern Norway, which leaves little room for additional buoyancy-effects below Moho. These observations do not seem consistent with the mechanisms of substantial buoyancy presently suggested to explain a significant Cenozoic uplift widely believed to be the cause of the high topography in present-day southern Norway.     

Curie-temperature depth estimation using a self-similar magnetization model

The Earth's crust is magnetized down to the Curie-temperature depth at about 10 to 50 km. This limited depth extent of the crustal magnetization is discernible in the power spectra of magnetic maps of South Africa and Central Asia. At short wavelengths, the power increases as rapidly towards longer wavelengths as expected for a self-similar magnetized crust with unlimited depth extent. Above wavelengths of about 100 km the power starts increasing less rapidly, indicating the absence of deep-seated sources. To quantify this effect we derive the theoretical power spectrum due to a slab carved out of a self-similar magnetization distribution. This model power spectrum matches the power spectra of South Africa and Central Asia for a self-similarity parameter of β = 4 and Curie temperature depths of 15 to 20 km.     

A geomagnetic study of the Great Glen Fault

Summary. An attempt is made to determine the range of two-dimensional current models consistent with the measured magnetovariational response, for periods from 5–30 min, near the Great Glen Fault in northern Scotland. All current models must be symmetric about the fault line but, because of uncertainty about the magnitude of the ocean effect, models ranging from a line current at 80 km depth to a uniform current sheet, 60 km wide, at 10 km depth are equally acceptable. Comparison with other geophysical studies of the same area suggests that a suitable conducting zone is unlikely to be present at shallow depths and interpretation in terms of a conducting zone in the 20–80 km depth range is favoured, although no such zone has been resolved by the other studies.     

On the anisotropy of sedimentary rocks from shear-wave analysis

Summary. The paper describes some results of experimental seismic prospecting investigations of the anisotropic properties of sedimentary rocks at depths of less than 2.5 km. Shear and converted PS -waves were recorded. Examples of SV and SH velocity distributions and shear-wave polarization are given. The main conclusion is that the medium essentially differs from that usually adopted for sedimentary rocks in that the symmetry axis is not normal to the layering.     

The crustal structure of the Reykjanes Ridge at 59° 30'N

Summary. In order to examine the development of the oceanic crust in the neighbourhood of a slowly spreading ridge, a seismic refraction experiment was carried out at 59° 30'N on the Reykjanes Ridge. Three 120 km long overlapped split profiles were shot parallel to the trend of the ridge, on the eastern flank, and recorded on up to five recording sonobuoys. The profiles were at distances of 0, 30 and 90km from the ridge axis, corresponding to approximate crustal ages of 0, 3 and 9 Myr. Data from the main profiles were supplemented by using a large chamber air gun during recovery of the buoys.
The analysis of the data combined standard travel-time interpretation, the 'tau' method of systematic travel-time inversion and detailed amplitude modelling using the Reflectivity Method to calculate synthetic seismograms. Detailed velocity-depth models were constructed for each of the profiles.
There is no indication of a significant magma chamber at the ridge crest, although a slight velocity inversion in layer 3 suggests a zone of elevated temperature. Away from the crest there was a slight positive velocity gradient in layer 3. Layer 2 was most effectively modelled by a region of varying velocity gradients, which thinned with age and the transition to layer 3 is marked by a sharp change in velocity gradient. The transition to mantle velocities is also best modelled by a high velocity gradient rather than an interface.
Although some lateral variation in properties is apparent along the profiles, the lateral velocity gradients were sufficiently weak to allow an effective analysis in terms of laterally uniform models.     

Upper mantle and lithospheric heterogeneities in central and eastern Europe as observed by teleseismic receiver functions

Data from 90 permanent broad-band stations spread over central and eastern Europe were analysed using Ps receiver functions to study the crustal and upper-mantle structure down to the mantle transition zone. Receiver functions provide valuable information on structural features, which are important for the resolution of European lithospheric dynamics. Moho depths vary from less than 25 km in extensional areas in central Europe to more than 50 km at stations in eastern Europe (Craton) and beneath the Alpine–Carpathian belt. A very shallow Moho depth can be observed at stations in the Upper Rhine Graben area ( ca. 25 km), whereas, for example, stations in the SW Bohemian Massif show a significantly deeper Moho interface at a depth of 38 km. V_p / V_s ratios vary between 1.60 and 1.96, and show no clear correlation to the major tectonic units, thus probably representing local variations in crustal composition. Delayed arrivals of converted phases from the mantle transition zone are observed at many stations in central Europe, whereas stations in the cratonic area show earlier arrivals compared with those calculated from the IASP91 Earth reference model. Differential delay times between the P₄₁₀s and P₆₆₀s phases indicate a thickened mantle transition zone beneath the eastern Alps, the Carpathians and the northern Balkan peninsula, whereas the transition zone thickness in eastern and central Europe agrees with the IASP91 value. The thickening of the mantle transition zone beneath the eastern Alps and the Carpathians could be caused by cold, deeply subducted oceanic slabs.