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.

The subduction zone in northern Chile is a well-identified seismic gap that last ruptured in 1877. The moment magnitude (Mw) 8.1 Iquique earthquake of 1 April 2014 broke a highly coupled portion of this gap. To understand the seismicity preceding this event, we studied the location and mechanisms of the foreshocks and computed Global Positioning System (GPS) time series at stations located on shore. Seismicity off the coast of Iquique started to increase in January 2014. After 16 March, several Mw > 6 events occurred near the low-coupled zone. These events migrated northward for ~50 kilometers until the 1 April earthquake occurred. On 16 March, on-shore continuous GPS stations detected a westward motion that we model as a slow slip event situated in the same area where the mainshock occurred.

The 1992 Landers, California, earthquake (magnitude 7.3) was modeled as the propagation of a spontaneous rupture controlled by a realistic prestress distribution with the use of a three-dimensional finite-difference method. The dynamic rupture reproduces the general slip pattern used to compute the initial stress distribution and generates near-fault ground motions at the surface similar to observations. The simulated rupture propagates on the fault along a complex path with highly variable speed and rise time, changing the magnitude and pattern of the stress significantly. The method provides the framework to estimate earthquake rupture parameters from recorded seismic and geodetic data.