A three-dimensional image of shallow subduction: crustal structure of the Raukumara Peninsula, New Zealand

Earthquake arrival time data from a 36-station deployment of portable seismographs on the Raukumara Peninsula have been used to determine the 3-D Vp and Vp/Vs structure of this region of shallow subduction. A series of inversions have been performed, starting with an inversion for 1-D structure, then 2-D, and finally 3-D. This procedure ensures a smooth regional model in places of low resolution. The subducted plate is imaged as a northwest-dipping feature, with Vp consistently greater than 8.5  km  s⁻¹ in the uppermost mantle of the plate. Structure in the overlying plate changes significantly along strike. In the northeast, there is an extensive low-velocity zone in the lower crust underlying the most rapidly rising part of the Raukumara Range. It is bounded on its arcward side by an upwarp of high velocity. A viable explanation for the low-velocity zone is that it represents an accumulation of underplated subducted sediment, while serpentinization of the uppermost mantle may be responsible for the adjacent high-velocity region. The low-velocity zone decreases and the adjacent high-velocity region is less extensive in the southwest. This change is interpreted to be related to a change in the thickness of the crust of the overlying plate. In the northeast the crust is thinner, and subducted sediment ponds against relatively strong uppermost mantle, while in the southwest the crust is thicker, and the relatively weak lower crust allows sediment subduction to greater depths. A narrow zone of high Vp/Vs parallels the shallow part of the plate interface. This suggests elevated fluid pressures, with the distribution of earthquakes about this zone further suggesting that these pressures may be close to lithostatic. The plate interface at 20  km depth beneath the Raukumara Peninsula may thus be a closed system for fluid flow, similar to that seen at much shallower depths in other subduction décollements.     

Crustal and uppermost mantle structure in southern Africa revealed from ambient noise and teleseismic tomography

Rayleigh wave phase velocity maps in southern Africa are obtained at periods from 6 to 40 s using seismic ambient noise tomography applied to data from the Southern Africa Seismic Experiment (SASE) deployed between 1997 and 1999. These phase velocity maps are combined with those from 45 to 143 s period which were determined previously using a two-plane-wave method by Li & Burke. In the period range of overlap (25–40 s), the ambient noise and two-plane-wave methods yield similar phase velocity maps. Dispersion curves from 6 to 143 s period were used to estimate the 3-D shear wave structure of the crust and uppermost mantle on an 1°× 1° grid beneath southern Africa to a depth of about 100 km. Average shear wave velocity in the crust is found to vary from 3.6 km s^–1 at 0–10 km depths to 3.86 km s^–1 from 20 to 40 km, and velocity anomalies in these layers correlate with known tectonic features. Shear wave velocity in the lower crust is on average low in the Kaapvaal and Zimbabwe cratons and higher in the surrounding Proterozoic terranes, such as the Limpopo and the Namaqua-Natal belts, which suggests that the lower crust underlying the Archean cratons is probably less mafic than beneath the Proterozoic terranes. Crustal thickness estimates agree well with a previous receiver function study of Nair et al. . Archean crust is relatively thin and light and underlain by a fast uppermost mantle, whereas the Proterozoic crust is thick and dense with a slower underlying mantle. These observations are consistent with the southern African Archean cratons having been formed by the accretion of island arcs with the convective removal of the dense lower crust, if the foundering process became less vigorous in arc environments during the Proterozoic.     

Structure and early evolution of the Arabian Sea and East Somali Basin

The Laxmi Ridge is a large-scale basement high buried beneath the sediments of the Indus Fan. The location of the ocean–continent transition (OCT) on this margin has previously been proposed at either the southern edge of the Laxmi Ridge or beyond it towards the India–Pakistan shelf. The former explains the margin-parallel Laxmi Basin as thinned continental crust, the latter as a failed rift of earlier seafloor spreading. To examine the structure of this margin, a reassessment of marine magnetic data has detailed seafloor-spreading magnetic anomalies prior to anomaly 24 in both the Arabian and East Somali basins. The previously identified anomaly 28 is not interpreted as a seafloor-spreading anomaly but as a magnetized basement feature adjacent to, and merging with, the ridge—the Laxmi Spur. New gravity models across the Laxmi Ridge and adjacent margin using ship and satellite data corroborate the existence of underplated crust beneath the Laxmi Ridge and Basin and the location of the OCT at the southern edge of the Ridge. The results are not compatible with the existence of a pre-anomaly 28 phase of seafloor spreading, although large-scale intrusions may be the origin of some of the basement features in the Laxmi Basin. The models also identify the Laxmi Spur as a low-density feature with a natural remanent magnetization (NRM) compatible with serpentinization. The Laxmi Ridge is mapped to the southeast, where it appears to terminate at a point coinciding with the appearance of E–W magnetic lineations and gravity anomalies at 15.5°N. Thereafter it becomes indistinct. This is interpreted as necessary in the reconstruction to the Mascarene Plateau to avoid continental overlap.     

LISPB – IV. Crustal structure of Northern Britain

Summary. This paper presents those results from the 1974 Lithospheric Seismic Profile in Britain (LISPB) which relate to the compressional velocity structure of the crust and uppermost mantle beneath Northern Britain. A combination of interpretation techniques suitable for modelling laterally inhomogeneous media, including two-dimensional ray-tracing and time-term analysis, has resulted in a detailed seismic cross-section across the Caledonian orogenic belt. The main features of this section are a possible horizontal discontinuity in the Pre-Caledonian'basement, a change in the relationship between the lower crust and the uppermost mantle from north to south and a considerable thickening of the crust beneath the Caledonian fold belt. These results place considerable constraints upon tectonic models for the evolution of the Caledonides in particular in their implication of differing crustal structures north and south of the Southern Uplands and their indication of the primary significance of the Southern Uplands Fault.     

Determination of the lower edges of magnetized bodies by using geothermal data

The determination of the lower edges of magnetized bodies in the Earth's crust is a complex geophysical problem, although these values can be estimated by using geothermal data. An analysis of the temperature regime and location of the lower edges of magnetized bodies has been carried out for the geosynclinal region of the southern Caucasus and the area joining the ancient platform with the Arabian Shield in Israel. Geothermal calculations for Israel have been performed for three models of the thermal regime for the Earth's crust and upper mantle. The process of ultrabasic rock serpentinization is accompanied by the transformation of iron suboxide to iron oxide. Both these processes run under identical thermodynamic conditions within an average temperature interval of 200°-400°C. The Curie surface controls the position of lower edges only in fault zones where oxidation conditions hold up to great depths.     

Upper mantle heterogeneities beneath Central Europe

Summary. An inversion of ISC travel-time data from selected earthquakes in the distance range 30°-90° to 53 stations in Central Europe has been used to model velocity down to 600 km depth. The model explains 0.1–0.2s of the residuals, as for other array studies, leaving 0.5 s unexplained as noise. The uppermost 100 km of the mantle and crust contains inhomogeneities that correlate remarkably well with the geology. This may be due to deep-seated thermal anomalies or, in some areas, to delays introduced by passage of the rays through sedimentary cover. The deeper anomalies are smaller and unrelated to those in the lithosphere, which suggests that the asthenosphere is decoupled from the rigid lithosphere. The structure at 600 km depth is again quite inhomogeneous and might be due to undulations of the 650 km discontinuity. The models show some suggestion of a high velocity slab trending from east to west beneath the Alps.     

Effect of the presence of a conducting channel between India and Sri Lanka Island on the features of the equatorial electrojet

Summary. Transient geomagnetic variations like SSCs, Bays, Sq and storm-time variations show anomalously large Z amplitudes at the three permanent magnetic observatories in India under the equatorial electrojet. Our earlier studies have shown that these anomalies cannot be explained in terms of the usual coastal effect. Another unique feature of this area is the small equatorial enhancement of all magnetic fluctuations, which still remains to be completely understood. It is also noticed that Z/H ratio at Annamalainagar is very large for night-time variations but becomes insignificant for day-time fluctuations, whereas no such large difference is seen at Trivandrum, the other coastal station. All the above features have been explained here by introducing the presence of a conducting channel in the lower crust or upper mantle between India and Sri Lanka Island. The anomalies of the equatorial electrojet in the Indian region thus do not necessitate the conducting surface of the mantle to be deeper in this area as has been suggested by earlier workers. A less conducting mantle in the Indian region otherwise could be difficult to incorporate in the present theories of mantle convection when this area is known to be tectonically active. The presence of the conducting channel is further confirmed by analogue model experiments of Papamastorakis. We further observe that the equatorial electrojet does have an associated internal part in the Indian region.     

Improved magnetic anomalies of the Antarctic lithosphere from satellite and near-surface data

The Antarctic magnetic anomaly map compiled marine and airborne surveys collected south of 60°S through 1999 and used Magsat data to help fill in the regional gaps between the surveys. Ørsted and CHAMP satellite magnetic observations with greatly improved measurement accuracies and temporal and spatial coverage of the Antarctic, have now supplanted the Magsat data. We combined the new satellite observations with the near-surface survey data for an improved magnetic anomaly map of the Antarctic lithosphere. Specifically, we separated the crustal from the core and external field components in the satellite data using crustal thickness variations estimated from the terrain and the satellite-derived free-air gravity observations. Regional gaps in the near-surface surveys were then filled with predictions from crustal magnetization models that jointly satisfied the near-surface and satellite crustal anomalies. Comparisons in some of the regional gaps that also considered newly acquired aeromagnetic data demonstrated the enhanced anomaly estimation capabilities of the predictions over those from conventional minimum curvature and spherical harmonic geomagnetic field models. We also noted that the growing number of regional and world magnetic survey compilations involve coverage gaps where these procedures can contribute effective near-surface crustal anomaly estimates.     

Seismic tomographic imaging of P- and S-waves velocity perturbations in the upper mantle beneath Iran

The inverse tomography method has been used to study the P - and S -waves velocity structure of the crust and upper mantle underneath Iran. The method, based on the principle of source–receiver reciprocity, allows for tomographic studies of regions with sparse distribution of seismic stations if the region has sufficient seismicity. The arrival times of body waves from earthquakes in the study area as reported in the ISC catalogue (1964–1996) at all available epicentral distances are used for calculation of residual arrival times. Prior to inversion we have relocated hypocentres based on a 1-D spherical earth's model taking into account variable crustal thickness and surface topography. During the inversion seismic sources are further relocated simultaneously with the calculation of velocity perturbations. With a series of synthetic tests we demonstrate the power of the algorithm and the data to reconstruct introduced anomalies using the ray paths of the real data set and taking into account the measurement errors and outliers. The velocity anomalies show that the crust and upper mantle beneath the Iranian Plateau comprises a low velocity domain between the Arabian Plate and the Caspian Block. This is in agreement with global tomographic models, and also tectonic models, in which active Iranian plateau is trapped between the stable Turan plate in the north and the Arabian shield in the south. Our results show clear evidence of the mainly aseismic subduction of the oceanic crust of the Oman Sea underneath the Iranian Plateau. However, along the Zagros suture zone, the subduction pattern is more complex than at Makran where the collision of the two plates is highly seismic.     

Deep seismic reflection studies in Israel - an update

Summary. The results of three deep crustal reflection lines are presently available from Israel. A 90 km line from near the Dead Sea rift to the Mediterranean coast was carried out for deep study. Two other lines in the Mediterranean coastal area were derived by recorrelation of oil exploration lines. The data shows a division between continental inner Israel and the coastal plain. In the first area a reflective lower crust is apparent with transparent upper crust and almost transparent upper mantle. Near the coast, in an area which was previously suggested as underlain by an ancient fossil oceanic crust, strong reflections characterize the uppermost mantle. Comparison between the reflection pattern and previous deep refraction and MT data indicates some agreement away from the coast and lack of correlation in the area of possible fossil oceanic crust near the coast.     

On the roughness of Mesozoic oceanic crust in the western North Atlantic

Seismic reflection profiles from Mesozoic oceanic crust around the Blake Spur Fracture Zone (BSFZ) in the western North Atlantic have been widely used in constraining tectonic models of slow-spreading mid-ocean ridges. These profiles have anomalously low basement relief compared to crust formed more recently at the Mid-Atlantic Ridge at the same spreading rate. Profiles from other regions of Mesozoic oceanic crust also have greater relief. The anomalous basement relief and slightly increased crustal thickness in the BSFZ survey area may be due to the presence of a mantle thermal anomaly close to the ridge axis at the time of crustal formation. If so, the intracrustal structures observed may be representative of an atypical tectonic regime.     

Travel-time anomalies in the mantle under the North Atlantic

Summary. Lateral heterogeneity exists in the Earth's mantle, and may result in seismic velocity anomalies up to several per cent. If convection cells and plumes extend down to the core, then these features may be associated with local inhomogeneities observed in the lower mantle.
Published data for direct and core-reflected P -wave residuals are used to delineate velocity anomalies in the middle—lower mantle under the North Atlantic. Differential ( PcP — P ) residuals indicate travel-time anomalies near the core—mantle transition, and may be due to core topography or lateral variations in velocity. It is assumed that the anomalies occur near the midpoints of the ray paths. The main source of error in the data set may arise from phases which have been identified incorrectly. Hence trend-surfaces are fitted to the residual data to show only the large-scale trends in anomaly values, with wavelengths of the order of 1000 km.
The Azores and Colorado hot spots occur in a region covered by the data. A possible interpretation of the trend maps is that an anomalous zone extends from a relatively fast region at the core boundary at 35° N, 50° W up to these hot spots, at about 30 degrees from the vertical. This may agree with the suggestion of Anderson that plumes are chemical rather than thermal in origin. If inclⁱned plumes do exist, the deviation from the ideal vertical plume or convection cell boundary may imply that lateral shear or other distortion effects exist in the mantle.     

Interpretation of aeromagnetic anomalies over the Lower and Middle Benue Trough of Nigeria

Summary. The results of an interpretation of aeromagnetic anomalies over the Lower and Middle Benue Trough of Nigeria are presented. Two-dimensional interpretation of the anomalies shows that the anomalies are best interpreted in terms of basic intrusive bodies occurring either within the Cretaceous sediments or within the metamorphic basement or both. The model suggests that the intrusive bodies are of variable thickness and magnetization polarity, and may have been emplaced during different epochs of polarity. The model is in good agreement with interpretations of gravity and magnetic anomalies over similar structures elsewhere in the world.     

The seismic structure of thinned continental crust in the northern Bay of Biscay

Summary. The stretching and thinning of the continental crust, which occurs during the formation of passive continental margins, may cause important changes in the velocity structure of such crust. Further, crust attenuated to a few kilometres' thickness, can be found underlying 'oceanic' water depths. This paper poses the question of whether thinned continental crust can be distinguished seismically from normal oceanic crust of about the same thickness. A single seismic refraction line shot over thinned continental crust as part of the North Biscay margin transect in 1979 was studied in detail. Tau— p inversion suggested that there are differences between oceanic and continental crust in the lower crustal structure. This was confirmed when synthetic seismograms were calculated. The thinned continental crust (β± 7.0) exhibits a two-gradient structure in the non-sedimentary crust with velocities between 5.9 and 7.4 km s⁻¹; an upper 0.8 s⁻¹ layer overlies a 0.4 s⁻¹ layer. No layer comparable to oceanic layer 3 was detected. The uppermost mantle also contains a low-velocity zone.     

Interpretation of aeromagnetic data across the central crystalline shield area of Nigeria

Summary. A residual map of the total magnetic field (above 25 000 nT base) is presented for a portion of the central crystalline shield area of Nigeria and overlapping small portions of the Chad basin and the Benue rift (8°30'−12° 00'lat, and 7°−10°30' long). The map (based on a dataset digitized from recently released aeromagnetic sheets of Nigeria) leads to four results. (1) A magnetic boundary, evident on the map, separates the Younger Granite complexes into two groups. The groups are petrologically different, and the boundary may be a fault line with uplift to the south. (2) South of the boundary the map is dominated by a system of sub-parallel anomalies striking NE–SW, possibly representing major tectonic trends, and a set of fractures through which the Younger Granite complexes were intruded. The trend of the system parallels the Benue rift and lineaments in the oceanic crust off West Africa. (3) Negative magnetic anomalies lie over most of the known ring complexes, and over some suspected buried ring complexes and other intrusions. (4) 2½-and 3-D modelling shows that the larger complexes extend to 12 km depth, and the smaller ones to 6 km. They have nearly vertical sides, and magnetization contrasts range from 0.3 to 0.5 A m⁻¹.     

Crustal structure of Atlantic fracture zones–II. The Vema fracture zone and transverse ridge

Summary. The crustal structure beneath the Vema fracture zone and its flanking transverse ridge was determined from seismic refraction profiles along the fracture zone valley and across the ridge. Relatively normal oceanic crust, but with an upwarped seismic Moho, was found under the transverse ridge. We suggest that the transverse ridge represents a portion of tectonically uplifted crust without a major root or zone of serpentinite diapirism beneath it. A region of anomalous crust associated with the fracture zone itself extends about 20 km to either side of the central fault, gradually decreasing in thickness as the fracture zone is approached. There is evidence to suggest that the thinnest crust is found beneath the edges of the 20 km wide fracture zone valley. Under the fracture zone valley the crust is generally thinner than normal oceanic crust and is also highly anomalous in its velocity structure. Seismic layer 3 is absent, and the seismic velocities are lower than normal. The absence of layer 3 indicates that normal magmatic accretionary processes are considerably modified in the vicinity of the transform fault. The low velocities are probably caused by the accumulation of rubble and talus and by the extensive faulting and fracturing associated with the transform fault. This same fracturing allows water to penetrate through the crust, and the apparently somewhat thicker crust beneath the central part of the fracture zone valley may be explained by the resultant serpentinization having depressed the seismic Moho below its original depth.     

Heterogeneous structure around the Jemez volcanic field, New Mexico, USA, as inferred from the envelope inversion of active-experiment seismic data

We analyse active-experiment seismic data obtained by the 1993 Jemez Tomography Experiment (JTEX) programme to elucidate the heterogeneous structure of the Jemez volcanic field, which is located at the boundary between the Colorado Plateau and the Rio Grande Rift. Using a single isotropic scattering assumption, we first calculate the envelope Green's functions for the upper and lower crust and the uppermost mantle. By fitting the theoretical envelopes with the observed three-component data, we estimate depth-dependent features of the scattering coefficients around Valles Caldera. We estimate the ratios of scattering coefficients, rather than scattering coefficients themselves, because of the uncertainty of the seismic efficiency of the explosive sources and knowledge of absolute site-amplification factors. The strongest scattering coefficients are observed at a shallow depth beneath the Valles Caldera. This is considered to be related to the complex structure caused by two episodes of caldera formation and the ensuing resurgent uplift in the caldera, etc. The depth-dependent characters of the scattering coefficients for the Colorado Plateau and the Rio Grande Rift are similar to each other: a transparent upper crust and a heterogeneous lower crust (small and large scattering coefficients, respectively). However, the scattering coefficients beneath the Rio Grande Rift are several times larger than those beneath the Colorado Plateau. Depths of the lower crust and the Moho boundary beneath the Rio Grande Rift are shallower than those of the Colorado Plateau. From their geological settings and other geophysical results around the region, we infer that the larger scattering coefficients of the rift are associated with rift formation and volcanic activity, such as magma ascent from the upper mantle to the crust.     

Regional three-dimensional gravity modelling of the NE Atlantic margin

A series of three‐dimensional models has been constructed for the structure of the crust and upper mantle over a large region spanning the NE Atlantic passive margin. These incorporate isostatic and flexural principles, together with gravity modelling and integration with seismic interpretations. An initial isostatic model was based on known bathymetric/topographic variations, an estimate of the thickness and density of the sedimentary cover, and upper mantle densities based on thermal modelling. The thickness of the crystalline crust in this model was adjusted to equalise the load at a compensation depth lying below the zone of lateral mantle density variations. Flexural backstripping was used to derive alternative models which tested the effect of varying the strength of the lithosphere during sediment loading. The models were analysed by comparing calculated and observed gravity fields and by calibrating the predicted geometries against independent (primarily seismic) evidence. Further models were generated in which the thickness of the sedimentary layer and the crystalline crust were modified in order to improve the fit to observed gravity anomalies. The potential effects of igneous underplating and variable upper mantle depletion were explored by a series of sensitivity trials. The results provide a new regional lithospheric framework for the margin and a means of setting more detailed, local investigations in their regional context. The flexural modelling suggests lateral variations in the strength of the lithosphere, with much of the margin being relatively weak but areas such as the Porcupine Basin and parts of the Rockall Basin having greater strength. Observed differences between the model Moho and seismic Moho along the continental margin can be interpreted in terms of underplating. A Moho discrepancy to the northwest of Scotland is ascribed to uplift caused by a region of upper mantle with anomalously low density, which may be associated with depletion or with a temperature anomaly.     

Crustal Structure and Origin of Peake and Freen Deeps, NE Atlantic

Summary Peake and Freen Deeps are elongate structures some 30 nautical miles long by 7 miles wide situated near 43° N 20° W on the lower flanks of the Mid-Atlantic Ridge. Seismic reflection records show that underneath about 400 fm of layered sediment the bedrock lies at a depth greater than 3600 fm in Peake Deep and 3300 fm in Freen Deep; the surrounding seafloor is at about 2100 fm. Freen Deep is the eastern end of King's Trough, a flat floored feature some 400 fms deeper than the adjacent seafloor. The Trough extends 220 miles west-north-westwards towards the crest of the Mid-Atlantic Ridge. The area is aseismic and heat flow is normal; there is no displacement of the crest of the mid-ocean ridge on the projected line of King's Trough. Gravity and magnetic surveys have been made. With minor exceptions, magnetic anomalies are not due to bodies elongated parallel with the structure, which, therefore, cannot be a volcanic collapse caldera. Seismic refraction results in the Peake-Freen area show that the crust is not thinned under the deeps although the Moho may be depressed by 2 km. Bouguer anomalies also suggest that the Moho is flat and does not rise to compensate the deeps. Models consistent with gravity and seismic information suggest there is a dense block in the upper mantle under the area. Since no reason to ascribe the origin of the structure to tear faulting has yet been acquired, it is interpreted in terms of over thrusting perpendicular to the deeps, followed by inversion of the lower part of the thickened basaltic crust to eclogite, and its subsequent sinking into the mantle.     

Numerical models of ductile rebound of crustal roots beneath mountain belts

Crustal roots formed beneath mountain belts are gravitationally unstable structures, which rebound when the lateral forces that created them cease or decrease significantly relative to gravity. Crustal roots do not rebound as a rigid body, but undergo intensive internal deformation during their rebound and cause intensive deformation within the ductile lower crust. 2-D numerical models are used to investigate the style and intensity of this deformation and the role that the viscosities of the upper crust and mantle lithosphere play in the process of root rebound. Numerical models of root rebound show three main features which may be of general application: first, with a low-viscosity lower crust, the rheology of the mantle lithosphere governs the rate of root rebound; second, the amount of dynamic uplift caused by root rebound depends strongly on the rheologies of both the upper crust and mantle lithosphere; and third, redistribution of the rebounding root mass causes pure and simple shear within the lower crust and produces subhorizontal planar fabrics which may give the lower crust its reflective character on many seismic images.