We use geodetic data to illuminate the distribution of horizontal flow beneath the California lithosphere. Major earthquakes in California have occurred above the regions of current plastic strain accumulation.

We explain the super-cycles of Mw 7 and Mw 9 earthquakes at the Japan trench with a numerical models of seismic cycles in a viscoelastic lithosphere.

We develop new numerical models of seismic cycles to account for Coulomb stress interaction between parallel faults.

We identify the non-dimensional parameters that control the rupture style and recurrence patterns of fault dynamics and illustrate the range of phenomena expected for a single fault.

We develop a constitutive framework for rate- and state-dependent friction whereby strength depend on the real area of contact, which is modulated by the flattening and erosion of microasperities.

We provide Green's functions to model the deformation caused by distributed deformation, such as viscoelastic flow or poroelastic rebound.

We expand the boundary integral method to incorporate the viscoelastic response of the asthenosphere at subduction zones.

We use numerical models of fault dynamics to simulate the irregular sequence of moment magnitude 6.0 Parkfield earthquakes on the San Andreas Fault.