Frictional afterslip following the 2005 Nias-Simeulue earthquake, Sumatra

Continuously recording Global Positioning System stations near the 28 March 2005 rupture of the Sunda megathrust [moment magnitude (Mw) 8.7] show that the earthquake triggered aseismic frictional afterslip on the subduction megathrust, with a major fraction of this slip in the up-dip direction from the main rupture. Eleven months after the main shock, afterslip continues at rates several times the average interseismic rate, resulting in deformation equivalent to at least a M(w) 8.2 earthquake. In general, along-strike variations in frictional behavior appear to persist over multiple earthquake cycles. Aftershocks cluster along the boundary between the region of coseismic slip and the up-dip creeping zone. We observe that the cumulative number of aftershocks increases linearly with postseismic displacements; this finding suggests that the temporal evolution of aftershocks is governed by afterslip.     

The 2010 Mw 8.8 Maule megathrust earthquake of Central Chile, monitored by GPS

Large earthquakes produce crustal deformation that can be quantified by geodetic measurements, allowing for the determination of the slip distribution on the fault. We used data from Global Positioning System (GPS) networks in Central Chile to infer the static deformation and the kinematics of the 2010 moment magnitude (M(w)) 8.8 Maule megathrust earthquake. From elastic modeling, we found a total rupture length of ~500 kilometers where slip (up to 15 meters) concentrated on two main asperities situated on both sides of the epicenter. We found that rupture reached shallow depths, probably extending up to the trench. Resolvable afterslip occurred in regions of low coseismic slip. The low-frequency hypocenter is relocated 40 kilometers southwest of initial estimates. Rupture propagated bilaterally at about 3.1 kilometers per second, with possible but not fully resolved velocity variations.     

Coseismic and Postseismic Fault Slip for the 17 August 1999, M = 7.5, Izmit, Turkey Earthquake

We use Global Positioning System (GPS) observations and elastic half-space models to estimate the distribution of coseismic and postseismic slip along the Izmit earthquake rupture. Our results indicate that large coseismic slip (reaching 5.7 meters) is confined to the upper 10 kilometers of the crust, correlates with structurally distinct fault segments, and is relatively low near the hypocenter. Continued surface deformation during the first 75 days after the earthquake indicates an aseismic fault slip of as much as 0.43 meters on and below the coseismic rupture. These observations are consistent with a transition from unstable (episodic large earthquakes) to stable (fault creep) sliding at the base of the seismogenic zone.     

Shallow dynamic overshoot and energetic deep rupture in the 2011 Mw 9.0 Tohoku-Oki earthquake

Strong spatial variation of rupture characteristics in the moment magnitude (M(w)) 9.0 Tohoku-Oki megathrust earthquake controlled both the strength of shaking and the size of the tsunami that followed. Finite-source imaging reveals that the rupture consisted of a small initial phase, deep rupture for up to 40 seconds, extensive shallow rupture at 60 to 70 seconds, and continuing deep rupture lasting more than 100 seconds. A combination of a shallow dipping fault and a compliant hanging wall may have enabled large shallow slip near the trench. Normal faulting aftershocks in the area of high slip suggest dynamic overshoot on the fault. Despite prodigious total slip, shallower parts of the rupture weakly radiated at high frequencies, whereas deeper parts of the rupture radiated strongly at high frequencies.     

Rupture process of the 2004 Sumatra-Andaman earthquake

The 26 December 2004 Sumatra-Andaman earthquake initiated slowly, with small slip and a slow rupture speed for the first 40 to 60 seconds. Then the rupture expanded at a speed of about 2.5 kilometers per second toward the north northwest, extending 1200 to 1300 kilometers along the Andaman trough. Peak displacements reached approximately 15 meters along a 600-kilometer segment of the plate boundary offshore of northwestern Sumatra and the southern Nicobar islands. Slip was less in the northern 400 to 500 kilometers of the aftershock zone, and at least some slip in that region may have occurred on a time scale beyond the seismic band.     

Slow earthquakes linked along dip in the Nankai subduction zone

We identified a strong temporal correlation between three distinct types of slow earthquakes distributed over 100 kilometers along the dip of the subducting oceanic plate at the western margin of the Nankai megathrust rupture zone, southwest Japan. In 2003 and 2010, shallow very-low-frequency earthquakes near the Nankai trough as well as nonvolcanic tremor at depths of 30 to 40 kilometers were triggered by the acceleration of a long-term slow slip event in between. This correlation suggests that the slow slip might extend along-dip between the source areas of deeper and shallower slow earthquakes and thus could modulate the stress buildup on the adjacent megathrust rupture zone.     

The 2011 magnitude 9.0 Tohoku-Oki earthquake: mosaicking the megathrust from seconds to centuries

Geophysical observations from the 2011 moment magnitude (M(w)) 9.0 Tohoku-Oki, Japan earthquake allow exploration of a rare large event along a subduction megathrust. Models for this event indicate that the distribution of coseismic fault slip exceeded 50 meters in places. Sources of high-frequency seismic waves delineate the edges of the deepest portions of coseismic slip and do not simply correlate with the locations of peak slip. Relative to the M(w) 8.8 2010 Maule, Chile earthquake, the Tohoku-Oki earthquake was deficient in high-frequency seismic radiation--a difference that we attribute to its relatively shallow depth. Estimates of total fault slip and surface secular strain accumulation on millennial time scales suggest the need to consider the potential for a future large earthquake just south of this event.     

A silent slip event on the deeper Cascadia subduction interface

Continuous Global Positioning System sites in southwestern British Columbia, Canada, and northwestern Washington state, USA, have been moving landward as a result of the locked state of the Cascadia subduction fault offshore. In the summer of 1999, a cluster of seven sites briefly reversed their direction of motion. No seismicity was associated with this event. The sudden displacements are best explained by approximately 2 centimeters of aseismic slip over a 50-kilometer-by-300-kilometer area on the subduction interface downdip from the seismogenic zone, a rupture equivalent to an earthquake of moment magnitude 6.7. This provides evidence that slip of the hotter, plastic part of the subduction interface, and hence stress loading of the megathrust earthquake zone, can occur in discrete pulses.     

Magnitude of late quaternary left-lateral displacements along the north edge of tibet

Images taken by the earth observation satellite SPOT of the Quaternary morphology at 18 sites on the 2000-kilometer-long Altyn Tagh fault at the north edge of Tibet demonstrate that it is outstandingly active. Long-term, left-lateral strike-slip offsets of stream channels, alluvial terrace edges, and glacial moraines along the fault cluster between 100 and 400 meters. The high elevation of the sites, mostly above 4000 meters in the periglacial zone, suggests that most offsets resulted from slip on the fault since the beginning of the Holocene. These data imply that slip rates are 2 to 3 centimeters per year along much of the fault length and support the hypothesis that the continuing penetration of India into Asia forces Tibet rapidly toward the east.     

Slow earthquakes coincident with episodic tremors and slow slip events

We report on the very-low-frequency earthquakes occurring in the transition zone of the subducting plate interface along the Nankai subduction zone in southwest Japan. Seismic waves generated by very-low-frequency earthquakes with seismic moment magnitudes of 3.1 to 3.5 predominantly show a long period of about 20 seconds. The seismicity of very-low-frequency earthquakes accompanies and migrates with the activity of deep low-frequency tremors and slow slip events. The coincidence of these three phenomena improves the detection and characterization of slow earthquakes, which are thought to increase the stress on updip megathrust earthquake rupture zones.     

Episodic tremor and slip on the Cascadia subduction zone: the chatter of silent slip

We found that repeated slow slip events observed on the deeper interface of the northern Cascadia subduction zone, which were at first thought to be silent, have unique nonearthquake seismic signatures. Tremorlike seismic signals were found to correlate temporally and spatially with slip events identified from crustal motion data spanning the past 6 years. During the period between slips, tremor activity is minor or nonexistent. We call this associated tremor and slip phenomenon episodic tremor and slip (ETS) and propose that ETS activity can be used as a real-time indicator of stress loading of the Cascadia megathrust earthquake zone.     

Earth's free oscillations excited by the 26 December 2004 Sumatra-Andaman earthquake

At periods greater than 1000 seconds, Earth's seismic free oscillations have anomalously large amplitude when referenced to the Harvard Centroid Moment Tensor fault mechanism, which is estimated from 300- to 500-second surface waves. By using more realistic rupture models on a steeper fault derived from seismic body and surface waves, we approximated free oscillation amplitudes with a seismic moment (6.5 x 10(22) Newton.meters) that corresponds to a moment magnitude of 9.15. With a rupture duration of 600 seconds, the fault-rupture models represent seismic observations adequately but underpredict geodetic displacements that argue for slow fault motion beneath the Nicobar and Andaman islands.     

A high-frequency secondary event during the 2004 Parkfield earthquake

By using seismic records of the 2004 magnitude 6.0 Parkfield earthquake, we identified a burst of high-frequency seismic radiation that occurred about 13 kilometers northwest of the hypocenter and 5 seconds after rupture initiation. We imaged this event in three dimensions by using a waveform back-projection method, as well as by timing distinct arrivals visible on many of the seismograms. The high-frequency event is located near the south edge of a large slip patch seen in most seismic and geodetic inversions, indicating that slip may have grown abruptly at this point. The time history obtained from full-waveform back projection suggests a rupture velocity of 2.5 kilometers per second. Energy estimates for the subevent, together with long-period slip inversions, indicate a lower average stress drop for the northern part of the Parkfield earthquake compared with that for the region near its hypocenter, which is in agreement with stress-drop estimates obtained from small-magnitude aftershocks.     

Fine structure of the landers fault zone: segmentation and the rupture process

Observations and modeling of 3- to 6-hertz seismic shear waves trapped within the fault zone of the 1992 Landers earthquake series allow the fine structure and continuity of the zone to be evaluated. The fault, to a depth of at least 12 kilometers, is marked by a zone 100 to 200 meters wide where shear velocity is reduced by 30 to 50 percent. This zone forms a seismic waveguide that extends along the southern 30 kilometers of the Landers rupture surface and ends at the fault bend about 18 kilometers north of the main shock epicenter. Another fault plane waveguide, disconnected from the first, exists along the northern rupture surface. These observations, in conjunction with surface slip, detailed seismicity patterns, and the progression of rupture along the fault, suggest that several simple rupture planes were involved in the Landers earthquake and that the inferred rupture front hesitated or slowed at the location where the rupture jumped from one to the next plane. Reduction in rupture velocity can tentatively be attributed to fault plane complexity, and variations in moment release can be attributed to variations in available energy.     

Three-dimensional splay fault geometry and implications for tsunami generation

Megasplay faults, very long thrust faults that rise from the subduction plate boundary megathrust and intersect the sea floor at the landward edge of the accretionary prism, are thought to play a role in tsunami genesis. We imaged a megasplay thrust system along the Nankai Trough in three dimensions, which allowed us to map the splay fault geometry and its lateral continuity. The megasplay is continuous from the main plate interface fault upwards to the sea floor, where it cuts older thrust slices of the frontal accretionary prism. The thrust geometry and evidence of large-scale slumping of surficial sediments show that the fault is active and that the activity has evolved toward the landward direction with time, contrary to the usual seaward progression of accretionary thrusts. The megasplay fault has progressively steepened, substantially increasing the potential for vertical uplift of the sea floor with slip. We conclude that slip on the megasplay fault most likely contributed to generating devastating historic tsunamis, such as the 1944 moment magnitude 8.1 Tonankai event, and it is this geometry that makes this margin and others like it particularly prone to tsunami genesis.     

Ancient ice islands in salt lakes of the central andes

Massive blocks of freshwater ice and frozen sediments protrude from shallow, saline lakes in the Andes of southwestern Bolivia and northeastern Chile. These ice islands range up to 1.5 kilometers long, stand up to 7 meters above the water surface, and may extend out tens of meters and more beneath the unfrozen lake sediments. The upper surfaces of the islands are covered with dry white sediments, mostly aragonite or calcite. The ice blocks may have formed by freezing of the fresh pore water of lake sediments during the "little ice age." The largest blocks are melting rapidly because of possibly recent increases in geothermal heat flux through the lake bottom and undercutting by warm saline lake water during the summer.     

Variation of interplate fault zone properties with depth in the japan subduction zone

The depth dependence of physical properties along the Japan subduction zone interface was explored using teleseismic recordings of earthquake signals. Broadband body waves were inverted to determine the duration of rupture and source depth for 40 interplate thrust earthquakes located offshore of Honshu between 1989 and 1995. After scaling for differences in seismic moment, there is a systematic decrease in rupture duration with increasing depth along the subducting plate interface. This indicates increases in rupture velocity or stress drop with depth, likely related to variation in rigidity of sediments on the megathrust.     

A mechanical model for intraplate earthquakes: application to the new madrid seismic zone

We present a time-dependent model for the generation of repeated intraplate earthquakes that incorporates a weak lower crustal zone within an elastic lithosphere. Relaxation of this weak zone after tectonic perturbations transfers stress to the overlying crust, generating a sequence of earthquakes that continues until the zone fully relaxes. Simulations predict large (5 to 10 meters) slip events with recurrence intervals of 250 to 4000 years and cumulative offsets of about 100 meters, depending on material parameters and far-field stress magnitude. Most are consistent with earthquake magnitude, coseismic slip, recurrence intervals, cumulative offset, and surface deformation rates in the New Madrid Seismic Zone. Computed interseismic strain rates may not be detectable with available geodetic data, implying that low observed rates of strain accumulation cannot be used to rule out future damaging earthquakes.     

Co-seismic strike-slip and rupture length produced by the 2001 Ms 8.1 Central Kunlun earthquake

Field investigations show that the surface wave magnitude (Ms) 8.1 Central Kunlun earthquake (Tibetan plateau) of 14 November 2001 produced a nearly 400-kilometer-long surface rupture zone, with as much as 16.3 meters of left-lateral strike-slip along the active Kunlun fault in northern Tibet. The rupture length and maximum displacement are the largest among the co-seismic surface rupture zones reported on so far. The strike-slip motion and the large rupture length generated by the earthquake indicate that the Kunlun fault partitions its deformation into an eastward extrusion of Tibet to accommodate the continuing penetration of the Indian plate into the Eurasian plate.     

The great Sumatra-Andaman earthquake of 26 December 2004

The two largest earthquakes of the past 40 years ruptured a 1600-kilometer-long portion of the fault boundary between the Indo-Australian and southeastern Eurasian plates on 26 December 2004 [seismic moment magnitude (Mw) = 9.1 to 9.3] and 28 March 2005 (Mw = 8.6). The first event generated a tsunami that caused more than 283,000 deaths. Fault slip of up to 15 meters occurred near Banda Aceh, Sumatra, but to the north, along the Nicobar and Andaman Islands, rapid slip was much smaller. Tsunami and geodetic observations indicate that additional slow slip occurred in the north over a time scale of 50 minutes or longer.