Introduction
BackgroundShallow groundwater supplies a majority of streamflow to most watersheds (Beck et al., 2013) and is the primary source of streamflow during seasonal low flows (Smakhtin, 2001). This mobile water resource provides an essential buffer to changes in temperature, nutrients, and precipitation (Ficklin et al., 2015) and sustains evapotranspiration Yang et al., 2011) during low flow periods. Even in snow-dominated climates, the hydrologic response to snowmelt is mediated by groundwater (Enzminger et al., 2019) with important implications for the response of these systems to global warming (Tague & Grant, 2009). Because shallow groundwater contributions to low flows are disproportionately sensitive to changes in near-term (i.e., years and decades) climate signals (Hare et al., 2021), it is essential that our hydrologic models are able to predict groundwater contributions to low flows accurately. Unfortunately, accurate simulation of low flows and groundwater contributions to them has proven to be notoriously difficult in distributed, physically based land surface models (LSMs;Clark et al., 2015, Holtzman et al., 2020 as well as lumped, conceptual "bucket" models (i.e., rainfall-runoff models;Fowler et al., 2020).Early LSMs described low flows as dependent on one-dimensional drainage below a soil column (Clark et al., 2015). This "free" drainage led to well-documented inaccuracies: too-fast drainage during wet periods, underestimates of seasonal storage, and the cessation of ET during short dry periods (