Travel times of Pn -waves in the Aegean and surrounding area

Summary. Based on accurately located 23 very shallow earthquakes ( h = 1–14 km) in northern and central Greece by portable networks of seismic stations and by the joint epicentre method, the travel times of the P_n -waves from the foci of these earthquakes to the sites of 54 permanent stations in the Balkan region have been determined. The travel times of P_n -waves in the central and eastern part of the area (eastern Greece, south-eastern Yugoslavia, the Aegean Sea, Bulgaria, southern Romania, western Turkey) fit a straight line very well with the P_n velocity equal to 7.9 ± 0.1 km s^-1. On the contrary, the travel times of P_n -waves to stations in the western part of the area (Albania, western Greece) do not fit this curve because the P_n -waves travelling to these stations are delayed by more than 1 s due to the thicker crust under the Dinarides–Hellenides mountain range. Time delays for P_n -waves have been calculated for each permanent station in the Balkan area with respect to the mean travel-time curve of these waves in the central and eastern part of the area. Corrections of the travel times for these delays contribute very much to the improvement of the accuracy in the location of the shallow earthquakes in the Aegean and surrounding area.     

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.     

Flexural uplift of the Stara Planina range, central Bulgaria

The Stara Planina is an E–W-trending range within the Balkan belt in central Bulgaria. This topographically high mountain range was the site of Mesozoic through early Cenozoic thrusting and convergence, and its high topography is generally thought to have resulted from crustal shortening associated with those events. However, uplift of this belt appears to be much younger than the age of thrusting and correlates instead with the age of Pliocene–Quaternary normal faulting along the southern side of the range. Flexural modelling indicates the morphology of the range is consistent with flexural uplift of footwall rocks during Pliocene–Quaternary displacement on S-dipping normal faults bounding the south side of the mountains, provided that the effective elastic plate thickness of 12  km under the Moesian platform is reduced to about 3  km under the Stara Planina. This small value of elastic plate thickness under the Stara Planina is similar to values observed in the Basin and Range Province of the western United States, and suggests that weakening of the lithosphere is due to heating of the lithosphere during extension, perhaps to the point that large-scale flow of material is possible within the lower crust. Because weakening is observed to affect the Moesian lithosphere for ≈10  km beyond (north of) the surface expression of extension, this study suggests that processes within the uppermost mantle, such as convection, play an active role in the extension process. The results of this study also suggest that much of the topographic relief in thrust belts where convergence is accompanied by coeval extension in the upper plate (or 'back arc'), such as in the Apennines, may be a flexural response to unloading during normal faulting, rather than a direct response to crustal shortening in the thrust belt.     

Evidence for underthrusting beneath the Queen Charlotte Margin, British Columbia, from teleseismic receiver function analysis

The Queen Charlotte Fault zone is the transpressive boundary between the North America and Pacific Plates along the northwestern margin of British Columbia. Two models have been suggested for the accommodation of the ∼20 mm yr⁻¹ of convergence along the fault boundary: (1) underthrusting; (2) internal crustal deformation. Strong evidence supporting an underthrusting model is provided by a detailed teleseismic receiver function analysis that defines the underthrusting slab. Forward and inverse modelling techniques were applied to receiver function data calculated at two permanent and four temporary seismic stations within the Queen Charlotte Islands. The modelling reveals a ∼10 km thick low-velocity zone dipping eastward at 28° interpreted to be underthrusting oceanic crust. The oceanic crust is located beneath a thin (28 km) eastward thickening (10°) continental crust.     

Basement faulting and the focal depths of the larger earthquakes in the Zagros mountains (Iran)

Summary. The present day seismicity of the Zagros seems to occur on high angle reverse faults distributed across the whole width of the belt. It does not indicate activity on a single inclined thrust surface and there do not seem to have been any well located events at intermediate depths. Modelling of the long period teleseismic body waves of seven large earthquakes presented here shows their focal depths to be in the range 8–15 km. This is thought to indicate faulting in the uppermost basement beneath the sedimentary cover, though the absence of published seismic refraction results renders the sediment thickness uncertain. Faulting of this type and depth may occur on inherited normal faults which have subsequently been reactivated as thrusts. Such reactivation allows considerable shortening to take place in the early stages of continental collision without the subduction or excessive thickening of continental crust.     

Lithospheric structure of the Arabian and Eurasian collision zone in eastern Turkey from S-wave receiver functions

Crustal and upper-mantle seismic discontinuities beneath eastern Turkey are imaged using teleseismic S -to- P converted phases. Three crustal phases are observed: the Moho with depth ranging between 30 and 55 km, indicating variable tectonic regimes within this continental collision zone; an upper-crustal discontinuity at approximately 10 km depth; and various crustal low-velocity zones, possibly associated with recent Quaternary volcanism. Imaging of the upper mantle is complicated by the 3-D geometry of the region, in particular due to the Bitlis–Zagros suture zone. However, several upper-mantle S -to- P converted phase are identified as being the signature of the lithosphere–asthenosphere boundary (LAB). The inferred LAB for the Eastern Anatolian Accretionary Complex indicates that eastern Turkey has an anomalously thin (between ∼60 and 80 km) lithosphere which is consistent with an oceanic slab detachment model. The observed LAB phases for the Arabian shield and Iranian plateau indicate that lithospheric thickness for these stable regions is on the order of 100 to 125 km thick, which is typical of continental margins.     

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.     

Basement depth and sedimentary velocity structure in the northern Arabian platform, eastern Syria

Basement depth in the Arabian plate beneath eastern Syria is found to be much deeper than previously supposed. Deep-seated faulting in the Euphrates fault system is also documented. Data from a detailed 300 km long reversed refraction profile, with offsets up to 54 km, are analysed and interpreted, yielding a velocity model for the upper 9 km of continental crust. The interpretation integrates the refraction data with seismic-reflection profiles, well logs and potential field data, such that the results are consistent with all available information. A model of sedimentary thicknesses and seismic velocities throughout the region is established. Basement depth on the north side of the Euphrates is interpreted to be around 6 km, whilst south of the Euphrates basement depth is at least 8.5 km. Consequently, the potentially hydrocarbon-rich pre-Mesozoic section is shown, in places, to be at least 7 km thick. The dramatic difference in basement depth on adjacent sides of the Euphrates graben system may suggest that the Euphrates system is a suture/shear zone, possibly inherited from Late Proterozoic accretion of the Arabian plate. Gravity modelling across the southeast Euphrates system tends to support this hypothesis. Incorporation of previous results allows us to establish the first-order trends in basement depth throughout Syria     

New insight into the crust and upper mantle structure under Alaska

To better understand the seismic structure of the subducting Pacific plate under Alaska, we determined the three-dimensional P-wave velocity structure to a depth of approximately 200 km beneath Alaska using 438,146 P-wave arrival times from 10,900 earthquakes. In this study an irregular grid parameterization was adopted to express the velocity structure under Alaska. The number of grid nodes increases from north to south in the study area so that the spacing between grid nodes is approximately the same in the longitude direction. Our results suggest that the subducting Pacific slab under Alaska can be divided into three different parts based on its geometry and velocity structure. The western part has features similar to those in other subduction zones. In the central part a thick low-velocity zone is imaged at the top of the subducting Pacific slab beneath north of the Kenai Peninsula, which is believed to be most likely the oceanic crust plus an overlying serpentinized zone and the coupled Yakutat terrane subducted with the Pacific slab. In the eastern part, significant high-velocity anomalies are visible to 60–90 km depth, suggesting that the Pacific slab has only subducted down to that depth.     

Crustal structure beneath exposed accreted terranes of Southern Alaska

Summary. The crustal structure beneath the exposed terranes of southern Alaska has been explored using coincident seismic refraction and reflection profiling. A wide-angle reflector at 8–9 km depth, at the base of an inferred low-velocity zone, underlies the Peninsular and Chugach terranes, appears to truncate their boundary, and may represent a horizontal decollement beneath the terranes. The crust beneath the Chugach terrane is characterized by a series of north-dipping paired layers having low and high velocities that may represent subducted slices of oceanic crust and mantle. This layered series may continue northward under the Peninsular terrane. Earthquake locations in the Wrangell Benioff zone indicate that at least the upper two low-high velocity layer pairs are tectonically inactive and that they appear to have been accreted to the base of the continental crust. The refraction data suggest that the Contact fault between two similar terranes, the Chugach and Prince William terranes, is a deeply penetrating feature that separates lower crust (deeper than 10 km) with paired dipping reflectors, from crust without such reflectors.     

Crustal thickness estimation beneath the southern central Andes at 30°S and 36°S from S wave receiver function analysis

We have used the S wave receiver function (SRF) technique to investigate the crustal thickness beneath two seismic profiles from the CHARGE project in the southern central Andes. A previous study employing the P wave receiver function method has observed the Moho interface beneath much of the profiles. They found, however, that the amplitude of the P to S conversion was diminished in the western part of the profiles and have attributed it to a reduction of the impedance contrast at the Moho due to lower crustal ecologitization. With SRF, we have successfully detected S to P converted waves from the Moho as well as possible conversions from other lithospheric boundaries. The continental South American crust reaches its maximum thickness of ∼70 km (along 30°S between 70°W and 68.5°W) beneath the Principal Cordillera and the Famatina system and becomes thinner towards the Sierras Pampeanas with a thickness of ∼40 km. Negative phases, possibly related to the base of the continental and oceanic lithosphere, can be recognized in the summation traces at different depths. By comparing our results with data obtained from previous investigations, we are able to further constrain the thickness of the crust and lithosphere beneath the central Andes.     

A seismological study of normal faulting in the Demirci, Alaşehir and Gediz earthquakes of 1969–70 in western Turkey: implications for the nature and geometry of deformation in the continental crust

Summary. In this study, seismological techniques are combined with surface observations to investigate the faulting associated with three large earthquakes in western Turkey. All involved normal faulting that nucleated at 6–10 km depth with dips in the range 30–50°. The two largest earthquakes, at Alaşehir (1969.3.28) and Gediz (1970.3.28), were clearly multiple events and their seismograms indicate that at least two discrete subevents were involved in producing the observed surface faulting. In addition, their seismograms contain later, longer-period signals that are likely to represent source, not structure or propagation, complexities. These later signals can be modelled by subevents with long time functions on almost flat detachment-type faults.
As a result of these observations, we propose a model for the deformation of the lower crust, in which brittle failure of the top part occurs when high strain rates are imposed during an earthquake that ruptures right through the upper, brittle crust. Under these special circumstances, seismic motion occurs on discrete faults in the lower crust, which otherwise normally deforms by distributed creep. In the case of the normal faults studied here, motion in the uppermost lower crust takes place on shallow dipping faults that are downward continuations of the steeper faults that break to the surface. The faults thus have an overall listric geometry, flattening into a weak zone below the brittle layer at a depth that is probably dependent on the termperature gradient. This interpretation explains why detachment-type mechanisms are not seen in first motion fault plane solutions of normal faulting earthquakes, and suggests an origin for the Metamorphic Core Complexes seen in the Basin and Range Province, which probably represent flat lower crustal faults, analogous to those postulated at Alaşehir and Gediz, that have been uplifted to the surface.     

Earthquake processes of the Himalayan collision zone in eastern Nepal and the southern Tibetan Plateau

Focal mechanisms determined from moment tensor inversion and first motion polarities of the Himalayan Nepal Tibet Seismic Experiment (HIMNT) coupled with previously published solutions show the Himalayan continental collision zone near eastern Nepal is deforming by a variety of styles of deformation. These styles include strike-slip, thrust and normal faulting in the upper and lower crust, but mostly strike-slip faulting near or below the crust–mantle boundary (Moho). One normal faulting earthquake from this experiment accommodates east–west extension beneath the Main Himalayan Thrust of the Lesser Himalaya while three upper crustal normal events on the southern Tibetan Plateau are consistent with east–west extension of the Tibetan crust. Strike-slip earthquakes near the Himalayan Moho at depths >60 km also absorb this continental collision. Shallow plunging P -axes and shallow plunging EW trending T -axes, proxies for the predominant strain orientations, show active shearing at focal depths ∼60–90 km beneath the High Himalaya and southern Tibetan Plateau. Beneath the southern Tibetan Plateau the plunge of the P -axes shift from vertical in the upper crust to mostly horizontal near the crust–mantle boundary, indicating that body forces may play larger role at shallower depths than at deeper depths where plate boundary forces may dominate.     

Seismic refraction profiles between Cyprus and Israel and their interpretation

Summary. A long seismic refraction profile was carried out between southern Israel and Cyprus. The seismic energy was generated by 33 sea shots each of 0.8 t explosives and was recorded by land stations in Israel and Cyprus and by ocean bottom seismographs deployed along the profile.
The results showed that the continental crust of southern Israel thins towards the Mediterranean underneath a northward thickening sedimentary cover. Cyprus is underlain by a 35 km thick continental crust thinning south-wards and extending to Mt Eratosthenes. Between Mt Eratosthenes and the Israel continental shelf the crystalline crust is composed of high velocity (6.5 km s^-1)material and is about 8 km thick. It is covered by 12–14 km of sediments and may represent a fossil oceanic crust.     

The continent‐ocean transition on the northwestern South China Sea

Rifted margins are created as a result of stretching and breakup of continental lithosphere that eventually leads to oceanic spreading and formation of a new oceanic basin. A cornerstone for understanding what processes control the final transition to seafloor spreading is the nature of the continent‐ocean transition (COT). We reprocessed multichannel seismic profiles and use available gravity data to study the structure and variability of the COT along the Northwest subbasin (NWSB) of the South China Sea. We have interpreted the seismic images to discern continental from oceanic domains. The continental‐crust domain is characterized by tilted fault blocks generally overlain by thick syn‐rift sedimentary units, and underlain by fairly continuous Moho reflections typically at 8–10 s twtt. The thickness of the continental crust changes greatly across the basin, from ~20 to 25 km under the shelf and uppermost slope, to ~9–6 km under the lower slope. The oceanic‐crust domain is characterized by a highly reflective top of basement, little faulting, no syntectonic strata and fairly constant thickness (over tens to hundreds of km) of typically 6 km, but ranging from 4 to 8 km. The COT is imaged as a ~5–10 km wide zone where oceanic‐type features directly abut or lap on continental‐type structures. The South China margin continental crust is cut by abundant normal faults. Seismic profiles show an along‐strike variation in the tectonic structure of the continental margin. The NE‐most lines display ~20–40 km wide segments of intense faulting under the slope and associated continental‐crust thinning, giving way to a narrow COT and oceanic crust. Towards the SW, faulting and thinning of the continental crust occurs across a ~100–110 km wide segment with a narrow COT and abutting oceanic crust. We interpret this 3D structural variability and the narrow COT as a consequence of the abrupt termination of continental rifting tectonics by the NE to SW propagation of a spreading centre. We suggest that breakup occurred abruptly by spreading centre propagation rather than by thinning during continental rifting. We propose a kinematic evolution for the oceanic domain of the NWSB consisting of a southward spreading centre propagation followed by a first narrow ridge jump to the north, and then a younger larger jump to the SE, to abandon the NWSB and create the East subbasin of the South China Sea.     

Seismic evidence for a sharp lithospheric base persisting to the lowermost mantle beneath the Caribbean

Broad-band data from South American earthquakes recorded by Californian seismic networks are analysed using a newly developed seismic wave migration method—the slowness backazimuth weighted migration (SBWM). Using the SBWM, out-of-plane seismic P -wave reflections have been observed. The reflection locations extend throughout the Earth's lower mantle, down to the core–mantle boundary (CMB) and coincide with the edges of tomographically mapped high seismic velocities. Modelling using synthetic seismograms suggests that a narrow (10–15 km) low- or high-velocity lamella with about 2 per cent velocity contrast can reproduce the observed reflected waveforms, but other explanations may exist. Considering the reflection locations and synthetic modelling, the observed out-of-plane energy is well explained by underside reflections off a sharp reflector at the base of the subducted lithosphere. We also detect weaker reflections corresponding to the tomographically mapped top of the slab, which may arise from the boundary between the Nazca plate and the overlying former basaltic oceanic crust. The joint interpretation of the waveform modelling and geodynamic considerations indicate mass flux of the former oceanic lithosphere and basaltic crust across the 660 km discontinuity, linking processes and structure at the top and bottom of the Earth's mantle, supporting the idea of whole mantle convection.     

Ccte d'Ivoire-Ghana margin: seismic imaging of passive rifted crust adjacent to a transform continental margin

During May 1990 and January-February 1991, an extensive geophysical data set was collected over the Côte d'Ivoire-Ghana continental margin, located along the equatorial coast of West Africa. The Ghana margin is a transform continental margin running subparallel to the Romanche Fracture Zone and its associated marginal ridge—the Côte d'Ivoire-Ghana Ridge. From this data set, an explosive refraction line running ∼ 150 km, ENE-WSW between 3°55'N, 3°21'W and 4°23'N, 2°4'W, has been modelled together with wide-angle airgun profiles, and seismic reflection and gravity data. This study is centred on the Côte d'Ivoire Basin located just to the north of the Côte d'Ivoire-Ghana Ridge, where bathymetric data suggest that a component of normal rifting occurred, rather than the transform motion observed along the majority of the equatorial West African margin.
Traveltime and amplitude modelling of the ocean-bottom seismometer data shows that the continental Moho beneath the margin rises in an oceanward direction, from ∼ 24 km below sea level to ∼ 17 km. In the centre of the line where the crust thins most rapidly, there exists a region of anomalously high velocity at the base of the crust, reaching some 8 km in thickness. This higher-velocity region is thought to represent an area of localized underplating related to rifting. Modelling of marine gravity data, collected coincident with the seismic line, has been used to test the best-fitting seismic model. This modelling has shown that the observed free-air anomaly is dominated by the effects of crustal thickness, and that a region of higher density is required at the base of the crust to fit the observed data. This higher-density region is consistent in size and location with the high velocities required to fit the seismic data.     

Axial Intrusion Zone beneath the Median Valley of the Mid-Atlantic Ridge at 37°N Detected by Explosion Seismology

Summary. Four seismic refraction lines, three of which had shots every 250 m, were shot across, along and parallel to the median valley of the Mid-Atlantic Ridge at 37° N. A method has been developed for calculating the effect on the travel times of the rough sea-floor relief beneath the profiles and has been used to correct all the travel times for this effect. Most arrivals were from a main refractor of apparent velocity 5·4 to 6·3 km s⁻¹; only beyond 35 km were faster arrivals observed from an 8·09 ± 36 km s⁻¹ refractor. The main refractor corresponds in depth, at least approximately, to the top of Layer 3 of the ocean basins but its velocity is significantly less than normal for Layer 3, perhaps due to dip. A study of time residuals along two profiles across the median valley indicates the presence of a 2 to 3 km wide low velocity zone (about 3·2 km s⁻¹) beneath the median valley floor. This zone extends over the upper 2·5 km of the crust and is believed to represent a zone of intrusion through which magma passes on its way to the sea floor.     

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.     

The lithosphere–asthenosphere boundary beneath the western United States

S receiver functions from 67 broad-band seismic stations in the western United States clearly reveal the existence of a mantle discontinuity with velocity reduction downward, which we interpret as the lithosphere–asthenosphere boundary (LAB). The average depth of the LAB is ∼70 km. The boundary is relatively sharp with an overall sharpness of less than 20 km. The boundary is more prominent south of the Mendocino Triple Junction, where the Farallon Plate has completely subducted. This may indicate partial melts at the base of the lithosphere caused by the upwelling of the asthenospheric flow through the slab window. A double low velocity zone is observed at base of the lithosphere beneath southern Sierra Nevada, implying a second melting zone at a depth of ∼100 km, well correlated with previous studies of lithospheric delamination in the area.