The growth of normal faults in mechanically layered sequences is numerically modelled using three-dimensional Distinct Element Method (DEM) models, in which rock comprises an assemblage of bonded spherical particles. Faulting is induced by movement on a pre-defined normal fault at the model base whilst a constant confining pressure is maintained by applying forces to particles lying at the model top. The structure of the modelled fault zones and its dependency on confining pressure, sequence (net:gross) and fault obliquity are assessed using various new techniques that allow (a) visualization of faulted horizons, (b) quantification of throw partitioning and (c) determination...