Finite Element (FE) modeling under plane stress condition is used to analyze the fault type variation with depth along and around the San Andreas Fault (SAF) zone. In this simulation elastic rheology was used and was thought justifiable as the variation in depth from 0.5 km to 20 km was considered. Series of calculations were performed with the variation in domain properties. Three types of models were created based on simple geological map of California, namely, 1) single domain model considering whole California as one homogeneous domain, 2) three domains model including the North American plate, Pacific plate, and SAF zone as separate domains, and 3) Four domains model including the three above plus the Garlock Fault zone. Mohr-Coulomb failure criterion and Byerlee’s law were used for the calculation of failure state. All the models were driven by displacement boundary condition imposing the fixed North American plate and Pacific plate motion along N34°W vector up to the northern terminus of SAF and N50°E vector motion for the subducting the Gorda and Juan de Fuca plates. Our simulated results revealed that as the depth increased, the fault types were generally normal, and at shallow depth greater strike slip and some thrust faults were formed. It is concluded that SAF may be terminated as normal fault at depth although the surface expression is clearly strike slip. Finite Element (FE) modeling under plane stress condition is used to analyze the fault type variation with depth along and around the San Andreas Fault (SAF) zone. In this simulation elastic rheology was used and was thought justifiable as the variation in depth from 0.5 Series of calculations were performed with the variation in domain properties. Three types of models were created based on simple geological map of California, namely, 1) single domain model considering whole California as one homogeneous domain, 2) Three domains model including the North American plate, Pacific plate, and SAF zone as separate domains, and 3) Four domains model including the three above plus the Garlock Fault zone. Mohr-Coulomb failure criterion and Byerlee’s law were used for the calculation of failure all the models were driven by displacement boundary condition imposing the fixed North American plate and Pacific plate motion along N34 ° W vector up to the northern t erminus of SAF and N50 ° E vector motion for the subducting the Gorda and Juan de Fuca plates. Our simulated results revealed that as the depth increased, the fault types were generally normal, and at shallow depth greater strike slip and some thrust faults were formed It is said that SAF may be terminated as normal fault at depth although the surface expression is clearly strike slip.