Lagged streamflow depletion due to pumping-induced stream drying: Incorporation into analytical streamflow depletion estimation methods

• Analytical depletion functions (ADFs) estimate streamflow depletion caused by pumping. • ADFs incorporating stream drying had strong agreement with observed streamflow. • Stream drying shifts timing of streamflow depletion due to hydrologic disconnection. • Modeling stream drying requires accurate estimates of non-depleted streamflow. • Strong ADF performance suggests potential use as a low-cost decision-support tool. Water management often requires accounting for reductions in streamflow caused by groundwater pumping (‘streamflow depletion’). Since streamflow depletion cannot be quantified from observational data, it is typically modeled. Analytical depletion functions (ADFs) are a low-cost, low-complexity approach for estimating streamflow depletion with utility for decision support, but ADFs adopt several simplifying assumptions, including an infinite supply of water within the stream. Here, we develop an approach to incorporate stream drying into ADFs to improve their estimation of streamflow and streamflow depletion. Using Scott Valley (California) as an example, we compare ADF results to observed streamflow data and the Scott Valley Integrated Hydrologic Model (SVIHM), a process-based numerical model. ADFs incorporating stream drying have strong agreement with observed streamflow and SVIHM results. Critically, ADFs with drying can simulate a temporal shift in streamflow depletion that occurs when summer stream drying causes stream network disconnections and a substantial fraction of streamflow depletion is lagged until the fall/winter, when the stream network rewets. Estimates of what streamflow would have been without groundwater pumping are required to incorporate stream drying into ADFs, and we evaluate the ability of a statewide statistical model of unimpaired monthly streamflow (the California Natural Flows Database [CNFD]) to meet this need. ADFs using CNFD data simulate appropriate temporal dynamics but overestimate streamflow. This suggests that regional unimpaired flow estimates combined with local bias-correction could provide a mechanism to apply ADFs in watersheds without local numerical models.