The spatially integrated model (SIM) which can evaluate temporal variation of groundwater quality is proposed in the stream-aquifer setting entered by nonpoint source contaminants. And the developed SIM included unsaturated soil zone and was tested against the spatially distributed model (SDM) of the coupled advection-dispersion and Richards equations for the various hydrologic and aquifer simulating conditions. The result of the comparison showed that the average concentration responses of saturated aquifer and groundwater outflow between the SIM and the SDM was in good agreement, except for the case of the large dispersivity ratio and thick aquifer system. And it is shown that for the cases of the large dispersivity ratio and thick aquifer system the performance of the nonlinear SIM is better than that of the linear SIM for evaluating the average concentration of groundwater outflow response. Keywords : aquifer-stream interaction, groundwater quality modeling, nonpoint source pollution, unsaturated flow