A MODFLOW package to linearize stream depletion analysis

The conventional numerical method is computationally intensive and prone to numerical noises for stream depletion analyses using MODFLOW. In this study, a new MODFLOW package has been developed to improve the computational efficiency and reduce the noises for each simulation. Using the assumption of unchanged flow coefficients between the baseline and scenario runs, the nonlinear groundwater flow system is linearized for solving the flow equations. The new package has been successfully applied to a regional groundwater model in Nebraska. The results show that the numerical noises, commonly identified in conventional approach, have been significantly reduced and a twenty-fold speedup has been achieved for a regional groundwater model in Nebraska. The results suggest this package can be adapted to be a component of optimization tools for water management scenario analyses especially when a large number of scenario model runs are involved.