Three-dimensional mechanical simulations of the San Andreas fault (SAF) within the Coachella Valley in California produce deformation that match geologic observations and demonstrate the impact of fault geometry on uplift patterns. Most models that include the Coachella Valley segment of the SAF have assumed a vertical orientation, but recent studies suggest that this segment dips 60-70° northeast. We compare models with varied fault geometry and evaluate how well they reproduce observed uplift patterns. Our model with a dipping SAF matches geologic observations, while models containing a vertical fault do not. This suggests that the active Coachella Valley segment of the SAF dips 60-70° northeast. Since ~1.5 Ma, the SAF in this region has undergone a major reorganization that entailed initiation of the San Jacinto fault and termination of slip on the West Salton detachment fault. The trace of the SAF itself has also evolved, with several shifts in activity through the San Gorgonio Pass. Despite a rich geologic record of these changes, the mechanisms that controlled abandonment of faults, initiation of new strands, and shifting loci of uplift are poorly understood. We model snapshots in time through the evolution of the fault system, and assess the mechanical viability of our snapshots by comparison with uplift patterns inferred from the stratigraphic record. Model results are compared with vertical axis rotation. We examine incipient faulting using maps of strain energy density, and explore changes to the mechanical efficiency of the system to better understand the evolution of this fault system.

Numerical simulations of a 15° restraining bend analog claybox experiment include considering the fault geometry, rheology, and boundary conditions. The numerical models show that a growing fault from an analog experiment propagates at depth rather than at the surface and is exposed in later stages of the experiment, and that the wet kaolin clay from the analog experiment is partially decoupled from the steel plate. The numerical models provide the stresses to predict accurate fault growth from the analog experiment and provide the evolution of external work within the fault system. The external work from the numerical models decrease as faults continue to grow, which agrees with the continuously increasing kinematic efficiency within the analog experiment. Three-dimensional mechanical models are used to simulate the southern San Andreas fault. These models show that incorporating fault interaction, time since last earthquake rupture, and nearby earthquakes affects the stress state along a fault. Absolute shear tractions are calculated by multiplying time since last earthquake rupture with the simulated interseismic stressing rates for each fault strand. From our multi-cycle model, fault interaction affects local normal stressing rates so that the stresses are not relieved in between earthquakes. We provide our absolute shear tractions and scale our multi-cycle normal stressing rates to be near to failure so that dynamic rupture modelers from University of California, Riverside use our results to simulate earthquake propagation for the complex fault region of the San Gorgonio Pass.