In addition to this new boundary condition, [@moench1995] included a finite radius pumping well with wellbore storage, conceptually similar to how @papadopulosandcooper1967 modified the solution of @theis1935. In all other respects, his definition of the problem was similar to [@neuman1974].

Moench’s solution resulted in improved fits to experimental data and produced more realistic estimates of specific yield [@moenchetal2001], including the use of multiple delay parameters $\alpha_m$ [@moench2003]. @moenchetal2001 used with $M = 3$ to estimate hydraulic parameters in the unconfined aquifer at Cape Cod. They showed that $M=3$ enabled a better fit to the observed drawdown data than obtained by $M<3$ or the model of [@neuman1974]. Similar to the parameter $\alpha $ in Boulton’s model, the physical meaning of $\alpha _m$ is not clear.

Unconfined Solutions Considering Unsaturated Flow
=================================================

As an alternative to linearizing the water table condition of @boulton1954a, the unsaturated zone can be explicitly included. The non-linearity of unsaturated flow is substituted for the non-linearity of . By considering the vadose zone, the water table is internal to the domain, rather than a boundary condition. The model-data misfit in Figure \[fig:capecod\] at “late intermediate” time is one of the motivations for considering the mechanisms of delayed yield 