Stress drops, wave spectra and recurrence intervals of great earthquakes — implications of the Etorofu earthquake of 1958 November 6

Summary . The great Etorofu earthquake of 1958 November 6 is characterized by a relatively small aftershock area (70 × 150 km²) and an extremely large felt area. The felt area is more extensive than those of any other large earthquakes which have occurred in the southern Kurile to northern Japan arc since the beginning of this century. The mechanism is a pure thrust fault typical of most great earthquakes in island arcs. A body wave magnitude of m _b = 8.2 is obtained at periods around 6 s using more than 40 observations, although an m _b value of only 7.6–7.7 would be expected empirically from the observed surface wave magnitude of M _s= 8.1–8.2. Both an unusually large felt area and a high m _b indicate a dominance of high-frequency components in the seismic waves. A seismic moment of M _o= 4.4 × 10²⁸ dyne cm is determined from long-period surface waves from which a high stress drop of Δσ = 78 bar is obtained using a relatively small aftershock area. Historic data indicate an anomalously long time interval between the 1958 event and any earlier great earthquake from the same source region. The observed high stress drop can be interpreted as a consequence of this long intervening period through which strain built up. The dominance of the high-frequency seismic waves can then be interpreted as a result of this high stress drop. Stress drops, seismic wave spectra and recurrence intervals of great earthquakes are in this way closely related to each other. The 1958 event may represent a high strength extreme of stochastic fluctuation of fracture strength relevant to great earthquakes.     

Coseismic well-level changes due to the 1992 Roermond earthquake compared to static deformation of half-space solutions

The M _w 5.4 Roermond earthquake of 1992 April 13 was one of the strongest events during the last 500 years in Central Europe. For the period March–May 1992, we collected records of 194 continuously operating well-level sensors, mostly located within 120  km of the epicentre. Nearly all wells penetrate unconfined or poorly confined Quaternary deposits with high hydraulic conductivities. 81 out of 194 raw data sets show a significant dynamic or step-like response of centimetre amplitude to the passage of seismic waves. Precursory anomalies are not obvious in these records. Coseismic well-level fluctuations could reflect a redistribution of stress and pore pressure in the brittle crust. Systematic analyses of such fluctuations may improve our knowledge of the role of pore fluids in crustal rheology and earthquake mechanics. The rather high number of individual observational records for a single event allows a regional correlation of the signs and amplitudes of the coseismic steps to changes in volume strain caused by the earthquake. The coseismic strain field at the surface was calculated for a homogeneous and a layered half-space. The results show reasonable agreement in the sign of the well-level steps but the amplitudes predicted from the wells' volumetric strain responses are much smaller than those that were recorded. Clearly, the coseismic well-level steps cannot be explained by volume strain changes, as derived from linear elastic models.     

Dynamic stress field of a kinematic earthquake source model with k-squared slip distribution

Source models such as the k -squared stochastic source model with k -dependent rise time are able to reproduce source complexity commonly observed in earthquake slip inversions. An analysis of the dynamic stress field associated with the slip history prescribed in these kinematic models can indicate possible inconsistencies with physics of faulting. The static stress drop, the strength excess, the breakdown stress drop and critical slip weakening distance D _c distributions are determined in this study for the kinematic k -squared source model with k -dependent rise time. Several studied k -squared models are found to be consistent with the slip weakening friction law along a substantial part of the fault. A new quantity, the stress delay, is introduced to map areas where the yielding criterion of the slip weakening friction is violated. Hisada's slip velocity function is found to be more consistent with the source dynamics than Boxcar, Brune's and Dirac's slip velocity functions. Constant rupture velocities close to the Rayleigh velocity are inconsistent with the k -squared model, because they break the yielding criterion of the slip weakening friction law. The bimodal character of D _c / D _tot frequency–magnitude distribution was found. D _c approaches the final slip D _tot near the edge of both the fault and asperity. We emphasize that both filtering and smoothing routinely applied in slip inversions may have a strong effect on the space–time pattern of the inferred stress field, leading potentially to an oversimplified view of earthquake source dynamics.     

Postseismic deformation at El Asnam (Algeria) in the seismotectonic context of northwestern Algeria

We use a combination of seismicity. tectonic features, focal mechanisms, seismic strain and postseismic movement to study the western part of North Algeria, the El Asnam region and its surrounding area in particular. A seismotectonic map of this part of Algeria, delimited by the Mediterranean Sea in the north and the Tellian mountains in the south, was built from available geological and seismological data. An examination of this map shows that the most significant earthquakes are concentrated along tectonic features and quaternary basins elongated in an east-west direction, suggesting NNW-SSE compressional movements. During the large El Asnam earthquake of 1980 October 10, M _w= 7.1, vertical movement was measured along a 40 km northeast-southwest thrust fault. These movements were determined geodetically in 1981 with reference to a basic network previously measured in 1976. In order to control postseismic movement and to ensure the monitoring of the seismic area, a dense geodetic network has been regularly measured since 1986, both in planemetry and altimetry. The results of the altimetric remeasurements show significant vertical movements. The elevation changes of the benchmarks have been deduced from precise levelling measurements: a remarkable uplift (5.1 ± 1.9 mm yr⁻¹) of the northwestern block, during the 1986-91 period has been observed, whereas the southeastern block is seen to be relatively stable. The Sar El Marouf anticline, situated along the central segment of the El Asnam surface breaks, appears to be growing with a maximum postseismic slip rate of (9.6 ± 1.4 mm yr⁻¹). The mean uplift rates computed for the northwestern block support the view that the 1954 earthquake did not occur on the same reverse fault as the 1980 event.     

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.     

Source parameters of large historical (1918–1962) earthquakes, South Island, New Zealand

We present the results of body waveform modelling studies for 17 earthquakes of M _w ≥5.7 occurring in the South Island, New Zealand region between 1918 and 1962, including the 1929 M _s = 7.8 Buller earthquake, the largest earthquake to have occurred in the South Island this century. These studies confirm the concept of slip partitioning in the northern South Island between strike-slip faulting in southwestern Marlborough and reverse and strike-slip faulting in the Buller region, but indicate that the zone of reverse faulting is quite localized. In the central South Island, all historical earthquakes appear to be associated with strike-slip faulting, although recent (post-1991) reverse faulting events suggest that slip partitioning also occurs within this region. The difference between historical and recent seismicity in the central South Island may also reflect stress readjustment occurring in response to the 1717 ad rupture along the Alpine fault. Within the Fiordland region (southwestern South Island) none of the historical earthquakes appears to have occurred along the Australian/Pacific plate interface, but rather they are associated with complex deformation of the subducting plate as well as with deformation of the upper (Pacific) plate. Two earthquakes in the Puysegur Bank region south of the South Island suggest that strike-slip deformation east of the Puysegur Trench is playing a major role in the tectonics of the region.     

Strain accumulation along an oblique plate boundary: the Reykjanes Peninsula, southwest Iceland

We use annual GPS observations on the Reykjanes Peninsula (RP) from 2000 to 2006 to generate maps of surface velocities and strain rates across the active plate boundary. We find that the surface deformation on the RP is consistent with oblique plate boundary motion on a regional scale, although considerable temporal and spatial strain rate variations are observed within the plate boundary zone. A small, but consistent increase in eastward velocity is observed at several stations on the southern part of the peninsula, compared to the 1993–1998 time period. The 2000–2006 velocities can be modelled by approximating the plate boundary as a series of vertical dislocations with left-lateral motion and opening. For the RP plate boundary we estimate left-lateral motion  18⁺⁴₋₃ mm yr⁻¹  and opening of  7⁺³₋₂ mm yr⁻¹  below a locking depth of  7⁺¹₋₂ km  . The resulting deep motion of  20⁺⁴₋₃ mm yr⁻¹  in the direction of  N(100⁺⁸₋₆)°E  agrees well with the predicted relative North America–Eurasia rate. We calculate the areal and shear strain rates using velocities from two periods: 1993–1998 and 2000–2006. The deep motion along the plate boundary results in left-lateral shear strain rates, which are perturbed by shallow deformation due to the 1994–1998 inflation and elevated seismicity in the Hengill–Hrómundartindur volcanic system, geothermal fluid extraction at the Svartsengi power plant, and possibly earthquake activity on the central part of the peninsula.     

Coseismic slip model of the 2007 August Pisco earthquake (Peru) as constrained by Wide Swath radar observations

The Pisco earthquake ( M _w 8.0; 2007 August 15) occurred offshore of Peru's southern coast at the subduction interface between the Nazca and South American plates. It ruptured a previously identified seismic gap along the Peruvian margin. We use Wide Swath InSAR observations acquired by the Envisat satellite in descending and ascending orbits to constrain coseismic slip distribution of this subduction earthquake. The data show movement of the coastal regions by as much as 85 cm in the line-of-sight of the satellite. Distributed-slip model indicates that the coseismic slip reaches values of about 5.5 m at a depth of ∼18–20 km. The slip is confined to less than 40 km depth, with most of the moment release located on the shallow parts of the interface above 30 km depth. The region with maximum coseismic slip in the InSAR model is located offshore, close to the seismic moment centroid location. The geodetic estimate of seismic moment is 1.23 × 10²¹ Nm ( M _w 8.06), consistent with seismic estimates. The slip model inferred from the InSAR observations suggests that the Pisco earthquake ruptured only a portion of the seismic gap zone in Peru between 13.5° S and 14.5° S, hence there is still a significant seismic gap to the south of the 2007 event that has not experienced a large earthquake since at least 1687.     

Seismic quiescence before the M 7, 1988, Spitak earthquake, Armenia

A detailed analysis of the 35  yr of seismicity between 1962 and 1997 using a gridding technique shows that the M 7, Spitak earthquake of 1988 December 7 was preceded by a quiescence anomaly that started at approximately 1984±0.5, and lasted about 5±0.5  yr, up to the main shock. This quiescence anomaly had a radius of about 20±3  km, estimated from circular areas with 75 per cent rate decrease, centred at the point of maximum significance of the anomaly. The quiescence was clearly present in the aftershock volume during the 5  yr before the 1988 main shock, but its statistically strongest expression was located 30  km NW of the epicentre. This anomaly fulfills the association rules between precursory quiescence anomalies and main shocks, even for a tight definition, and is therefore proposed as a case of precursory quiescence. The largest value of the standard deviate Z , found by random selection of samples by gridding, was Z =14 for a time window of T _w=3  yr, using a sample size of N =300 events. This makes this anomaly the strongest observed so far, and it is the first documented in an environment of continental collision. There are no false alarms exceeding in significance the precursor. The Armenian earthquake catalogue used for this study had 4600 earthquakes with M ≥ M _min=2.2 in the area bounded by 39.5° to 42°N/42.5° to 47°E. From the point of view of homogeneous reporting this is the best catalogue we have analysed so far. The limits of the data used and the density of the grid are dictated by the data, and have no influence on the results. The choice of free parameters does not influence the results significantly within the following limits: 100≤ N ≤500, 2≤ T _w≤7, 2.2≤ M _min≤2.8.     

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.     

Imaging an asperity of the 2003 Tokachi-oki earthquake using a dense strong-motion seismograph network

The 2003 Tokachi-oki earthquake ruptured a large area of approximately 100 km². The location of the largest asperity was estimated to be several dozen kilometres offshore of Hokkaido, Japan. The magnitude measured 8.0 on the Japan Meteorological Agency scale, and several studies used waveform inversion analysis to estimate the moment-magnitude as M _w 8.0–8.2. Several studies reported that there was a minor asperity at the northeastern edge of the fault plane, and that the rupture velocity towards the minor asperity was less than that towards the main asperity. One of them illustrated that the location and timing of the minor asperity were poorly constrained. In this paper, we introduce a procedure based on semblance analysis to image the location of the minor northeastern asperity with improved resolution. We group 15 strong-motion seismographs into three arrays, and we perform semblance analysis on impulsive waves that were possibly generated from the minor asperity and were conspicuously observed at stations in eastern Hokkaido. By projecting the semblance values onto the fault plane, we estimate the location of the minor asperity. We find it to be shallower and farther from the coast than the previous results indicated. The average rupture velocity towards the asperity is estimated to be 2.5 km s⁻¹, which is slower than the 3.6 km s⁻¹ obtained by waveform inversion analysis.     

Shear-wave anisotropy: spatial and temporal variations in time delays at Parkfield, Central California

Shear-wave splitting is analysed on data recorded by the High Resolution Seismic Network (HRSN) at Parkfield on the San Andreas fault, Central California, during the three-year period 1988-1990. Shear-wave polarizations either side of the fault are generally aligned in directions consistent with the regional horizontal maximum compressive stress, at some 70° to the fault strike, whereas at station MM in the immediate fault zone, shear-wave polarizations are aligned approximately parallel to the fault. Normalized time delays at this station are found to be about twice as large as those in the rock mass either side. This suggests that fluid-filled cracks and fractures within the fault zone are elastically or seismically different from those in the surrounding rocks, and that the alignment of fault-parallel shear-wave polarizations are associated with some fault-specific phenomenon.
Temporal variations in time delays between the two split shear-waves before and after a M_L = 4 earthquake can be identified at two stations with sufficient data: MM within the fault zone and VC outside the immediate fault zone. Time delays between faster and slower split shear waves increase before the M_L = 4 earthquake and decrease near the time of the event. The temporal variations are statistically significant at 68 per cent confidence levels. Earthquake doublets and multiplets also show similar temporal variations, consistent with those predicted by anisotropic poroelasticity theory for stress modifications to the microcrack geometry pervading the rock mass. This study is broadly consistent with the behaviour observed before three other earthquakes, suggesting that the build-up of stress before earthquakes may be monitored and interpreted by the analysis of shear-wave splitting.     

Source process and tectonic implications of the great 1975 North Atlantic earthquake

Summary. The Atlantic segment of the Africa–Europe plate boundary has usually been interpreted as a transform boundary on the basis of the bathymetric expression of the Gloria fault and dextral strike-slip first-motion mechanisms aligned along the Azores–Gibraltar line of seismicity. The 1975 May 26 earthquake ( M _s=7.9) was assumed to fit into this framework because it occurred in the general area of this line and has a similar first-motion focal mechanism (strike=288°, dip=72°, slip angle=184°). However, several anomalies cast doubt on this picture: the event is abnormally large for an oceanic transform event; a sizeable tsunami was excited; the aftershock area is unusually small for such a large event; and most significantly, the epicentre is 200 km south of the presumed plate boundary. The Rayleigh wave radiation pattern indicates a change in focal mechanism to one with a significant dip-slip component. The short duration of the source time history (20 s, as deconvolved from long-period P -waves), the lack of directivity in the Rayleigh waves, and the small one-day aftershock area suggest a fault length less than 80 km. One nodal plane of the earthquake is approximately aligned with the trace of an ancient fracture zone.
We have compared the Pasadena 1-90 record of the 1975 earthquake to that of the 1941 North Atlantic strike-slip earthquake (200 km to the NNW) and confirmed the large size of the 1941 event ( M =8.2). The non-colinear relationship of the 1975 and 1941 events suggests that there is no well-defined plate boundary between the Azores and Gibraltar. This interpretation is supported by the intraplate nature of both the 1975 event and the large 1969 thrust event 650 km to the east. This study also implies that the largest oceanic strike-slip earthquakes occur in old lithosphere in a transitional tectonic regime.     

New techniques for the analysis of earthquake sources from local array data with an application to the 1993 Scotts Mills, Oregon, aftershock sequence

We analysed aftershocks recorded by a temporary digital seismic network following the moderate M _w = 5.5 1993, Scotts Mills, Oregon, earthquake. A technique to retrieve source moment tensors from local waveforms was developed, tested, and applied to 41 small earthquakes ( M _w ranging from 1.6 to 3.2). The derived focal mechanisms, although well resolved, are highly variable and do not share a common nodal plane. In contrast, the majority of the events, relocated with a joint hypocentre determination algorithm, collapse to a well-focused plane. The incompatibility of the nodal planes of most events with the plane defined by their locations suggests that the aftershocks did not occur on the fault plane, but tightly around it, outlining the rupture area rather than defining it. Furthermore, the moment tensors reveal stable P -axes, whereas T  -axes plunges are highly dispersed. We detect a rotation of average T  -axis plunge with depth, indicating a change from shallower, predominantly dip-slip mechanisms to deeper strike-slip mechanisms. These characteristics are difficult to explain by remnant stress concentrations on the main-shock rupture plane or asperity- and barrier-type models. We suggest that the aftershocks occurred under the ambient regional stress, triggered by a sudden weakening of the region surrounding the main-shock slip, rather than from a shear stress increase due to the main shock.     

Seismic hazard in central southern Africa

Seismic hazard maps of central-southern Africa where hazard has been expressed in terms of peak ground acceleration for an annual probability in excess of 10^-1 show relatively high values that distinguish the seismic hazard potential of the Deka fault zone, the mid-Zambezi basin-Luangwa rift and western central Mozambique. In areas such as central-southern Africa where little is known about the geology of the region and the fault systems have not been fully mapped, seismic hazard potential may be estimated from seismicity and broad-scale fault features. For this region, such potential is based on earthquake magnitude M_s ≥ 6. Events of such magnitude have recently occurred in the mid-Zambezi basin, southern Zimbabwe and western-central Mozambique. This paper follows the conventional probabilistic hazard analysis procedure, defining seismic source zones from seismicity based on instrumental records from a cataloque that spans a period of 83 years. Geological and geomorphological features in the region are described on the mesoscale and are correlated with the seismicity as broad fault zones. The scarcity of strong-motion accelerogram data necessitated the formulation of attenuation values based on random vibration theory (RVT).     

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).     

A new derivation of the ratio Q(wobble)/Qμ (mantle)

Summary An extension of the Love-Larmor theory to a low-loss unelastic earth model, leads to the surprisingly simple approximation
   
where τ_s= 447.4 sidereal day is the static wobble period, τ_R= 306 sidereal day is the rigid-earth wobble period and τ_w= 433 sidereal day is the observed Chandler period. Q _W, Q _μ are the respective average Q values of the wobble and the Earth's mantle at τ_W. The known numerical factor F is only slightly dependent on the Earth structure.     

A conjugate strike-slip fault system within the extensional tectonics of Western Turkey

Three main shocks M-1, M-2 and M-3 (17 October 2005 at 05:45 UTC, M _w 5.4; 17 October at 09:46 UTC, M _w 5.8 and 20 October at 21:40 UTC, M _w 5.9) and their associated aftershocks within the Gulf of S i ğac i k, 50 km southwest of Izmir, Turkey were studied in detail. A temporary seismic network deployed during the activity allowed the hypocentre of M-3 and subsequent aftershocks to be determined with high accuracy. A relative relocation technique was used to improve the epicentres of M-1 and M-2. All three main shocks have strike-slip mechanisms which agree with the linear trends of the aftershock locations. Two distinct zones were illuminated by the aftershock locations. The zones contain clear echelon patterns with slightly different orientations from the trend of the aftershock distribution. M-2 and M-3 ruptured along of the eastern rupture zone which aligns N45°E. However the strike direction of M-1 is not clearly identified. The alignment of the two rupture zones intersect at their southern terminus at an angle of 90°. The fault zones form conjugate pair system and static triggering is considered as a probable mechanism for the sequential west to east occurrence of M-1, M-2 and M-3. This earthquake sequence provides seismological evidence for conjugate strike-slip faulting co-existing within a region dominated by north–south extension and well-developed east–west trending normal faults.     

On the use of satellite altimetry to infer the earthquake rupture characteristics: application to the 2004 Sumatra event

Standard data and methods, such as the inversion of seismic and GPS data, have been used extensively to infer the details of the 2004 December 26 earthquake. The unprecedented large size of this event gave the opportunity to modern altimeters to provide the first clear records of a tsunami in deep ocean, therefore allowing us to study the rupture history from an independent perspective. We invert the Jason-1 and Topex–Poseidon altimetry records, considering the new constraints available on the geometry of the fault plane, and taking them into account in a 3-D rupture model. The data are corrected for the non-negligible effect of satellite motion during measurements. Our results show that the rupture propagated over the 1500 km of subduction zone initially identified by the aftershock distribution, with a magnitude of   M _w= 9.1  . Our solution compares well with the latitudinal distribution of slip inferred from other data sets, with a maximum of energy release north of Sumatra, and two other slip patches near the Nicobar and Andaman islands. Based on waveform comparison, we assert that the shallow portion of the megathrust offshore Banda Aceh had slip amplitudes of more than 20 m. Also, we find that significant amounts of slip (about 10 m) concentrated below the Andaman islands and did not propagate on the shallow portion of the interface. Although synthetic tests tend to show less resolution in the northern part of the rupture, this solution is compatible with the near-field data (GPS, coral heads and imagery), and would allow one to explain the apparent paradox between the large local displacements and the moderate tsunami observed locally. Finally, we demonstrate the rapidly dominating effect of propagation and slip distribution over the rupture velocity, and how it precludes the direct estimate of this latter parameter.     

Comments on the Chandler wobble Q

Summary. The Chandler wobble Q _w, as obtained from the astronomical data cannot be equated with the Q _m of the source of damping, as an examination of Chandler wobble energetics reveals. We find that if dissipation occurs in the mantle then Q _w≃ 9 Q _m, implying that either the mantle Q is frequency dependent or the wobble Q is much larger than 100. If the dissipation is in the oceans then Q _w≃ 20 Q _o, and the pole tide must be far from equilibrium.