The deep seismicity of the Tyrrhenian Sea

The study reappraises the deep seismicity of the Tyrrhenian Sea. Careful examination of the quality of reported hypocentres shows that the earthquakes define a zone dipping NW, about 200 km along strike, 50 km thick, and reaching a depth of about 500 km. The zone is slightly concave to the NW at a depth of 300 km, but, contrary to many previous reports, is not tightly concave, nor are there significant spatial gaps in the seismicity, which is effectively continuous with depth. Seismicity is, however, concentrated in the depth interval 250–300 km, where the dip of the seismic zone changes from 70° (above 250 km) to a more gentle dip of 45° at greater depths. Seven fault-plane solutions are available for the largest earthquakes in this depth interval, all of them consistent with a P -axis down the dip of the seismic zone, and all of them requiring movement on faults out of the plane of the subducting slab.
Two deep earthquakes near Naples lie well outside the main zone of activity; for one of which a fault-plane solution is available that has a P -axis not aligned with the dip of the seismic zone. The tightly concave slab-geometry favoured by other reports is supported mainly by the location of these events near Naples, which we think may represent deformation in a separate, probably shallower dipping, piece of subducted lithosphere.
The lack of shallow seismicity, and particularly of thrust faulting earthquakes, at the surface projection of the Benioff zone suggests that active subduction has ceased. Estimates of the convergence rate responsible for subduction in the last 10 Myr far exceed the present convergence rate of Africa and Eurasia, suggesting that the subduction was related instead to the stretching and thinning of the crust in the Tyrrhenian Sea.     

Microearthquakes near the eastern end of St Paul's Fracture Zone

Summary . Four ocean-bottom seismographs were deployed near the eastern end of the St Paul's Fracture Zone in 1974 December. Microearthquakes were observed both along the fracture zone and in the median valley of the Mid-Atlantic Ridge. Seventy-six of them have been located and reliable depths obtained for 51. The range of depths observed suggests that the thickness of the lithosphere close to the ridge axis is 7 km. The absence of earthquakes on the ridge axis between 1 and 5 km depth may be the result of a highly cracked crust and thus indicates the depth to which hydrothermal fluids penetrate.     

A catalogue of seismicity in Greece and adjacent areas

Summary. A new earthquake catalogue for Greece has been formed to cover the instrumental period 1901–78, in particular 605 earthquakes for the period 1917–63 inclusive are relocated using first arrival data from the International Seismological Summary. These relocations incorporate macro-seismically and other well-controlled master events into an ensuing joint epicentre determination technique. The largest annual average shift is 165 km for 78 earthquakes during the decade after 1917, decreasing to 17 km for the later years 1957–63. Magnitudes are redetermined mainly using readings from the Swedish network and Uppsala since as early as 1908. Catalogue completeness exists for magnitudes around 5.5 for at least the last 60 years, but magnitudes of about 4.7 are completely reported only during the most recent 15 years.
The new epicentral positions lead to a better delineation of the seismic zones than has previously been achieved. The majority of shallow earthquakes are contained in a belt parallel to the Hellenic Arc which extends north into Albania, in the south-east the epicentral locations give a more diffuse extension of this zone into the west coast of Turkey. Depths of intermediate earthquakes along the Hellenic Arc tend to relocate to shallower depths, in the south-west part of Crete no earthquake focus deeper than 100 km is found, although in the south-eastern section of the arc a tendency to increased depths is observed. A second zone starts at Leukas Island in the west and extends through central Greece to near Volos on the east coast, where it divides into two branches, which are less well defined, but which eventually join the seismicity of western Turkey. A third zone follows the Saronikos and Corinth gulfs.     

The Nootka Fault Zone — a new plate boundary off western Canada

Summary. Relative motion across a boundary between the main Juan de Fuca plate and its northern extension, the Explorer plate, had earlier been suggested from sea-floor magnetic anomaly analysis and from earthquakes recorded on the western Canada land seismic network. The location of the boundary, called the Nootka fault zone, and the motion across it have been examined through seismic reflection profiles, accurate location of earthquakes with an array of ocean bottom seismometers and through analysis of magnetic, gravity and bathymetric data. The fault zone extends from a ridge-fault—fault triple point at the northern end of the Juan de Fuca ridge to a fault—trench—trench triple junction at the margin off north-central Vancouver Island. The active portion of the fault zone is about 20 km wide, and has produced extensive disturbance in the 0.5 to 1 km of overlying sediments. Magnetic anomaly analysis suggests present left-lateral strike slip motion of about 3 cm/yr, with convergence at the margin being more rapid to the south than to the north of the fault zone. Because of rapidly changing spreading parameters on the Explorer and Juan de Fuca ridges over the past 5 Myr the Nootka fault zone has had a very complex history.     

Regional mid-ocean stress and a proposed focal mechanism of 'stress-discordant' earthquakes

Summary. Rock stress measurements in Iceland show maximum horizontal compression perpendicular to the trend of Reykjanes Ridge crest and of its extension, the active volcanic zone of Iceland. Fault-plane solutions of dormant stage earthquakes are consistent with the measured stress orientations, but strike—slip earthquakes associated with volcanic surges and some earthquake swarms in active geothermal areas exhibit apparent reversals of mechanism and are here defined as 'stress-discordant' in the sense that they yield deduced stress orientations 90° from the regional stress field as determined by hydrofracturing and strain relief methods. It is proposed, supported by comparison with the pore-pressure induced Denver earthquakes, that the 'stress-discordant' volcanic earthquakes are triggered by increased pore pressure and probably involve stick-slip motion similar to that reported for some laboratory tests of the pore pressure effect, characterized by gradual onset and sudden stopping of each slip episode. The question is raised as to whether stress-discordant earthquakes are dominated by a stopping phase or terminal shock with consequent reversal of the deduced shear couple. A possible stopping mechanism is suggested: the dilatant stiffening of fault gouge during shear.
It is proposed that direct measurements of stress orientation be made by hydrofracturing tests at other places along the mid-ocean ridge crest and on the margins of the Red Sea and East African rifts. The Icelandic stress data indicate the need for sceptical re-examination of some fundamentals of plate tectonics theory.     

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.     

Seismotectonics and analysis of earthquakes from the Hindukush region

Summary. Over 80 earthquakes, exclusively from the Hindukush focal region, which were recorded at the Gauribidanur seismic array (GBA) have been used in this study. These events have similar epicentral distances and a narrow azimuthal range from GBA but varying focal depths from 10 to 240 km. A fault plane dipping steeply (75°) in the north-west direction and striking N 66° E has been investigated on the basis of the spatial distribution of earthquakes in two vertical planes through 68° E and 32° N. Short period P -wave recordings up to 30 s were processed using the adaptive cross-correlation filtering technique. Slowness and azimuthal anomalies were obtained for first arrivals. These anomalies show positive as well as negative bias and are attributed to a steep velocity gradient in the upper mantle between the 400–700 km depth range where the seismic rays have their maximum penetration. Relative time residuals between the stations of GBA owe their origin very near to the surface beneath the array. A search of the signals across the array revealed that most of the events occurring at shallower depths had complex signatures as compared to the deeper events. The structure near the source region, complicated source functions and the scattering confined to the crust—upper mantle near source are mainly responsible for the complexity of the Hindukush earthquakes as the transmission zone of the ray tubes from turning point to the recording station is practically the same.     

Imaging the Hikurangi subduction zone, New Zealand, using teleseismic receiver functions: crustal fluids above the forearc mantle wedge

We image the Hikurangi subduction zone using receiver functions derived from teleseismic earthquakes. Migrated receiver functions show a northwest dipping low shear wave feature down to 60 km depth, which we associate with the crust of the subducted Pacific Plate. Receiver functions (RF) at several stations also show a pair of negative and positive polarity phases with associated conversion depths of ∼20–26 km, where the subducted Pacific Plate is at a depth of ∼40–50 km beneath the overlying Australian Plate. RF inversion solutions model these phases with a thin low S -wave velocity zone less than 4 km thick, and an S -wave velocity contrast of more than ∼0.5 km s⁻¹ with the overlying crust. We interpret this phase pair as representing fluids near the base of the lower crust of the Australian Plate, directly overlying the forearc mantle wedge.     

Upper mantle shear structure of North America

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

Aseismic fault movement before the 1995 Kobe earthquake detected by a GPS survey: implication for preseismic stress localization?

By inversion analysis of the baseline changes and horizontal displacements observed with GPS (Global Positioning System) during 1990–1994, a high-angle reverse fault was detected in the Shikoku-Kinki region, southwest Japan. The active blind fault is characterized by reverse dip-slip (0.7±0.2  m yr⁻¹ within a layer 17–26  km deep) with a length of 208±5  km, a (down-dip) width of 9±2  km, a dip-angle of 51°±2° and a strike direction of 40°±2° (NE). Evidence from the geological investigation of subfaults close to the southwestern portion of the fault, two historical earthquakes ( M _L=7.0, 1789 and 6.4, 1955) near the centre of the fault, and an additional inversion analysis of the baseline changes recorded by the nationwide permanent GPS array from 18 January to 31 December 1995 partially demonstrates the existence of the fault, and suggests that it might be a reactivation of a pre-existing fault in this region. The fact that hardly any earthquakes ( M _L>2.0) occurred at depth on the inferred fault plane suggests that the fault activity was largely aseismic. Based on the parameters of the blind fault estimated in this study, we evaluated stress changes in this region. It is found that shear stress concentrated and increased by up to 2.1 bar yr⁻¹ at a depth of about 20  km around the epicentral area of the 1995 January 17  Kobe earthquake ( M _L=7.2, Japan), and that the earthquake hypocentre received a Coulomb failure stress of about 5.6 bar yr⁻¹ during 1990–1994. The results suggest that the 1995  Kobe earthquake could have been induced or triggered by aseismic fault movement.     

Temporal and three-dimensional spatial analyses of the frequency–magnitude distribution near Long Valley Caldera, California

The 3-D distribution of the b value of the frequency–magnitude distribution is analysed in the seismically active parts of the crust near Long Valley Caldera, California. The seismicity is sampled in spherical volumes, containing N =150 earthquakes and centred at nodes of a grid separated by 0.3  km. Significant variations in the b value are detected, with b ranging from b ≈0.6 to b ≈2.0. High b -value volumes are located near the resurgent dome, and at depths below 5  km at Mammoth Mountain. b values are found to be much lower south of the Long Valley Caldera. We interpret this to indicate that an active magma body has advanced from depths below 8  km to depths of 4 to 5  km beneath Mammoth Mountain in 1989, and that anomalous crust, either highly fractured or containing unusually high pore pressure, such as is the case in the vicinity of active magma bodies, exists north of the seismically active area beneath the resurgent dome at all depths. We also investigate the spatial distribution of temporal variations of the frequency–magnitude distribution by introducing differential b -value maps. b values increased from b ≈0.8 to b ≈1.5 underneath Mammoth Mountain at the onset of the 1989 earthquake swarm and remained high thereafter. This suggests that an intrusion permanently altered the average distribution of cracks at 5–10  km depth, or that the pore pressure permanently increased. We propose that high b values are a necessary (but not sufficient) condition near a magmatic body, and therefore spatial b -value mapping can be used to aid in the identification of active magma bodies.     

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.     

The evolution of the Gulf of Corinth (Greece): an aftershock study of the 1981 earthquakes

Summary. A preliminary study of the aftershocks of three earthquakes that occurred near to Corinth (Greece) in 1981 is combined with observations of the morphology and faulting to understand the evolution of the Eastern Gulf of Corinth. The well located aftershocks form a zone 60km long and 20km wide. They do not lie on the main fault planes and are mostly located between the north-dipping faulting on which the first two earthquakes occurred and the south-dipping faulting associated with the third event. A cluster of aftershocks also lies in the footwall of the eastern end of the south-dipping fault of the third event.
Morphologically, it is observed that in the evolution of the Eastern Gulf of Corinth, antithetic faulting apparently predates the appearance of the main faulting at the surface. This evolution can be explained by motion on a deep seated, shallow angle, aseismic normal fault. A model based on such a fault also accounts for the aftershock distribution of the 1981 earthquakes.     

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.     

An earthquake sequence studied with ocean-bottom seismographs

Summary. A tripartite ocean-bottom seismograph array at the junction of the East Pacific Rise and Rivera Fracture Zone recorded an eathquake sequence, consisting of three main shocks ( m _B= 4.3, 4.3 and 4.8) and numerous aftershocks from the fracture zone, in the distance range 35–50 km. Delineation of the rupture zones by aftershocks indicates that the first two main shocks took place on overlapping fault areas, while the third occurred over a fault area separated from the first by several kilometres. Both rupture zones were about 4 km long. Surface wave spectra indicate a shallow (about 3 km below the sea floor) source, as does OBS array phase velocity data. The seismic moments, obtained from teleseismic surface wave data, of 1.3, 2.1 and 2.8 × 10²³ dyn cm, with the fault areas as delineated by aftershocks, imply a stress drop of about 8 bars for the main shocks. Aftershock sequences of each of the main shocks are similar, with a b -value of about 0.65. Teleseismic P travel times are similar to those from near-surface sources in Nevada.     

Glacial rebound of the British Isles—II. A high-resolution, high-precision model

Observations of ice movements across the British Isles and of sea-level changes around the shorelines during Late Devensian time (after about 25 000 yr BP) have been used to establish a high spatial and temporal resolution model for the rebound of Great Britain and associated sea-level change. The sea-level observations include sites within the margins of the former ice sheet as well as observations outside the glaciated regions such that it has been possible to separate unknown earth model parameters from some ice-sheet model parameters in the inversion of the glacio-hydro-isostatic equations. The mantle viscosity profile is approximated by a number of radially symmetric layers representing the lithosphere, the upper mantle as two layers from the base of the lithosphere to the phase transition boundary at 400 km, the transition zone down to 670 km depth, and the lower mantle. No evidence is found to support a strong layering in viscosity above 670 km other than the high-viscosity lithospheric layer. Models with a low-viscosity zone in the upper mantle or models with a marked higher viscosity in the transition zone are less satisfactory than models in which the viscosity is constant from the base of the lithosphere to the 670 km boundary. In contrast, a marked increase in viscosity is required across this latter boundary. The optimum effective parameters for the mantle beneath Great Britain are: a lithospheric thickness of about 65 km, a mantle viscosity above 670 km of about (4-5) 10²⁰ Pa s, and a viscosity below 670 km greater than 4 × 10²¹ Pa s.     

A hypothesis for the seismogenesis of a double seismic zone

The seismogenesis of a double seismic zone, in particular the lower layer of a double seismic zone, has not been adequately explained in the literature. On the basis of seismic data and geothermal structures along three well-studied cross-sections in the Kuril-Kamchatka and Japan subduction zones, we investigate the temperature/pressure conditions associated with seismogenic structures of the double seismic zones. the corresponding T/P loci seem to suggest that earthquakes observed in the lower layer and in the lower part (below approximately 130 ± 20 km) of the top layer of a double seismic zone were caused by metastable phase transition-a mechanism similar to that responsible for deep-focus earthquakes only at lower temperature/pressure conditions. Under this hypothesis, the wedge-shaped configuration of a double seismic zone is interpreted to represent the loci of the kinetic boundary of the phase transition. According to theoretical/experimental studies and the constraints imposed by our observations, a likely candidate for such a phase transition is the metastable Al-rich enstatite decomposing into the assemblage of Al-poor enstatite plus garnet. Earthquakes in the upper part of the top layer were most probably due to conventional mechanisms such as dehydration of subducted materials and/or facies change from basalt to eclogite. That the top layer involves more than one seismogenic mechanism is also implied by the distinct behaviour of seismicity in the vicinity of 130 ± 20 km. Because the presence of deviatoric stress is critical to the reaction rate of a metastable phase transition, it is inferred that single seismic zones are also caused by the same mechanisms, except that the implicit layer of a supposed double seismic zone is missing, due to the insufficient amount of appropriate metastable minerals or to the lack of appropriate deviatoric stresses in the source region.     

The transition from subduction to continental collision: crustal structure in the North Canterbury region, New Zealand

The North Canterbury region marks the transition from Pacific plate subduction to continental collision in the South Island of New Zealand. Details of the seismicity, structure and tectonics of this region have been revealed by an 11-week microearthquake survey using 24 portable digital seismographs. Arrival time data from a well-recorded subset of microearthquakes have been combined with those from three explosions at the corners of the microearthquake network in a simultaneous inversion for both hypocentres and velocity structure. The velocity structure is consistent with the crust in North Canterbury being an extension of the converging Chatham Rise. The crust is about 27 km thick, and consists of an 11 km thick seismic upper crust and 7 km thick seismic lower crust, with the middle part of the crust being relatively aseismic. Seismic velocities are consistent with the upper and middle crust being composed of greywacke and schist respectively, while several lines of evidence suggest that the lower crust is the lower part of the old oceanic crust on which the overlying rocks were originally deposited.
The distribution of relocated earthquakes deeper than 15 km indicates that the seismic lower crust changes dip markedly near 43S. To the south-west it is subhorizontal, while to the north-east it dips north-west at about 10. Fault-plane solutions for these earthquakes also change near 43S. For events to the south, P -axes trend approximately normal to the plate boundary (reflecting continental collision), while for events to the north, T -axes are aligned down the dip of the subducted plate (reflecting slab pull). While lithospheric subduction is continuous across the transition, it is not clear whether the lower crust near 43S is flexed or torn.     

Amplification of ground motion and waveform complexity in fault zones: examples from the San Andreas and Calaveras Faults

Summary. P -wave seismograms at ranges less than 10 km are synthesized by asymptotic ray theory and by summation of Gaussian beams for point sources located in a low-velocity wedge surrounding a fault. The computations are performed using models of the wedge inferred from the analysis of reflection and refraction experiments across the San Andreas and Hayward-Calaveras faults. Calculations in these models show that the 10–20Hz vertical displacements of earthquakes located at 3–10km depth are amplified by up to an order of magnitude in a 1–2km wide region centred on the fault trace compared to displacements predicted by laterally homogeneous models of the crust. This amplification is not cancelled by high attentuation in the fault zone and compensates for the reduction in amplitudes directly above the source predicted from the radiation pattern of a strike-slip earthquake. Depending on the source depth of the earthquake and the structure and velocity contrast of the wedge, multiple triplications in the travel-time curve of direct P - and S -waves will occur at stations in the fault zone. A wedge model successfully predicts the triplications observed in the P waveforms of aftershocks of the Coyote Lake earthquake recorded in the fault zone, showing that body waves from microearthquakes can be used to determine the three-dimensional velocity structure of the fault zone. The amplification, waveform complexity, and distortion of ray paths introduced by the low- velocity wedge suggest that its effects should be included in the interpretation of strong ground motions and travel times observed in the fault zone. For realistic models of the wedge, asymptotically approximate methods of calculating the body waveforms are strictly valid for frequencies greater than 20Hz. Numerical methods may be necessary to calculate accurately the wavefield at lower frequencies.     

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.