Earthquake swarms commonly occur in upper-crustal hydrothermal-magmatic systems and activate mesh-like fault-fracture networks at zone of fault complexity. How these networks develop through space and time along seismic faults is poorly constrained in the geological record. Here, we describe a spatially dense array of small-displacement (< 1.5 m) epidote-rich fault-veins within granitoids, occurring at the intersections of subsidiary faults with the exhumed seismogenic Bolfin Fault Zone (Atacama Fault System, Northern Chile). Epidote faulting and veining occurred at 3-7 km depth and 200-300 °C ambient temperature. At distance ≤ 1 cm to fault-veins, the magmatic quartz of the wall-rock shows (i) thin (<10- µm-thick) interlaced deformation lamellae, and (ii) crosscutting quartz-healed veinlets. The epidote-rich fault-veins (i) include clasts of deformed magmatic quartz, with deformation lamellae and quartz-healed veinlets, and (ii) record cyclic events of extensional-to-hybrid veining and either aseismic and seismic shearing. Deformation of the wall-rock quartz is interpreted to record the large stress perturbations associated with the rupture propagation of small earthquakes. Instead, dilation and shearing forming the epidote-rich fault-veins are interpreted to record the later development of a mature and hydraulically-connected fault-fracture system. In this latter stage, the fault-fracture system cyclically ruptured due to fluid pressure fluctuations, possibly correlated with swarm-like earthquake sequences.

How major crustal-scale seismogenic faults nucleate and evolve in the crystalline basement represents a long-standing, but poorly understood, issue in structural geology and fault mechanics. Here, we address the spatio-temporal evolution of the Bolfin Fault Zone (BFZ), a >40-km-long exhumed seismogenic splay fault of the 1000-km-long strike-slip Atacama Fault System. The BFZ has a sinuous fault trace across the Mesozoic magmatic arc of the Coastal Cordillera (Northern Chile) and formed during the oblique subduction of the Aluk plate beneath the South American plate. Seismic faulting occurred at 5-7 km depth and ≤ 310 °C in a fluid-rich environment as recorded by extensive propylitic alteration and epidote-chlorite veining. Ancient (125-118 Ma) seismicity is attested by the widespread occurrence of pseudotachylytes. Field geological surveys indicate nucleation of the BFZ on precursory geometrical anisotropies represented by magmatic foliation of plutons (northern and central segments) and andesitic dyke swarms (southern segment) within the heterogeneous crystalline basement. Seismic faulting exploited the segments of precursory anisotropies that were favorably oriented with respect to the long-term far-stress field associated with the oblique ancient subduction. The large-scale sinuous geometry of the BFZ resulted from hard linkage of these anisotropy-pinned segments during fault growth.