The HYDRUS-based flow package for MODFLOW (the HPM or the HYDRUS pack-age) is an existing unsaturated zone flow package for MODFLOW. In MODFLOW with the HPM, the groundwater modeling domain is discretized into regular grids that can be combined into multiple zones based on similarities in soil hydrology, topographical characteristics, and the depth to the groundwater. Each of these zones is assigned one unsaturated soil profile (the HPM profile). In this model, after every MODFLOW time step, the flux at the bottom of the HPM profile is given as an input recharge flux to MODFLOW. MODFLOW simulates groundwater flow, and the water table depth at the end of the MODFLOW time step is assigned as the bottom boundary condition in the HPM profile. The current coupling algorithm assumes that the groundwater table in the HPM profile remains constant throughout the entire MODFLOW time step. This results in unrealistic sudden inflow and/or outflow fluxes at the bottom of the HPM profile after every time step. The objective of this study was to develop a methodology to eliminate the error in the determination of the recharge flux at the bottom of the HPM profile. This was achieved by updating or modifying the pressure head profile in the HPM profile after every MODFLOW time step. The effectiveness and the applicability of the new coupling algorithm were evaluated using different case studies. The new coupling algorithm is effective in eliminating unrealistic sudden variations in the bottom flux in the HPM profiles
