“…Numerical models solving the Richards equation are often limited to small‐scale simulations (i.e., < 100 m 2 ) and to one‐dimensional (1D) flow in the shallow unsaturated zone (Scanlon et al., 2002), whereas other types of models may be used to assess the spatial and temporal variability of groundwater recharge at larger scales if the main hydrological variables, water table level and discharge, are known with sufficient accuracy (Scanlon et al., 2002). Commonly, large‐scale models (i.e., watershed, aquifer, or regional scale) rely on a spatial representation within a GIS platform, which is used to generate the various spatial inputs to the model (e.g., land cover, slope, and soil type), and then on water balance calculations to create the temporal dimension of recharge (Croteau, Nastev, & Lefebvre, 2010; Dripps & Bradbury, 2007; Jyrkama& Sykes, 2007; King, Raiber, Cox, & Cendón, 2017; Melo & Wendland, 2017; Sophocleous, 2000; Westenbroeck, Kelson, Dripps, Hunt, & Bradbury, 2009). To be calibrated, even the simplest model requires at least sufficiently long‐term discharge monitoring (e.g., GR [Génie Rural] models; Perrin, Michel, & Andréassian, 2003).…”