Despite two decades of major advances in the field of thermochronological modeling, state-of-the-art numerical implementations still rely mostly on burial and exhumation processes to explain radiometric measurements. Even though such an approach has proved valuable, failing to account for other first-order geological variables has led to misinterpretations and therefore, calls for a refinement. In this study a new version of the Fetkin (finite element temperature kinematics, Ecopetrol) program is presented. Its new algorithm couples time-dependent hydrological and thermal calculations, thus rendering thermochronological ages that, instead of being solely dependent on the kinematical evolution of a system, conditioning by the fluid flow is also present. In contrast with previous thermochronological models, this work considers the influence of effective stress on rock properties (porosity and permeability) and therefore, in thermal conductivity. Sensitivity analyses addressing relevant geological questions show not only the versatility of the code but also, new perspectives on forward low-temperature thermochronological modeling. Groundwater circulation through pure-sandstone settings produce colder thermal architectures than those obtained in impermeable domains. Differences in cooling ages from models with and without fluid circulation are up to 5 Myr. A 4-fold variation in thrusting rates (0.5 km/Myr to 2 km/Myr) produces a 15-Myr difference in cooling ages in models with fluid flow, which contrasts to much lower differences, only 2 Myr, in domains without (or minimal) fluid circulation. 2D thermal solutions in fold-bend-fold thrust belts composed of sandstones remain static despite substantial relief development by kinematic folding. A case-study from Western Argentina, in the Andean Precordillera, confirms the plausibility of the numerical algorithm here posed and raises new questions on the first-order thermal controls in settings under deformation.

尽管在热年代学建模领域已取得了二十年的重大进展，但最新的数字实现仍主要依靠掩埋和挖掘过程来解释辐射测量。即使已证明这种方法很有价值，但不能考虑其他一阶地质变量也会导致误解，因此需要改进。在这项研究中，提出了Fetkin（有限元温度运动学，Ecopetrol）程序的新版本。它的新算法结合了随时间变化的水文和热力计算，从而提供了热年代学年龄，而不是仅取决于系统的运动学演化，还存在通过流体流动进行调节的情况。与以前的热年代学模型相比，这项工作考虑了有效应力对岩石性质（孔隙度和渗透率）的影响，并因此影响了热导率。针对相关地质问题的敏感性分析不仅显示了规范的多功能性，而且还展示了正向低温热年代学建模的新观点。通过纯砂岩环境进行的地下水循环产生的热结构比在不可渗透区域获得的热结构更冷。带和不带流体循环的型号的冷却时间差异最大为5 Myr。在有流体流动的模型中，推力速率的4倍变化（0.5 km / Myr至2 km / Myr）在冷却年龄方面产生了15 Myr的差异，而在没有（或最小）的流体循环。尽管通过运动学折叠产生了明显的浮雕，但由砂岩构成的折弯-折叠推力带中的二维热解仍保持静态。在安第斯山脉的安第斯山脉中的阿根廷西部进行了一项案例研究，证实了此处提出的数值算法的合理性，并对变形环境中的一阶热控制提出了新的问题。