High-resolution structures of the Landers fault zone inferred from aftershock waveform data

High-frequency body waves recorded by a temporary seismic array across the surface rupture trace of the 1992 Landers, California, earthquake were used to determine fault-zone structures down to the seismogenic depth. We first developed a technique to use generalized ray theory to compute synthetic seismograms for arbitrarily oriented tabular low-velocity fault-zone models. We then generated synthetic waveform record sections of a linear array across a vertical fault zone. They show that both arrival times and waveforms of P and S waves vary systematically across the fault due to transmissions and reflections from boundaries of the low-velocity fault zone. The waveform characteristics and arrival-time patterns in the record sections allow us to locate the boundaries of the fault zone and to determine its P - and S -wave velocities independently as well as its depth extent. Therefore, the trade-off between the fault-zone width and velocities can be avoided. Applying the method to the Landers waveform data reveals a low-velocity zone with a width of 270–360 m and a 35–60 per cent reduction in P and S velocities relative to the host rock. The analysis suggests that the low-velocity zone extends to a depth of ∼7 km. The western boundary of the low-velocity zone coincides with the observed main surface rupture trace.     

Anomalous high-frequency wave propagation from the Tonga–Kermadec seismic zone to New Zealand

Summary. Bulletins of the International Seismological Centre (ISC) show very large residuals, up to 15 s early, for arrivals from events in the Tonga–Kermadec subduction zone to the New Zealand network of seismometers. The very early arrivals are confined to events south of about 22°S, and shallower than about 350 km. The waveforms show two distinct phases: an early, emergent, first phase with energy in the high-frequency band 2–10 Hz, and a distinct second phase, containing lower frequency energy, arriving at about the time predicted by JB tables.
The residuals are attributed to propagation through the cold, subducted lithosphere, which has a seismic velocity 5 per cent faster, on average, than normal. Ray tracing shows that the ray paths lie very close to the slab for events south of 22°S, but pass well beneath the slab for events further north, corresponding to the change in residual pattern. This characteristic of the ray paths is due to the curved shape of the seismic zone, and in particular to the bend in the zone where the Louisville ridge intersects the trench at 25°S.
The residuals can only be explained if the high velocity anomaly extends to a depth of 450 km in the region of the gap in deep seismicity from 32 to 36°S. The very high-frequency character of the first phase requires the path from the bottom of the slab to the stations to be of high Q , and to transmit 2–10 Hz energy with little attenuation.
The absence of low-frequency energy in the first phase is due to the narrowness of the high-velocity slab, which transmits only short-wavelength waves. The second phase, which contains low frequencies, is identified as a P -wave travelling beneath the subducted slab in normal mantle. There is no need to invoke any special structures, such as low-velocity waveguides or reflectors, to explain any of the observations. The S -wave arrivals show similar effects.     

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.     

Core-log integration studies in hole-A of Taiwan Chelungpu-fault Drilling Project

Taiwan Chelungpu-fault Drilling Project (TCDP) was initiated to understand the physical mechanisms involved in the large displacements of the 1999 Taiwan Chi-Chi earthquake. Continuous measurements of cores (including laboratory work) and a suite of geophysical downhole logs, including P - and S -wave sonic velocity, gamma ray, electrical resistivity, density, temperature, electrical borehole images and dipole-shear sonic imager, were acquired in Hole-A over the depth of 500–2003 m. Integrated studies of cores and logs facilitate qualitative and quantitative comparison of subsurface structures and physical properties of rocks. A total of 10 subunits were divided on the basis of geophysical characteristics. Generally, formation velocity and temperature increase with depth as a result of the overburden and thermal gradient, respectively. Gamma ray, resistivity, formation density, shear velocity anisotropy and density-derived porosity are primarily dependent on the lithology. Zones with changes of percentage of shear wave anisotropy and the fast shear polarization azimuth deduced from Dipole Shear-Imager (DSI) are associated with the appearance of fractures, steep bedding and shear zones. The fast shear wave azimuth is in good agreement with overall dip of the bedding (approximately 30° towards SE) and maximum horizontal compressional direction, particularly in the Kueichulin Formation showing strong shear wave velocity anisotropy. Bedding-parallel fractures are prevalent within cores, whereas minor sets of high-angle, NNW–SSE trending with N- and S-dipping fractures are sporadically distributed. The fault zone at depth 1111 m (FZA1111) is the Chi-Chi earthquake slip zone and could be a fluid conduit after the earthquake. The drastic change in fast shear wave polarization direction across the underlying, non-active Sanyi thrust at depth 1710 m reflects changes in stratigraphy, physical properties and structural geometry.     

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.     

The reflection seismic survey of project TIPTEQ—the inventory of the Chilean subduction zone at 38.2° S

We describe results of an active-source seismology experiment across the Chilean subduction zone at 38.2°S. The seismic sections clearly show the subducted Nazca plate with varying reflectivity. Below the coast the plate interface occurs at 25 km depth as the sharp lower boundary of a 2–5 km thick, highly reflective region, which we interpret as the subduction channel, that is, a zone of subducted material with a velocity gradient with respect to the upper and lower plate. Further downdip along the seismogenic coupling zone the reflectivity decreases in the area of the presumed 1960 Valdivia hypocentre. The plate interface itself can be traced further down to depths of 50–60 km below the Central Valley. We observe strong reflectivity at the plate interface as well as in the continental mantle wedge. The sections also show a segmented forearc crust in the overriding South American plate. Major features in the accretionary wedge, such as the Lanalhue fault zone, can be identified. At the eastern end of the profile a bright west-dipping reflector lies perpendicular to the plate interface and may be linked to the volcanic arc.     

Upper crustal shear velocity models from higher mode Rayleigh wave dispersion in Scotland

Summary. Group velocities for first and second higher mode Rayleigh waves, in the frequency range 0.8–4.8 Hz, generated from a local earthquake of magnitude 3.7 M _L in western Scotland, are measured at stations along the 1974 LISPB line. These provide detailed information about the crustal structure west of the line. The data divide the region into seven apparently homogeneous provinces. Averaged higher mode velocity dispersion curves for each province are analysed simultaneously using a linearized inversion technique, yielding regionalized shear velocity profiles down to a depth of 17 km into the upper crust. Shear wave velocity is between 3.0 and 3.4 km s⁻¹ in the upper 2 km, with a slow increase to around 3.8 km s⁻¹. P -wave models computed using these results agree with profiles from the LISPB and LUST refraction experiments.     

Coseismic deformation revealed by inversion of strong motion and GPS data: the 2003 Chengkung earthquake in eastern Taiwan

A moderate earthquake of   M _w= 6.8  occurred on 2003 December 10. It ruptured the Chihshang Fault in eastern Taiwan which is the most active segment of the Longitudinal fault as a plate suture fault between the Luzon arc of the Philippine Sea plate and the Eurasian plate. The largest coseismic displacements were 13 cm (horizontal) and 26 cm (vertical). We analyse 40 strong motion and 91 GPS data to model the fault geometry and coseismic dislocations. The most realistic shape of the Chihshang fault surface is listric in type. The dipping angle of the seismic zone is steep (about 60°–70°) at depths shallower than 10 km and then gradually decreases to 40°–50° at depths of 20–30 km. Thus the polygonal elements in Poly3D are well suited for modelling complex surfaces with curving boundaries. Using the strong motion data, the displacement reaches 1.2 m dip-slip on the Chihshang Fault and decreases to 0.1 m near surface. The slip averages 0.34 m, releasing a scalar moment of 1.6E26 dyne-cm. For GPS data, our model reveals that the maximal dislocation is 1.8 m dip-slip. The dislocations decrease to 0.1 m near the surface. The average slip is 0.48 m, giving a scalar moment of 2.2E26 dyne-cm. Regarding post-seismic deformation, a displacements of 0.5 m were observed near the Chihshang Fault, indicating the strain had not been totally released, as a probable result of near-surface locking of the fault zone.     

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

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

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.     

The three-dimensional shear wave structure in the mantle by overtone waveform inversion - I. Radial seismogram inversion

Summary. The three-dimensional (3-D) shear wave structure of the mantle, down to the depth of about 900 km, is obtained by inverting waveforms of radial component seismograms. Radial component seismograms contain large amplitude overtone signals which circle the Earth as wave packets and are sometimes called X1, X2, X3, … We use data which contain R1, X1 and X2 and filtered between 2 and 10mHz. It is shown that, unless each seismogram is weighted, all seismograms are not fitted uniformly. Only data from large earthquakes are fitted and the final velocity anomalies are biased by the small number of large earthquake data. Resolution is good at shallow depths, becomes worse in the intermediate depth range between about 400 and 500 km and then becomes better at greater depth ranges (600–900km). Even though we use only spheroidal mode data, velocity anomalies in the shallow structure show excellent correlation with the age of the surface rocks of the Earth. In the deeper regions, between about 600 and 900km, South America shows a fast velocity anomaly which may indicate the slab penetration beyond 700 km there. Another region which shows a fast velocity anomaly is the Mariana trench, but other subduction regions do not show such features.     

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.     

Viscoelastic relaxation and long-lasting after-slip following the 1997 Umbria-Marche (Central Italy) earthquakes

We combine Global Positioning System (GPS) measurements with forward modelling of viscoelastic relaxation and after-slip to study the post-seismic deformation of the 1997 Umbria-Marche (Central Apennines) moderate shallow earthquake sequence. Campaign GPS measurements spanning the time period 1999–2003 are depicting a clear post-seismic deformation signal. Our results favour a normal faulting rupture model where most of the slip is located in the lower part of the seismogenic upper crust, consistent with the rupture models obtained from the inversion of strong motion data. The preferred rheological model, obtained from viscoelastic relaxation modelling, consists of an elastic upper crust, underlain by a transition zone with a viscosity of 10¹⁸ Pa s, while the rheology of deeper layers is not relevant for the observed time-span. Shallow fault creep and after-slip at the base of the seismogenic upper crust are the first order processes behind the observed post-seismic deformation. The deep after-slip, below the fault zone at about 8 km depth, acting as a basal shear through localized time-dependent deformation, identifies a rheological discontinuity decoupling the seismogenic upper crust from the low-viscosity transition zone.     

Small-scale heterogeneities below the Gräfenberg array, Germany from seismic wavefield fluctuations of Hindu Kush events

Small-scale elastic heterogeneities (<5  km) are found in the upper lithosphere underneath the Gräfenberg array, southeast Germany. The results are based on the analysis of broadband recordings of 17 intermediate-depth (201–272  km) events from the Hindu Kush region. The wavefront of the first P arrival and the following 40  s coda are separated into coherent and incoherent (scattered) parts in the frequency range from 0.05 to 5  Hz. The frequency-dependent intensities of the mean and fluctuation wavefields are used to describe the scattering characteristics of the lithosphere underneath the receivers. It is possible to discriminate a weak-fluctuation regime of the wavefield in the frequency range below approximately 1.5–2.5  Hz and a strong-fluctuation regime starting at 2.0–2.5  Hz and continuing to higher frequencies. In order to explain the observed wavefield fluctuations, an approach with seismic scattering at random media-type structures is proposed. The preferred model contains heterogeneities with 3–7 per cent perturbations in seismic velocity and correlation lengths of 0.6–4.8  km in the crust. This is compatible with models from active seismic experiments. Scattering in the lithospheric mantle is not required, but cannot be excluded at weak velocity contrasts (<3 per cent).     

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.     

Size and duration of the high-frequency radiator in the source of the 2004 December 26 Sumatra earthquake

We recover the gross space–time characteristics of high-frequency (HF) radiator of the great Sumatra-Andaman islands earthquake of 2004 December 26 ( M _w= 9.1–9.3) using the time histories of the power of radiated HF P waves. To determine these time histories we process teleseismic P waves at 36 BB stations, using, in sequence: (1) bandpass filtering (four bands: 0.4–1.2, 1.2–2, 2–3 and 3–4 Hz); (2) squaring wave amplitudes, making 'power signals' for each band and (3) stripping the propagation-related distortion ( P coda, etc.) from the power signal and thus recovering source time function for HF power. In step (3) we employ an inverse filter constructed from an empirical Green's function, which is estimated as the power signal from an aftershock. For each ray we thus obtain signals with relatively well-defined end and no coda. From these signals we extract: total duration (joint estimate for all four bands) and temporal centroid of signal power for each band. Through linear inversion, the set of duration values for a set of rays delivers estimates of the rupture stopping point and stopping time. Similarly, the set of temporal centroids can be inverted to obtain the position of the space–time centroid of HF energy radiator. The quality of inversion for centroid is acceptable for lower-frequency bands but deteriorates for higher-frequency bands where only a fraction of stations provide useful data. For the source length and duration the following joint estimates were obtained: 1241 ± 224 km, 550 ± 10 s. The estimated stopping point position corresponds to the northern extremity of the aftershock zone. Spatial HF radiation centroids are located at distances 350–700 km from the epicentre, in a systematic way: the higher is the frequency, the farther is the centroid from the epicentre. Average rupture propagation velocity is estimated as 2.25 km s^–1.     

Comparison of the Upper Mantle Shear Wave Velocity Structure of the Indian Shield and the Tibetan Plateau and Tectonic Implications

Summary. Travel times and waveforms of long-period SH -waves recorded at distances of 10–30° and some SS waveforms are used to constrain the upper mantle velocities down to a depth of 400km beneath both the Indian Shield and the Tibetan Plateau. the shear velocity in the uppermost mantle beneath both the Indian Shield and the Tibetan Plateau is high and close to 4.7 km s⁻¹. the Indian Shield has a fairly thick high velocity lid, and the mean velocity between 40 and 250 km is between 4.58 and 4.68 km s⁻¹. In contrast, S -wave travel times and waveforms of S -waves, as well as a few for SS , show that the mean velocity between 70 and 250km beneath the central and northern part of the Tibetan Plateau is slower by 4 per cent or more than that beneath the Indian Shield and probably is between 4.4 and 4.5km s⁻¹. No large differences in the structure of the two areas below 250 km are required to explain both the arrival times and the waveforms of SH phases crossing Tibet or the Indian Shield. These results show that the structure of Tibet is not that of a shield and imply that the Indian plate is not underthrusting the whole of the Tibetan Plateau at the present time.     

Thrust faulting and crust—upper mantle structure in East Australia

Summary. The mid-crustal earthquake of 1973 March 9 (m_b= 5.5, h ≤ 20 km) located 60 km south-west of Sydney, Australia, provides unambiguous evidence of contemporary thrust faulting in South-eastern Australia — a region of high heat flow and Cenozoic basaltic volcanism. Aftershock locations suggest a steeply dipping fault in the depth range from 8 to 24 km with a lateral extent of about 8 km. The mechanism solution is consistent with a tectonic stress field that is dominated by east—west horizontal compression. A seismic moment of 5.7 ± 10²³± 20 per cent dyne-cm was computed from surface-wave amplitudes. Minimum values of slip and stress drop, 2 cm and 1 bar respectively, were estimated from the moment and a fault size taken' from aftershock locations.
Refinement modelling by a controlled Monte Carlo technique was used to provide unbiased models directly from multimode group velocities. The dispersion of fundamental and higher mode surface waves recorded at the digital high-gain station at Charters Towers, Queensland, and the WWSSN station at Adelaide, South Australia, is satisfied by crust- and upper-mantle models which have neither pronounced S-wave low-velocity zones nor thick high-velocity lids within 140 km of the Earth's surface. These models have subcrustal shear velocities of 4.20–4.32 km/s which are 0.4–0.5 km/s slower than Canadian shield shear velocities (CANSD).     

Scattering parameters of the lithosphere below the Massif Central, France, from teleseismic wavefield records

Teleseismic P -wave recordings are analysed in the frequency range 0.3–6  Hz to derive structural (statistical) parameters of the lithosphere underneath the French Massif Central. For this we analyse differences in frequency-dependent intensities of the mean wavefield and the fluctuation wavefield. It is possible to discriminate a weak fluctuation regime of the wavefield in the frequency range below 1  Hz and a strong fluctuation regime starting above 1  Hz and continuing to higher frequencies. The observed wavefield fluctuations in the frequency range 0.3–3  Hz can be explained by scattering of the teleseismic P wave front at elastic inhomogeneities in the lithosphere. A statistical distribution of the inhomogeneities is assumed and the concept of random media is applied. The lithospheric structure under the Massif Central can be described as a 70  km thick heterogeneous layer with velocity fluctuations of 3–7 per cent and correlation lengths of the heterogeneities of 1–16  km.