“…The model domain size and mesh discretization should be carefully chosen to best represent the system scale and to allow accurate reproduction of heterogeneities and hydraulic and thermal gradients while keeping in mind the study's objectives. For example, recent publications have shown the importance of shallow supra‐permafrost groundwater on cold‐regions hydrology (Lamontagne‐Hallé et al, 2018; O'Connor, Cardenas, Neilson, Nicholaides, & Kling, 2019; Rey, Walvoord, Minsley, Rover, & Singha, 2019; Toohey, Herman‐Mercer, Schuster, Mutter, & Koch, 2016). To simulate these local processes included in a larger groundwater system, while limiting computational time, a finer mesh size near the surface and a coarser mesh size at greater depths can be used (e.g., McKenzie & Voss, 2013).…”