In this dissertation I used Interferometric Synthetic Aperture Radar (InSAR) and Global Positioning System (GPS) to recover crustal deformation caused by earthquake cycle processes. The studied areas span three different types of tectonic boundaries: a continental thrust earthquake (M7.9 Wenchuan, China) at the eastern margin of the Tibet plateau, a mega-thrust earthquake (M8.8 Maule, Chile) at the Chile subduction zone, and the interseismic deformation of the San Andreas Fault System (SAFS). A new L-band radar onboard a Japanese satellite ALOS allows us to image high-resolution surface deformation in vegetated areas, which is not possible with older C-band radar systems. In particular, both the Wenchuan and Maule InSAR analyses involved L-band ScanSAR interferometry which had not been attempted before. I integrated a large InSAR dataset with dense GPS networks over the entire SAFS. The integration approach features combining the long-wavelength deformation from GPS with the short-wavelength deformation from InSAR through a physical model. The recovered fine-scale surface deformation leads us to better understand the underlying earthquake cycle processes. The geodetic slip inversion reveals that the fault slip of the Wenchuan earthquake is maximum near the surface and decreases with depth. The coseismic slip model of the Maule earthquake constrains the down-dip extent of the fault slip to be at 45 km depth, similar to the Moho depth. I inverted for the slip rate on 51 major faults of the SAFS using Green's functions for a 3-dimensional earthquake cycle model that includes kinematically prescribed slip events for the past earthquakes since the year 1000. A 60 km thick plate model with effective viscosity of 1019 Pa · s is preferred based on the geodetic and geological observations. The slip rates recovered from the plate models are compared to the half-space model. The InSAR observation reveals that the creeping section of the SAFS is partially locked. This high-resolution deformation model will refine the moment accumulation rates and shear strain rates, which are not well resolved by previous models.

These proceedings contain a selection of peer-reviewed papers presented at the International Symposium on Geodesy for Earthquake and Natural Hazards (GENAH), Matsushima, Japan, 22-26 July, 2014. The scientific sessions focused on monitoring temporal and spatial changes in Earth's lithosphere and atmosphere using geodetic satellite systems, high rate GNSS as well as high resolution imaging (InSAR, Lidar). Researchers in various fields of geodesy discussed the role of geodesy in disaster mitigation and how groups with different techniques can collaborate toward such a goal.

Geodetic observations provide kinematic constraints on the behavior of tectonically active fault systems. Estimates of earthquake cycle activity derived from these constraints may depend on modeling assumptions and/or regularization of a geodetic inverse problem, which is often poorly conditioned. Common model assumptions may affect kinematic solutions and conclusions about physical properties of faults and fault zones. For example, within a geometrically complex fault system, parameterization of nearby faults may affect slip estimates on an individual fault. In addition, fault slip models are often regularized by assuming that slip varies smoothly in space, which may artificially smear slip estimates beyond physical boundaries. As an alternative to smooth regularization, the applied mathematics field of compressed sensing provides a suite of methods for recovering sparse solutions. Applied to GPS observations of the 2011 Tohoku earthquake, compressed sensing algorithms enable imaging of spatially localized slip during and following the earthquake, and identification of a sharp boundary between coseismic and postseismic slip. Similar algorithms recover quantized solutions and may be applied to models of plate boundary deformation. Beginning with a dense array of tectonic micro-plates bounded by mapped faults in North America, these methods can be used to detect coherent motions of groups of micro-plates behaving as larger active blocks, effectively quantifying the complexity of North America plate boundary deformation. By improving our ability to identify and compare kinematic constraints on earthquake cycle processes, we are able to characterize the spectrum of earthquake cycle behaviors and gain a deeper understanding of earthquake phenomenology and physics.

In the last decade of the 20th century, there has been great progress in the physics of earthquake generation; that is, the introduction of laboratory-based fault constitutive laws as a basic equation governing earthquake rupture, quantitative description of tectonic loading driven by plate motion, and a microscopic approach to study fault zone processes. The fault constitutive law plays the role of an interface between microscopic processes in fault zones and macroscopic processes of a fault system, and the plate motion connects diverse crustal activities with mantle dynamics. An ambitious challenge for us is to develop realistic computer simulation models for the complete earthquake process on the basis of microphysics in fault zones and macro-dynamics in the crust-mantle system. Recent advances in high performance computer technology and numerical simulation methodology are bringing this vision within reach. The book consists of two parts and presents a cross-section of cutting-edge research in the field of computational earthquake physics. Part I includes works on microphysics of rupture and fault constitutive laws, and dynamic rupture, wave propagation and strong ground motion. Part II covers earthquake cycles, crustal deformation, plate dynamics, and seismicity change and its physical interpretation. Topics in Part II range from the 3-D simulations of earthquake generation cycles and interseismic crustal deformation associated with plate subduction to the development of new methods for analyzing geophysical and geodetical data and new simulation algorithms for large amplitude folding and mantle convection with viscoelastic/brittle lithosphere, as well as a theoretical study of accelerated seismic release on heterogeneous faults, simulation of long-range automaton models of earthquakes, and various approaches to earthquake predicition based on underlying physical and/or statistical models for seismicity change.