Re-examination of the earth's free oxcillations excited by the Kamchatka earthquake of November 4, 1952

Benioff's suggestion that the 58-min period sinusoidal oscillation found on a Pasadena strain seismogram after the Kamchatka earthquake of November 4, 1952 may represent the earth's gravest normal mode is re-examined in terms of a slow large-scale post-seismic deformation. The mechanism and the seismic moment of the main shock of the Kamchatka earthquake are determined by using the amplitude and the initial phase of G₂ and R₂ recorded at Pasadena and R₆ recorded at Palisades. By constraining the dip angle and the strike of the fault at 30° (towards NW) and N34°E, respectively, on the basis of the geometry of the Benioff zone, the slip angle is determined as 110° which represents 74% thrust and 26% right-lateral faulting. The direction of the slip angle agrees with the slip direction of the Pacific plate. A seismic moment of 3.5 · 10²⁹ dyn cm is obtained. If a fault area of 650 · 200 km² is assumed, an average dislocation of 5 m is obtained. Spectral analyses of the Pasadena strain records show that the 58-min sinusoidal oscillation in fact consists of a spectral peak near 54 min which is very close to the ₀S₂ mode and other high-frequency peaks which can be correlated to the earth's normal modes. The records from two independent recording galvanometers correlate with each other very well, indicating that the recorded oscillation represents a real strain and not instrumental noise. The phase relation between the NS and EW components is consistent with the strain field associated with ₀S₂ mode. Although these results provide positive evidence for a slow post-seismic deformation, the cause of the abrupt termination of the oscillation and the excitation mechanism remain unresolved.     

A large normal-fault earthquake at the junction of the Tonga trench and the Louisville ridge

The source mechanism of a large (M_s ? 7.2) earthquake that occurred in the oceanic plate at the junction of the Tonga—Kermadec trench systems with the aseismic Louisville ridge is found by inverting long-period vertical-component Rayleigh waves recorded by the IDA network. The solution is an almost-pure normal fault, on a plane striking roughly parallel to the trench axis, with seismic moment of 1.7 × 10²⁷ dyn cm, and thus is among the ten largest documented shallow normal-fault earthquakes. A point-source depth of 20 km for the event is resolved by modeling teleseismic body waves; the actual rupture may have extended deeper, to 30 or 40 km. The earthquake was a multiple event, consisting of two sources separated by 16 s. A rupture velocity of 3.5 km s^?1 is inferred. The earthquake can be interpreted as tensional failure in the shallow portion of the downgoing plate caused by the gravitational pull of the slab. The Louisville ridge may be creating a local degree of decoupling of the oceanic plate from the overriding plate, and/or a zone of extension within the slab, which could enhance the effect of the gravitational forces in the shallower part of the downgoing plate. In particular, the earthquake could be associated with the break-up of the leading seamount of the ridge, which is currently right at the trench. Alternatively, the earthquake may have been caused by stresses associated with the bending of the plate prior to subduction.     

Spectral analysis of body waves for earthquakes and their source parameters in the Himalaya and nearby regions

The source characteristics of 33 earthquakes with magnitude m_b between 4.4 and 6.0, which occurred in the Himalayan and nearby regions, are investigated using the records of the Hyderabad seismograph station. The P- and S-wave spectra of these events are interpreted in terms of Brune's seismic source model for estimating the source parameters, i.e., seismic moment, source dimension, stress drop, average dislocation, apparent stress and the radiated energy. Seismic moments, M₀, vary between 0.3 × 10²⁴ and 9.0 × 10²⁶ dyne cm; source dimensions, r, between 4.3 and 18.6 km; stress-drops, Δσ between 0.3 and 151.6 bar; average dislocations, $\text{u}$ between 0.6 and 381 cm; apparent stresses, $\text{η}\text{σ}$ between 0.1 and 73.2 bar. The radiated energy, E_R is estimated by the spectrum integration method and is found to vary between 0.2 × 10¹⁸ and 9.3 × 10²² erg. In general, the stress drop and apparent stress are found to be high, indicating high stresses in these regions.     

Source process and tectonic implications of the Spanish deep-focus earthquake of March 29, 1954

The source process of the deep-focus Spanish earthquake of March 29, 1954 (m_b = 7.1, h = 630 km) has been studied by using seismograms recorded at teleseismic distances. Because of its unusual location, this earthquake is considered to be one of the most important earthquakes that merit detailed studies. Long-period body-wave records reveal that the earthquake is a complicated multiple event whose wave form is quite different from that of usual deep earthquakes. The total duration of P phases at teleseismic distances is as long as 40 s. This long duration may explain the considerable property damage in Granada and Malaga, Spain, which is rather rare for deep earthquakes. Using the azimuthal distribution of the differences between the arrival times of the first, the second and later P phases, the hypocenters of the later events are determined with respect to the first event. The focus of the second event is located on the vertical nodal plane of the first shock suggesting that this vertical plane is the fault plane. This fault plane which strikes in N2°E and dips 89.1°E defines a nearly vertical dip-slip fault, the block to the west moving downwards. The time interval and spatial separation between the first and the second events are 4.3 s and 19 km respectively, giving an apparent rupture velocity of 4.3 km/s which is about 74% of the S-wave velocity at the source. A third event occurred about 8.8 s after the first event and about 35.6 km from it. At least six to ten events can be identified during the whole sequence. The mechanism of some of the later events, however, seems to differ from the first two events. Synthetic seismograms are generated by superposition of a number of point sources and are matched with the observed signals to determine the seismic moment. The seismic moments of the later events are comparable to, or even larger than, that of the first. The total seismic moment is determined to be 7 · 10²⁷ dyn cm while the moments of the first and the second shocks are 2.1 · 10²⁶ dyn cm and 5.1 · 10²⁶ dyn cm, respectively. The earthquake may represent a series of fractures in a detached piece of the lithosphere which sank rapidly into the deep mantle preserving the heterogeneity of material property at shallow depths.     

Use of the mantle magnitudeM m for the reassessment of the moment of historical earthquakes

The mantle magnitudeM _m is used on a dataset of more than 180 wavetrains from 44 large shallow historical earthquakes to reassess their moments, which in many cases had been previously estimated only on the basis of the earthquake's rupture area. We provide 27 new or revised values ofM _o, based on the spectral amplitudes of surface waves recorded at a number of stations, principally Uppsala and Pasadena. Among them, and most significantly, we document a large low-frequency component to the source of the 1923 Kanto earthquake: the low-frequency seismic moment is 2.9×10²⁸ dyn-cm, in accord with geodetic observations. On the other hand, we revise downwards the seismic moment of the 1906 Ecuador event, which did not exceed 6×10²⁸ dyn-cm.Finally, the study of the 1960 Chilean and 1964 Alaskan earthquakes whose exceptionally large moments are properly retrieved throughM _m measurements, serves proof that this approach performs flawlessly even for the very greatest earthquakes, and is therefore successful in its goal to avoid the saturation effects plaguing any magnitude scale measured at a fixed period.     

Use of long-period surface waves for rapid determination of earthquake-source parameters

A method for rapid retrieval of earthquake-source parameters from long-period surface waves is developed. With this method, the fault geometry and seismic moment can be determined immediately after the surface wave records have been retrieved. Hence, it may be utilized for warning of tsunamis in real time. The surface wave spectra are inverted to produce either a seismic moment tensor (linear) or a fault model (nonlinear). The method has been tested by using the IDA (International Deployment of Accelerographs) records. With these records the method works well for the events larger than M_s = 6, and is useful for investigating the nature of slow earthquakes.For events deeper than 30 km, all of the five moment tensor elements can be determined. For very shallow events (d ? 30 km) the inversion becomes ill-conditioned and two of the five source moment tensor elements become unresolvable. This difficulty is circumvented by a two-step inversion. In the first step, the unresolvable elements are constrained to be zero to yield a first approximation. In the second step, additional geological and geophysical data are incorporated to improve the first approximation. The effect of the source finiteness is also included.     

On the observability of isotropic seismic sources: The July 31, 1970 Colombian earthquake

Theoretical calculations are made to study the observability of isotropic components of seismic sources. In particular we consider the 1970 deep Colombian earthquake, for which a precursory isotropic component was previously reported by Gilbert and Dziewonski.We compare an ultra-long period vertical record at Pasadena of the 1970 event to synthetic seismograms calculated both for Gilbert and Dziewonski's source model and for the pure double-couple source of Furumoto and Fukao, and obtain better overall agreement for the latter. The amplitude of the long-period synthetic for the isotropic source is about 5–15 times smaller than the synthetic for the deviatoric source, suggesting that the data may be relatively insensitive to the presence of a small isotropic source. When this possibility was tested, the overall agreement was found to be almost completely insensitive to the presence of even a reasonably large isotropic component.However, the isotropic source was derived from multi-station moment tensor inversion, rather than from single-station studies. A numerical experiment on the effect of lateral heterogeneity of eigenfrequencies and of Q on the inversion for the moment tensor shows that even relatively small amounts of heterogeneity can produce spurious isotropic sources from moment tensor inversion.     

Re-evaluation of the large deep earthquake of January 21, 1906

The large deep earthquake of January 21, 1906 is re-evaluated using old seismogram data and updated analysis techniques. From the P and pP-P time data the hypocentre parameters are determined as follows: origin time, 13h 49min 35s; latitude, 33.8°N; longitude, 137.5°E; depth, 340 km. The body-wave magnitude m_B is re-evaluated from the amplitude and periods of P, PP and S waves. The average value of 7.4 is obtained. This value is the smallest among any values assigned previously to this shock, and it is denied that the earthquake is the world's largest deep shock in this century. The focal mechanism is estimated from the P-wave first motions and amplitude distribution of P and S waves. Synthetic body waves are used to constrain the mechanism and to determine the seismic moment. The mechanism solution suggests the down-dip compression typical of this region. A seismic moment of 1.5 × 10²⁷ dyn · cm is obtained. This value and the re-evaluated value of m_B are consistent with the moment-_B relation obtained for other deep earthquakes.     

Source Characteristics of the 12 November 1996 Mw 7.7 Peru Subduction Zone Earthquake

—The 12 November 1996 M _w 7.7 Peru subduction zone earthquake occurred off the coast of southern Peru, near the intersection of the South American trench and the highest topographical point of the subducting Nazca Ridge. We model the broadband teleseismic P-waveforms from stations in the Global Seismic Network to constrain the source characteristics of this subduction zone earthquake. We have analyzed the vertical component P-waves for this earthquake to constrain the depth, source complexity, seismic moment and rupture characteristics. The seismic moment determined from the nondiffracted P-waves is 3–5 × 10²⁰ N·m, corresponding to a moment magnitude M _w of 7.6–7.7. The source time function for the 1996 Peru event has three pulses of seismic moment release with a total duration of approximately 45–50 seconds. The largest moment release occurs at approximately 35–40 seconds and is located ～90km southeast of the rupture initiation. Approximately 70% of the seismic moment was released in the third pulse.¶We find that the 1996 event reruptured part of the rupture area of the previous event in 1942. The location of the 1996 earthquake corresponds to a region along the Peru coast with the highest uplift rates of marine terraces. This suggests that the uplift may be due to repeated earthquakes such as the 1996 and 1942 events.     

Source process of the 1990 Gonghe,China, earthquake and tectonic stress field in the northeastern Qinghai-Xizang (Tibetan) plateau

TheM _s =6.9 Gonghe, China, earthquake of April 26, 1990 is the largest earthquake to have been documented historically as well as recorded instrumentally in the northeastern Qinghai-Xizang (Tibetan) plateau. The source process of this earthquake and the tectonic stress field in the northeastern Qinghai-Xizang plateau are investigated using geodetic and seismic data. The leveling data are used to invert the focal mechanism, the shape of the slipped region and the slip distribution on the fault plane. It is obtained through inversion of the leveling data that this earthquake was caused by a mainly reverse dip-slipping buried fault with strike 102°, dip 46° to SSW, rake 86° and a seismic moment of 9,4×10¹⁸ Nm. The stress drop, strain and energy released for this earthquake are estimated to be 4.9 MPa, 7.4×10^–5 and 7.0×10¹⁴ J, respectively. The slip distributes in a region slightly deep from NWW to SEE, with two nuclei, i.e., knots with highly concentrated slip, located in a shallower depth in the NWW and a deeper depth in the SEE, respectively.Broadband body waves data recorded by the China Digital Seismograph Network (CDSN) for the Gonghe earthquake are used to retrieve the source process of the earthquakes. It is found through moment-tensor inversion that theM _s =6.9 main shock is a complex rupture process dominated by shear faulting with scalar seismic moment of the best double-couple of 9.4×10¹⁸ Nm, which is identical to the seismic moment determined from leveling data. The moment rate tensor functions reveal that this earthquake consists of three consecutive events. The first event, with a scalar seismic moment of 4.7×10¹⁸ Nm, occurred between 0–12 s, and has a focal mechanism similar to that inverted from leveling data. The second event, with a smaller seismic moment of 2.1×10¹⁸ Nm, occurred between 12–31 s, and has a variable focal mechanism. The third event, with a sealar seismic moment of 2.5×10¹⁸ Nm, occurred between 31–41 s, and has a focal mechanism similar to that inverted from leveling data. The strike of the 1990 Gonghe earthquake, and the significantly reverse dip-slip with minor left-lateral strike-slip motion suggest that the pressure axis of the tectonic stress field in the northeastern Qinghai-Xizang plateau is close to horizontal and oriented NNE to SSW, consistent with the relative collision motion between the Indian and Eurasian plates. The predominant thrust mechanism and the complexity in the tempo-spatial rupture process of the Gonghe earthquake, as revealed by the geodetic and seismic data, is generally consistent with the overall distribution of isoseismals, aftershock seismicity and the geometry of intersecting faults structure in the Gonghe basin of the northeastern Qinghai-Xizang plateau.Contribution No. 96 B0006 Institute of Geophysics, State Seismological Bureau, Beijing, China.     

Comments on: “The partitioning of nickel between olivine and silicate melt” by S.R. Hart and K.E. Davis

We study a set of very high-quality records of first-order overtone Rayleigh waves from the deep-focus earthquake of September 29, 1973, in the Japan Sea. Standard surface wave techniques are used with these overtones, treated as individual seismic phases, to retrieve radiation pattern, Q, moment and phase velocity. A figure of M₀ = (6.7 ± 1.4) × 10²⁷dyn-cm is obtained, in total agreement with published values computed from either P waves, or fundamental Rayleigh waves. We also demonstrate the feasibility of using overtones as individual seismic phases in order to investigate their dispersion and attenuation properties.     

Analysis of seismological and tsunami data from the 1993 Guam earthquake

The fault parameters of the Guam earthquake of August 8, 1993 are estimated from seismological analyses, and the possibility of identifying the actual fault plane from tsunami waveforms is tested. The Centroid Moment Tensor solution of long-period surface waves shows one nodal plane shallowly dipping to the north and the other nodal plane steeply dipping to the south. The seismic moment is 3.5×10²⁰ Nm and the corresponding moment magnitude is 7.7. The Moment Tensor Rate Function inversion ofP waves also yields a similar focal mechanism and seismic moment. The point source depth is estimated as 40–50 km.This earthquake generated tsunamis that propagated toward the Japanese coast along the Izu-Bonin-Mariana ridge system. The tsunamis are recorded on ocean bottom pressure gauges and tide gauges. Numerical computation of tsunamis shows that the computed waveforms from the two possible fault planes match well with the observed tsunami waveforms. The numerical computation also shows that the tsunami waveforms at Guam Island, just above the fault, should contain useful information regarding the identification of the actual fault plane. However, the current sampling rate of the tide gauges is so small that the records cannot help the identification.     

A Study of Source Parameters, Site Amplification Functions and Average Effective Shear Wave Quality Factor Qseff from Analysis of Accelerograms of the 1999 Chamoli Earthquake, Himalaya

The accelerograms of the 1999 Chamoli earthquake and nine of its aftershocks, which occurred in Uttaranchal Himalaya, have been analyzed to investigate their source parameters, the site amplification functions and the average effective shear-wave quality factor Q_seff in the region. The fault plane solution of the main shock is obtained using the spectral amplitudes of SH waves (approximated by transverse components of accelerograms) of the high-energy packets observed in the accelerograms of the main shock. It is found to be comparable with the reported solutions in other studies. Similarly the other source parameters (viz., seismic moment = (5.03±1.7) × 10²⁵ dyne-cm, stress drop = 65 bars, source duration = 5.2 s and moment magnitude = 6.4) estimated for the main shock are consistent with the values obtained in other studies. The stress drops estimated for the aftershocks vary from 23 bars to 153 bars and the seismic moment from 1.4 × 10²³ dyne-cm to 2.9 × 10²³ dyne-cm. The average estimated values of the effective shear-wave quality factor Q_seff vary from 655±359 in the Uttaranchal sector of Himalaya and 1475±130 in the Delhi region. In general, the Q_seff value increases with an increase in the epicentral distance reflecting the penetration of the waves into deeper layers of the crust as the epicentral distance of the observation point increases. These values of Q_seff indicate that in general the curst is at low temperatures that will promote brittle behavior and conditions for episodic failure as compared to creep, under the accumulated strains from plate collision at the Himalaya plate boundary. The site amplification characteristics at sites have been identified from the frequency bands of significant amplification observed in the spectral ratios of the horizontal to the vertical component records. The decay of peak ground acceleration (PGA) values with distance has been investigated using the empirical regression curves vis-à-vis the site amplification factors.     

Estimation of Source Parameters for the Aftershocks of the 2001 Mw 7.7 Bhuj Earthquake, India

The source parameters for 213 Bhuj aftershocks of moment magnitude varying from 2.16 to 5.74 have been estimated using the spectral analysis of the SH- waveform on the transverse component of the three-componnet digital seismograms as well as accelerograms. The estimated stress drop values for Bhuj aftershocks show more scatter (Mo^{0.5 to 1} ∞ Δσ) toward the larger seismic moment values (log Mo ≥ 10^14.5 N-m, larger aftershocks), whereas, they show a more systematic nature (Mo³ ∞ Δσ) for smaller seismic moment (log Mo < 10^14.5 N-m, smaller aftershocks) values. This size dependency of stress drop has also been seen from the relation between our estimated seismic moment and source radius, however, this size-dependent stress drop is not observed for the source parameter estimates for the other stable continental region earthquakes in India and around the world. The estimated seismic moment (Mo), source radius (r) and stress drop (Δσ) for aftershocks of moment magnitude 2.16 to 5.74 range from 1.95 × 10¹² to 4.5 × 10¹⁷ N-m, 239 to 2835 m and 0.63 to 20.7 MPa, respectively. The near-surface attenuation factor (k) is found to be large of the order of 0.03 for the Kachchh region, suggesting thick low velocity sediments beneath the region. The estimated stress drop values show an increasing trend with the depth indicating the base of seismogenic layer (as characterized by larger stress drop values (>15 MPa)) lying in 22–26km depth range beneath the region. We suggest that the concentration of large stress drop values at 10–36km depth may be related to the large stress/strain associvated with a brittle, competent intrusive body of mafic nature.     

Applicability of some infiltration formulae to rangeland infiltrometer data

Three algebraic infiltration equations (Kostiakov's, Horton's and Philip's) were examined to determine which one would best fit infiltrometer data collected from a variety of mostly semi-arid rangeland plant communities from both Australia and the United States. Approximately 1,100 infiltrometer plots were included in the analysis. Results indicated that, in every instance, Horton's equation best fit the infiltrometer data. Variability of “point” measures of short-term infiltration rates were never satisfactorily accounted for by using either Kostiakov's or Philip's equation. Though Horton's equation provided a best fit to the overall infiltration data, R² values indicated a potential usefulness of this equation only under the certain conditions that were sampled in several rangeland plant communities in the Northern Territory, Australia. The equation could not be considered consistently useful under conditions sampled on rangelands in the United States.     

Iterative deconvolution of complex body waves from great earthquakes—the Tokachi-Oki earthquake of 1968

A method of body-wave inversion is developed in an attempt to extract the information about asperities or barriers in a fault zone. A sequence of point sources, each being characterized with the seismic moment, the onset time and the location, are iteratively derived from observed records at multi-stations, where the two-dimensional extent of the source location is taken into account. A modification is made of the iterative method of Kikuchi and Kanamori on the formulation of inversion procedure to facilitate the computation.Using this method, we analyse long period P waves of the Tokachi-Oki earthquake of 1968 (M_w = 8.2) and obtain several significant subevents with time durations of ～ 10 s. Their spatio-temporal distribution shows that the rupture process consists of three characteristic stages: (A) a stage of introductory rupture, (B) a stage of main rupture and (C) a stage of aftershocks. The main rupture takes place in the form of clustering around a few sites of the fault plane. The largest subevent occurs in the northwestern corner. The stress drop associated with this event is estimated to be ～ 200 bars, one order of magnitude higher than the stress drop averaged over the entire fault plane. The sum of the seismic moments of the individual subevents amounts to 2.3 × 10²⁸ dyn. cm which approximately coincides with the one estimated from the analysis of long-period surface waves. This implies that the source of the Tokachi-Oki earthquake consists of several major subevents with time durations of ～ 10 s in addition to other minor subevents.     

The 1957 great Aleutian earthquake

The 9 March 1957 Aleutian earthquake has been estimated as the third largest earthquake this century and has the longest aftershock zone of any earthquake ever recorded—1200 km. However, due to a lack of high-quality seismic data, the actual source parameters for this earthquake have been poorly determined. We have examined all the available waveform data to determine the seismic moment, rupture area, and slip distribution. These data include body, surface and tsunami waves. Using body waves, we have estimated the duration of significant moment release as 4 min. From surface wave analysis, we have determined that significant moment release occurred only in the western half of the aftershock zone and that the best estimate for the seismic moment is 50–100×10²⁰ Nm. Using the tsunami waveforms, we estimated the source area of the 1957 tsunami by backward propagation. The tsunami source area is smaller than the aftershock zone and is about 850 km long. This does not include the Unalaska Island area in the eastern end of the aftershock zone, making this area a possible seismic gap and a possible site of a future large or great earthquake. We also inverted the tsunami waveforms for the slip distribution. Slip on the 1957 rupture zone was highest in the western half near the epicenter. Little slip occurred in the eastern half. The moment is estimated as 88×10²⁰ Nm, orM _w =8.6, making it the seventh largest earthquake during the period 1900 to 1993. We also compare the 1957 earthquake to the 1986 Andreanof Islands earthquake, which occurred within a segment of the 1957 rupture area. The 1986 earthquake represents a rerupturing of the major 1957 asperity.     

The 15 August 2007 Peru Earthquake and Tsunami: Influence of the Source Characteristics on the Tsunami Heights

The tsunami caused by the 2007 Peru earthquake (M_w 8.0) provoked less damage than by the seismic shaking itself (numerous casualties due to the earthquake in the vicinity of Pisco). However, it propagated across the Pacific Ocean and small waves were observed on one tide gauge in Taiohae Bay (Nuku Hiva, Marquesas, French Polynesia). We invert seismological data to recover the rupture pattern in two steps. The first step uses surface waves to find a solution for the moment tensor, and the second step uses body waves to compute the slip distribution in the source area. We find the slip distribution to consist of two main slip patches in the source area. The inversion of surface waves yields a scalar moment of 8.9 10²⁰ Nm, and body-wave inversion gives 1.4 10²¹ Nm. The inversion of tsunami data recorded on a single deep ocean sensor also can be used to compute a fault slip pattern (yielding a scalar moment of 1.1 10²¹ Nm). We then use these different sources to model the tsunami propagation across the Pacific Ocean, especially towards Nuku Hiva. While the source model taken from the body-wave inversion yields computed tsunami waves systematically too low with respect to observations (on the central Pacific Ocean DART buoy as on the Polynesian tide gauge), the source model established from the surface-wave inversion is more efficient to fit the observations, confirming that the tsunami is sensitive to the low frequency component of the source. Finally we also discuss the modeling of the late tsunami arrivals in Taiohae Bay using several friction coefficients for the sea bottom.     

Sensitivity of the strong ground motion time histories to a finite source model: A case study for the January 12, 2010 Haiti earthquake (Mw=7.0)

We analyze the waveforms generated by the January 12, 2010 Haiti earthquake (M_w=7.0) for its source characteristics. A 60 to 25 km source model is retrieved by the Kikuchi and Kanamori finite source inversion technique that uses broadband teleseismic body wave records. The derived rupture model points out unilateral rupture propagation commenced at the eastern side of the fault plane where the major seismic moment release occurred. The rupture front propagated westward and terminated at a site where the largest aftershocks occurred. Our estimates yield a seismic moment of M_o=8.17×10¹⁹ N m released on a 60 km-long fault plane. A patch at the eastern side of the ruptured fault plane inferred as a region of maximum moment release.     

From 3-Hz P Waves to 0 S 2: No Evidence of A Slow Component to the Source of the 2011 Tohoku Earthquake

In the hours following the 2011 Honshu event, and as part of tsunami warning procedures at the Laboratoire de Géophysique in Papeete, Tahiti, the seismic source of the event was analyzed using a number of real-time procedures. The ultra-long period mantle magnitude algorithm suggests a static moment of 4.1 × 10²⁹ dyn cm, not significantly different from the National Earthquake Information Center (NEIC) value obtained by W-phase inversion. The slowness parameter, $\Uptheta = -5.65, $ is slightly deficient, but characteristic of other large subduction events such as Nias (2005) or Peru (2001); it remains significantly larger than for slow earthquakes such as Sumatra (2004) or Mentawai (2010). Similarly, the duration of high-frequency (2–4 Hz) P waves in relation to seismic moment or estimated energy, fails to document any slowness in the seismic source. These results were confirmed in the ensuing weeks by the analysis of the lowest-frequency spheroidal modes of the Earth. A dataset of 117 fits for eight modes (including the gravest one, ₀ S ₂, and the breathing mode, ₀ S ₀) yields a remarkably flat spectrum, with an average moment of 3.5 × 10²⁹ dyn cm (*/1.07). This behavior of the Tohoku earthquake explains the generally successful real-time modeling of its teleseismic tsunami, based on available seismic source scaling laws. On the other hand, it confirms the dichotomy, among mega-quakes (M ₀ > 10²⁹ dyn cm) between regular events (Nias, 2005; Chile, 2010; Sendai, 2011) and slow ones (Chile, 1960; Alaska, 1964; Sumatra, 2004; and probably Rat Island, 1965), whose origin remains unexplained.