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.

In an active volcanic arc, magmatically sourced fluids are channeled through the brittle crust by structural features. This interaction is observed in the Andean volcanic mountain belt, where volcanoes, geothermal springs and the locations of major mineral deposits coincide with NNE-striking, convergent margin-parallel faults and margin-oblique, NW/SE-striking Andean Transverse Faults (ATF). The Tinguiririca and Planchón-Peteroa volcanoes in the Andean Southern Volcanic Zone (SVZ) demonstrate this relationship, as both volcanic complexes and their spatially associated thermal springs show strike alignment to the outcropping NNE oriented El Fierro Thrust Fault System. This study aims to constrain the 3D architecture of this fault system and its interaction with volcanically sourced hydrothermal fluids from a combined magnetotelluric (MT) and seismicity survey. The 3D conductivity model and seismic hypocenter locations show correlations between strong conductivity contrasts and seismic clusters in the top 10km of the crust. This includes a distinct WNW-striking seismogenic feature which has characteristics of the ATF domains. As the surveyed region is characterized by high heat flow regimes, volcanic activity and hydrothermal systems related to the volcanic arc, the conductivity contrast suggests that magmatically derived fluids meet an impenetrable barrier, most likely the sealed core of the fault. The resulting increase in hydrostatic fluid pressure facilitates seismic activity on this WNW oriented structure. These results provides the first observation of the mechanism behind the reactivation and seismogenesis of ATF. The study also uncovers the role of the ATF the compartmentalization of magmatic-derived fluids that accumulate to form hydrothermal reservoirs in the SVZ.