Normal faults growing in the Earth's crust are subject to the effects of an increasing frictional resistance to slip caused by the increasing lithostatic load with depth. We use three-dimensional (3-D) boundary element method numerical models to evaluate these effects on planar normal faults with variable elliptical tip line shapes in an elastic solid. As a result of increasing friction with depth, normal fault slip maxima for a single slip event are skewed away from the fault center toward the upper fault tip. There is a correspondingly greater propagation tendency at the upper tip. However, the tall faults that would result from such a propagation tendency are generally not observed in nature. We show how mechanical interaction between laterally stepping fault segments significantly competes with the lithostatic loading effect in the evolution of a normal fault system, promoting lateral propagation and possibly segment linkage. Resultant composite faults are wider than they are tall, resembling both 3-D seismic data interpretations and previously documented characteristics of normal fault systems. However, this effect may be greatly complemented by the influence of a heterogeneous stratigraphy, which can control fault nucleation depth and inhibit fault propagation across the mechanical layering. Our models demonstrate that although lithostatic loading may be an important control on fault evolution in relatively homogeneous rocks, the contribution of lithologic influences and mechanical interaction between closely spaced, laterally stepping faults may predominate in determining the slip behavior and propagation tendency of normal faults in the Earth's crust. ¿ 1999 American Geophysical Union