“…However, it has been shown that a two-way coupling approach is more applicable for temporally protracted perturbations (Neuzil, 2003;Rutqvist et al, 2002), where subsurface strain affects hydraulic rock properties (e.g., κ, ϕ) which in turn govern the flow behaviour and in turn stresses and strains. Similar to studies by e.g., Hutnak et al (2009); Cur- Neither of our simulations account for the temperature dependence of parameters such as permeability (Ikard & Revil, 2014), fluid properties (Arens et al, 2020) or elastic parameters (Head et al, 2021), all of which have an effect on geophysical anomalies modelled in our study; a dedicated analysis is required to assess this influence. Although we neglect thermoelectric processes caused by strong thermal gradients (Corwin & Hoover, 1979;Fitterman & Corwin, 1982) in the HTU simulations, we find that for a maximum temperature change of 0.18 • C (unrest III) at the plateau, the thermoelectric potential (TEP) is ± 0.3 mV and 0.1 mV using a thermoelectric coupling coefficient of ± 1.5 mV/ • C and ± 0.5 mV/ • C (Revil & Mahardika, 2013;Ikard & Revil, 2014), respectively.…”