We computed the global wavefield excited by Model III using the spectral element method (SEM) (1) for a 3D Earth model composed of mantle model s20rts (1), model Crust 2.0 (2), and topography from ETOPO5. We show the predicted velocities and displacements on the Earth's surface as movies. The 3D simulations were used to calibrate the effect of 3D structure on the 1D waveforms used in the inversion for Model III. By comparing the fits of synthetics computed for different models to data that were not used in the inversions these synthetics can be used to distinguish between models.This was done qualitatively when developing model III to estimate the improvements between successive versions of the model. The waveforms for model III match the overall amplitude and directivity in observed seismograms (Fig. S1) over a wide

We determine coseismic and the first-month postseismic deformation associated with the Sumatra–Andaman earthquake of 26 December 2004 from near- field Global Positioning System (gps) surveys in northwestern Sumatra and along the Nicobar-Andaman islands, continuous and campaign gps measurements from Thailand and Malaysia, and in situ and remotely sensed observations of the vertical motion of coral reefs. The coseismic model shows that the Sunda subduction megathrust ruptured over a distance of about 1500 km and a width of less than 150 km, releasing a total moment of 6.7–7.0 × 1022 N m, equivalent to a magnitude Mw ∼9.15. The latitudinal distribution of released moment in our model has three distinct peaks at about 4° N, 7° N, and 9° N, which compares well to the latitudinal variations seen in the seismic inversion and of the analysis of radiated T waves. Our coseismic model is also consistent with interpretation of normal modes and with the amplitude of very-long-period surface waves. The tsunami predicted from this model fits relatively well the altimetric measurements made by the jason and topex satellites. Neither slow nor delayed slip is needed to explain the normal modes and the tsunami wave. The near-field geodetic data that encompass both coseismic deformation and up to 40 days of postseismic deformation require that slip must have continued on the plate interface after the 500-sec-long seismic rupture. The postseismic geodetic moment of about 2.4 × 1022 N m (Mw ∼8.8) is equal to about 30 ± 5% of the coseismic moment release. Evolution of postseismic deformation is consistent with rate-strengthening frictional afterslip. Online material: Summary of geodetic data used in this study.

40Large volcanic eruptions on Earth commonly occur with collapse of the roof of a crustal magma 41 reservoir, forming a caldera. Only a few such collapses occur per century and lack of detailed 42 observations has obscured insight on mechanical interplay between collapse and eruption. We use Calderas are 1 -100 km diameter depressions found in volcanic regions of Earth and other planets. basaltic andesite) intrusive activity and eruptions (2,(9)(10)(11)(12). 59The consensus from field and modelling studies is that caldera collapse progresses from initial 60 surface downsag to fault-controlled subsidence (1, 8, 13, 14). The pre-collapse topography is obtained by subtracting the subsidence observed at the surface. As we recorded the caldera subsidence mainly on the ice (Fig. 1, Fig. S1), we made corrections and (Fig. 3A). We therefore conclude that suggestions of a large increase in ice flow out of the caldera 147 during these events (25) cannot be fitted with our data. 148Bedrock subsidence exceeding 1 m occurred within an area of 110 km 2 that extended beyond the 149 pre-existing caldera (Fig. 1, Fig. S1). After termination of collapse the total subsidence at the pre-150 existing caldera rims amounted to 3 to 11 meters ( Fig. 1D and 1E). Using subglacial radio-echo GPS station in the center of the caldera (Fig. 1A), including the rate of vertical rate of ice surface Cumulative number of M>4 caldera earthquakes, with magnitude evolution colored in red, blue and 176 grey representing clusters on the southern rim, the northern rim and smaller clusters, respectively 177 (see Fig. S5). E) Cumulative seismic moment for M>4 caldera earthquakes. from analysis of subaerial gas measurements (Fig. 4). This depth concurs with our regional on FTIR and Multi-GAS measurements (24). 194Seismicity and subsurface structure 195 We used seismic data and Distinct Element Method (DEM) numerical modelling (24), to 196 characterize the deeper collapse structure as the reactivation of a steeply-inclined ring fault (Fig. 5). 197We mostly observed seismicity at depths of 0-9 km beneath the northern and southern caldera rims 198( Fig. 5B), with earthquakes being more numerous on the northern rim. This spatial pattern of 199 seismicity is consistent with fracturing above a deflating magma reservoir that was elliptical in (Fig. 5C, D). Our best fitting models had preexisting faults dipping out at 80-85¡ from the caldera 207 center on the north side and at 85-90¡ toward the caldera center on the south side. The modeled pre- 208existing faults lay at 1-2 km below the surface on the north side and 3-4 km on the south side. 209Modeling of a more complex fault geometry or the inclusion of greater material heterogeneity may 210 further improve the data fit, but presently lacks robust geophysical constraints. components of the observed earthquakes at B ‡rdarbunga. We, however, narrowed down on 222 plausible solutions by using the micro-earthquakes (Fig. 5A). The moment tensor solutions are well 223 constrained, but the inferred d...

The 2005 M w 8.6 Nias-Simeulue earthquake was caused by rupture of a portion of the Sunda megathrust offshore northern Sumatra. This event occurred within an array of continuous Global Positioning System (GPS) stations and produced measurable vertical displacement of the fringing coral reefs above the fault rupture. Thus, this earthquake provides a unique opportunity to assess the source characteristics of a megathrust event from the joint analysis of seismic data and near-field static co-seismic displacements. Based on the excitation of the normal mode data and geodetic data we put relatively tight constraints on the seismic moment and the fault dip, where the dip is determined to be 8Њ to 10Њ with corresponding moments of 1.24 ‫ן‬ 10 22 to 1.00 ‫ן‬ 10 22 N m, respectively. The geodetic constraints on slip distribution help to eliminate the trade-off between rupture velocity and slip kinematics. Source models obtained from the inversion of various combinations of the teleseismic body waves and geodetic data are evaluated by comparing predicted and observed long-period seismic waveforms (100-500 sec). Our results indicate a relatively slow average rupture velocity of 1.5 to 2.5 km/sec and long average rise time of up to 20 sec. The earthquake nucleated between two separate slip patches, one beneath Nias and the other beneath Simeulue Island. The gap between the two patches and the hypocentral location appears to be coincident with a local geological disruption of the forearc. Coseismic slip clearly tapers to zero before it reaches the trench probably because the rupture propagation was inhibited when it reached the accretionary prism. Using the models from joint inversions, we estimate the peak ground velocity on Nias Island to be about 30 cm/sec, an order of magnitude slower than for thrust events in continental areas. This study emphasizes the importance of utilizing multiple datasets in imaging seismic ruptures.

No abstract

No abstract

International audienc

S U M M A R YPrevious analyses of geological and geodetic data suggest that the obliquely compressive relative motion across the Caribbean-North America plate boundary in Hispaniola is accommodated through strain partitioning between near-vertical transcurrent faults on land and low-angle thrust faults offshore. In the Dominican Republic, earthquake focal-mechanism geometries generally support this interpretation. Little information has been available about patterns of seismic strain release in Haiti, however, due to the small numbers of moderateto-large earthquakes occurring in western Hispaniola during the modern instrumental era. Here, we analyse the damaging M W = 7.0 earthquake that occurred near Port au Prince on 2010 January 12 and aftershocks occurring in the four months following this event, to obtain centroid-moment-tensor (CMT) solutions for 50 earthquakes with magnitudes as small as M W = 4.0. While the 2010 January main shock exhibited primarily strike-slip motion on a steeply dipping nodal plane (strike=250 • , dip=71 • and rake=22 • ), we find that nearly all of the aftershocks show reverse-faulting motion, typically on high-angle (30 • -45 • ) nodal planes. Two small aftershocks (M W 4.5 and 4.6), located very close to the main shock epicentre, show strike-slip faulting with geometries similar to the main shock. One aftershock located off the south coast of Haiti shows low-angle thrust faulting. We also examine earthquakes occurring in this region from 1977-2009; successful analysis of four such events provides evidence for both strike-slip and reverse faulting. The pattern of seismic strain release in southern Haiti thus indicates that partitioning of plate motion between transcurrent and reverse structures extends far west within Hispaniola. While we see limited evidence for low-angle underthrusting offshore, most reverse motion appears to occur on high-angle fault structures adjacent to the Enriquillo fault. Our results highlight the need to incorporate seismogenic slip on compressional structures into hazard assessments for southern Haiti.