“…Heat transfer will occur through thermal conduction and convection of fluids circulating in porous and fractured host rocks (e.g., Norton & Taylor, 1979; Parmentier & Schedl, 1981; Turcotte & Schubert, 2002). The occurrence of convection is more evident in metallogenic districts where fluids rising along faults and fracture zones or highly permeable rocks lead to localized mineralization processes (e.g., Fu et al., 2010; Luijendijk, 2019; Sánchez et al., 2021). To explore and quantify the transient thermal effect of magmatic intrusions on host rocks (i.e., the spatial distribution, magnitude, and timing of the perturbation with and without fluids), we designed a forward, 3D numerical thermal model using the TOUGH2 code (Transport Of Unsaturated Groundwater and Heat; W. M. Kissling, 1995; Pruess, 1991, https://tough.lbl.gov/).…”