Seismic parameters controlling far-field tsunami amplitudes: A review

We present a review of the influence of various parameters of the sources of major oceanic earthquakes on the amplitude of tsunamis at transoceanic distances. We base our computations on the normal mode formalism, applied to realistic Earth models, but interpret our principal results in the simpler framework of Haskell theory in the case of a water layer over a Poisson half-space. Our results show that source depth and focal geometry play only a limited role in controlling the amplitude of the tsunami; their combined influence reaches at most 1 order of magnitude down to a depth of 150 km into the hard rock. More important are the effects of directivity due to rupture propagation along the fault, which for large earthquakes can result in a ten-fold decrease in tsunami amplitude by destructive interference, and the possibility of enhanced tsunami excitation in material with weaker elastic properties, such as sedimentary layers. Modelling of the so-called tsunami earthquakes suggests that an event for which 10% of the moment release takes place in sediments generates a tsunami 10 times larger than its seismic moment would suggest. We also investigate the properties of non-double couple sources and find that their relative excitation of tsunamis and Rayleigh waves is in general comparable to that of regular seismic sources. In particular, landslides involving weak sediments could result in very large tsunamis. Finally, we emphasize that the final amplitude at a receiving shore can be strongly affected by focusing and defocusing effects, due to variations in bathymetry along the path of the tsunami.     

Subduction and back-arc activity at the Hikurangi convergent Margin,New Zealand

The Hikurangi Margin is a region of oblique subduction with northwest-dipping intermediate depth seismicity extending southwest from the Kermadec system to about 42°S. The current episode of subduction is at least 16–20 Ma old. The plate convergence rate varies along the margin from about 60 mm/a at the south end of the Kermadec Trench to about 45 mm/a at 42°S. The age of the Pacific lithosphere adjacent to the Hikurangi Trench is not known.The margin divides at about latitude 39°S into two quite dissimilar parts. The northern part has experienced andesitic volcanism for about 18 Ma, and back-arc extension in the last 4 Ma that has produced a back-arc basin onshore with high heaflow, thin crust and low upper-mantle seismic velocities. The extension appears to have arisen from a seawards migration of the Hikurangi Trench north of 39°S. Here the plate interface is thought to be currently uncoupled, as geodetic data indicate extension of the fore-arc basin, and historic earthquakes have not exceededM _s=7.South of 39°S there is no volcanism and a back-arc basin has been produced by downward flexure of the lithosphere due to strong coupling with the subducting plate. Heatflow in the basin is normal. Evidence for strong coupling comes from historic earthquakes of up to aboutM _s=8 and high rates of uplift on the southeast coast of the North Island.The reason for this division of the margin is not known but may be related to an inferred increase, from northeast to southwest, in the buoyancy of the Pacific lithosphere.     

Active tectonics of the Himalaya

The Himalayan mountains are a product of the collision between India and Eurasia which began in the Eocene. In the early stage of continental collision the development of a suture zone between two colliding plates took place. The continued convergence is accommodated along the suture zone and in the back-arc region. Further convergence results in intracrustal megathrust within the leading edge of the advancing Indian plate. In the Himalaya this stage is characterized by the intense uplift of the High Himalaya, the development of the Tibetan Plateau and the breaking-up of the central and eastern Asian continent. Although numerous models for the evolution of the Himalaya have been proposed, the available geological and geophysical data are consistent with an underthrusting model in which the Indian continental lithosphere underthrusts beneath the Himalaya and southern Tibet. Reflection profiles across the entire Himalaya and Tibet are needed to prove the existence of such underthrusting. Geodetic surveys across the High Himalaya are needed to determine the present state of the MCT as well as the rate of uplift and shortening within the Himalaya. Paleoseismicity studies are necessary to resolve the temporal and spatial patterns of major earthquake faulting along the segmented Himalayan mountains.     

Comparison of amplitude decay rates in reflection,refraction, and local earthquake records

Three types of seismic data recorded near Coalinga, California were analyzed to study the behavior of scattered waves: 1) aftershocks of the May 2, 1983 earthquake, recorded on verticalcomponent seismometers deployed by the USGS; 2) regional refraction profiles using large explosive sources recorded on essentially the same arrays above; 3) three common-midpoint (CMP) reflection surveys recorded with vibrator sources over the same area. Records from each data set were bandpassed filtered into 5 Hz wide passbands (over the range of 1–25 Hz), corrected for geometric spreading, and fit with an exponential model of amplitude decay. Decay rates were expressed in terms of inverse codaQ (Q _c ^–1 ).Q _c ^–1 values for earthquake and refraction data are generally comparable and show a slight decrease with increasing frequency. Decay rates for different source types recorded on proximate receivers show similar results, with one notable exception. One set of aftershocks shows an increase ofQ _c ^–1 with frequency.Where the amplitude decay rates of surface and buried sources are similar, the coda decay results are consistent with other studies suggesting the importance of upper crustal scattering in the formation of coda. Differences in the variation ofQ _c ^–1 with frequency can be correlated with differences in geologic structure near the source region, as revealed by CMP-stacked reflection data. A more detailed assessment of effects such as the depth dependence of scattered contributions to the coda and the role of intrinsic attenuation requires precise control of source-receiver field geometry and the study of synthetic seismic data calculated for velocity models developed from CMP reflection data.     

Statistical analysis and long-term prediction of seismicity for linear zones

The model of the Poisson point process is too vague for earthquake locations in space and time: earthquakes tend to cluster in middle distances and to repulse in large ones. The Poisson point model with variable density makes it possible to describe the tendency for clustering but does not reveal the periodicity of clusters. The author proposes the point-process model where locations of points are determined not by densities of point distribution, but by densities of interpoint differences distribution. In the model, a latent periodicity is revealed and used for prediction of a point process. In 1983, the point-process model prediction was made for the Kuril Islands for 1983–1987 and two signs of danger in time and location were determined. Then they were confirmed by strong earth-quakes. In 1989, a similar prediction was made for North Armenia. The Spitak earthquake in 1988 is clearly seen from the data of previous earthquakes.     

Dynamic faulting under rate-dependent friction

We discuss the effects of rate-dependent friction on the propagation of seismic rupture on active faults. Several physicists using Burridge and Knopoff's box and spring model of faulting have proposed that fault complexity may arise from the spontaneous development of a self-similar stress distribution on the fault plane. If this model proves to be correct, it has important consequences for the origin of the complexity of seismic sources. In order to test these ideas on a more realistic earthquake model, we developed a new boundary integral equation method for studying rupture propagation along an antiplane fault in the presence of nonlinear rate-dependent friction. We study rupture dynamics of models with single and twin asperities. In our models, asperities are places on the fault with a higher value of prestress. Othewise all fault parameters are homogeneous. We show that for models with such asperities, a slip velocity weakening friction leads to the propagation of supersonic healing phases and to the spontaneous arrest of fracture if the prestress outside the asperities is low enough. For models with asperities, we can also observe narrow slip velocity pulses, qualitatively similar to the so-called Heaton pulses observed in some earthquake accelerograms. We also observe a complex distribution of stress after the rupture that depends on details of the initial distribution of asperities and on the details of the friction law.     

Seismic activity related to the 1991 eruption of Colima Volcano,Mexico

Ten years after the last effusive eruption and at least 15 years of seismic quiescence, volcanic seismic activity started at Colima volcano on 14 February 1991, with a seismic crisis which reached counts of more than 100 per day and showed a diversity of earthquake types. Four other distinct seismic crises followed, before a mild effusive eruption in April 1991. The second crisis preceded the extrusion of an andesitic scoriaceous lava lobe, first reported on 1 March; during this crisis an interesting temporary concentration of seismic foci below the crater was observed shortly before the extrusion was detected. The third crisis was constituted by shallow seismicity, featuring possible mild degassing explosion-induced activity in the form of hiccups (episodes of simple wavelets that repeat with diminishing amplitude), and accompanied by increased fumarolic activity. The growth of the new lava dome was accompanied by changing seismicity. On 16 April during the fifth crisis which consisted of some relatively large, shallow, volcanic earthquakes and numerous avalanches of older dome material, part of the newly extruded dome, which had grown towards the edge of the old dome, collapsed, producing the largest avalanches and ash flows. Afterwards, block lava began to flow slowly along the SW flank of the volcano, generating frequent small incandescent avalanches. The seismicity associated with the stages of this eruptive activity shows some interesting features: most earthquake foci were located north of the summit, some of them relatively deep (7–11 km below the summit level), underneath the saddle between the Colima and the older Nevado volcanoes. An apparently seismic quiet region appears between 4 and 7 km below the summit level. In June, harmonic tremors were detected for the first time, but no changes in the eruptive activity could be correlated with them. After June, the seismicity decreasing trend was established, and the effusive activity stopped on September 1991.     

The frictional behavior of lizardite and antigorite serpentinites: Experiments,constitutive models,and implications for natural faults

Laboratory studies of the frictional behavior of rocks can provide important information about the strength and sliding stability of natural faults. We have conducted friction experiments on antigorite and lizardite serpentinites, rocks common to both continental and oceanic crustal faults. We conducted both velocity-step tests and timed-hold tests on bare surfaces and gouge layers of serpentinite at room temperature. We find that the coefficient of friction of lizardite serpentinite is quite low (0.15–0.35) and could explain the apparent low stresses observed on crustal transform faults, while that of antigorite serpentinite is comparable to other crustal rocks (0.50–0.85). The frictional behavior of both types of serpentinite is well described by a two-mechanism model combining state-variable-dominated behavior at high slip velocities and flow-dominated behavior at low velocities. The two-mechanism model is supported by data from velocity-step tests and timed-hold tests. The low velocity behavior of serpentinite is strongly rate strengthening and should result in stable fault creep on natural faults containing either antigorite or lizardite serpentinite.