Abstract. Field observations of two overlapping normal faults and associated deformation document features common to many normal-fault relay zones: a topographic ramp between the fault segments, tapering slip on the faults as they enter the overlap zone, and associated fracturing, especially at the top of the ramp. These observations motivate numerical modeling of the development of a relay zone. A three-dimensional boundary element method numerical model, using simple fault-plane geometries, material properties, and boundary conditions, reproduces the principal characteristics of the observed fault scarps. The model, with overlapping, semicircular fault segments under orthogonal extension, produces a region of high Coulomb shear stress in the relay zone that would favor fault linkage at the center to upper relay ramp. If the fault height is increased, the magnitude of the stresses in the relay zone increases, but the position of the anticipated linkage does not change. The amount of fault overlap changes the magnitude of the Coulomb stress in the relay zone: the greatest potential for fault linkage occurs with the closest underlapping fault tips. Ultimately, the mechanical interaction between segments of a developing normal-fault system promote the development of connected, zigzagging fault scarps.